SERVICE DE GEODESIE

ET NIVELLEMENT NOTES TECHNIQUES

NT/G 71

PROJECTION CARTOGRAPHIQUE
CONIQUE CONFORME DE LAMBERT
S
G
N

2
Algorithmes
7
8
1
0

1ère édition
Janvier 1995

 I N S T I T U T G E O G R A P H I Q U E N A T I O N A L
2-4, AVENUE PASTEUR - 94165 SAINT MANDE CEDEX
 ALGORITHMES NECESS AIRES

A LA

PROJECTION CARTOGRAPHIQUE

CONIQUE CONFORME DE LAMBERT

SOMMAIRE NOMBRE de PAGES

ALG0001 2

ALG0002 3

ALG0003 3

ALG0004 3

ALG0019 3

ALG0021 2

ALG0054 4

APPLICATION 1
ALG0001 1/2

CALCUL DE LA LATITUDE ISOMETRIQUE.

Numéro : ALG0001.

Description :

Calcul de la latitude isométrique sur un ellipsoïde de première excentricité e

au point de latitude ϕ.

Variables :

- paramètres en entrée :

ϕ : latitude.
e : première excentricité de l’ellipsoïde.

- paramètre en sortie :

L : latitude isométrique.

Schéma séquentiel :

E : ϕ , e.
S : L

e
π ϕ 1 − e ⋅ sin ϕ
L = ln( tan(
4
+
2
)⋅(
1 + e ⋅ sin ϕ
)2 )

S
ALG0001 2/2

CALCUL DE LA LATITUDE ISOMETRIQUE.

Jeux d’essai :

ϕ (rad) 0,872 664 626 00 -0,300 000 000 00 0,199 989 033 70

e 0,081 991 889 98 0,081 991 889 98 0,081 991 889 98

L 1,005 526 536 49 -0,302 616 900 63 0,200 000 000 009

Remarque :

On notera L(ϕϕ,e) la valeur de la latitude isométrique sur l’ellipsoïde de

première excentricité e au point de latitude ϕ.
ALG0002 1/3

CALCUL DE LA LATITUDE A PARTIR DE LA LATITUDE ISOMETRIQUE.

Numéro : ALG0002.

Description :

Calcul de la latitude ϕ à partir de la latitude isométrique L.

Variables :

- paramètres en entrée :

L : latitude isométrique.
e : première excentricité de l’ellipsoïde.
ε : tolérance de convergence.

- paramètre en sortie :

ϕ : latitude en radian.
ALG0002 2/3

CALCUL DE LA LATITUDE A PARTIR DE LA LATITUDE ISOMETRIQUE.

Schéma séquentiel :

E : L , e , ε.
S : ϕ.

ϕ0 = 2 ⋅ arctan ( exp ( L ) )−
π
2

i ← 0

i ← i + 1

e
ϕi = 2 ⋅ arctan ( (
1 + e ⋅ sin ϕ i − 1
1 − e ⋅ sin ϕ i − 1
) ⋅ exp (
2
L ) )−
π
2

non
ϕ i − ϕ i− 1 〈ε

oui
ALG0002 3/3

CALCUL DE LA LATITUDE A PARTIR DE LA LATITUDE ISOMETRIQUE.

Schéma séquentiel (suite) :

ϕ = ϕi

Jeux d’essai :

L 1,005 526 536 48 -0,302 616 900 60 0,200 000 000 0

e 0,081 991 889 98 0,081 991 889 98 0,081 991 889 98

ε 1.10-11 1.10-11 1.10-11

ϕ (rad) 0,872 664 626 00 -0,299 999 999 97 0,199 989 033 69

Remarque :

On notera L-1(L,e) la valeur de la latitude à partir de la latitude isométrique

L pour un ellipsoïde de première excentricité e.
ALG0003 1/3

TRAN SFO RMA TIO N DE COO RDO NNE ES

λ , ϕ →
→ X , Y Lambert.

Numéro : ALG0003.

Description :

Transformation de coordonnées géographiques en coordonnées en projection conique

conforme de Lambert.

Variables :

- paramètres en entrée :

λ : longitude par rapport au méridien origine.

