Le portail francophone de la géomatique

Bon, je viens de retrouver ma question de 2015 sur le même sujet (mais avec des TIF) :
http://www.forumsig.org/showthread.php/ … ement-TFWX

J'ai bien peur que la seule solution soit de recréer des Raster rectifiés depuis ArcGis...

Problème : il y a plusieurs dizaines de RASTER, et ça va être long...
Peut-être y aurait-il moyen de scripter ce workflow ?
(je devrais peut-être poster dans la section ArcGis, du coup ...)

Si vous avez des pistes !

Sujet intéressant ! si je comprends bien le j2wx contient la transformation affine (degré 1) approximative, le fichier aux.xml définit les points de contrôle pour une transformation de degré 2 ou 3. La solution proposée dans le lien indiqué par Gene a l'air pas mal : gdal_translate pour assigner les points de contrôle puis gdal_warp pour faire la rectification avec l'option -order pour définir le degré de la transformation.

Si ça marche, un petit script Python pour parser le aux.xml et lancer les commandes gdal fera l'affaire !

Voici une image si quelqu'un peut tester : https://we.tl/t-svARYdrGFf

Pardon Dominique pour le délai de réponse, mais j'étais absent ces derniers jours.

OK, donc c'est récupérable.
Mais j'avoue que ça fait pas mal d'étapes pour mes 30 fichiers, et je ne sais pas trop comment les scripter.
Tu penses que GDAL pourrait faire ça ? (et du coup, aurais tu un exemple pour un fichier ?)

Salut Sylvain

Désolé pour la réponse tardive je n''avais pas vu que tu avais remonté le sujet. J'ai essayé avec ton nouveau fichier, en le convertissant en tiff j'arrive à lui faire lire le aux.xml mais le calage n'est pas correcte. Je pense qu'il y avait déjà ce problème avec le précédent fichier mais je n'avais pas vérifié en détail. En cherchant un peu j'ai pu déterminer comment GDAL lisait le aux.xml et c'est assez bizarre : il considère que les coordonnées images des points de contrôle sont exprimées en pouces et que l'origine de l'image est au point inférieur gauche. J'ai bien vérifié le aux.xml et les coordonnées sont définitivement exprimées en pixel avec un origine fixée à l'angle supérieur gauche.

J'ai écris à la mailing liste gdal-dev pour essayer d'avoir plus d'info :
http://osgeo-org.1560.x6.nabble.com/gda … 98985.html

On verra ce qu'il en ressort mais je pense qu'il va falloir revenir à l'idée du script. J'ai un peu de congés cette semaine j'aurai certainement un peu de temps pour commencer qqchose ça ne devrait pas être trop chronophage.

a+

Dernière modification par dominique.lys (Wed 20 March 2019 12:45)

T'embête pas à tout uploader, le mieux c'est que tu fasses tourner le script de ton côté. J'ai fais un premier jet ici. Il faut fuste modifier le dossier en entrée, les rasters peuvent être en jpg ou jp2, on peut ajouter d'autres extension en modifiant la variable targetExt au début du code. Attention ici tous les rasters sont considérés comme étant en Lambert 2 étendu, mais dans les faits je ne sais pas si ce sera le cas. Donc le code parse le aux.xml pour extraire les points de contrôle et écrire un fichier script système qui sera a exécuter indépendamment par la suite. Ce script lance successivement les commandes gdal_translate et gdalwarp nécessaires pour chaque raster. Ce qui implique que le chemin vers les binaires gdal soit connu. Sous Windows si QGIS est installé alors ils doivent être dans C:/OSGeo4w/bin (ou quelque-chose dans ce goût là), il faut donc ajouter ce chemin dans les variables d'environnement du système.

Par ailleurs le script écrit un fichier *.points qui peut être réutilisé dans le géoréférenceur de QGIS, afin de reprendre les géoréférencements si besoin !

Dernière modification par dominique.lys (Thu 21 March 2019 14:25)