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
_______________________________________________
Talk-fr mailing list
Talk-fr@openstreetmap.org
http://lists.openstreetmap.org/listinfo/talk-fr
