Hi Kim, On Tue, Jun 29, 2010 at 6:37 PM, Kim Besson <[email protected]> wrote: > Hello Markus > Thanks for the reply > To do what you suggest, I would have to retrieve LOCATION's CRS EPSG code > and I'm not able to get it automatically. I only can get that PROJ4 string > by using g.proj. Is there a way to automatically get a LOCATIONS's EPSG > code?
That not possible (yet?) - also GDAL/OGR does not have a reverse lookup - too complicated or too many exceptions. The string +proj=utm +south +no_defs +zone=24 +a=6378137 +rf=298.257223563 +towgs84=0.000,0.000,0.000 +to_meter=1 looks like http://spatialreference.org/ref/epsg/32724/ to me. cs2cs -le | grep 'WGS 84' WGS84 a=6378137.0 rf=298.257223563 WGS 84 Markus _______________________________________________ grass-user mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/grass-user
