Zaimportuj shp do Postgis przy użyciu Pythona i ogr

14

Właśnie eksportuję tabelę Postgis do shp, korzystając z tych wskazówek, ale nie jestem w stanie zaimportować shp do Postgis przy użyciu tej samej biblioteki (ogr). Dowolny pomysł? Wielkie dzięki

python postgis gdal ogr franco_b
źródło
1
Czy naprawdę potrzebujesz do tego Pythona, czy po prostu próbujesz zaimportować plik? Zgadzam się, że musisz zaimportować plik, a następnie użyć wier2ogr w wierszu polecenia ogr2ogr -f "PostgreSQL" PG:”dbname=DBNAME host=localhost" file.shp -nln TABLENAME
Jesse Crocker

Odpowiedzi:

29

W czystym Pythonie, bez użycia modułu podprocesu (os.system jest przestarzały) do wywołania ogr2ogrlub shp2pgsql, na przykład):

1) z ogrodem

2) z ogrodem i psycopg2 z książki Python Geospatial Development (Eric Westra), rozdział 7, s.219

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

3) tylko z psycopg2

4) z psycopg2 i innymi bibliotekami przestrzennymi

gen
źródło