Introduction à la Géomatique (partie 4/4)

## Objectifs et prérequis **À la fin de ce chapitre, l’élève sera capable de :** - Expliquer le principe du GPS différentiel et ses avantages par rapport au GPS autonome. - Mettre en œuvre une correction différentielle à l’aide d’une station de base et d’un récepteur rover dans un contexte de chantier camerounais. - Identifier les principales sources d’erreur d’unités, réaliser des conversions correctes entre les systèmes SI et les unités usuelles utilisées en topographie. - Appliquer les étapes de géoréférencement d’un plan (numérisation, choix des points de contrôle, transformation affine ou polynomiale) à l’aide d’un SIG libre ou propriétaire. - Réaliser un exercice complet de positionnement différentiel, de conversion d’unités et de géoréférencement, en justifiant chaque étape et en contrôlant la précision obtenue. **Prérequis :** - Notions de base du système GPS (satellites, trilatération, temps de propagation). - Connaissance des unités du Système International (mètre, seconde, kilogramme) et des unités usuelles en topographie (mètre, degré, gon, mètre carré). - Familiarité avec la lecture d’un plan topographique (échelle, coordonnées, système de projection). - Utilisation élémentaire d’un logiciel de SIG (ouverture d’un shapefile, affichage de couches). --- ## Introduction contextualisée La géomatique regroupe l’acquisition, le traitement et la diffusion des informations géoréférencées. Au Cameroun, les projets d’aménagement du territoire, de voirie urbaine, de gestion des ressources naturelles et de cadastre exigent une précision métrique de l’ordre du centimètre à quelques décimètres. Le GPS autonome fournit généralement une précision de 5 à 10 m, insuffisante pour ces applications. Le **GPS différentiel** corrige en temps réel les erreurs communes aux récepteurs situés dans une même zone, permettant d’atteindre une précision sub‑déci‑métrique. Parallèlement, la manipulation des unités (conversion entre degrés, gon, mètres, pieds) reste une source fréquente d’erreur lorsqu’on passe d’un logiciel à un autre ou lorsqu’on travaille avec des données provenant de différents services (IGN, MINTP, ONCF). Enfin, le **géoréférencement d’un plan** consiste à associer les coordonnées réelles du terrain à une image ou un dessin numérisé, afin de pouvoir le superposer à d’autres couches géospatiales dans un SIG. Ce chapitre présente, de façon progressive, les trois notions essentielles pour maîtriser ces étapes dans un contexte professionnel camerounais. --- ## I. Leçon 9 : GPS différentiel ### 9.1 Principe du GPS Le système GPS repose sur une constellation de 24 satellites émettant des signaux contenant l’heure de transmission et leur **éphéméride**. Un récepteur mesure le temps de propagation du signal depuis chaque satellite visible et en déduit une pseudo‑distance. Avec au moins quatre satellites, on résout le système d’équations pour obtenir les trois coordonnées spatiales (X, Y, Z) et l’horloge du récepteur. > Le principe de trilatération nécessite au moins quatre satellites pour résoudre les inconnues de position et d’horloge. ### 9.2 Sources d’erreur Les erreurs affectant le GPS se classent en deux catégories : - **Erreurs satellitaires** (éphémérides, horloge des satellites) – communes à tous les récepteurs. - **Erreurs locales** (délai ionosphérique, délai troposphérique, multipath, bruit du récepteur) – partiellement communes selon la distance entre récepteurs. ### 9.3 Fonctionnement du GPS différentiel Une **station de base** (référent) installée sur un point connu (coordonnées précises provenant d’un réseau géodésique tel que le Réseau Géodésique du Cameroun (RGC)) reçoit les mêmes signaux que le récepteur **rover** (mobile). Elle calcule, pour chaque satellite, la différence entre la pseudo‑distance mesurée et la vraie distance connue (correction). Ces corrections sont transmises en temps réel (via radio UHF, liaison cellulaire ou Internet) au rover, qui les applique à ses propres mesures. > La correction différentielle élimine essentiellement les erreurs d’éphéméride et d’horloge satellitaire, réduisant l’erreur résiduelle à quelques centimètres. ### 9.4 Mise en œuvre au Cameroun 1. **Choix de la station de base** : installer le récepteur sur un point du réseau géodésique national (ex. : point IGN‑Yaoundé 1). 2. **Configuration** : définir le taux d’échantillonnage (1 Hz), le masque d’élévation (15°) et le protocole de diffusion (RTCM 3.2). 3. **Transmission** : utiliser une radio VHF 25 km de portée ou un modem 3G/4G pour envoyer les corrections au rover. 4. **Réception au rover** : le récepteur applique les corrections en temps réel et affiche la position corrigée. 5. **Contrôle de qualité** : vérifier le PDOP (Position Dilution of Precision) < 3 et le nombre de satellites utilisés ≥ 5. ### 9.5 Exemple chiffré Supposons une station de base située aux coordonnées connues : X = ‑162 450 m, Y = ‑4 850 100 m, Z = 4 200 050 m (système WGS84, projection UTM 32N). Le rover mesure une pseudo‑distance au satellite PRN 08 de 20 250 000 m alors que la vraie distance, calculée à partir des coordonnées connues du satellite et de la base, est de 20 249 950 m. La correction pour ce satellite est donc ‑ 50 m (vraie − mesurée). Après application de l’ensemble des corrections (moyenne sur 8 satellites), la position du rover passe de (‑162 406,67 m, ‑4 850 056,67 m, 4 200 093,33 m) à (‑162 450 m, ‑4 850 100 m, 4 200 050 m), soit une erreur résiduelle de moins de 0,01 m. --- ## II. Leçon 10 : Les problèmes d’unités ### 10.1 Système d’unités SI Le Système International d’unités (SI) définit sept unités de base : mètre (m) pour la longueur, kilogramme (kg) pour la masse, seconde (s) pour le temps, ampère (A) pour le courant électrique, kelvin (K) pour la température thermodynamique, mole (mol) pour la quantité de substance et candela (cd) pour l’intensité lumineuse. En topographie, les unités dérivées les plus utilisées sont le mètre carré (m²) pour les surfaces et le mètre cube (m³) pour les volumes. > Le mètre est défini comme la distance parcourue par la lumière dans le vide pendant 1/299 792 458 seconde. ### 10.2 Conversions d’unités fréquentes | Grandeur | Unité SI | Unité usuelle | Facteur de conversion | |----------|----------|---------------|-----------------------| | Longueur | mètre (m) | pied (ft) | 1 m = 3,28084 ft | | Longueur | mètre (m) | degré (décimal) | dépend de la latitude : 1° ≈ 111 km à l’équateur | | Angle | radian (rad) | gon (grad) | 1 rad = 63,662 gon | | Angle | degré (°) | gon | 1° = 1,111… gon | | Surface | m² | hectare (ha) | 1 ha = 10 000 m² | | Pression | pascal (Pa) | bar | 1 bar = 100 000 Pa | ### 10.3 Erreurs d’unités fréquentes - Confondre **degrés décimaux** et **degrés‑minutes‑secondes** (DMS) lors de la saisie de coordonnées. - Utiliser le **gon** alors que le logiciel attend des **degrés** (ou inversement) sans appliquer le facteur 0,9. - Oublier de convertir les **pieds** en mètres lorsqu’on importe un DAO provenant d’un projet américain. - Mélanger **mètres carrés** et **hectares** dans le calcul d’une superficie de parcelle, entraînant une erreur de facteur 100. ### 10.4 Tableau de conversion rapide (à garder à portée de main) | De | Vers | Multiplier par | |----|------|----------------| | mètre | pied | 3,28084 | | pied | mètre | 0,3048 | | degré décimal | gon | 1,111111 | | gon | degré décimal | 0,9 | | hectare | m² | 10 000 | | m² | hectare | 0,0001 | | radian | degré | 57,2958 | | degré | radian | 0,0174533 | ### 10.5 Exemple d’application Un levé topographique fournit la longueur d’un tronçon de route : 1250 ft. Le projet exige la longueur en mètres pour le calcul du volume de remblai. Conversion : 1250 ft × 0,3048 m/ft = 381,0 m. Le volume de remblai prévu est de 0,8 m de hauteur sur une largeur moyenne de 6 m. Volume = Longueur × largeur × hauteur = 381,0 m × 6 m × 0,8 m = 1 828,8 m³. Si l’on avait oublié de convertir les pieds, on aurait obtenu 1250 ft × 6 ft × 0,8 ft = 6 000 ft³ ≈ 169,9 m³, soit une sous‑estimation de plus de 90 %. --- ## III. Leçon 11 : Géoréférencement d’un plan ### 11.1 Définition Le géoréférencement consiste à associer à chaque point d’une image raster ou d’un dessin vectoriel des coordonnées géographiques réelles (système de projection connu) afin que le document puisse être superposé à d’autres données spatiales dans un SIG. > Un plan non géoréférencé ne possède aucune information de localisation ; il ne peut être utilisé que comme simple illustration. ### 11.2 Méthodes de géoréférencement 1. **Points de contrôle terrestres (GCP)** : relever sur le terrain des points facilement identifiables (coin de bâtiment, borne géodésique) avec un récepteur GPS différentiel. 2. **Transformation affine** : supposée valable lorsque le plan subit uniquement des translations, rotations et changements d’échelle uniformes. Elle nécessite au moins 3 points de contrôle non alignés. 3. **Transformation polynomiale d’ordre 2** : utilisée lorsqu’il y a des déformations locales (étirement, courbure) ; nécessite au moins 6 points de contrôle. 4. **Transformation par caoutchouc (rubber‑sheeting)** : ajuste localement l’image en déplaçant chaque point selon un maillage de points de contrôle ; utile pour des scans de plans anciens fortement déformés. ### 11.3 Étapes de géoréférencement (procédure pas à pas) 1. **Préparation du plan** : numériser le document à 300 dpi minimum, enregistrer au format TIFF ou PNG sans compression destructive. 2. **Choix du système de projection** : déterminer la projection utilisée par les données de référence (souvent UTM 32N, WGS84 pour le Cameroun). 3. **Saisie des points de contrôle** : - Ouvrir le plan dans le SIG (QGIS, ArcGIS). - Activer l’outil de géoréférencement. - Cliquer sur un point connu du plan, entrer ses coordonnées X,Y (obtenues par GPS différentiel ou provenant d’une base de données géodésique). - Répéter pour au moins 3 points (affine) ou 6 points (polynomiale d’ordre 2). 4. **Choix du type de transformation** : sélectionner affine ou polynomiale selon la qualité du plan et le nombre de GCP. 5. **Calcul du résidu** : le SIG affiche l’erreur quadratique moyenne (RMSE) ; accepter si RMSE < 0,5 m pour un projet de voirie urbaine. 6. **Application de la transformation** : créer le raster géoréférencé (GeoTIFF) ou mettre à jour le shapefile vectoriel. 7. **Contrôle final** : superposer le plan géoréférencé à une couche de référence (orthophoto, tracé routier) et vérifier visuellement l’alignement. ### 11.4 Exemple de chantier au Cameroun Un bureau d’études doit intégrer le plan de lotissement d’un quartier de Douala (échelle 1/500) dans le SIG communal. - **Étape 1** : le plan est scanné à 400 dpi, enregistré en TIFF. - **Étape 2** : la projection choisie est UTM 32N, datum WGS84. - **Étape 3** : quatre points de contrôle sont relevés avec un récepteur GPS différentiel : - P1 : coin Sud‑Ouest du lot 1 → X = 578 000 m, Y = 451 800 m - P2 : coin Nord‑Est du lot 5 → X = 578 350 m, Y = 452 200 m - P3 : borne géodésique IGN‑Douala‑03 → X = 578 250 m, Y = 451 500 m - P4 : intersection rue principale / rue secondaire → X = 578 250 m, Y = 452 050 m - **Étape 4** : les mêmes points sont identifiés sur le plan numérique (en pixels) aux coordonnées suivantes : - P1 : Xpixel = 120, Ypixel = 340 - P2 : Xpixel = 480, Ypixel = 300 - P3 : Xpixel = 300, Ypixel = 180 - P4 : Xpixel = 420, Ypixel = 120 - **Étape 5** : transformation affine sélectionnée (le SIG fournit les paramètres a,b,c,d,e,f et un RMSE = 0,12 m). - **Étape 6** : le raster géoréférencé est exporté en GeoTIFF et ajouté au SIG communal. - **Résultat** : le lotissement peut désormais être croisé avec le réseau d’assainissement et le plan d’urbanisme pour vérifier les servitudes. ### 11.5 Utilisation de logiciels - **QGIS (gratuit)** : barre d’outils Géoréférencer → ajouter points → choisir transformation → démarrage. - **ArcGIS Pro** : fenêtre Géoréférencer → ajouter points de contrôle → choisir transformation → mise à jour. Dans les deux cas, il est possible d’enregistrer les points de contrôle sous forme de fichier CSV pour une réutilisation ultérieure. --- ## Synthèse du chapitre - **GPS différentiel** : correction en temps réel des erreurs satellitaires grâce à une station de base connue ; précision sub‑déci‑métrique requise pour les travaux de voirie, de cadastre et de suivi d’ouvrages. - **Problèmes d’unités** : maîtriser les conversions entre SI et unités usuelles (pied, degré, gon, hectare) ; éviter les confusions DMS/DD et les erreurs de facteur (ex. : 1 ha = 10 000 m²). - **Géoréférencement d’un plan** : associer des coordonnées réelles à une image ou un dessin via des points de contrôle terrestres et une transformation (affine, polynomiale ou rubber‑sheeting) ; contrôler la qualité par le RMSE et superposer aux références SIG. - **Procédure type** : préparation du relevé → acquisition des points de contrôle (GPS différentiel) → choix de la transformation → calcul et application → vérification finale. - **Outils** : récepteurs GPS différentiels (ex. : Trimble R10, Leica GS18), radios UHF ou liaisons cellulaires pour les corrections, SIG libres (QGIS) ou propriétaires (ArcGIS) pour le géoréférencement. --- ## Glossaire | Terme | Définition | |-------|------------| | GPS différentiel | Technique qui améliore la précision du GPS en transmettant les erreurs mesurées par une station de base connue à un récepteur mobile. | | Station de base | Récepteur GPS fixe installé sur un point de coordonnées connues, chargé de calculer les corrections. | | Rover | Récepteur GPS mobile qui reçoit les corrections de la station de base pour déterminer sa position en temps réel. | | RTCM | Protocole de transmission des corrections différentielles (Radio Technical Commission for Maritime services). | | PDOP | Position Dilution of Precision, indicateur de la qualité géométrique de la configuration des satellites utilisés. | | GCP (Ground Control Point) | Point de contrôle terrestre dont les coordonnées réelles sont connues, utilisé pour le géoréférencement. | | Transformation affine | Modèle mathématique qui traduit, rotate et échelle uniformément un objet géométrique (nécessite ≥ 3 GCP non alignés). | | RMSE (Root Mean Square Error) | Racine carrée de la moyenne des carrés des écarts entre les coordonnées estimées et les coordonnées de référence ; mesure de la précision du géoréférencement. | | UTM | Universal Transverse Mercator, système de projection cartographique divisant la Terre en zones de 6° de longitude. | | Datum | Modèle mathématique de la forme de la Terre servant de référence pour les coordonnées géodésiques (ex. : WGS84). | | Rubber‑sheeting | Méthode de déformation locale d’une image où chaque point est déplacé selon un maillage de points de contrôle. | | SIG | Système d’Information Géographique, logiciel permettant la saisie, le stockage, l’analyse et la visualisation de données géoréférencées. | --- ## Exercices d'application corrigés ### Exercice 1 – GPS différentiel **Situation** : Une équipe de topographie doit implanter un poteau électrique sur la route Yaoundé‑Bafia. La station de base est installée sur le point géodésique IGN‑Yaoundé‑02 (X = ‑162 450 m, Y = ‑4 850 100 m, Z = 4 200 050 m). Le rover mesure les pseudo‑distances suivantes (en mètres) pour trois satellites : | Satellite | Pseudo‑distance mesurée (m) | Vraie distance (calculée à partir des coordonnées du satellite et de la base) (m) | |-----------|----------------------------|---------------------------------------------------| | PRN 05 | 20 120 350 | 20 120 300 | | PRN 12 | 20 098 750 | 20 098 720 | | PRN 20 | 20 150 100 | 20 150 050 | **Question** : Calculer les corrections pour chaque satellite, puis déterminer la position corrigée du rover sachant que sa position brute (avant correction) est (‑162 406,67 m, ‑4 850 056,67 m, 4 200 093,33 m). **Corrigé** - Correction PRN 05 = Vraie – Mesurée = 20 120 300 – 20 120 350 = ‑ 50 m - Correction PRN 12 = 20 098 720 – 20 098 750 = ‑ 30 m - Correction PRN 20 = 20 150 050 – 20 150 100 = ‑ 50 m Moyenne des corrections = (‑50 ‑ 30 ‑ 50)/3 = ‑ 43,33 m (on applique la même correction aux trois coordonnées, hypothèse d’erreur commune). Position corrigée = position brute + correction = (‑162 406,67 ‑ 43,33, ‑4 850 056,67 ‑ 43,33, 4 200 093,33 ‑ 43,33) = (‑162 450 m, ‑4 850 100 m, 4 200 050 m). Cette position coïncide pratiquement avec les coordonnées connues de la base (erreur résiduelle < 0,01 m), montrant l’efficacité du différentiel. --- ### Exercice 2 – Conversion d’unités **Situation** : Un relevé fournit la superficie d’une parcelle agricole : 3,75 ha. Le projet d’irrigation nécessite la superficie en mètres carrés pour calculer le débit nécessaire (débit = 5 L s⁻¹ ha⁻¹). **Question** : a) Convertir la superficie en m². b) Calculer le débit total nécessaire en L s⁻¹. **Corrigé** a) 1 ha = 10 000 m² → 3,75 ha × 10 000 m²/ha = 37 500 m². b) Débit par hectare = 5 L s⁻¹ ha⁻¹ → Débit total = 3,75 ha × 5 L s⁻¹ ha⁻¹ = 18,75 L s⁻¹. --- ### Exercice 3 – Géoréférencement d’un plan **Situation** : Un plan papier d’un lotissement de Ngaoundéré (échelle 1/1000) doit être intégré au SIG communal. Quatre points de contrôle sont relevés avec un GPS différentiel : | Point | Coordonnées X (m) | Coordonnées Y (m) | |-------|-------------------|-------------------| | A | 578 000 | 451 800 | | B | 578 350 | 452 200 | | C | 578 250 | 451 500 | | D | 578 250 | 452 050 | Les mêmes points sont identifiés sur le plan numérique (en pixels) aux coordonnées suivantes : | Point | Xpixel | Ypixel | |-------|--------|--------| | A | 120 | 340 | | B | 480 | 300 | | C | 300 | 180 | | D | 420 | 120 | **Question** : Déterminer les paramètres de la transformation affine qui passe des coordonnées pixel aux coordonnées réelles (X = a·Xpixel + b·Ypixel + c, Y = d·Xpixel + e·Ypixel + f). **Corrigé** : Le SIG calcule les paramètres de la transformation affine (a,b,c,d,e,f) à partir des points de contrôle. Pour cet exemple, le RMSE obtenu est de 0,12 m, inférieur au seuil de 0,5 m exigé pour un projet de voirie urbaine, donc la transformation est acceptable. --- ## Erreurs fréquentes et règles de sécurité | N° | Erreur fréquente | Conséquence | Règle de prévention / bonne pratique | |----|------------------|-------------|--------------------------------------| | 1 | Oublier d’activer le masque d’élévation (angle < 15°) sur le récepteur GPS | Inclusion de satellites à faible élévation → augmentation du multipath et de l’erreur ionosphérique | Toujours régler le masque d’élévation à 15° ou plus avant le démarrage d’une session. | | 2 | Confondre degrés décimaux et degrés‑minutes‑secondes lors de la saisie de coordonnées | Décalage de plusieurs kilomètres voire dizaines de kilomètres | Vérifier l’unité affichée dans le logiciel ; convertir systématiquement en degrés décimaux avant toute opération. | | 3 | Utiliser un facteur de conversion incorrect (pied → mètre) lors du import d’un DAO | Erreur d’échelle sur tout le projet (ex. : voirie trop longue ou trop courte) | Conserver un tableau de conversion à portée de main ; vérifier une longueur connue (ex. : largeur d’une voie) après import. | | 4 | Ne pas réaliser de contrôle de qualité (RMSE) après géoréférencement | Superposition inexacte des couches → décisions d’urbanisme erronées | Exiger un RMSE < 0,5 m pour les projets de voirie urbaine ; refuser le travail si la valeur dépasse ce seuil. | | 5 | Négliger la vérification de la batterie du récepteur GPS différentiel en milieu isolé | Perte de liaison base‑rover → arrêt du travail | Toujours emporter une batterie de rechange et tester l’autonomie avant de quitter la base. | | 6 | Utiliser un radio UHF non homologué ou une puissance excessive | Interférences avec les services de communication, risque de sanction légale | Respecter les fréquences autorisées par l’ANRT (Agence Nationale des Réglementations des Télécommunications) et limiter la puissance à celle spécifiée dans le manuel du matériel. | --- *Fin du cours.*