Get srid from a shapefile using ogr

Question

I'm writing a script to import shp to postgis, using python, osgeo and psycopg2 ( see my previous question: Import shp to Postgis using Python and ogr )

I'd know if there is a way to get by osgeo library the SRID (for example 4326) of the shp.

So I can confirm the SRID or transform it using script directly.

Answer 1

Go to the Python GDAL/OGR Cookbook 1.0 documentation and you'll have the answers to all your questions:

    from osgeo import ogr
    driver = ogr.GetDriverByName('ESRI Shapefile')
    shape = driver.Open('my.shp')
    layer= shape.GetLayer()
    # the crs
    crs = layer.GetSpatialRef()

and you can also create a projection file if the shapefile does not have one or reproject the shapefile with osgeo.osr

Answer 2

ogrinfo layer.shp layer -so will give you general information about the shapefile.

Answer 3

for display crs finally you can use crs.proj4

from osgeo import ogr
driver = ogr.GetDriverByName('ESRI Shapefile')
shape = driver.Open('my.shp')
layer= shape.GetLayer()
# the crs

crs = layer.GetSpatialRef()

crs.proj4()

Answer 4

If you don't want to use Python, gdalsrsinfo (which you probalby have installed if you have ogrinfo) will also work to display the SRS from a shapefile:

gdalsrsinfo shapefile.shp

Answer 5

You could user ogrinfo

ogrinfo -so -al shapefile.shp

-so : Summary Only: suppress listing of individual features and show only summary information like projection, schema, feature count and extents.

-al : List all features of all layers (used instead of having to give layer names as arguments).