[gdal-dev] gdal 1.6 doesn't recognize his own projection definition (LL -> LAEA)

[Markus Neteler]

Hi,

I have taken a couple of CGIAR SRTM Geotiffs (LatLong, unfortunately
they missed to insert that into the metadata in V4) and gdalwarp'ed them
to the EU LAEA projection on the fly. The definition originates from
the CORINE landcover mapset available at the European Environmental
Agency, .prj file is attached:

gdalwarp -s_srs epsg:4326 -t_srs /geo/laea_etrs89.prj
/geo/cgiar_LL/srtm_41_0[1-3].tif srtmV4_cgiar_laea.tif
Creating output file that is 8277P x 22800L.

But recognition then fails:

gdalinfo europe_all_srtmV4_cgiar_250m_LAEA_EU.tif
Driver: GTiff/GeoTIFF
Files: europe_all_srtmV4_cgiar_250m_LAEA_EU.tif
Size is 24213, 17341
Coordinate System is:
PROJCS["ETRS_1989_LAEA_L52_M10",
    GEOGCS["GCS_ETRS_1989",
        DATUM["unknown",
            SPHEROID["unnamed",6378137,298.2572221010002]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433]],
    PROJECTION["Lambert_Azimuthal_Equal_Area"],
    PARAMETER["latitude_of_center",0],
    PARAMETER["longitude_of_center",0],
    PARAMETER["false_easting",4321000],
    PARAMETER["false_northing",3210000],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]]]
Origin = (1772847.359583478420973,4574867.517244365997612)
Pixel Size = (250.000000000000000,-250.000000000000000)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  ( 1772847.360, 4574867.517) ( 23d26'46.28"W, 12d 6'34.61"N)
Lower Left  ( 1772847.360,  239617.517) ( 25d 6'35.90"W, 26d28'25.00"S)
Upper Right ( 7826097.360, 4574867.517) ( 32d26'40.92"E, 11d52'36.96"N)
Lower Right ( 7826097.360,  239617.517) ( 34d42'27.86"E, 25d54'34.36"S)
Center      ( 4799472.360, 2407242.517) (  4d19'30.13"E,  7d15'33.16"S)
Band 1 Block=24213x1 Type=Int16, ColorInterp=Gray
  Overviews: 12107x8671, 6054x4336, 3027x2168, 1514x1084, 757x542

Apparently the datum isn't properly morphed? I was using 1.6-branch to generate
the map and used both 1.6-branch and SVN trunk to read it, no way.

Bug or feature?

Ah, generating it in SVN trunk and reading it in SVN trunk appears to work:
...
Coordinate System is:
PROJCS["ETRS_1989_LAEA_L52_M10",
    GEOGCS["GCS_ETRS_1989",
        DATUM["D_ETRS_1989",
            SPHEROID["GRS_1980",6378137,298.257222101]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433]],
    PROJECTION["Lambert_Azimuthal_Equal_Area"],
...

Since I prefer to use for the production machine the 1.6-branch version,
could the relevant fix be backported?

thanks,
Markus
-------------- next part --------------
A non-text attachment was scrubbed...
Name: laea_etrs89.prj
Type: application/octet-stream
Size: 388 bytes
Desc: not available
Url : http://lists.osgeo.org/pipermail/gdal-dev/attachments/20090205/441be19c/laea_etrs89.obj