GÉODÉSIE

Élaboration d’une grille de conversion

de coordonnées géographiques
Application à la région d’Oran
Sid Ahmed MEDJAHED - Ali ZEGGAI

L’utilisation du système GPS conduit à une solution globale, les coordonnées géodésiques MOTS-CLÉS
(longitude GPS & latitude GPS) sont exprimées par rapport au système WGS84. Généralement,
Transformation, WGS84,
chaque pays ou chaque territoire possède sa propre référence géodésique. Afin d’exploiter Nord Sahara, Grille
localement les avantages de l’outil GPS, il est nécessaire de transformer les coordonnées GPS en planimétrique.
coordonnées exprimées dans le système local.
En géodésie, il existe plusieurs modèles de transformation qui permettent le passage d’un système géodésique à un autre.
L’emploi de ces modèles passe par des étapes de traitement plus ou moins complexes (ex. : le choix de l’approche de
transformation, le choix des points d’appui et la sélection des points de contrôle…), d’où la nécessité de mettre en place
une méthodologie de transformation destinée directement à l’utilisateur. Cette approche est appelée conventionnellement
transformation par grille. Elle est basée sur le principe de construction et l’utilisation d’une grille planimétrique (2D)
pour la transformation de la longitude et la latitude entre deux datum géodésiques (WGS84 et Nord Sahara).

E
n Algérie, le système géodésique Nord Sahara 1959 nées obtenues après le traitement et l’ajustement des
constitue l’infrastructure sur laquelle s’appuient tous les observations GPS sont exprimées dans un système de
travaux géodésiques et cartographiques. Ce système a référence géodésique, géocentrique, tridimensionnel et
été déduit du système Europe 1950 (ED50) par la transforma- mondial appelé WGS84 basé sur l’ellipsoïde WGS84
tion suivante : caractérisé par :
λNS =λED50 -0Gr.0004 Demi grand axe (a) : 6378137.00 m
£NS=£ED50-0Gr.00484 Excentricité (e) : 0.0818191980426

Le système Nord Sahara 1959 est basé sur l’ellipsoïde de L’utilisation du GPS en Algérie est souvent liée à la transfor-
Clarke 1880A (recommandé à la conférence de Bukavu au mation des coordonnées. Les résultats obtenus sont exprimés
Zaïre) et dont les paramètres sont [8] : dans le système en vigueur en Algérie (Nord Sahara).
Demi grand axe (a) : 6378249.145 m Le présent article a pour but de présenter l’approche de trans-
Aplatissement (f) : 293.46500 formation par grilles. L’exploitation de la grille est générée et
utilisée par un programme (grand public) pour l’exploitation
Outre sa vieillesse, le système géodésique national Nord
des coordonnées dans une cartographie existante basée sur
Sahara connaît de nombreux problèmes dus essentiellement
le système géodésique national.
aux techniques utilisées pour sa réalisation, parmi ces pro-
blèmes, nous citons [1] :
• la mauvaise définition de l’origine de ce système car il est Méthodes de transformation
déduit du système ED50 ;
entre systèmes géodésiques
• la non-homogénéité du réseau qui est due à la combinaison
entre le réseau du nord et le canevas du sud, issus de deux Il existe plusieurs modèles de transformation qui permettent
sources d’observations différentes (géodésique et astrono- le passage d’un système géodésique à un autre :
mique) ; • modèles globaux : le modèle de Bursa Wolf, le modèle de
• le manque d’information de la hauteur du géoïde relative à Molodensky-Badekas, le modèle de Veis, le modèle géogra-
l’ellipsoïde de Clarke 1880A (Nord Sahara) ; phique de Molodensky...
• la couverture partielle du réseau sur le territoire national. • modèles locaux : le modèle des lignes géodésiques, la trans-
formation polynomiale par la multiple régression (MREs)…
Aujourd’hui, l’ensemble des réseaux géodésiques natio-
naux, continentaux et mondiaux sont réalisés par les tech- Le choix d’un modèle de transformation adéquat est influencé
niques de positionnement spatial, parmi ces systèmes, par [6] :
l’outil GPS est le plus utilisé actuellement, les coordon- • le domaine pour lequel le modèle doit être appliqué q
Revue XYZ • N° 131 – 2e trimestre 2012 47
 GÉODÉSIE

