Comment

peut-on effectuer
un changement de système géodésique ?

Processus de changement de système

On entend par changement de système géodésique la transformation qui permet d'exprimer

les coordonnées cartésiennes, géographique ou planes d'un point dans un autre système
géodésique. Il peut être détaillé selon ce modèle :

modèle
système A système B
de transformation

similitude 3D
cartésiennes ← à 7 paramètres → cartésiennes

↕ ↕
formules
géographiques ← de Molodensky → géographiques

↕ ↕
transformation
planes ← polynomiale → planes

La transformation la plus utilisée est la similitude 3D car elle présente l'avantage de pouvoir
être utilisée dans les deux sens : les mêmes paramètres servent à transformer des coordonnées
exprimées dans le système A vers le système B ou du système B vers le système A.
Les formules de Molodensky sont des développements limités dont l'ordre influe sur la
précision finale : le passage inverse nécessite l'application de formules différentes.
La transformation polynomiale ne s'applique que sur des zones restreintes (pour conserver
une précision comparable à celle obtenue par l'emploi d'une similitude).

IGN • Service de Géodésie et de Métrologie • 73 avenue de Paris, 94165 Saint-Mandé Cedex • 01 43 98 83 25 • [email protected] 1/7
Similitude 3D à sept paramètres (également appelée transformation de Helmert)

Les formules permettant de passer du système A au système B font intervenir trois

translations (TX, TY, TZ) en mètres, trois rotations exprimées selon la convention de
l’IERS 1 (RX, RY, RZ) en radians et un facteur d’échelle (Δ) en ppm. Compte tenu des
faibles rotations utilisées (quelques secondes d’arc), les formules peuvent s’écrire :

XB XA TX XA 0 -RZ RY XA
YB = YA + TY +Δ YA + RZ 0 -RX YA
ZB ZA TZ ZA -RY RX 0 ZA

Transformation standard NTFf↔ WGS84

A partir des éléments de géodésie spatiale (Doppler) des années 1970 et 1980, il a pu être
établi pour la France un modèle entre la NTF et le système WGS84. Aucune rotation ou
défaut d'échelle globaux n'ont pu être mis en évidence. C'est donc un modèle ainsi réduit à
trois paramètres (translations) qui a été retenu et largement utilisé.

TX – 168
TY = – 60
TZ + 320

Les valeurs de ces paramètres correspondent à la valeur officielle définie par l'IGN pour la
France dans le document Nouveaux systèmes géodésiques utilisables en France (WGS84,
ED87) (RT/G n°14).
A l'époque de sa détermination, la précision de ce modèle était estimée à environ 2 mètres.
Les déterminations effectuées au travers de l'établissement du RBF (environ 1000 points dans
les systèmes RGF93 et NTF) confirment ces valeurs à 20 cm près, ainsi que la précision
annoncée (soit 2 mètres, avec des valeurs extrêmes inférieures à 5 mètres).
Cette transformation standard a été adoptée dans la quasi totalité des traitements jusqu'à la
mise à disposition de la transformation par grille de paramètres (voir page 4).

1
International Earth Rotation and Reference Systems Service
2/7 IGN • Service de Géodésie et de Métrologie • 73 avenue de Paris, 94165 Saint-Mandé Cedex • 01 43 98 83 25 • [email protected]
Transformation standard NTF↔ ED50

TX – 84
TY = + 37
TZ + 437

Il existe plusieurs jeux de paramètres de transformation entre NTF et ED50. Celui-ci

correspond à la valeur officielle défini par l'IGN pour la France dans le document
Définition des systèmes géodésiques utilisés en France (NTF, ED50, WGS72) (RT/G n°7).

Transformation standard ED50↔ WGS84

TX – 84
TY = – 97
TZ – 117

Cette transformation est issue de la composition des transformations ED50 vers WGS72 et
WGS72 vers WGS84. Les valeurs de ces paramètres correspondent à la valeur officielle
définie par l'IGN pour la France dans les documents Nouveaux systèmes géodésiques
utilisables en France (WGS84, ED87) et Définition des systèmes géodésiques utilisés en
France (NTF, ED50, WGS72) (RT/G n°14 et n°7f).

IGN • Service de Géodésie et de Métrologie • 73 avenue de Paris, 94165 Saint-Mandé Cedex • 01 43 98 83 25 • [email protected] 3/7
Transformations par grille de paramètres NTF↔ RGF93

Dans le contexte de l'introduction d'un nouveau système géodésique de référence pour la

France métropolitaine, l'Institut national de l’information géographique et forestière a
développé un nouveau processus de transformation de coordonnées. Il facilite le passage de
l'ancien système (NTF – Nouvelle Triangulation de la France) au nouveau (RGF93).
Le principe de ce processus de transformation est l'interpolation de paramètres
tridimensionnels de translation entre systèmes dans un semis de points régulièrement répartis
(appelé grille).
Ces points constituent la grille de paramètres GR3DF97A au pas régulier de 0.1° en
longitude et latitude.
La grille de paramètres de transformation de coordonnées NTF ↔ RGF93 est fournie sous la
forme du fichier texte (ASCII) :

