Re: [OSM-talk-fr] Projet d'import WMS du cadastre avec le s rasters
Salut Fabien, Le script de Fabien est dispo ici: http://crapouillou.net/~jocelyn/ressources/osm/cadastre_tools_180608.tar.bz2 J'ai modifié le Wiki vers ce nouvel emplacement C'est du Python donc c'est Bon ;-) Petite remarque, tu peux dans ce script comme dans l'outil d'import de Fred extraire juste le calque qui t'intéresse grace aux arguments layer et style. Au max tu as ça: layer=CDIF:LS3,CDIF:LS2,CDIF:LS1,CDIF:PARCELLE,CDIF:NUMERO,CDIF:PT3,CDIF:PT2,CDIF:PT1,CDIF:LIEUDIT,CDIF:SUBSECTION,CDIF:SECTION,CDIF:COMMUNE style=LS3_90,LS2_90,LS1_90,PARCELLE_90,NUMERO_90,PT3_90,PT2_90,PT1_90,LIEUDIT_90,SUBSECTION_90,SECTION_90,COMMUNE_90 Si tu veux juste le contour de commune tu peux faire: layer=CDIF:COMMUNE style=COMMUNE_90 Pour le détail de chaque valeur voir Wiki Cadastre Techniques: http://wiki.openstreetmap.org/wiki/WikiProject_Cadastre_Fran%C3% A7ais/Technique#Les_diff.C3.A9rents_calques Attention aussi aux communes qui sont aux limites de deux zones Lambert PS: Je fais mon petit retour sur la liste après un LONG congée paternité, salut Pierre ! ;-) signature.asc Description: This is a digitally signed message part ___ Talk-fr mailing list Talk-fr@openstreetmap.org http://lists.openstreetmap.org/listinfo/talk-fr
[OSM-talk-fr] Regarde mon profil Facebook
Bonjour Liste, Je vous ai invité à rejoindre Facebook récemment et je voulais vous rappeler que dès que vous serez enregistré, vous pourrez vous connecter en ligne, partager des photos, organiser des groupes et des événements, et plus. Merci, Ratzillas Pour vous inscrire à Facebook, suivez le lien ci-dessous : http://www.facebook.com/p.php?i=636867474k=SZE35YRX4W5M5J1CWGWYSVr Ratzillas Da Rats#039; Prince a invité talk-fr@openstreetmap.org à rejoindre Facebook. Si vous ne voulez plus recevoir ce type de messages de la part de Facebook, veuillez cliquer sur le lien ci-dessous pour annuler votre inscription. http://www.facebook.com/o.php?k=07787bu=1843658522mid=885a13G6de3ff1aG0G8 Les bureaux de Facebook se trouvent à : 1601 S. California Ave., Palo Alto, CA 94304. ___ Talk-fr mailing list Talk-fr@openstreetmap.org http://lists.openstreetmap.org/listinfo/talk-fr
[OSM-talk-fr] Résolution décalage Cadastre DOM-TOM
Bonjour à tous je vous livre les résultats de mes investigations pour le
cadastre dans les DOM. J'ai fais des vérifications graphiques sous
Quantum Gis d'un fichier échantillons de points extrait du cadastre en
ligne. La bibliothèque libre proj4 (dispo chez nos amis de l'OsGeo)
permet de configurer les paramètres de projections.
Dans cette bibliothèque on retrouve 'cs2cs' qui permet de convertir des
coordonnées entre deux systèmes géodésiques différents. L'avantage c'est
que cet utilitaire peut être très finement configuré (Ellipsoïde de
référence, Système géodésique etc..). Géoconv semble utiliser par défaut
l'ellisoïde du WGS84. Or lorsque l'on regarde les paramètres de
projection pour les DOM (Gpe et Mqe) c'est l'ellisoïde de Hayford qui
est utilisé. (Je n'ai pas réussi à avoir exactement le système IGN, il y
a toujours des décalages)
Mais ce n'est pas la seule explication au décalage des fichiers OSM de
l'outil de Fred
Régardez cette fiche géodésique de l'IGN:
http://geodesie.ign.fr/fiche_point_OM.asp?num_site=9710106no_ptg=03
On va pouvoir tester grâce à ces coordonnées fiables les résultats de
geoconv et de cs2cs
-
Coordonnées Degrés Minutes Secondes
61 ° 29 ' 10.05520 ''Ouest 16 ° 15 ' 40.53420 ''Nord
Coordonnées Degrés Décimaux
61.48612644 16.261259500
Coordonnées XY WGS84
661 775,546 1 798 433,350
Coordonnées XY Système Sainte Anne - Cadastre
662 198,180 1 798 736,709
-
Résultats Obtenus par transformation des Coord XY WGS84
code ligne de commande: echo 661775.546 1798433.350 | cs2cs -f %.16f
+proj=utm +zone=20N +to +proj=latlon -
-61.48612644(45438920)16.261259(4232077470) - OK
Résultats Obtenus par transfo des Coord XY Ste Anne (Syst du Cadastre)
code ligne de commande: echo 662198.180 1798736.709 | cs2cs -f %.16f
+proj=utm +zone=20N +to +proj=latlon -
--Décalage--
-61.482151435039242016.2639725768512022
-
Résultats de GEOCONV sur système WGS84 XY
WGS84 -61.486126444(54062) 16.261259(424403377) OK
--Décalage--
Résultats de GEOCONV sur système Ste Anne XY - CADASTRE
WGS84 -61.48215143503598 16.26397257804702
-
Les outils sont donc fiables puisque l'on a jusqu'à la 7ème décimale les
mêmes résultats.Mais il faut partir des coordonnées XY WGS84.J'ai donc
effectué la différence entre les deux systèmes géodésiques et effectué
les transfos après le calcul suivant:
/
| X_wgs84=Xcadastre-422.634m
| Y_wgs84=Ycadastre-303.359m
\
Piti Schéma : Il faut donc décaler la grille cadastrale de 422.634m à
l'Ouest et de 303.359m au Sud pour avoir la WGS84
- -422m ___C_C: Système Cadastrale
| |W: Système WGS84
-303m __W_| |
|| || |
v| ||__|
||
||
Tous les points géodésiques que j'ai pu controler ont cette
particularité à +-20cm ce qui me fait dire que sur un territoire comme
la Guadeloupe cette transformation est homogène.
J'ai dans la foulée pris des points de la côte Guadeloupéenne dans OSM
pour voir si par transformation inverse:
WGS84 (lat/lon) - UTM20 Nord (X Y)
Je tombais dans l'Océan ou la Mer des Caraïbes ;-( et ... et ... Non ;-)
Donc la côte dans OSM à peu de chose est bien calée, les outils de
conversions sont fiables, le cadastre est fiable ;-$
J'ai donc modifié les scripts de l'Outil d'Import suivants comme suit:
conv-lambert.sh
java -jar tools/geoconv.jar -o WGS84 -deg -in Lambert $lambertZone \
${1} \${2} 0 -out 'trkpt lat=${2.3} lon=${2.2}/trkpt'
${inputFile}
remplacé par:
awk '{printf %.3f %.3f\n,$1-303.4,$2-422.6}' ${inputFile} | cs2cs
-f %.16f +proj=utm +zone=20N +to +proj=latlon - | awk '{print trkpt
lat=\ $2 \ lon=\ $1 \/trkpt }'
conv-lambert-gnuplot.sh
java -jar tools/geoconv.jar -o WGS84 -deg -in UTM 20 N \${1} \${2} 0
-out 'trkpt lat=${2.3} lon=${2.2}/' ${inputFile}
remplacé par
awk '{printf %.3f %.3f\n,$1-303.4,$2-422.6}' ${inputFile} | cs2cs
-f %.16f +proj=utm +zone=20N +to +proj=latlon - | awk '{print trkpt
lat=\ $2 \ lon=\ $1 \/ }'
---
Explications des arguments:
Décalage des Systèmes géodésique
Colonne1($1 X) - 303.4 et Colonne2($2 Y)-422.6:
awk '{printf %.3f %.3f\n,$1-303.4,$2-422.6}' ${inputFile}
-f %.16f : sortie en degrée décimaux 16 décimales
+proj=utm +zone=20N +to +proj=latlon: Système de départ UTM20N
+to +proj=latlon : Système
