Ich bin sehr neu in GIS.
Ich benutze gdal
, um eine Landnutzungs- / Landbedeckungskarte einzulesen, und ich muss die Länge bestimmter Landbedeckungstypen auswählen, um sie in einem anderen Datensatz zu indizieren, der nur in Länge ausgedrückt wird. Unfortuantely, ich verstehe nicht , die Form der x und y , die ich gegeben Koordinaten aus dem Geotransformation, insbesondere die originX
und originY
unten:
geotransform = dataset.GetGeoTransform()
originX = geotransform[0]
originY = geotransform[3]
Das Drucken dieser Werte gibt mir Koordinaten wie (447466.693808, 4952570.40529)
. Wie hängen diese mit dem ursprünglichen Breiten- und Längengrad zusammen?
Bearbeiten:
Hier ist ein einfaches Python-Beispiel, mit dem ich gefunden habe, wonach ich gesucht habe:
srs = osr.SpatialReference()
srs.ImportFromWkt(dataset.GetProjection())
srsLatLong = srs.CloneGeogCS()
ct = osr.CoordinateTransformation(srs,srsLatLong)
print ct.TransformPoint(originX,originY)
Gestohlen von: tolatlong.py
dataset.GetProjectionRef()
und herausfinden kann, dass ich "UTM Zone 10" verwende, aber was dann? Ich google nach Methoden wie "unproject", komme aber auf null.