Systèmes de référence cartographique

Rappels

Les systèmes de référence sont aussi appelés systèmes géodésiques ou en

anglais Datum ou Spatial reference.
En géométrie, on vous a appris à représenter des points aux coordonnées connues
dans un repère orthonormé bidimensionnel ou tridimensionnel défini auparavant.
On commence d'abord par définir ce repère bidimensionnel ou tridimensionnel,
ensuite on y représente des points grâce à leurs coordonnées (X,Y,Z). En reliant
ces points, on peut retrouver des formes géométriques comme celles d'un carré,
d'un rectangle, etc.
De la même manière, en topographie où l'on doit lever des détails et en
cartographie, où l'on doit représenter ces détails sous la meilleure forme, on doit
également définir un système de coordonnées de référence pour y rattacher tous
ces détails.

Pour cela, il faut d'abord définir le système géocentrique cartésien:

2- 1 Système géocentrique cartésien

Les détails du terrain sont en 3D, le système de référence est alors toujours
tridimensionnel. Lorsqu'on construit mathématiquement un système de référence
⃗⃗ ), on parle forcément de coordonnées X,Y,Z (Fig.1). Les coordonnées de
(o, 𝑖⃗, 𝑗⃗, 𝑘
base d'un système géodésique (o, 𝑖⃗, 𝑗⃗, 𝑘 ⃗⃗ ) sont des coordonnées cartésiennes
géocentriques (liées à la Terre) X, Y, Z tels que:

O est proche du centre de la Terre. Z

(O, 𝑖⃗, 𝑗⃗): plan de l'équateur
(o, 𝑘⃗⃗ ): axe de rotation de la terre
(o, 𝑖⃗, 𝑘⃗⃗ ): plan du méridien origine
M

O Y

Figure 1: Référentiel cartésien

géocentrique
X
Le système cartésien géocentrique est utilisé par les scientifiques seulement. Cela
dit, les coordonnées mesurées par le GPS sont des coordonnées géocentriques
cartésiennes.

Exemple:
Un point au niveau de l'aéroport:
X = 5108608,479, Y= 268484,128, Z=3797047.925

ENSSMAL Benzouaï S. 2020/2021

 2-2 Systèmes de coordonnées géographiques :
Pour passer à des coordonnées géographiques, le choix d'une surface
mathématique qui représente la Terre s'impose. Rappelez-vous que depuis le
primaire, on vous a toujours appris que les longitudes se mesurent par rapport au
méridien de Greenwich et que les latitudes se mesurent par rapport à l'équateur,
tout cela sur une sphère qui serait la Terre.
Et bien maintenant, il faut savoir que représenter la Terre par une sphère est tout
simplement une approximation (pour simplifier les choses).
En réalité, la forme mathématique qui représente au mieux la Terre est l'ellipsoïde
de révolution.

Par conséquent, on munit le référentiel géocentrique cartésien d'un ellipsoïde

(Fig.2). Un ellipsoïde est défini par 2 paramètres (a,b) ou (a,e), (a, é), (a, f), etc.

- a: demi-grand axe (~6378 km)

b
- b: demi-petit axe (~6356 km)

𝑎2 − 𝑏 2 a
- 1ère excentricité : 𝑒 2 = ,
𝑎2
𝑎2 − 𝑏2
- 2ème excentricité: 𝑒′2 = 𝑏2
Figure 2: Paramètres d'un ellipsoïde
𝑎−𝑏
- f: aplatissement: 𝑓 = 𝑎

On le munit également d'un :

- méridien origine qui est le méridien de Greenwich (décision de l'association
internationale de géodésie).
- d'une unité angulaire: DMS (Degré Minute Secondes) ou dans la norme NMEA
(Degré Minute Décimale)
- d'une unité de longueur: mètre

De cette manière, on peut définir des coordonnées géographiques (l,f,h) avec h

: hauteur sur l'ellipsoïde (Fig.3).
Z

M
h

Y
O

X
Figure 3: Coordonnées géographiques

On définit alors:

ENSSMAL Benzouaï S. 2020/2021

 l: angle entre le plan du méridien origine et le plan contenant M compté de O, 180 E
et 0, 180 W ou (O,360)
f: angle entre le plan de l'équateur et la normale à l'ellipsoïde contenant le point M
h: hauteur au-dessus de l'ellipsoïde.

Exemple:
Voici les coordonnées géographiques du même point de l'aéroport:
l= 3°00'30.33173 f= 36°46'07.11280 h = 318.892

X,Y,Z -----------> l,f,h

Il est possible de passer des coordonnées géocentriques cartésiennes aux

coordonnées géographiques et vis-versa.

2-3 Systèmes de coordonnées planes :

Maintenant, nous voulons représenter ces détails 3D connus en coordonnées
géographique l,f,h sur un plan. Le seul moyen pour y arriver est d'utiliser une
fonction mathématique qui n'est autre qu'une projection (Fig. 4).

N(y)
Z
F
PROJECTION

O
Y

E(x)
X
Surface courbe (ellipsoïde) Surface plane
Figure 4: Passage du 3D à une représentation plane 2D

F est appelée "Projection", une fonction mathématique qui permet de passer d'une
surface courbe à une surface plane.
Nous venons ainsi de définir les 3 éléments indispensables et indissociables
rentrant dans la définition d'un système géodésique ou système cartographique
de référence.

ENSSMAL Benzouaï S. 2020/2021

 Systèmes de coordonnées géocentriques cartésiennes
-----> Coordonnées cartésiennes géocentriques
X,Y,Z

Systèmes de coordonnées géocentriques cartésiennes +

Ellipsoïde
l,f,h -----> Coordonnées géographiques

Systèmes de coordonnées
géocentriques cartésiennes +
Ellipsoïde + Projection
E1,N1 E2,N2 -----> Coordonnées planes

Il est à noter que les projections peuvent conserver:

- les angles, on les appelle projections conformes: ce type de projection est utilisé
dans les cartes topographiques, géologiques et marines. En effet, dans ce type de
cartes la notion d'angle et de forme est très importante.
- les surfaces, on les appelle projections équivalentes: utile dans certaines cartes
thématiques où il est important de montrer la préservation des superficies (carte
agricole, carte du monde, etc.)

Une projection ne conserve jamais les distances (altération linéaire) raison pour
laquelle on construit plusieurs projections de façon à minimiser ces altérations en
s'arrangeant toujours à avoir des coordonnées positives.

Pour cela une surface mathématique de projection (Fig.5) est choisie : cylindre, cône,
plan, etc.

ENSSMAL Benzouaï S. 2020/2021

 Figure 5: Surfaces de projection cylindrique, conique et azimutale

La surface de projection est soit sécante ou tangente. Chaque projection a ses propres
paramètres (méridien(s) central(aux), parallèle(s) de référence(s), point tangent, etc.)

Ces projections portent souvent des noms particuliers permettant de les distinguer.
Exemple d'une projection conique: Projection Lambert. Exemple d'une projection
cylindrique : Projection UTM (Universal Transverse Mercator).
La projection UTM est celle utilisée en Algérie, la surface de projection est un cylindre
disposé horizontalement. La Terre est subdivisée en 60 fuseaux ou zones d'amplitude
6°. L'Algérie est couverte par le fuseau 29, 30, 31 et 32 sachant que le méridien de
Greenwich 0° passe par Stidia (Mostaganem).
Fuseau 29 (de 12°W à 6°W) avec un méridien origine à 9° W;
Fuseau 30 (de 6°W à 0°) avec un méridien origine à 3° W;
Fuseau 31 (de 0° à 6°E) avec un méridien origine à 3° E;
Fuseau 32 (de 6° à 12°E) avec un méridien origine à 9° E.

Les paramètres de la projection UTM :

- Le méridien central dont la projection sur le cylindre va représenter l'axe des Y et
l'équateur qui va porter l'axe des X. Dans un souci d'avoir des coordonnées positives,
l'axe des Y est décalé de 500 km vers l'ouest (False X= 500000m). Il n’y a pas besoin
de décaler l'axe des X puisque l'Algérie se trouve bien au-dessus de l'équateur (False
Y=0).
- Un facteur d'échelle égal à E=0.9996.

