[GRASS-dev] [GRASS GIS] #3048: r.out.gdal does not write EPSG code of the projection

[GRASS GIS]

#3048: r.out.gdal does not write EPSG code of the projection
------------------------+-------------------------
 Reporter:  sbl         |      Owner:  grass-dev@…
     Type:  defect      |     Status:  new
 Priority:  normal      |  Milestone:  7.0.5
Component:  Raster      |    Version:  svn-trunk
 Keywords:  r.out.gdal  |        CPU:  Unspecified
 Platform:  All         |
------------------------+-------------------------
 GeoTiffs produced with GRASS GIS 6.4, 7.0.4 and 7.1.svn do not contain the
 AUTHORITY parameter for the projection, which is thus not properly
 recognized by other tools (e.g. QGIS, but esp. GeoServer).

 Related discussion on the GDAL mailinglist can be found here:
 http://thread.gmane.org/gmane.comp.gis.gdal.devel/43019

 g.proj -w gave the projection information as in the GeoTIFF, though
 without any AUTHORITY parameters (as it is supposed to according to the
 manual). The AUTHORITY tags seem to be added by r.out.gdal, except for the
 authority parameter for the entire projection. Even if g.proj -g returns
 epsg=25832.

 Compared to GDALs GeoTIFFs also AXIS parameters are missing (see below).

 {{{

 GRASS 7.0.4svn (ETRS_32N):~ >g.proj -g
 name=ETRS89 / UTM zone 32N
 datum=etrs89
 ellps=grs80
 proj=utm
 zone=32
 towgs84=0,0,0,0,0,0,0
 no_defs=defined
 epsg=25832
 unit=meter
 units=meters
 meters=1
 GRASS 7.0.4svn (ETRS_32N):~ >
 GRASS 7.0.4svn (ETRS_32N):~ >r.out.gdal input=test output=$HOME/test.tif
 GRASS 7.0.4svn (ETRS_32N):~ >gdalinfo $HOME/test.tif
 (...)
 Coordinate System is:
 PROJCS["UTM Zone 32, Northern Hemisphere",
     GEOGCS["grs80",
         DATUM["European_Terrestrial_Reference_System_1989",
 SPHEROID["Geodetic_Reference_System_1980",6378137,298.257222101,
                 AUTHORITY["EPSG","7019"]],
             TOWGS84[0,0,0,0,0,0,0],
             AUTHORITY["EPSG","6258"]],
         PRIMEM["Greenwich",0],
         UNIT["degree",0.0174532925199433]],
     PROJECTION["Transverse_Mercator"],
     PARAMETER["latitude_of_origin",0],
     PARAMETER["central_meridian",9],
     PARAMETER["scale_factor",0.9996],
     PARAMETER["false_easting",500000],
     PARAMETER["false_northing",0],
     UNIT["metre",1,
         AUTHORITY["EPSG","9001"]]]
 Origin = (0.000000000000000,1.000000000000000)
 Pixel Size = (1.000000000000000,-1.000000000000000)
 (...)
 GRASS 7.0.4svn (ETRS_32N):~ >
 GRASS 7.0.4svn (ETRS_32N):~ >gdal_translate -a_srs EPSG:25832
 $HOME/test.tif $HOME/test2.tif
 GRASS 7.0.4svn (ETRS_32N):~ >gdalinfo $HOME/test2.tif
 (...)
 Coordinate System is:
 PROJCS["ETRS89 / UTM zone 32N",
     GEOGCS["ETRS89",
         DATUM["European_Terrestrial_Reference_System_1989",
             SPHEROID["GRS 1980",6378137,298.257222101,
                 AUTHORITY["EPSG","7019"]],
             TOWGS84[0,0,0,0,0,0,0],
             AUTHORITY["EPSG","6258"]],
         PRIMEM["Greenwich",0,
             AUTHORITY["EPSG","8901"]],
         UNIT["degree",0.0174532925199433,
             AUTHORITY["EPSG","9122"]],
         AUTHORITY["EPSG","4258"]],
     PROJECTION["Transverse_Mercator"],
     PARAMETER["latitude_of_origin",0],
     PARAMETER["central_meridian",9],
     PARAMETER["scale_factor",0.9996],
     PARAMETER["false_easting",500000],
     PARAMETER["false_northing",0],
     UNIT["metre",1,
         AUTHORITY["EPSG","9001"]],
     AXIS["Easting",EAST],
     AXIS["Northing",NORTH],
     AUTHORITY["EPSG","25832"]]
 Origin = (0.000000000000000,1.000000000000000)
 Pixel Size = (1.000000000000000,-1.000000000000000)
 (...)
 }}}

--
Ticket URL: <https://trac.osgeo.org/grass/ticket/3048>
GRASS GIS <https://grass.osgeo.org>