Le portail francophone de la géomatique

Bonjour,
Il suffit d'afficher les informations sue ces dalles (Raster -> Divers -> Informations) pour voir que l'une des dalles ne contient pas d'information de SRS, et l'autre un SRS qui semble erroné, étant données les latitudes affichées :
Corner Coordinates:
Upper Left  (   0.0000000,   0.0000000) (  6d 8'21.03"W, 24d58' 3.34"N)
Lower Left  (       0.000,   -4000.000) (  6d 8'11.52"W, 24d55'46.78"N)
Upper Right (    4000.000,       0.000) (  6d 5'54.10"W, 24d58'12.15"N)
Lower Right (    4000.000,   -4000.000) (  6d 5'44.65"W, 24d55'55.59"N)
Center      (    2000.000,   -2000.000) (  6d 7' 2.83"W, 24d56'59.47"N)

Après plusieurs test, il se trouve que mon fichier TFW est correcte, j'ai réussi à mettre à la bonne place des tiff que j'ai créé.
Du coup je me dit qu'il suffit que je supprime toutes les données de géoréférencement lié à mon tiff, que je mette dans un même dossier mon tif et mon tfw et rien d'autre, j'ouvre qgis et hop tout va aller pour le mieux.
Sauf que non. Bien que dans le dossier il y ai juste les tfw correctement géoréférencé et le tif, bah il me remet mon tiff à la mauvaise place ... alors qu'avec une autre image ça se place bien ...

Comment supprimer ces mauvaise infos pour repartir avec un beau scan25 tout neuf?

Bonjour,

voila un zip avec mes deux .tiff et les deux tfw correspondant:

http://www.mediafire.com/?eoqpgrj4n87dud0

alors j'ai mis dans un dossier mon tiff (le F085_31) tout seul sans le fichier tfw, lorsque je regarde les détails avec gdalinfo, voila ce que j'ai:

Driver: GTiff/GeoTIFF
Files: C:\Users\user21\Desktop\Nouveau dossier (2)\F085_031.TIF
Size is 4000, 4000
Coordinate System is:
PROJCS["NTF (Paris) / France II",
    GEOGCS["NTF (Paris)",
        DATUM["Nouvelle_Triangulation_Francaise_Paris",
            SPHEROID["Clarke 1880 (IGN)",6378249.2,293.4660212936265,
                AUTHORITY["EPSG","7011"]],
            AUTHORITY["EPSG","6807"]],
        PRIMEM["Paris",2.5969213],
        UNIT["grad",0.01570796326794897],
        AUTHORITY["EPSG","4807"]],
    PROJECTION["Lambert_Conformal_Conic_1SP"],
    PARAMETER["latitude_of_origin",46.8],
    PARAMETER["central_meridian",0],
    PARAMETER["scale_factor",0.99987742],
    PARAMETER["false_easting",600000],
    PARAMETER["false_northing",2200000],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","27582"]]
Origin = (0.000000000000000,0.000000000000000)
Pixel Size = (1.000000000000000,-1.000000000000000)
Metadata:
  AREA_OR_POINT=Area
  TIFFTAG_DOCUMENTNAME=F085_031
  TIFFTAG_IMAGEDESCRIPTION=GEOCODE:LAMBERT2
  TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
  TIFFTAG_XRESOLUTION=72
  TIFFTAG_YRESOLUTION=72
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (   0.0000000,   0.0000000) (  6d 8'21.03"W, 24d58' 3.34"N)
Lower Left  (       0.000,   -4000.000) (  6d 8'11.52"W, 24d55'46.78"N)
Upper Right (    4000.000,       0.000) (  6d 5'54.10"W, 24d58'12.15"N)
Lower Right (    4000.000,   -4000.000) (  6d 5'44.65"W, 24d55'55.59"N)
Center      (    2000.000,   -2000.000) (  6d 7' 2.83"W, 24d56'59.47"N)
Band 1 Block=4000x10 Type=Byte, ColorInterp=Palette
  Color Table (RGB with 256 entries)
    0: 255,0,51,255

puis les 255 autres couleurs

du coup on voit bien que le tiff tout seul garde en mémoire le mauvais géoréférencement. comment dire à qgis de supprimer ces infos et de lui dire de regarder le tfw qui contient:

2.500000
-0.000000
-0.000000
-2.500000
850001.250000
2369998.750000

tout en spécifiant qu'on est en lambert II étendu bien sur.

j'ai fait la même chose avec l'autre tiff qui se positionne bien (le F084_31) lorsqu'il y a le fichier tfw dans le meme dossier, voila ce que donne gdalinfo (quand le tiff est tout seul dans un dossier sans le tfw):

Driver: GTiff/GeoTIFF
Files: C:\Users\user21\Desktop\Nouveau dossier (2)\F084_031.TIF
Size is 4000, 4000
Coordinate System is `'
Metadata:
  TIFFTAG_SOFTWARE=Handmade Software, Inc. Image Alchemy v1.9.1
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0, 4000.0)
Upper Right ( 4000.0,    0.0)
Lower Right ( 4000.0, 4000.0)
Center      ( 2000.0, 2000.0)
Band 1 Block=4000x2 Type=Byte, ColorInterp=Palette
  Color Table (RGB with 256 entries)
    0: 0,0,0,255

et voila le contenu du fichier tfw:

2.500
0
0
-2.500
840001.250
2369998.750

pour conclure:
pour le F085_31.tiff, avec ou sans le tfw, le tiff se place à une origine 0,0
pour le F084_31.tiff, sans le tfw, le tiff se place à une origine 0,0 et avec le tfw, le tiff se place à une origine 840000, 2370000.

Bonjour
Le plus simple est d'intégrer les coordonnées dans le fichier GeoTiff, et de ne plus utiliser les fichiers .TFW.
Exemple :
gdal_edit -a_srs epsg:27572 -a_ullr 840000 2370000 850000 2360000 F085_031.TIF

La commande gdal_edit est décrite là : http://www.gdal.org/gdal_edit.html

Voici le résultat de gdalinfo :

Driver: GTiff/GeoTIFF
Files: F085_031.TIF
       F085_031.TFW
Size is 4000, 4000
Coordinate System is:
PROJCS["NTF (Paris) / Lambert zone II",
    GEOGCS["NTF (Paris)",
        DATUM["Nouvelle_Triangulation_Francaise_Paris",
            SPHEROID["Clarke 1880 (IGN)",6378249.2,293.4660212936265,
                AUTHORITY["EPSG","7011"]],
            TOWGS84[-168,-60,320,0,0,0,0],
            AUTHORITY["EPSG","6807"]],
        PRIMEM["Paris",2.5969213],
        UNIT["grad",0.01570796326794897],
        AUTHORITY["EPSG","4807"]],
    PROJECTION["Lambert_Conformal_Conic_1SP"],
    PARAMETER["latitude_of_origin",46.8],
    PARAMETER["central_meridian",0],
    PARAMETER["scale_factor",0.99987742],
    PARAMETER["false_easting",600000],
    PARAMETER["false_northing",2200000],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AUTHORITY["EPSG","27572"]]
Origin = (840000.000000000000000,2370000.000000000000000)
Pixel Size = (2.500000000000000,-2.500000000000000)
Metadata:
  AREA_OR_POINT=Area
  TIFFTAG_DOCUMENTNAME=F085_031
  TIFFTAG_IMAGEDESCRIPTION=GEOCODE:LAMBERT2
  TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)
  TIFFTAG_XRESOLUTION=72
  TIFFTAG_YRESOLUTION=72
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (  840000.000, 2370000.000) (  3d18'11.36"E, 48d27'30.93"N)
Lower Left  (  840000.000, 2360000.000) (  3d17'54.15"E, 48d21'31.19"N)
Upper Right (  850000.000, 2370000.000) (  3d26'26.41"E, 48d27'18.15"N)
Lower Right (  850000.000, 2360000.000) (  3d26' 8.48"E, 48d21'18.43"N)
Center      (  845000.000, 2365000.000) (  3d22'10.10"E, 48d24'24.74"N)
Band 1 Block=4000x10 Type=Byte, ColorInterp=Palette
  Color Table (RGB with 256 entries)

Le gros avantage de gdal_edit par rapport à gdal_translate (méthode proposée par Yves) est qu'il n'y a pas de phase de décompression/recompression de l'image, le résultat est immédiat et l'image n'est pas altérée.

Dernière modification par jcr83 (Fri 08 February 2013 13:56)

Ou pour une seule image on peut éditer manuellement l'en-tête GeoTIFF du .tif avec un éditeur hexa. Si le tfw est OK, virer l'en-tête ou y mettre les bonnes valeurs de calage à l'image du Scan25 qui a une en-tête GeoTIFF correcte.
http://portableapps.com/apps/development/frhed_portable

La commande gdal_edit n'existe qu'à partir de GDAL 1.9, et encore elle était dans le dossier "samples". A partir de la version 1.10, elle est devenue une commande officielle.
Si tu ne la trouves pas sur ton PC, tu peux la télécharger à l'adresse : http://trac.osgeo.org/gdal/browser/trun … al_edit.py
Mais elle ne fonctionnera que si tu utilises la version 1.9 de GDAL, ou une version ultérieure.

Si tu ne veux pas avoir de problème d'indentation, tu dois enregistrer directement le fichier gdal_edit.py , sans passer par un copier-coller (trop risqué).

Et pour l'autre problème, je pense que tu n'as pas choisi la bonne version. Lis la page http://cartometric.com/blog/2011/10/17/ … n-windows/ pour des instructions détaillées.

Cela veut simplement dire que gdal_edit ne trouve pas le fichier .TIF que tu lui indiques. Une piste : ça peut provenir des parenthèses, ce n'est pas un caractère conseillé dans les noms de fichiers ou de dossiers (de même que les espaces et les lettres accentuées).