Я сейчас работаю над проектом спутниковых данных. Основой следующей части является открытие данных спутника Sentinel-2 и сохранение визуального спектра в виде обычного изображения. Для этой задачи я использовал satpy.
Теперь я хочу обрезать плитку в меньшие изображения заданного размера. Чтобы обрезать, нужно
- , чтобы убедиться, что координаты (long, lat) для обрезки находятся в пределах плитки
- , обрезать изображение с заданными координатами long, lat
Я не смог найти первую часть. У кого-нибудь есть решение для этого. Я думаю, что этот ридер (safe_msi) не является функцией для этого?
Во-вторых, я хотел вызвать функцию обрезки класса Scene с примером:
(11.99900880451099, 48.96301655297934 , 12.080991195489009, 49.016968756890996)
Но я получаю ошибку:
NotImplementedError: Исходная проекция должна быть 'geos', если исходная / целевая проекции не равны.
, чего я не понял относительно:
https://satpy.readthedocs.io/en/latest/api/satpy.html#satpy .scene.Scene.crop
, в котором говорится:
ll_bbox (кортеж, список) - кортеж из 4 элементов, значения в градусах по долготе и долготе. Элементы (xmin, ymin, xmax, ymax), где X - долгота, а Y - широта.
Понимаю ли я что-то неправильно? Что я могу сделать, чтобы обрезать сцену?
Большое спасибо!