Я просто експортую таблицю Postgis до shp, використовуючи ці поради, але я не в змозі імпортувати shp до Postgis за допомогою тієї ж бібліотеки (ogr). Будь-яка ідея? Велике спасибі f.
Я просто експортую таблицю Postgis до shp, використовуючи ці поради, але я не в змозі імпортувати shp до Postgis за допомогою тієї ж бібліотеки (ogr). Будь-яка ідея? Велике спасибі f.
Відповіді:
У чистому Python, не використовуючи модуль підпроцесу (os.system застаріло) для виклику ogr2ogr
або shp2pgsql
, наприклад):
import os.path
import psycopg2
import osgeo.ogr
connection = psycopg2.connect("dbname=... user=...")
cursor = connection.cursor()
cursor.execute("DELETE FROM countries")
srcFile = os.path.join("DISTAL-data", "TM_WORLD_BORDERS-0.3","TM_WORLD_BORDERS-0.3.shp")
shapefile = osgeo.ogr.Open(srcFile)
layer = shapefile.GetLayer(0)
for i in range(layer.GetFeatureCount()):
feature = layer.GetFeature(i)
name = feature.GetField("NAME").decode("Latin-1")
wkt = feature.GetGeometryRef().ExportToWkt()
cursor.execute("INSERT INTO countries (name,outline) " +"VALUES (%s, ST_GeometryFromText(%s, " +"4326))", (name.encode("utf8"), wkt))
connection.commit()
ogr2ogr -f "PostgreSQL" PG:”dbname=DBNAME host=localhost" file.shp -nln TABLENAME