Even, I've opened #4880 on this - it is problematic in the R-spatial setting because the CRS (coordinate reference system) object recorded in each Spatial object uses the PROJ.4 string to represent spatial reference. A user can create a CRS, write out a raster (which expands the description to include , read it back, and the CRS are not the same. The same does not appear to happen with OGR write/read (ESRI Shapefile driver). Is the conclusion that calling programs must first find "NAD83" in the raster input WKT, if it is there set the environment variable if the default GDAL behaviour is not desired, export to Proj4, then remove the environment variable if set? Does this only affect NAD83 - it does seem to be limited to this datum, as WGS84 does not show the same behaviour.
-- View this message in context: http://osgeo-org.1560.n6.nabble.com/gdal-dev-Fwd-OSRExportToProj4-losing-datum-tp5013658p5013681.html Sent from the GDAL - Dev mailing list archive at Nabble.com. _______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev