[GRASS-user] EPSG 3035

[Moritz Lennert]

[please keep the discussion on the mailing list]

On 26/03/09 10:34, Laura Poggio wrote:
> I used the "create from georeferenced file" and selected the file I 
> posted the gdalinfo in this mail.

Ok, as Paul pointed out the projection information in that file is 
insufficient, so it is logical that GRASS cannot create the correct 
location, but it should at least say so during the process and probably 
refuse to create the location.

> 
> I have the same output also from the location created using EPSG code 
> and the value 3035.

This is strange, you should get below results. Which versions of GRASS, 
proj and gdal are you using ?

Moritz

> 
> Thank you
> 
> Laura
> 
> 
> 
> 2009/3/26 Moritz Lennert <mlennert at club.worldonline.be 
> <mailto:mlennert at club.worldonline.be>>
> 
>     On 26/03/09 10:03, Laura Poggio wrote:
> 
>         The output of g.proj - p
> 
>         -PROJ_INFO-------------------------------------------------
>         datum      : etrs89
>         ellps      : grs80
>         -PROJ_UNITS------------------------------------------------
>         unit       : metre
>         units      : metres
>         meters     : 1
> 
>         That looks better
> 
> 
>     No, you should have something like this:
> 
>     GRASS 7.0.svn (test3035):~ > g.proj -p
>     -PROJ_INFO-------------------------------------------------
>     name       : Lambert Azimuthal Equal Area
>     proj       : laea
>     datum      : etrs89
>     ellps      : grs80
>     lat_0      : 52
>     lon_0      : 10
>     x_0        : 4321000
>     y_0        : 3210000
>     no_defs    : defined
> 
>     -PROJ_UNITS------------------------------------------------
>     unit       : metre
>     units      : metres
>     meters     : 1
>     GRASS 7.0.svn (test3035):~ > 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["metre",1]]
> 
>     How did you create this location ?
> 
>     Moritz
> 
>