GR3D 002024 024 20370201

GR3D1 -5.5000 10.0000 41.0000 52.0000 .1000 .1000
GR3D2 INTERPOLATION BILINEAIRE
GR3D3 PREC CM 01:5 02:10 03:20 04:50 99>100
-5.500000000 41.000000000 -165.027 -67.100 315.813 99 -0158
-5.500000000 41.100000000 -165.169 -66.948 316.007 99 -0157
-5.500000000 41.200000000 -165.312 -66.796 316.200 99 -0157
.....
-4.900000000 47.600000000 -170.859 -58.180 322.743 99 -0321
-4.900000000 47.700000000 -170.815 -58.186 322.693 03 L0320
-4.900000000 47.800000000 -170.773 -58.204 322.644 03 L0320
-4.900000000 47.900000000 -170.732 -58.233 322.596 02 0319
…

4/7 IGN • Service de Géodésie et de Métrologie • 73 avenue de Paris, 94165 Saint-Mandé Cedex • 01 43 98 83 25 • [email protected]
Entête
Champ Définition Valeurs Définition
GR3D codes IGN 002024 002 (NTF) vers 024 (RGF93)
024 RGF93
2 coordonnées géographiques
037 ellipsoïde IAG-GRS80
02 degrés décimaux
01 méridien international (Greenwich)
GR3D1 valeurs en degrés -5.5000 longitude minimale
10.0000 longitude maximale
41.0000 latitude minimale
52.0000 latitude maximale
.1000 pas en longitude (0.1°)
.1000 pas en latitude (0.1°)
GR3D2 texte Interpolation bilinéaire
GR3D3 Codes IGN PREC CM… précision de la transformation

Corps
Valeur Définition
-5.500000000 Longitude dans le système géodésique RGF93
41.000000000 Latitude dans le système géodésique RGF93
-165.027 TX : Translation en X dans le sens NTF → RGF93 2
-67.100 TY : Translation en Y dans le sens NTF → RGF93
315.813 TZ : Translation en Z dans le sens NTF → RGF93
99 code précision
« », L, - caractère f50 :
« », feuille existante
L, feuille fictive en limite de zone d'application de la grille
-, feuille fictive hors zone d'application de la grille
0158 n de feuille 1: 50000

Le processus de transformation de coordonnées est dépendant du sens de la transformation.

L'interpolation est directe à partir des coordonnées géographiques dans le sens RGF93 →
NTF. Elle nécessite un premier calcul approché dans le sens NTF → RGF93.

2
Les 3 paramètres correspondent aux coordonnées, dans le système RGF93, de l'origine du système NTF.

IGN • Service de Géodésie et de Métrologie • 73 avenue de Paris, 94165 Saint-Mandé Cedex • 01 43 98 83 25 • [email protected] 5/7
Schéma du processus de transformation de coordonnées utilisant la grille de paramètres

Rappel : Conversions de coordonnées dans un même système de référence géodésique

Géométrie de l'ellipsoïde
À partir des valeurs du demi-grand axe a et de l’aplatissement f, il est possible de calculer le
demi-petit axe b de l'ellipsoïde, ainsi que son aplatissement e.
b = a (1 − f )

a 2 − b2
e=
a2

Transformation de coordonnées géographiques ↔ cartésiennes

X = [ N + h ] cos φ cos λ
Y = [ N + h ] cos φ sin λ
Z = [ N (1 – e²) + h ] sin φ

e2 = ( a2 – b2) / a2

W = �𝟏𝟏 – 𝒆𝒆² 𝒔𝒔𝒔𝒔𝒔𝒔² 𝝋𝝋

N=a/W
ρ = a (1 – e²) / W3
r = N cos φ

6/7 IGN • Service de Géodésie et de Métrologie • 73 avenue de Paris, 94165 Saint-Mandé Cedex • 01 43 98 83 25 • [email protected]
Transformation de coordonnées cartésiennes ↔ géographiques 3

1 1 − e2
f =−

R= X 2 +Y2 + Z2

Y 
λ = arctg  
X

 Z   e2 a   
=µ arctg  ⋅  (1 − f ) +   
 X 2 + Y 2   R   

 
 Z (1 − f ) + e 2 a sin 3 µ 
ϕ = arctg 
 
 (1 − f )  X + Y − e a cos µ  
2 2 2 3

h  X 2 + Y 2 ⋅ cos ϕ  + [ Z sin ϕ ] −  a 1 − e 2 sin 2 ϕ 

=
   

Transformation de coordonnées géographiques ↔ planes

Ces transformations, qui dépendent de formules spécifiques à chaque projection, ne sont pas
détaillées dans le présent document.

3
Source : Bowring, 1985, « The accuracy of geodetic latitude and height equations », Survey Review, 28
IGN • Service de Géodésie et de Métrologie • 73 avenue de Paris, 94165 Saint-Mandé Cedex • 01 43 98 83 25 • [email protected] 7/7