Géoréférencement de carte

[gibtoul]

bonjour,

On m'a fourni des fichiers *.png avec un fichier text contenant les coordonnées géographique de chaque fichier.
Ex le fichier paris_132778_90181_18.png je retrouve sur le fichier texte :
ImageFileName Left_Edge_Longitude Right_Edge_Longitude Top_Edge_Latitude Bottom_Edge_Latitude
paris_132778_90181_18.png: 2.34283447265625 2.34420776367188 48.86019753088974 48.85929403653233
paris_132778_90182_18.png: 2.34283447265625 2.34420776367188 48.85929403653233 48.85839052586621
paris_132778_90183_18.png: 2.34283447265625 2.34420776367188 48.85839052586621 48.85748699889132

Est- il possible avec qgis de georeferencé correctement ces images

merci pour votre aide

[SkåL]

La réponse est "oui".

http://download.osgeo.org/qgis/doc/m...r_guide_fr.pdf

Page 83 : charger un raster dans QGIS,
Page 167 : utiliser l'extension de géoréférencement.

Les coordonnées fournies sont (très certainement) en WGS84 et représentent, comme l'indique la première ligne, les coordonnées gauche, droite, haute et basse de l'image.

Du coup, quand vous souhaiterez géoréférencer votre image "paris_132778_90181_18.png" en WGS 84, il faudra faire 4 points, correspondant aux 4 angles de l'image.

Le coin haut gauche aura les coordonnées X de 2.34283447265625 et Y de 48.86019753088974.
Et ainsi de suite pour les 3 autres angles de cette image.

[Barbapapa]

Bonjour,
C'est possible de tout faire automatiquement en générant un fichier de géoréférencement par image (.pgw pour du png). Mais pour cela, il faut faire un peu de programmation en python.
Vous avez tout ce qu'il faut pour le faire dans l'excellent article de gene :
http://www.portailsig.org/content/py...es-tfw-jgwpngw

[gibtoul]

Bonsoir

Merci pour vos réponses

Effectivement une solution "automatique" m'aurait grandement aidé
N'ayant pas les compétences en python , ... je suis complétement perdu
sur le lien fourni

Si éventuellement une petite routine existe dans la communauté , je saurais preneur

merci

a+

[Barbapapa]

La routine se fait en récupérant des éléments de l'article de Gene plus en se référant à la structure de votre fichier texte.

Le but est donc de générer un fichier de géoréférencement standardisé par image png (donc un fichier .pngw de même nom par fichier png contenant 6 lignes comme expliqué dans l'article de Gene).

- ouverture en lecture de votre fichier texte contenant les infos de géoréférencement.
- dans cette boucle récupération des 5 paramètres (nom_fichier, Xmin, Xmax, Ymin et Ymax) à l'aide de la méthode split
- calcul des paramètres du wordile (cf article) grâce aux paramètres de l'image (largeur et hauteur en pixels)
- écriture de ces 6 paramètres dans un nouveau fichier sur 6 lignes.

[gibtoul]

bonsoir barbapapa

Je comprends bien que ma solution est dans le descriptif que tu me communiques.
Mais ..... mes compétences s’arrête la.
J'arrive a faire fonctionner les logiciels mais la programmation ... c'est un autre monde

Je vais essayer de voir si dans mon entourage quelqu'un connait python


Je te remercie pour ton aide

a+

[gene]

Pour un cas comme le votre, il n'y a pas besoin de Python, cela peut être fait à la main (comme je faisais avant). Vous avez les coordonnées des 4 coins et tout ce qu'il manque c'est la taille en pixels de vos fichiers png (en x et y). Après, c'est un simple règle de trois.

[gibtoul]

Bonsoir GENE

Oui mais j'ai plus de 1000 fichiers cela risque d'être long

a+

[Barbapapa]

Avec 1000 fichiers, pas d'autre choix que de créer un script, le plus simple étant en python.

Ceci dit, n'importe quel langage fera l'affaire : du C, du php...
Il s'agit juste de lire un fichier texte, de récupérer des largeurs et hauteurs d'images puis de faire quelques petits calculs et d'écrire le résultat dans des fichiers textes. Ce n'est vraiment pas compliqué pour quelqu’un qui sait un peu programmer.

Il faut juste comprendre le fonctionnement des fichiers tfw, pngw... Et l'article de Gene est là pour çà.

[Lud]

En parallèle de l'article de gene, tu peux peut être te baser sur un script existant (que l'on a développé avec gene) faisant un peu ce genre de chose, les sources étant accessibles ici
http://www.portailsig.org/content/pci2geotiff
plus particulièrement le fichier pci2geotiff_class.py

[jcr83]

Ce qui m'étonne, c'est que personne n'a demandé dans quelle projection était la carte. Il y a de fortes chances pour que ce soit du pseudo-Mercator, et dans ce cas, les méthodes ci-dessus ne conviennent pas du tout.

En effet, le fichier world doit contenir les coordonnées projetées, et non les coordonnées géographiques.

En partant du nom du fichier, qui contient les numéros des dalles, on trouve, pour le coin supérieur gauche du premier fichier :
X = 20037508.343 * (132778 / 2^17 -0.5 ) = 260803.140511764
Y = 20037508.343 * (0.5 - 90181/ 2^17 ) = 6251173.047284035

Si on convertit ces coordonnées pseudo-Mercator en coordonnées géographiques, on retrouve les valeurs indiquées dans le fichier de calibrage :

$ cs2cs -f %12.10f +init=epsg:3857 +to +init=epsg:4326
260803.140511764 6251173.047284035
2.3428344727 48.8601975350 0.0000000000

Attention, ce ne sont pas ces valeurs qu'il faut mettre dans le fichier world, ce sont les coordonnées du centre du pixel supérieur gauche. Il faut donc ajouter (pour l'axe X) ou retrancher (pour l'axe Y) à ces valeurs la moitié de la taille d'un pixel.

Ensuite, un petit script qui extrait les numéros des dalles du nom du fichier et applique les formules ci-dessus fera l'affaire.

[gene]

moi je veux bien mais c'est ce qui est déjà expliqué en grand détail dans le lien fourni pour les scripts en Python...
Le problème pour giptoul c'est Python et non la compréhension du principe

[jcr83]

Ça, je n'en suis pas du tout certain, d'autant plus que la réponse de SkåL et la tienne d'hier à 19h58 laissaient entendre qu'il fallait mettre les coordonnées géographiques dans le fichier world.

[Lud]

Ce qui m'étonne, c'est que personne n'a demandé dans quelle projection était la carte. Il y a de fortes chances pour que ce soit du pseudo-Mercator, et dans ce cas, les méthodes ci-dessus ne conviennent pas du tout.

Juste pour calmer un peu le jeu, là tu suposes des choses et tu as peut être raison mais ça c'est à gibtoul de le dire, on ne peut pas le deviner (ou difficilement) donc à mon avis les méthodes fonctionent, reste à le faire correctement, nous n'avons pas les fichiers, il faut donc faire avec les informations initiales.

[gene]

car, de toutes manières, il manque la taille en pixels des images