[gdal-dev] Getting SRID from ESRI shapefile

[Etienne Tourigny]

Alessandro,

unfortunately there is not a reliable tool on GDAL at the moment, but
I have been working on experimental EPSG lookup, see ticket #4345

http://trac.osgeo.org/gdal/ticket/4345

It requires gdalsrsinfo from gdal-1.9 and some extra data files
installed in the gdal data directory (see the ticket)

usage:
$gdalsrsinfo -e 102067.prj

However, it does not work with your example because it is not an EPSG
definition, but an ESRI one.

Since SRID 102067 is in esri_extra.wkt,  is recognized by gdal (and
spatialreferenc.org), but this has not been added in my experimental
EPSG lookup.

$gdalsrsinfo EPSG:102067
PROJ.4 : '+proj=krovak +lat_0=49.5 +lon_0=24.83333333333333
+alpha=30.28813975277778 +k=0.9999 +x_0=0 +y_0=0 +ellps=bessel
+units=m +no_defs '

OGC WKT :
PROJCS["S-JTSK_Krovak_East_North",
    GEOGCS["GCS_S_JTSK",
        DATUM["Jednotne_Trigonometricke_Site_Katastralni",
            SPHEROID["Bessel_1841",6377397.155,299.1528128]],
        PRIMEM["Greenwich",0],
        UNIT["Degree",0.017453292519943295]],
    PROJECTION["Krovak"],
    PARAMETER["False_Easting",0],
    PARAMETER["False_Northing",0],
    PARAMETER["Pseudo_Standard_Parallel_1",78.5],
    PARAMETER["Scale_Factor",0.9999],
    PARAMETER["Azimuth",30.28813975277778],
    PARAMETER["Longitude_Of_Center",24.83333333333333],
    PARAMETER["Latitude_Of_Center",49.5],
    PARAMETER["X_Scale",-1],
    PARAMETER["Y_Scale",1],
    PARAMETER["XY_Plane_Rotation",90],
    UNIT["Meter",1],
    AUTHORITY["EPSG","102067"]]


regards, Etienne


On Mon, Feb 6, 2012 at 11:11 AM, Alessandro Candini <candini at meeo.it> wrote:
> Is it possible to get a SRID code from an ESRI shapefile?
> For example, I've got the following:
>
> ~ $ ogrinfo -al LC_FIN_smooth.shp | head -n 30
> INFO: Open of `LC_FIN_smooth.shp'
>      using driver `ESRI Shapefile' successful.
>
> Layer name: LC_FIN_smooth
> Geometry: Polygon
> Feature Count: 33308
> Extent: (-598426.000000, -1172874.000000) - (-583748.000000,
> -1158198.000000)
> Layer SRS WKT:
> PROJCS["S-JTSK_Krovak_East_North",
>    GEOGCS["GCS_S_JTSK",
>        DATUM["Jednotne_Trigonometricke_Site_Katastralni",
>            SPHEROID["Bessel_1841",6377397.155,299.1528128]],
>        PRIMEM["Greenwich",0.0],
>        UNIT["Degree",0.0174532925199433]],
>    PROJECTION["Krovak"],
>    PARAMETER["False_Easting",0.0],
>    PARAMETER["False_Northing",0.0],
>    PARAMETER["Pseudo_Standard_Parallel_1",78.5],
>    PARAMETER["Scale_Factor",0.9999],
>    PARAMETER["Azimuth",30.28813975277778],
>    PARAMETER["Longitude_Of_Center",24.83333333333333],
>    PARAMETER["Latitude_Of_Center",49.5],
>    PARAMETER["X_Scale",-1.0],
>    PARAMETER["Y_Scale",1.0],
>    PARAMETER["XY_Plane_Rotation",90.0],
>    UNIT["Meter",1.0]]
> Class_name: String (254.0)
> LC: Real (32.16)
> OGRFeature(LC_FIN_smooth):0
>  Class_name (String) = out_of_AOI
>
> At spatialreference.org I found that this correspond to SRID 102067.
> But is there a way to obtain this number with gdal tools?
>
> Thanks in advance.
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/gdal-dev