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.
Respostas:
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:
fonte
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.
fonte
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.
fonte
Uma explicação detalhada sobre como substituir a
ogr2ogr
ferramenta de linha de comando pela API OGR Python pode ser encontrada aqui .fonte