[gdal-dev] non-fully conform output of gdalwarp

[Even Rouault]

Matthieu,

I can also reproduce it with GDAL 1.6.1, but no longer with GDAL 1.7.0dev. 
Significatn changes have been done in trunk concerning preservation of 
citations in the GeoTIFF driver. The most significant for your use case is 
probably http://trac.osgeo.org/gdal/changeset/15993 that adds Imagine 
citation support. Not sure how other software will like it.

Here's the output with GDAL 1.7.0dev

$ 
gdalwarp -t_srs 'PROJCS["ETRS89_GRS80_3",GEOGCS["LUA_ETRS89",DATUM["D_GRS_1980",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",3500000.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",15.0],PARAMETER["Scale_Factor",0.9996],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0]]' 
byte.tif byte_etrs.tif

$ gdalinfo byte_etrs.tif 
Driver: GTiff/GeoTIFF
Files: byte_etrs.tif
Size is 23, 44
Coordinate System is:
PROJCS["ETRS89_GRS80_3",
    GEOGCS["LUA_ETRS89",
        DATUM["D_GRS_1980",
            SPHEROID["GRS_1980",6378137,298.257222101]],
        PRIMEM["Greenwich",0],
        UNIT["degree",0.0174532925199433]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",15],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",3500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]]]
Origin = (5754460.888622452504933,35335406.682907812297344)
Pixel Size = (999.222226262158188,-999.222226262158188)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  ( 5754460.889,35335406.683) ( 15d 0'0.00"E, 90d 0'0.00"N)
Lower Left  ( 5754460.889,35291440.905) ( 15d 0'0.00"E, 90d 0'0.00"N)
Upper Right ( 5777443.000,35335406.683) ( 15d 0'0.00"E, 90d 0'0.00"N)
Lower Right ( 5777443.000,35291440.905) ( 15d 0'0.00"E, 90d 0'0.00"N)
Center      ( 5765951.944,35313423.794) ( 15d 0'0.00"E, 90d 0'0.00"N)
Band 1 Block=23x44 Type=Byte, ColorInterp=Gray

$ tiffdump -m 200 byte_etrs.tif
byte_etrs.tif:
Magic: 0x4949 <little-endian> Version: 0x2a
Directory 0: offset 1158 (0x486) next 0 (0)
ImageWidth (256) SHORT (3) 1<23>
ImageLength (257) SHORT (3) 1<44>
BitsPerSample (258) SHORT (3) 1<8>
Compression (259) SHORT (3) 1<1>
Photometric (262) SHORT (3) 1<1>
StripOffsets (273) LONG (4) 1<146>
SamplesPerPixel (277) SHORT (3) 1<1>
RowsPerStrip (278) SHORT (3) 1<44>
StripByteCounts (279) LONG (4) 1<1012>
PlanarConfig (284) SHORT (3) 1<1>
SampleFormat (339) SHORT (3) 1<1>
33550 (0x830e) DOUBLE (12) 3<999.222 999.222 0>
33922 (0x8482) DOUBLE (12) 6<0 0 0 5.75446e+06 3.53354e+07 0>
34735 (0x87af) SHORT (3) 88<1 1 0 20 1024 0 1 1 1025 0 1 1 1026 34737 15 0 
2048 0 1 32767 2049 34737 82 15 2050 0 1 32767 2054 0 1 9102 2056 0 1 32767 
2057 34736 1 5 2059 34736 1 6 2061 34736 1 7 3072 0 1 32767 3074 0 1 32767 
3075 0 1 1 3076 0 1 9001 3080 34736 1 1 3081 34736 1 0 3082 34736 1 3 3083 
34736 1 4 3092 34736 1 2 0 0 0 0>
34736 (0x87b0) DOUBLE (12) 8<0 15 0.9996 3.5e+06 0 6.37814e+06 298.257 0>
34737 (0x87b1) ASCII (2) 98<ETRS89_GRS80_3|GCS Name = LUA_ETRS89|Datum = 
D_GRS_1980|Ellipsoid = GRS_1980|Primem = Greenwich||\0>


Le Friday 29 May 2009 12:11:53 Matthieu Rigal, vous avez écrit :
> 'PROJCS["ETRS89_GRS80_3",GEOGCS["LUA_ETRS89",DATUM["D_GRS_1980",SPHEROID["G
>RS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0
>174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting
>",3500000.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",15
>.0],PARAMETER["Scale_Factor",0.9996],PARAMETER["Latitude_Of_Origin",0.0],UNI
>T["Meter",1.0]]' 8_WGS84.tif 8_ETRS89.tif