[GRASS-dev] [GRASS GIS] #3193: r.out.gdal: add AUTHORITY node to the srs info of the exported raster

[GRASS GIS]

#3193: r.out.gdal: add AUTHORITY node to the srs info of the exported raster
--------------------------+---------------------------------------
  Reporter:  hellik       |      Owner:  grass-dev@…
      Type:  enhancement  |     Status:  new
  Priority:  major        |  Milestone:  7.4.0
 Component:  Raster       |    Version:  svn-trunk
Resolution:               |   Keywords:  r.out.gdal,r.external.out
       CPU:  All          |   Platform:  All
--------------------------+---------------------------------------

Comment (by sbl):

 Hmm... This is still an issue for me both in GRASS 7.4 (r72108) and GRASS
 7.5 (r72157):\\
 See also (sort of duplicate): #3048

 {{{
 g.proj -p
 -PROJ_INFO-------------------------------------------------
 name       : ETRS89 / UTM zone 33N
 datum      : etrs89
 ellps      : grs80
 proj       : utm
 zone       : 33
 no_defs    : defined
 -PROJ_EPSG-------------------------------------------------
 epsg       : 25833
 -PROJ_UNITS------------------------------------------------
 unit       : meter
 units      : meters
 meters     : 1

 r.external.out --verbose directory=/tmp extension=.tif format=GTiff
 options=COMPRESS=LZW
 r.mapcalc expression="test=1"
 gdalinfo /tmp/test.tif

 r.external.out -r
 g.remove type=raster name=test -f
 rm /tmp/test.tif
 r.mapcalc expression="test=1"
 r.out.gdal -c input=test output=/tmp/test.tif createopt="COMPRESS=LZW" --o
 gdalinfo /tmp/test.tif
 }}}

 In both cases gdalinfo shows:

 {{{
 Coordinate System is:
 PROJCS["UTM Zone 33, Northern Hemisphere",
     GEOGCS["grs80",
         DATUM["European_Terrestrial_Reference_System_1989",
 SPHEROID["Geodetic_Reference_System_1980",6378137,298.257222101,
                 AUTHORITY["EPSG","7019"]],
             AUTHORITY["EPSG","6258"]],
         PRIMEM["Greenwich",0],
         UNIT["degree",0.0174532925199433]],
     PROJECTION["Transverse_Mercator"],
     PARAMETER["latitude_of_origin",0],
     PARAMETER["central_meridian",15],
     PARAMETER["scale_factor",0.9996],
     PARAMETER["false_easting",500000],
     PARAMETER["false_northing",0],
     UNIT["metre",1,
         AUTHORITY["EPSG","9001"]]]
 Origin = (256495.000000000000000,6653695.000000000000000)
 Pixel Size = (10.000000000000000,-10.000000000000000)
 }}}

 That means:

 {{{
     AXIS["Easting",EAST],
     AXIS["Northing",NORTH],
     AUTHORITY["EPSG","25833"]]
 }}}

 is still missing at the bottom of the CRS information, compared to CRS
 info with:

 {{{
 gdal_edit.py -a_srs EPSG:25833 /tmp/test.tif
 gdalinfo /tmp/test.tif
 }}}

--
Ticket URL: <https://trac.osgeo.org/grass/ticket/3193#comment:7>
GRASS GIS <https://grass.osgeo.org>