[gdal-dev] Question about Ticket #2706

Jean-Claude REPETTO jrepetto at free.fr
Wed Apr 28 08:48:37 EDT 2010


I have a problem that seems related to ticket #2706. Has this bug been 
fixed in GDAL 1.7.2 ?

I have created a Tiff file with this WKT :
PROJCS["Geoportail - France metropolitaine",GEOGCS["Reseau geodesique 
francais 1993",DATUM["Reseau geodesique francais 1993",SPHEROID["IAG GRS 
1980",6378137,298.2572221010000,AUTHORITY["IGNF","ELG037"]],TOWGS84[0,0,0,0,0,0,0],AUTHORITY["IGNF","REG024"]],PRIMEM["Greenwich",0,AUTHORITY["IGNF","LGO01"]],UNIT["degree",0.01745329251994330],AXIS["Longitude",EAST],AXIS["Latitude",NORTH],AUTHORITY["IGNF","RGF93G"]],PROJECTION["Equirectangular",AUTHORITY["IGNF","PRC9001"]],PARAMETER["semi_major",6378137],PARAMETER["semi_minor",6378137],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",0],PARAMETER["standard_parallel_1",46.5],PARAMETER["false_easting",0.000],PARAMETER["false_northing",0.000],UNIT["metre",1],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["IGNF","GEOPORTALFXX"]]

The true_scale latitude (46.5°) seems to be present in the Tiff file :
$ tiffinfo /tmp/test.tif
   Tag 34735: 
1,1,0,20,1024,0,1,1,1025,0,1,1,1026,34737,35,0,2048,0,1,32767,2049,34737,128,35,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,17,3076,0,1,9001,3078,34736,1,2,3082,34736,1,3,3083,34736,1,4,3088,34736,1,1,3089,34736,1,0,0,0,0,0
   Tag 34736: 
0.000000,0.000000,46.500000,0.000000,0.000000,6378137.000000,298.257222,0.000000

but not in gdalinfo :

Coordinate System is:
PROJCS["Geoportail - France metropolitaine",
     GEOGCS["Reseau geodesique francais 1993",
         DATUM["Reseau_geodesique_francais_1993",
             SPHEROID["IAG GRS 1980",6378137,298.257222101]],
         PRIMEM["Greenwich",0],
         UNIT["degree",0.0174532925199433]],
     PROJECTION["Equirectangular"],
     PARAMETER["latitude_of_origin",0],
     PARAMETER["central_meridian",0],
     PARAMETER["standard_parallel_1",0],
     PARAMETER["false_easting",0],
     PARAMETER["false_northing",0],
     UNIT["metre",1,
         AUTHORITY["EPSG","9001"]]]

FYI, I am using libtiff 3.9.2, not the internal Tiff library.

Regards,
Jean-Claude



More information about the gdal-dev mailing list