[OSM-talk-fr] décalage de projection orthophoto

[Vincent de Chateau-Thierry]

Bonsoir,

Le 10/12/2010 23:39, jean navarro a écrit :
>
> -----BEGIN PGP SIGNED MESSAGE-----
> Hash: SHA1
>
> bonjour
>
> j'ai obtenu une orthophoto de ma commune (avec les droits pour
> l'utilisation pour un site web communal)
>
> c'est un fichier .ecw en lambert II centre.
> j'utilise la commande "gdal2tiles.py nom_fichier.ecw" qui me produit des
> dalles au format png et un fichier openlayers.html...
>
> dans un navigateur, je vois bien l'orthophoto qui peut être superposée à
> une carte google... mais cette orthophoto apparait décalé d'environ 40m
> vers l'est !!
>
> dans le fichier html il y a :
>    projection: new OpenLayers.Projection("EPSG:900913"),
>    displayProjection: new OpenLayers.Projection("EPSG:4326"),
>
>
> besoin d'un peu d'aide pour trouver l'origine du pb :
> 	-manque une ou des options sur la commande gdal2tiles ?
> 	-un pb d'origine de la carte mal définie ?
>

La distance (40m) me fait penser à un problème de système géodésique 
("datum"). Ton ortho d'origine en Lambert II est normalement basée sur 
le système NTF (Nouvelle Triangulation Française), et non sur WGS84, que 
tu utilises au travers de la projection 900913. Il faudrait qu'au moment 
où tes données changent de projection, ou avant dans une étape à part, 
tu appliques aussi un changement de datum, pour passer de NTF à WGS84.
Je ne connais pas gdal2tiles.py autrement que par cette page [1] à 
l'instant, où je lis :
"Les fichiers world et les références spatiales incluses sont utilisés 
durant la génération des tuiles"
donc voir si la définition de la projection Lambert II de ton ECW 
indique bien le NTF, et voir ensuite si cet aspect (le changement de 
datum) est bien géré par gdal.

vincent

[1] : http://georezo.net/wiki/main/logiciels/gdal/couteau_suisse/gdal2tiles