[gdal-dev] UTM -> LAEA

[Thomas Becker]

Hello listers,

using gdalwarp to transform a Landsat image from UTM to lambert Azimuthal Equal Area, the output is showing "unknown" for the Datum and "unnamed" for the Spheroid.

I was using the following command:
gdalwarp -t_srs '+proj=laea +lat_0=52 +lon_0=10 +x_0=4321000 +y_0=3210000 +ellps=GRS80' p193r002_7t19990711_z33_nn40.tif test.tif

The outputof gdalinfo:

$ gdalinfo.exe /g/etm+/utm/test.tif
Driver: GTiff/GeoTIFF
Files: g:/etm+/utm/test.tif
Size is 9601, 8837
Coordinate System is:
PROJCS["unnamed",
    GEOGCS["GRS 1980(IUGG, 1980)",
        DATUM["unknown",
            SPHEROID["unnamed",6378137,298.2572221010042]],
        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,
        AUTHORITY["EPSG","9001"]]]
Origin = (4451567.686025693100000,3623709.044894170500000)
Pixel Size = (28.522398941536249,-28.522398941536249)
Metadata:
  AREA_OR_POINT=Area
Corner Coordinates:
Upper Left  ( 4451567.686, 3623709.045)
Lower Left  ( 4451567.686, 3371656.605)
Upper Right ( 4725411.238, 3623709.045)
Lower Right ( 4725411.238, 3371656.605)
Center      ( 4588489.462, 3497682.825)
Band 1 Block=9601x1 Type=Byte, ColorInterp=Gray

Can you tell me how i have to set the parameters for the DATUM and SPHEROID? Wenn i am using the EPSG:3035 ther will be no adjustments for lat_0 and so on. Is there a possibility to use these parameters as well with the EPSG Code?

Cheers,
Thomas
-- 
Psssst! Schon vom neuen GMX MultiMessenger gehört?
Der kann`s mit allen: http://www.gmx.net/de/go/multimessenger