Usando ogr2ogr para converter GML em shapefile em Python?

Estou tentando converter um GML em um shapefile ESRI usando o utilitário ogr2ogr em um script python.

Instalei com sucesso o pacote GDAL / OGR via osgeo, mas agora estou lutando para encontrar / entender quaisquer detalhes sobre a sintaxe para usar ogr2ogr em python.

Depois de importar ogr, tudo o que encontrei é ogr2ogr -f "ESRI Shapefile" output.shp input.gml

Como eu uso isso em python? É tão simples quanto atribuir os arquivos de entrada e saída?

Tudo o que eu tento recebe apenas erros de sintaxe. Qualquer ponteiro na direção certa para começar seria ótimo.

Pegue uma cópia da porta Python ogr2ogr, que é distribuída com o download do código-fonte GDAL ou pode ser encontrada aqui: http://svn.osgeo.org/gdal/trunk/gdal/swig/python/samples/ogr2ogr.py

Depois de importá-lo para o seu código, você pode usá-lo assim:

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"])

Você provavelmente encontrará respostas muito melhores do que posso fornecer em http://www.gis.usu.edu/~chrisg/python/2009/lectures/ospy_slides1.pdf e http://trac.osgeo.org/gdal/ wiki / GdalOgrInPython .

Tente lá ou outro tutorial vinculado a partir do segundo URL e talvez faça outra pergunta mais específica se você estiver preso.

Leia sobre os.system ou os.popen ou subprocesso . Ou apenas use ogr2ogr e seu shell. Até o cmd.exe (Windows) permite fazer um loop sobre os arquivos de entrada.

Uma explicação detalhada sobre como substituir a ogr2ogrferramenta de linha de comando pela API OGR Python pode ser encontrada aqui .