LIVRE EN LIGNECartographie humanitaire avec ArcGIS Desktop
📖 Livre ⌂ Site SMX
Module 2 – Systèmes de coordonnées et projectionsPage 4 / 53
MODULE 2

Systèmes de coordonnées et projections

C'est le module le plus important du livret. La plupart des erreurs de cartographie en contexte humanitaire viennent d'un système de coordonnées mal compris.

À LA FIN DE CE MODULE, VOUS SEREZ CAPABLE DE :
  • expliquer la différence entre géoïde, ellipsoïde et datum ;
  • lire et convertir des coordonnées en degrés décimaux et en degrés-minutes-secondes ;
  • expliquer pourquoi on projette, et ce qu'une projection déforme ;
  • choisir la zone UTM adaptée à Goma, Bukavu ou Kisangani, et gérer le cas de l'équateur ;
  • vérifier, définir et changer le système de coordonnées d'une couche dans ArcMap sans corrompre les données.

2.1 La forme de la Terre : géoïde, ellipsoïde, datum

Pour donner la position d'un point, il faut d'abord un modèle de la Terre. La Terre réelle a un relief irrégulier ; on utilise donc des surfaces de référence plus simples.

Surface topographique, géoïde et ellipsoïde (schéma, écarts exagérés).
FIGURE 2.1 — Surface topographique, géoïde et ellipsoïde (schéma, écarts exagérés).

2.2 Le système de coordonnées géographiques

Un système de coordonnées géographiques (en anglais Geographic Coordinate System, GCS) repère un point par deux angles mesurés depuis le centre de l'ellipsoïde : la latitude et la longitude. L'unité est le degré.

Latitude, longitude et convention de signe utilisée par les SIG.
FIGURE 2.2 — Latitude, longitude et convention de signe utilisée par les SIG.
ÉlémentValeur pour WGS 84
Nom dans ArcMapGCS_WGS_1984
Code EPSG4326
Unitédegré
Méridien d'origineGreenwich (0°)
Plage de latitude−90 à +90 (sud négatif)
Plage de longitude−180 à +180 (ouest négatif)

2.3 Degrés décimaux et degrés-minutes-secondes

Une même position peut s'écrire de deux façons. Les SIG travaillent en degrés décimaux (DD). Les cartes anciennes et certains GPS affichent des degrés, minutes, secondes (DMS). Un degré compte 60 minutes, et une minute 60 secondes.

Conversion de DMS en DD

DD = degrés + minutes / 60 + secondes / 3600, puis on met un signe moins si le point est au sud (S) ou à l'ouest (O).

LieuDMSCalculDD (Y ; X)
Goma1°40′48″ S ; 29°13′12″ E1 + 40/60 + 48/3600 = 1,68 ; 29 + 13/60 + 12/3600 = 29,22−1,68 ; 29,22
Bukavu2°30′00″ S ; 28°52′00″ E2 + 30/60 = 2,5 ; 28 + 52/60 = 28,8667−2,5 ; 28,8667
Kisangani0°31′00″ N ; 25°12′00″ E31/60 = 0,5167 ; 25 + 12/60 = 25,20,5167 ; 25,2

Coordonnées arrondies, données à titre d'exemple pour les centres-villes.

Conversion de DD en DMS

  1. Les degrés sont la partie entière : 29,2228 → 29°.
  2. Multipliez la partie décimale par 60 : 0,2228 × 60 = 13,368 → 13′.
  3. Multipliez la nouvelle partie décimale par 60 : 0,368 × 60 = 22,08 → 22″. Résultat : 29°13′22″ E.

Combien de décimales garder ?

Le nombre de décimales indique la précision. Près de l'équateur, un degré de latitude mesure environ 111 km au sol :

DécimalesExemplePrécision au sol (env.)Usage
1−1,711 kmLocaliser une ville
2−1,681,1 kmLocaliser un quartier, un village
3−1,679110 mLocaliser un site de déplacés
4−1,679211 mLocaliser un bâtiment, un forage
5−1,679211,1 mLimite de précision d'un GPS de smartphone

2.4 Les systèmes de coordonnées projetés

Une carte est plane, la Terre ne l'est pas. Une projection est une méthode mathématique qui transforme les latitudes et longitudes en coordonnées planes X et Y, exprimées en mètres. Un système de coordonnées projeté (Projected Coordinate System, PCS) combine un système géographique (datum) et une projection.

Aplatir une surface courbe entraîne toujours des déformations. Une projection ne peut pas conserver à la fois les formes, les surfaces, les distances et les directions. On choisit donc la projection selon ce que l'on veut mesurer :

Propriété conservéeNomUtile pourExemple
Les angles et les formes localesConformeNavigation, cartes topographiques, mesures localesMercator, UTM, Lambert conforme conique
Les surfacesÉquivalenteComparer des superficies, densités, statistiques par zoneAlbers conique équivalente
Les distances depuis un point ou le long de certaines lignesÉquidistanteRayons d'action, cartes de distanceÉquidistante azimutale
Les trois grandes familles de projections, selon la surface de projection.
FIGURE 2.3 — Les trois grandes familles de projections, selon la surface de projection.

