Heim  >  Artikel  >  Backend-Entwicklung  >  Python-GDAL-Tutorial: Geometrie und Projektion

Python-GDAL-Tutorial: Geometrie und Projektion

黄舟
黄舟Original
2016-12-24 16:01:174800Durchsuche

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(, , []). Die Z-Koordinate wird im Allgemeinen weggelassen und der Standardwert ist 0

Zum Beispiel:

point = ogr.Geometry(ogr.wkbPoint)

point.AddPoint(10, 20 )

Neue Zeile

Verwenden Sie AddPoint(, , []), um Punkte hinzuzufügen

Verwenden Sie SetPoint(< index> , , , []) Ändern Sie die Koordinaten des Punktes

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)!


Stellungnahme:
Der Inhalt dieses Artikels wird freiwillig von Internetnutzern beigesteuert und das Urheberrecht liegt beim ursprünglichen Autor. Diese Website übernimmt keine entsprechende rechtliche Verantwortung. Wenn Sie Inhalte finden, bei denen der Verdacht eines Plagiats oder einer Rechtsverletzung besteht, wenden Sie sich bitte an admin@php.cn