[GRASS-user] EPSG 3035

[Markus Neteler]

On Wed, Mar 25, 2009 at 11:29 PM, Paul Kelly
<paul-grass at stjohnspoint.co.uk> wrote:
> On Wed, 25 Mar 2009, Moritz Lennert wrote:
>
>> On 25/03/09 18:05, Laura Poggio wrote:
>>>
>>> The output of gdalinfo on the file is
>>>
>>> Driver: GTiff/GeoTIFF                                   Files:
>>> data/JRC/map/DTM/ast1.tif                        Size is 360, 364 Coordinate
>>> System is:                                   PROJCS["ETRS89 / ETRS-LAEA",
>>>                              GEOGCS["ETRS89",
>>> DATUM["European_Terrestrial_Reference_System_1989",
>>>            SPHEROID["GRS 1980",6378137,298.2572221010042,
>>>                AUTHORITY["EPSG","7019"]], AUTHORITY["EPSG","6258"]],
>>> PRIMEM["Greenwich",0], UNIT["degree",0.0174532925199433],
>>> AUTHORITY["EPSG","4258"]],                            UNIT["metre",1,
>>> AUTHORITY["EPSG","9001"]], AUTHORITY["EPSG","3035"]]
>>
>> Mmmh, so coordinate system info is included and GRASS should recognize it.
>
> The WKT is invalid though as it says it's projected, but doesn't include
> either the name of the projection nor any parameters for it.

But if I used OGR's tool, I can generate it, too:


PYTHONPATH=/usr/local/lib64/python2.5/site-packages/osgeo epsg_tr.py 3035
PROJCS["ETRS89 / ETRS-LAEA",
    GEOGCS["ETRS89",
        DATUM["European_Terrestrial_Reference_System_1989",
            SPHEROID["GRS 1980",6378137,298.257222101,
                AUTHORITY["EPSG","7019"]],
            AUTHORITY["EPSG","6258"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.01745329251994328,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4258"]],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    PROJECTION["Lambert_Azimuthal_Equal_Area"],
    PARAMETER["latitude_of_center",52],
    PARAMETER["longitude_of_center",10],
    PARAMETER["false_easting",4321000],
    PARAMETER["false_northing",3210000],
    AUTHORITY["EPSG","3035"],
    AXIS["X",EAST],
    AXIS["Y",NORTH]]

grep 'ETRS89 / ETRS-LAEA' lib/proj/pcs.csv
3035,"ETRS89 / ETRS-LAEA",9001,4258,19986,9820,1,0,8802,10,9102,8801,52,9102,8807,3210000,9001,8806,4321000,9001,,,,,,,,,

?
Markus