Impossible de trouver l'emplacement UTM de l'OTAN dans Sentinel-2


10

Voir les coordonnées 31.96212, -103.004715

Les convertisseurs UTM donnent ses coordonnées UTM 13/R/FR.

Un exemple de convertisseur se trouve ici: http://www.rcn.montana.edu/resources/converter.aspx

Mais ils sont nombreux et donnent des réponses similaires pour ces coordonnées.

Simultanément, dans l'ensemble de données Sentinel-2 ici http://sentinel-s2-l1c.s3-website.eu-central-1.amazonaws.com/#tiles/13/R/

Je ne trouve pas de FRsous-répertoire.

Dans google, cet emplacement est ici:

entrez la description de l'image ici

Et trouver le même endroit dans le navigateur d'images Sentinel que je vois, cette tuile est différente

entrez la description de l'image ici

qui signifie par13/S/FR exemple la même UTMbande carrée mais différente.

Comment est-ce possible?

METTRE À JOUR

KML avec les tuiles Sentinel-2 signale également la Stuile à un emplacement donné

entrez la description de l'image ici

MISE À JOUR 2

D'après cette photo

entrez la description de l'image ici

prise d'ici , la FRplace est située à moitié en Szone UTM et à moitié en Rzone. De toute évidence, la plupart des convertisseurs automatiques attribuent ce carré à la Rzone, tandis que Sentinel-2 en tient compte S.

Y a-t-il une vérité ici?

MISE À JOUR 3

Le code Python simple, extrait d'ici 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)

revient également 13R.

Est-ce une erreur dans les données Sentinel-2 ou quoi?



C'est le cas S/FR, alors que les convertisseurs UTM donnent R/FR. Comment calculer l'emplacement si les convertisseurs UTM ne fonctionnent pas correctement?
Dims

La latitude est un peu moins de 32 degrés nord. Cela le place carrément dans la "bande" de latitude R. Sentinel-2 peut avoir une tuile en utilisant le point central de la tuile qui pourrait être dans la bande "S" à la place.
mkennedy

@mkennedy comment simuler cet algorithme à partir de coordonnées?
Dims

2
Vous pouvez également envisager de signaler cela à eosupport@copernicus.esa.int, car cela ressemble en effet à un comportement inattendu.
Kersten

Réponses:


1

En réponse à votre question de commentaire "comment simuler cet algorithme":

C'est une solution assez brute, mais facile à implémenter et devrait donner de bonnes performances:

  1. Utilisez l'un des convertisseurs UTM qui fonctionnent "comme prévu", en plaçant les coordonnées dans 13R.
  2. Vérifiez ensuite si le dossier existe dans la structure de données Sentinel 2. Si oui, vous avez terminé, hourra.

  3. Sinon, vérifiez les grilles UTM voisines et voyez si la tuile / dossier "FR" existe dedans. Étant donné qu'il y a des chevauchements partout, vous devez vérifier les 8 grilles environnantes.
    L'ordre de vérification le plus probable serait 13S, 13Q, 12R, 14R, 12S, 14S, 12Q, 14Q.
    Les quatre derniers pourraient être pertinents si vos coordonnées se trouvent dans les coins d'une zone UTM, mais sont hautement improbables.

Étant donné la façon dont Sentinel2 étiquette les tuiles, un seul des voisins devrait jamais avoir un tel dossier, ce qui vous garantit d'obtenir le bon fichier.

Toute autre solution géographiquement plus "correcte" impliquerait beaucoup plus de frais de calcul que je ne le pense justifié ici.

Et certainement, indiquez-le à l'équipe de l'ESA, comme l'a suggéré Kersten dans les commentaires. Je ne comprends vraiment pas pourquoi ils ont choisi un système organisationnel aussi inutilement alambiqué.


0

Article connexe ici

Ce qui a fonctionné pour moi, c'est d'utiliser le S2 KML fourni par l'ESA pour calculer toutes les tuiles qui se croisent avec mon AOI, puis rechercher ces tuiles dans AWS.

Ce KML semble fonctionner comme une définition de tous les identifiants de tuiles possibles générés par S2, éliminant de nombreuses options qui se chevauchent.

En regardant le KML (inspection visuelle uniquement, pas sûr à 100%), il me semble que dans le pire des cas, il faudrait rechercher 4 tuiles.

Ce serait bien d'avoir l'algorithme que l'ESA a utilisé pour définir le KML pour le rendre plus efficace.

En utilisant notre site, vous reconnaissez avoir lu et compris notre politique liée aux cookies et notre politique de confidentialité.
Licensed under cc by-sa 3.0 with attribution required.