ϕ : latitude.
n : exposant de la projection.
c : constante de la projection.
e : première excentricité de l’ellipsoïde.
λc : longitude de l’origine par rapport au méridien origine.
Xs, Ys : coordonnées en projection du pôle.

- paramètres en sortie :

X, Y : coordonnées en projection du point.

Autre algorithme utilisé :

ALG0001 : calcul de la latitude isométrique

Algorithmes dont les résultats sont utilisés en entrée :

ALG0019 : détermination des paramètres de calcul n, c, λc, Xs, Ys en

fonction des paramètres de définition usuels, dans le cas tangent.

ALG0054 : détermination des paramètres de calcul n, c, λc, Xs, Ys en

fonction des paramètres de définition usuels, dans le cas sécant.
ALG0003 2/3

TRAN SFO RMA TIO N DE COO RDO NNE ES

λ , ϕ →
→ X , Y Lambert.

Schéma séquentiel :

E : e , n , c , λc , Xs , Ys , λ , ϕ.
S : X , Y.

E
↓
L = L (ϕ, e ) ALG0001

↓
X = X s + c ⋅ exp(− n ⋅ L )⋅ sin (n (λ − λ c ))
Y = Ys − c ⋅ exp(− n ⋅ L )⋅ cos(n (λ − λ c ))

↓
S

Notation utilisée :


(ϕ ) : latitude isométrique croissante sur l’ellipsoïde

ALG0003 3/3

TRAN SFO RMA TIO N DE COO RDO NNE ES

λ , ϕ →
→ X , Y Lambert.

Jeu d’essai :

e 0,082 483 256 8

n 0,760 405 966

c (m) 11 603 796,976 7

λc (rad) 0,040 792 344 33

Xs (m) 600 000,000 0

Ys (m) 5 657 616,674 0

λ (rad) 0,145 512 099 00

ϕ (rad) 0,872 664 626 00

X (m) 1 029 705,081 8

Y (m) 272 723,851 0

ALG0004 1/3

TRANSFORMATION DE COORDONNEES

X , Y Lambert 
→ λ , ϕ.

Numéro : ALG0004.

Description :

Transformation de coordonnées en projection conique conforme de Lambert, en

coordonnées géographiques.

Variables :

- paramètres en entrée :

X, Y : coordonnées en projection conique conforme de Lambert du

point.
n : exposant de la projection.
c : constante de la projection.
e : première excentricité de l’ellipsoïde.
λc : longitude de l’origine par rapport au méridien origine.
Xs, Ys : coordonnées en projection du pôle.
ε : tolérance de convergence

- paramètres en sortie :

λ : longitude par rapport au méridien origine.

ϕ : latitude.

Autre algorithme utilisé :

ALG0002 : calcul de la latitude à partir de la latitude isométrique.

Algorithmes dont les résultats sont utilisés en entrée:

ALG0019 : détermination des paramètres de calcul n, c, λc, Xs, Ys

dans le cas d’une projection Lambert tangente en fonction
des paramètres de définition usuels.
ALG0054 : détermination des paramètres de calcul n, c, λc, Xs, Ys en
fonction des paramètres de définition usuels, dans le cas
sécant.
ALG0004 2/3

TRANSFORMATION DE COORDONNEES

X , Y Lambert 
→ λ , ϕ.

Schéma séquentiel :

E : n , e , c , λc , Xs , Ys , X , Y, ε.
S : λ , ϕ.

R = ( X − X s )2 + ( Y − Y s )2

X − Xs
γ = arctan ( )
Ys − Y

γ
λ = λ + n
c

1 R
L = − n ⋅
ln
c

ϕ = L (L,e)
-1
ALG0002

Notation utilisée :

L-1(L,e) : latitude à partir de la latitude isométrique L, calculée avec la

tolérance ε.
ALG0004 3/3

TRANSFORMATION DE COORDONNEES

X , Y Lambert 
→ λ , ϕ.

Jeux d’essai :

X (m) 1 029 705,083 0

Y (m) 272 723,849 0

n 0,760 405 966

c (m) 11 603 796,976 7

Xs (m) 600 000,000 0

Ys (m) 5 657 616,674 0

λc (rad) 0,040 792 344 33

