Projection conique conforme de Lambert

La directive européenne INSPIRE préconise pour les échelles inférieures ou égales à 1:500000, l'utilisation de la projection ETRS89-LCC (https://inspire.ec.europa.eu).

Cette nouvelle projection présente un méridien central et deux parallèles automécoïques. Il utilise sur l'ellipsoïde global GRS80, associé à ETRS89.

Depuis que le programme EUREF permet la construction du réseau géodésique international ETRS89, précis et de surcroît compatible avec le système WGS84, la nouvelle projection Lambert 2008 (http://www.ngi.be/) a été adopté pour les projections cartographiques nationales. Le Lambert 2008 ( EPSG:3812) remplace le Lambert 72 (EPSG:31370).

Cette nouvelle projection présente un méridien central et deux parallèles automécoïques. Il utilise sur l'ellipsoïde global GRS80, associé à ETRS89.

Ce nouveau Lambert est basé sur l'ETRS89 et est conforme INSPIRE contrairement au précédent qui était sur DB 1972.

Lambert 93 et WGS84

La projection Lambert93 (projection officielle pour les cartes de France métropolitaine depuis le décret 2000-1276 du 26 décembre 2000) est la projection liée au système géodésique RGF93.

Son principal intérêt réside dans son référentiel RGF93, qui est d'une part commun aux voisins européens de la France car fondé sur ETRS89. Ce système a en commun avec le WGS84 - utilisé notamment par le système GPS de positionnement par satellite - l'ellipsoïde de référence IAG GRS80.

La projection RGF93/Lambert93 a pour référence EPSG:2154 [3] ou IGNF:LAMB93.

« Définition précise de la projection Lambert93 par l'IGN »

En pratique, la projection Lambert 93 est peu utilisée, en partie du fait des altérations linéaires importantes qui y sont associées (2,3 m/km à Dunkerque, 0,60 m/km à Marseille et 2,95 m/km à Bonifacio, légèrement inférieures à celles de la projection Lambert 2 étendue). Pour y remédier, le décret 2006-272 a entériné la création de 9 projections coniques conformes sécantes, couvrant 9 zones du nord au sud, telles que proposées par un rapport du CNIG et adoptées notamment par les géomètres-experts et le service du Cadastre. Elles ont en commun avec le Lambert93 le système géodésique RGF93 et le méridien de référence 3°E (Méridien de Greenwich)[4]^,[5].

Alors que le nombre 93 désigne l'année 1993 dans « L93 », le nombre désigne le parallèle dans les zones CC (version de 2009). Ainsi le « CC43 » est centré sur le parallèle 43°N, entre les latitudes 42.25°N et 43.75°N, en passant par 43.00°N. Une numérotation par zone existe (par exemple Lambert 93 zone 2 pour le CC43), mais est inutile et peut amener des confusions avec les anciennes zones. Dans cet exemple le Lambert CC zone 2 (« RGF-CC2 ») risque d'être confondu avec le Lambert Carto Zone Centre (« NTF-LII ») et le Lambert zone II (« NTF-L02 »).

Il existe une « Recommandation d'usage des différentes projections CC en fonction du département » créée par le Certu et reprise telle quelle par la DGFip (« le cadastre »). Cette recommandation n'est pas contraignante, par exemple au nord de la ville de Troyes, dans le département de l'Aube, le CC48 est recommandé (car la latitude est 48.5°N) mais on peut utiliser le CC49 voisin (celui de la capitale Paris) par commodité, les déformations seront les mêmes entre le CC48 et le CC49 lorsqu'on est à mi-chemin entre 48°N et 49°N. Cependant, utiliser un CC qui n'est pas voisin comme le CC50 (celui de Lille) est autorisé mais n'est pas rigoureux[6]^,[7].

L'Institut national de l'information géographique et forestière et un autre organisme public français, appelé Service de l'information aéronautique commercialisent au format papier des cartes destinées au vol à vue. Elles sont aux normes de l'Organisation de l'aviation civile internationale[8]. Leur projection est Lambert 93 mais leurs coordonnées sont WGS84, coordonnées utilisées notamment par le système GPS de positionnement par satellite, ce qui fait que leurs coordonnées ne doivent pas être confondues avec le Lambert 93 tout court (qui est en coordonnées RGF93), et finalement cette projection est appelée Lambert OACI[9].

Cinq cartes sont commercialisées pour la France métropolitaine, aucune pour l'outre-mer. Sur les 4 cartes au 1/500000 le navigant mesure 10 cm sur le papier qui représentent 50 km (10 x 500000 = 50000000 cm = 50 km). Sur la carte au 1/1000000 le navigant mesure 10 cm sur le papier qui représentent 100 km (10 x 1000000 = 10000000 cm = 100 km).

En anglais le vol à vue est appelé VFR, ce qui fait que l'abréviation Lambert OACI-VFR désigne la carte elle-même et non la projection.

Dans les faits, les projections dites Lambert 4 Zones définies en 1920, associées au système NTF (Nouvelle triangulation de la France), restent encore utilisées.

Leur ellipsoïde de référence est l'ellipsoïde de Clarke 1880 IGN, et leur méridien de référence le méridien de Paris (qui passe au milieu de la façade de l'Observatoire de Paris, soit 2°20´14.025" Est du méridien de Greenwich). Les coordonnées angulaires associées au système NTF sont exprimées en grades (gon) et non en degrés, pour profiter de l'équivalence 1 grade de latitude = 100 km. La France est ainsi découpée en quatre zones qui s'étendent à 150 km au sud et au nord de chaque parallèle d'origine, sauf au nord de la zone 1 (200 km), et au sud de la zone 3 (200 km).

Les projections Lambert associées à la NTF sont traditionnellement définies comme tangentes avec facteur d'échelle.

À cela s'ajoute la projection Lambert Grand Champ (très peu utilisée, citée pour mémoire), qui elle se définit avec des coordonnées angulaires en degrés.

La projection Lambert 93 est une projection « carto » dans le sens où elle s'applique à toute la France métropolitaine et pas seulement à des zones restreintes. Il existait déjà une projection carto avant 1993.

Les projections « Lambert zone Carto » sont constituées de 4 zones couvrant toute la France métropolitaine. Ces projections ont été utilisées par l'IGN depuis 1972, alors que le « Lambert zone » « standard » (zones numérotées de 1 à 4) a été utilisé par les services du Cadastre ou autres utilisateurs locaux (collectivités et géomètres notamment). Ainsi, en Lambert 3, on pouvait trouver Y=3 256 315 sur les cartes IGN et Y=256 315 au Cadastre (le même nombre, mais sans le préfixe 3 indiquant la zone). Cette notation dite « standard » est une survivance d'anciennes définitions IGN, et ne devrait plus être utilisée depuis 1972. Des tentatives de distinguer 1920 et 1972 à l'aide des chiffres romains (Lambert zone III, sans le préfixe) ou arabe (Lambert zone 3, avec le préfixe) ont amené des confusions car elles ne sont pas officielles.

La « Lambert 2 étendue » est tout simplement la « Lambert 2 Carto » : une zone couvrant toute la France métropolitaine. Elle est « dépréciée », c'est-à-dire que pour les nouveaux projets la projection Lambert93 remplace la Lambert 2 étendue, car cette dernière ne présente plus de compatibilité et de précision suffisantes avec le RGF93.

De 1993 à 2006, la projection Lambert93 est en phase de transition et n'est pas obligatoire (l'obligation ne s'applique qu'aux agents de l'état).

En 2015, cette projection est encore utilisée dans de nombreux SIG français, qui ont démarré avant 2006. Il en est exactement de même pour les « Zones » (1920, 1972) utilisée à la place des « Zones CC » (2006), tant que cet usage ne pose aucun problème. Par exemple, l'ensemble des réseaux SFR-Numéricâble sont dessinés en Lambert zones (1972) ou Lambert 2 étendu, car tant que les erreurs de longueur de câble sont seulement de quelques centimètres, cela ne pose pas de problème.

En base de données et en gestion documentaire elle est codée indifféremment « LE », « L2E » (notamment dans Autocad Map 3D), ou « L09 » (pour la distinguer de L02, zone II non étendue). Son EPSG est « 27572 ».

Partant des coordonnées (latitude, longitude) = (${\displaystyle \varphi }$, ${\displaystyle \lambda }$) d'un point du globe supposé être un ellipsoïde de révolution de demi-grand axe ${\displaystyle a}$ et de demi-petit axe ${\displaystyle b}$, calculons ses coordonnées ${\displaystyle (X,Y)}$ sur la carte Lambert. Pour cela, nous allons passer par les coordonnées (${\displaystyle \rho }$, ${\displaystyle \theta }$) du point projeté sur le cône. L'axe des ${\displaystyle X}$ est croissant vers l'est et l'axe des ${\displaystyle Y}$ croissant vers le nord.

L'origine correspond aux valeurs ${\displaystyle \lambda =\lambda _{0}}$ et ${\displaystyle \varphi =\varphi _{0}}$. On a bien pour ce point ${\displaystyle X=X_{0}}$ et ${\displaystyle Y=Y_{0}}$.

Les parallèles sont définis par ${\displaystyle \varphi }$ = constante, d'où ${\displaystyle \rho }$ = constante. Ils sont représentés par des cercles concentriques.

Les méridiens sont définis par ${\displaystyle \lambda }$ = constante, d'où ${\displaystyle \theta }$ = constante. Ils sont représentés par des droites passant par le centre commun des cercles précédents, de coordonnées : ${\displaystyle X=X_{0}}$ et ${\displaystyle Y=Y_{0}+\rho _{0}}$.

Les deux constantes ${\displaystyle n}$ et ${\displaystyle \rho (0)}$ sont calculées en utilisant les deux parallèles sécants de référence. Sur ces parallèles, appelés automécoïques, de latitude ${\displaystyle \varphi _{1}\,}$ et ${\displaystyle \varphi _{2}\,}$, l'échelle est conservée :

On peut donc calculer ${\displaystyle n}$ et ${\displaystyle \rho (0)}$ :

De même que ci-dessus nous allons passer par les coordonnées (${\displaystyle \rho }$, ${\displaystyle \theta }$) du point projeté sur le cône :

• ${\displaystyle \rho ={\sqrt {(X-X_{0})^{2}+(Y_{0}-Y+\rho _{0})^{2}}}}$
• ${\displaystyle \theta =2\tan ^{-1}({\frac {X-X_{0}}{Y_{0}-Y+\rho _{0}+\rho }})}$

Comme on a vu plus haut, on a alors :

• ${\displaystyle \lambda ={\frac {\theta }{n}}+\lambda _{0}\,}$
• ${\displaystyle F(\phi )=\tan({\frac {\phi }{2}}+{\frac {\pi }{4}})\left({\frac {1-e\sin(\phi )}{1+e\sin(\phi )}}\right)^{\frac {e}{2}}=\left({\frac {\rho (0)}{\rho }}\right)^{\frac {1}{n}}}$

Si on appelle F la fonction de gauche, puisque F est strictement décroissante :

${\displaystyle \phi =F^{-1}\left(\left({\frac {\rho (0)}{\rho }}\right)^{\frac {1}{n}}\right)}$

Mais le calcul de ${\displaystyle F^{-1}}$ n'est pas trivial. Pour calculer ${\displaystyle \phi }$ avec une précision arbitraire, puisque F est continue sur son intervalle de définition, on peut utiliser la méthode de dichotomie. On peut aussi, pour plus de rapidité, utiliser l'algorithme proposé par l'IGN :

• ${\displaystyle \phi _{0}=2\tan ^{-1}\left(\left({\frac {\rho (0)}{\rho }}\right)^{\frac {1}{n}}\right)-{\frac {\pi }{2}}}$
• ${\displaystyle \phi _{i}=2\tan ^{-1}\left(\left({\frac {\rho (0)}{\rho }}\right)^{\frac {1}{n}}\left({\frac {1+e\sin(\phi _{i-1})}{1-e\sin(\phi _{i-1})}}\right)^{\frac {e}{2}}\right)-{\frac {\pi }{2}}}$
• jusqu'à ce que ${\displaystyle \left|\phi _{i}-\phi _{i-1}\right|<\epsilon }$

• Logiciel Circé
• Liste de projections cartographiques

• Cours de l'ENSG
• Historique des projections françaises
• Base de données des EPSG
• Projection Lambert dans proj
• Algorithmes de conversion IGN
• « Lambert et la projection cartographique », Physiciens en Alsace, témoins de leurs temps, Université de Strasbourg, p. 34-35 (exposition)

• Portail de la cartographie
• Portail de la géométrie