tutoriel python gdal : géométrie et projection

Créer un objet géométrique vide : ogr.Geometry

Les méthodes utilisées pour définir les différentes géométries sont différentes (point, ligne, polygone, etc)

Créer un nouveau point, Utiliser la méthode AddPoint(, , []). La coordonnée z est généralement omise et la valeur par défaut est 0

Par exemple :

point = ogr.Geometry(ogr.wkbPoint)

point.AddPoint(10, 20)

Nouvelle ligne

Utilisez AddPoint(, , []) pour ajouter des points

Utilisez SetPoint(< index> , , , []) Changer les coordonnées du point

Par exemple, le code suivant modifie les coordonnées du point 0 :

line = ogr.Geometry(ogr.wkbLineString)

line.AddPoint(10,10)

line.AddPoint(20,20)

line.SetPoint(0 ,30, 30) #(10,10) -> (30,30)

Comptez le nombre de tous les points

imprimer la ligne.GetPointCount()

Lire point 0 La coordonnée x et la coordonnée y

print line.GetX(0)

print line.GetY(0)

Pour créer un nouveau polygone, créez d'abord un nouveau anneau, puis Ajouter des anneaux aux objets polygonaux.

Comment créer une bague ? Créez d’abord un nouvel objet anneau, puis ajoutez-y des points un par un.

ring = ogr.Geometry(ogr.wkbLinearRing)

ring.AddPoint(0,0)

ring.AddPoint(100,0)

ring.AddPoint(100,100)

ring.AddPoint(0,100)

À la fin, utilisez CloseRings pour fermer l'anneau, ou définissez les coordonnées du dernier point pour qu'elles soient les mêmes que celles du premier indiquer.

ring.CloseRings()

ring.AddPoint(0,0)

Ce qui suit est un exemple pour créer une boîte. Il s'agit d'un objet polygonal composé de deux couches d'anneaux.

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)

Les trois dernières phrases sont plus importantes, c'est-à-dire créez d'abord un objet polygone, puis ajoutez l'anneau extérieur et l'anneau intérieur

La phrase suivante peut vous aider à compter le nombre d'anneaux que votre polygone peut avoir

imprimer polygon.GetGeometryCount()

Lire les anneaux du polygone, l'ordre d'index et les ajouter lors de la création du polygone. l'ordre des anneaux est le même

outring = polygon.GetGeometryRef(0)

inring = polygon.GetGeometryRef(1)

Créer une géométrie composite multi-géométrie

Par exemple, MultiPoint, MultiLineString, MultiPolygon

Utilisez AddGeometry pour ajouter des formes géométriques ordinaires aux formes géométriques composites, par exemple :

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)

Lire la géométrie dans MultiGeometry La méthode est la même que la lecture de l'anneau de Polygon. a dit que Polygon Est une MultiGeometry intégrée.

Ne supprimez pas la géométrie d'une fonctionnalité existante, car cela ferait planter python

Vous ne pouvez supprimer que la géométrie créée lors de l'exécution du script, par exemple, créée manuellement ou automatiquement par appelant d'autres fonctions. Même si cette géométrie a été utilisée pour créer d'autres fonctionnalités, vous pouvez toujours la supprimer.

Par exemple : Polygon.Destroy()

Concernant les projections de projection, utilisez l'objet SpatialReference

Une variété de projections, GDAL prend en charge WKT, PROJ.4, ESPG, USGS, ESRI.prj

peut lire les projections de la couche et de la géométrie, par exemple :

spatialRef = layer.GetSpatialRef()

spatialRef = geom.GetSpatialReference()

Les informations de projection sont généralement stockées dans un fichier .prj. S'il n'existe pas de fichier de ce type, la fonction ci-dessus renvoie Aucun

Créez une nouvelle projection :

Importez d'abord la bibliothèque osr, puis utilisez osr.SpatialReference () Après avoir créé l'objet SpatialReference

, utilisez les instructions suivantes pour importer les informations de projection dans l'objet SpatialReference

•ImportFromWkt()

•ImportFromEPSG()

•ImportFromProj4()

•ImportFromESRI()

•ImportFromPCI( , ,

  )

 •ImportFromUSGS(, )

•ImportFromXML()

Exporter la projection, utiliser Les instructions suivantes peuvent être exportées sous forme de chaînes

•ExportToWkt()

•ExportToPrettyWkt()

•ExportToProj4()

•ExportToPCI()

•ExportToUSGS()

•ExportToXML()

Pour effectuer une transformation de projection sur une forme géométrique Geometry, vous devez d'abord initialiser deux Projections, puis créer un objet CooperativeTransformation et l'utiliser pour la transformation

sourceSR = osr.SpatialReference()

sourceSR .ImportFromEPSG(32612) #UTM 12N WGS84

targetSR = osr.SpatialReference()

targetSR.ImportFromEPSG(4326) #Geo WGS84

coordTrans = osr.CoordonnéeTransformation(sourceSR , targetSR)

geom.Transform(coordTrans)

Mais ce code est très pénible ! Cela ne fonctionne pas sous Windows. Il y a une discussion sur le forum des étrangers, disant qu'il n'y a pas de problème sous Linux, mais que Windows est mort, hé. . .

Il y a quelques autres points auxquels il faut faire attention :

Modifiez la géométrie au moment approprié. Il est préférable de ne pas la déplacer après la transformation de la projection.

Pour effectuer une transformation de projection sur toute la géométrie d'une DataSource, vous devez le faire une par une. Utilisez une boucle.

Écrire votre projection dans un fichier .prj est en fait très simple. Tout d'abord, MorphToESRI(), convertissez-le en chaîne, puis ouvrez un fichier texte et écrivez-le dedans. Par exemple :

targetSR.MorphToESRI()

file = open('test.prj', 'w')

file.write(targetSR.ExportToWkt())

ffile.close()

Ce qui précède est le contenu du didacticiel python gdal : géométrie et projection. Pour plus de contenu connexe, veuillez faire attention au site Web PHP chinois (www.php. cn) !