2.5 Le système UTM et les zones de la RDC

Le système UTM (Universal Transverse Mercator) découpe la Terre en 60 zones de 6° de longitude, numérotées de 1 à 60 vers l'est à partir de 180° O. Chaque zone utilise sa propre projection transverse de Mercator, centrée sur un méridien central. Les déformations restent faibles à l'intérieur d'une zone, ce qui rend l'UTM très pratique pour mesurer des distances et des surfaces à l'échelle d'une province.

ParamètreValeur
Largeur d'une zone6° de longitude
Numéro de zonepartie entière de (longitude + 180) / 6, plus 1
Facteur d'échelle sur le méridien central0,9996
Fausse abscisse (False Easting)500 000 m
Fausse ordonnée (False Northing)0 m au nord de l'équateur ; 10 000 000 m au sud
Unitémètre

Exemple pour Goma : (29,22 + 180) / 6 = 34,87 ; partie entière 34, plus 1 : zone 35. La zone 35 couvre les longitudes de 24° E à 30° E, avec un méridien central à 27° E. Goma, Bukavu et Kisangani sont toutes les trois dans la zone 35.

Zones UTM couvrant la RDC et position des trois chefs-lieux. La RDC s'étend sur les zones 33 à 36.
FIGURE 2.4 — Zones UTM couvrant la RDC et position des trois chefs-lieux. La RDC s'étend sur les zones 33 à 36.
ZoneLongitudesHémisphère sud (EPSG)Hémisphère nord (EPSG)Exemples
3312° E – 18° E3273332633Kinshasa, Matadi
3418° E – 24° E3273432634Kananga, Mbandaka
3524° E – 30° E3273532635Goma, Bukavu, Kisangani, Lubumbashi
3630° E – 36° E3273632636Bunia, extrême est de l'Ituri

2.6 Cas particulier : le Nord-Kivu et la Tshopo à cheval sur l'équateur

Goma (1,68° S) et Bukavu (2,51° S) sont au sud de l'équateur. Kisangani (0,52° N), Beni (0,49° N) et Butembo (0,13° N) sont au nord. Le Nord-Kivu et la Tshopo sont donc traversés par l'équateur. En UTM, les hémisphères nord et sud n'ont pas la même fausse ordonnée : une même province se retrouve « coupée » entre 35N et 35S.

Trois stratégies sont possibles. Elles sont correctes à condition d'être appliquées à toute l'équipe et écrites dans les méthodes du bureau.

BesoinSystème recommandéPourquoi
Stocker et échanger les donnéesWGS 84 géographique (EPSG:4326)Format attendu par Kobo, GPS, HDX et les partenaires.
Mesurer distances et zones tampons dans une provinceUne seule zone UTM pour toute la province, par exemple WGS 1984 UTM Zone 35S (EPSG:32735), y compris pour les points juste au nord de l'équateurLa projection reste valide de part et d'autre de l'équateur ; les points au nord ont simplement une ordonnée supérieure à 10 000 000 m.
Comparer des superficies entre provincesProjection équivalente : Africa Albers Equal Area Conic (ESRI:102022)Conserve les surfaces à l'échelle du continent.

2.7 Pourquoi c'est important : quatre erreurs réelles

ErreurCe qui se passeConséquence
Latitude et longitude inverséesGoma (Y = −1,68 ; X = 29,22) importé avec X = −1,68 et Y = 29,22Le point apparaît dans le nord de l'Algérie, à plus de 3 500 km.
Signe moins oubliéY = 1,68 au lieu de −1,68Le point est déplacé d'environ 370 km vers le nord, de l'autre côté de l'équateur.
Distance calculée en degrésDistance Goma–Bukavu mesurée dans un système géographique : environ 0,9 « degré »Un résultat sans unité au sol ; la bonne valeur est d'environ 100 km à vol d'oiseau, calculée dans un système projeté ou par une formule géodésique.
Mauvais outil de projectionDefine Projection utilisé pour « convertir » une couche UTM en WGS 84Les coordonnées ne changent pas, mais l'étiquette est fausse : la couche disparaît de la zone ou se retrouve à des milliers de kilomètres.

Superficies : une maille de 0,1° × 0,1° autour de Bukavu couvre environ 123 km². Dans un système géographique, ArcMap ne peut pas donner directement cette surface en km², car l'unité est le degré et la taille d'un degré de longitude varie avec la latitude.

La projection « à la volée »

Le bloc de données (data frame) d'ArcMap possède son propre système de coordonnées. Il prend par défaut celui de la première couche ajoutée. Les couches suivantes sont reprojetées à l'affichage pour se superposer correctement : c'est la projection à la volée. Les fichiers sur le disque ne sont pas modifiés. Cette fonction ne marche que si chaque couche a un système de coordonnées défini.

2.8 Dans ArcMap : vérifier, définir, projeter

a) Connaître le système de coordonnées d'une couche

b) Le système du bloc de données

