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 <[email protected]> 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 > [email protected] > http://lists.osgeo.org/mailman/listinfo/gdal-dev _______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
