[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.1
 Component:  Raster       |    Version:  svn-trunk
Resolution:               |   Keywords:  r.out.gdal,r.external.out
       CPU:  All          |   Platform:  All
--------------------------+---------------------------------------

Comment (by hellik):

 Replying to [comment:10 mmetz]:
 > Replying to [comment:9 hellik]:
 > > Replying to [comment:6 mmetz]:
 > > > Replying to [comment:5 martinl]:
 > > > > This functionality was recently implemented in trunk, right?
 Reference is missing here.
 > > >
 > > > It was implemented by martinl with r71303
 > >
 > > tested with the example above:
 > > [...]
 > >
 > > {{{
 > > AUTHORITY["EPSG","31255"]
 > > }}}
 > >
 > > seems to be still missing
 >
 > The functions were implemented in r71303. An EPSG code was only exported
 with v.out.ogr. As of trunk r72160, an existing EPSG code is exported with
 r.external.out, r.out.gdal, v.external.out, v.out.ogr.

 tested with

 {{{
 System Info
 GRASS version: 7.5.svn
 GRASS SVN revision: r72194
 Build date: 2018-01-30
 Build platform: x86_64-w64-mingw32
 GDAL: 2.2.3
 PROJ.4: 4.9.3
 GEOS: 3.5.0
 SQLite: 3.17.0
 Python: 2.7.5
 wxPython: 2.8.12.1
 Platform: Windows-8-6.2.9200 (OSGeo4W)
 }}}

 {{{
 g.proj -p
 -PROJ_INFO-------------------------------------------------
 name       : MGI / Austria GK Central
 datum      : hermannskogel
 ellps      : bessel
 proj       : tmerc
 lat_0      : 0
 lon_0      : 13.33333333333333
 k          : 1
 x_0        : 0
 y_0        : -5000000
 no_defs    : defined
 towgs84    : 577.326,90.129,463.919,5.1366,1.4742,5.2970,2.4232
 -PROJ_EPSG-------------------------------------------------
 epsg       : 31255
 -PROJ_UNITS------------------------------------------------
 unit       : meter
 units      : meters
 meters     : 1
 }}}


 {{{
 gdalinfo test_raster_authority.tif
 Driver: GTiff/GeoTIFF
 Files: test_raster_authority.tif
 Size is 170, 141
 Coordinate System is:
 PROJCS["MGI / Austria GK Central",
     GEOGCS["MGI",
         DATUM["Militar_Geographische_Institute",
             SPHEROID["Bessel 1841",6377397.155,299.1528128,
                 AUTHORITY["EPSG","7004"]],
             TOWGS84[577.326,90.129,463.919,5.1366,1.4742,5.297,2.4232],
             AUTHORITY["EPSG","6312"]],
         PRIMEM["Greenwich",0,
             AUTHORITY["EPSG","8901"]],
         UNIT["degree",0.0174532925199433,
             AUTHORITY["EPSG","9122"]],
         AUTHORITY["EPSG","4312"]],
     PROJECTION["Transverse_Mercator"],
     PARAMETER["latitude_of_origin",0],
     PARAMETER["central_meridian",13.33333333333333],
     PARAMETER["scale_factor",1],
     PARAMETER["false_easting",0],
     PARAMETER["false_northing",-5000000],
     UNIT["metre",1,
         AUTHORITY["EPSG","9001"]],
     AUTHORITY["EPSG","31255"]] <=
 Origin = (-161500.000000000000000,273250.000000000000000)
 Pixel Size = (250.000000000000000,-250.000000000000000)
 Metadata:
   AREA_OR_POINT=Area
   TIFFTAG_SOFTWARE=GRASS GIS 7.5.svn with GDAL 2.2.3
 Image Structure Metadata:
   COMPRESSION=LZW
   INTERLEAVE=BAND
 Corner Coordinates:
 Upper Left  ( -161500.000,  273250.000) ( 11d11'10.74"E, 47d34'42.28"N)
 Lower Left  ( -161500.000,  238000.000) ( 11d11'56.99"E, 47d15'41.55"N)
 Upper Right ( -119000.000,  273250.000) ( 11d45' 3.81"E, 47d35'15.35"N)
 Lower Right ( -119000.000,  238000.000) ( 11d45'37.92"E, 47d16'14.25"N)
 Center      ( -140250.000,  255625.000) ( 11d28'27.33"E, 47d25'29.61"N)
 }}}

  AUTHORITY["EPSG","XXXX"] is now exported in rasters too.

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