Heim > Artikel > Backend-Entwicklung > Python-GDAL-Tutorial: Geometrie und Projektion
Erstellen Sie ein leeres Geometrieobjekt: ogr.Geometry
Die zum Definieren verschiedener Geometrien verwendeten Methoden sind unterschiedlich (Punkt, Linie, Polygon usw.)
Erstellen Sie einen neuen Punkt. Verwenden Sie die Methode AddPoint(
Zum Beispiel:
point = ogr.Geometry(ogr.wkbPoint)
point.AddPoint(10, 20 )
Neue Zeile
Verwenden Sie AddPoint(
Verwenden Sie SetPoint(< index> ,
Der folgende Code ändert beispielsweise die Koordinaten von Punkt 0:
line = ogr.Geometry(ogr.wkbLineString)
line.AddPoint(10,10)
line.AddPoint(20,20)
line.SetPoint(0 ,30, 30) #(10,10) -> (30,30)
Zähle die Anzahl aller Punkte
print line.GetPointCount()
Lesen Punkt 0 Die x-Koordinate und die y-Koordinate
print line.GetX(0)
print line.GetY(0)
Um ein neues Polygon zu erstellen, erstellen Sie zunächst ein neues Klicken Sie auf „Ring“ und dann auf „Ringe zum Polygonobjekt hinzufügen“.
Wie erstelle ich einen Ring? Erstellen Sie zunächst ein neues Ringobjekt und fügen Sie ihm dann nacheinander Punkte hinzu.
ring = ogr.Geometry(ogr.wkbLinearRing)
ring.AddPoint(0,0)
ring.AddPoint(100,0)
ring.AddPoint(100,100)
ring.AddPoint(0,100)
Verwenden Sie am Ende CloseRings, um den Ring zu schließen, oder legen Sie die Koordinaten des letzten Punkts so fest, dass sie mit denen des ersten übereinstimmen Punkt.
ring.CloseRings()
ring.AddPoint(0,0)
Das Folgende ist ein Beispiel zum Erstellen einer Box. Dies ist ein Polygonobjekt, das aus zwei Ringschichten besteht.
outring = ogr.Geometry(ogr.wkbLinearRing)
outring.AddPoint(0,0)
outring.AddPoint(100,0)
outring.AddPoint(100,100)
outring.AddPoint(0,100)
outring.AddPoint(0,0)
inring = ogr.Geometry(ogr.wkbLinearRing)inring = ogr.Geometry(ogr.wkbLinearRing)
inring.AddPoint(25,25)
inring.AddPoint(75,25)
inring.AddPoint(75,75)
inring.AddPoint(25,75)
inring.CloseRings()
polygon = ogr.Geometry(ogr.wkbPolygon)
polygon.AddGeometry (Outring)
polygon.AddGeometry(inring)
Die letzten drei Sätze sind wichtiger, dh erstellen Sie zuerst ein Polygonobjekt und fügen Sie dann den Außenring und den Innenring hinzu
Der folgende Satz kann Ihnen dabei helfen, zu zählen, wie viele Ringe Ihr Polygon haben kann
polygon drucken.GetGeometryCount()
Ringe aus dem Polygon lesen, indizieren und beim Erstellen des Polygons hinzufügen Reihenfolge der Ringe ist die gleiche
outring = polygon.GetGeometryRef(0)
inring = polygon.GetGeometryRef(1)
Erstellen Sie eine zusammengesetzte Geometrie mit mehreren Geometrien
Zum Beispiel MultiPoint, MultiLineString, MultiPolygon
Verwenden Sie AddGeometry, um gewöhnliche geometrische Formen zu zusammengesetzten geometrischen Formen hinzuzufügen, zum Beispiel:
multipoint = ogr.Geometry(ogr.wkbMultiPoint)
point = ogr.Geometry(ogr.wkbPoint)point = ogr.Geometry(ogr.wkbPoint)
point.AddPoint(10,10)
multipoint.AddGeometry(point)
point.AddPoint(20,20)
multipoint.AddGeometry(point)
Geometrie lesen in MultiGeometry sagte, dass Polygon eine integrierte MultiGeometrie ist.
Löschen Sie nicht die Geometrie eines vorhandenen Features, da dies zum Absturz von Python führt.
Sie können nur die Geometrie löschen, die während der Ausführung des Skripts erstellt wurde, z. B. manuell erstellt oder automatisch erstellt Aufruf anderer Funktionen. Auch wenn diese Geometrie zum Erstellen anderer Features verwendet wurde, können Sie sie dennoch löschen.
Zum Beispiel: Polygon.Destroy()
Für Projektionsprojektionen verwenden Sie das SpatialReference-Objekt
Eine Vielzahl von Projektionen, GDAL unterstützt WKT, PROJ.4, ESPG, USGS , ESRI.prj
kann Projektionen aus Layer und Geometrie lesen, zum Beispiel:
spatialRef = layer.GetSpatialRef()
spatialRef = geom.GetSpatialReference()
Projektionsinformationen werden im Allgemeinen in einer .prj-Datei gespeichert, wenn keine solche Datei vorhanden ist.
Erstellen Sie eine neue Projektion:
Importieren Sie zunächst die OSR-Bibliothek und Verwenden Sie dann osr.SpatialReference (). Nachdem Sie das SpatialReference-Objekt
erstellt haben, verwenden Sie die folgenden Anweisungen, um die Projektionsinformationen in das SpatialReference-Objekt
zu importieren: •ImportFromWkt(
•ImportFromEPSG(
•ImportFromProj4(
•ImportFromESRI(
•ImportFromPCI(
•ImportFromUSGS(
•ImportFromXML(
Projektion exportieren, verwenden Die folgenden Anweisungen können als Zeichenfolgen exportiert werden
•ExportToWkt()
•ExportToPrettyWkt()
•ExportToProj4()
•ExportToPCI()
•ExportToUSGS()
•ExportToXML()
Um eine Projektionstransformation für eine geometrische Formgeometrie durchzuführen, müssen Sie zuerst zwei Projektionen initialisieren, dann ein CoordinateTransformation-Objekt erstellen und es für die Transformation verwenden
sourceSR = osr.SpatialReference()
sourceSR .ImportFromEPSG(32612) #UTM 12N WGS84
targetSR = osr.SpatialReference()
targetSR.ImportFromEPSG(4326) #Geo WGS84
coordTrans = osr.CoordinateTransformation(sourceSR , targetSR)
geom.Transform(coordTrans)
Aber dieser Code ist sehr schmerzhaft! Unter Windows funktioniert es nicht. Im Ausländerforum gibt es eine Diskussion darüber, dass es unter Linux kein Problem gibt, aber Windows ist tot, hey. . .
Es gibt noch ein paar andere Dinge, auf die Sie achten sollten:
Bearbeiten Sie die Geometrie zum richtigen Zeitpunkt. Es ist am besten, sie nach der Projektionstransformation nicht zu verschieben.
Um eine Projektionstransformation für alle Geometrien in einer Datenquelle durchzuführen, müssen Sie diese einzeln durchführen. Verwenden Sie eine Schleife.
Das Schreiben Ihrer Projektion in eine .prj-Datei ist eigentlich sehr einfach. Zuerst MorphToESRI(), wandeln Sie es in einen String um, öffnen Sie dann eine Textdatei und schreiben Sie es hinein. Zum Beispiel:
targetSR.MorphToESRI()
file = open('test.prj', 'w')
file.write(targetSR.ExportToWkt())
ffile.close()
Das Obige ist der Inhalt des Python-Gdal-Tutorials: Geometrie und Projektion. Weitere verwandte Inhalte finden Sie auf der chinesischen PHP-Website (www.php. cn)!