ENSSMAL Benzouaï S. 2020/2021

 Figure 6: Paramètres de la projection UTM

Il est à noter que le S (Fig. 7) qui apparaît sur les GPS de navigation indique un
identifiant de zone et non pas la zone sud, ceci est valable aussi pour Google Earth.
De même les coordonnées, à 3 chiffres en x et à 4 chiffres en y, lues sur une carte
topographique sont des coordonnées planes (Fig.9) dans le Nord Sahara 1959,
Clarke 1880, UTM (31 pour l'extrait de carte que vous avez géoréférencé en TD). Ce
sont des coordonnées kilométriques, il faut les convertir en mètre (x1000), unité du
référentiel. Sur la même carte, vous trouverez aussi une indication sur les
coordonnées géographiques

ENSSMAL Benzouaï S. 2020/2021

 Figure 7: Grille de numérotation des zones UTM

Coordonnées planes Nord Sahara 1959, Clarke 1880,UTM31

Coordonnées géographiques Nord Sahara 1959, Clarke 1880

En Algérie, il existe 2 Figure 9: Grille

systèmes de numérotation
géodésiques des zones UTM
officiel:

 Le Nord Sahara 1959, projection UTM associée à l’ellipsoïde Clarke 1880, fuseau
29,30,31 et 32.
 Dès l’an 2000, le datum WGS84 (utilisé dans le positionnement spatial GPS) est
entré en vigueur en Algérie.

ENSSMAL Benzouaï S. 2020/2021

L'ancienne édition des cartes topographiques de l'Algérie, héritage de la cartographie
française en Algérie, étaient dans le référentiel Voirol 1875 Nord Algérie qui est basé
sur une projection Lambert (conique).

Ces Systèmes de référence sont répertoriés dans tous les logiciels de SIG et de
Cartographie.

Exemple 1: ARCGIS
Nord_Sahara_1959_UTM_Zone_31N
WKID : 30731 Autorité : EPSG

Projection: Transverse_Mercator
False_Easting: 500000,0
False_Northing: 0,0
Central_Meridian: 3,0
Scale_Factor: 0,9996
Latitude_Of_Origin: 0,0
Linear Unit: Meter (1,0)

Geographic Coordinate System: GCS_Nord_Sahara_1959

Angular Unit: Degree (0,0174532925199433)
Prime Meridian: Greenwich (0,0)
Datum: D_Nord_Sahara_1959
Spheroid: Clarke_1880_RGS
Semimajor Axis: 6378249,145
Semiminor Axis: 6356514,869549776
Inverse Flattening: 293,465

Exemple 2: QGis

ENSSMAL Benzouaï S. 2020/2021

Exemple 3: ARCGIS (Système de coordonnées Voirol Nord Algérie)

Voirol_1875_Nord_Algerie_ancienne
Autorité : personnalisée

Projection: Lambert_Conformal_Conic
false_easting: 500000,0
false_northing: 300000,0
central_meridian: 2,7
standard_parallel_1: 36,0
scale_factor: 0,999625544
latitude_of_origin: 36,0
Linear Unit: Meter (1,0)

Geographic Coordinate System: GCS_Voirol_1875

Angular Unit: Degree (0,0174532925199433)
Prime Meridian: Greenwich (0,0)
Datum: D_Voirol_1875
Spheroid: Clarke_1880_IGN
Semimajor Axis: 6378249,2
Semiminor Axis: 6356515,0
Inverse Flattening: 293,4660212936265

2-4 Système de référence local et global:

Au début, les systèmes de référence étaient définis localement de manière à les
adapter au mieux localement à la région à cartographier (Fig. 9). Le Nord Sahara 1959
en est un exemple.
Avec certaines techniques spatiales (laser, VLBI, etc.) et surtout avec l'apparition des
GNSS, les systèmes de référence globaux sont devenus plus populaires, notamment
le WGS84, utilisé par le GPS.
Ces systèmes de référence s'appuient sur des ellipsoïdes locaux et globaux.
L'ellipsoïde local (Exemple: Clarke 1880) épouse au mieux la forme d'une région
donnée (Exemple: Algérie). Pour cartographier une autre région du monde, on choisit
un autre ellipsoïde.
En revanche, l'ellipsoïde global comme le WGS84 est adapté globalement à toutes les
régions du monde.

Figure 9: Ellipsoïde local, global et géoïde

ENSSMAL Benzouaï S. 2020/2021

2-5 Système de référence des altitudes:
Pour les altitudes, il existe 2 surfaces de référence: le géoïde et l'ellipsoïde.
Le géoïde est la forme dynamique ou physique de la Terre. Elle correspond par
définition au niveau moyen des mers prolongé sous les continents. C'est aussi une
surface équipotentielle (même valeur de la gravité). Les altitudes mesurées par rapport
à cette surface sont dites orthométriques ou verticales. La technique qui permet
l'obtention de ce type d'altitude est le nivellement (voir cours de topographie). Les
repères de nivellement fournissent des altitudes orthométriques, les altitudes, lues sur
une carte topographique, des points cotés et des courbes de niveaux sont des altitudes
orthométriques.
L'ellipsoïde a été défini précédemment; c'est une surface mathématique qui représente
le mieux la forme de la Terre. Les altitudes mesurées par rapport à cette surface sont
dites ellipsoïdales ou élévations.
L'ondulation du géoïde (Fig. 10) est la différence entre les 2 altitudes N= h-H.

Figure 10: Systèmes de référence des altitudes

2-6 Passage d'un référentiel géodésique à un autre

Le passage d'un référentiel géodésique à un autre est une transformation
mathématique connu, la similitude; 3 translations (de l'origine), 3 rotations (des axes)
et un facteur d'échelle.

Les paramètres de cette transformation sont livrés par l'INCT (voir site) mais aussi par
la majorité des logiciels de SIG ou de cartographie.

Tx = 209.36 m
Ty = 87.82 m
Tz = - 404.62 m

Rx = -0".0046 (sexagésimales)
Ry = -3".4784
Rz = - 0".5805

Facteur d'échelle : 1.454722

ENSSMAL Benzouaï S. 2020/2021