[GRASS-user] Projection problems with EU soil data

[Francesco P. Lovergine]

The ESDB data here available:

http://eusoils.jrc.ec.europa.eu/ESDB_Archive/ESDB_Data_Distribution/ESDB_data.html

still have a peculiar LAEA projection with spherical ellipsoid as 
explained here:

http://eusoils.jrc.ec.europa.eu/gisco_dbm/dbm/home.htm 
(see The Lambert Azimuthal Equal Area projection)
and indeed this is the output of ogrinfo about them:

---8<----------------------------------------------------------------------------
INFO: Open of `sgdbe4_0.shp'
      using driver `ESRI Shapefile' successful.

Layer name: sgdbe4_0
Geometry: Polygon
Feature Count: 51729
Extent: (-1594087.250000, -1350119.375000) - (5100808.000000, 7673966.500000)
Layer SRS WKT:
PROJCS["User_Defined_Lambert_Azimuthal_Equal_Area",
    GEOGCS["GCS_User_Defined",
        DATUM["User_Defined",
            SPHEROID["User_Defined_Spheroid",6378388.0,0.0]],
        PRIMEM["Greenwich",0.0],
        UNIT["Degree",0.0174532925199433]],
    PROJECTION["Lambert_Azimuthal_Equal_Area"],
    PARAMETER["False_Easting",0.0],
    PARAMETER["False_Northing",0.0],
    PARAMETER["Central_Meridian",9.0],
    PARAMETER["Latitude_Of_Origin",48.0],
    UNIT["Meter",1.0]]
AREA: Real (13.3)
PERIMETER: Real (13.3)
SOIL_: Real (11.0)
SOIL_ID: Real (11.0)
SMU: Integer (7.0)
STU_DOM: Integer (7.0)
PCAREA: Real (20.6)
---8<----------------------------------------------------------------------------

I created a custom ellipsoid and a correct location to host them as follows:
(g.proj -j output:)

+proj=laea
+lat_0=48.0000000000
+lon_0=9.0000000000
+x_0=0.0000000000
+y_0=0.0000000000
+a=6378388
+b=6378388
+no_defs
+towgs84=0,0,0

I guess, I had to use the same center of WGS84 (maybe not?)

The problem is that the resulting data are south apparently shifted of about 22 km,
which is quite bad for a 1x1 km grid like that.

Anyone had the same experience with those data? Any advice? 


-- 
Francesco P. Lovergine