#1 Thu 08 February 2024 17:53
- BegBlev
- Juste Inscrit !
- Date d'inscription: 8 Feb 2024
- Messages: 3
Probl de conversion RGF93 vers WGS84
Bonjour,
Dans le cadre de mes activités SNSM, je travaille sur un projet de conversion des fichiers Litto3D du SHOM en images puis en fichier ".kap" afin de disposer d'une cartographie fine du littoral utilisable avec des logiciels de navigation type OpenCPN.
Dans ce but j'ai créé un projet open source sur Github:
https://github.com/BegBlev/litto3d-tools
La partie conversion en images donne de très bon résultats mais je rencontre des difficultés concernant les coordonnées.
Je ne parviens pas à convertir correctement les positions RGF93 en positions WGS84.
Si je prends par exemple le fichier fournis par le SHOM:
BZH-MAR_FRA_0218_6877_MNT_20220731_L93_RGF93_IGN69.asc
l'en-tête du fichier indique:
xllcenter 218000.000000
yllcenter 6876001.000000
Si je converti ces valeurs RGF en WGS84 en utilisant différents moyens:
- le programme sous Linux "cs2cs"
- sites de conversion en ligne
Je trouve toujours les mêmes coordonnées
N48.803511
W3.569999
Pour ces conversions j'utilise le système EPSG:2154 (RFG93) en entrée et EPSG:4326 (WGS84) en sortie.
Lorsque je reporte ces coordonnées sur Google Maps ou sur OpenCPN je me retrouve quelques centaines de mètres au sud ouest de la position attendue.
Sauriez-vous d'où peut provenir le problème?
Merci d'avance, Cordialement
Hors ligne
#2 Thu 08 February 2024 19:00
- Yves Egels
- Participant assidu
- Lieu: Paris
- Date d'inscription: 29 Sep 2011
- Messages: 256
- Site web
Re: Probl de conversion RGF93 vers WGS84
Bonjour,
Le problème, c'est qu'il n'y a pas de problème de transformation ... La transfo est correcte (vérifiée avec un outil officiel certifié : Circe développé par mes camarades géodésiens de l'IGN). Le positionnement du point est identique dans Google earth et dans Géoportail (lui aussi officiel et certifié), que ce soit en géographiques ou en Lambert 93. Le problème se situe donc en amont, probablement au niveau du piquage des coordonnées ? Est-ce que le décalage est constant, et uniquement en Y ?
Ingénieur géographe honoraire
École nationale des sciences géographiques
Société française de photogrammétrie et télédétection
Hors ligne
#3 Tue 13 February 2024 10:20
- BegBlev
- Juste Inscrit !
- Date d'inscription: 8 Feb 2024
- Messages: 3
Re: Probl de conversion RGF93 vers WGS84
Bonjour Yves,
Merci pour ta réponse.
Je crois que je viens de comprendre le problème.
Le fichier de données vient du SHOM et me semblait donc, à priori, fiable.
Pour info, le fichier qui contient le maillage d'altitude Lidar vient des données Litto3D du SHOM que l'on récupère ici:
https://diffusion.shom.fr/partie-mariti … -2021.html
Le fichier que j'utilise pour l'exemple ici est:
BZH-MAR_FRA_0218_6877_MNT_20220731_L93_RGF93_IGN69.asc
L'en-tête du fichier contient les coordonnées:
xllcenter 218000.000000
yllcenter 6876001.000000
La spécification technique du fichier peut-être récupérée ici:
https://services.data.shom.fr/static/sp … c_v1_5.pdf
Cette spécification précise les éléments suivants:
- Système géodésique: RGF93
- Projection: Lambert 93
- Les coordonnées de l'en-tête sont celles de l'angle NO (§6.1.2)
Lorsque je transforme le relevé en carte (ci-jointe) je constate que le coin NO est à peut près là:
https://maps.app.goo.gl/RSkkbwuMZ5uieJCe9
Alors que les coordonnées obtenues par conversion m'envoient ici:
https://maps.app.goo.gl/xiXaEK18P9UuGfcL8
En réalité il s'agit du point SO de la carte.
Je viens de faire le test avec un autre fichier et j'obtiens le même résultat.
Il semble donc qu'il y ait une erreur dans la spécification du SHOM.
Merci beaucoup pour ton aide.
Cordialement.
Vincent
Dernière modification par BegBlev (Tue 13 February 2024 10:21)
Hors ligne
#4 Tue 13 February 2024 10:51
- fbecir
- Participant assidu
- Lieu: Saint-Mandé
- Date d'inscription: 16 Sep 2008
- Messages: 512
Re: Probl de conversion RGF93 vers WGS84
Bonjour
Dans le format ASC, xllcenter et yllcenter signifient X Lower Left Center et Y Lower Left Center, donc les coordonnées du centre du pixel le plus au Sud-Ouest.
La spécification du SHOM n'est pas fausse : dans le nom du fichier, on a bien le X ; Y du coin Nord-Ouest (donc le coin du pixel) et dans le fichier, l'indication du coin Sud-Ouest, mais cela c'est dû au format ASC.
(https://desktop.arcgis.com/fr/arcmap/la … format.htm)
Cordialement
Hors ligne
#5 Tue 13 February 2024 10:56
- fbecir
- Participant assidu
- Lieu: Saint-Mandé
- Date d'inscription: 16 Sep 2008
- Messages: 512
Re: Probl de conversion RGF93 vers WGS84
Pour compléter, le SHOM a adopté la nomenclature de nom de fichier utilisée par l'IGN pour ses produits comme la BDORTHO où l'on nomme le fichier à partir des coordonnées du pixel (0 ; 0) de l'image, donc le pixel supérieur gauche.
Dans le cas des fichiers raster, les fichiers de géoréférencement TFW (donc qui viennent d'ESRI comme le ASC), référencent les coordonnées terrain du centre du pixel (0 ; 0), donc le pixel le plus au Nord-Ouest.
Et ce même ESRI, pour les MNT et le format ASC, choisit une autre convention, le pixel en bas à gauche, donc Sud-Ouest ...
Rien n'est simple dans ce bas monde !
Hors ligne
#6 Tue 13 February 2024 11:59
- BegBlev
- Juste Inscrit !
- Date d'inscription: 8 Feb 2024
- Messages: 3
Re: Probl de conversion RGF93 vers WGS84
C'est maintenant clair comme de l'eau de roche
Merci à tous
Vince
Hors ligne