systèmes de positionnement satellitaires (partie 1/2)

## Objectifs et prérequis **À la fin de ce chapitre, l’élève sera capable de :** - Définir un système de positionnement satellitaire (GNSS) et en citer les principaux constituants. - Expliquer le principe de mesure du temps de propagation et l’équation de base du pseudorange. - Décrire la configuration du GPS (orbites, constellation, segments) et les corrections nécessaires au positionnement. - Réaliser un calcul simple de position à partir de pseudoranges mesurés par un récepteur. - Identifier les applications professionnelles du GPS au Cameroun (topographie, agriculture de précision, suivi de chantiers, gestion des risques). - Appliquer les bonnes pratiques de sécurité et reconnaître les erreurs fréquentes en utilisation du GNSS. **Prérequis :** - Notions de base en physique (vitesse de la lumière, référentiels inertiels). - Connaissances en géodésie élémentaire (systèmes de coordonnées, ellipsoïde de référence). - Maîtrise du calcul matriciel de base (résolution d’un système linéaire 4×4). - Utilisation d’un récepteur GPS de terrain (mode lecture de coordonnées). --- ## Introduction contextualisée Au Cameroun, les projets d’aménagement du territoire, les travaux routiers (comme l’autoroute Yaoundé‑Douala) et les programmes d’agriculture de précision nécessitent une localisation fiable au centimètre près. Les systèmes de positionnement satellitaire, en particulier le GPS (Global Positioning System), offrent cette précision lorsqu’ils sont utilisés avec des corrections différentielles ou en temps réel (DGPS, RTK). Ce chapitre présente d’abord les concepts généraux des GNSS, puis se concentre sur le GPS, son fonctionnement, ses corrections et ses usages professionnels dans le contexte camerounais. --- ## I. Systèmes de positionnement satellitaires : généralités, présentation ### A. Définition et principe de base Un **système de positionnement satellitaire (GNSS)** est une constellation de satellites émettant des signaux radiofréquence permettant à un récepteur terrestre de déterminer sa position (latitude, longitude, altitude) et l’heure universelle en mesurant le temps de propagation des signaux. > **> Principe fondamental** : le pseudorange $\rho$ est égal à la vitesse de la lumière multipliée par le délai de propagation mesuré ($\rho = c\,\Delta t$). ### B. Historique et systèmes majeurs | Système | Pays/Opérateur | Année de mise en service opérationnelle | Nombre de satellites nominale | |---------|----------------|----------------------------------------|------------------------------| | GPS | États‑Unis (DoD) | 1995 | 24 | | GLONASS | Russie | 1995 | 24 | | Galileo | Union européenne | 2016 (services initiaux) | 24 | | BeiDou | Chine | 2020 (constellation complète) | 35 | > **> Fait à retenir** : le GPS reste le système le plus déployé mondialement et constitue la référence pour les applications civiles au Cameroun. ### C. Composants d’un système GNSS 1. **Segment spatial** : satellites sur orbite moyenne terrestre (≈20 200 km pour le GPS) portant des horloges atomiques et des émetteurs L‑band. 2. **Segment de contrôle** : stations de suivi au sol (ex. stations du segment de contrôle du GPS, du GLONASS, de Galileo, de BeiDou) qui mesurent les orbites et les horloges, puis envoient les corrections aux satellites. 3. **Segment utilisateur** : récepteurs GNSS (smartphones, récepteurs topographiques, modules embarqués) qui décodent les signaux et calculent la position. ### D. Applications au Cameroun - **Topographie et levé de plans** : création de cartes cadastrales, suivi des limites de parcelles agricoles. - **Agriculture de précision** : guidage de tracteurs pour semis, épandage d’engrais, réduction des intrants. - **Gestion des chantiers routiers** : contrôle du tracé, suivi du compactage, prévention des glissements de terrain. - **Navigation maritime et fluviale** : positionnement des navires sur le fleuve Sanaga et le lac Tchad. - **Gestion des catastrophes** : localisation rapide des zones inondées ou touchées par des glissements de terrain après des pluies abondantes. --- ## II. Systèmes de positionnement satellitaires : configuration, fonctionnement et utilisations du GPS ### A. Architecture du GPS - **Constellation** : 24 satellites opérationnels répartis sur 6 plans orbitaux inclinés à 55°, période orbitale de 11 h 58 min. - **Signaux** : deux fréquences civiles L1 (1575,42 MHz) et L2 (1227,60 MHz) ; code C/A (coarse/acquisition) sur L1 et code P (precision) sur L1 et L2. - **Horloges** : chaque satellite embarque une horloge à rubidium ou à césium, stabilisée à mieux que $1\times10^{-13}\ \text{s/s}$ (stabilité de l’ordre de $10^{-13}$). > **> Formule du période orbitale** : $T = 2\pi\sqrt{\frac{a^{3}}{\mu}}$ où $a$ est le demi‑grand axe et $\mu = 3,986\times10^{14}\ \text{m}^{3}\text{s}^{-2}$ (constante gravitationnelle terrestre). ### B. Principe de mesure : pseudorange Le récepteur mesure le temps de transmission $t_{tr}$ entre l’émission du signal par le satellite et sa réception. Le **pseudorange** $\rho$ est donné par : \[ \rho = c\,(t_{r} - t_{s}) + \text{erreurs} \] où $c = 299\,792\,458\ \text{m/s}$ est la vitesse de la lumière, $t_{r}$ l’heure de réception (horloge du récepteur) et $t_{s}$ l’heure d’émission (horloge du satellite). > **> À retenir** : le pseudorange diffère de la vraie distance géométrique du fait des erreurs d’horloge (satellite et récepteur), d’orbite, d’ionosphère, de troposphère et de relativité. ### C. Corrections nécessaires | Source d’erreur | Typique (m) | Modèle de correction | |-----------------|-------------|----------------------| | Erreur d’horloge du satellite | ±10 | Données d’horloge transmises dans le message de navigation | | Erreur d’orbite (éphéméride) | ±2,5 | Éphémérides diffusées ; corrections précises (IGS) en post‑traitement | | Délai ionosphérique | ±5 (jour) | Modèle Klobuchar (L1) ou mesure dual‑fréquence (L1‑L2) | | Délai troposphérique | ±2,3 | Modèle Saastamoinen ou modèle humide basé sur pression/température | | Effet relativiste (Einstein) | +38 µs/jour ≈ +7,2 km | Correction intégrée dans la fréquence d’émission du satellite | | Brouillage/multipath | variable | Antenne à gain élevé, traitement de signaux, sites ouverts | > **> Règle de sécurité** : en zone urbaine dense, privilégier un récepteur à antenne choke‑ring pour réduire le multipath. ### D. Procédure de positionnement 1. **Acquisition** : recherche des codes C/A des satellites visibles. 2. **Suivi** : verrouillage de la porteuse et mesure du pseudorange et de la phase. 3. **Résolution d’équation** : pour $n$ satellites, on obtient $n$ équations : \[ \rho_i = \sqrt{(x-x_i)^2+(y-y_i)^2+(z-z_i)^2} + c\,\delta t_r - c\,\delta t^s + \varepsilon_i \] où $(x,y,z)$ sont les coordonnées inconnues du récepteur, $(x_i,y_i,z_i)$ celles du satellite $i$, $\delta t_r$ l’erreur d’horloge du récepteur, $\delta t^s$ l’erreur d’horloge du satellite, $\varepsilon_i$ les erreurs résiduelles (ionosphère, troposphère, multipath, bruit). 4. **Solution** : linéarisation autour d’une position approximative et résolution par moindres carrés (généralement 4 satellites minimum pour résoudre les 4 inconnues $x,y,z,\delta t_r$). **Différentiel GPS (DGPS)** : une station de référence connue émet les corrections des pseudoranges en temps réel (via radio VHF ou liaison cellulaire). Le récepteur applique ces corrections, réduisant l’erreur à ≈0,5 m. **RTK (Real‑Time Kinematic)** : utilisation de la phase de porteuse avec double différence, permettant une précision centimétrique (1‑2 cm) sous condition de liaison radio stable entre base et rover. > **> À retenir** : le RTK nécessite une base fixe connue à moins de 20 km du rover pour garantir la corrélation des erreurs atmosphériques. ### E. Utilisations professionnelles du GPS au Cameroun | Domaine | Exemple concret | Précision visée | Équipement typique | |---------|----------------|----------------|--------------------| | Levés topographiques | Bornage d’un lot de 5 ha dans la région de l’Ouest | 1‑2 cm (RTK) | Récepteur GNSS double fréquence + antenne fixe | | Suivi de chantier routier | Contrôle du tracé de l’autoroute Yaoundé‑Douala (section 12 km) | ≈0,5 m (DGPS) | Rover GPS + liaison radio UHF | | Agriculture de précision | Guidage d’un tracteur pour semis de maïs dans le Nord‑Ouest | 2‑3 cm (RTK) | Tracteur équipé d’un récepteur RTK + écran de guidage | | Foresterie | Inventaire des essences dans la réserve du Dja | 1‑2 m (GPS autonome + post‑traitement PPP) | Récepteur handheld avec possibilité de corrections PPP ou SBAS si disponible | | Gestion des risques | Cartographie des zones inondables après les pluies de septembre 2023 | ≈1 m (DGPS) | Rover GPS + poste de référence temporaire | ### F. Exemple chiffré : calcul de position à partir de 4 pseudoranges **Données** (unités en mètres) : | Satellite | $x_i$ | $y_i$ | $z_i$ | Pseudorange $\rho_i$ | |-----------|-------|-------|-------|----------------------| | 1 | 20 000 000 | 0 | 0 | 20 000 002 | | 2 | 0 | 20 000 000 | 0 | 19 999 999 | | 3 | 0 | 0 | 20 000 000 | 20 000 003 | | 4 | -20 000 000 | 0 | 0 | 19 999 998 | On suppose une position approximative $(x_0,y_0,z_0) = (0,0,0)$ et une erreur d’horloge nulle. On linéarise chaque équation : \[ \rho_i \approx \rho_i^0 + \frac{(x_0-x_i)}{\rho_i^0}\Delta x + \frac{(y_0-y_i)}{\rho_i^0}\Delta y + \frac{(z_0-z_i)}{\rho_i^0}\Delta z + c\,\Delta t_r \] où $\rho_i^0 = \sqrt{x_i^2+y_i^2+z_i^2}$ est la distance géométrique à l’approximation. Calcul des $\rho_i^0$ : - $\rho_1^0 = \rho_2^0 = \rho_3^0 = \rho_4^0 = 20\,000\,000\ \text{m}$ On forme la matrice $A$ (4×4) dont les lignes sont $\big[-(x_i-x_0)/\rho_i^0,\ -(y_i-y_0)/\rho_i^0,\ -(z_i-z_0)/\rho_i^0,\ 1\big]$ et le vecteur $b = [\rho_i-\rho_i^0]$. Ainsi : \[ A = \begin{bmatrix} -1 & 0 & 0 & 1\\ 0 & -1 & 0 & 1\\ 0 & 0 & -1 & 1\\ 1 & 0 & 0 & 1 \end{bmatrix}, \qquad b = \begin{bmatrix} 2\\ -1\\ 3\\ -2 \end{bmatrix} \] En résolvant $A\,\Delta\mathbf{x}= \mathbf{b}$ on obtient : \[ \Delta x = -2\ \text{m},\quad \Delta y = 1\ \text{m},\quad \Delta z = -3\ \text{m},\quad c\,\Delta t_r = 0\ \text{m} \] Ainsi la position estimée du récepteur est : \[ \boxed{(x,y,z) = (-2\ \text{m},\ 1\ \text{m},\ -3\ \text{m})} \] et l’erreur d’horloge du récepteur vaut $\Delta t_r = 0 / c = 0\ \text{s}$. > **> Point clé** : cet exemple montre que, avec une bonne approximation initiale, la résolution du système linéaire fournit une position au mètre près ; l’utilisation de mesures de phase et de corrections différentielles permet d’atteindre le centimètre. --- ## Synthèse du chapitre - Les GNSS se composent d’un segment spatial (satellites), de contrôle (stations au sol) et utilisateur (récepteurs). - Le GPS, opérationnel depuis 1995, utilise 24 satellites sur des orbites moyennes terrestres à 20 200 km d’altitude. - Le principe de mesure repose sur le pseudorange : $\rho = c\Delta t + \text{erreurs}$. - Les erreurs principales (horloge, orbite, ionosphère, troposphère, relativité) doivent être modélisées ou corrigées (DGPS, RTK, modèles). - Le positionnement nécessite au moins 4 satellites pour résoudre les 4 inconnues $(x,y,z,\delta t_r)$. - Au Cameroun, le GPS est employé en topographie, agriculture de précision, suivi de chantiers routiers, foresterie et gestion des risques, avec des précisions allant du mètre (GPS autonome) au centimètre (RTK). - Les bonnes pratiques d’utilisation comprennent l’antenne adaptée, la vérification des corrections, la séparation de la base et du rover (<20 km en RTK) et le respect des consignes de sécurité sur les chantiers. --- ## Glossaire - **GNSS** : Global Navigation Satellite System, ensemble des systèmes de positionnement satellitaire (GPS, GLONASS, Galileo, BeiDou). - **Pseudorange** : mesure de distance entre satellite et récepteur incluant l’erreur d’horloge du récepteur. - **Segment spatial** : partie d’un GNSS constituée des satellites en orbite. - **Segment de contrôle** : réseau de stations au sol qui surveillent les satellites et envoient les corrections. - **Segment utilisateur** : récepteurs GNSS utilisés pour déterminer la position. - **DGPS** : Differential GPS, technique de correction en temps réel utilisant une station de référence connue. - **RTK** : Real‑Time Kinematic, méthode de positionnement centimétrique basée sur la phase de porteuse. - **Éphéméride** : données décrivant l’orbite d’un satellite à un instant donné. - **Horloge atomique** : dispositif de mesure du temps d’une stabilité exceptionnelle utilisé dans les satellites GNSS. - **Multipath** : phénomène où le signal atteint l’antenne après réflexion sur des objets proches, provoquant des erreurs de mesure. - **SBAS** : Satellite‑Based Augmentation System (ex. EGNOS, WAAS), fournissant des corrections d’intégrité et de précision via satellite. - **Relativité** : effet de la théorie d’Einstein sur les horloges des satellites (dérive de +38 µs/jour). - **Klobuchar** : modèle ionosphérique diffusé dans le message de navigation GPS pour estimer le délai ionosphérique. - **Saastamoinen** : modèle troposphérique utilisé pour corriger le délai dû à l’humidité et à la pression atmosphérique. --- ## Exercices d'application corrigés ### Exercice 1 – Calcul de pseudorange Un récepteur situé à Yaoundé (approx. 3°52′ N, 11°31′ E) reçoit le signal du satellite GPS PRN 08 à l’heure $t_r = 512\,345,678\,901\,234\ \text{s}$ (temps GPS). L’heure d’émission du satellite, tirée du message de navigation, est $t_s = 512\,345,678,900,000\ \text{s}$. **a)** Calculer le pseudorange $\rho$. **b)** En supposant que l’erreur d’horloge du satellite est nulle, quelle est la distance géométrique approximative entre le récepteur et le satellite ? **Corrigé** a) $\Delta t = t_r - t_s = 0,000,901,234\ \text{s} = 9,01234\times10^{-4}\ \text{s}$. $\rho = c\Delta t = 299\,792\,458 \times 9,01234\times10^{-4} \approx 270\,200\ \text{m}$. b) Puisque l’erreur d’horloge du satellite est négligeable, le pseudorange approximatif vaut la distance géométrique : $\approx 270,2\ \text{km}$. --- ### Exercice 2 – Positionnement différentiel (DGPS) Une station de référence située à Bertoua (coordonnées connues : $X_R = 1\,200\,000\ \text{m},\ Y_R = 4\,500\,000\ \text{m},\ Z_R = 1\,200\,000\ \text{m}$) mesure les pseudoranges aux mêmes quatre satellites que dans l’exemple du chapitre et obtient les corrections suivantes (en mètres) : $\Delta\rho_1 = +0,8$, $\Delta\rho_2 = -0,5$, $\Delta\rho_3 = +0,3$, $\Delta\rho_4 = -0,2$. Un rover situé à proximité mesure les pseudoranges bruts : $\rho_1^{raw}=21\,450\,000$, $\rho_2^{raw}=20\,980\,000$, $\rho_3^{raw}=22\,110\,000$, $\rho_4^{raw}=20\,760\,000$. **a)** Déterminer les pseudoranges corrigés. **b)** En utilisant la même linéarisation que dans l’exemple du chapitre, estimer la correction de position $(\Delta x,\Delta y,\Delta z)$ apportée par le DGPS (on peut supposer que l’erreur d’horloge du récepteur reste nulle). **Corrigé** a) $\rho_i^{corr} = \rho_i^{raw} + \Delta\rho_i$ : - $\rho_1^{corr}=21\,450\,000,8$ m - $\rho_2^{corr}=20\,979\,999,5$ m - $\rho_3^{corr}=22\,110\,000,3$ m - $\rho_4^{corr}=20\,759\,999,8$ m b) En recomputant les $\rho_i^0$ (distances géométriques à l’approximation $(0,0,0)$) on obtient les mêmes valeurs que précédemment (environ 20,0 Mm pour chaque satellite dans l’exemple simplifié). Les résidus $b_i = \rho_i^{corr} - \rho_i^0$ sont alors : - $b_1 = +0,8$ m - $b_2 = -0,5$ m - $b_3 = +0,3$ m - $b_4 = -0,2$ m En résolvant $A\,\Delta\mathbf{x}= \mathbf{b}$ (même matrice $A$ que précédemment) on trouve : $\Delta x \approx +0,4\ \text{m},\ \Delta y \approx -0,2\ \text{m},\ \Delta z \approx +0,1\ \text{m}$. Ainsi le DGPS corrige la position du rover de quelques décimètres dans chaque direction, montrant l’efficacité de la différentiation. --- ### Exercice 3 – Choix d’une solution RTK pour un levé topographique Un ingénieur doit réaliser le bornage d’une parcelle de 3 ha située à proximité de la ville de Bafoussam. La précision requise est de 2 cm. Il dispose d’un récepteur GNSS double fréquence (L1/L2) et d’une radio UHF pouvant transmettre à 5 km. **a)** Quels sont les éléments indispensables pour mettre en place une solution RTK fiable sur ce site ? **b)** Si la seule base disponible se trouve à 18 km du chantier, quelle alternative proposer pour atteindre la précision demandée ? **Corrigé** a) Pour du RTK il faut : - Une base fixe installée sur un point connu (coordonnées précises à l’échelle centimétrique, idéalement issue d’un réseau géodésique national). - Un lien de communication en temps réel (radio UHF, modem cellulaire ou liaison Internet) entre la base et le rover, avec une portée suffisante pour garantir la réception des corrections de phase. - Un récepteur rover capable de suivre les signaux L1 et L2 en mode double fréquence et de réaliser les doubles différences de phase. - Un environnement avec un masque d’horizon ouvert (angle d’élévation >15°) afin de minimiser le multipath et les pertes de signal. b) Avec une base à 18 km, la corrélation des erreurs ionosphériques et troposphériques devient insuffisante pour du RTK classique. Deux alternatives : 1. **Utiliser un réseau de référence permanente (CORS)** : si le Cameroun dispose d’un réseau de stations GNSS publiques (ex. stations de l’IGNG), le rover peut se connecter à la correction réseau (NTRIP) qui fournit des corrections interpolées, valables jusqu’à 50 km. 2. **Passer à une solution PPP (Precise Point Positioning)** avec corrections précises d’orbite et d’horloge (ex. service IGS) ; la précision atteinte est généralement de 3‑5 cm en statique et peut être améliorée à 2 cm avec des temps d’observation de 30 min et un filtrage adapté. Dans les deux cas, il faut vérifier la disponibilité du service (couverture NTRIP ou accès à Internet) et prévoir une batterie suffisante pour la durée de la mesure. --- ## Erreurs fréquentes et règles de sécurité | # | Erreur fréquente | Conséquence | Règle de prévention / bonne pratique | |---|------------------|-------------|--------------------------------------| | 1 | Négliger le masque d’horizon (antenne placée près d’un bâtiment ou d’un arbre) | Augmentation du multipath, pertes de lock, erreurs de plusieurs mètres | Installer l’antenne sur un trépied dégagé, angle d’élévation minimal de 15°, vérifier visibilité du ciel avec un outil de planification (ex. GPSTest). | | 2 | Utiliser une batterie faible pendant une mesure RTK en mouvement | Coupure du lien de données, perte de la correction, dérive de la position | Vérifier le niveau de charge avant le départ, prévoir une batterie de secours ou un panneau solaire portable pour les longues journées. | | 3 | Oublier d’appliquer les corrections ionosphériques (modèle Klobuchar) en zone équatoriale pendant les heures de forte activité solaire | Erreurs pouvant dépasser 5 m, surtout en période de tempête solaire | Activer le modèle ionosphérique du récepteur ou utiliser une correction dual‑fréquence (L1‑L2) ; consulter les prévisions d’indice Kp avant les campagnes critiques. | | 4 | Confondre les unités lors de la saisie des coordonnées (degrés décimaux vs degrés‑minutes‑secondes) | Décalage de plusieurs centaines de mètres à kilomètres | Toujours vérifier l’unité affichée dans le logiciel de saisie ; convertir systématiquement en degrés décimaux avant tout calcul. | | 5 | Ne pas sécuriser le périmètre autour de la station de base (ex. chantier routier) | Risque de choc avec engins de chantier, détérioration de l’antenne | Délimiter la zone de la base avec des rubans de signalisation et des panneaux « Zone de travail – Antenne GNSS – Ne pas toucher ». | | 6 | Utiliser un récepteur non homologué pour des travaux de topographie officielle | Résultats non conformes aux normes nationales (IGN, MINTP) pouvant entraîner le rejet du dossier | Sélectionner du matériel agréé par le MINTP/IGNG, conserver le certificat d’étalonnage et réaliser un contrôle de précision sur un point connu avant chaque campagne. | --- *Fin du cours.*