[GRASS-dev] ETRS_1989_LAEA issue

[Paul Kelly]

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
>