[GRASS-dev] Projection name lookup table needed?

[Markus Neteler]

Hi,

I have generated a new location from the EU DEM 25m which comes with
these metadata:

gdalinfo eudem_dem_3035_europe.tif
Driver: GTiff/GeoTIFF
Files: eudem_dem_3035_europe.tif
       eudem_dem_3035_europe.tif.ovr
Size is 240000, 200000
Coordinate System is:
PROJCS["ETRS89 / LAEA Europe",
    GEOGCS["ETRS89",
        DATUM["European_Terrestrial_Reference_System_1989",
            SPHEROID["GRS 1980",6378137,298.2572221010002,
                AUTHORITY["EPSG","7019"]],
            AUTHORITY["EPSG","6258"]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433],
        AUTHORITY["EPSG","4258"]],
    PROJECTION["Lambert_Azimuthal_Equal_Area"],
    PARAMETER["latitude_of_center",52],
    PARAMETER["longitude_of_center",10],
    PARAMETER["false_easting",4321000],
    PARAMETER["false_northing",3210000],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","3035"]]
Origin = (2000000.000000000000000,6000000.000000000000000)
Pixel Size = (25.000000000000000,-25.000000000000000)
[...]

I have then imported the file, cut out a portion with r.mapcalc and
packed up the result with r.pack since I have to deliver the data to
another GRASS user.

For a test, I created a new location with

grass70 -c epsg:3035 ~/grassdata/eu_laea_extra

and tried to import it again. Surprisingly, different spelling created an issue:

r.unpack subregion_eudem_25m.pack
WARNING: Difference between PROJ_INFO file of packed map and of current
         location:
         - name: ETRS89 / LAEA Europe
         ? -------
         + name: ETRS89 / ETRS-LAEA
         ? +++++
         datum: etrs89
         ellps: grs80
         proj: laea
         lat_0: 52
         lon_0: 10
         x_0: 4321000
         y_0: 3210000
         - towgs84: 0,0,0,0,0,0,0
         no_defs: defined
ERROR: Projection of dataset does not appear to match current location. In
        case of no significant differences in the projection definitions,
        use the -o flag to ignore them and use current location definition.

Question:
in lib/proj/convert.c we have a function papszDatumEquiv(). Do we also
need a similar function for names?

thanks,
Markus

PS: If it matters, here the GDAL version (a Linux Mint box):
gdalinfo --version
GDAL 1.10.1, released 2013/08/26