[GRASS-dev] [GRASS GIS] #2456: read CSV from GDAL data directory

[GRASS GIS]

#2456: read CSV from GDAL data directory
--------------------------------+-------------------------------------------
 Reporter:  martinl             |       Owner:  grass-dev@…              
     Type:  task                |      Status:  new                      
 Priority:  critical            |   Milestone:  7.1.0                    
Component:  Projections/Datums  |     Version:  unspecified              
 Keywords:  gdal                |    Platform:  Unspecified              
      Cpu:  Unspecified         |  
--------------------------------+-------------------------------------------

Comment(by martinl):

 Replying to [comment:3 martinl]:
 > For record list of affected files:
 >  * pcs.csv

 The attached patch attachment:gdal_data.diff forces GRASS to use CSV files
 from GDAL. First found problem:

 GRASS defines 5514 as

 {{{
 pcs.csv:5514,"S-JTSK / Krovak East
 North",9001,4818,19952,9819,1,0,6501,8806,0,9001,8807,0,9001,8811,49.3,9110,8813,30.1717303,9110,8818,78.3,9110,8819,0.9999,9201,8833,42.3,9110
 }}}

 ->

 {{{
 g.proj epsg=5514 -p
 -PROJ_INFO-------------------------------------------------
 name       : Krovak
 proj       : krovak
 datum      : hermannskogel
 ellps      : bessel
 lat_0      : 49.5
 lon_0      : 42.5
 alpha      : 30.28813972222222
 k          : 0.9999
 x_0        : 0
 y_0        : 0
 pm         : ferro
 no_defs    : defined
 -PROJ_UNITS------------------------------------------------
 unit       : meter
 units      : meters
 meters     : 1
 }}}

 GDAL as

 {{{
 pcs.csv:5514,"S-JTSK / Krovak East
 North",9001,4156,5510,1041,1,0,4499,1036,30.1717303,9110,8806,0,9001,8807,0,9001,8811,49.3,9110,8818,78.3,9110,8819,0.9999,9201,8833,24.5,9110
 }}}

 ->

 {{{
 .proj epsg=5514 -p
 -PROJ_INFO-------------------------------------------------
 name       : Krovak
 proj       : krovak
 ellps      : bessel
 lat_0      : 49.5
 lon_0      : 24.83333333333333
 alpha      : 30.28813972222222
 k          : 0.9999
 x_0        : 0
 y_0        : 0
 towgs84    : 589,76,480,0,0,0,0
 no_defs    : defined
 -PROJ_UNITS------------------------------------------------
 unit       : meter
 units      : meters
 meters     : 1
 }}}

 The `wgs84` causes that

 {{{
 g.proj epsg=5514 datum_trans=-1
 ERROR: No output format specified, define one of flags -p, -g, -j, or -w
 }}}

 will not work.

 Should be:

 {{{
 g.proj epsg=5514 datum_trans=-1
 ---
 1
 Used in whole hermannskogel region
 towgs84=653.000,-212.000,449.000
 Default 3-Parameter Transformation (May not be optimum for older datums;
 use this only if no more appropriate options are available.)
 ---
 2
 Used in Austria
 towgs84=577.326,90.129,463.919,5.1366,1.4742,5.2970,2.4232
 Accuracy approx. 1.5m
 ---
 3
 Used in Czech Republic
 towgs84=570.8,85.7,462.8,4.998,1.587,5.261,3.56

 ---
 4
 Used in Slovakia
 towgs84=485.021,169.465,483.839,7.786342,4.397554,4.102655,0

 ---
 5
 Used in Slovenia
 towgs84=426.9,142.6,460.1,4.91,4.49,-12.42,17.1
 }}}

-- 
Ticket URL: <http://trac.osgeo.org/grass/ticket/2456#comment:6>
GRASS GIS <http://grass.osgeo.org>