[gdal-dev] gdal_translate changes coordinate system

[Hermann Peifer]

I am trying to create a blank raster file, following the hint at 
http://lists.osgeo.org/pipermail/gdal-dev/2008-February/016061.html

I might be doing something terribly wrong, but after a simple 
gdal_translate, longitude/latitude of origin are changed from 10/52 to 
0/0, see at the end of this post.

I tried to fix the issue with -a_srs EPSG:3035, but things are getting 
worse. Now the coordinate system doesn't have any projection parameters:

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"]]


Any hint?

Hermann



[peifer at whitefish:lceugr100_06]> gdalinfo lceugr100_06_pct.tif
Driver: GTiff/GeoTIFF
Files: lceugr100_06_pct.tif
        lceugr100_06_pct.aux
Size is 67000, 58000
Coordinate System is:
PROJCS["ETRS_1989_LAEA_L52_M10",
     GEOGCS["GCS_ETRS_1989",
         DATUM["European_Terrestrial_Reference_System_1989",
             SPHEROID["GRS_1980",6378137.0,298.257222101]],
         PRIMEM["Greenwich",0.0],
         UNIT["Degree",0.0174532925199433]],
     PROJECTION["Lambert_Azimuthal_Equal_Area"],
     PARAMETER["False_Easting",4321000.0],
     PARAMETER["False_Northing",3210000.0],
     PARAMETER["Central_Meridian",10.0],   <---
     PARAMETER["Latitude_Of_Origin",52.0], <---
     UNIT["Meter",1.0]]
Origin = (800000.000000000000000,6500000.000000000000000)
Pixel Size = (100.000000000000000,-100.000000000000000)

(...)

[peifer at whitefish:lceugr100_06]> gdal_translate lceugr100_06_pct.tif 
-projwin 3900000 3200000 4000000 3100000 -scale 0 255 0 0  template.tif
Input file size is 67000, 58000
Computed -srcwin 31000 33000 1000 1000 from projected window.
0...10...20...30...40...50...60...70...80...90...100 - done.

[peifer at whitefish:lceugr100_06]> gdalinfo template.tif
Driver: GTiff/GeoTIFF
Files: template.tif
Size is 1000, 1000
Coordinate System is:
PROJCS["ETRS_1989_LAEA_L52_M10",
     GEOGCS["GCS_ETRS_1989",
         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]],
     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 = (3900000.000000000000000,3200000.000000000000000)
Pixel Size = (100.000000000000000,-100.000000000000000)
(...)