PythonでPostGISテーブルをShapefileに変換しますか？

PostGISテーブルをシェープファイルに変換したい（pgsql2shpを使用しない）。

シェープファイルでジオメトリを作成するには、Xmin、YminとXmax、Ymaxを指定する必要があり、PostGISテーブルにあるジオメトリは不規則な形状です（境界ボックスを使用して外部を取得できますが、これには以下が含まれます）不穏な私の領域よりもいくつかの余分な領域）。タスクを完了する方法はありますか？

プログラムで、Pythonを使用して実行したいと考えています。

プログラムでPythonを使用して実行する場合は、おそらくGDAL / OGRが最適です。PostgreSQLからSHPファイルにテーブルをコピーするサンプルコードの例を見てください。例は完璧ではありませんが、ニーズに合わせて簡単に変更できます。

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()

プログラミング環境がWindowsの場合は、ここからGDAL / OGRをダウンロードできます。あなたがここで見つけることができるいくつかの良い出発材料。それが役に立てば幸い。

GDAL / OGRライブラリをご覧になることをお勧めします：http : //www.gdal.org/ogr/index.html