Data Frame Properties › Coordinate System : (1) onglet, (2) choix de WGS 1984 UTM Zone 35S, (3) détail des paramètres.
FIGURE 2.5 — Data Frame Properties › Coordinate System : (1) onglet, (2) choix de WGS 1984 UTM Zone 35S, (3) détail des paramètres.

c) Définir une projection (Define Projection)

Define Projection sert uniquement à écrire l'étiquette d'un jeu de données dont le système est inconnu mais connu de vous (par exemple un shapefile créé depuis Kobo, en WGS 84, sans fichier .prj). L'outil ne modifie pas les coordonnées.

Outil Define Projection (ArcToolbox › Data Management Tools › Projections and Transformations).
FIGURE 2.6 — Outil Define Projection (ArcToolbox › Data Management Tools › Projections and Transformations).

d) Projeter (Project)

Project crée un nouveau jeu de données dont les coordonnées sont recalculées dans le système de sortie. C'est l'outil à utiliser pour obtenir une copie en UTM 35S d'une couche en WGS 84.

Outil Project : (1) couche d'entrée, (2) couche de sortie, (3) système de sortie.
FIGURE 2.7 — Outil Project : (1) couche d'entrée, (2) couche de sortie, (3) système de sortie.

e) L'avertissement de transformation

Avertissement affiché quand une couche utilise un autre datum que le bloc de données.
FIGURE 2.8 — Avertissement affiché quand une couche utilise un autre datum que le bloc de données.

Cet avertissement signale que la couche et le bloc de données n'utilisent pas le même datum. Cliquez sur Transformations… et choisissez la transformation adaptée à la région. Ignorer ce message peut créer des décalages de plusieurs dizaines à plusieurs centaines de mètres selon les datums.

f) Calculer une superficie ou une longueur

g) Importer des coordonnées depuis Kobo ou Excel

2.9 Exercices

EXERCICE 2.ASuperficie des zones de santé du Sud-Kivu45 min

Objectif. Comparer des superficies calculées dans plusieurs systèmes et comprendre pourquoi les degrés ne conviennent pas.

Données. Zones de santé du Sud-Kivu (jeu du livret, dérivé des limites du Référentiel Géographique Commun publiées sur HDX).

Étapes

  1. Ouvrez un document vierge et ajoutez 01_Limites_admin/SudKivu_ZonesSante.shp. Vérifiez son système (attendu : GCS_WGS_1984).
  2. Dans la table attributaire, ajoutez trois champs Double : Sup_UTM, Sup_ALB, Ecart_pct.
  3. Réglez le bloc de données sur WGS 1984 UTM Zone 35S, puis calculez Sup_UTM en km² avec Calculate Geometry.
  4. Réglez le bloc sur Africa Albers Equal Area Conic, puis calculez Sup_ALB en km².
  5. Calculez Ecart_pct = (Sup_UTM − Sup_ALB) / Sup_ALB × 100 avec le Field Calculator.
  6. Remettez le bloc en GCS_WGS_1984 et essayez de calculer la superficie : notez ce qu'ArcMap propose.
  7. Enfin, projetez la couche en UTM 35S avec l'outil Project et comparez le champ automatique Shape_Area à Sup_UTM.

Questions

  1. Quel est l'écart maximal entre UTM et Albers ? Est-il important pour un rapport de situation ?
  2. Pourquoi ArcMap refuse-t-il de calculer une surface en km² dans un système géographique ?
  3. Quel système choisiriez-vous pour un tableau de densité de population par zone de santé ? Justifiez.
EXERCICE 2.BChasse aux erreurs dans un export Kobo25 min

Objectif. Repérer et corriger les erreurs de coordonnées les plus courantes avant de cartographier.

Étapes

  1. Lisez le tableau ci-dessous, extrait fictif d'une évaluation de sites autour de Goma.
  2. Pour chaque ligne, dites si les coordonnées sont plausibles. Si non, nommez l'erreur et proposez une correction.
  3. Importez ensuite le fichier complet (04_Exercices/M02/kobo_sites.csv) avec Display XY Data et vérifiez vos réponses sur la carte.
Site_gps_latitude_gps_longitudePrécision (m)
Site A – Lushagala-1.635229.19876
Site B – Bulengo29.1495-1.65388
Site C – Kanyaruchinya1.574129.25345
Site D – Rusayo-1.601229.13111 450
Site E – Mugunga00—
Site F – Kibati-1°33′40″29°16′05″7

Noms de lieux réels, coordonnées et précisions fictives, construites pour l'exercice.

EXERCICE 2.CConversions et zones UTM15 min

Objectif. S'entraîner aux conversions et au choix de zone.

Étapes

  1. Convertissez en DD : Uvira 3°24′00″ S ; 29°08′00″ E.
  2. Convertissez en DMS : Butembo 0,1333 ; 29,2833.
  3. Calculez la zone UTM et l'hémisphère de : Uvira, Butembo, Isangi (0,77° N ; 24,25° E).
  4. Donnez le code EPSG de la zone UTM utilisée pour Isangi.