Lambert-93

Lambert-93

Projection conique conforme de Lambert

Page d'aide sur l'homonymie Pour les articles homonymes, voir Lambert.
Page d'aide sur l'homonymie Ne doit pas être confondu avec Projection azimutale équivalente de Lambert.

La projection conique conforme de Lambert, ou plus simplement, la projection de Lambert est l'une des projections cartographiques présentées par le mathématicien mulhousien Johann Heinrich Lambert en 1772[1].

Dans ce système de projection conforme, les méridiens sont des droites concourantes, et les parallèles des arcs de cercle centrés sur le point de convergence des méridiens.

Le système a été adopté par l'artillerie française au cours de la première guerre mondiale pour les cartes à grande échelle (1/20 000 et au-dessus), une carte conforme étant nécessaire pour la préparation des tirs, ce que ne permettait pas la projection de Bonne alors en usage.

Elle est depuis cette époque la projection officielle utilisée pour représenter la France métropolitaine, avec différents paramètres successifs suivant les époques. C'est aussi la projection officielle en Belgique et en Estonie, ainsi que pour les cartes couvrant toute l'Europe à des échelles inférieures ou égales au 1:500 000[2].

Sommaire

Définition

Projection conique.jpg

Cette projection est une projection conique conforme (qui conserve les angles).

Le sommet du cône appartient à l'axe des pôles et est tangent à un ellipsoïde de référence en un point défini par un méridien de référence et un parallèle de référence de latitude φ0 qui est aussi l'angle au sommet du cône. Le cône peut également être sécant à l'ellipsoïde au niveau de deux parallèles, dits alors parallèles automécoïques φ1 et φ2. La définition d'une projection de Lambert peut ainsi se faire soit par un parallèle tangent et un facteur d'échelle, soit par deux parallèles sécants et isométriques.

Le cône est ensuite développé sur un plan, sans déformation.

Dans cette projection on retrouve :

  • les parallèles (latitude constante) sont des cercles concentriques autour du point P, projection du pôle Nord et sommet du cône.
  • les méridiens (longitude constante) sont des droites concourantes en P.
  • l'axe des ordonnées est la projection du méridien de référence.
  • l'origine se trouve au point de référence.
  • le cercle, projection du point de référence, est appelé isomètre ou isomètre de référence. En effet, c'est selon ce cercle que l'on définit l'échelle de la carte.
  • les angles sont conservés.

La déformation linéaire est relativement importante lorsque l'on s'éloigne du parallèle d'origine, il est donc courant de définir des zones différant par leur angle φ0.

Enfin cette projection étant utilisée en cartographie, pour éviter d'avoir des abscisses et des ordonnées négatives, on effectue une translation en affectant des coordonnées positives au point de référence.

Projections officielles en France métropolitaine

La projection Lambert93 (projection officielle pour les cartes de France métropolitaine depuis le décret du 26 décembre 2000) est la projection liée au système géodésique RGF93. Elle utilise deux parallèles sécants : 44°N et 49°N, le méridien de référence 3°E (Greenwich), le parallèle d'origine (qui n'est pas de référence puisqu'on utilise des parallèles sécants) 46°30', et les coordonnées d'origine (700 000 m, 6 600 000 m).

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, et d'autre part compatible avec le WGS84 utilisé notamment par le système GPS de positionnement par satellite (le passage de RGF93 à WGS84 se fait par une translation métrique, et des rotations et mises à l'échelle de l'ordre de 10-9, les rendant virtuellement identiques pour la plupart des applications pratiques). Ces trois systèmes ont de plus en commun l'ellipsoïde de référence IAG GRS80.

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

En pratique, cette projection est encore 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 entérine 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. Elles ont en commun avec le Lambert93 le système géodésique RGF93 et le méridien de référence 3°E (Greenwich).

Projection φ0 φ1 φ2 X0 Y0
CC42 42° 41.25° 42.75° 1 700 000 m 1 200 000 m
CC43 43° 42.25° 43.75° 1 700 000 m 2 200 000 m
CC44 44° 43.25° 44.75° 1 700 000 m 3 200 000 m
CC45 45° 44.25° 45.75° 1 700 000 m 4 200 000 m
CC46 46° 45.25° 46.75° 1 700 000 m 5 200 000 m
CC47 47° 46.25° 47.75° 1 700 000 m 6 200 000 m
CC48 48° 47.25° 48.75° 1 700 000 m 7 200 000 m
CC49 49° 48.25° 49.75° 1 700 000 m 8 200 000 m
CC50 50° 49.25° 50.75° 1 700 000 m 9 200 000 m

