NATO UTM-Standort in Sentinel-2 kann nicht gefunden werden

10

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 FRUnterverzeichnis finden .

In Google ist dieser Ort hier:

Geben Sie hier die Bildbeschreibung ein

Und wenn ich den gleichen Ort im Sentinel-Bildbrowser finde, den ich sehe, ist diese Kachel anders

Geben Sie hier die Bildbeschreibung ein

was für13/S/FR dh gleiche UTMund quadratische, aber unterschiedliche Band steht.

Wie ist das möglich?

AKTUALISIEREN

KML mit Sentinel-2-Kacheln meldet auch SKacheln an einem bestimmten Ort

Geben Sie hier die Bildbeschreibung ein

UPDATE 2

Nach diesem Bild

Geben Sie hier die Bildbeschreibung ein

genommen von hier , das FRist Platz befindet sich die Hälfte in Sin UTM - Zone und eine halbe RZone. Offensichtlich weisen die meisten automatischen Konverter dieses Quadrat der RZone zu, während Sentinel-2 es für die SZone 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?

Dims
quelle
Es ist S/FR, während UTM-Konverter geben R/FR. Wie berechnet man den Standort, wenn UTM-Konverter nicht richtig funktionieren?
Dims
Der Breitengrad liegt bei knapp 32 Grad Nord. Damit liegt es genau im R-Breitengrad "Band". Sentinel-2 hat möglicherweise gekachelt, indem der Mittelpunkt der Kachel verwendet wurde, der sich stattdessen im "S" -Band befinden könnte.
Mkennedy
@mkennedy Wie simuliere ich diesen Algorithmus ausgehend von Koordinaten?
Dims
2
Sie können dies auch an [email protected] melden, da dies tatsächlich nach unerwartetem Verhalten aussieht.
Kersten

Antworten:

1

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:

  1. Verwenden Sie einen der UTM-Konverter, die "wie erwartet" funktionieren, und platzieren Sie die Koordinaten in 13R.
  2. Überprüfen Sie anschließend, ob der Ordner in der Sentinel 2-Datenstruktur vorhanden ist. Wenn ja, sind Sie fertig, Hurra.

  3. 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.

Senshi
quelle
0

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.

Fred Liporace
quelle