Koordinaten berücksichtigen 31.96212, -103.004715
UTM-Konverter geben die UTM-Koordinaten an 13/R/FR
.
Beispielkonverter finden Sie hier: http://www.rcn.montana.edu/resources/converter.aspx
Aber es gibt viele von ihnen und sie geben ähnliche Antworten für diese Koordinaten.
Gleichzeitig im Sentinel-2-Datensatz hier http://sentinel-s2-l1c.s3-website.eu-central-1.amazonaws.com/#tiles/13/R/
Ich kann kein FR
Unterverzeichnis finden .
In Google ist dieser Ort hier:
Und wenn ich den gleichen Ort im Sentinel-Bildbrowser finde, den ich sehe, ist diese Kachel anders
was für13/S/FR
dh gleiche UTM
und quadratische, aber unterschiedliche Band steht.
Wie ist das möglich?
AKTUALISIEREN
KML mit Sentinel-2-Kacheln meldet auch S
Kacheln an einem bestimmten Ort
UPDATE 2
Nach diesem Bild
genommen von hier , das FR
ist Platz befindet sich die Hälfte in S
in UTM - Zone und eine halbe R
Zone. Offensichtlich weisen die meisten automatischen Konverter dieses Quadrat der R
Zone zu, während Sentinel-2 es für die S
Zone berücksichtigt .
Gibt es hier eine Wahrheit?
UPDATE 3
Der einfache Python-Code von hier /gis//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)
kehrt auch zurück 13R
.
Ist dieser Fehler in Sentinel-2-Daten oder was?
S/FR
, während UTM-Konverter gebenR/FR
. Wie berechnet man den Standort, wenn UTM-Konverter nicht richtig funktionieren?Antworten:
Als Antwort auf Ihre Kommentarfrage "Wie simuliere ich diesen Algorithmus?":
Dies ist eine ziemlich brutale Lösung, aber einfach zu implementieren und sollte eine gute Leistung bringen:
Überprüfen Sie anschließend, ob der Ordner in der Sentinel 2-Datenstruktur vorhanden ist. Wenn ja, sind Sie fertig, Hurra.
Wenn nicht, überprüfen Sie die benachbarten UTM-Gitter und prüfen Sie, ob die Kachel / der Ordner "FR" in ihnen vorhanden ist. Da es überall Überlappungen gibt, müssten Sie alle umliegenden 8 Gitter überprüfen.
Die wahrscheinlichste zu überprüfende Reihenfolge wäre 13S, 13Q, 12R, 14R, 12S, 14S, 12Q, 14Q.
Die letzten vier können relevant sein, wenn Ihre Koordinaten in den Ecken einer UTM-Zone liegen, aber höchst unwahrscheinlich sind.
Angesichts der Art und Weise, wie Sentinel2 Kacheln beschriftet, sollte nur einer der Nachbarn jemals einen solchen Ordner haben, um sicherzustellen, dass Sie die richtige Datei erhalten.
Jede andere, geografisch "korrektere" Lösung würde viel mehr Rechenaufwand bedeuten, als ich hier für gerechtfertigt halte.
Und definitiv, definitiv, melde dies dem ESA-Team, wie von Kersten in den Kommentaren vorgeschlagen. Ich verstehe wirklich nicht, warum sie sich für ein so unnötig kompliziertes Organisationssystem entschieden haben.
quelle
In Verbindung stehender Beitrag hier
Was für mich funktioniert hat, ist, die von der ESA bereitgestellte S2 KML zu verwenden, um alle Kacheln dort zu berechnen, die sich mit meiner AOI überschneiden, und dann in AWS nach diesen Kacheln zu suchen.
Diese KML scheint als Definition aller möglichen von S2 generierten Kachel-IDs zu funktionieren, wodurch viele überlappende Optionen eliminiert werden.
Wenn ich mir die KML ansehe (nur visuelle Inspektion, nicht 100% sicher), scheint es mir, dass Sie im schlimmsten Fall nach 4 Kacheln suchen müssten.
Es wäre schön, den Algorithmus zu haben, mit dem die ESA die KML definiert, um dies effizienter zu gestalten.
quelle