e 0,082 483 256 8

ε 1. 10-11

λ (rad) 0,145 512 099 25

ϕ (rad) 0,872 664 625 67

ALG0019 1/3

PARAMETRES DE PROJECTION

Projection Lambert conique conforme dans le cas tangent

Numéro : ALG0019.

Description :

Détermination des paramètres de calcul d’une projection Lambert conique

conforme dans le cas tangent, avec ou sans facteur d'échelle en fonction des
paramètres de définition usuels.

Variables :

- paramètres en entrée :

a : demi-grand axe de l'ellipsoïde.

e : 1 ère excentricité de l'ellipsoïde.
λ0 : longitude origine par rapport au méridien origine.
ϕ0 : latitude origine.
k0 : facteur d’échelle à l’origine.
X0, Y0 : coordonnées en projection du point origine.

- paramètres en sortie :

e : 1 ère excentricité de l'ellipsoïde.

λc : longitude origine par rapport au méridien origine.
n : exposant de la projection.
C : constante de la projection.
Xs, Ys : coordonnées du pôle en projection.

Autres algorithmes utilisés :

ALG0001 : Calcul de la latitude isométrique.

ALG0021 : Calcul de la grande normale.
ALG0019 2/3

PARAMETRES DE PROJECTION

Projection Lambert conique conforme dans le cas tangent

Schéma séquentiel :

E : a , e , λ0 , ϕ0 , k0 , X0 , Y0.
S : e , n , C , λ c , Xs , Ys.

λc = λ0

n = sin ϕ 0

C = k 0 ⋅ N ( ϕ 0 , a , e ) ⋅ cot an ϕ0 ⋅ exp ( n ⋅ L ( ϕ0 ,e ) ) ALG0021/ALG0001

XS = X0

Y S = Y 0 + k 0 ⋅ N ( ϕ 0 , a , e ) ⋅ cot an ϕ0
ALG0021

Notation utilisée :

L(ϕ,e) : latitude isométrique croissante sur l'ellipsoïde de première

excentricité au point de latitude ϕ.
ALG0019 3/3

PARAMETRES DE PROJECTION

Projection Lambert conique conforme dans le cas tangent

Jeux d'essai:

λ0 (rad) 0,181 128 088 00 0,040 792 344 33

ϕ0 (rad) 0,977 384 381 00 0,863 937 980 00

k0 1,000 000 000 0 0,999 877 340 0

X0 (m) 0,000 0 600 000,000 0

Y0 (m) 0,000 0 200 000,000 0

a (m) 6 378 388,000 0 6 378 249,200 0

e 0,081 991 890 0,082 483 256 8

λc (rad) 0,181 128 088 00 0,040 792 344 33

e 0,081 991 890 0,082 483 256 8

n 0,829 037 572 5 0,760 405 965 8

C (m) 11 464 828,219 2 11 603 796,976 0

Xs (m) 0,000 0 600 000,000 0

Ys (m) 4 312 250,971 8 5 657 616,671 2

ALG0021 1/2

CALCUL DE LA GRANDE NORMALE

Numéro : ALG0021.

Description :

Calcul de la grande normale de l’ellipsoïde.

Variables :

- paramètres en entrée :

ϕ : latitude.
a : demi-grand axe de l’ellipsoïde.
e : première excentricité de l’ellipsoïde.

- paramètre en sortie :

N : grande normale.

Schéma séquentiel :

E : ϕ , a , e.
S : N.

a
N = 2 2
1 − e ⋅ sin ϕ

S
ALG0021 2/2

CALCUL DE LA GRANDE NORMALE

Jeux d’essai :

0,977 384 381 00

ϕ(rad)

a(m) 6 378 388,000 0

e 0,081 991 890

N(m) 6 393 174,975 5

Remarque :

On notera N(ϕ,e,a) la valeur de la grande normale d’un ellipsoïde donné (a,e)

en un point de latitude ϕ.
ALG0054 1/4

PARAMETRES DE PROJECTION

Projection Lambert sécante.

Numéro : ALG0054.

Description :

Calcul des constantes d'une projection Lambert conique conforme dans le cas
sécant.

Variables :

- paramètres en entrée :

