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

## Objectifs et prérequis **Objectifs** : À la fin de ce chapitre, l’élève sera capable de : - Identifier et décrire les principaux outils et documents utilisés en télédétection (satellites, capteurs, images orthorectifiées, métadonnées). - Expliquer la notion de base de données géographique (BDG) et de base de données spatiale (BDS), leurs modèles de données et leurs formats courants. - Concevoir pas à pas une BDS adaptée à un besoin professionnel camerounais (ex. : gestion du parcellaire agricole, suivi de la couverture forestière) en suivant les étapes de modélisation conceptuelle, logique et physique. - Réaliser un petit projet de création d’une BDS avec un SIG libre (QGIS) et produire les documents de spécification associés. **Prérequis** : - Notions de base en géographie (systèmes de coordonnées, projections, échelle). - Familiarité avec l’utilisation d’un ordinateur (système d’exploitation Windows/Linux, gestion de fichiers). - Connaissances élémentaires en informatique (notions de fichier, de tableur). --- ## Introduction contextualisée La géomatique regroupe l’ensemble des techniques permettant d’acquérir, de stocker, d’analyser et de diffuser des informations géoréférencées. Au Cameroun, les défis liés à la planification urbaine (Yaoundé, Douala), à la gestion des ressources naturelles (forêts du Congo, bassin du Logone) et à la sécurité alimentaire (suivi des rendements agricoles) imposent une maîtrise des outils de télédétection et des bases de données spatiales. Ce chapitre introduit les fondements indispensables pour passer de l’image brute à une information exploitable dans un SIG professionnel. --- ## I. Outils et documents de la télédétection ### 1. Définition et champ d’application > La télédétection est l’acquisition d’informations sur un objet ou une surface sans contact physique, généralement par des capteurs embarqués sur des plateformes aériennes ou spatiales. ### 2. Plateformes et capteurs principaux | Plateforme | Altitude typique | Exemple de mission camerounaise | Type de capteur | |------------|------------------|--------------------------------|-----------------| | Satellite optique | 600‑800 km | Landsat 8 (NASA/USGS) | Multispectral (VNIR‑SWIR) | | Satellite radar | 500‑700 km | Sentinel‑1A (ESA) | SAR (C‑band) | | Drone (UAV) | 50‑150 m | Missions de cartographie parcellaire dans le Nord‑Ouest | RGB, multispectral, LiDAR léger | | Avion léger | 1 000‑3 000 m | Surveys aériens de l’IGN Cameroun (Institut Géographique National du Cameroun) | Photogrammétrie, LiDAR | ### 3. Niveaux de traitement des images - **Niveau 0** : données brutes (DN – Digital Numbers). - **Niveau 1** : correction radiométrique et géométrique de base (géoréférencement). - **Niveau 2** : orthorectification et produits dérivés (indices de végétation, classification d’usage des sols). > Un produit de niveau 2 couramment utilisé au Cameroun est l’indice NDVI = (NIR − RED)/(NIR + RED), permettant de suivre la vigueur végétale des cultures de maïs et de cacao. ### 4. Documents associés - **Métadonnées** : informations sur la date d’acquisition, le capteur, la résolution spatiale (ex. : 30 m pour Landsat), la projection (UTM zone 32N, WGS84). - **Rapport de qualité** : taux de couverture nuageuse, erreur géométrique (RMSE). - **Catalogue d’images** : base de recherche (ex. : USGS EarthExplorer, Copernicus Open Access Hub). ### 5. Exemple d’atelier professionnel Dans un bureau d’études de Yaoundé, un technicien télécharge une scène Landsat 8 du 15 mai 2024 couvrant la région du Littoral, applique une correction atmosphérique (méthode DOS), calcule l’NDVI, puis classe les pixels en cinq catégories d’usage des sols (forêt, savane, culture, eau, urbain) à l’aide d’un algorithme de maximum de vraisemblance. Le résultat est exporté en GeoTIFF et intégré à la BDS du projet de suivi de la déforestation. --- ## II. Base de données géographiques et spatiales ### 1. Concepts de base > Une base de données géographique (BDG) stocke des données attributaires liées à des entités géographiques (points, lignes, polygones). > Une base de données spatiale (BDS) étend la BDG avec des capacités de requête topologique (intersection, distance, containment) et d’indexation spatiale (R‑tree, Quadtree). ### 2. Modèles de données | Modèle | Représentation | Avantages | Inconvénients | |--------|----------------|-----------|---------------| | Vecteur | Points, lignes, polygones (coordonnées XYZ) | Précision géométrique, bon pour les limites administratives | Moins adapté aux phénomènes continus | | Raster | Grille de cellules (valeur attributaire par cellule) | Adapté aux images satellitaires, modèles numériques de terrain | Redondance de stockage, perte de détail vectoriel | | Objet‑orienté | Entités avec méthodes et héritage | Modélisation complexe (réseaux, 3D) | Nécessite un SGBD objet‑orienté spécialisé | ### 3. Formats d’échange courants - **Shapefile (.shp)** : format ESRI, limité à 2 Go, champs d’attributs en dBase. - **GeoPackage (.gpkg)** : conteneur SQLite, supporte vecteur et raster, norme OGC. - **File Geodatabase (.gdb)** : ESRI, performant, propriétaire. - **PostGIS** : extension spatiale de PostgreSQL, idéal pour les applications web et les gros volumes. ### 4. Exemple de schéma attributaire pour une parcelle agricole | Champ | Type | Description | Unité | |-------|------|-------------|-------| | `ID_Parcelle` | Integer (PK) | Identifiant unique | – | | `Culture` | Text | Type de culture (maïs, cacao, etc.) | – | | `Surface` | Real | Aire de la parcelle | ha | | `Date_Semis` | Date | Date de plantation | – | | `Rendement` | Real | Production estimée | t/ha | | `Geom` | Geometry (Polygon) | Limite de la parcelle | – | > La contrainte d’intégrité `Surface > 0` garantit qu’aucune parcelle ne possède une aire nulle ou négative. ### 5. Notions de normalisation et d’indexation - **Première forme normale (1FN)** : chaque cellule contient une valeur atomique. - **Deuxième forme normale (2FN)** : élimination des dépendances partielles par rapport à la clé primaire. - **Index spatial R‑tree** : structure d’arbre qui regroupe les enveloppes englobantes minima (MBR) pour accélérer les recherches d’intersection. --- ## III. Conception et Implémentation d’une base donnée spatiale ### 1. Méthodologie de conception (cycle de vie) 1. **Analyse des besoins** : recueil des exigences fonctionnelles (ex. : suivi des parcelles agricoles, génération de rapports de rendement). 2. **Modélisation conceptuelle** : diagramme Entité‑Relation (ER) ou UML, identification des entités (Parcelle, Culture, Exploitant) et de leurs relations. 3. **Modélisation logique** : traduction en schéma relationnel, définition des clés primaires/étrangères, normalisation. 4. **Modélisation physique** : choix du SGBD (PostGIS), définition des types de données géométriques, création des index spatiaux. 5. **Implémentation** : écriture des scripts SQL de création de tables, insertion de données de test. 6. **Validation** : requêtes de contrôle (ex. : `SELECT ST_Area(geom) FROM parcelle WHERE culture='maïs';`). 7. **Documentation** : dictionnaire de données, métadonnées ISO 19115, manuel d’utilisation. ### 2. Exemple complet : création d’une BDS pour le suivi du parcellaire agricole dans la région de l’Ouest #### Étape 1 – Script SQL de création (PostGIS) ```sql -- Activation de l'extension PostGIS CREATE EXTENSION IF NOT EXISTS postgis; -- Table des exploitants CREATE TABLE exploitant ( id_exploitant SERIAL PRIMARY KEY, nom TEXT NOT NULL, contact TEXT ); -- Table des cultures CREATE TABLE culture ( id_culture SERIAL PRIMARY KEY, libelle TEXT NOT NULL UNIQUE ); -- Table des parcelles CREATE TABLE parcelle ( id_parcelle SERIAL PRIMARY KEY, id_exploitant INTEGER REFERENCES exploitant(id_exploitant), id_culture INTEGER REFERENCES culture(id_culture), surface REAL CHECK (surface > 0), date_semis DATE, rendement REAL, geom GEOMETRY(Polygon, 32632) -- SRID UTM zone 32N, WGS84 ); -- Index spatial sur la géométrie CREATE INDEX idx_parcelle_geom ON parcelle USING GIST (geom); ``` #### Étape 2 – Insertion de données de référence ```sql INSERT INTO exploitant (nom, contact) VALUES ('Jean Kamga', 'jean.kamga@email.cm'), ('Aïssa Nkoum', 'aissa.nkoum@email.cm'); INSERT INTO culture (libelle) VALUES ('Maïs'), ('Cacao'); INSERT INTO parcelle (id_exploitant, id_culture, surface, date_semis, rendement, geom) VALUES (1, 1, 1.5, '2024-03-10', 3.2, ST_GeomFromText('POLYGON((610900 555555, 666350 555555, 666350 611111, 610900 611111, 610900 555555))', 32632)), (2, 2, 0.8, '2024-02-20', 0.6, ST_GeomFromText('POLYGON((633032 577777, 688462 577777, 688462 633333, 633032 633333, 633032 577777))', 32632)); ``` #### Étape 3 – Requête d’analyse > Calcul du rendement moyen par culture pour les parcelles supérieures à 1 ha. ```sql SELECT c.libelle, AVG(p.rendement) AS rendement_moyen FROM parcelle p JOIN culture c ON p.id_culture = c.id_culture WHERE p.surface > 1.0 GROUP BY c.libelle; ``` #### Étape 4 – Export et partage - Export du tableau attributaire en CSV (`\copy (SELECT * FROM parcelle) TO 'parcelles.csv' CSV HEADER;`). - Export de la couche géographique en GeoPackage (`ogr2ogr -f GPKG parcelles.gpkg PG:"host=localhost dbname=geodb user=postgres" -sql "SELECT * FROM parcelle"`). - Génération d’un fichier de métadonnées ISO 19115 à l’aide de l’outil `GeoNetwork` ou `MDEditor`. ### 3. Bonnes pratiques de sécurité et d’hygiène des données - Sauvegarde quotidienne du cluster PostgreSQL (`pg_dump -Fc geodb > geodb.backup`). - Contrôle d’accès basé sur les rôles (`GRANT SELECT ON parcelle TO lecteur;`). - Journalisation des modifications (`CREATE TRIGGER log_parcelle AFTER INSERT OR UPDATE ON parcelle FOR EACH ROW EXECUTE FUNCTION audit_log();`). - Vérification de l’intégrité topologique avec `ST_IsValid(geom)` avant toute insertion massive. --- ## Synthèse du chapitre - **Télédétection** : acquisition d’images satellitaires, aériennes ou drones ; niveaux de traitement (0‑2) ; métadonnées essentielles (date, résolution, projection). - **Base de données géographique** : stockage d’entités géométriques avec attributs ; modèles vecteur/raster/objet‑orienté ; formats d’échange (Shapefile, GeoPackage, PostGIS). - **Base de données spatiale** : ajout de capacités topologiques et d’index spatial (R‑tree) ; workflow de conception (analyse → conceptuel → logique → physique → implémentation → validation → documentation). - **Exemple camerounais** : suivi du parcellaire agricole en région de l’Ouest avec PostGIS, création de tables, insertion de données, requêtes de rendement, export en GeoPackage et métadonnées ISO 19115. - **Sécurité** : sauvegardes, contrôle d’accès, journalisation, validation topologique. > **À retenir** : la résolution spatiale d’une image Landsat 8 est de 30 m × 30 m. > **À retenir** : un indice NDVI varie entre –1 et +1 ; valeurs proches de +1 indiquent une végétation dense. > **À retenir** : la clé primaire d’une table doit être unique et non nulle. > **À retenir** : l’index GIST sur une colonne `geom` accélère les recherches d’intersection dans PostGIS. --- ## Glossaire | Terme | Définition | |-------|------------| | Télédétection | Acquisition d’informations sur la Terre sans contact physique, via des capteurs embarqués. | | Orthorectification | Processus de correction géométrique d’une image pour éliminer les distorsions dues au relief et à l’inclinaison du capteur. | | NDVI (Normalized Difference Vegetation Index) | Indice spectral = (NIR − RED)/(NIR + RED) utilisé pour estimer la vigueur végétale. | | Base de données géographique (BDG) | Système organisé de stockage de données attributaires liées à des entités géoréférencées. | | Base de données spatiale (BDS) | BDG enrichie de fonctions de requête topologique et d’index spatial. | | Shapefile | Format d’échange vectoriel propriétaire d’ESRI, composé de plusieurs fichiers (.shp, .shx, .dbf, .prj). | | GeoPackage | Conteneur SQLite ouvert, norme OGC, pouvant stocker vecteur et raster. | | PostGIS | Extension spatiale de PostgreSQL permettant le stockage, l’indexation et l’analyse de données géométriques. | | SRID (Spatial Reference System Identifier) | Code identifiant un système de coordonnées de projection (ex. : 32632 pour UTM zone 32N, WGS84). | | Index R‑tree | Structure d’arbre qui regroupe les enveloppes englobantes minima (MBR) pour accélérer les recherches spatiales. | | Métadonnées ISO 19115 | Norme internationale décrivant les caractéristiques des données géographiques (titre, date, résolution, ligne de base, etc.). | | Géotraitement | Ensemble d’opérations (buffer, intersection, union, etc.) appliquées aux données géométriques dans un SIG. | | Normalisation (1FN, 2FN, 3FN) | Processus d’organisation des tables d’une base relationnelle pour éliminer les redondances et assurer l’intégrité des données. | | Gestion des droits (RBAC) | Modèle de contrôle d’accès où les permissions sont associées à des rôles plutôt qu’à des utilisateurs individuels. | --- ## Exercices d'application corrigés ### Exercice 1 – Calcul d’NDVI à partir de valeurs de réflectance **Données** : Une parcelle de maïs observée par Sentinel‑2 présente les réflectances suivantes (en pourcentage) : - Bande rouge (B4) = 12 % - Bande proche infrarouge (B8) = 45 % **Question** : Calculer l’indice NDVI et interpréter le résultat. **Corrigé** NDVI = (NIR − RED)/(NIR + RED) Convertir les pourcentages en valeurs décimales : RED = 0,12 ; NIR = 0,45 NDVI = (0,45 − 0,12)/(0,45 + 0,12) = 0,33/0,57 ≈ 0,579 Interprétation : Un NDVI de 0,58 indique une végétation saine et relativement dense, typique d’un maïs en bonne phase de croissance. --- ### Exercice 2 – Création d’une table parcellaire et requête de superficie totale **Contexte** : Vous devez concevoir une BDS pour le suivi des parcelles de cacao dans la région du Sud‑Ouest. **Étapes** : 1. Écrire le script SQL de création de la table `parcelle_cacao` avec les champs suivants : `id_parcelle` (clé primaire), `id_exploitant` (clé étrangère vers table `exploitant`), `surface` (real, > 0), `date_plantation` (date), `rendement` (real), `geom` (geometry(Polygon, 32632)). 2. Créer un index spatial GIST sur `geom`. 3. Insérer deux parcelles : - Parcelle A : superficie 0,9 ha, date plantation 2023‑11‑15, rendement 0,4 t/ha, géométrie POLYGON((10.0 4.0,10.5 4.0,10.5 4.5,10.0 4.5,10.0 4.0)). - Parcelle B : superficie 1,3 ha, date plantation 2023‑10‑01, rendement 0,6 t/ha, géométrie POLYGON((10.2 4.2,10.7 4.2,10.7 4.7,10.2 4.7,10.2 4.2)). 4. Réaliser une requête qui retourne la superficie totale des parcelles de cacao. **Corrigé** ```sql -- 1. Création de la table CREATE TABLE exploitant ( id_exploitant SERIAL PRIMARY KEY, nom TEXT NOT NULL ); CREATE TABLE parcelle_cacao ( id_parcelle SERIAL PRIMARY KEY, id_exploitant INTEGER REFERENCES exploitant(id_exploitant), surface REAL CHECK (surface > 0), date_plantation DATE, rendement REAL, geom GEOMETRY(Polygon, 32632) ); -- 2. Index spatial CREATE INDEX idx_parcelle_cacao_geom ON parcelle_cacao USING GIST (geom); -- 3. Insertion des données de référence INSERT INTO exploitant (nom) VALUES ('Samuel Mvondo'), ('Madeleine Tchinda'); INSERT INTO parcelle_cacao (id_exploitant, surface, date_plantation, rendement, geom) VALUES (1, 0.9, '2023-11-15', 0.4, ST_GeomFromText('POLYGON((611050 444444, 666575 444444, 666575 500000, 611050 500000, 611050 444444))', 32632)), (2, 1.3, '2023-10-01', 0.6, ST_GeomFromText('POLYGON((633200 466666, 688700 466666, 688700 522222, 633200 522222, 633200 466666))', 32632)); -- 4. Requête de superficie totale SELECT SUM(surface) AS superficie_totale_ha FROM parcelle_cacao; ``` **Résultat attendu** : `superficie_totale_ha = 2.2` ha. --- ## Erreurs fréquentes et règles de sécurité | # | Erreur fréquente | Conséquence | Règle de prévention / bonne pratique | |---|------------------|-------------|--------------------------------------| | 1 | Oublier de définir un SRID lors de la création d’une colonne `geom`. | Les géométries sont considérées comme inconnues ; les opérations de distance ou d’intersection retournent NULL ou provoquent des erreurs. | Toujours spécifier le SRID (ex. : `GEOMETRY(Polygon, 32632)`) et vérifier avec `SELECT Find_SRID('public','parcelle','geom');`. | | 2 | Insérer une valeur de surface négative ou nulle. | Fausse estimation de superficie, biais dans les rapports de rendement. | Ajouter une contrainte `CHECK (surface > 0)` lors de la création de la table. | | 3 | Utiliser un shapefile pour stocker plus de 2 Go de données. | Échec de l’écriture, perte de données. | Privilégier le GeoPackage ou PostGIS pour les jeux de données volumineux. | | 4 | Négliger la sauvegarde régulière de la base PostGIS. | Perte irréversible en cas de panne matérielle ou de corruption. | Mettre en place une sauvegarde automatisée quotidienne (`pg_dump -Fc`) et tester la restauration mensuellement. | | 5 | Effectuer une requête spatiale sans utiliser l’index GIST. | Temps de réponse très lent sur de gros tableaux (secondes à minutes). | Toujours créer un index `USING GIST (geom)` sur les colonnes géométriques fréquemment interrogées. | | 6 | Ne pas contrôler la validité topologique des géométries avant insertion. | Géométries auto‑intersectantes ou non fermées provoquant des erreurs d’analyse. | Utiliser `SELECT ST_IsValid(geom) FROM parcelle WHERE NOT ST_IsValid(geom);` et corriger avec `ST_MakeValid`. | --- *Fin du cours.*