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