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
laea_etrs89.prj
Description: Binary data
_______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
