Importación de shp a Postgis usando Python y ogr

Question

Acabo de exportar una Tabla Postgis para shp usando este consejos pero yo no soy capaz de importar un shp a Postgis usar la misma biblioteca(ogr). Alguna idea? Muchas gracias f.

Answer 1

En Python puro, sin necesidad de utilizar el módulo subprocess (os.el sistema está en desuso) para llamar a ogr2ogr o shp2pgsql, por ejemplo):

1) con ogr

• Anexar un archivo shapefile a un postgis tabla con el GDAL/OGR interfaz de Python

2) con la ogr y psycopg2 del libro de Python Geoespacial de Desarrollo (Eric Westra), Capítulo 7, pág.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) con psycopg2 sólo

• Cómo guardar un nuevo shapefile en PostGIS?

4) con psycopg2 y otros espacial de las bibliotecas

• shp-postgis-esquema con pyshp (shapefile)
• shapefile_importer.pytambién con pyshp