q • la présence des distorsions dans l’un ou l’autre système géo- Application de la MREs
désique ou bien les deux à la fois. dans la région d’Oran
• le type de transformation souhaitée (2D ou 3D), et le sys-
L’objectif de l’application est de déterminer les paramètres de
tème de coordonnées utilisé (géographique, cartésien ou
transformation bidimensionnelle entre les coordonnées géo-
plan)
désiques rapportées à l’ellipsoïde de révolution WGS84 et
• la précision recherchée suivant le type d’application.
celles rapportées à l’ellipsoïde de Clarke 1880A associé au
système Nord Sahara. Ces paramètres sont les coefficients Ai
Transformation polynomiale et Bi calculés par les équations de la transformation polyno-
par la multiple régression (MREs) miale.
Les données utilisées sont collectées de la manière suivante :
L’emploi d’une méthode de transformation pour le calcul données GPS : issues des observations menées lors des cam-
des paramètres de passage entre deux référentiels géodé- pagnes GPS à travers la région d’Oran
siques exige l’homogénéité entre ces deux référentiels. données Nord Sahara : issues des observations de la trian-
Certains référentiels géodésiques locaux réalisés par la gulation et de la géodésie classique en Algérie.
méthode classique, présentent un manque d’information
L’équation de la régression sur λ s’écrit :
en termes d’orientation et d’échelle et par conséquent le
∆λ = −0,000000006−0,000025582. U + 0,000020706. V
problème de transformation se pose.
−0,000153838. U2 − 0,000036369 U.V + 0,000040856. V2+
Diverses approches ont été proposées, parmi ces
0,000822261. U3 − 0,001138648. U2.V − 0,001457295 . U.V2
méthodes, la transformation polynomiale par multiple
−0,000546643. V3
régression ou MREs [4]. En termes simples ce sont des
fonctions de polynômes qui représentent les distorsions L’équation de la régression sur ϕ s’écrit :
entre les deux systèmes géodésiques, en fonction de la ∆ϕ = 0,000000001 − 0,000009969.U + 0,000046024.V +
position en latitude et en longitude [6] : 0,000067966.U2 − 0,000120882.U.V + 0,000006113. V2 −
λ2 = λ1 + ∆λ 0,000321816. U3 + 0,000986396.U2.V − 0,000583960. U.V2 +
ϕ2 = ϕ1 + ∆ϕ 0.000350958. V3
Avec :
Avec : U = λi − 0°508837967
λ1, ϕ1 : Les coordonnées géodésiques dans le premier sys-
V=ϕi+ 35°594102044
tème (WGS84 )
λ2, ϕ2 : Les coordonnées géodésiques dans le deuxième Les résultats de la transformation sont résumés dans les
système ( Nord Sahara 1959 ) tableaux suivants :
∆λ : La différence en longitude entre les deux systèmes
Sur les quatorze (14) points d’appui
géodésiques.
∆ϕ : La différence en latitude entre les deux systèmes géo- Composante Min Max RMS RMS Calcul
désiques.
dλ (Sec) -0.0043 0.0051 0.0021
Pour ∆λ, la formule de la régression multiple s’écrit :
dϕ (Sec) -0.0044 0.0063 0.0023
∆λ =A0+ Ordre 0 0.088 m
A1U+A2V+ Ordre 1 dE (m) -0.1074 0.1284 0.0536
A3U2+A4UV+A5V2+ Ordre 2
dN(m) -0.1351 0.1943 0.0702
A6U3+A7U2V+A8UV2+A9V3+ Ordre 3
A10U4+A11U3V+A12U2V2+A13UV3+A14V4+ Ordre 4
Tableau 1. Résultats sur les points d’appui
A92U13 + A93 U12V+…..+ A105V13 Ordre 13

Pour ∆ϕ, la formule de la régression multiple s’écrit : Sur les quatre (04) points de contrôle
∆ϕ = B0+ Ordre 0
Composante Min Max RMS RMS Calcul
B1U+B2V+ Ordre 1
B3U2+B4UV+B5V2+ Ordre 2 dλ (Sec) -0.0023 0.0009 0.0016
B6U3+B7U2V+B8UV2+B9V3+ Ordre 3
dϕ (Sec) -0.0055 0.0050 0.0041
B10U4+B11U3V+B12U2V2+B13UV3+B14V4+ Ordre 4 0.133 m
B92U13 + B93 U12V+….……..+ B105V13 Ordre 13 dE (m) -0.059 0.024 0.041
Où dN(m) -0.172 0.155 0.127
Ai et Bi : coefficients de la MREs.
U=K.(ϕ − ϕ0) : latitude géodésique normalisée du point de Tableau 2. Résultats sur les points de contrôle
calcul.
V=K.(λ − λ0) : longitude géodésique normalisée du point de Les écarts représentent les différences entre les coordonnées
calcul. réelles obtenues dans le cadre de la géodésie classique (sys-
(λ0, ϕ0) : coordonnées origine du premier système. tème Nord Sahara) et celles calculées par l’approche de trans-
K : facteur d’échelle. formation polynomiale. L’écart moyen obtenu sur les points

