Перетворення таблиці PostGIS в Shapefile в Python?


10

Я хочу перетворити таблицю PostGIS у файл форми (не використовуючи pgsql2shp).

Для того, щоб створити геометрію у формі форми, я повинен надати Xmin, Ymin та Xmax, Ymax, а геометрія, яку я маю у своїй таблиці PostGIS, має неправильну форму (я можу отримати зовнішній вигляд за допомогою обмежувального вікна, але воно буде включати якусь зайву площу більше, ніж моя область неінтересних). Чи є якийсь метод, за допомогою якого я можу виконати завдання?

Я хочу зробити цю справу програмно і за допомогою Python.

Відповіді:


12

Якщо ви хочете це робити програмно з 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 звідси . Деякі хороші вихідні матеріали ви можете знайти тут . Сподіваюся, це допомагає.


1

Можу порадити поглянути на бібліотеку GDAL / OGR: http://www.gdal.org/ogr/index.html

Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.