Utiliser ogr2ogr pour convertir du GML en shapefile en Python?

J'essaie de convertir un GML en un fichier de formes ESRI en utilisant l'utilitaire ogr2ogr dans un script python.

J'ai installé avec succès le paquet GDAL / OGR via osgeo mais j'ai du mal à trouver / comprendre tous les détails sur la syntaxe pour utiliser ogr2ogr en python.

Après avoir importé ogr, tout ce que j'ai trouvé est ogr2ogr -f "ESRI Shapefile" output.shp input.gml

Comment utiliser cela en python? Est-ce aussi simple que d'attribuer les fichiers d'entrée et de sortie?

Tout ce que j'essaie, je reçois juste des erreurs de syntaxe. Tout pointeur dans la bonne direction pour me lancer serait formidable.

Récupérez une copie du port Python ogr2ogr, qui est distribué avec le téléchargement du code source GDAL ou peut être trouvé ici: http://svn.osgeo.org/gdal/trunk/gdal/swig/python/samples/ogr2ogr.py

Une fois que vous avez importé cela dans votre code, vous pouvez l'utiliser comme ceci:

import ogr2ogr

def main():
  #note: main is expecting sys.argv, where the first argument is the script name
  #so, the argument indices in the array need to be offset by 1
  ogr2ogr.main(["","-f", "KML", "out.kml", "data/san_andres_y_providencia_administrative.shp"])

Vous trouverez probablement de bien meilleures réponses que celles que je peux vous donner sur http://www.gis.usu.edu/~chrisg/python/2009/lectures/ospy_slides1.pdf et http://trac.osgeo.org/gdal/ wiki / GdalOgrInPython .

Essayez-le ou un autre didacticiel lié à partir de la deuxième URL et posez peut-être une autre question plus spécifique si vous êtes bloqué.

Lisez sur os.system ou os.popen ou sous - processus . Ou utilisez simplement ogr2ogr et votre shell. Même cmd.exe (Windows) vous permet de parcourir les fichiers d'entrée.

Une explication détaillée sur la façon de remplacer l' ogr2ogroutil de ligne de commande par l'API OGR Python peut être trouvée ici .