Le portail francophone de la géomatique

ftsig31 · Tue 31 January 2012 10:13

Bonjour,

Pour nos projets QGis en Lambert 93, dans préférences, propriétés du projet, j'utilisais jusqu'à présent avec bonheur et sans me poser la moindre question philosophique, le système de coordonnées EPSG 2154, ID 45, RGF93 / Lambert-93.

Puis suite à la lecture du site de l'IGN http://lambert93.ign.fr/index.php?id=29, un petit doute a commencé à s'instiller sournoisement. J'ai ainsi découvert le système IGNF:LAMB93, ID 10089, Lambert 93.

QGis nous propose les 2.
En y regardant plus en détail, il y a tout de même des différences entre les 2 définitions...
(lat_1 et lat_2 inversés, +towgs84 avec 2 paramètres +a et +rf)

Pourquoi le lambert93 a il 2 définitions différentes ??
Lequel faut il prendre?

EPSG 2154, ID 45, RGF93 / Lambert-93 :
+proj=lcc
+lat_1=49
+lat_2=44
+lat_0=46.5
+lon_0=3
+x_0=700000
+y_0=6600000
+ellps=GRS80
+towgs84=0,0,0,0,0,0,0
+units=m
+no_defs

IGNF:LAMB93, ID 10089, Lambert 93 :
+proj=lcc
+towgs84=0.0000,0.0000,0.0000 +a=6378137.0000 +rf=298.2572221010000
+lat_0=46.500000000
+lon_0=3.000000000
+lat_1=44.000000000
+lat_2=49.000000000
+x_0=700000.000
+y_0=6600000.000
+units=m
+no_defs

Je vous remercie

Fabrice

Fichier(s) joint(s) :
Pour accéder aux fichiers vous devez vous inscrire.

Yves · Tue 31 January 2012 10:51

Bonjour,

Le code IGNF:LAMB93 utilise la grille de déformation pour reprojeter les données avec une précision centimétrique. Pour des donnés à petite échelle cela ne pose pas de problème pour des grandes échelles la différence commence à se voir lors des reprojections.

Y.

mich · Tue 31 January 2012 11:08

Bonjour,
Pourriez vous apporter quelques compléments quand je vous cite :

Yves a écrit:

Le code IGNF:LAMB93 utilise la grille de déformation pour reprojeter les données avec une précision centimétrique. Pour des donnés à petite échelle cela ne pose pas de problème pour des grandes échelles la différence commence à se voir lors des reprojections.

Y.

De quel échelle de grandeur parlez vous? Qu'est ce que des données à petites échelles ? et des grandes échelles?

En vous remerçiant

David Michallet

Yves · Tue 31 January 2012 11:16

Bonjour,

une carte à grande échelle présente une petite portion de territoire avec beaucoup de détails;
une carte à petite échelle présente une grande portion de territoire avec peu de détails.

Source : http://www.quebecgeographique.gouv.qc.c … chelle.asp

par exemple :
* grand échelle : 1 / 1 000
* petite échelle : 1 / 50 000 000

Y.

dominique.lys · Tue 31 January 2012 13:55

Bonjour,

Il me semblait que le paramètre proj4 qui permet d'utiliser la grille de transformation est +nadgrids=ntf_r93.gsb or on ne trouve ce paramètre que pour les projections IGNF dans le datum NTF (il est absent pour le datum RGF93). Du coup j'aurais tendance à penser qu'une reprojection de EPSG2154 vers IGNF:LAMBE par exemple prendra en compte la grille NTv2; à vérifier.

Bref, je trouve que ça n'explique pas vraiment les différence constatées.

Vu ici :
+lat_1 Latitude of first standard parallel
+lat_2 Latitude of second standard parallel
+a Semimajor radius of the ellipsoid axis
+rf Reciprocal of the ellipsoid flattening term

Personnellement ça me parle pas beaucoup !

Yves · Tue 31 January 2012 14:24

Bonjour,

Pour reprojeter des données lambert 2 étendue en lambert 93 il faut utiliser le code IGNF du lambert 2 étendue et éventuellement EPSG:2154 pour la cible.

La grille est utilisée pour l'ellipsoide du lambert2étendue (NTF) vers WGS84 utilisé en interne par proj4.