Icône pour souligner l'importance du texte Dans l'attente d'une définition officielle de ces projections, les Y0 ci-dessus sont donnés sous toutes réserves et selon cette publication du CERTU, des valeurs de 42 200 000 m à 50 200 000 m ayant également été publiées!

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).

A 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.

Projection φ0 φ1 φ2 Translation X Translation Y
Lambert I (Nord) 55 grades 48°35’54,682” 50°23’45,282” 600 000 m 200 000 m
Lambert II (Centre) 52 grades 45°53’56,108” 47°41’45,652” 600 000 m 200 000 m
Lambert III (Sud) 49 grades 43°11’57,449” 44°59’45,938” 600 000 m 200 000 m
Lambert IV (Corse) 46,85 grades 41°33’37,396” 42°46’03,588” 234,358 m 185 861,369 m
Lambert I Carto (Nord) 55 grades 48°35’54,682” 50°23’45,282” 600 000 m 1 200 000 m
Lambert II Carto (Centre) 52 grades 45°53’56,108” 47°41’45,652” 600 000 m 2 200 000 m
Lambert III Carto (Sud) 49 grades 43°11’57,449” 44°59’45,938” 600 000 m 3 200 000 m
Lambert IV Carto (Corse) 46,85 grades 41°33’37,396” 42°46’03,588” 234,358 m 4 185 861,369 m
Lambert Grand Champ 47° 45° 49° 600 000 m 600 000 m

On signalera que avec les coordonnées Lambert, Le Lambert Zone "Carto" est utilisé par l'IGN depuis 1972, alors que le Lambert zone "standard" reste utilisé par les services du Cadastre ou autres utilisateurs locaux (géomètres notamment) : ainsi, en Lambert 3, on pourra trouver Y=3 256 315 sur les cartes IGN et Y=256 315 au cadastre. Cette notation dite "standard" est une survivance d'anciennes définitions IGN, et ne devrait plus être utilisée depuis 1972.

Pour les représentations devant couvrir la France entière en une projection, la projection utilisée est la "Lambert 2 étendue", qui reprend tous les paramètres de la Lambert 2 Carto. C'est à ce jour la projection généralement utilisée dans les SIG français. À terme la projection Lambert93 remplacera le lambert 2 étendu, qui ne présente plus de compatibilité et de précision suffisantes.

Formules mathématiques

Pour l'application numérique, on rappelle ici les paramètres pour la projection Lambert 93 :

  • a = 6 378 137 m
  • b = a.(1 - f)
  • f = 1 / 298.257222101
  • φ0 = 46°30
  • φ1 = 44°
  • φ2 = 49°
  • λ0 = 3°
  • X0 = 700 000 m
  • Y0 = 6 600 000 m

Des coordonnées géographiques (lat, lon) aux cartographiques (X, Y)

Partant d'un point (φ, λ) du globe supposé être une ellipsoïde de demi-grand axe a et de demi-petit axe b, calculons ses coordonnées (X, Y) sur la carte. Pour cela, nous allons passer par les coordonnées (ρ, θ) du point projeté sur le cône. L'axe des X est croissant vers l'Est et l'axe des Y croissant vers le Nord.

