Не удалось найти местоположение НАТО UTM в Сентинел-2

10

Считать координаты 31.96212, -103.004715

UTM-конвертеры дают его UTM-координаты 13/R/FR.

Пример конвертера находится здесь: http://www.rcn.montana.edu/resources/converter.aspx

Но их много, и они дают похожие ответы для этих координат.

Одновременно в наборе данных Sentinel-2 здесь http://sentinel-s2-l1c.s3-website.eu-central-1.amazonaws.com/#tiles/13/R/

Я не могу найти FRподкаталог.

В Google это место находится здесь:

введите описание изображения здесь

И найдя то же самое место в браузере изображений Sentinel, я вижу, что плитка отличается

введите описание изображения здесь

который выступает за13/S/FR то есть то же самое , UTMи площади, но разные группы.

Как это возможно?

ОБНОВИТЬ

KML с тайлами Sentinel-2 также сообщает о Sтайле в указанном месте

введите описание изображения здесь

ОБНОВЛЕНИЕ 2

По этой картине

введите описание изображения здесь

взятая отсюда , FRплощадь находится наполовину в Sзоне UTM и наполовину в Rзоне. Очевидно, что большинство автоматических преобразователей присваивают этот квадрат Rзоне, а Sentinel-2 учитывает его как Sзону.

Здесь есть правда?

ОБНОВЛЕНИЕ 3

Простой код Python, взятый отсюда https://gis.stackexchange.com/a/224994/32207

bandVals = "CDEFGHJKLMNPQRSTUVWXX"

lon = 31.96212
lat = -103.004715

zone = int(lat + 186.0) / 6

if (lon >= 84.0):
    band = 'Y' if (lat < 0.0) else 'Z'
elif (lon <= -80.0):
    band = 'A' if (lat < 0.0) else 'B'
else:
    band = bandVals[int(lon + 80.0) / 8]

print '{:02d}{:s}'.format(zone,band)

также возвращается 13R.

Эта ошибка в данных Sentinel-2 или как?

Dims
источник
Это S/FRпока UTM конвертеры дают R/FR. Как рассчитать местоположение, если преобразователи UTM работают некорректно?
Димс
Значение широты чуть меньше 32 градусов северной широты. Это помещает это прямо в R полосу широты. Sentinel-2 может быть мозаичным, используя центральную точку плитки, которая вместо этого может находиться в полосе «S».
Mkennedy
@mkennedy, как смоделировать этот алгоритм, начиная с координат?
Димс
2
Вы также можете сообщить об этом по адресу eosupport@copernicus.esa.int, поскольку это действительно выглядит как неожиданное поведение.
Керстен

Ответы:

1

В ответ на ваш комментарий на вопрос «как смоделировать этот алгоритм»:

Это довольно грубое решение, но простое в реализации и должно дать хорошую производительность:

  1. Используйте любой из преобразователей UTM, которые работают «как положено», поместив координаты в 13R.
  2. Затем проверьте, существует ли папка в структуре данных Sentinel 2. Если да, то все готово, ура.

  3. Если нет, проверьте соседние сетки UTM и посмотрите, существует ли в них плитка / папка "FR". Поскольку повсюду есть перекрытия, вам придется проверить все окружающие 8 сеток.
    Наиболее вероятный порядок проверки будет 13S, 13Q, 12R, 14R, 12S, 14S, 12Q, 14Q.
    Последние четыре могут иметь значение, если ваши координаты лежат в углах зоны UTM, но маловероятно.

Учитывая способ, которым Sentinel2 маркирует плитки, такая папка должна быть только у одного из соседей, гарантируя, что вы получите правильный файл.

Любое другое, более географически более «правильное» решение потребовало бы намного больше вычислительных затрат, чем я считаю здесь оправданным.

И, безусловно, обязательно сообщите об этом команде ESA, как это было предложено Керстеном в комментариях. Я действительно не понимаю, почему они выбрали такую ​​излишне запутанную организационную систему.

Senshi
источник
0

Связанный пост здесь

Для меня было полезно использовать S2 KML, предоставленный ESA, для вычисления всех плиток, которые пересекаются с моим AOI, и затем искать эти плитки в AWS.

Этот KML, кажется, работает как определение всех возможных идентификаторов тайлов, генерируемых S2, устраняя множество перекрывающихся опций.

Смотря на KML (только визуальный осмотр, не уверен на 100%), мне кажется, что в худшем случае вам придется искать 4 плитки.

Было бы неплохо иметь алгоритм, который ESA использовал для определения KML, чтобы сделать это более эффективным.

Фред Липорас
источник