Y.

MathieuR · Tue 31 January 2012 14:34

la définition du lambert 93 par l'IGN d'une part, et par l'EPSG d'autre part n'aboutit à aucune déformation notable lors d'une reprojection. Par ex. lambert 2 etendu -> Lambert 93 IGN ou EPSG => résultats identiques
Par contre, celle du lambert 2 étendu diffère.
Lambert 2 etendu IGN ou EPSG-> Lambert 93 => résultats différents, décalages parfois de plusieurs mètres

Dernière modification par baobazz (Tue 31 January 2012 14:37)

dominique.lys · Tue 31 January 2012 14:37

Donc, le choix entre EPSG 2154 et IGNF LAMB93 n'est pas conditionné par l'utilisation de la grille, alors pourquoi ces différences?

Edit

Bon en fait je crois que c'est assez simple les paramètres +a=6378137.0000 et +rf=298.2572221010000 de IGNF:LAMB93 permettent de définir l’ellipsoïde, ils sont résumés par le paramètre +ellps=GRS80 dans EPSG:2154. Après j'imagine que l'ordre entre lat_1 et lat_2 a peu d'importance.
Bref c'est strictement les même.

Dernière modification par dominique.lys (Tue 31 January 2012 15:19)

ftsig31 · Tue 31 January 2012 16:44

Rebonjour,

Je viens de faire un test concret avec QGis 1.7.3 pour tenter de voir ce que tout ça donne.

J'ai un projet avec l'ortho du RGE d'IGN (ortho fournie par IGN en Lambert 93).

Test1 : reprojection d'un point en WGS84
Le point est dans un fichier csv en WGS84 (EPSG:4326) du genre 43.519236;1.496761.
Préférences, propriétés du projet, activer la projection à la volée est coché.

Je mets le projet en EPSG:2154
J'ajoute la couche du point en WGS84 et je fais une capture écran au 1/50ème (1 carré = 50 cm, cf image jointe)
Je mets le projet en IGNF:LAMB93
J'ajoute la couche du point en WGS84 et je fais une capture écran au 1/50ème (1 carré = 50 cm, cf image jointe)
-> Pas de différence, on dira peut être que l'ellipsoide du WGS84 est compatible avec celui du Lambert 93.
-> Intellectuellement, ça n'est pas satisfaisant, avec des paramètres complètement différents on obtient le même résultat...

Pris dans un élan de curiosité, je fais un second test pour voir si ça jour sur les reprojections qui demandent un changement d'éllipsoïde NTF->Lambert93