a : demi-grand axe de l'ellipsoïde.

e : première excentricité de l'ellipsoïde.
λ0 : longitude origine en radian par rapport au méridien origine.
ϕ0 : latitude origine.
ϕ1 : latitude en radian du 1er parallèle automécoïque.
ϕ2 : latitude en radian du 2ème parallèle automécoïque.
X0, Y0 : coordonnées en projection du point origine.

- paramètres en sortie :

e : première excentricité de l'ellipsoïde.

λc : longitude origine par rapport au méridien origine.
n : exposant de la projection.
c : constante de la projection.
Xs, Ys : coordonnées du pôle en projection.

Autres algorithmes utilisés :

ALG0001 : Calcul de la latitude isométrique.

ALG0021 : Calcul de la grande normale.
ALG0054 2/4

PARAMETRES DE PROJECTION

Projection Lambert sécante.

Schéma séquentiel :

E : a , e , λ0 , ϕ0 , ϕ1 , ϕ2 , k0 , X0 , Y0.

S : e , n , c , λ c , Xs , Ys.

λc = λ 0

N ( ϕ 2 , a , e ) ⋅ cos ϕ 2
ln ( ) ALG0021
N ( ϕ 1 , a , e ) ⋅ cos ϕ 1
n = ( )
L (ϕ 1 , e) − L (ϕ 2 , e)
ALG0001

C =
N ( ϕ 1 , a , e ) ⋅ cos ϕ 1
n
⋅ exp ( n ⋅ L (ϕ 1 , e ))
ALG0001 / ALG0021

Xs = X 0

Ys = Y 0
ALG0054 3/4

PARAMETRES DE PROJECTION

Projection Lambert sécante.

Schéma séquentiel (suite) :

oui π non
ϕ0 = 2

XS X0 Xs X0
= =
YS Y0
= Ys
=
Y0
+
c
⋅
exp ( n
− ⋅ L ( 0 , e) )
ϕ

Notation utilisée :

L(ϕ,e) : latitude isométrique croissante sur l'ellipsoïde

ALG0054 4/4

PARAMETRES DE PROJECTION

Projection Lambert sécante.

Jeux d'essai:

λ0 (rad) 0,000 000 000 00 0,076 235 545 39

ϕ0 (rad) 0,000 000 000 00 1,570 796 327 00

X0 (m) 0,000 0 150 000,000 0

Y0 (m) 0,000 0 5 400 000,000 0

ϕ1 (rad) -0,575 958 653 00 0,869 755 744 00

ϕ2 (rad) -0,785 398 163 00 0,893 026 801 00

a (m) 6 378 388,000 0 6 378 388,000 0

e 0,081 991 890 0,081 991 890

λc (rad) 0,000 000 000 00 0,076 235 545 39

e 0,081 991 890 0,081 991 890

n -0,630 496 330 0 0,771 642 186 7

c (m) -12 453 174,179 5 11 565 915,829 4

Xs (m) 0,000 0 150 000,000 0

Ys (m) -12 453 174,179 5 5 400 000,000 0

 Lambert France 1/1

CONSTANTES DE PROJECTION

Projection Lambert France.

Description :

Valeurs des constantes n , c , Xs , Ys , λ0 des 5 projections de types Lambert

conique conforme en usage en France et du Lambert II étendu.

Valeurs des constantes Lambert France :

Lambert I Lambert II Lambert III Lambert IV Lambert-93

n 0,760 405 965 6 0,728 968 627 4 0,695 912 796 6 0,671 267 932 2 0,725 607 765 0

c (m) 11 603 796,98 11 745 793,39 11 947 992,52 12 136 281,99 11 754 255,426

Xs (m) 600 000,0 600 000,0 600 000,0 234,358 700 000,0

Ys (m) 5 657 616,674 6 199 695,768 6 791 905,085 7 239 161,542 12 655 612,050

λ0 = 2°20’14,025” E par rapport au méridien de Greenwich

= 0 grades par rapport au méridien de Paris.

e = 0,082 483 256 76 (première excentricité de l’ellipsoïde Clarke 1880

français).

Les constantes en usage pour le Lambert II étendu sont celles du

Lambert II avec Ys = 8 199 695,768.