0% ont trouvé ce document utile (0 vote)
107 vues23 pages

NTG 71

Ce document décrit plusieurs algorithmes nécessaires pour effectuer des transformations entre coordonnées géographiques et coordonnées projetées en projection conique conforme de Lambert. Les algorithmes calculent des valeurs telles que la latitude isométrique, les coordonnées projetées à partir des coordonnées géographiques, et inversement.

Transféré par

ilyesin
Copyright
© © All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats PDF, TXT ou lisez en ligne sur Scribd
0% ont trouvé ce document utile (0 vote)
107 vues23 pages

NTG 71

Ce document décrit plusieurs algorithmes nécessaires pour effectuer des transformations entre coordonnées géographiques et coordonnées projetées en projection conique conforme de Lambert. Les algorithmes calculent des valeurs telles que la latitude isométrique, les coordonnées projetées à partir des coordonnées géographiques, et inversement.

Transféré par

ilyesin
Copyright
© © All Rights Reserved
Nous prenons très au sérieux les droits relatifs au contenu. Si vous pensez qu’il s’agit de votre contenu, signalez une atteinte au droit d’auteur ici.
Formats disponibles
Téléchargez aux formats PDF, TXT ou lisez en ligne sur Scribd
Vous êtes sur la page 1/ 23

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

1re 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.

Numro : ALG0001.

Description :

Calcul de la latitude isomtrique sur un ellipsode de premire excentricit e


au point de latitude .

Variables :

- paramtres en entre :

: latitude.
e : premire excentricit de lellipsode.

- paramtre en sortie :

L : latitude isomtrique.

Schma squentiel :

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 dessai :

