#1 Mon 30 December 2013 15:21
- benjamin_68
- Participant assidu
- Lieu: Lapoutroie
- Date d'inscription: 3 May 2010
- Messages: 200
décalage de 40/50 mètres en L93
Bonjour,
Cela fait plusieurs fois que je rencontre ce problème sans y trouver de solution.
Lorsque je reprojette une couche en L93 (EPSG 2154), il arrive que celle-ci, sur certains projets, soit décalée de 40 à 50 mètres par rapport aux autres couches. Cela arrive sur ArcGis, Géoconcept...alors que ces projets sont en L93 !
Cette fois-ci, j'ai le problème sur Géoconcept. Quelqu'un pourrait-il m'éclairer ??
Merci d'avance
Benjamin
Hors ligne
#2 Thu 02 January 2014 14:15
- Nicolas Granier
- Participant assidu
- Date d'inscription: 19 Apr 2007
- Messages: 271
Re: décalage de 40/50 mètres en L93
Bonjour,
Ce décalage est classiquement dû à la méthode de projection entre les Datum RGF et NTF (ou entre le Lambert 93 et le Lambert 2 Etendu).
Il existe deux méthode pour reprojeter les données
- méthode par transformation à 3 paramètres
- méthode par transformation à l'aide de la grille IGN à 7 paramètres (http://geodesie.ign.fr/index.php?page=algorithmes)
La transformation à 7 paramètres est de bien meilleure qualité que celle à 3 paramètres. Le décalage suivant les zones géographiques peut être de plusieurs dizaines de mètre.
Arcgis permet lors de la reprojection de passer en paramètre la grille IGN comme cela est détaillé ici:
http://ressources.esrifrance.fr/transfo … gis9x.aspx
A+
Nicolas GRANIER
Hors ligne
#3 Thu 02 January 2014 16:00
- benjamin_68
- Participant assidu
- Lieu: Lapoutroie
- Date d'inscription: 3 May 2010
- Messages: 200
Re: décalage de 40/50 mètres en L93
Merci beaucoup pour la réponse,
Malheureusement, je n'ai pas ArcGis mais QGIS. Or, QGIS ne permet pas de reprojeter des couches vecteur (sauf erreur de ma part ?), du coup j'utilise IGNMap. Mais IGNMap ne permet pas de choisir le type de transformation... Comment faire dans mon cas ?
Hors ligne
#4 Thu 02 January 2014 16:41
- Nicolas Granier
- Participant assidu
- Date d'inscription: 19 Apr 2007
- Messages: 271
Re: décalage de 40/50 mètres en L93
IGN Map (si je ne me trompe pas) fait la projection à 7 paramètres par défaut, IGN oblige.
C'est confirmé ici (http://ignmap.ign.fr/spip.php?article34)
Qgis fait bien de la reprojection vectorielle.
Quid du résultat avec IGNMAP ?
Dernière modification par Nicolas Granier (Thu 02 January 2014 16:46)
Hors ligne
#5 Thu 02 January 2014 18:07
- FERNANE
- Participant actif
- Date d'inscription: 14 Sep 2005
- Messages: 139
Re: décalage de 40/50 mètres en L93
Bonjour,
QGIS permet de reprojeter une couche vecteur, il faut "sauvegarder sous" et choisir le SCR de destination/ liste déroulante/ bouton "parcourir".
Au préalable la couche vecteur doit dans ses propriétés être bien déclaré dans son système de coordonnées.
Meilleurs voeux 2014
cordialement,
Béatrice FERNANE
Hors ligne
#6 Fri 03 January 2014 12:30
- benjamin_68
- Participant assidu
- Lieu: Lapoutroie
- Date d'inscription: 3 May 2010
- Messages: 200
Re: décalage de 40/50 mètres en L93
Nicolas : Je ne suis pas encore en mesure de connaître le résultat donné à l'aide d'IGNMAP : je le saurai la semaine prochaine.
Béatrice : Merci, c'était tout simple, je pensais que cette commande permettait seulement d'affecter un SCR!
Merci encore pour vos réponses.
Meilleurs voeux également.
Hors ligne
#7 Sun 05 January 2014 14:28
- GAILLOT Alexis
- Participant actif
- Lieu: REMIREMONT
- Date d'inscription: 26 Mar 2009
- Messages: 53
Re: décalage de 40/50 mètres en L93
Bonjour,
Si quelqu'un peut confirmer que sous QGIS, lors de la manip suivante :
"sauvegarder sous" et choisir le SCR de destination/ liste déroulante/ bouton "parcourir"
dans le cas du passage de la NTF au RGF93, QGIS utilise bien la grille ...
Sinon, il y a toujours la possibilité de faire la transformation avec ogr2ogr (faire une recherche sur ce forum), et de spécifier explicitement la transformation avec la grille.
Exemple ci-dessous de la transformation d'un shape de Lambert 1 vers Lambert 93 :
ogr2ogr -t_srs "+init=IGNF:LAMB93" -s_srs "+init=IGNF:LAMB1 + wktext" C:\test\outl93.shp C:\test\LIMLPX.shp
On peut tout faire avec ogr2ogr, simplement la prise en main n'est pas très aisée, mais une fois qu'on a dépassé ce cap, après c'est un outil vraiment très efficace (et indépendant de tout logiciel SIG) pour toutes les questions de conversions et/ou transformations, par fichier ou par lot.
Alexis GAILLOT
Hors ligne
#8 Sun 05 January 2014 21:00
- Sylvain_M
- Membre
- Lieu: Lyon
- Date d'inscription: 3 Aug 2009
- Messages: 298
Re: décalage de 40/50 mètres en L93
QGIS prend en compte la grille si on utilise une projection du registre IGNF pour la couche d'origine (par exemple IGNF:LAMBE pour du lambert 2 etendu), pour la couche de destination peut importe le registre IGNF ou EPSG.
c'est facile à vérifier dans la liste des paramètres de la projection d'origine, on doit y trouver "+nadgrids=ntf_r93.gsb,null"
Hors ligne
#9 Mon 06 January 2014 10:23
- benjamin_68
- Participant assidu
- Lieu: Lapoutroie
- Date d'inscription: 3 May 2010
- Messages: 200
Re: décalage de 40/50 mètres en L93
J'ai tenté la reprojection par QGIS et celle-ci ne fonctionne pas. L'erreur suivante apparaît (et ce peut importe les données) :
L'export du fichier vectoriel a échoué.
Erreur : Impossible de transformer un point pendant le dessin de l'entité avec l'id '0'. L'écriture est stoppée. (Exception: Poursuivre la transformation de
(0.120548, 0.840293)
PROJ.4: +proj=lcc +lat_1=47.25 +lat_2=48.75 +lat_0=48 +lon_0=3 +x_0=1700000 +y_0=7200000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs +to +proj=lcc +lat_1=46.8 +lat_0=46.8 +lon_0=0 +k_0=0.99987742 +x_0=600000 +y_0=2200000 +a=6378249.2 +b=6356514.999978254 +nadgrids=ntf_r93.gsb,null +pm=2.337229167 +units=m +no_defs
Erreur: failed to load datum shift file)
Hors ligne
#10 Mon 06 January 2014 17:38
- Sylvain_M
- Membre
- Lieu: Lyon
- Date d'inscription: 3 Aug 2009
- Messages: 298
Re: décalage de 40/50 mètres en L93
L'erreur indiquée sur la dernière ligne (Erreur: failed to load datum shift file = impossible de charger le fichier de la grille de conversion) fait penser à un problème lors de l'installation de QGIS ou un problème de droits d'accès au fichier ...
Hors ligne