\left\{\begin{matrix} X& = &\rho \sin(\theta) + X_0\\ Y& =& Y_0 + \rho_0 - \rho \cos(\theta)\end{matrix}\right. où \left\{\begin{matrix}\theta & =& n(\lambda - \lambda_0)\\\rho & = &\rho(\phi)\end{matrix}\right.


Démonstration

Puisque c'est une projection conforme pour un point M à la surface de la terre :

\frac{ \frac{\partial X}{\partial\lambda}}{\frac{\partial Y}{\partial\phi}} = \frac{\left| \frac{\partial \vec M}{\partial \lambda} \right|}{\left| \frac{\partial \vec M}{\partial \phi} \right|}\ \left(1\right)

On peut donc calculer ρ le long du méridien d'origine (λ = λ0)

  • \frac{\partial X}{\partial \lambda} = \rho\;n\ \left(2\right)
  • \frac{\partial Y}{\partial \phi} = - \frac{\partial \rho}{\partial \phi}\ \left(3\right)

Le méridien est une ellipse d'équation :

\left\{\begin{matrix} r & = a\cos(t) \\ z & = b\sin(t)\end{matrix}\right.\forall t \in [0;2\pi]

La latitude étant l'angle que fait le méridien du lieu avec le plan équatorial, on trouve :

\tan(\phi) = -\frac{\left.\frac{\partial \vec M}{\partial t} \right\rfloor_r}{\left.\frac{\partial \vec M}{\partial t} \right\rfloor_z} = \frac{a}{b}\tan(t)\ \left(4\right)

\left\|\frac{\partial \vec M}{\partial \lambda}\right\| = a\cos(t)\ \left(5\right)

\left\|\frac{\partial \vec M}{\partial \phi}\right\| = \left\|\frac{\partial \vec M}{\partial t}\right\| \left|\frac{\partial t}{\partial \phi}\right|

\left\|\frac{\partial \vec M}{\partial t}\right\| = \sqrt{a^2 \sin(t)^2 + b^2 \cos(t)^2} = \left|b\cos(t)\right|\sqrt{ \frac{a^2}{b^2}tan(t)^2 + 1} = \left|b\cos(t)\right|\sqrt{\tan(\phi)^2 + 1} = \left|\frac{b\cos(t)}{\cos(\phi)}\right| = \frac{b\cos(t)}{\cos(\phi)}

On peut dériver (4) par rapport à φ :

\frac{1}{\cos(\phi)^2} = \frac{a}{b}\frac{\partial t}{\partial \phi}\frac{1}{\cos(t)^2}

Ce qui donne en regroupant les trois dernières équations :

\left\|\frac{\partial \vec M}{\partial \phi}\right\| = \frac{b\cos(t)}{\cos(\phi)}\frac{b}{a}\frac{\cos(t)^2}{\cos(\phi)^2}\ \left(6\right)

Et en regroupant (1), (2), (3), (5) et (6) :

\frac{\partial \rho}{\partial \phi} = - \rho\;n \frac{b\cos(t)}{\cos(\phi)}\frac{b}{a}\frac{\cos(t)^2}{\cos(\phi)^2} \frac{1}{a\cos(t)} = -\rho\;n\frac{b^2}{a^2}\frac{\cos(t)^2}{\cos(\phi)^3}\ \left(7\right)

\frac{b^2}{a^2}\frac{\cos(t)^2}{\cos(\phi)^3} = \frac{b^2}{a^2}\frac{1}{(1 + \tan(t)^2)\cos(\phi)^3} = \frac{b^2}{a^2}\frac{1}{(1 + \frac{b^2}{a^2}\tan(\phi)^2)\cos(\phi)^3} = \frac{b^2}{(a^2\cos(\phi)^2 + b^2\sin(\phi)^2)\cos(\phi)}\,

\frac{b^2}{a^2}\frac{\cos(t)^2}{\cos(\phi)^3} = \frac{1}{\cos(\phi)} + \frac{b^2 - b^2\sin(\phi)^2 - a^2\cos(\phi)^2}{(a^2\cos(\phi)^2 + b^2\sin(\phi)^2)\cos(\phi)} = \frac{1}{\cos(\phi)} + \frac{(b^2-a^2)\cos(\phi)}{a^2 - (a^2 - b^2)\sin(\phi)^2}\,

Or e\,, l'excentricité d'une ellipse vérifie : e^2 = \frac{a^2 - b^2}{a^2}. Ce qui nous donne :

\frac{b^2}{a^2}\frac{\cos(t)^2}{\cos(\phi)^3} = \frac{1}{\cos(\phi)} - \frac{e^2\cos(\phi)}{1 - e^2\sin(\phi)^2} = \frac{1}{\sin(\phi + \frac{\pi}{2})} - \frac{e}{2}\frac{2e\cos(\phi)}{(1 + e\sin(\phi))(1 - e\sin(\phi))} = \frac{1}{2\sin(\frac{\phi}{2} + \frac{\pi}{4})\cos(\frac{\phi}{2} + \frac{\pi}{4})} - \frac{e}{2}\frac{e\cos(\phi)(1 + e\sin(\phi)) + e\cos(\phi)(1 - e\sin(\phi))}{(1 + e\sin(\phi))(1 - e\sin(\phi))}

Et finalement en utilisant les dérivées de la fonction \ln\, :

\frac{b^2}{a^2}\frac{\cos(t)^2}{\cos(\phi)^3} = \frac{\partial}{\partial \phi}\left(\ln(\tan(\frac{\phi}{2} + \frac{\pi}{4})) - \frac{e}{2}\ln(\frac{1 + e\sin(\phi)}{1 - e\sin(\phi)})\right)

On peut donc réécrire l'équation (7) en :

\frac{\partial \ln(\rho)}{\partial \phi} = -n\frac{\partial}{\partial \phi}\left(\ln(\tan(\frac{\phi}{2} + \frac{\pi}{4})) - \frac{e}{2}\ln(\frac{1 + e\sin(\phi)}{1 - e\sin(\phi)})\right)

\rho(\phi) = \rho(0)\left(\cot(\frac{\phi}{2} + \frac{\pi}{4})(\frac{1 + e\sin(\phi)}{1 - e\sin(\phi)})^{\frac{e}{2}}\right)^n

Nous allons travailler ici avec deux parallèles sécants de référence. Sur ces parallèles, appelés automécoïques, de latitude \phi_1\, et \phi_2\,, l'échelle est conservée : n\;\rho_i = a\cos(t_i) = \frac{a\cos(\phi_i)}{\sqrt{1 - e^2\sin(\phi_i)^2}}\, On peut donc calculer n et ρ(0) :

n = \left.\frac{\ln\left(\frac{\cos(\phi_2)}{\cos(\phi_1)}\right)+\frac{1}{2}\ln\left(\frac{1-e^2\sin(\phi_1)^2}{1-e^2\sin(\phi_2)^2}\right)}{\ln\left(\frac{\tan(\frac{\phi_1}{2}+\frac{\pi}{4})(1 - e\sin(\phi_1))^{\frac{e}{2}}(1 + e\sin(\phi_2))^{\frac{e}{2}} }{\tan(\frac{\phi_2}{2}+\frac{\pi}{4})(1 + e\sin(\phi_1))^{\frac{e}{2}}(1 - e\sin(\phi_2))^{\frac{e}{2}} }\right)}\right.

\rho(0) = \frac{a\cos(\phi_1)}{n\sqrt{1 - e^2\sin(\phi_1)^2}}\left(\tan(\frac{\phi_1}{2} + \frac{\pi}{4})\left(\frac{1-e\sin(\phi_1)}{1+e\sin(\phi_1)}\right)^{\frac{e}{2}}\right)^n

Des coordonnées cartographiques (X,Y) aux géographiques (lat, lon)

De même que ci-dessus nous allons passer par les coordonnées (ρ, θ) du point projeté sur le cone :

  • \rho = \sqrt{(X - X_0)^2 + (Y_0 - Y + \rho_0)^2}
  • \theta = 2\tan^{-1}(\frac{X - X_0}{Y_0 - Y + \rho_0 + \rho})

Comme on a vu plus haut, on a alors :

  • \lambda = \frac{\theta}{n} + \lambda_0\,
  • \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 :

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

Mais le calcul de F − 1 n'est pas trivial. Pour calculer φ avec une précision arbitraire, puisque F est continue sur son intervalle de définition, on peut utiliser un algorithme dichotomique. On peut aussi, pour plus de rapidité, utiliser l'algorithme proposé par l'IGN :

  • \phi_0 = 2\tan^{-1}\left(\left(\frac{\rho(0)}{\rho}\right)^\frac{1}{n}\right) - \frac{\pi}{2}
  • \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 \left|\phi_i - \phi_{i-1}\right| < \epsilon

Bibliographie

  1. ↑ Karen Mulcahy, « Cylindrical Projections », City University of New York. Consulté le 30 mars 2007
  2. ↑ (en) [pdf] Institute for Environment and Sustainability, European Communities,Map Projections for Europe, 2001, page 43.

Voir aussi

Ce document provient de « Projection conique conforme de Lambert#Projections officielles en France m.C3.A9tropolitaine ».

Wikimedia Foundation. 2010.

Contenu soumis à la licence CC-BY-SA. Source : Article Lambert-93 de Wikipédia en français (auteurs)

Игры ⚽ Поможем сделать НИР

Regardez d'autres dictionnaires:

  • Lambert — may refer to*Lambert of Maastricht, bishop, saint, and martyr *Lambert of St Bertin or Lambert of St Omer, medieval encyclopedist *Lambert Mieszkowic, son of Mieszko I of Poland *Lambert McKenna, Irish scholar, Editor and Lexicographer. (1870… …   Wikipedia

  • Lambert — ist ein Vor und Familienname, siehe Lambert (Name) Orte in den Vereinigten Staaten: Lambert (Arkansas), Stadt in Arkansas Lambert (Mississippi), Stadt in Mississippi Lambert (Missouri), Stadt in Missouri Lambert (Montana), Stadt in Montana… …   Deutsch Wikipedia

  • LAMBERT (J. H.) — LAMBERT JOHANN HEINRICH (1728 1777) Mathématicien, astronome, physicien et philosophe suisse et allemand d’ascendance française. Né à Mulhouse, qui faisait alors partie de la Suisse, Johann Heinrich Lambert, dès l’âge de douze ans, quitte l’école …   Encyclopédie Universelle

  • LAMBERT (M.) — LAMBERT MICHEL (1610 env. 1696) Musicien français, compositeur, luthiste et chanteur célèbre. Le rôle de Lambert fut déterminant dans l’élaboration du récitatif à la française. La méthode de chant de cet excellent pédagogue exerça une influence… …   Encyclopédie Universelle

  • lambert — LAMBÉRT, lamberţi, s.m. (fiz.) Unitate (tolerată) de măsură a luminanţei. [abr.: L] – Din fr. lambert, germ. Lambert. Trimis de LauraGellner, 16.05.2004. Sursa: DEX 98  lambért s. m., pl. lambérţi; simb. la …   Dicționar Român

  • Lambert — Lambert, MO U.S. village in Missouri Population (2000): 49 Housing Units (2000): 17 Land area (2000): 0.052007 sq. miles (0.134698 sq. km) Water area (2000): 0.000000 sq. miles (0.000000 sq. km) Total area (2000): 0.052007 sq. miles (0.134698 sq …   StarDict's U.S. Gazetteer Places

  • LAMBERT (J.) — LAMBERT JOHN (1619 1684) Issu de la gentry du Yorkshire, John Lambert reçoit une formation universitaire et juridique. Au début de la guerre civile, il est capitaine sous les ordres de Thomas Fairfax, mais, dès 1643, il est à la tête de son… …   Encyclopédie Universelle

  • lambert — (saint ) (sin lan bêr) C est aujourd hui la Saint Lambert, Qui quitte sa place la perd ; locution qui vient, dit on, de ce que la ville de Troyes, prise le jour de la Saint Lambert en 1590, fut reprise presque aussitôt après. Il est aujourd hui… …   Dictionnaire de la Langue Française d'Émile Littré

  • Lambert, MO — U.S. village in Missouri Population (2000): 49 Housing Units (2000): 17 Land area (2000): 0.052007 sq. miles (0.134698 sq. km) Water area (2000): 0.000000 sq. miles (0.000000 sq. km) Total area (2000): 0.052007 sq. miles (0.134698 sq. km) FIPS… …   StarDict's U.S. Gazetteer Places

  • Lambert, MS — U.S. town in Mississippi Population (2000): 1967 Housing Units (2000): 676 Land area (2000): 0.813085 sq. miles (2.105881 sq. km) Water area (2000): 0.000000 sq. miles (0.000000 sq. km) Total area (2000): 0.813085 sq. miles (2.105881 sq. km) FIPS …   StarDict's U.S. Gazetteer Places

  • Lambert, OK — U.S. town in Oklahoma Population (2000): 9 Housing Units (2000): 6 Land area (2000): 0.114863 sq. miles (0.297494 sq. km) Water area (2000): 0.000000 sq. miles (0.000000 sq. km) Total area (2000): 0.114863 sq. miles (0.297494 sq. km) FIPS code:… …   StarDict's U.S. Gazetteer Places

Share the article and excerpts

Direct link
… Do a right-click on the link above
and select “Copy Link”