(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 isomtrique sur lellipsode de


premire excentricit e au point de latitude .
ALG0002 1/3

CALCUL DE LA LATITUDE A PARTIR DE LA LATITUDE ISOMETRIQUE.

Numro : ALG0002.

Description :

Calcul de la latitude partir de la latitude isomtrique L.

Variables :

- paramtres en entre :

L : latitude isomtrique.
e : premire excentricit de lellipsode.
: tolrance de convergence.

- paramtre en sortie :

: latitude en radian.
ALG0002 2/3

CALCUL DE LA LATITUDE A PARTIR DE LA LATITUDE ISOMETRIQUE.

Schma squentiel :

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.

Schma squentiel (suite) :

= i

Jeux dessai :

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 isomtrique


L pour un ellipsode de premire excentricit e.
ALG0003 1/3

TRAN SFO RMA TIO N DE COO RDO NNE ES

,
X , Y Lambert.

Numro : ALG0003.

Description :

Transformation de coordonnes gographiques en coordonnes en projection conique


conforme de Lambert.

Variables :

- paramtres en entre :

: longitude par rapport au mridien origine.


: latitude.
n : exposant de la projection.
c : constante de la projection.
e : premire excentricit de lellipsode.
c : longitude de lorigine par rapport au mridien origine.
Xs, Ys : coordonnes en projection du ple.

- paramtres en sortie :

X, Y : coordonnes en projection du point.

Autre algorithme utilis :

ALG0001 : calcul de la latitude isomtrique

Algorithmes dont les rsultats sont utiliss en entre :

ALG0019 : dtermination des paramtres de calcul n, c, c, Xs, Ys en


fonction des paramtres de dfinition usuels, dans le cas tangent.

ALG0054 : dtermination des paramtres de calcul n, c, c, Xs, Ys en


fonction des paramtres de dfinition usuels, dans le cas scant.
ALG0003 2/3

TRAN SFO RMA TIO N DE COO RDO NNE ES

,
X , Y Lambert.

Schma squentiel :

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

Notation utilise :


( ) : latitude isomtrique croissante sur lellipsode



ALG0003 3/3

TRAN SFO RMA TIO N DE COO RDO NNE ES

,
X , Y Lambert.

Jeu dessai :

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

Numro : ALG0004.

Description :

Transformation de coordonnes en projection conique conforme de Lambert, en


coordonnes gographiques.

Variables :

- paramtres en entre :

X, Y : coordonnes en projection conique conforme de Lambert du


point.
n : exposant de la projection.
c : constante de la projection.
e : premire excentricit de lellipsode.
c : longitude de lorigine par rapport au mridien origine.
Xs, Ys : coordonnes en projection du ple.
: tolrance de convergence

- paramtres en sortie :

: longitude par rapport au mridien origine.


: latitude.

Autre algorithme utilis :

ALG0002 : calcul de la latitude partir de la latitude isomtrique.

Algorithmes dont les rsultats sont utiliss en entre:

ALG0019 : dtermination des paramtres de calcul n, c, c, Xs, Ys


dans le cas dune projection Lambert tangente en fonction
des paramtres de dfinition usuels.
ALG0054 : dtermination des paramtres de calcul n, c, c, Xs, Ys en
fonction des paramtres de dfinition usuels, dans le cas
scant.
ALG0004 2/3

TRANSFORMATION DE COORDONNEES

X , Y Lambert
, .

Schma squentiel :

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 utilise :

L-1(L,e) : latitude partir de la latitude isomtrique L, calcule avec la


tolrance .
ALG0004 3/3

TRANSFORMATION DE COORDONNEES

X , Y Lambert
, .

Jeux dessai :

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

Numro : ALG0019.

Description :

Dtermination des paramtres de calcul dune projection Lambert conique


conforme dans le cas tangent, avec ou sans facteur d'chelle en fonction des
paramtres de dfinition usuels.

Variables :

- paramtres en entre :

a : demi-grand axe de l'ellipsode.


e : 1 re excentricit de l'ellipsode.
0 : longitude origine par rapport au mridien origine.
0 : latitude origine.
k0 : facteur dchelle lorigine.
X0, Y0 : coordonnes en projection du point origine.

- paramtres en sortie :

e : 1 re excentricit de l'ellipsode.
c : longitude origine par rapport au mridien origine.
n : exposant de la projection.
C : constante de la projection.
Xs, Ys : coordonnes du ple en projection.

Autres algorithmes utiliss :

ALG0001 : Calcul de la latitude isomtrique.


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

PARAMETRES DE PROJECTION

Projection Lambert conique conforme dans le cas tangent

Schma squentiel :

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 utilise :

L(,e) : latitude isomtrique croissante sur l'ellipsode de premire


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

Numro : ALG0021.

Description :

Calcul de la grande normale de lellipsode.

Variables :

- paramtres en entre :

: latitude.
a : demi-grand axe de lellipsode.
e : premire excentricit de lellipsode.

- paramtre en sortie :

N : grande normale.

Schma squentiel :

E : , a , e.
S : N.

a
N = 2 2
1 e sin

S
ALG0021 2/2

CALCUL DE LA GRANDE NORMALE

Jeux dessai :

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 dun ellipsode donn (a,e)


en un point de latitude .
ALG0054 1/4

PARAMETRES DE PROJECTION

Projection Lambert scante.

Numro : ALG0054.

Description :

Calcul des constantes d'une projection Lambert conique conforme dans le cas
scant.

Variables :

- paramtres en entre :

a : demi-grand axe de l'ellipsode.


e : premire excentricit de l'ellipsode.
0 : longitude origine en radian par rapport au mridien origine.
0 : latitude origine.
1 : latitude en radian du 1er parallle automcoque.
2 : latitude en radian du 2me parallle automcoque.
X0, Y0 : coordonnes en projection du point origine.

- paramtres en sortie :

e : premire excentricit de l'ellipsode.


c : longitude origine par rapport au mridien origine.
n : exposant de la projection.
c : constante de la projection.
Xs, Ys : coordonnes du ple en projection.

Autres algorithmes utiliss :

ALG0001 : Calcul de la latitude isomtrique.


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

PARAMETRES DE PROJECTION

Projection Lambert scante.

Schma squentiel :

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

Schma squentiel (suite) :

oui non
0 = 2

XS X0 Xs X0
= =
YS Y0
= Ys
=
Y0
+
c

exp ( n
L ( 0 , e) )

Notation utilise :

L(,e) : latitude isomtrique croissante sur l'ellipsode


ALG0054 4/4

PARAMETRES DE PROJECTION

Projection Lambert scante.

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 = 22014,025 E par rapport au mridien de Greenwich


= 0 grades par rapport au mridien de Paris.

e = 0,082 483 256 76 (premire excentricit de lellipsode Clarke 1880


franais).

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


Lambert II avec Ys = 8 199 695,768.

Vous aimerez peut-être aussi