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
[pei...@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)
(...)
[pei...@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.
[pei...@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)
(...)
_______________________________________________
gdal-dev mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/gdal-dev