e
48 Revue XYZ • N° 131 – 2 trimestre 2012
 Étude statistique et validation
des résultats
Dans le but d’étudier les résultats de la transformation des
coordonnées par la MREs deux tests statistiques seront
examinés.

1- Test de Khi deux : Ce test a pour but d’éliminer les grosses

erreurs. Les observations issues d’une population sont accep-
tées si : [5]

Où est la valeur de la distribution de Khi-deux à un

degré de liberté et un niveau de signification α.

La probabilité est donc :

Le test statistique de khi-deux est donné par la formule sui-

vante :

σo : c’est l’écart-type a priori prise égale à ± 0. 005 seconde en

référence à la précision des points d’appui Nord Sahara sur la
Figure 1. Le repère géodésique local
zone d’Oran (±15 cm) [2]
de contrôle est de l’ordre de 14 cm. La précision de la trans-
2-Test de Student’s : Le deuxième test est le test de Student’s,
formation est acceptable (< 20 cm) pour les applications géo-
l’observation est acceptée si : [5]
désiques à moyenne précision. L’approche de la transforma-
tion par les équations de la multiple régression peut être
appliquée pour transformer les coordonnées WGS84 vers le Où t : c’est la valeur tabulée de la distribution de Student’s
système national Nord Sahara. —
m : la moyenne des résidus
Remarque : Les différences dλ et dϕ sont respectivement les
σ
^ : c’est l’écart-type du vecteur des résidus
écarts en longitude et latitude. Ces écarts sont exprimés (cf. :
tous les tableaux) en seconde sexagésimale [1]. Les diffé- Pour le degré de liberté (ddl) = 4, les résultats des tests statis-
rences dE et dN correspondent aux valeurs de dλ et dϕ sur tiques sont résumés dans les deux tableaux ci-dessous.
la surface de la terre suivant l’Est et le Nord (E : Easting,
Les résultats de l’étude statistique indiquent que le test de
N : Northing). Les écarts correspondent à :
khi-deux est positif, ceci signifie que les résidus suivent la
dE = dX = Ncosϕ.dλ
loi normale et le test de Student’s indique qu’aucune des
dN = dY = ρdϕ
observations n’a été rejetée.
Avec N et ρ : Respectivement la normale et le rayon de cour-
Enfin et d’après les résultats obtenus, on conclut que l’ap-
bure et comme dE et dN sont petits, nous pouvons appliquer
proche de la transformation par les équations de la régres-
les formules approchées suivantes :
sion multiple employée pour la transformation des coor-
dE = dX = Rcosϕ.dλ
données GPS dans notre application (Région d’Oran) est
dN = dY = Rdϕ
suffisante et elle sert à élaborer la grille planimétrique de
Où R est le rayon moyen de la terre égal à 6 371 km. conversion de coordonnées.

sur la longitude
M
σ0 a priori σ
^0 a postériori t-distribution de Student’s Test de khi-deux Test de Student
R
E
0. 005 0.0021 95% 9.488 0.7301 95% 2.776 Positif O obs suspecte
S

Tableau 3. Test de Student et de khi-deux sur la longitude

sur la latitude
M
σ0 a priori σ
^0 a postériori t-distribution de Student’s Test de khi-deux Test de Student
R
E
0. 005 0.0023 95% 9.488 0.8276 95% 2.776 Positif O obs suspecte
S

Tableau 4. Test de Student’s et khi-deux sur latitude

Revue XYZ • N° 131 – 2e trimestre 2012 49

 GÉODÉSIE

