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

Считать координаты 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 или как?

coordinate-system utm sentinel-2 mgrs Dims
источник

это не? sentinel-s2-l1c.s3.amazonaws.com/tiles/13/S/FR/2017/5/6/0/… разные группы - sentinel-s2-l1c.s3-website.eu-central-1.amazonaws .com / # плитки / 13 /…

Mapperz

Это S/FRпока UTM конвертеры дают R/FR. Как рассчитать местоположение, если преобразователи UTM работают некорректно?

Димс

Значение широты чуть меньше 32 градусов северной широты. Это помещает это прямо в R полосу широты. Sentinel-2 может быть мозаичным, используя центральную точку плитки, которая вместо этого может находиться в полосе «S».

Mkennedy

@mkennedy, как смоделировать этот алгоритм, начиная с координат?

Димс

Вы также можете сообщить об этом по адресу eosupport@copernicus.esa.int, поскольку это действительно выглядит как неожиданное поведение.

Керстен

Ответы:

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

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

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

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

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

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

Senshi
источник

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

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

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

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

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

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