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 _______________________________________________ grass-user mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/grass-user
