Convertendo a tabela PostGIS para Shapefile em Python?

10

Eu quero converter uma tabela PostGIS em um shapefile (sem usar o pgsql2shp).

Para criar uma geometria no arquivo shapefile, preciso fornecer o Xmin, Ymin e Xmax, Ymax, e a geometria que tenho na minha tabela PostGIS é de formato irregular (posso obter o exterior usando a caixa delimitadora, mas isso incluirá alguma área extra mais do que a minha área de ineterest). Existe algum método pelo qual eu possa executar a tarefa?

Eu quero fazer a coisa de forma programática e usando Python.

python postgis Vicky
fonte

Respostas:

12

Se você quiser fazer isso programaticamente com Python, provavelmente o GDAL / OGR é a melhor maneira de fazê-lo. Veja o exemplo de código de exemplo que copia a tabela do PostgreSQL para o arquivo SHP. O exemplo não é perfeito, mas você pode modificá-lo facilmente para atender às suas necessidades.

import os
os.environ['PATH'] = "c:\\Program Files\\GDAL\\bin" + ';' + os.environ['PATH']
os.environ['GDAL_DRIVER_PATH'] = "c:\\Program Files\\GDAL\\bin\\gdal\\plugins-optional"
os.environ['GDAL_DATA'] = "c:\\Program Files\\GDAL\\bin\\gdal-data"
import ogr

conn=ogr.Open("PG: host=192.168.5.3 dbname=some_database user=postgres password=xxxx")
if conn is None:
    print 'Could not open a database or GDAL is not correctly installed!'
    sys.exit(1)

output = "d:\\points.shp"

# Schema definition of SHP file
out_driver = ogr.GetDriverByName( 'ESRI Shapefile' )
out_ds = out_driver.CreateDataSource(output)
out_srs = None
out_layer = out_ds.CreateLayer("point", out_srs, ogr.wkbPoint)
fd = ogr.FieldDefn('name',ogr.OFTString)
out_layer.CreateField(fd)


layer = conn.GetLayerByName("point_data")
#layer = conn.ExecuteSQL(sql)

feat = layer.GetNextFeature()
while feat is not None:
    featDef = ogr.Feature(out_layer.GetLayerDefn())
    featDef.SetGeometry(feat.GetGeometryRef())
    featDef.SetField('name',feat.TITLE)
    out_layer.CreateFeature(featDef)
    feat.Destroy()
    feat = layer.GetNextFeature()

conn.Destroy()
out_ds.Destroy()

Se o seu ambiente de programação for o Windows, você poderá fazer o download do GDAL / OGR aqui . Alguns bons materiais de partida você pode encontrar aqui . Espero que ajude.

Mario Miler
fonte
1

Eu posso aconselhá-lo a dar uma olhada na biblioteca GDAL / OGR: http://www.gdal.org/ogr/index.html

Alex Markov
fonte