Le portail francophone de la géomatique

Bonjour,

Si vous avez installé OGR/GDAL et que votre répertoire contient AUBENAS.png et AUBENAS.tab (en admettant que vous n'ayez pas d'autres fichiers portant le même nom et ayant une autre extension), vous aurez le moyen expliqué ci-dessous.

Démarche :
Les commandes ci-dessous se passent dans le répertoire contenant les fichiers.
Le Tiff gère les TAB avec géoréférencement, on va donc passer par son intermédiaire (cf http://georezo.net/wiki/main:logiciels: … rencement)
Le contenu du WLD (World file) est indépendant du format d'image : c'est l'extension qui change généralement.

Cas pratique

Passer votre png en tif

gdal_translate AUBENAS.png AUBENAS.tif

Remplacer le contenu de votre AUBENAS.tab et manipuler les fichiers

sed -e "s/AUBENAS.png/AUBENAS.tif/g" AUBENAS.tab > tmp.txt
rm AUBENAS.tab
mv tmp.txt AUBENAS.tab

Rebasculer en png pour récupérer les coordonnées depuis le tab

gdal_translate -co "WORLDFILE=YES" AUBENAS.tif AUBENAS_bis.png
rm AUBENAS_bis.png

mv AUBENAS_bis.wld AUBENAS.pngw

Il ne vous reste plus qu'à récupérer le AUBENAS.pngw et le AUBENAS.png d'origine et c'est calé

Merci de confirmer que cela marche (j'ai testé avec un jeu perso mais dans le doute...)

Cordialement

ThomasG

Dernière modification par ThomasG (Wed 07 July 2010 22:09)

Bonjour,

Tout d'abord merci beaucoup pour vos réponses!
J'ai adapté un poil la procédure pour automatiser un peu plus (j'ai plus d'un raster dans mon répertoire) j'ai donc fait.
Par ailleurs j'ai eu quelques difficultés car ma version de GDAL ne semble pas supporter la génération de worldfile pour les tiffs (bizarre non? il est possible que j'ai fait une boulette à la compilation).

Warning 6: Driver GTiff does not support WORLDFILE creation option

Bref, je m'en suis sorti en faisant des pirouettes assez simples (après pas mal de tâtonnements):



La conversion de png vers tif

find ./ -iwholename "*.png" -execdir gdal_translate {} {}.tif \;

Ça marche parfaitement!


Ensuite j'ai choisi de renommer les fichiers tif obtenus plutôt que de remplacer la chaîne de caractères du .tab:

rename -f -v 's/\.png//' *.tif

Je lance la création des fichiers wld:

find ./ -iwholename "*.tif" -execdir gdal_translate -co "WORLDFILE=YES" {} {}_bis.png \;

Je renomme les fichier world:

rename -f -v 's/tif_bis\.wld/pngw/' *.wld

Et enfin je fini le ménage pour ne garder que les png et les pngw:

rm -vrf *.TAB *.tif*

Sinon dans mon tâtonement j'ai aussi essayé ça (c'est moins pratique je trouve):

find ./ -iwholename "*.png" -execdir gdal_translate {} {}.tiff \; #conversion png vers tiff
for i in `find ./ -iwholename "*.tab"|xargs`; do sed -e "s/png/png.tiff/" $i > $i.txt; done #remplacement de la chaine de caractères de .tab
rename -v -f 's/TAB\.txt/TAB/' * #rempalcement des anciens .tab
find ./ -iwholename "*.tiff" -execdir gdal_translate -co "WORLDFILE=YES" {} {}_bis.png \; #création des wld

Le renommage des fichiers et le ménage reste identique sur la fin.

C'est encore un peu embrouillé dans ma tête.
Par contre j'aurais une autre question:

Comment fais tu le calcul de la taille des pixels?

Salut tous,

Du nouveau coté forumsig.org avec un petit éclaircissement sur la manière dont est structurée un fichier tab (il devient alors possible de faire un script python plutot que de passer par GDAL/bash):
http://www.forumsig.org/showthread.php? … post242890

Ma question aurait était suivante sans avoir compris comment était fait un tab: "comment devines tu que ton image fait 4000px pour 10000m de coté?" Maintenant que j'ai compris comment est structuré un tab je comprend comment tu peux calculer ces infos.

Donc là si j'ai bien tout suivi, on a 10000m de coté car change angle est séparé par 10000m (jusque là tout va bien!) car 77000 - 76000 = 10000 idem pour les y: 19600000 - 1950000 = 10000 soit un carré de 10 000m de coté.
L'unité est précisée dans la dernière ligne de mon tab si je ne me trompe pas.

Idem pour les 4000 px (Le nombre pixel est inscrit dans le tab mais je ne savais extraire cette info).
(0,0) --> coin inférieur gauche
(4000,0) --> coin inférieur droit
(0,4000) --> coin supérieur gauche
(4000,4000) --> coin supérieur droit

Peux tu me confirmer que c'était le raisonnement attendu stp (histoire d'avoir l'air moins c*n la prochaine fois!)?

EDIT:
La commande gdalinfo donne pas mal de ces infos notamment: "Size is 4000, 4000".
Pour la synthaxe: gdalinfo monraster.truc

Avec les explications de gene (voir  le lien précédent vers le forumsig) oui j'ai compris que

(760000.0,1960000.0) (0,0) Label "Pt1"

se rapporte à l'origine, et que 

(770000.0,1950000.0) (4000,4000) Label "Pt4"

se rapporte au coin en haut à droite. Comme les worldfile ne calculent pas de la même manière, il faut compenser le fait que l'origine est au centre du pixel supérieur gauche, c'est bien ça?
Tout ça crééé un décalage c'est ce qui est expliqué ici : http://www.portailsig.org/content/pytho … fw-jgwpngw (paragraphe taille d'un pixel).

Si tu peux m'indiquer si mon raisonnement tiens la route, ça m'aidera à progresser.

Merci d'avance.

PS: Je suis débutant et je n'ai suivi aucune formation en géomatique. Étant autodidacte il est possible que j'ai des lacunes majeures (exemple: ne pas savoir "lire" un fichier .tab) qui m'empêchent de voir des "évidences", donc soyez indulgents svp!

Dernière modification par Matmat (Fri 09 July 2010 11:54)