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

## Objectifs et prérequis **Objectifs** : à la fin de ce chapitre l’élève sera capable de - concevoir et mettre en œuvre une base de données spatiales adaptée à un besoin local (ex. suivi des réseaux d’eau) ; - initialiser un récepteur GPS de navigation, acquérir les satellites et évaluer la qualité de la position ; - expliquer le principe du GPS absolu, formuler l’équation de pseudo‑distance et identifier les principales sources d’erreur. **Prérequis** : notions de SIG (coordonnées géographiques, systèmes de projection), bases de données relationnelles (modèle SQL), notions élémentaires d’électricité et de propagation d’ondes radioélectriques. --- ## Introduction contextualisée Au Cameroun, la géomatique soutient les projets d’aménagement du territoire, la gestion des ressources naturelles et la planification urbaine. Les bases de données spatiales permettent de centraliser les informations topographiques, les réseaux d’infrastructure et les données environnementales. Le GPS, quant à lui, fournit la référence géodésique indispensable au relevé de terrain, à la navigation maritime le long du littoral atlantique et à la précision des levés cadastraux. Maîtriser l’implémentation d’une base de données spatiales et l’utilisation du GPS constitue donc une compétence clé pour les techniciens supérieurs en géomatique, très recherchée dans les bureaux d’études, les collectivités territoriales et les entreprises de travaux publics. --- ## I. Leçon 6 : Implémentation des bases de données spatiales ### 6.1 Concepts de base de données spatiales Une **base de données spatiales** (SDB) stocke des objets géométriques (points, lignes, polygones) associés à leurs attributs alphanumériques. Elle permet des requêtes spatiales (intersection, distance, voisinage) grâce à des index spécifiques (R‑tree, quadtree). > **Définition** : Une SDB est un SGBD étendu par des types de données géométriques et des opérateurs topologiques conformes aux normes OGC (Open Geospatial Consortium). ### 6.2 Modèles de données | Modèle | Représentation | Avantages | Inconvénients | |--------|----------------|-----------|---------------| | Vecteur | Points, lignes, polygones (coordonnées discrètes) | Précision élevée, bon pour les réseaux | Moins adapté aux phénomènes continus | | Raster | Grille de cellules (valeur par pixel) | Simple pour les images satellitaires, MNT | Redondance mémoire, perte de précision | | Topologie | Relations d’incidence (sommet‑arete‑face) | Intégrité géométrique, détection d’erreurs | Complexité de mise en place | ### 6.3 Systèmes de gestion de bases de données spatiales (SGBD‑S) Les principaux SGBD‑S utilisés au Cameroun sont : - **PostGIS** (extension de PostgreSQL) : open‑source, support complet des normes OGC, bonnes performances avec index GIST. - **Spatialite** (extension de SQLite) : léger, adapté aux projets de terrain sans serveur dédié. - **Oracle Spatial** : solution commerciale, haute disponibilité, utilisé dans les grandes administrations. ### 6.4 Étapes d’implémentation 1. **Analyse des besoins** : identifier les thématiques (hydrographie, voirie, cadastre) et les niveaux de précision requis. 2. **Conception logique** : modéliser les entités (ex. `Hydrant`, `Pipe`, `Valve`) et leurs attributs (diamètre, matériau, date de pose). 3. **Création du schéma** : écrire les requêtes SQL `CREATE TABLE` avec colonne `geometry` de type PostGIS. 4. **Importation des données** : utiliser `shp2pgsql` ou `ogr2ogr` pour convertir shapefiles ou GeoJSON en tables spatiales. 5. **Indexation spatiale** : créer un index GIST (`CREATE INDEX ON table USING GIST (geom);`). 6. **Validation** : vérifier les contraintes topologiques (ex. `ST_IsValid`) et exécuter des requêtes test (ex. recherche de tous les tuyaux situés à moins de 50 m d’un point d’incendie). ### 6.5 Exemple d’atelier : suivi du réseau d’eau potable à Yaoundé **Contexte** : la société camerounaise des eaux (SCE) souhaite gérer les interventions de maintenance sur son réseau de distribution. **Données** : shapefile contenant les tuyaux (ligne), les vannes (point) et les réservoirs (polygone). **Procédure** : 1. Charger le shapefile dans QGIS, vérifier le système de projection (UTM zone 32N, WGS84). 2. Exporter en PostGIS : `shp2pgsql -s 32632 -I réseau_eau.shp public.reseau | psql -U postgres -d geomatics`. 3. Créer une vue matérialisée pour les sections critiques : ```sql CREATE MATERIALIZED VIEW vue_critique AS SELECT * FROM reseau WHERE diamètre < 100 AND matériau = 'PVC'; ``` 4. Mettre en place un déclencheur qui enregistre la date de dernière intervention dans une table `maintenance` lorsqu’une vanne est mise à jour. **Résultat** : les techniciens peuvent interroger la base pour connaître immédiatement les tronçons à remplacer et planifier les interventions avec un gain de temps estimé à 30 %. > **Règle de métier** : toujours sauvegarder le schéma SQL et les fichiers de projection avant toute importation massive afin de garantir la traçabilité. --- ## II. Leçon 7 : Initialisation à l’utilisation du GPS de navigation ### 7.1 Principe du GPS Le système GPS se compose de 24 satellites (plus des réserves) répartis sur six plans orbitaux, émettant deux fréquences porteuses : L1 (1575,42 MHz) et L2 (1227,60 MHz). Le récepteur mesure le temps de propagation du signal depuis chaque satellite afin de déterminer sa distance (pseudo‑distance). ### 7.2 Étapes d’initialisation 1. **Mise sous tension** : le récepteur charge l’almanach (orbites approximatives) et l’éphéméride (orbites précises) stockés en mémoire. 2. **Recherche des satellites** : le récepteur scanne le ciel, identifie les satellites visibles grâce à leurs codes PRN. 3. **Calcul de la dilution de précision (PDOP)** : valeur inférieure à 5 indique une bonne géométrie. 4. **Application des corrections** : si disponible, le récepteur utilise le signal SBAS (EGNOS en Afrique) pour améliorer la précision à < 3 m. 5. **Obtention du fix** : lorsque au moins quatre satellites sont suivis avec un PDOP acceptable, le récepteur fournit une position 3D (latitude, longitude, altitude) et une heure précise. ### 7.3 Procédure pas à pas sur un récepteur Garmin eTrex 30x (modèle fréquent au Cameroun) | Étape | Action | Affichage attendu | |-------|--------|-------------------| | 1 | Allumer l’appareil (bouton Power) | Écran d’accueil, recherche satellites | | 2 | Attendre l’acquisition des premières satellites | Icône satellite clignotante, nombre de satellites > 4 | | 3 | Vérifier le PDOP | Menu → Setup → System → PDOP < 5 | | 4 | Activer EGNOS si disponible | Menu → Setup → System → SBAS → ON | | 5 | Confirmer le fix 3D | icône fixe apparaît, coordonnées affichées | | 6 | Sauvegarder le waypoint actuel | Mark Waypoint → Nom → OK | ### 7.4 Précautions de sécurité et d’hygiène - Éviter les zones urbaines denses (immeubles hauts) qui provoquent du masquage et du multipath. - Vérifier l’état de la batterie avant toute sortie terrain ; une faible tension peut entraîner une perte de synchronisation. - Utiliser une antenne externe magnétique sur le véhicule lorsqu’on travaille en milieu forestier pour améliorer la réception. - Ne pas laisser le récepteur exposé à la pluie directe ; utiliser une housse étanche. ### 7.5 Exemple chiffré : temps de première fixation (TTFF) Le TTFF dépend du nombre de satellites visibles (N) et de la qualité de l’almanach. Une approximation empirique : \[ \text{TTFF (s)} \approx \frac{30}{N} + 10 \] Si le récepteur voit 6 satellites : \[ \text{TTFF} \approx \frac{30}{6}+10 = 5+10 = 15\text{ s} \] > **Définition** : TTFF (Time To First Fix) est le délai entre la mise sous tension du récepteur et l’obtention d’une première position valide. --- ## III. Leçon 8 : GPS absolu, mesure de pseudo‑distance ### 8.1 Concept de pseudo‑distance La pseudo‑distance ρ est la mesure du temps de propagation du signal multiplié par la vitesse de la lumière dans le vide (c ≈ 299 792 458 m s⁻¹). Elle comprend l’erreur d’horloge du récepteur et les retardations atmosphériques. ### 8.2 Équation de mesure Pour un satellite s : \[ \rho_s = \bigl\| \mathbf{X} - \mathbf{X}_s \bigr\| + c\,\delta t + I_s + T_s + \varepsilon_s \] où : - \(\mathbf{X}\) = position inconnue du récepteur (ECEF) - \(\mathbf{X}_s\) = position connue du satellite (éphéméride) - \(\delta t\) = erreur d’horloge du récepteur (secondes) - \(I_s\) = retard ionosphérique (mètres) - \(T_s\) = retard troposphérique (mètres) - \(\varepsilon_s\) = bruit de mesure et multipath (mètres) > **Formule fondamentale** : la pseudo‑distance est la somme de la distance géométrique vraie, du terme d’horloge et des corrections atmosphériques. ### 8.3 Sources d’erreur principales | Source | Effet typique (mètres) | Méthode d’atténuation | |--------|------------------------|-----------------------| | Erreur d’horloge du récepteur | 10 – 30 (équivalent à 30‑90 ns) | Estimation simultanée de \(\delta t\) dans la résolution | | Retard ionosphérique | 5 – 15 (jour) | Modèle Klobuchar ou utilisation de fréquences duales (L1/L2) | | Retard troposphérique | 2 – 4 | Modèle Saastamoinen ou modèle météorologique local | | Multipath (réflexion) | 0,5 – 2 | Antenne à gain élevé, traitement de signal (smoothing) | | Erreur d’orbite satellite | 0,5 – 1,5 | Éphémérides précises (IGS) ou corrections DGPS | ### 8.4 Méthode de résolution On résout le système non linéaire d’équations de pseudo‑distance par la méthode des moindres carrés pondérés (WLS). Après linéarisation autour d’une position approximative \(\mathbf{X}_0\), on obtient : \[ \Delta \boldsymbol{\rho} = \mathbf{H}\,\Delta \mathbf{X} + \mathbf{v} \] où \(\mathbf{H}\) est la matrice de géométrie (lignes de vue unitaires) et \(\mathbf{v}\) le vecteur des résidus. La solution est : \[ \Delta \mathbf{X} = (\mathbf{H}^T \mathbf{W} \mathbf{H})^{-1} \mathbf{H}^T \mathbf{W} \Delta \boldsymbol{\rho} \] \(\mathbf{W}\) étant la matrice de poids liée à la variance des mesures. Un filtrage de Kalman peut être appliqué pour suivre en temps réel la position et l’horloge. ### 8.5 Exemple de calcul de pseudo‑distance Un récepteur reçoit un signal du satellite PRN 08 avec un délai de propagation mesuré de **70 ms**. \[ \rho = c \times \Delta t = 299\,792\,458 \times 0,070 = 20\,985\,472,06\ \text{m} \approx 20\,985\ \text{km} \] Si la distance géométrique vraie entre le récepteur et le satellite est de 20 200 km, l’excédent de 785 km provient de l’erreur d’horloge du récepteur (\(\delta t\)) et des retards atmosphériques. > **Encadré à retenir** : la pseudo‑distance mesure le temps de propagation du signal GPS multiplié par la vitesse de la lumière ; elle n’est pas une distance vraie tant que les corrections d’horloge et atmosphériques ne sont pas appliquées. --- ## Synthèse du chapitre - **Implémentation des bases de données spatiales** : conception logique, schéma SQL/PostGIS, importation shapefile, index GIST, requêtes spatiales. - **Initialisation du GPS de navigation** : acquisition almanach, recherche satellites, contrôle PDOP, correction SBAS, procédure pas à pas sur récepteur Garmin eTrex 30x. - **GPS absolu et pseudo‑distance** : équation de mesure, sources d’erreur (horloge, ionosphère, troposphère, multipath, orbite), résolution par moindres carrés pondérés, exemple de calcul. - **Compétences professionnelles** : capacité à mettre en place une base de données géospatiale pour un projet d’infrastructure, à préparer et utiliser un récepteur GPS en terrain camerounais, à interpréter et corriger les mesures de pseudo‑distance pour obtenir une position précise. --- ## Glossaire (6 termes) - **Base de données spatiales (SDB)** : SGBD étendu par des types de données géométriques et des opérateurs topologiques conformes aux normes OGC. - **Index GIST** : méthode d’indexation équilibrée permettant des recherches rapides sur les objets géométriques (intersection, distance). - **PDOP (Position Dilution of Precision)** : facteur géométrique qui quantifie l’effet de la configuration satellitaire sur la précision de la position. - **Pseudo‑distance (ρ)** : mesure du temps de propagation du signal GPS multiplié par la vitesse de la lumière, incluant les erreurs d’horloge et atmosphériques. - **SBAS (Satellite Based Augmentation System)** : système d’amélioration du signal GPS (ex. EGNOS) fournissant des corrections d’orbite et d’horloge pour augmenter la précision. - **WLS (Weighted Least Squares)** : méthode d’estimation qui minimise la somme pondérée des carrés des résidus, utilisée pour résoudre les équations de pseudo‑distance. --- ## Exercices d'application corrigés ### Exercice 1 – Création d’une base de données pour le suivi des pistes rurales **Consigne** : Vous devez concevoir une base de données PostGIS pour inventorier les pistes rurales d’un département (attributs : `id_piste`, `nom`, `longueur_km`, `revêtement`, `géométrie`). 1. Écrire la requête SQL de création de la table `pistes_rurales` avec une colonne géométrique de type `LINESTRING`. 2. Indiquer la commande permettant de créer un index GIST sur la colonne géométrique. 3. Donner la requête permettant de sélectionner toutes les pistes dont le revêtement est « laterite » et dont la longueur dépasse 5 km. **Corrigé** : 1. ```sql CREATE TABLE pistes_rurales ( id_piste SERIAL PRIMARY KEY, nom VARCHAR(100), longueur_km NUMERIC(5,2), revêtement VARCHAR(30), geom GEOMETRY(LINESTRING, 32632) -- UTM zone 32N ); ``` 2. ```sql CREATE INDEX idx_pistes_geom ON pistes_rurales USING GIST (geom); ``` 3. ```sql SELECT * FROM pistes_rurales WHERE revêtement = 'laterite' AND longueur_km > 5; ``` ### Exercice 2 – Évaluation de la précision d’un récepteur GPS **Consigne** : Un récepteur Garmin eTrex 30x affiche un PDOP de 4,2 après l’acquisition de 7 satellites. En utilisant l’approximation \(\sigma_{pos} \approx PDOP \times \sigma_{URE}\) où \(\sigma_{URE}=2,5\) m (erreur moyenne de l’orbite et de l’horloge satellite), calculez l’erreur théorique de position en mètres. **Corrigé** : \[ \sigma_{pos} = 4,2 \times 2,5 = 10,5\ \text{m} \] L’erreur théorique de position est donc d’environ **10,5 m**. --- ## Erreurs fréquentes et règles de sécurité | N° | Erreur fréquente | Conséquence | Règle de prévention / bonne pratique | |----|------------------|-------------|--------------------------------------| | 1 | Oublier de spécifier le SRID lors de la création de la colonne géométrique | Les coordonnées sont interprétées dans un système de projection inconnu, entraînant des décalages de plusieurs kilomètres | Toujours indiquer le SRID correspondant (ex. 32632 pour UTM zone 32N) lors de `CREATE TABLE … GEOMETRY(..., SRID)`. | | 2 | Laisser le récepteur GPS en mode « économie d’énergie » pendant l’acquisition initiale | Le récepteur ne suit pas suffisamment de satellites, le fix reste en 2D ou échoue | Désactiver le mode économie d’énergie pendant les premières minutes d’utilisation ou jusqu’à l’obtention d’un fix 3D. | | 3 | Négliger la vérification de l’état de la batterie avant une sortie terrain | Coupure soudaine du récepteur, perte de données de suivi | Mesurer la tension de la batterie avec un multimètre ; remplacer ou recharger si < 3,7 V pour un Li‑ion standard. | | 4 | Utiliser une antenne interne en milieu urbain dense sans tenir compte du masquage | Augmentation du multipath et du bruit, dégradation de la précision | Privilégier une antenne externe à gain élevé placée sur le toit du véhicule ou sur un trépied lorsqu’on travaille en ville. | | 5 | Appliquer directement les pseudo‑distances brutes sans correction d’horloge | Position erronée de plusieurs dizaines à centaines de mètres | Effectuer toujours une étape de résolution (moindres carrés ou filtre de Kalman) qui estime \(\delta t\) avant d’interpréter les coordonnées. | | 6 | Ne pas sauvegarder le fichier de projection (.prj) lors de l’export d’un shapefile | Lors de la réimportation, les données peuvent être décalées ou mal alignées | Toujours conserver le fichier .prj associé au shapefile ; lors de l’importation dans PostGIS, spécifier explicitement le SRID. | --- *Fin du cours.*