Я создаю инструмент, который загружает изображения из "Sentinel-2 на AWS"
У меня проблема с переводом координат Long / Lat в конкретную плитку Sentinel. Я знаю, что Sentinel 2 использует MGRS, и я попытался просто конвертировать Long / Lat в MGRS, и с низкой точностью я получил идентификатор плитки S2A. И это работает, но не всегда.
Например, Long / Lat 34.665,31.625 преобразуется в «36R XA» в MGRS, но плитка Sentinel 36RXA не существует.
Каков будет правильный способ определения плитки Sentinel с использованием координат Long / Lat?
SELECT * FROM [bigquery-public-data:cloud_storage_geo_index.sentinel_2_index] where north_lat < 50 and north_lat > 49 order by sensing_time LIMIT 10
, но вы упомянули AWS, поэтому, думаю, это не поможет.Ответы:
Разделение S2 в MGRS особенно сложно, когда данные находятся между зонами UTM и перекрываются тайлы, см. Рисунок ниже. Плитки слева направо: 20_M_QE, 20_M_RE, 21_M_SV и 21_M_TV:
Поиск, выполненный с помощью opensearch для региона, возвращает одно приобретение S2, разбитое на две плитки: 20_M_RE и 21_M_TV:
Поэтому, если вы заинтересованы в данных в представленной области тайла 21_M_SV, вам следует использовать 20_M_RE. Обратите внимание, что представление, используемое в примерах (шейп-файл, полученный из NGA), перекрытие между тайлами не представлено.
По-видимому, самым безопасным способом получения плиток является поиск в вашем AOI с использованием opensearch, а затем получение возвращенных плиток или использование инструмента sentinel2-util, как упомянуто bugmenot123.
Если вы решите отобразить Long / Lat на плитки, я рекомендую использовать KML, предоставленный ESA, в качестве определения сетки листов:
https://sentinel.esa.int/web/sentinel/missions/sentinel-2/data-products
источник
Такой инструмент уже существует. Он называется Sentinelsat, а источник доступен на GitHub . Он предлагает интерфейс командной строки и Python API. Он работает с Sentinel 1 и 2. Пространственный запрос основан на многоугольнике, а не на точке, но в остальном это именно то, что вам нужно.
РЕДАКТИРОВАТЬ: 1) вы можете вернуть идентификатор продукта (или список идентификаторов продуктов) с помощью функции запроса. А из командной строки вы можете опустить опцию -d, чтобы не запускать загрузку (и использовать опцию -f, чтобы создать след каждого изображения). 2) Что касается геоджона, вы можете увидеть ниже, что wkt может быть используется непосредственно вместо Как упомянул @ThingumaBob, вы можете создать крошечный многоугольник на основе ваших координат широты и долготы. Но на самом деле я протестировал его с точкой GeoJason, и он также работал нормально. Таким образом, нет необходимости конвертировать вашу точку в многоугольник.
которые могут быть скорректированы в
источник