Transformation par la grille planimétrique Si : correspondent à la distance angulaire entre le point cible
et le point source défini par :
La transformation de λ et ϕ entre les deux datum consiste en
la construction d’une surface de correction de coordonnées,
cette surface représente le décalage en λ et ϕ entre les deux
datum. Chaque nœud de la grille contient les informations
Avec : (λ,ϕ) sont des coordonnées géographiques relatives
nécessaires pour transformer les coordonnées. Il s’agit de la
aux données sources (λ0,ϕ0) et à la grille cible (λi,ϕi)
distorsion en longitude (dλGrille) et en latitude (dϕGrille).
3e Étape : Transformer les coordonnées λ1 et ϕ1 exprimées
Le principe de l’utilisation de la grille de conversion planimé-
dans le 1er datum vers le 2e datum (λ2 et ϕ2) en utilisant les for-
trique est le suivant :
mules de la MREs suivantes :
1re Étape : Positionner le point à transformer sur une maille de
la grille et les 4 nœuds de la maille sont à rechercher.
2e Étape : Calculer la distorsion en longitude et en latitude
par une méthode d’interpolation (deux interpolateurs
seront utilisés) :

1 - Interpolation bilinéaire : Soit une maille élémentaire consti-

Élaboration de la grille planimétrique
tuée de quatre nœuds notés 1,2,3 4 délimitée par [2] :
dans la région d’Oran
L’élaboration de la grille de conversion planimétrique dans la
Les longitudes : λ1, λ2, λ3, λ4
région d’Oran a été effectuée en utilisant les mêmes coeffi-
Les latitudes : ϕ1, ϕ2, ϕ3, ϕ4
cients calculés par la MREs dans la transformation 2D dans la
Avec : λ1= λ2, λ3= λ4 et ϕ1= ϕ3, ϕ2 =ϕ4
région d’Oran. Chaque grille plane est fournie sous forme
d’un fichier (Texte/Binaire) selon la configuration suivante :
• En tête : λ Min, λMax, ϕ Min, ϕMax, ∆λ , ∆ϕ , Nc , Nl

• Corps de la grille : dλ , dϕ
La grille plane calculée sert à la transformation de la longitude
et la latitude obtenue par GPS vers le datum lié au système
Nord Sahara. Les nœuds de la grille planimétrique sont expri-
més en coordonnées géographiques.
Figure 2. Interpolation de dλ et dϕ
Afin de tester la validité de la grille 2D pour la transformation
des coordonnées GPS en coordonnées Nord Sahara, une
Soit un point ’’P’’ appartenant à cette maille de coordonnées
étude comparative a été réalisée sur le point (305B), ce point
(λ,ϕ), on obtient la valeur de la variation de la longitude (dλ) et
a été mesuré par observations classiques et par observations
la valeur de la variation de latitude (dϕ) en fonction des
GPS, avec :
valeurs aux nœuds (dλ1,dλ2,dλ3,dλ4) et (dϕ1,dϕ2,dϕ3,dϕ4) par
l’interpolation bilinéaire de la façon suivante : WGS84 (GPS ) NS59 (Triangulation)
La distorsion en longitude : dλP -0°313807883 -0°312976469
λ
ϕ 35°858269158 35°858163164

Tableau 5. Les coordonnées du point 305

La distorsion en latitude : dϕP
La position du point 305B et les limites de la grille 2D calculée
dans la région d’Oran sont représentés dans le schéma suivant

Avec : et

2- Interpolation par le plus proche voisin (PPV) : la méthode

d’interpolation du PPV consiste à attribuer à chaque site Si
cible dans la sélection un poids inversement proportionnel à
la distance entre ce site et le point à estimer S0, dans notre cas
on obtient la distorsion en longitude et en latitude par la for-
mule suivante [3] :

et

Figure 3. Position du point 305B

Avec : Les distorsions en longitude et en latitude sur les Quatre
n0 : Le nombre de données prises en compte (le nombre de nœuds (maille) avoisinant le point 305B sont représentées
voisins utilisés) dans le tableau suivant :

e
50 Revue XYZ • N° 131 – 2 trimestre 2012
 Nœud N° 01 Nœud N° 02 à l’utilisateur de transformer directement leurs coordonnées
