[GRASS-dev] ETRS_1989_LAEA issue

[Markus Neteler]

Hello Paul,

all the funny non-standard datum names, sigh (official EU data file!).

The original .prj content is:
PROJCS["ETRS_1989_LAEA",GEOGCS["GCS_ETRS_1989",DATUM["D_ETRS_1989",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Lambert_Azimuthal_Equal_Area"],PARAMETER["False_Easting",4321000.0],PARAMETER["False_Northing",3210000.0],PARAMETER["Central_Meridian",10.0],PARAMETER["Latitude_Of_Origin",52.0],UNIT["Meter",1.0]]

Using your recent patch, it seems to work:
g.proj -w
PROJCS["Lambert Azimuthal Equal Area",
    GEOGCS["grs80",
        DATUM["European_Terrestrial_Reference_System_1989",
            SPHEROID["Geodetic_Reference_System_1980",6378137,298.257222101]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433]],
    PROJECTION["Lambert_Azimuthal_Equal_Area"],
    PARAMETER["latitude_of_center",52],
    PARAMETER["longitude_of_center",10],
    PARAMETER["false_easting",4321000],
    PARAMETER["false_northing",3210000],
    UNIT["Meter",1]]

"grs80" was well detected.
Thanks for the quick fix, I have merged it into the 6.2-relbranch for 6.2.1.

Markus

On Tue, Nov 07, 2006 at 11:45:29AM +0000, Paul Kelly wrote:
> Hello Markus
> ETRS_1989 is a non-standard datum name (i.e. not what's in the EPSG 
> database) - it should be European_Terrestrial_Reference_System_1989. I've 
> committed a patch to CVS to handle the non-standard version so you should 
> try it now.
> 
> BTW (not relevant after the above patch, but interesting anyway) with my 
> updates from Sunday night I would think it should at least recognise the 
> ellipsoid, but it will probably put in WGS_1984 instead of GRS_1980 as 
> they are the same within the tolerances we can check - I'm not sure which 
> I should have first in the ellipsoid table (the first matching one will be 
> reported)!
> 
> Paul
> 
> On Tue, 7 Nov 2006, Markus Neteler wrote:
> 
> >Hi,
> >
> >I tried to create a new location from the CLC2000 data from
> >EEA, distributed on DVD. The SHAPE file comes with a .proj file:
> >
> >ogrinfo -so 100kmE58N17.shp 100kmE58N17
> >INFO: Open of `100kmE58N17.shp'
> >     using driver `ESRI Shapefile' successful.
> >
> >Layer name: 100kmE58N17
> >Geometry: Polygon
> >Feature Count: 392
> >Extent: (5800000.000180, 1699999.999117) - (5867022.238090, 1799999.999127)
> >Layer SRS WKT:
> >PROJCS["ETRS_1989_LAEA",
> >   GEOGCS["GCS_ETRS_1989",
> >       DATUM["ETRS_1989",
> >           SPHEROID["GRS_1980",6378137.0,298.257222101]],
> >       PRIMEM["Greenwich",0.0],
> >       UNIT["Degree",0.0174532925199433]],
> >   PROJECTION["Lambert_Azimuthal_Equal_Area"],
> >   PARAMETER["False_Easting",4321000.0],
> >   PARAMETER["False_Northing",3210000.0],
> >   PARAMETER["Central_Meridian",10.0],
> >   PARAMETER["Latitude_Of_Origin",52.0],
> >   UNIT["Meter",1.0]]
> >CODE_00: String (3.0)
> >AREA_HA: Real (17.5)
> >SHAPE_area: Real (19.11)
> >SHAPE_len: Real (19.11)
> >
> >
> >But GRASS creates the location like this:
> >
> >GRASS 6.3.cvs (clc2000):/tmp > g.proj -w
> >PROJCS["Lambert Azimuthal Equal Area",
> >   GEOGCS["unnamed",
> >       DATUM["unknown",
> >           SPHEROID["unnamed",6378137,298.257222101]],
> >       PRIMEM["Greenwich",0],
> >       UNIT["degree",0.0174532925199433]],
> >   PROJECTION["Lambert_Azimuthal_Equal_Area"],
> >   PARAMETER["latitude_of_center",52],
> >   PARAMETER["longitude_of_center",10],
> >   PARAMETER["false_easting",4321000],
> >   PARAMETER["false_northing",3210000],
> >   UNIT["Meter",1]]
> >
> >g.proj -p
> >-PROJ_INFO-------------------------------------------------
> >name       : Lambert Azimuthal Equal Area
> >proj       : laea
> >a          : 6378137.0
> >es         : 0.006694380022900787
> >lat_0      : 52
> >lon_0      : 10
> >x_0        : 4321000
> >y_0        : 3210000
> >no_defs    : defined
> >-PROJ_UNITS------------------------------------------------
> >unit       : Meter
> >units      : Meters
> >meters     : 1
> >
> >
> >What's the trick to make GRASS recognize it completely?
> >
> >thanks,
> >Markus
> >
> >_______________________________________________
> >grass-dev mailing list
> >grass-dev at grass.itc.it
> >http://grass.itc.it/mailman/listinfo/grass-dev
> >

-- 
Markus Neteler  <neteler itc it>  http://mpa.itc.it/markus/
ITC-irst -  Centro per la Ricerca Scientifica e Tecnologica
MPBA - Predictive Models for Biol. & Environ. Data Analysis
Via Sommarive, 18        -       38050 Povo (Trento), Italy