Test2 : reprojection d'un point en Lambert 3 Sud
Le point est contenu dans un vieux shapefile en Lambert III Sud (les chefs lieux, reçus de la BDTopo IGN)
Préférences, propriétés du projet, activer la projection à la volée est coché.
Je mets le projet en EPSG:2154
J'ajoute la couche du point en LIII et je fais une capture écran au 1/50ème (1 carré = 50 cm, cf image jointe)
Je mets le projet en IGNF:LAMB93
J'ajoute la couche du point en LIII et je fais une capture écran au 1/50ème (1 carré = 50 cm, cf image jointe)
-> Je ne vois pas du tout la différence d'utilisation de ces 2 projections avec QGIS.
-> Je ne vois pas l'intérêt de le projection IGNF:LAMB93 par rapport à EPSG:2154.
(peut être que d'autres logiciels ne donnent pas le même résultat !)

Ma manip vous semble-t-elle bonne?
Constatez vous la même chose?

Cordialement

Fichier(s) joint(s) :
Pour accéder aux fichiers vous devez vous inscrire.

Yves · Tue 31 January 2012 17:21

Bonjour,

-> Pas de différence, on dira peut être que l'ellipsoide du WGS84 est compatible avec celui du Lambert 93.

C'est en partie l'intérêt du lambert93, en plus d'être compatible avec nos voisins.

1 carré = 50 cm, cf image jointe

Normalement 1 / 50 signifie 1 cm sur la carte = 50 cm sur le terrain. Es tu sure que ton pixel fait 1 cm ou 50 cm (je ne comprend pas trop ce que tu as voulu dire) ?

Test2 : reprojection d'un point en Lambert 3 Sud

Tu as ton projet en EPSG:2154 : ta couche est en lambert 3 sud (façon EPSG) elle est reprojeté en epsg:2154
Tu as ton projet en IGNF:LAMB93 : ta couche est en lambert 3 sud (façon EPSG) elle est reprojeté en IGNF:LAMB93

Comme dit plus haut cela ne changera rien. Pour utiliser la grille il faut définir la projection de ta couche en IGNF:LAMB3 ou IGNF:LAMB3C (je ne sais pas laquelle est la bonne). Ton projet peut être en EPSG ou en IGNF.

l'IGN a fournit un ensemble de code propre aux projections françaises, le code LAMB93 doit y être mais effectivement il n'y a pas de différence (mais je peux me tromper).

Y.

dominique.lys · Tue 31 January 2012 17:25

Attention pour que la grille soit prise en compte lors d'une reprojection à la volée dans QGIS il faut utiliser le plugin Transformation Tools

Voilà une comparaison un peu + simple mais + précise:

Je prend un point quelconque en Lambert 2 étendu : x655000 y2052940
J'utilise Circé pour convertir la coordonnée en Lambert 93 en utilisant la grille (précision centimétrique) et je trouve x703112.347 y6486136.042

Maintenant via FWTools (distribution de GDAL) j'utilise la commande cs2cs pour faire des conversions avec les différentes définitions proj4

cs2cs +init=EPSG:27572 +to +init=EPSG:2154
--> 703112.91 6486136.43

cs2cs +init=EPSG:27572 +to +init=IGNF:LAMB93
--> 703112.91 6486136.43

cs2cs +init=IGNF:LAMBE +to +init=EPSG:2154
--> 703112.35 6486136.04

cs2cs +init=IGNF:LAMBE +to +init=IGNF:LAMB93
--> 703112.35 6486136.04

Conclusion: EPSG:2154 ou IGNF:LAMB93 cela n'a pas d'influence sur le résultat par contre en utilisant IGNF:LAMBE j'ai un résultat identique à Circé alors que sans j'ai pratiquement 50cm d'écart. Selon la localisation du point on doit pouvoir trouver des écart bien + grand.

ftsig31 · Wed 01 February 2012 10:21

OK, merci pour ces infos, les choses sont clarifées, voyons ce que ça donne.

Que je mette le projet QGis en IGNF:LAMB93 ou EPSG:2154 est strictement identique.
(@Yves : la BD Ortho a un résolution de 50 cm, 1 pixel = 1 carré = 50 cm)

Je termine mes tests uniquement à l'aide de QGIS 1.7.3 et la propriété du projet "reprojection à la volée" cochée.
Pas de pb pour l'affichage de données en WGS84.
Je vois bien une différence si j'affiche une couche en Lambert III Sud (-> changement NTF->RGF93).

Mon shapefile en Lambert III provient directement d'une livraison IGN de la BD Topo.
Il a déjà un .PRJ associé qui contient :
PROJCS["NTF_Lambert_Zone_III",GEOGCS["GCS_NTF",DATUM["D_NTF",SPHEROID["Clarke_1880_IGN",6378249.2,293.46602]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Lambert_Conformal_Conic"],PARAMETER["False_Easting",600000.0],PARAMETER["False_Northing",200000.0],PARAMETER["Central_Meridian",2.3372291667],PARAMETER["Standard_Parallel_1",43.1992913889],PARAMETER["Standard_Parallel_2",44.9960938889],PARAMETER["Scale_Factor",1.0],PARAMETER["Latitude_Of_Origin",44.1],UNIT["Meter",1.0]]

Je reprojete ce shapefile avec IGNMap pour avoir une reprojection de référence (Lambert3 -> Lambert93).

J'ouvre le shapefile BDTopo initial (en Lambert III Sud avec son .PRJ associé), le point tombe à environ 50m du point reprojeté avec IGN Map (c'est énorme, bien plus que l'altération linéaire attendue... cf écran joint...).

Je supprime le .PRJ du shapefile en Lambert III Sud.
Dans les options SRC du projet, je coche "Demander le SCR" pour les couches qui n'en ont pas.

J'ouvre le Shapefile et je choisis IGNF:LAMB3 -> le point tombe à 3 mm du point reprojeté avec IGN Map (parfait !).
J'ouvre le shapefile et je choisis EPSG:27573 -> le point n'est pas reprojeté à la volée et s'affiche 3000km au sud.

Bilan de ces tests :
1) Que le projet soit en IGNF:LAMB93 ou EPSG:2154 c'est pareil
2) reprojection à la volée : pas de pb avec des données en WGS84.
3) pour la reprojection à la volée de données en Lambert NTF, prendre les IGNF:LAMBxxx et se méfier des .PRJ proposés par défaut dans certaines livraisons...

Cordialement

Fabrice

Fichier(s) joint(s) :
Pour accéder aux fichiers vous devez vous inscrire.

Yves · Wed 01 February 2012 11:34

Bonjour Fabrice,

Quelques précisions :

Les fichiers prj ne contiennent pas la définition pour utiliser la grille, cela revient à utiliser un code EPSG.

Tu n'es pas obligé de supprimer ton fichier prj si tu veux écraser la projection utilisé par QGIS, change la juste à partir des propriétés de la couche.

Le décalage de 3 000 km est dû à l'utilisation d'un mauvais code. 3 000 km pour du lambert 3 => ce ne serai pas du lambert 3 carto ?

Y.

ftsig31 · Wed 01 February 2012 13:16

Bonjour,

Non pas de 3000000, le shapefile est bien en Lambert 3.
Ex 527000, 145000.

A+

Communauté	Ressources	L'emploi	Le marché	L'association	Fils RSS
Forum Blogs Wiki Planet Sites amis	Agenda Biblio Fiches AMO-CNIG Paris PCGIS Téléchargements	Offres d'emploi Banque de CV Formations Démarche Métiers	Geo-Entreprises Geo-Communiqué Appels d'offres	L'association Qui sommes-nous ? Nous aider Bienfaiteurs Partenaires Historique Statuts	Tout le site Job Géomatique Services GeoBlogs ...	Faire un don Mentions légales

Le portail francophone de la géomatique

Annonce

#1 Tue 31 January 2012 10:13

QGIS 1.7.0 : Lambert 93 EPSG:2154 ou IGNF:LAMB93 ?

#2 Tue 31 January 2012 10:51

Re: QGIS 1.7.0 : Lambert 93 EPSG:2154 ou IGNF:LAMB93 ?

#3 Tue 31 January 2012 11:08

Re: QGIS 1.7.0 : Lambert 93 EPSG:2154 ou IGNF:LAMB93 ?

#4 Tue 31 January 2012 11:16

Re: QGIS 1.7.0 : Lambert 93 EPSG:2154 ou IGNF:LAMB93 ?

#5 Tue 31 January 2012 13:55

Re: QGIS 1.7.0 : Lambert 93 EPSG:2154 ou IGNF:LAMB93 ?

#6 Tue 31 January 2012 14:24

Re: QGIS 1.7.0 : Lambert 93 EPSG:2154 ou IGNF:LAMB93 ?

#7 Tue 31 January 2012 14:34

Re: QGIS 1.7.0 : Lambert 93 EPSG:2154 ou IGNF:LAMB93 ?

#8 Tue 31 January 2012 14:37

Re: QGIS 1.7.0 : Lambert 93 EPSG:2154 ou IGNF:LAMB93 ?

#9 Tue 31 January 2012 16:44

Re: QGIS 1.7.0 : Lambert 93 EPSG:2154 ou IGNF:LAMB93 ?

#10 Tue 31 January 2012 17:21

Re: QGIS 1.7.0 : Lambert 93 EPSG:2154 ou IGNF:LAMB93 ?

#11 Tue 31 January 2012 17:25

Re: QGIS 1.7.0 : Lambert 93 EPSG:2154 ou IGNF:LAMB93 ?

#12 Wed 01 February 2012 10:21

Re: QGIS 1.7.0 : Lambert 93 EPSG:2154 ou IGNF:LAMB93 ?

#13 Wed 01 February 2012 11:34

Re: QGIS 1.7.0 : Lambert 93 EPSG:2154 ou IGNF:LAMB93 ?

#14 Wed 01 February 2012 13:16

Re: QGIS 1.7.0 : Lambert 93 EPSG:2154 ou IGNF:LAMB93 ?

Pied de page des forums