GPS sans passer par les étapes de traitement des logiciels de
λ -0°314 λ -0°314
transformations disponibles en Algérie et elle permet aussi
ϕ 35°858 ϕ 35°859
d’éviter les travaux de rattachement en planimétrie sur les
dλ 0°000574972 dλ 0°000572918 points géodésiques. ●
dϕ -0°000063711 dϕ -0°000063288
Nœud N° 03 Nœud N° 04
Contacts
λ -0°313 λ -0°313
ϕ 35°858 ϕ 35°859 Sid Ahmed MEDJAHED - Ali ZEGGAI
dλ 0°000573838 dλ 0°000571776 [email protected]
dϕ -0°000063637 dϕ -0°000063215 Centre des Techniques Spatiales –Division de Géodésie Spatiale
(Arzew- Oran/Algérie)
Tableau 6. Calculs sur les nœuds avoisinant le point 305B
Les distorsions calculées sur le point 305B par l’emploi de l’in- Bibliographie
terpolation bilinéaire (IBL) et l’interpolation du plus proche
[1] S.A.Medjahed (Février 2006) Intégration des données GPS, Nord
voisin ( PPV) sont :
Sahara et de Nivellement en vue d’une transformation
dλIBL : 0°000574201
tridimensionnelle. Cas de l’Algérie. Mémoire de Magister CNTS
dϕ IBL :-0°000063583
[2] S.A Medjahed (2008) Génération de grilles de conversions de
dλPPV : 0°000573887
coordonnées géodésiques - Rapport de Recherche –CTS
dϕPPV :-0°000063540
[3] K.Kokoche Omar (2007) Élaboration d’une grille de conversion
Le tableau suivant illustre l’écart en longitude et en latitude
altimétrique pour la Mans Métropole-mémoire d’ingénieur – École
calculé entre la position réelle déterminée par mesures géo- Supérieure des Géomètres et Topographes (ESGT).
désiques et la position calculée par la grille de transformation
[4] European Petroleum Suvey Group (EPSG) - Guidance Note
planimétrique : Number 7 – Coordinate Conversions and Transformations including
Formulas.
Position Différence [5 ] A.Zeggai, S.A Benahmed Daho, B. Ghezali, H.Taibi , R. Ait
Position calculée par grille
mesurée (Sec) Ahmed Lamara (2006) Conversion altimétrique des hauteurs
λIBL -0°313233682 2.0671
ellipsoïdales par GPS. Revu XYZ N° 109 – 4e trimestre 2006
0°312976469 [6] Merrin Pearse and Chris Crook Crook (Report: GS 1997/11)
λPPV -0°313233996 2.0660 Recommended transformation parameters from WGS84 to NZGD49-
ϕIBL 35°858205575 -0.2289 Geodetic System Technical Report - Toitu te Land Whenua
35°858163164 Information New Zealand.
ϕPPV 35°858205618 -0.2287
[7] NATIONAL IMAGERY AND MAPPING AGENCY (2000) – Its Definition
Tableau 7. Résultats de la transformation du point 305B and relations ships with local geodetic systems –TECHNICAL REPORT.
[8] A .Reyt (1960) Les principales étapes de la géodésie en Afrique
On remarque qu’il existe une petite différence entre la posi- du Nord et leur incidence cartographique. éditeur. Imprimerie IGN
tion réelle et la position calculée en utilisant la grille de conver-
sion planimétrique. ABSTRACT
La précision de la grille 2D est liée à la précision de l’approche Words Key: Transformation, WGS84, Nord Sahara, two
de transformation employée pour sa construction, il s’agit de dimensional grid.
la MREs dans notre application. The use of the GPS system, conducted us to a global solution,
the geodetic co-ordinates (GPS longitude & GPS latitude) are
expressed in relation to the WGS84 system. Generally, every
Conclusion country or each territory possesses its geodetic reference. In
Les modèles de transformation à trois dimensions restent mal order to exploit advantages of the GPS tool locally, it is
adaptés pour le cas de l’Algérie à cause du manque d’informa- necessary to transform the GPS co-ordinates in the local system.
tion sur la hauteur du géoïde relatif à l’ellipsoïde de Clarke 1880A. In geodesy there are several models of geodetic transformations
Afin de résoudre le problème nous avons transformé les deux that allow the passage of a geodetic system another. The using
composantes (longitude et latitude) par l’emploi d’une transfor- of these models through the processing steps more or less
mation polynomiale type régression multiple (MREs). Les résul- complex (e.g. the choice of the transformation approach, the
tats obtenus dans la région d’Oran montrent que cette approche choice of support points, selection controls points etc.).From
peut être appliquée pour transformer les observations GPS vers where the necessity to put an intended transformation
le système national Nord Sahara avec une précision de 14 cm, methodology in place directly to the user. This approach is called
une grille planimétrique (2D) pour la transformation des coor- transformation conventionally by grid and it based on the
donnés géodésiques dans la région d’Oran a été élaborée. principle of construction and the use of a two dimensional grid
La grille élaborée dans la région d’Oran représente la surface (2D) for the transformation of the longitude and the latitude
de correction de cordonnées entre les deux datum. Elle permet between two geodetic datum (WGS84 and North Sahara).

Revue XYZ • N° 131 – 2e trimestre 2012 51