Travaux Pratiques de Teledetection Spatiale II Uliege
Travaux Pratiques de Teledetection Spatiale II Uliege
Travaux Pratiques de Teledetection Spatiale II Uliege
Pratiques de
Télédétection
Spatiale II
Antoine DENIS
1
Avant-propos
Ce manuel présente une série d’exercices de télédétection spatiale, prêts à l’emploi, avec
données et logiciels libres. Les exercices vont du « niveau découverte » à un « niveau plus
avancé » tout en restant très abordables pour des étudiants de Master 1 et 2. Les exercices
peuvent se réaliser en autonomie grâce aux instructions méthodologiques détaillées et aux
nombreuses illustrations. En plus des exercices en tant que tels de nombreuses ressources
complémentaires sont proposées (sites web, applications).
Ce manuel et les données y afférentes sont libres de droit et redistribuables moyennant une
mention de la source*. Ils sont téléchargeables gratuitement à l’adresse web :
https://orbi.uliege.be/handle/2268/240835.
Merci de signaler toute erreur rencontrée dans ce manuel par email à l’adresse ci-dessus.
Les mises à jour de ce manuel seront disponibles via l’adresse web ci-dessus.
Images de couverture : images extraites du livre « EARTH as ART », NASA disponible à l’adresse
http://www.nasa.gov/connect/ebooks/earth_art_detail.html et vue d’artiste du satellite « Sentinel 2 ».
* Comment citer ce manuel ? : « Travaux Pratiques de télédétection spatiale II - Avec données et logiciels libres, Antoine
DENIS, 2021, ULIEGE, Belgique, https://orbi.uliege.be/handle/2268/240835 ».
** N’imprimez pas ces notes ! Si vous devez les imprimer, utilisez le verso de feuilles déjà utilisées, imprimez recto-verso
et/ou imprimez en 2 pages par feuilles.
2
Du même auteur
Sont également disponibles, du même auteur, les manuels de Travaux Pratiques suivants:
Adresse de
Manuel de TP Année Pages Commentaires
téléchargement
3
Table des matières
Avant-propos ...................................................................................................................... 2
Du même auteur ................................................................................................................. 3
Table des matières ............................................................................................................. 4
1. Introduction ............................................................................ 17
1.1. Objectifs de ce manuel ............................................................................................. 17
1.2. Accès aux données et ressources relatives à ce manuel ......................................... 18
1.3. Conseil d’utilisation de ce manuel ........................................................................... 19
1.4. Liens avec les cours ULIEGE ...................................................................................... 20
4
4.5. Ajout des bandes spectrales individuelles d’une image satellite dans un projet QGIS
et visualisation des bandes en « dégradé de gris » ................................................. 40
4.6. Accès aux propriétés des bandes spectrales d’une image satellite ......................... 41
4.7. Fusion de bandes spectrales individuelles d’une image satellite en un seul fichier
unique multi-bandes ................................................................................................ 43
4.8. Visualisation de compositions colorées d’une image satellite ................................ 45
4.8.1. Composition « vraies couleurs » et composition « fausses couleurs »..... 45
4.8.2. Amélioration de la visualisation d’une image satellite ............................. 48
4.8.2.1. Symbologie ....................................................................................... 48
4.8.2.2. Transparence.................................................................................... 51
4.8.2.3. Rendu ............................................................................................... 51
4.8.2.4. Pyramides ......................................................................................... 51
4.9. Autres outils d’interrogation des images ................................................................. 52
4.9.1. L’outil « Identifier les entités » .................................................................. 52
4.9.2. Outil de mesures de distances, superficies et angles ................................ 53
6
8.2.3. Choix du logiciel : QGIS 3.10 (ou version ultérieure)............................... 126
8.2.4. Contextualisation ..................................................................................... 126
8.3. Préparation du projet QGIS .................................................................................... 129
8.3.1. Création d’un nouveau projet QGIS ou ouverture du projet QGIS
préexistant ............................................................................................... 129
8.3.2. Identification du système de coordonnées à utiliser .............................. 131
8.3.3. Présentation et installation du plugin « Semi-Automatic Classification
Plugin » (SCP) ........................................................................................... 132
8.4. Données auxiliaires : présentation, création, acquisition et prétraitements ........ 134
8.4.1. Données auxiliaires vectorielles .............................................................. 134
8.4.2. Données auxiliaires raster ....................................................................... 135
8.4.2.1. Raster d’altitude ............................................................................. 135
8.4.2.1.1. Téléchargement du Modèle Numérique de Terrain (MNT) ............................ 136
8.4.2.1.2. Ajout des 4 MNT dans QGIS ............................................................................ 138
8.4.2.1.3. Fusion des 4 MNT ........................................................................................... 138
8.4.2.1.4. Reprojection du MNT fusionné dans le bon système de coordonnées .......... 138
8.4.2.1.5. Découpage du MNT sur la zone d’étude ......................................................... 139
8.4.2.1.6. Rééchantillonner le MNT à une taille de pixel plus faible ............................... 140
7
8.5.2.3.1. Recherche de l’image Landsat-5 via le site web « Earth Explorer » ................ 150
8.5.2.3.2. Téléchargement de l’image Landsat-5 via le site web « Earth Explorer »....... 154
9
8.7.5. Observation multi-temporelle ................................................................. 216
8.7.6. Observation des images en 3D via une représentation du relief naturel 217
8.7.7. Utilisation de données complémentaires pour la photo-interprétation
d’une image ............................................................................................. 220
8.8. Classification supervisée par maximum de vraisemblance ................................... 222
8.8.1. La classification supervisée par maximum de vraisemblance, une méthode
de classification parmi d’autres ............................................................... 222
8.8.2. La classification supervisée : un processus itératif.................................. 224
8.8.3. Sélection des « ROI – Region Of Interest » de calibration ou aires
d’entrainement pour chaque classe : bonnes pratiques ........................ 224
8.8.4. Création des ROI de calibration dans SCP ............................................... 227
8.8.4.1. Préalables à la création de ROI ...................................................... 227
8.8.4.1.1. Définir le « Jeu de bandes » à considérer pour la création des ROI ................ 227
8.8.4.1.2. Paramétrage des préférences SCR de QGIS .................................................... 228
10
8.11.1. Considérations générales sur l’évaluation de la qualité d’une classification
249
8.11.1.1. Evaluation de la qualité d’une classification et notions de calibration
vs validation ................................................................................... 249
8.11.1.2. Stratégies d’échantillonnage des zones de validation (ROI de
validation)....................................................................................... 252
8.11.1.3. De la qualité d’une classification.................................................... 253
8.11.2. Validation par échantillonnage stratifié aléatoire ................................... 254
8.11.2.1. Ressources complémentaires ........................................................ 254
8.11.2.2. Les grandes étapes ......................................................................... 254
8.11.2.3. Calcul du nombre de points de validation par classe .................... 254
8.11.2.4. Création des points de validation .................................................. 256
8.11.2.5. Attribution des classes aux ROI de validation via photo-
interprétation (ou autre technique) .............................................. 258
8.11.2.6. Création de la matrice de confusion (avant et après regroupement
des classes et/ou post-traitements) .............................................. 259
8.11.2.7. Matrice de confusion : construction et interprétation, un exemple
illustré 262
8.12. Matrice de transition .............................................................................................. 266
8.12.1. Définition de la matrice de transition ..................................................... 266
8.12.2. Création de la matrice de transition........................................................ 266
8.12.3. Interprétation de la matrice de transition .............................................. 267
8.12.4. Mise en relation des précisions de classification et de l’ampleur des
changements détectés ............................................................................ 267
8.13. Graphiques exprimant l’évolution de l’occupation du sol ..................................... 268
8.13.1. Histogramme par classe .......................................................................... 269
8.13.2. Camembert par classe ............................................................................. 269
8.14. Réalisation de cartes d’occupation du sol - Mise en page cartographique ........... 270
8.15. Discussion et limitations de l’approche ................................................................. 271
11
9.2. Objective (video timing: 0 min 00 sec) ................................................................... 273
9.3. Technical description of the 2 Sentinel-1 images used.......................................... 273
9.4. Download of the 2 Sentinel-1 images (not in the video) ....................................... 275
9.5. SNAP software - Download and installation (not in the video) ............................. 276
9.6. Opening of the 2 Sentinel-1 images in SNAP (video timing: 0 min 30 sec)............ 277
9.7. Main steps of the workflow described in this section (video timing: 2 min 18 sec)
............................................................................................................................... 278
9.8. Visual analysis of the 2 images side by side (video timing: 3 min 09 sec) ............. 279
9.9. Analysis of images Metadata (video timing: 3 min 15 sec) .................................... 280
9.10. Spatial subset of the images, on an area common to the 2 images, where the main
floods occur (video timing: 3 min 53 sec) .............................................................. 281
9.11. Intermediate steps: closing/opening images (video timing: 6 min 00 sec) ........... 282
9.12. Multilooking (video timing: 6 min 30 sec) .............................................................. 283
9.13. Radiometric Calibration (video timing: 7 min 35 sec) ............................................ 284
9.14. Intermediate steps: visualization (video timing: 8 min 22 sec) ............................. 285
9.15. Terrain Correction (video timing: 10 min 22 sec) .................................................. 287
9.16. Visualization of the terrain corrected images (video timing: 11 min 32 sec) ........ 288
9.17. Combination of the Archive and Crisis images to produce a RGB composite of the
flooded area, preliminary step: images stack (video timing: 13 min 25 sec) ........ 290
9.18. Overlay of the stacked images for visual comparison of the 2 images (optional) . 292
9.19. RGB composite image (video timing: 16 min 50 sec)............................................. 293
9.20. Interpretation of the RGB composite image (video timing: 17 min 28 sec) .......... 294
9.21. Exportation of the RGB composite of the flood map in another format and
visualization of the result in GOOGLE EARTH PRO (video timing: 19 min 22 sec) 295
9.22. Autre tutoriel: Satellite Remote Sensing of Flood Monitoring and Management
(NASA ARSET) ......................................................................................................... 297
12
10.4.1. Commande des images NDVI via le « Copernicus Global Land Service » 301
10.4.2. Téléchargement via « FileZilla Client » des images NDVI commandée... 308
10.4.3. Dé-compression des images téléchargées .............................................. 309
10.4.4. Téléchargement des fichiers vectoriels shapefile ................................... 312
10.5. Aperçu des données dans QGIS (facultatif) ........................................................... 313
10.6. Utilisation de SPIRITS.............................................................................................. 313
10.6.1. Téléchargement et installation de SPIRITS .............................................. 313
10.6.2. Démarrage de SPIRITS et découverte de l’interface ............................... 313
10.6.3. Création d’un projet SPIRITS.................................................................... 314
10.6.4. Importation des images dans SPIRITS et conversion de leur format ...... 317
10.6.5. Visualisation d’une série d’images NDVI ................................................. 322
10.6.5.1. Création d’un « map template » ou template de carte (QNQ)...... 323
10.6.5.2. Production d’une série de cartes à partir d’un « map template » 323
10.6.6. Extraction de statistiques à partir des images NDVI ............................... 326
10.6.6.1. Création de la structure de la base de données ............................ 326
10.6.6.2. Création de l’image raster des zones d’intérêt pour l’extraction des
statistiques ..................................................................................... 333
10.6.6.3. Extraction des statistiques sur les valeurs de NDVI à partir des
images 335
10.6.6.4. Visualisation des profils temporels de NDVI par région ................ 339
10.6.7. Lissage temporel des images NDVI et extraction des statistiques .......... 346
10.6.7.1. Concept du lissage temporel d’images satellites ........................... 346
10.6.7.2. Lissage temporel des images NDVI ................................................ 346
10.6.7.3. Extraction des statistiques sur les images NDVI lissées
temporellement ............................................................................. 348
10.6.8. Création de carte d’anomalie .................................................................. 351
10.6.8.1. Calcul de « Long Term Average » (LTA) des images NDVI ............. 351
10.6.8.2. Calcul d’anomalie du NDVI via la méthode « Relative difference »
353
10.6.8.3. Visualisation des images d’anomalie de NDVI ............................... 355
10.6.9. Analyse des variables phénologiques ...................................................... 359
10.6.9.1. Présentation du concept d’extraction de variables phénologiques à
partir de courbes d’évolution du NDVI .......................................... 359
13
10.6.9.2. Calcul des variables phénologiques à partir des images NDVI
temporellement lissées (résultat = des images correspondant à des
paramètres phénologiques) ........................................................... 360
10.6.9.3. Visualisation des images correspondant aux paramètres
phénologiques ................................................................................ 363
10.6.9.4. Création de la structure de la base de données pour accueillir les
variables phénologiques ................................................................ 370
10.6.9.5. Extraction des statistiques par zone (RUM) pour les variables
phénologiques ................................................................................ 371
10.6.9.6. Visualisation des statistiques par zone (RUM) pour les variables
phénologiques ................................................................................ 373
14
13. Introduction à GOOGLE EARTH ENGINE (GEE) ........................ 391
13.1. Préalable !! ............................................................................................................. 392
13.2. Présentation de Google Earth Engine (GEE) (30 min) ............................................ 392
13.2.1. Page web principale................................................................................. 392
13.2.2. Vidéo d’introduction (2 minutes) ............................................................ 393
13.2.3. Timelapse................................................................................................. 393
13.2.4. Exemples de cas d’étude avec GEE.......................................................... 394
13.2.5. Les datasets ............................................................................................. 394
13.2.6. 2 Interfaces: EXPLORER et CODE EDITOR ................................................ 395
13.2.6.1. EXPLORER ....................................................................................... 395
13.2.6.2. CODE EDITOR: ................................................................................ 395
13.2.7. Documentations ...................................................................................... 397
13.3. Exercices ................................................................................................................. 398
13.3.1. GEE EXPLORER ......................................................................................... 398
13.3.1.1. Découverte de l’EXPLORER et classification supervisée (20 minutes)
398
13.3.1.2. Advanced - Visualizing change over time (5 min) .......................... 399
13.3.1.3. Raster de classification (2 minutes) ............................................... 400
13.3.1.4. Advanced - Viewing Classified Rasters (5 min) .............................. 400
13.3.2. GEE CODE EDITOR.................................................................................... 400
13.3.2.1. Découvrir et visualiser une collection via le code editor (45 minutes)
400
13.3.2.2. Tutoriel officiel « The Earth Engine API »....................................... 404
13.3.2.2.1. Visualizing Images and Image Bands (20 minutes) ......................................... 404
13.3.2.2.2. Computations using Images (20 min) ............................................................. 405
13.3.2.2.3. Image Collections (20 minutes) ...................................................................... 406
13.3.2.2.4. Compositing, Masking, and Mosaicking (20 minutes) .................................... 407
13.3.2.2.5. NDVI, Mapping a Function over a Collection, Quality Mosaicking .................. 409
13.3.2.2.6. Exporting Charts and Images .......................................................................... 409
15
14.2. Traduction des commentaires de la vidéo HOLOGLOBE_v3 .................................. 411
16
1. Introduction
Ce manuel propose :
17
1.2. Accès aux données et ressources relatives à ce manuel
Ce manuel et les données nécessaires pour faire les exercices sont disponibles gratuitement
via le site web ORBI de l’ULIEGE à l’adresse :
https://orbi.uliege.be/handle/2268/240835
18
1.3. Conseil d’utilisation de ce manuel
Aidez-vous, pour naviguer dans ce document, de la table des matières et des « signets » du
document PDF disponibles via le panneau latéral de gauche de ce document (Figure 1).
Vous remarquerez que de nombreux hyperliens renvoyant à des sections parfois éloignées
dans le document sont utilisés. Pour revenir facilement à l’endroit à partir duquel l’hyperlien
a été cliqué, utilisez, dans le logiciel Adobe Acrobat Reader, les touches « alt flèche
gauche ».
Figure 1 : Utilisez les signets pour naviguer plus facilement dans ce manuel !
19
1.4. Liens avec les cours ULIEGE
Ce document est le manuel des travaux pratiques de télédétection spatiale utilisé pour les
cours ULIEGE (Belgique) suivants :
Les parties de ce manuel vues lors de ces différents cours sont sélectionnées en fonction des
objectifs spécifiques de chaque cours.
Ce manuel présente, en plus des exercices réalisés dans le cadre de ces cours avec les
étudiants, une série de ressources supplémentaires.
20
2. Notions de télédétection - Quelques
rappels théoriques
2.1. La télédétection
La télédétection (détection à distance, Remote Sensing (RS) en anglais) désigne, dans son
acception la plus large, la mesure ou l'acquisition d'informations sur un objet ou un
phénomène, par l'intermédiaire d'un instrument de mesure n'ayant pas de contact avec
l'objet étudié (Figure 2).
C'est l'utilisation à distance (par exemple, d’un drone, d'un avion, d'un engin spatial, d'un
satellite ou encore d'un bateau) de n'importe quel type d'instrument permettant
l'acquisition d'informations sur l'environnement. On fait souvent appel à des instruments
tels qu’appareils photographiques, lasers, radars, sonars, sismographes ou gravimètres.
21
L’on distingue en général les moyens de télédétection « actif » et « passif » (Figure 3):
22
Certains instruments de télédétection utilisent des ondes sonores de façon similaire, et
d'autres mesurent des variations dans des champs magnétiques ou gravitaires.
23
2.3. La signature spectrale et les capteurs satellites
La signature spectrale d’un élément/d’une matière correspond à la représentation
graphique de la variation de sa réflectance (axe Y) en fonction des longueurs d’onde du
spectre auquel le capteur utilisé est sensible (axe X) (Figure 5). Différents éléments/matières
(ex : terre vs végétation) ou des éléments/matières présentant différents états (ex : plus ou
moins humide/sec) présenteront en général une signature spectrale différente. La signature
spectrale peut donc être utilisée pour identifier/reconnaitre/caractériser différents
éléments/matières et/ou leur état.
Les capteurs multispectraux embarqués sur la plupart des satellites d’observation de la terre
n’enregistrent pas la signature spectrale de la surface terrestre avec autant de précision que
présenté à la Figure 5. Ils enregistrent en général l’énergie réfléchie par la surface de la terre
dans quelques gammes assez larges de longueurs d’ondes tel que présenté pour quelques
capteurs satellites bien connus à la Figure 6. La « signature spectrale » dérivée de ce type de
capteur ne sera donc composée que de quelques valeurs, une par bande spectrale.
Cependant certains capteurs satellites (rares), dit « hyperspectraux », peuvent enregistrer
une signature spectrale très précise grâce à un très grand nombre de bandes spectrales
étroites, comme pour le capteur CHRIS embarqué sur le satellite PROBA 1 (Figure 6).
24
Figure 6 : Signature spectrale typique de la végétation verte et résolutions spectrales et spatiales
de différents capteurs satellites multispectraux bien connus et d’un capteur hyperspectral (capteur
CHRIS sur le satellite PROBA 1). Les rectangles de couleurs correspondent aux différentes bandes
spectrales disponibles dans une image satellite. (Source : adapté de DENIS A, 2018, Can satellites
help organic crops certification? https://orbi.uliege.be/handle/2268/226209)
25
2.4. Exemples d’applications de la télédétection
Les exemples d’applications de la télédétection spatiale sont très nombreux :
26
2.5. Les différents types de résolutions des images
satellites
Série temporelle de 36 images NDVI décadaires (36 images * 10 jours = 360 jours = 1 an)
issues du capteur SPOT VEGETATION pour le Pakistan (gauche) et évolution temporelle du
NDVI en un pixel présentant une alternance marquée de l’état de la végétation selon les
saisons (droite). La résolution temporelle de cette série est de « 1 image / 10 jours ».
27
2.5.3. Résolution spectrale
o Aptitude d'un système de détection à distinguer des rayonnements
électromagnétiques de fréquences (/de longueurs d’onde) différentes.
o Le nombre et la largeur (gamme de longueurs d’ondes) des bandes spectrales d’une
image (ou auxquelles un capteur est sensible).
o Plus la résolution spectrale est fine, plus les gammes de longueurs d’ondes des
différents canaux du capteur sont étroites.
28
2.5.4. Résolution radiométrique
o La résolution radiométrique d'un système de télédétection décrit sa capacité à
reconnaître de petites différences dans l'énergie électromagnétique. Plus la
résolution radiométrique d'un capteur est fine, plus le capteur est sensible à de
petites différences dans l'intensité de l'énergie reçue.
o La résolution radiométrique est définie par le nombre maximum de niveaux
d'intensité d'énergie électromagnétique que le capteur peut enregistrer et dépend
du nombre de bits utilisés pour représenter l'intensité enregistrée, par exemple :
o Si 8 bits sont utilisés, 28 soit 256 valeurs numériques allant de 0 à 255 sont
disponibles
o Si 4 bits sont utilisés, 24 soit 16 valeurs numériques allant de 0 à 15 sont
disponibles
4 bits 8 bits
2 bits 8 bits
29
2.6. Longueurs d’ondes et fréquences du spectre
électromagnétique
(Source: http://en.wikipedia.org/wiki/Electromagnetic_spectrum)
30
3. Découverte de « Google Earth Pro »
Objectifs
Découverte de Google Earth Pro et de ses principales fonctionnalités.
pédagogiques
Temps
~ 30 minutes.
approximatif
Prérequis Aucun.
31
« Google Earth Pro » est probablement la plus belle application de la télédétection spatiale
réalisée à ce jour. C’est une représentation de la planète terre en 3D, interactive et
disponible en ligne, construites à partir d’un très grand nombre (milliers/millions ?) d’images
satellites issues de différents capteurs, de différentes résolutions spatiales et enregistrées à
différents moments. Il ne s’agit donc pas (encore) d’une visualisation en directe de la terre !
Google Earth Pro intègre une série de fonctionnalités et de très nombreuses données
supplémentaires (limites administratives, photos, lieux, météo, relief, bâtiments 3D, etc)
(Figure 8 et Figure 9).
https://www.google.com/intl/fr_ALL/earth/
Google Earth
32
Figure 8 : Aperçu de quelques fonctionnalités de Google Earth Pro
33
Figure 9 : Aperçu de la fonctionnalité « Images d’archive » de Google Earth Pro pour la Mer d’Aral
34
Dézoomez et positionnez-vous au-dessus de quelques montagnes. Via l’activation du
bouton « Lumière du soleil en fonction de l’heure du jour » dans la partie
supérieure de l’interface et l’utilisation de la ligne du temps qui apparait, constatez
l’effet de la position du soleil sur le réflectance des images, effet particulièrement
visible en zone montagneuse (cet effet d’ombrage est modélisé).
Tapez ensuite « Manhattan, New-York » dans l’onglet « Recherche ». Vous devriez
arriver à « Central Park ». Déplacez-vous ensuite vers la pointe sud de Manhattan et
cochez l’option « Bâtiments 3D » dans le panneau « Calques ». Inclinez à nouveau la
fenêtre de visualisation dans le plan horizontal et visitez Manhattan en 3D. Déplacez-
vous ensuite vers la « Statue of Liberty ». (Remarque : tout aussi impressionnant, le
campus d’Arlon au Sud de la Belgique est également modélisé en bâtiment 3D).
Déplacez-vous ensuite vers Paris et explorez « l’Avenue des Champs Elysées » et
Google Earth permet également d’explorer l’espace, la Lune et Mars via le bouton .
Rem : Attention ! Il est possible de récupérer les images de Google Earth :
Via une « capture d’écran » (« print screen » en anglais) ou
Via le menu de Google Earth « Fichier > Enregistrer > Enregistrer l’image… ».
Temps
~ 45 minutes
approximatif
Il est préférable, mais non obligatoire, que l’utilisateur ait déjà quelques
Prérequis
notions de base concernant la gestion d’un projet QGIS.
36
Vraies couleurs RGB = B3-B2-B1 Vraies couleurs RGB = B3-B2-B1 Fausses couleurs RGB = B4-B7-B3
SANS amélioration de contraste AVEC amélioration de contraste AVEC amélioration de contraste
Figure 10 : Image Landsat-5 TM utilisée dans le cadre de cet exercice avec différents effets de
visualisation : amélioration de contraste et compositions colorées.
Plus d’information sur cette image et la manière de la rechercher et télécharger sur le net
sont données à la section « 8.5.2.3 Recherche et téléchargement de l’image Landsat-5 via le
site web « EarthExplorer » page 150. Ces informations ne sont cependant pas utiles pour cet
exercice-ci.
Les caractéristiques techniques les plus importantes de cette image, dont la résolution
spatiale et la résolution spectrale, sont reprises dans les Figure 12 et Figure 11 ci-dessous.
37
Nom du fichier téléchargé LT05_L1TP_173060_19870807_20170211_01_T1
Nom du satellite LANDSAT 5 TM
Lieu dans le système LANDSAT (Path, Row) WRS_PATH = 173
et « Tile » de Sentinel-2 WRS_ROW = 60
Système de coordonnées EPSG 32635
Etendue spatiale approximative d’une 170 km north-south by 183 km east-west
image/tile
Date d’acquisition 07 08 1987
Résolution spatiale 30 m
Nombre de bandes spectrales 7
Couverture nuageuse Partielle
Format original des images lors du .tif
téléchargement
Poids numérique du dossier 138 Mo / 460 Mo
compressé/décompressé contenant une
image
Site web de référence pour le https://earthexplorer.usgs.gov/
téléchargement
Page web reprenant les caractéristiques https://www.usgs.gov/land-resources/nli/landsat/landsat-5?qt-
science_support_page_related_con=0#qt-
techniques des images science_support_page_related_con
Figure 11 : Caractéristiques principales de l’image Landsat-5 Thematic Mapper ™ utilisée dans le
cadre de cet exercice.
Numéro de la bande Longueur d’onde (µm) Nom de la bande Résolution spatiale (m)
Band 1 0.45-0.52 Blue 30
Band 2 0.52-0.60 Green 30
Band 3 0.63-0.69 Red 30
Band 4 0.76-0.90 Near Infrared (NIR) 30
Band 5 1.55-1.75 Shortwave Infrared (SWIR) 1 30
Band 6 10.40-12.50 Thermal 120 (30)*
Band 7 2.08-2.35 Shortwave Infrared (SWIR) 2 30
Figure 12 : Description des 7 bandes spectrales des images Landsat-5 Thematic Mapper ™. (Source :
https ://www.usgs.gov/faqs/what-are-band-designations-landsat-satellites ?qt-news_science_products=0#qt-
news_science_products)
* TM Band 6 was acquired at 120-meter spatial resolution, but products are resampled to
30-meter pixels.
38
4.2. Visualisation des fichiers correspondant à l’image
satellite dans l’explorateur de l’ordinateur
Via l’explorateur de votre ordinateur (explorateur Windows ou autre) l’image satellite
Landsat-5 se présente dans le dossier « …\1_VISUALISATION\... » comme illustré à la Figure
13 :
Figure 13 : Visualisation des fichiers correspondant à l’image satellite Landsat-5 dans l’explorateur
de l’ordinateur (explorateur Windows dans ce cas-ci).
39
4.3. Ouverture du projet QGIS à utiliser
Ouvrez le projet QGIS préalablement créé pour vous et disponible dans le répertoire
« …\1_VISUALISATION\VISUALISATION.qgz » en double-cliquant dessus
Ce projet QGIS est presque vide mais contient une couche shapefile des pays du monde qui
vous permettra de facilement vous localiser.
Vous pouvez éventuellement ajouter une couche de fond (« base map ») qui correspond par
exemple aux images satellites « Google Satellite » via l’utilisation de l’extension QGIS
« QuickMapServices » (à installer et paramétrer correctement si ce n’est pas encore fait
comme indiqué dans le manuel de TP QGIS disponible ici
https://orbi.uliege.be/handle/2268/190559). Cette couche de fond en très haute résolution
spatiale peut parfois être utile pour mieux comprendre le paysage de la zone étudiée.
Cliquez sur le menu « Projet > Propriétés… > SCR > et tapez « 32635 » dans le filtre de
recherche de SCR afin de trouver plus facilement ce SCR > Sélectionnez-le > Cliquez
sur OK
Les images satellites sont des données de type raster et doivent donc s’ouvrir
Via le bouton « Gestionnaire des sources Open Data » > Raster > Naviguer vers le
répertoire
« 1_VISUALISATION\DATA\LT05_L1TP_173060_19870807_20170211_01_T1\ » >
Sélectionner les bandes B1 à B7.tif > Open > Ajouter > Fermer
40
Pour plus de facilité de manipulation, placez ces bandes dans le groupe de couche
«LANDSAT BRUT » du projet QGIS préexistant, par simple glissé-déposé dans le
panneau « Couches » de QGIS
Chacune des 7 bandes est visible en « dégradé de gris » (en anglais : « gray scale »).
Les 7 bandes étant reprises comme 7 couches indépendantes, il n’est pas possible de réaliser
directement une « composition colorée » multi-bandes (une image couleur RGB). Pour
pallier à ce problème, il est possible de « fusionner » les 7 bandes en un seul fichier unique
multi-bandes (confer section 4.7 page 43 ci-dessous).
Les propriétés de chaque couche d’information ajoutée dans QGIS sont disponibles via :
Ces informations renseignent notamment sur le SCR dans lequel est exprimée la couche
d’information, le nombre de pixels constituant l’image en ligne (hauteur) et colonne
(largeur), le format de l’image (GeoTIFF dans ce cas-ci), quelques statistiques sur les bandes
spectrales, la taille du pixel, etc.
41
Figure 14 : Accès aux propriétés d’une couche raster (bande 1 d’une image Landsat-5 TM ici) dans
QGIS
Ces informations et bien d’autres sont également reprises dans les fichiers de type texte
(avec l’extension « .txt ») qui accompagnent l’image au moment du téléchargement depuis
le site web source. Ces fichiers textes, que vous pouvez ouvrir par curiosité avec le logiciel
« bloc note », sont les fichiers (dans le répertoire
« …\1_VISUALISATION\DATA\LT05_L1TP_173060_19870807_20170211_01_T1\ ») :
42
4.7. Fusion de bandes spectrales individuelles d’une image
satellite en un seul fichier unique multi-bandes
A l’aide des indications ci-dessous, fusionnez les 7 bandes individuelles de l’image Landsat-5
TM en un seul fichier unique multi-bandes.
Pour fusionner les différentes bandes spectrales en une seule image multi-bandes (Figure
9) :
Utilisez la fonction QGIS disponible dans le menu « Raster > Divers > Fusionner… » et
paramétrez-la comme indiqué ci-dessous.
Couches en entrée : sélectionner les 7 bandes B1 à B7 de l’image Landsat-5 à
fusionner
Cocher l’option « Placer chaque fichier en entrée dans une bande séparée », ce qui
permettra d’obtenir un raster multi-bandes. Sinon le raster résultant ne contient
qu’une seule bande.
Fusionné : indiquer le répertoire de sortie et nommer le raster en sortie, par exemple
« LT05_L1TP_173060_19870807_20170211_01_T1_FUSION_B1_B7.TIF »
Décocher « Ouvrir le fichier en sortie après l’exécution de l’algorithme » car cela
ajouterait automatiquement le raster fusionné dans le projet QGIS mais avec le nom
par défaut « Fusionné », ce qui n’est pas idéal.
Cliquer sur « Exécuter »
Après exécution de la fonction (quelques secondes sont nécessaires) et vérification
de l’absence de message d’erreur dans le journal de la fonction « Fusionner », fermer
la fenêtre de cette fonction.
Ajouter le raster fusionné produit manuellement dans votre projet QGIS via le bouton
d’ajout de raster. Pour plus de facilité de manipulation, placez-le dans le groupe de
couche du panneau « Couches » «LANDSAT FUSION »
43
Figure 15 : Fusion des 7 bandes spectrales d’une image Landsat-5 TM pour produire un raster
multi-bandes, dans QGIS via la fonction « Fusion »
44
4.8. Visualisation de compositions colorées d’une image
satellite
A l’aide des indications ci-dessous, réalisez une composition vraies couleurs et 2-3
compositions fausses couleurs de qualité de l’image satellite Landsat-5 TM.
Prenez ensuite quelques secondes pour photo-interpréter l’image satellite, c’est-à-dire
reconnaitre les différentes classes d’occupation du sol (+ nuages) présentes sur cette
image.
45
Vraies couleurs RGB = B3-B2-B1 Fausses couleurs RGB = B4-B7-B3
Figure 16 : Image Landsat-5 TM en composition « vraies couleurs » (à gauche) et en composition
« fausses couleurs » (un exemple parmi les nombreuses possibles) (à droite), toutes deux avec
amélioration de contraste.
Pour afficher une « composition vraies couleurs » d’une image satellite dans QGIS (Figure
17):
Cliquez-droit sur l’image satellite multi-bandes fusionnée (la fusion de l’image a été
faite à la section 4.7 page 43) > Propriétés… > Symbologie
Dans le menu « Rendu des bandes raster » choisissez
o « Type de rendu » = « Couleur à bandes multiples »
o « Bande rouge » = Bande 3 de l’image Landsat-5 TM qui correspond donc à la
longueur d’onde du rouge
o « Bande verte » = Bande 2 de l’image Landsat-5 TM qui correspond donc à la
longueur d’onde du vert
o « Bande bleue » = Bande 1 de l’image Landsat-5 TM qui correspond donc à la
longueur d’onde du bleu
ATTENTION !
L’ordre des bandes RGB=3-2-1 correspond à une composition colorée vraies couleurs
dans ce cas-ci uniquement parce que les bandes spectrales numéro 3-2-1 des images
Landsat-5 TM correspondent aux longueurs d’onde rouge-vert-bleu respectivement.
Si un autre capteur est utilisé, il faut se renseigner sur la correspondance entre
numéro de bandes et longueurs d’onde et adapter les 3 bandes à affecter aux 3
canaux RGB pour une composition vraies couleurs. Par exemple, pour les images de
type Landsat-8, une composition vraies couleurs sera obtenue via la composition RGB
= 4-3-2 (confer ici https://www.usgs.gov/faqs/what-are-band-designations-landsat-
satellites?qt-news_science_products=0#qt-news_science_products)
46
Si un capteur ne possède pas les 3 bandes spectrales du rouge, vert et bleu, alors il
devient impossible de réaliser une composition vraies couleurs pour ce capteur.
Le rendu visuel par défaut d’une image satellite dans QGIS est souvent de mauvaise qualité.
A vous donc d’améliorer la visualisation de l’image via les différentes possibilités reprises ci-
dessous.
47
4.8.2. Amélioration de la visualisation d’une image satellite
Le rendu visuel par défaut d’une image satellite dans QGIS est souvent de mauvaise qualité.
A vous donc d’améliorer la visualisation de l’image via les différentes possibilités reprises ci-
dessous.
Les possibilités d’amélioration de la visualisation d’une image satellite sont disponibles via
les menus:
4.8.2.1. Symbologie
Ce menu est le plus important et permet notamment de (Figure 17):
48
données sources (la valeur des pixels). Elle ne fait qu’attribuer différemment
la palette de couleurs disponibles à votre image afin de faire ressortir certains
éléments plus clairement.
o Dans la Figure 18 - droite, l’amélioration de contraste appliquée est celle de
la Figure 17. La présence de nuages dans la partie nord de l’image résulte en
des pixels qui saturent et présentent la valeur maximale possible pour ce type
d’image (soit 255 dans ce cas-ci vu que l’image est codée sur 8 bit (28= 256
valeurs possibles allant de 0 à 255)). La distribution de la fréquence des
valeurs de pixels pour chacune des bandes spectrales de l’image Landsat-5
TM (Figure 19) montre très clairement que la gamme de valeurs de la plupart
des pixels pour une bande spectrale donnée est assez étroite (entre 50 et 100
pour la bande 1 par exemple) mais que certains pixels (ceux identifiant les
nuages) présentent une valeur saturant à 255. Ceci a pour conséquence que
la gamme de couleurs possibles est mal distribuée entre d’une part des
pixels présentant de relativement faibles valeurs (pour les diverses
occupations du sol) et d’autres part des pixels présentant la valeur maximale
possible (255 pour les nuages), et résulte en une image très contrastée avec
des pixels très sombres et des pixels très clairs, ce qui n’est pas idéale pour
une bonne analyse visuelle de l’image.
o L’amélioration de contraste peut être calculée soit, via le menu,
« Statistiques de l’emprise > Raster entier », sur l’ensemble du raster.
Cette option présente le grand avantage de résulter en une
symbologie fixe et donc comparable indépendamment de la partie du
raster visualisée.
« Statistiques de l’emprise > Emprise actualisée », sur la partie du
raster qui est affichée à l’écran. Cette option présente le grand
avantage de résulter en une symbologie qui s’adapte dynamiquement
à la partie du raster affichée à l’écran, ce qui est très pratique pour
visualiser via « zoom in » des parties de raster qui sont soit trop
claires, soit trop sombres lorsque le raster est affiché dans son
entièreté.
Modifier d’autres paramètres de visualisation (luminosité, contraste, etc)
o Via les menus « Rendu de la couleur » et « Ré-échantillonage »
49
Figure 19 : Histogramme de la distribution de la fréquence des valeurs de pixels pour chacune des
bandes spectrales de l’image Landsat-5 TM
50
4.8.2.2. Transparence
Diverses options de transparence de l’image (Figure 20) sont disponibles via le menu :
4.8.2.3. Rendu
Permet de définir la gamme d’échelle de visualisation pour laquelle l’image sera affichée à
l’écran
Clic-droit sur la couche raster > Propriétés… > Rendu > Visibilité dépendante de
l’échelle
4.8.2.4. Pyramides
« Des couches raster de haute résolution peuvent ralentir la navigation dans QGIS car leur
taille empêche leur mise en mémoire complète. La création de copies de résolution
inférieure (que l'on nomme pyramides) peut améliorer la performance d'affichage. »
Attention ! Altération possible de l’image d’origine ! (Source = menu QGIS de construction
des pyramides).
51
4.9. Autres outils d’interrogation des images
Plusieurs outils, présentés ci-dessous, sont disponibles pour accéder aux informations
caractérisant une image satellite.
L’outil « Identifier les entités » permet d’afficher les valeurs de pixels des différentes
bandes spectrales d’une image et d’obtenir les coordonnées géographiques du point
consulté.
Pour l’utiliser :
Figure 21 : Fenêtre de résultat de l’identification d’un pixel de l’image Landsat-5 TM fusionnée via
l’outil « Identifier les entités »
52
4.9.2. Outil de mesures de distances, superficies et angles
L’outil « Mesurer » permet de mesurer des distances, superficies et angles sur une image :
53
5. Analyse de changement simple dans QGIS
54
QGIS 3.12 ou supérieur
(https://www.qgis.org/fr/site/forusers/download.html).
Logiciel ou Attention ! QGIS 3.10 a présenté un bug lors de l’utilisation de la «
application Calculatrice raster » (section 5.3 page 63). QGIS 3.12 par contre a
fonctionné sans problème. Utilisez donc QGIS 3.12 ou une version
ultérieure de QGIS pour cet exercice.
Temps
~ 2 heures
approximatif
Prérequis Aucun.
55
5.1. Accès aux images satellites basse résolution de type
NDVI du capteur SPOT-VEGETATION: recherche et
téléchargement
A l’aide des indications ci-dessous, téléchargez 2 images satellites basses résolution de type
NDVI sur l’Afrique pour janvier et septembre 2006. ~ 30 minutes.
https://land.copernicus.eu/global/
56
Choisissez le menu « NDVI » : https://land.copernicus.eu/global/products/ndvi
Jetez un œil à :
o La définition du NDVI (aussi développée dans la Figure 25)
o La section « NDVI characteristics » pour le NDVI 1km, et en particulier les sous-
sections
Access
Algorithm
Technical (Figure 24)
Figure 24 : Caractéristiques des images NDVI distribuées sur le site web « Copernicus Global Land
Service » accessibles via le menu « NDVI characteristics > NDVI 1km > Technical » avec en
particulier les paramètres « scaling » et « offset » utilisés pour transformer les Digital Number (DN)
en valeurs de NDVI.
57
NDVI – Normalized Difference Vegetation Index
Le NDVI est un indice de végétation qui se définit comme la différence normalisée des mesures de
réflectance spectrale acquises dans les zones de longueurs d’onde « Proche Infra-Rouge » (« PIR ») et
« Rouge ».
( PIR Rouge)
NDVI
( PIR Rouge)
Sa valeur théorique varie entre -1 et 1. En pratique, une surface d’eau libre (océan, lac,…) prendra des
valeurs de NDVI proches de 0, un sol nu prendra des valeurs de 0.1 à 0.2, alors qu’une végétation
dense aura des valeurs de 0.5 à 0.8.
Les plantes vertes absorbent une grande partie des radiations solaires dans une zone de longueurs
d’onde appelée « PAR » (« Photosynthetically Active Radiation »), dont fait partie la zone du
« Rouge ». Par contre, le « PIR » est fortement diffusé (non absorbé : transmis et réfléchi) par la
plante. Par conséquent, une plante verte vive apparait sombre dans le PAR (dont le rouge) car ces
longueurs d’ondes sont absorbées par la plante (petite valeur de réflectance), et « brillante » dans le
« PIR » car ces longueurs d’ondes sont réfléchies en partie (grande valeur de réflectance).
58
Les valeurs de NDVI, qui varient théoriquement entre -1 et 1 (confer Figure 25), sont
encodées dans les pixels des images à l’aide de nombres entiers. Ces valeurs entières
réellement encodées dans les pixels sont appelées les « Digital Number (DN) ». Dans le cas
d’une image codée sur 8 bits, comme c’est le cas pour les images NDVI SPOT-VEGETATION,
les DN de chaque pixel peuvent prendre 28= 256 valeurs possibles, allant de 0 à 255. Les
valeurs de NDVI correspondant à ces DN peuvent être obtenues par simple transformation
mathématique (Figure 24 et Figure 26), dans ce cas-ci via la formule
Figure 26 : Principe de la conversion des DN (Digital Numbers) en valeurs NDVI via une
fonction mathématique du premier degré présentant les 2 coefficients pente (ou scaling)
(a) et intercepte (ou offset) (b).
Dans l’onglet « NDVI 1 km > Access », cliquez sur « Near-real time products »
59
Vous arrivez sur la page web
https://land.copernicus.vgt.vito.be/PDF/portal/Application.html#Browse;Root=5131
86;Collection=1000085;Time=NORMAL,NORMAL,-1,,,-1,, (Figure 27)
Figure 27 : Portail web « Copernicus Global Land Service » pour accéder aux images NDVI basse
résolution
60
Une fois inscrit, loguez-vous via le bouton « Login » en haut à droite de l’interface
(Figure 27)
Une fois logué, paramétrez votre recherche de produits comme indiqué dans la
Figure 29 en sélectionnant
o La collection « NDVI 1km V1 ». Comprendre de quel type de produit il s’agit
avec le bouton “XML” : par quelle méthode ces images sont-elles produites ?
(« 30-days composite with a sliding window of 10 days »)
o La période du 01/01/2006 au 31/12/2006
o La zone : dessinez un petit carré sur le centre de l’Afrique avec l’outil « Draw
Region Of Interest (ROI) »
o Le « Coverage » : « BioPar_NDVI_V1_CONTINENTS »
Lancez la recherche via le bouton « Search »
Identifiez dans la liste de résultats les images d’intérêt, soit les 2 images suivantes
dans ce cas-ci :
o 1ère décade de janvier 2006 : NDVI_200601030000_AFRI_VGT_V1.3
o 1ère décade de septembre 2006 : NDVI_200609030000_AFRI_VGT_V1.3
o ATTENTION les dates, reprises dans les noms des images, sont indiquées dans
le format « YYYYMMDD » ou « Année_Mois_Jours_dans_la_décade ».
61
Pré-visualisez ces 2 images à l’aide de l’icône dans la colonne « Thumbnail ». (Pour
l’évaluation, assurez-vous qu’il y a un contraste suffisant entre les 2 dates choisies
pour que votre analyse diachronique soit facile à interpréter).
Téléchargez ces 2 images en
o Cliquant sur la flèche bleue de la ligne correspondant au produit à télécharger
o Dans la fenêtre s’étant ouverte, cliquez sur la flèche bleue du « Data Files »
pour télécharger le fichier « .zip » (par exemple le fichier
« g2_BIOPAR_NDVI_200601030000_AFRI_VGT_V1.3.zip » de 46.4 Mo)
Les 2 dossiers se téléchargent en quelques secondes.
Après téléchargement, déplacez ces 2 dossiers dans votre répertoire de travail,
c’est-à-dire dans le répertoire « \2_ANALYSE_DE_CHANGEMENT_SIMPLE\ ».
Décompressez ces 2 dossiers. Ceci peut se faire par exemple avec le logiciel gratuit 7-
zip disponible ici https://www.7-zip.org/ (à n’installer que si un logiciel de ce type
n’est pas déjà installé sur votre ordinateur).
Le contenu d’un dossier « .zip » est présenté à la Figure 30. Dans ce dossier :
Les images « .tiff » sont les images qui seront utilisées dans cet exercice
Les images « .h5 » ne seront pas utilisées et peuvent être supprimées pour alléger
vos discs durs (441 Mo par image).
Figure 30 : Fichiers contenus dans les dossiers « .zip » des images NDVI basse résolution
Les images sont maintenant prêtes à être importées dans QGIS pour analyse (confer section
suivante).
62
5.2. Analyse visuelle des images satellites NDVI basse
résolution
A l’aide des indications ci-dessous, réaliser une analyse visuelle de l’évolution spatio-
temporelle des conditions de végétation entre janvier et septembre 2006 sur l’Afrique. ~
30 minutes.
Figure 31 : 2 images NDVI basse résolution du capteur SPOT-VEGETATION sur l’Afrique pour janvier
et septembre 2006.
Un dossier contenant l’ensemble de l’exercice déjà réalisé (projet QGIS et données utilisées)
est disponible dans le répertoire :
« …\2_ANALYSE_DE_CHANGEMENT_SIMPLE\»
Cependant il vous est conseillé de faire l’ensemble des manipulations par vous-même, pour
un bon apprentissage, dans un dossier séparé.
Les étapes à suivre pour réaliser une analyse visuelle des images satellites NDVI basse
résolution sont les suivantes :
Identifier le lien entre les couleurs de l’image, les valeurs de pixels (DN) et le NDVI.
o Les Figure 24 et Figure 26 nous renseignent que il y a un lien directement
proportionnel entre les DN (variant entre [0-255]) et le NDVI (variant dans la
pratique entre 0 et 1), c’est-à-dire qu’au plus les DN sont élevés au plus le
NDVI sera élevé.
o L’analyse de la palette de couleurs nous indique que :
La couleur brune correspond à de faibles DN, c’est-à-dire de faibles
valeurs de NDVI, c’est-à-dire à un développement de la végétation
verte faible
64
La couleur jaune correspond à des valeurs de DN intermédiaires, c’est-
à-dire des valeurs de NDVI intermédiaires, c’est-à-dire à un
développement de la végétation verte intermédiaire
La couleur verte correspond à des valeurs de DN élevées, c’est-à-dire
à des valeurs de NDVI élevées, c’est-à-dire à un développement de la
végétation verte plus élevé
La couleur blanche correspond aux DN 251 à 255, qui sont ce que l’on
appelle des « Flags », c’est-à-dire des « valeurs code » pouvant
signifier plusieurs choses, par exemple qu’il n’y a pas de données
valide, qu’il y a des nuages, que c’est une zone non continentale
(océan, mer), etc. Ces pixels sont donc à ignorer.
o ATTENTION ! Les couleurs utilisées sont donc des couleurs artificielles et pas
des « couleurs naturelles ». Il ne faut donc les interpréter QUE sur base de
leur signification précise et pas sur base d’une impression visuelle subjective
qui pourrait se révéler totalement fausse si la palette n’est pas intuitivement
interprétable.
L’analyse visuelle des 2 images peut être facilitée dans certains cas par l’utilisation
d’autres outils :
o La valeur des pixels de l’image est consultable également de manière
interactive via l’outil « Identifier des entités » en sélectionnant d’abord
l’image à interroger dans le panneau « Couches » puis en cliquant avec cet
outil sur l’image à interroger.
o L’utilisation de l’extension (à installer) « MapSwipe Tool » qui permet de
réaliser un swipe vertical ou horizontal sur les 2 images à comparer
o L’utilisation de plusieurs fenêtres de visualisation spatialement
synchronisées lors de la navigation mais à différentes échelles (par exemple :
une fenêtre sur l’ensemble de l’Afrique et une fenêtre de zoom à une échelle
plus grande). De nouvelles fenêtres peuvent être ajoutées via le menu « Vue
> Nouvelle vue cartographique » et ensuite être paramétrées selon les
souhaits de l’utilisateur.
Notez enfin la visibilité sur les images NDVI de quelques « poches de végétation » notables :
66
Figure 32 : Aperçu de l’Afrique telle que présentée dans Google Earth avec localisation des 2 villes
de Mopti au Mali et Gaborone au Bostwana. Le vert correspond à des zones avec végétation et le
jaune-orange aux zones plus désertiques
67
5.3. Analyse de changement mathématique
A l’aide des indications ci-dessous, réaliser une analyse de changement mathématique de
l’évolution spatio-temporelle des conditions de végétation entre janvier et septembre
2006 sur l’Afrique. ~ 60 minutes.
68
Figure 35 : Illustration de l’opération de soustraction algébrique de 2 images satellites (codées sur
8 bits 28 = 256 valeurs possibles dans la gamme [0-255]) avec des exemples de valeurs de pixels
(Digital Number (DN)) proportionnel au NDVI.
Pour calculer une simple différence mathématique entre 2 images via la « Calculatrice
Raster » de QGIS (Figure 36) :
Cliquez sur le menu « Raster > Raster Calculator… », ce qui ouvrira la « Calculatrice
Raster »
Paramétrez l’interface de la « Calculatrice Raster » comme indiqué à la Figure 36 en
nommant votre raster en sortie « DIFFERENCE_SEPTEMBRE_JANVIER.tif »
Cliquez sur « OK »
69
Ajoutez manuellement le raster de changement «
DIFFERENCE_SEPTEMBRE_JANVIER.tif » dans votre projet QGIS s’il ne s’est pas ajouté
automatiquement.
Visualisez-le brièvement (Figure 37)
Via la gestion des pixels ne correspondant pas à du NDVI Confer section 5.3.2
page 70.
Via l’utilisation d’une symbologie qui permet de plus facilement
interpréter/comprendre les changements de NDVI identifiés Confer section 5.3.3
page 76.
70
Dans l’analyse ci-dessus, ces pixels n’ont pas été exclus de l’analyse alors qu’ils auraient dû
l’être. En particulier, si un pixel présente une valeur comprise entre 251 et 255 dans l’une
des 2 images de NDVI ou dans les 2 images NDVI, il convient de lui attribuer une valeur
spécifique (nous lui attribuerons la valeur « Nodata » dans ce cas-ci) dans l’image de
« changement de NDVI » afin que l’on puisse facilement identifier ces pixels qui ne
traduisent pas un changement de NDVI.
Pour ce faire la méthode ci-dessous est proposée. Nombre d’autres méthodes pourraient
être appliquées. Cette méthode bug avec QGIS 3.10.0 mais fonctionne avec QGIS 3.12.
1 identifie les pixels pour lesquels les 2 images de janvier et septembre présentent
des valeurs proportionnelles au NDVI (DN < 251)
0 identifie les pixels pour lesquels une des 2 images de janvier et septembre, ou les 2
images, ne présentent pas des valeurs proportionnelles au NDVI (DN ≥ 251)
Ce raster binaire de masque peut être facilement obtenu via la « Calculatrice raster » de
QGIS :
Cliquez sur le menu « Raster > Raster Calculator… », ce qui ouvrira la « Calculatrice
Raster »
Paramétrez l’interface de la « Calculatrice Raster » comme indiqué à la Figure 39 en
nommant votre raster en sortie « RASTER_MASQUE_BINAIRE_0_1.tif ». La formule
utilisée est la suivante :
o "g2_BIOPAR_NDVI_QL_200601030000_AFRI_VGT_V1.3@1" < 251 AND
"g2_BIOPAR_NDVI_QL_200609030000_AFRI_VGT_V1.3@1" < 251
71
Le « AND » (en français : « ET ») est un opérateur logique booléen permettant de
réaliser une comparaison entre 2 valeurs (0 ou 1) selon la logique reprise dans la
Figure 38.
Le tableau ci-dessous reprend les 3 cas de figures possibles dans notre cas avec des
exemples de valeurs de pixels et les résultats qui en découlent.
72
Figure 39 : Création d’un raster de masque binaire 0 1, via la « Calculatrice raster » de QGIS, avec la
valeur de 1 qui identifiera les pixels présentant une valeur de NDVI valide dans les 2 images de
janvier et septembre, et la valeur de 0 qui identifiera les autres pixels
Figure 40 : Raster de masque binaire 0 1 identifiant en blanc (1) les pixels présentant une valeur de
NDVI valide dans les 2 images de janvier et septembre et en noir (0) les autres pixels
73
5.3.2.2. Masquage de l’image de changement de NDVI par le raster de
masque binaire 0 1
La deuxième étape consiste en l’application du raster de masque binaire 0 1 créé à l’étape
précédente sur le raster de « changement de NDVI » pour produire un nouveau raster de
changement de NDVI pour lequel :
Les pixels correspondant à du NDVI dans les 2 images de départ conserveront leur
valeur de changement de NDVI
Les pixels ne correspondant pas à du NDVI dans une des deux, ou dans les deux
images NDVI de départ seront identifiés par le code « Nodata ».
Cette étape peut être vue comme un masquage du raster de « changement de NDVI » par le
raster de masque binaire 0 1.
Cliquez sur le menu « Raster > Raster Calculator… », ce qui ouvrira la « Calculatrice
Raster »
Paramétrez l’interface de la « Calculatrice Raster » comme indiqué à la Figure 41 en
nommant votre raster en sortie « DIFFERENCE_SEPTEMBRE_JANVIER_MASQUE.tif »
La formule utilisée est la suivante : "DIFFERENCE_SEPTEMBRE_JANVIER@1" /
"RASTER_MASQUE_BINAIRE_0_1@1"
74
Figure 41 : Division du raster de changement par le raster binaire 0 1 dans la « Calculatrice raster »
de QGIS
Attention ! Dans le raster masqué, les pixels ne traduisant pas un changement de NDVI
prennent :
75
la valeur « Nodata » (en français : « sans donnée ») : vous pouvez le vérifier via
l’outil « Identifier des entités »
une couleur transparente : attention ! Si une autre couche est
activée/cochée/visible et positionnée en-dessous du raster masqué dans le panneau
« Couches » de QGIS, c’est cette autre couche qui sera visible à l’endroit des pixels
transparents du raster masqué. Pour visualiser correctement le raster masqué, il
vaut donc mieux décocher toutes les autres couches de votre projet QGIS.
Pour à la fois conserver l’image en dégradé de gris dans le projet QGIS et produire une
représentation colorée de cette image, il est possible de simplement dupliquer la couche de
changement dans le panneau « Couches », comme indiqué ci-dessous :
76
conditions de végétation. Typiquement un dégradé bicolor (noir à blanc par exemple)
ne conviendra pas.
Cliquez-droit sur l’image de changement dupliquée > Propriétés… > Symbologie >
Rendu des bandes raster > Type de rendu > Pseudo-couleur à bande unique
Choisissez un dégradé de couleur adapté pour représenter les changements de NDVI,
comme indiqué par exemple à la Figure 43 (dégradé de couleur « RdYlGn » par
exemple).
Figure 43 : Choix d’une symbologie adaptée (dégradé de couleur « RdYlGn ») pour visualiser
l’image de changement de NDVI.
Modifiez éventuellement cette palette de couleur selon votre souhait. Dans ce cas-ci
la couleur centrale par défaut a été remplacée par un jaune un peu plus vif afin de
mieux différencier les pixels présentant un changement de NDVI nul ou faible des
pixels ne représentant pas un changement de NDVI (« Nodata » transparent
apparaissant en blanc si le background est blanc). Pour ce faire :
77
o Cliquez-gauche sur la « Palette de couleurs » sélectionnée dans la fenêtre des
propriétés de la couche > Symbologie de la couche
Amélioration de contraste
Enfin, il peut être souhaitable de réaliser une amélioration de contraste (quelques exemples
sont donnés à la Figure 45) afin de mieux faire ressortir les différentes classes de
changement de NDVI. Ceci peut se faire via :
Clic-droit sur l’image de changement dupliquée > Propriétés… > Symbologie > Rendu
des bandes raster > Paramètres « Min » et « Max » (Figure 43) à définir de manière
symétrique autour de 0 si vous voulez que la couleur centrale de la palette
corresponde à un changement de NDVI nul.
78
Au plus les valeurs des bornes inférieure et supérieure (paramètres « Min » et « Max ») de la
« Palette de couleurs » seront resserrées autour de 0, au moins la représentation colorée
sera nuancée étant donné que tous les pixels ayant des valeurs inférieures/supérieures aux
paramètres « Min »/« Max » prendront la couleur correspondant à la borne
inférieure/supérieure de la « Palette de couleurs » respectivement, comme illustré à la
Figure 45. A vous de choisir la représentation qui vous semble être la plus adaptée à vos
objectifs de représentation et d ‘analyse.
Figure 45 : Image de changement du NDVI entre janvier 2006 et septembre 2006, en dégradé de
gris et dégradés colorés utilisant différentes améliorations de contraste via l’utilisation des
paramètres « Min » et « Max » des propriétés de symbologie de la couche. L’image en bas à droite
utilise comme background le raster de masque binaire 0 1 avec la couleur noire correspondant aux
pixels pour lesquels un changement de NDVI ne peut être calculé (nuages, etc).
L’image de changement affichée avec une palette de couleurs appropriée (Figure 45, images
colorées) est facilement interprétable et le lien avec les conditions climatiques évident,
comme déjà interprété dans la section 5.2 page 63.
79
6. Analyse de séries temporelles d’images
sur bases d’applications web
Temps
~ 1h30
approximatif
Prérequis Aucun.
Consignes de
Confer la « Contextualisation » reprise ci-après
départ
80
Contextualisation*
Production de médias mettant en évidence la déforestation à partir de séries
temporelles d’images satellites disponibles sur des applications web basées sur la
télédétection.
Un peuple d’indiens d’Amazonie, les Arumbayas, dans l’Etat de Rondônia au Brésil, voit
chaque jour leur forêt détruite un peu plus, destruction qui s’est accélérée récemment
suite à l’accession à la présidence du Brésil de Jair Bolsonaro le 1er janvier 2019.
Les Arumbayas aimeraient sensibiliser l’opinion nationale et internationale à la
problématique rencontrée sur leur territoire afin d’en stopper la destruction.
En tant que géographe fraichement diplômé, et sensibilisé à leur cause, vous proposez de
les aider via la création d’un ou plusieurs médias qui mettront en évidence la gravité et
l’urgence de la situation à partir de l’analyse d’images satellites.
Les Arumbayas ont identifié le périmètre de leur territoire grâce à un GPS reçu de la part
d’activistes de l’ONG « Treepeace ». Ce territoire est repris dans le shapefile
« Paraiso_verde.shp » et dans le fichier KML « PARAISO_VERDE.KML ».
Face à l’urgence de la situation, le délai de réalisation est assez court et vous décidez
d’utiliser les applications web existantes plutôt que de vous lancez dans un travail
fastidieux de téléchargement et traitement en local des images satellites.
Vous testerez les différentes applications présentées ci-dessous et comparerez leurs
qualités, défauts, avantages et inconvénients.
Vous pourrez montrer à la fois l’évolution historique (sur plusieurs années ou décennies)
et l’évolution très récente (derniers mois ou années récentes)!
Parmi les médias potentiels à produire :
Images satellites de la zone montrant la situation (capture d’écran,…)
Animation de type « Timelapse » montrant l’évolution de la situation
Graphique de l’évolution temporelle de la végétation
* Cette contextualisation s’inspire de faits réels mais ne se veut pas véridique.
81
6.1. « QuickMapServices » (QMS) via QGIS
A l’aide des indications ci-dessous, installez « QGIS » sur votre ordinateur si ce n’est pas
encore fait et utilisez l’extension « QuickMapServices » (QMS) pour consulter une image
satellite de très haute résolution spatiale sur la zone d’étude de « Paraiso Verde ». (~ 10
minutes)
QGIS n’est pas une application web mais un logiciel à installer sur votre ordinateur. Une fois
installé, il permet d’accéder facilement, dans une interface SIG, à une série de données en
lignes, dont des images satellites très hautes résolutions spatiales, via différents services
web dont le service « QuickMapServices » qui est utilisé dans cette section.
! QGIS doit être installé sur votre ordinateur pour cette section. QGIS est disponible
ici https://qgis.org/en/site/.
Si ce n’est pas encore fait, dans votre interface QGIS, installez l’extension
« QuickMapServices » via le Menu Extension > Installer / Gérer les extensions
Après l’installation de base de QMS, ajoutez les services supplémentaires de QMS
via le menu « Internet > QuickMapServices > Settings > More > Get contributed
pack »
Ajoutez le fond de carte “Google Satellite » à votre projet QGIS via le menu « Internet
> QuickMapServices > Google > Google Satellite »
Analysez visuellement le territoire de « Paraiso Verde ». Pensez-vous pouvoir
produire un média intéressant via l’interface QGIS ?
Avantages de QMS :
82
Accessible directement dans QGIS et donc facilement comparable avec des données
géographiques personnelles préalablement disponibles dans un projet QGIS
Désavantages de QMS :
Vue statique sans information sur la date des images satellites disponibles
Pas d’accès à une série temporelle d’images qui pourraient montrer l’évolution de
l’occupation du sol dans la zone
Figure 46 : Ajout du fond de carte « Google Satellite » via l’extension « QuickMapServices » dans
un projet QGIS et visualisation de la zone d’étude de « Paraiso Verde » dans l’Etat de Rondônia au
Brésil. Les « couloirs » de déforestation sont très visibles sous la forme de motifs linéaires gris-
brun, le vert correspondant à la forêt.
83
6.2. Google Earth Pro Desktop
A l’aide des indications ci-dessous, installez « Google Earth Pro » sur votre ordinateur si ce
n’est pas encore fait et utilisez cette application pour consulter les images satellites
d’archives disponibles sur la zone d’étude de « Paraiso Verde ». (~ 10 minutes)
« Google Earth Pro » permet d’afficher des images satellites d’archive, ce qui permet de se
faire une idée de la dynamique de l’occupation du sol dans une région.
! Google Earth Pro doit être installé sur votre ordinateur pour cette section : Google
Earth Pro est disponible ici :
https://www.google.be/intl/fr_ALL/earth/versions/#earth-pro.
Google Earth Pro nécessite une connexion à internet pour fonctionner correctement
et accéder aux images satellites d’archives notamment.
Par défaut, Google Earth Pro présente le même fond de carte que celui utilisé via l’extension
« QuickMapServices » dans QGIS ci-dessus.
Affichez l’outil « Afficher des images d’archive » est disponible comme illustré à la
Figure 47 :
Notez que les images d’archives disponibles à l’affichage varient en fonction du niveau de
zoom. A partir d’une échelle d’approximativement 1/80 000 en dé-zoomant (zoom out), des
images annuelles sans nuage sont disponibles de 1984 à 2016. Ce sont probablement
essentiellement des images LANDSAT (information non vérifiée). Sur la zone, la plus
ancienne image semble remonter à 1970. En zoomant à une échelle supérieure à 1/80 000
(zoom in), d’autres images, souvent de plus hautes résolutions spatiales, apparaissent dans
la liste des images d’archive.
84
Passez en revue les images disponibles sur la zone d’étude afin de vous faire une idée
de la dynamique d’occupation du sol dans la région (Figure 48).
Réponse : attention ! Les images affichées pour 1970 et 1975 datent en fait de 2001, même
lorsque l’on zoom sur une toute petite zone ! Il semble donc y avoir une erreur dans Google
Earth à ce niveau ! Confer série d’images ci-dessous !
1984 2016
Figure 48 : Images satellites d’archive sélectionnées pour quelques dates et affichées dans Google
Earth Pro sur la zone d’étude de « Paraiso Verde » dans l’Etat de Rondônia au Brésil. Les « couloirs
» de déforestation sont très visibles sous la forme de motifs linéaires gris-brun, le vert
correspondant à la forêt.
Pour information, une animation (image animée GIF) peut facilement être réalisée dans le
logiciel gratuit GIMP (https://www.gimp.org/) à partir d’une série d’images. Un exemple
d’animation a été créé avec 4 images de la zone (1984, 2001, 2008, 2016) et est disponible
dans le dossier « \3_SERIE_TEMPORELLE_APPLI_WEB\GIMP_ANIMATION_GIF\ » pour créer
ce type d’animation :
85
Ajouter et ordonner chronologiquement les images dans GIMP (la plus récente au-
dessus)
File > Export as > et donner un nom de fichier avec l’extension .GIF
Export > As animation + loop forever + paramétrer le temps d’affichage d’une image
de la série (par exemple : 500 ms) > Export
Le résultat est le fichier GIF
« \3_SERIE_TEMPORELLE_APPLI_WEB\GIMP_ANIMATION_GIF\GOOGLE EATH PRO
PARAISO VERDE 1984 2001 2008 2016.gif » qui s’ouvrira comme animation dans un
navigateur web.
86
6.3. Google EarthEngine Timelapse
“Google EarthEngine Timelapse” est une « vidéo mondiale zoomable » qui permet de
visualiser l’évolution de l’occupation du sol de la planète pendant ces 35 dernières années
(1984 - 2018) à l’aide d’une couverture mondiale annuelle sans nuage d’images satellites
(Landsat, Sentinel 2, ~15 millions d’images satellites utilisées). Cette application a été
réalisée dans le « Google Earth Engine’s cloud computing model ».
https://earthengine.google.com/timelapse/
Une série d’exemples bien connus d’évolution impressionnante de l’occupation du sol sont
présentés dans le panneau latéral de gauche, dont la déforestation massive au Brésil dans
l’Etat de Rondônia où se situe Paraiso Verde.
Explorez quelques exemples à votre guise, selon votre intérêt (déforestation, mine,
développement urbain, fonte de glacier, assèchement de la mer d’Aral, irrigation,
feux, modification du lit de fleuves, etc). Il faut parfois un peu de temps avant que les
images ne soient nettes.
En plus de ces exemples bien connus, cette vidéo est disponible pour tout endroit de la
planète.
87
Figure 49 : Aperçu de l’application “Google EarthEngine Timelapse” pour la zone d’étude de
« Paraiso Verde » dans l’Etat de Rondônia au Brésil. Les « couloirs » de déforestation sont très
visibles sous la forme de motifs linéaires gris-brun, le vert correspondant à la forêt.
88
6.4. Global Forest Change par « University of MARYLAND »
Un autre portail notable, aussi basé sur l’outil Google Earth Engine, est le portail « Global
Forest Change » développé par la « University of MARYLAND » (Figure 50), disponible ici :
http://earthenginepartners.appspot.com/science-2013-global-forest
Cette application résulte d’analyses de séries temporelles d’images satellites Landsat entre
2000 et 2018 et met en évidence l’étendue et les changements des forêts sur l’ensemble de
la planète.
Consultez les différentes couches d’information disponible via le menu « Data Products ».
Quelle(s) couche(s) d’information vous semblent la (les) plus pertinentes par rapport à la
problématique de cet exercice ?
Remarquez que la couche « Forest Loss Year (2018 Highlight) » permet de bien mettre en
évidence la déforestation 2018 qui semble inaugurer assez clairement un nouveau front de
déforestation dans la partie Nord de « Paraiso Verde ».
89
Figure 50 : Portail web « Global Forest Change » développé par « University of MARYLAND » pour
la zone d’étude de « Paraiso Verde » dans l’Etat de Rondônia au Brésil. Donnée affichée : « perte
en forêt » annuelle avec une mise en évidence de l’année 2018 en bleu clair.
90
6.5. EO Browser
EOBrowser est une application plus récente, développée par l’ESA (European Space Agency),
disponible ici (Figure 51):
https://www.sentinel-hub.com/explore/eobrowser
Accédez au portail via le bouton “Start exploring” (via le bouton « sign in », vous
risquez de rencontrer un problème de connexion)
91
“EO Browser combines a complete archive of Sentinel-1, Sentinel-2, Sentinel-3, ESA’s archive
of Landsat 5, 7 and 8, global coverage of Landsat 8, Envisat Meris, Proba-V and MODIS
products in one place.”
Enregistrez-vous pour avoir accès à toutes les fonctionnalités de cette interface (Sign
in > Sign up).
!! Pour faire une recherche d’images, attention, respectez les étapes suivantes :
92
Pour une image Sentinel-2 et une image Sentinel-1, explorez rapidement les
différents modes de visualisation disponibles (indices de végétation, compositions
fausse couleur, etc) (exemple à la Figure 52).
Réalisation d’un graphique de l’évolution temporelle du signal satellite sur la zone (Figure
53 et Figure 54).
Pour réaliser un graphique de l’évolution temporelle d’un produit, il faut utiliser une
« ROI » (Region Of Interest), soit dessiné manuellement (bouton « Draw area of interest » à
droite de l’interface), soit importé à partir d’un fichier KML par exemple.
Après avoir obtenu les résultats de recherche pour les images de type Sentinel-2 2A et choisi
le rendu « NDVI »:
Timelapse
93
Figure 52 : Aperçu de la zone d’étude avec le NDVI d’une image Sentinel-2
Figure 53 : Graphique de l’évolution temporelle (1 an) du NDVI de Sentinel-2 L2A sur la zone
d’étude, sans filtre de nuages : signal très bruyant en dents de scie.
94
Figure 54 : Graphique de l’évolution temporelle (1 an) du NDVI de Sentinel-2 L2A sur la zone
d’étude, AVEC filtre de nuages à 52 % : signal moins bruyant qu’à la figure précédente.
95
Figure 55 : Aperçu de la zone d’étude avec une image Sentinel-1
L’exemple repris dans les figures ci-dessous correspond aux coordonnées LAT -8.33 et LON -
62.88 entre le 15 et le 20/7/2019 avec le NDVI de Sentinel-2. Un graphique de l’évolution
temporelle du NDVI met clairement en évidence une cassure (en anglais : « break »)
correspondant à une forte diminution du signale NDVI entre ces 2 dates.
96
Figure 56 : Parcelle déforestée dans la zone d’étude telle qu’elle apparait sur des images NDVI de
Sentinel-2, avant (en haut) et après (en bas) déforestation.
97
Figure 57 : Graphique de l’évolution temporelle (6 mois) du NDVI de Sentinel-2 L2A sur la parcelle
déforestée, SANS (en haut) et AVEC (en bas) filtre de nuages: baisse soudaine du NDVI due à la
déforestation bien visible en juillet 2019.
98
6.6. Time series viewer SPOT-VGT & PROBA-V
A l’aide des indications ci-dessous, découvrez l’application « Time series viewer » et utilisez-
la sur la zone d’étude de « Paraiso Verde ». (~ 15 minutes)
Analyse d’une série temporelle d’images NDVI basse résolution via une application web
(temps : 15 minutes).
A l’aide du « Time series viewer » développé par le VITO et disponible ici (il faut copier-coller
le lien manuellement dans une page web pour accéder à l’application):
https://proba-v-mep.esa.int/applications/time-series-viewer/app/app.html,
Produits disponibles :
NDVI : issus des capteurs PROBAV et SPOT-VGT (S10 = synthèse décadaire (sur une
période de 10 jours) à partir d’images acquises quotidiennement)
4 produits CGLS (Copernicus Global Land Service) : FAPAR, LAI, FCOVER, Dry Matter
Productivity
Météo
99
Remarquez la possibilité de travailler avec des masques d’occupation du sol (« landcover »).
Figure 58 : Double saison de végétation annuelle observée pour le Punjab au Pakistan sur base du
NDVI PROBA-V (courbe verte) et précipitations (bleu)
100
Figure 59 : Différents types de saisons de végétation pour quelques régions du monde sur base du
NDVI PROBA-V (courbe verte) et précipitations
101
Figure 60 : Evolution latitudinale des types de saisons de végétation (courbes vertes) pour 3
régions du Bénin (Afrique de l’Ouest), depuis une double saison sur la côte dans le Sud jusqu’à une
saison unique dans le Nord en passant par une saison de type intermédiaire (double petit pic) dans
le centre. Précipitations en bleu.
102
Remarquez que la qualité des données disponibles via ce portail ne semble pas toujours être
garantie : confer la Figure 61.
Figure 61 : Exemple d’évolution du FAPAR pour une région du Ghana : il semble y avoir une erreur
dans les données puisque la moyenne historique (Long Term Average - LTA) passe sous le
minimum.
103
Figure 62 : Evolution du NDVI PROBA-V pour l’Etat de Rondônia (courbe verte) et pour la zone
d’étude Paraiso Verde (courbe mauve).
L’écart de NDVI entre les 2 capteurs SPOT VEGETATION et PROBA V pour 2014
La différence d’amplitude et de bruit de NDVI entre SPOT VGT ET PROBAV
104
6.7. Terrascope
Enregistrez-vous pour avoir accès à toutes les fonctionnalités de cette interface (Register
and Log in).
105
Figure 64 : Aperçu du « viewer » de la plateforme « Terrascope » avec les images Sentinel-2
disponibles pour le 28 mars 2020.
https://proba-v-mep.esa.int/applications/geo-viewer/app/
106
6.9. Système mondial d’information et d’alerte rapide sur
l’alimentation et l’agriculture (SMIAR)
Site web :
http://www.fao.org/giews/earthobservation/index.jsp?lang=fr
107
6.10. Fonte des glaces du pôle Nord
Exemple supplémentaire de monitoring qui peut être fait à partir d’une série temporelle
d’images satellites.
108
7. Manipulations de données
hyperspectrales dans QGIS
Prérequis Aucun.
109
7.1. Introduction
Dans cette section vous aurez brièvement l’occasion de manipuler quelques données de
type « hyperspectral ». Un signal est dit « hyperspectral » lorsqu’il est composé de
nombreuses bandes spectrales (généralement 100 ou plus), étroites (1 à 500 nm de largeur
par exemple) et contiguës (signal continu dans la gamme spectrale considérée) (exemple à la
Figure 68).
110
7.2. Manipulation de données hyperspectrales ponctuelles
(ASD) dans « ViewSpecPro »
Les données que vous allez utiliser ont été acquises sur des échantillons de sol agricoles
étudiés en laboratoire pour différents niveaux de teneurs en eau et d’ombrages.
(Remarque : pour les sols ombragés, l’énergie incidente correspond à la pleine lumière et
l’énergie réfléchie correspond à la lumière réfléchie par un sol partiellement ombragé.)
Le logiciel à utiliser est « ViewSpecPro ». Ce logiciel gratuit peut être téléchargé à l’adresse
web suivante : https://www.malvernpanalytical.com/fr/support/product-
support/software/ViewSpecProSoftwareInstall.html après inscription sur ce site.
Sélectionnez tous les fichiers et cliquez sur « View > Graph Data ».
Modifiez éventuellement la symbologie des courbes via un double-clic gauche sur le
graphique > Style
Les spectres correspondants aux fichiers sélectionnés s’affichent à l’écran (Figure 68).
112
7.3. Manipulation de données hyperspectrales images
(CHRIS-PROBA)
Figure 69 : Vue d’artiste du satellite PROBA (crédit image: Verhaert Space, ESA) (gauche) et
illustration d’une séquence de prises de vue multi-angulaire de PROBA-1 (crédit image : SSTL, ESA)
(droite). Source: https://directory.eoportal.org/web/eoportal/satellite-missions/p/proba-1.
Plus d’information sur le couple CHRIS-PROBA est disponible aux pages web suivantes :
https://earth.esa.int/web/guest/missions/esa-operational-eo-missions/proba, ESA
earthnet online (ESA = European Space Agency)
https://directory.eoportal.org/web/eoportal/satellite-missions/p/proba-1
L’image que vous allez utiliser se trouve dans le dossier
« \4_HYPERSPECTRAL\HYPERSPECTRAL_IMAGE_CHRIS_PROBA\» et est nommée «
CH100728_N.BIN ». Elle a été enregistrée le 28 juillet 2010 au-dessus d’une région agricole
au Sud-Ouest de Leipzig (Allemagne) en position « NADIR » (« N ») (capteur à la verticale du
lieu enregistré). L’extension « .BIN » signifie que l’image est codée en binaire.
A l’aide des indications ci-dessous, découvrez l’image CHRIS-PROBA via son ouverture dans
QGIS et la réalisation d’une ou plusieurs compositions vraies couleurs et éventuellement
une ou plusieurs compositions fausses couleurs.
114
o Remarquez que plusieurs combinaisons sont possibles étant donné que
plusieurs bandes spectrales correspondent à chacune des 3 couleurs Rouge,
Vert et Bleu (Figure 70).
Identifiez rapidement quelques classes d’occupation du sol dont de la végétation
cultivée (champs cultivés), de la terre (champs nus), de l’eau (lac) et remarquez la
présence importante de nuages sur l’image.
A l’aide des indications ci-dessous, constatez la variation des signatures spectrales pour
différentes occupation du sol et pour les nuages.
Plusieurs extensions de QGIS permettent d’extraire la signature spectrale d’un pixel ou d’un
groupe de pixels (via un polygone) d’une image, de manière interactive et son affichage sous
la forme d’un graphique (Figure 71).
Afin de consulter de manière interactive la signature spectrale (composée par les 62 bandes
spectrales de l’image CHRIS-PROBA) de différentes occupations du sol présentes au sein de
l’image :
115
Installez l’extension QGIS « Temporal/Spectral Profile Tool » qui, après installation,
est disponible sous la forme du bouton ou dans le menu « Extension > Profile
Tool > Temporal/Spectral Profile » (Menu Extension > Installer/Gérer les extensions >
Recherche et installation de l’extension)
(Remarque : une autre extension permettant également de réaliser un graphique de
Sol nu Végétation
Eau Nuage
Figure 71 : Signatures spectrales obtenues pour 3 occupations du sol et pour les nuages à partir de
l’image hyperspectrale CHRIS-PROBA composée de 62 bandes spectrales.
116
7.3.2.3. Graphique 2D de bandes spectrales 2 à 2
Le nuage de points du graphique des bandes spectrales 1 et 2 est très « aplati » sur la
« diagonale X = Y » car les valeurs de pixels de ces 2 bandes sont extrêmement proches ou
« extrêmement corrélées » (Figure 72, gauche). 2 bandes spectrales proches sont quelque
part assez redondantes étant donné qu’elles portent une information spectrale similaire.
Le nuage de points du graphique des bandes spectrales 1 et 42 est très « éclaté » car les
valeurs de pixels de ces 2 bandes sont « très peu corrélées » (la variation des valeurs de
pixels de la bande 1 ne correspond pas à une variation similaire des valeurs de pixels de la
bande 42) (Figure 72, droite). 2 bandes spectrales distantes ne sont pas redondantes car
elles portent chacune une information spectrale différente/particulière. 2 bandes spectrales
distantes seront donc plus riche en information et donc plus intéressantes pour caractériser
une occupation du sol donnée que 2 bandes spectrales proches.
117
Figure 72 : Graphiques de type « nuage de points » (anglais : « scatter plot ») pour l’ensemble de
l’image CHRIS-PROBA pour les bandes spectrales 1 et 2 (gauche) et 1 et 42 (droite), réalisés avec
l’extension « Raster Data Plotting » de QGIS
118
7.3.2.4. Cube hyperspectral
QGIS ne permet pas de réaliser un cube hyperspectral. Un exemple est simplement présenté
à la Figure 73.
ème
Figure 73 : Cube hyperspectral. La 3 dimension du cube (la profondeur) correspond aux longueurs d’ondes
du spectre électromagnétique. Illustration adaptée depuis : https://www.dronelogx.com/geospatial-
imagery-analysis.
119
7.3.2.5. Importance du « Fichier entête » ou « Header file »
Le « Fichier entête » (anglais « Header file ») est un petit fichier texte accompagnant une
image, qui contient une série d’information décrivant l’image, et qui permet à un logiciel de
lire et d’afficher correctement cette image.
120
8. Classification supervisée d’images
satellites multispectrales et analyse
diachronique : méthode simple et
méthode avancée (dans QGIS)
Il est préférable mais non obligatoire que l’utilisateur ait déjà quelques
Prérequis
notions de base en SIG (QGIS) et télédétection.
121
8.1. INFORMATIONS PREALABLES IMPORTANTES !
Cette section numéro 8 est assez longue (~140 pages au total) et peut être parcourue
selon 2 CHEMINS METHODOLOGIQUES correspondant à une :
METHODE SIMPLE de CLASSIFICATION SUPERVISÉE D’UNE SEULE image satellite (de
type Landsat-5) (+acquisition et prétraitement de l’image + représentation 3D +
analyse de séparabilité spectrale des classes + validation + mise en page
cartographique), SANS ANALYSE DIACHRONIQUE
o ATTENTION ! POUR LA METHODE SIMPLE :
Les titres de sections utiles à sa réalisation sont surlignés en vert
Les étapes méthodologiques utiles à sa réalisation sont encadrées
d’un halo vert dans le schéma méthodologique global de la Figure 75.
o Cette méthode simple est présentée sur un total approximatif de 90 pages
METHODE AVANCEE de CLASSIFICATION SUPERVISÉE DE PLUSIEURS images
satellites (de type Landsat-5 et Sentinel-2) reprenant les étapes de la méthode
simple ET EN ADDITION la réalisation de traitements supplémentaires (masques de
nuages, construction d’une image désennuagée, reclassification conditionnelle) et la
réalisation d’une ANALYSE DIACHRONIQUE (matrice de transition, graphiques
d’évolution du paysage, etc)
o L’ensemble des sous-sections de cette section 8 sont à considérer pour mener
à bien la méthode avancée (~140 pages) à l’exception de quelques sections
qui sont données à titre purement informatif.
La METHODE AVANCEE est volontairement complexe, ceci afin de montrer une large
gamme (mais très loin d’être exhaustive) des méthodes potentiellement applicables à ce
type d’analyse. Sachez cependant que nombre de manipulations peuvent se révéler
facultatives dans certains cas et que leur utilisation devra donc être adaptée en fonction de
vos données d’entrée, des objectifs de l’étude (classes d’occupation du sol d’intérêt,
précision souhaitée), de la difficulté de la classification de la zone d’intérêt, etc. Dans
certains cas, une classification supervisée d’images satellites ainsi qu’une analyse
diachronique peuvent donc être conduites beaucoup plus directement si désiré.
122
8.1.2. Traitements déjà réalisés
Certains traitements de données ont déjà été effectués pour vous et les données résultantes
sont mises à disposition dans les sous-dossiers du répertoire
«…\5_CLASSIFICATION_SUPERVISEE_ANALYSE_DIACHRONIQUE\».
Ces étapes, reprises sous la forme d’encadrés oranges dans la Figure 75, sont clairement
renseignées tout au long de la procédure et l’utilisateur pourra décider de :
Soit, pour gagner du temps : repartir des données « prétraitées » mises à disposition
Soit, pour pratiquer l’ensemble de la procédure et un meilleur apprentissage :
refaire l’ensemble des étapes, indépendamment des données mises à disposition.
8.2.1. Objectifs
Cet exercice a pour objectif principal de donner à l’utilisateur les outils méthodologiques
nécessaires que pour réaliser correctement :
Nombre de commentaires sont fournis sur les éléments auxquels il faut faire attention pour
réaliser une analyse méthodologiquement valide, sur les limitations de l’approche et sur les
pistes d’amélioration possible.
123
Etapes supplémentaires pour la METHODE
METHODE SIMPLE
AVANCEE
Acquisition des données à analyser via Acquisition des données à analyser via
recherche sur internet et téléchargement : recherche sur internet et téléchargement :
o 1 image satellite multispectrale o 2 images satellites multispectrales
Landsat-5 Sentinel-2
o Modèle Numérique de Terrain
(MNT) (pour la représentation 3D)
o Données auxiliaires telles que
données climatiques, shapefile et
informations contextuelles
Prétraitements sur l’image Landsat-5 Prétraitements sur les images Sentinel-2
o Conversion en réflectance, o Conversion en réflectance,
o Sélection des bandes spectrales, sélection des bandes, découpage
o Découpage spatiale, spatiale, rééchantillonnage
Prétraitements sur le raster d’altitude o Gestion des nuages : masque de
o Fusion, nuages, mosaïquage, etc
o Reprojection,
o Découpage spatiale,
o Rééchantillonnage
Observation des images en 2D et en 3D et
identification des classes d’occupation du
sol
Sélection des zones d’entrainement (ROI
de calibration)
Classification supervisée
Séparabilité spectrale des classes
Amélioration de la classification par
reclassification conditionnelle en fonction
de l’altitude et de la proximité aux rivières
Post-traitements (filtre spatiale,…)
Réalisation d’une matrice transition
Réalisation de graphiques (camembert et
histogramme) exprimant l’évolution de
l’occupation du sol
124
Figure 75 : Etapes méthodologiques de la section 8 de ce manuel concernant la classification supervisée
d’images satellites multispectrales et l’analyse diachronique.
125
8.2.3. Choix du logiciel : QGIS 3.10 (ou version ultérieure)
Le logiciel QGIS 3.10 (ou version ultérieure) a été sélectionné car il est gratuit et
permet de réaliser toutes les opérations de la méthode développée dans cette
section dans un seul environnement, ce qui est assez confortable.
En particulier, le plugin « Semi-Automatic Classification Plugin (SCP) » de QGIS
permet de réaliser un grand nombre d’étapes clefs.
QGIS 3.10 (et versions ultérieures) est disponible ici
https://www.qgis.org/fr/site/forusers/download.html
A télécharger et installer si ce n’est pas encore fait sur votre ordinateur
8.2.4. Contextualisation
La zone étudiée se situe dans la région de la ville de Kasese en Ouganda (Afrique de l’Est), à
la frontière Est de la République Démocratique du Congo (RDC) et contient de nombreuses
zones protégées (Figure 76), dont le « Parc national Rwenzori Mountains » dans la partie
Nord-Ouest (Figure 76 et Figure 77).
Figure 76 : Localisation de la zone d’étude (rectangle noir) : région de la ville de Kasese en Ouganda
(Afrique de l’Est)
Un aperçu rapide de la topographie de la zone d’étude peut être obtenu facilement via
l’application « Google Earth Pro » (Figure 77) déjà présentée à la section 3 page 31
(installation et utilisation de « Google Earth Pro »).
126
Figure 77 : Aperçu de la zone d’étude dans Google Earth Pro
Une idée de l’évolution de la couverture forestière dans la zone d’étude (Figure 78) est
disponible via l’application « Global Forest Change » de la « University of Maryland »
disponible via le lien ci-dessous et déjà présentée à la section 6.4 page 89.
https://earthenginepartners.appspot.com/science-2013-global-forest
127
Figure 78 : Aperçu de l’étendue, des pertes et des gains de la couverture forestière entre 2000 et
2018 dans la zone d’étude via l’application « Global Forest Change » de la « University of
Maryland ».
128
8.3. Préparation du projet QGIS
METHOD SIMPLE QGIS Desktop METHODE AVANCEE QGIS Desktop with GRASS
Figure 79 : Démarrage de « QGIS Desktop » ou « QGIS Desktop with GRASS » (Il est fortement
conseillé d’utiliser la dernière version de QGIS disponible et non pas la version 3.10).
Après avoir lancé la version de QGIS désirée et l’interface de QGIS ouverte, choisissez une
des 2 méthodes proposées ci-dessous:
129
« SRTM_FUSION_EPSG32635_CLIP_90m.tif » (raster)
« RIVIERES.shp » (vecteur) (Méthode avancée seulement)
« RIVIERES_TAMPON_60_m.shp » (vecteur) (Méthode avancée
seulement)
« RIVIERES_TAMPON_60_m_RASTER.tif » (raster) (Méthode
avancée seulement)
o Si ce n’est pas encore fait, installez l’extension QGIS
« QuickMapServices » et ajoutez le service « Google Satellite » (fond de
carte correspondant aux images Google satellite) dans votre projet QGIS,
comme décrit à la section 0 page 81.
o Organisez les différentes données disponibles correctement dans le
panneau « Couches » et adaptez éventuellement leur symbologie (« clic-
droit sur la couche > Propriétés… > Symbologie ».
o Sauvegardez ce projet QGIS dans votre répertoire de travail via le menu
« Projet > Enregistrer sous… » et nommez ce projet (exemple :
« Classification_votre_nom.qgz ».
130
Figure 80 : Aperçu du projet QGIS mis à disposition « CLASSIF_SUP_ANALYSE_DIA.qgz » mettant en
évidence la zone d’étude
Attention ! Les images satellites Landsat-5 et Sentinel-2 utilisées dans le cadre de cet
exercice pour la zone d’intérêt (ville de Kasese en Ouganda) sont toutes fournies, lors du
téléchargement sur internet, dans le SCR « EPSG32635 - WGS84 UTM Zone 35N - projeté ».
Cette information est disponible (vous pourrez le vérifier une fois les images satellites
acquises, confer section 8.5 page 142):
Soit via les fichiers de métadonnées accompagnant les images lorsqu‘elles sont
téléchargées de manière classique
131
o Par exemple, pour une image Landsat-5, on retrouve cette information dans
le fichier de type texte qui se termine par « MTL.txt » en fin de texte («
DATUM = "WGS84", ELLIPSOID = "WGS84", UTM_ZONE = 35 »).
Soit via un clic-droit sur les bandes spectrales des images une fois celles-ci ajoutées
dans QGIS > Propriétés… > Information > SCR
Nous adopterons donc, dans ce cas-ci, le SCR mentionné ci-dessus pour l’ensemble du
projet QGIS:
Pour le projet QGIS en tant que tel. Pour définir le SCR du projet QGIS :
o Cliquez sur le menu QGIS « Projet > Propriétés…> SCR » et choisissez le SCR
d’intérêt.
Pour les données géographiques images et vecteurs ajoutées dans ce projet. Note :
si vous deviez changer le SCR d’une donnée géographique qui présente un autre SCR
que le SCR d’intérêt, utilisez les fonctions de Reprojection suivantes dans QGIS :
o Pour les vecteurs : menu « Vecteur > Outils de gestion de données >
Reprojeter une couche… »
o Pour les rasters (images): menu « Raster > Projections > Projection (warp)… »
La plupart des indications concernant spécifiquement l’utilisation du plugin SCP seront ici
assez sommaires, dans un but de concision. Le lecteur intéressé pourra se référer au manuel
officiel de SCP pour plus d’explications (disponible, en anglais seulement, via les liens ci-
dessous) et en particulier aux tutoriaux officiels de SCP numéro 1 et 3.
https://fromgistors.blogspot.com/
https://fromgistors.blogspot.com/p/user-manual.html?spref=scp
Ce manuel contient notamment 3 tutoriaux dont :
o “Tutorial 1: Your First Land Cover Classification” qui explique comment
réaliser une classification supervisée d’une image satellite Landsat
o “Tutorial : Accuracy Assessment of a Land Cover Classification” qui explique
comment évaluer la précision d’une classification d’une image via la
132
production d’une matrice de confusion) (ce tutoriel n’est plus accessible
depuis la page web principale du manuel).
https://plugins.qgis.org/plugins/SemiAutomaticClassificationPlugin/
Après installation, les fonctionnalités du plugin sont disponibles via (Figure 81):
Le menu principal de QGIS « SCP » qui, lors de l’utilisation d’une fonction, ouvre
automatiquement une nouvelle interface appelée « Semi-Automatic Classification
Plugin »
2 barres d’outils : une principale et une secondaire
Un panneau latéral appelé « Menu SCP », qui peut être ajouté comme panneau près
du panneau « Couches » de QGIS pour plus de facilité, comme illustré à la Figure 81
(déplacez le panneau jusqu’à ce qu’il vienne se coller automatiquement au panneau
« Couches »).
Figure 81 : Menu, barres d’outils, panneau latéral et fenêtre principale du plugin SCP
133
8.4. Données auxiliaires : présentation, création,
acquisition et prétraitements
Shapefile « ZONE_ETUDE_KASESE.shp »
o Shapefile délimitant la zone d’étude, exprimé dans le Système de
Coordonnées de Référence (SCR) EPSG 32635
o Shapefile créé pour cet exercice
Shapefile « ne_50m_admin_0_countries_EPSG32635.shp »
o Shapefile des limites administratives des pays du monde, utile pour faciliter la
localisation relative de la zone d’étude,
o Shapefile téléchargé via le site web « Natural Earth »
https://www.naturalearthdata.com/downloads/ > Large scale data, 1:10m
(meilleure résolution parmi les 3 disponibles sur ce site web : 1 cm = 100 km)
> Cultural > Admin 0 – Countries > Download countries. A
dézipper/décompresser après téléchargement et avant l’ajout dans QGIS (si
vous vouliez le re-télécharger par vous-même).
Shapefile « RIVIERES.shp » (Figure 82) (méthode avancée seulement)
o Shapefile de quelques rivières principales dans la zone d’étude
o Shapefile créé par numérisation rapide dans QGIS sur base d’une photo-
interprétation du fond de carte « Google Satellite »
Shapefile « RIVIERES_TAMPON_60_m.shp » (Figure 82) (méthode avancée
seulement)
o Shapefile d’une zone tampon de 60 mètres de part et d’autre des rivières du
fichier « RIVIERES.shp »
o Shapefile créé à l’aide de la fonction QGIS « Vecteur > Outils de géotraitement
> Tampon… »
134
Figure 82 : Aperçu des rivières numérisées et de leurs zones tampon de 60 m
Un MNT est un raster (image) dont la valeur des cellules ou pixels correspond à une
information sur l’altitude.
Un MNT à 90 mètres de résolution spatiale découpé sur la zone d’étude et exprimé dans le
SCR EPSG32635 est disponible dans le dossier
« \5_CLASSIFICATION_SUPERVISEE_ANALYSE_DIACHRONIQUE\DATA\DONNEES_AUXILIARES
» sous le nom de
« SRTM_FUSION_EPSG32635_CLIP_90m.tif ».
135
Vous pouvez :
ATTENTION ! Des informations détaillées décrivant comment télécharger des données sur ce
site web sont données dans la section 8.5.2.3 page 150 ci-dessous, avec pour exemple le
téléchargement d’une image Landsat-5. Le principe est le même pour le téléchargement du
MNT.
Le MNT mondial (donnée très lourde) est mis à disposition au téléchargement sous la forme
de très nombreux petits morceaux de MNT d’approximativement 111 km * 111 km sur la
zone d’étude, appelés « tiles » (données plus légères de ~25 Mo chacune (Figure 85)), afin de
faciliter le téléchargement de MNT pour des régions particulières. La zone d’étude de cet
exercice-ci est couverte par 4 « tiles » comme illustré aux Figure 83 et Figure 84. La
résolution spatiale de ce MNT, c’est-à-dire la taille des pixels, est de 1 arc-seconde de côté (1
seconde d’arc = 1/36000 de degré = 0.00028 degré) ce qui correspond à ~30 mètres.
136
Figure 83 : Aperçu du positionnement des 4 « tiles » SRTM (MNT) couvrant la zone d’étude sur le
site web https://earthexplorer.usgs.gov/
Figure 84 : Aperçu des 4 « tiles » SRTM (MNT) couvrant la zone d’étude sur le site web
https://earthexplorer.usgs.gov/
137
Figure 85 : Aperçu des 4 « tiles » SRTM (MNT) couvrant la zone d’étude, dans l’explorateur
Windows après téléchargement
Après téléchargement, les opérations reprises ci-dessous peuvent être réalisées afin
d’obtenir un seul fichier raster d’altitude couvrant toute la zone d’étude et exprimé dans le
même système de coordonnées que les images satellites utilisées dans le cadre de cet
exercice-ci.
Ajoutez les 4 « tiles » MNT dans votre projet QGIS (via le menu d’ajout d’une coucher
raster)
Vous remarquerez que le MNT au Nord-Ouest de la zone d’étude (Figure 84) semble
présenter quelques mauvaises données en montagne. En effet, au milieu du massif
montagneux du « Parc national Rwenzori Mountains » (confer section 8.2.4 page 126), qui
est identifié par des couleurs blanches (le blanc correspondant aux altitudes les plus
élevées), quelques taches noires présentant une altitude de 0 sont présentes (le noir
correspondant aux altitudes les plus faibles). Nous ignorerons ce problème potentiel dans le
cadre de cet exercice. Remarque : ce problème persiste (un peu amoindri) dans la version
« 3Arc void filled » du MNT.
Afin d’obtenir un seul fichier MNT plus facilement manipulable que 4 « tiles » séparées, il
faut fusionner les 4 « tiles ». Pour ce faire :
Afin d’obtenir un MNT dans le même Système de Coordonnées de Référence (SCR) que celui
qui a été choisi comme SCR à utiliser dans le cadre de cet exercice (pour rappel, c’est le SCR
138
EPSG 32635, confer section 8.3.2 page 131), il faut reprojeter le MNT fusionné car celui-ci se
trouve par défaut dans le SCR « EPSG:4326 - WGS 84 - Géographique » (pour le vérifier :
cliquez-droit sur la couche MNT fusionnée > Propriétés > Information > SCR).
Afin d’obtenir un MNT qui ne couvre que la zone d’étude et qui sera donc plus léger et plus
facilement manipulable qu’un MNT qui couvrirait une plus grande zone, il convient de
découper le MNT sur la zone d’étude (Figure 86). Pour ce faire :
Menu QGIS: Raster > Extraction > Découper un raster selon une emprise…
Couche source : le MNT reprojeté résultat de l’étape précédente
Etendue de découpage (xmin, xmax, ymin, ymax) :
o > Utiliser l’emprise de la couche… : utilisez le shapefile identifiant la zone
d’étude « ZONE_ETUDE_KASESE.shp », exprimé lui-aussi dans le SCR EPSG
32635. Attention, il faut que ce shapefile ait une information « emprise »
dans ses métadonnées :
Cette information sur l’emprise est disponible via un clic droit sur
cette couche dans le panneau Couches > Propriétés… > Information >
Emprise.
Si cette information n’est pas présente, possibilité de créer cette
métadonnée via un clic droit sur cette couche dans le panneau
Couches > Propriétés… > Source > Mettre à jour l’emprise.
o Si vous n’avez pas de fichier shapefile identifiant votre zone d’étude :
Soit créez un fichier shapefile de type polygone dans le SCR d’intérêt
et éditez ce shapefile en y ajoutant un polygone qui couvre votre zone
d’étude (méthode non détaillée dans ce manuel mais disponible dans
le manuel QGIS disponible ici :
https://orbi.uliege.be/handle/2268/190559)
Soit spécifiez manuellement les valeurs des coordonnées
géographiques définissant l’étendue de la zone d’étude comme
indiqué : xmin, xmax, ymin, ymax où x fait référence à la longitude et y
à la latitude, coordonnées exprimées dans le SCR et ses unités de la
couche à découper, c’est-à-dire du MNT reprojeté dans ce cas-ci. Pour
139
connaitre ces coordonnées, autrement que via le calcul de l’emprise
proposé ci-dessus, positionnez le curseur de la souris aux extrémités x
et y de la zone d’étude : les coordonnées du pointeur de la souris
s’affichent en bas de l’interface QGIS dans la boite « Coordonnée ».
Ces coordonnées s’affichent dans le SCR du projet QGIS (modifiable via
le menu Projet > Propriétés… > SCR, choisir dans ce cas-ci EPSG 32635)
et dans les unités choisies via le menu QGIS « Projet > Propriétés… >
Général > Affichage des coordonnées > Afficher les coordonnées
avec » : choisir dans ce cas-ci « Unités de carte (mètres) ». De cette
manière vous obtenez les coordonnées du pointeur de la souris dans
le SCR du projet QGIS, qui est le même que le SCR du MNT reprojeté,
et dans ses unités (des mètres).
Découpé (étendue) : choisir l’option « Enregistrer vers un fichier… » puis le répertoire
de sortie (votre répertoire de travail) et donner un nom (par exemple : «
SRTM_FUSION_EPSG32635_CLIP.tif »)
Figure 86 : 4 « tiles » SRTM (MNT) (avec symbologie « RdYlGn ») et étendue de la zone d’étude
(rectangle noir) (gauche) et MNT après fusion et découpage sur la zone d’étude (avec symbologie
« wiki-schwarzwald-cont »*) (droite)
Afin d’alléger encore plus le MNT à utiliser, et donc de le rendre encore plus
facilement/rapidement manipulable, il est possible de le rééchantillonner à une taille de
140
pixel plus faible. Dans ce cas-ci le MNT d’origine à une taille de pixel de 30 mètres de côté, ce
qui n’est pas nécessaire dans le cadre de cet exercice. Nous allons donc rééchantillonner le
MNT à une taille de pixel de 90 mètres de côté, ce qui devrait diviser par 9 le poids
numérique du MNT (il y a 9 cellules de 30*30 m dans une cellule de 90*90m). Pour ce faire :
Pour ce faire, il faut un raster indiquant la proximité aux rivières. C’est le rôle du raster
« RIVIERES_TAMPON_60_m_RASTER.tif »
o Ce raster prend une valeur de 1 dans les zones tampon du fichier
« RIVIERES_TAMPON_60_m.shp » (à 60m de part et d’autre des rivières
numérisées, confer section 8.4.1 page 134) et une valeur de « NoData »
ailleurs.
o Ce raster a été créé via la fonction « Raster > Conversion > Rastérisation
(vecteur vers raster)… » avec le paramétrage suivant :
Couche en entrée : « RIVIERES_TAMPON_60_m.shp »
Valeur fixe : 1
Unité du raster résultat : Unités géoréférencée
Largeur/Résolution horizontale : 30
Hauteur/Résolution verticale : 30
Emprise du résultat (xmin, xmax, ymin, ymax) : utiliser l’emprise de la
couche : sélectionner le shapefile délimitant la zone d’étude
Rasterisé : choisir l’option « Enregistrer vers un fichier… » puis le
répertoire de sortie (votre répertoire de travail) et donner un nom
141
8.5. Acquisition des images satellites
8.5.1. Choix des dates des images satellites, notamment pour une analyse
diachronique
Afin d’augmenter la comparabilité des images analysées, il faudra donc veiller à ce que
celles-ci aient été acquises :
Ces considérations sont bien sûr à nuancer en fonction de vos objectifs (intéressé ou pas par
la classe végétation,…) et de la zone étudiée (impact saisonnier plus ou moins important :
régions tempérées ou régions sahéliennes vs régions tropicales).
142
distinguables à certaines périodes de l’année alors qu’elles le seront parfaitement à d’autres
périodes. Par exemple :
Imaginons que nous voulons réaliser l’étude en comparant les périodes suivantes :
Ces périodes sont ici très larges (5-10 ans), ceci afin de faciliter notre recherche d’images
satellites de qualité (sans ou avec peu de nuages). En fonction des caractéristiques
météorologiques de votre zone d’étude, vous pourrez adapter la largeur de la fenêtre
temporelle de recherche.
143
8.5.2.1.1. Recherche et identification des capteurs satellites
Dans un premier temps il faut trouver une / des images satellites correspondant à la période
d’intérêt / à chacune des 2 périodes d’intérêt. Pour ce faire il est indiqué de se renseigner
sur :
les capteurs satellites actifs durant ces périodes et leurs caractéristiques techniques,
l’accessibilité des images produites : gratuité ?, prix ?, restriction d’utilisation ?.
Une recherche exhaustive ne sera pas menée ici et notre choix se portera sur 1 / 2 capteurs
multi-spectraux de haute résolution spatiale parmi les plus connus et disponibles pour la
période / chacune des 2 périodes. Nous choisirons :
Une image Landsat-5 (satellite américain) comme image « historique ». Les images
Landsat-5 sont disponibles entre 1984 et 2013. (Méthode simple et méthode
avancée)
Une image Sentinel-2 (satellite européen) comme image récente. Les images
Sentinel-2 sont disponibles depuis 2015 jusqu’actuellement. (Méthode avancée
uniquement)
Un aperçu des périodes couvertes par les différents satellites Landsat 1 à 9 est disponible
ici :
https://landsat.gsfc.nasa.gov/a-landsat-timeline/
o Landsat-5 :
o https://www.usgs.gov/land-resources/nli/landsat/landsat-5?qt-
science_support_page_related_con=0#qt-science_support_page_related_con
o Sentinel-2 :
o https://sentinels.copernicus.eu/web/sentinel/missions/sentinel-2
D’autres capteurs auraient bien sûr pu être mobilisés. Pour information, des listes de
capteurs satellites sont disponibles ici :
https://webapps.itc.utwente.nl/sensor/default.aspx?view=allsensors
https://www.satimagingcorp.com/satellite-sensors/
https://www.indexdatabase.de/db/s.php
144
Le diagramme climatique de la région de Kasese est repris à la Figure 87. Cette figure nous
apprend que la température est constante tout au long de l’année mais qu’il y a 2 saisons
de précipitations plus importantes, entre mars et mai et entre août et novembre, ce qui
peut avoir un impact sur les conditions de végétation. Une analyse diachronique
interannuelle évitera donc par exemple de se baser sur une comparaison d’une image de
janvier (saison plus sèche) avec une image d’octobre (fin de saison pluvieuse) car ce type de
comparaison serait fortement impacté par l’évolution saisonnière des conditions de
végétations elle-même liée aux conditions de précipitations.
Une série temporelle multi-annuelle des précipitations pour la région de Kasese sur la
période 1981-2019 est présentée à la Figure 88 sous la forme de déviations par rapport à la
moyenne des précipitations totales annuelles. Ce graphique a été obtenu très facilement à
partir de l’application « Climate Engine » (https://clim-engine.appspot.com/climateEngine)
en utilisant les données de précipitations « CHIRPS - Pentad ». Ce graphique permet
d’identifier les années plutôt sèches, normales ou pluvieuses, ce qui peut constituer une
première observation intéressante. Notez cependant que, dans ce cas-ci, étant donné la
double saison de précipitations annuelles dans la zone d’étude, ce graphique, présentant des
totaux annuels, n’est pas extrêmement pertinent car il ne reflète en rien la possible variation
de la répartition des précipitations au cours de l’année et en particulier pendant les
quelques mois précédant l’acquisition des images satellites qui seront utilisées pour réaliser
l’analyse diachronique. Par exemple, une année pourrait être considérée comme plutôt
sèche alors que le déficit pluviométrique serait dû uniquement aux mois postérieurs à la
date d’acquisition de l’image satellite à utiliser. Remarquez que généralement les
températures conditionnent moins le développement de la végétation que les
précipitations. Notez enfin qu’il est toujours bon de questionner la qualité et la fiabilité des
données météorologiques utilisées.
145
Figure 88 : Déviations par rapport à la moyenne des précipitations totales annuelles dans la région
de Kasese, obtenues à partir de l’application « Climate Engine » utilisant les données de
précipitations « CHIRPS - Pentad » sur la période 1981-2019. Source : https://clim-
engine.appspot.com/climateEngine
Un graphique des précipitations mensuelles des mois précédant la date d’acquisition des 2
images satellites à comparer (Figure 89) est donc plus pertinent pour s’assurer que les
conditions météorologiques de ces années ne soient pas un facteur de variation trop
important de l’occupation du sol. Par exemple, dans le cas de la Figure 89-gauche, une
analyse diachronique basée sur des images satellites acquises aux alentours de la période
correspondant à la mesure 35 (confer les graphiques) pour les années 1984 et 1996 pourrait
présenter un biais de l’occupation du sol lié au fait que les semaines précédant cette période
(cercle noir sur le graphique) présente des précipitations relativement différentes entre
1984 (plus sec) et 1996 (plus humide), bien que cet exemple-ci ne soit pas un exemple très
marqué.
146
Figure 89 : Comparaisons des précipitations dans la région de Kasese pour différents couples
d’années : à gauche, 2 années relativement contrastées en termes de total de précipitations
annuels, et à droite, les 2 années considérées pour l’analyse diachronique menée dans le cadre de
cet exercice-ci. Source : données extraites de « Climate Engine »
(https://app.climateengine.org/climateEngine) et mises en graphique dans Excel.
Enfin, l’analyse d’une série temporelle d’un indice de végétation, tel quel le NDVI, le FAPAR,
etc, peut nous renseigner sur la saisonnalité des conditions de végétation, et donc, de
manière indirecte, sur les conditions météorologiques, que ce soit de manière générale
(comportement moyen) ou pour des années particulières que nous voudrions comparer. La
Figure 90 montre l’évolution de l’indice de végétation fAPAR tel que calculé à partir des
capteurs SPOT-VEGETATION et PROBA-V pour la région de Kasese entre 1999 et 2020. Cette
figure a été produite à partir de l’application web présentée à la section « 6.6 Time series
viewer SPOT-VGT & PROBA-V » page 99. Cette figure nous apprend que la région de Kasese
est caractérisée par une double saison de végétation annuelle avec la saison d’avril-mai
présentant généralement un pic de fAPAR très légèrement supérieur à celui de la saison de
septembre à novembre. Cette évolution trouve son origine dans la variation intra-annuelle
des conditions de précipitations comme discuté ci-dessus et illustré à la Figure 87. Cette
application ne permet cependant pas de remonter avant 1999, ce qui ne nous permet pas
d’évaluer la comparabilité météorologique (indirectement) d’une année dans la période
1985-1995 avec une année dans la période 2015-2019. Notez enfin les 2 limitations
suivantes à l’utilisation de ce type d’information:
147
Figure 90 : Aperçu de l’évolution de l’indice de végétation fAPAR (des capteurs SPOT-VEGETATION
et PROBA-V) pour la région de Kasese entre 1999 et 2020. Source : https://proba-v-
mep.esa.int/applications/time-series-viewer/app/app.html
l’impact potentiel que ces conditions peuvent avoir OU PAS sur les classes
d’occupation du sol dont vous voulez étudier l’évolution,
la disponibilité des images satellite pour les périodes que vous désirez comparer.
148
8.5.2.2. Remarque méthodologique sur la recherche et le téléchargement
des images satellites : site web vs plugin SCP de QGIS
Notez que pour certains prétraitements à faire sur les images avec le plugin SCP (correction
atmosphérique, création d’un jeu de bandes, sélection des bandes d’intérêt,
rééchantillonnage automatique, etc)
Si le téléchargement est effectué via SCP dans QGIS, ces prétraitements peuvent se
faire automatiquement dans la foulée du téléchargement (confer section 8.5.4 ci-
dessous)
Si le téléchargement est effectué indépendamment de SCP, via les sites web de
référence, ces prétraitements devront constituer une étape de traitement à part
entière (confer section 8.6 page 176).
Attention ! Dans le cadre de ce manuel, voici la méthode proposée, afin de ne pas perdre
trop de temps :
Recherche des images Landsat-5 et Sentinel-2 d’intérêt: cette recherche se fera via
les sites web de référence qui sont plus conviviaux que le plugin SCP (confer les 2
sections 8.5.2.3 et 8.5.2.4 ci-dessous), puis pour
Landsat-5 (méthode simple et méthode avancée): téléchargement et prétraitements
se feront via le plugin SCP dans QGIS, car, l’image n’étant pas trop lourde, cette étape
peut se faire en quelques minutes et cela permet d’explorer ces fonctionnalités de
SCP
Sentinel-2 (méthode avancée seulement): les images déjà téléchargées et déjà
prétraitées sont directement mises à disposition des étudiants car ces étapes peuvent
prendre un certain temps (10-15 minutes par image Sentinel-2). Ces images sont
disponibles dans le dossier « …\DATA\IMAGES\SENTINEL_2\ » et commencent par le
préfixe « 30m_clip_RT… ».
o Les prétraitements préalablement faits sur les images Sentinel-2 sont:
Téléchargement
Correction DOS1 dans SCP
Découpage spatial sur la zone d’étude
Sélection des bandes 2, 3, 4, 5, 8, 11 (430 Mo pour 6 bandes)
Rééchantillonnage des bandes à 30 m (108 Mo pour 6 bandes)
Dans les sections suivantes, des consignes reprises sous le titre « ! Méthode proposée dans le
cadre de ce manuel » vous indiqueront la marche à suivre.
149
téléchargement et prétraitements des images satellites, via site web de référence OU via le
plugin SCP, sont données, à toutes fins utiles, dans les sections suivantes :
Recherche et téléchargement des images satellites via les sites web de référence :
sections 8.5.2.3 et 8.5.2.4 pages 150 et 158 respectivement
Téléchargement des images avec prétraitements dans la foulée via le plugin SCP
dans QGIS : section 8.5.4 page 166
Prétraitements des images satellites non couplés au téléchargement via le plugin
SCP dans QGIS : section 8.6 page 176
Le site web de référence pour rechercher et télécharger une image Landsat est le suivant :
150
Figure 91 : Site web https://earthexplorer.usgs.gov/ utilisé ici pour la recherche des images
Landsat
La recherche d’une image satellite sur ce site peut se faire librement, sans création de
compte.
Le téléchargement d’une image satellite à partir de ce site, que ce soit directement à partir
du site web OU via le plugin SCP, nécessite par contre que vous :
Vous créiez un compte via le « EROS Registration System (ERS) » disponible ici :
https://ers.cr.usgs.gov/, et ensuite, que vous (Figure 92)
Vous loguiez sur le site https://earthexplorer.usgs.gov/ à l’aide de vos identifiant et
mot de passe (Figure 91).
151
Figure 92 : Page du site web https://ers.cr.usgs.gov/ permettant de créer un compte personnel
nécessaire pour pouvoir télécharger les images satellites Landsat
Les critères à utiliser pour la recherche d’une image satellite Landsat-5 sur le site web «
Earth Explorer » sont dans ce cas-ci (Figure 93):
152
WRS_Row : 60
o Scene Cloud cover : < 30 %. Ceci permettra de n’obtenir que les images
présentant une couverture nuageuse de moins de 30% et de cette manière ne
pas avoir trop d’images trop ennuagées, et donc inutilisables, dans les
résultats de la requête.
Ces critères de recherche résultent en 7 images Landsat-5 dont la plus intéressante est
l’image du 7 août 1987 qui présente très peu de nuage sur la zone d’étude (Figure 94). Un
aperçu géolocalisé des images est disponible via le bouton « Show Browse Overlay »
153
Figure 94 : Résultat de la recherche d’une image Landsat-5 sur la zone d’étude sur le site
https://earthexplorer.usgs.gov/
! Méthode proposée dans le cadre de ce manuel : une fois l’image d’intérêt trouvée,
prenez note de ses références (type de capteur (Landsat4-5 TM), date (7 août 1987), lieu
(Kasese)) et passez à la section 8.5.4 page 166 de ce manuel afin de faire le téléchargement
de l’image via le plugin SCP dans QGIS.
! Méthode proposée dans le cadre de ce manuel : cette section est donnée à titre
d’information. A ne pas faire.
Le téléchargement d’une image satellite ne sera possible que si vous vous loguez sur le site
https://earthexplorer.usgs.gov/ après vous être créé un compte via le « EROS Registration
System (ERS) » disponible ici https://ers.cr.usgs.gov/, comme déjà indiqué plus haut en
début de section 8.5.2.3.1 page 150.
Une fois l’image d’intérêt trouvée, le téléchargement direct s’effectue via le bouton
154
Figure 95 : Options de téléchargement d’une image Landsat-5 sur le site web
https://earthexplorer.usgs.gov/.
Les images sont téléchargées sous la forme d’un dossier doublement compressé « .tar.gz »
(Figure 96, 2ème ligne) qu’il faudra donc décompresser 2 fois avant de pouvoir les ouvrir dans
QGIS :
La 1ère décompression produira le dossier avec l’extension « .tar » (Figure 96, 1ère
ligne)
La 2ème décompression rendra disponibles tous les fichiers contenus dans ce dossier,
à savoir :
o Les bandes spectrales B1 à B7
o La bande « Quality Band » « BQA » renseignant sur la qualité de chaque pixel
(confer section 8.6.12.2 page 199)
o Le fichier de métadonnées « …MTL.txt »
o Un aperçu de l’image « …VER.jpg »
155
Figure 96 : Dossier « .tar.gz » d’une image Landsat-5 tel que téléchargé (2ème ligne) ET dossier
« .tar » obtenu après la première décompression (1ère ligne) ET fichiers contenus dans ce dossier
visibles après la 2ème décompression. Parmi ces fichiers, notamment, les bandes spectrales B1 à B7,
la « Quality Band » « BQA », les métadonnées « …MTL.txt » et un aperçu de l’image « …VER.jpg »
Remarques :
156
Figure 97 : Fenêtre « Options des dossiers » d’une fenêtre d’exploration Windows permettant de
faire apparaitre les extensions de fichiers de types connus en décochant l’option « Cacher les
extensions de types de fichier connu »
157
8.5.2.4. Recherche et téléchargement des images Sentinel-2 via le site web
« SciHub »
Le site web de référence pour rechercher et télécharger une image Sentinel-2 est le
suivant (Figure 98) :
https://scihub.copernicus.eu/dhus/#/home
158
Figure 98 : Page d’accueil du site web « Copernicus Open Access Hub » disponible ici
https://scihub.copernicus.eu/dhus/#/home » et permettant l’accès aux images Sentinel-2
notamment.
La recherche d’une image satellite sur ce site peut se faire librement, sans création de
compte.
Le téléchargement d’une image satellite à partir de ce site, que ce soit directement à partir
du site web OU via le plugin SCP, nécessite par contre que vous :
Vous créiez un compte via le bouton de login > Sign up, et ensuite que vous
Vous loguiez via le bouton de login > LOGIN à l’aide de vos identifiant et mot de
passe.
Les critères à utiliser pour la recherche d’une image satellite Sentinel-2 sur le site web «
SciHub » sont dans ce cas-ci (Figure 99):
Lieu : le lieu d’intérêt doit être identifié en se déplaçant dans l’interface via la
159
La période et le type de produit (satellite) doivent être paramétrés via le panneau
latéral de gauche affichable via le bouton . Dans ce cas-ci, les critères suivants ont
été utilisés :
o Période : entre le 1er juin 2018 et le 31 août 2018
o Type de produit : Mission : Sentinel-2
Figure 99 : Site web « Copernicus Open Access Hub » paramétré pour rechercher une image
Sentinel-2 dans la zone d’étude pour la période d’intérêt. Panneau latéral de gauche = critères de
recherche d’une image (période, produit (satellite, etc), couverture nuageuse, etc)
160
Figure 100 : Résultats de la recherche d’image Sentinel-2 dans le « Copernicus Open Access Hub »
La liste des images résultantes peut être explorée à la recherche de la meilleure image.
Des informations détaillées sur une image sont disponibles via le bouton
(Figure 101).
161
Figure 101 : Aperçu de l’image Sentinel-2 du 07 juillet 2018
Face à la difficulté de trouver des images satellites présentant une couverture nuageuse
suffisamment faible que pour permettre une analyse de qualité suffisante, il peut être
judicieux de procéder à un mosaïquage de 2 ou plusieurs images, chacune partiellement
ennuagées, et acquises dans un intervalle de temps relativement restreint (à moduler en
fonction des objectifs de l’étude) (Figure 102), afin de reconstruire une image
« désennuagée » à partir des zones non-ennuagées de chacune des images. Cette méthode
de « mosaïquage temporel » est proposée dans la section 8.6.9 page 187.
Les 2 images Sentinel-2 qui seront utilisées à cette fin dans le cadre de cet exercice (Figure
102) sont les images du :
12 juin 2018
07 juillet 2018
162
Image Sentinel-2 du 12 juin 2018 Image Sentinel-2 du 7 juillet 2018
Figure 102 : Aperçu des 2 images Sentinel-2, toutes 2 partiellement ennuagées, utilisées dans le
cadre de cet exercice
Une fois identifiée, une image peut être téléchargée via le bouton .
Les images sont téléchargées sous la forme de dossiers compressés « .zip » qu’il convient de
décompresser avant l’intégration de l’image dans QGIS. Ce dossier contient un nombre
important de sous-dossiers et fichiers et il ne sera pas détaillé ici.
163
Landsat-5, 7 août 1987 Sentinel-2, 12 juin 2018 Sentinel-2, 7 juillet 2018
Figure 103 : les 3 images satellites retenues dans le cadre de l’analyse proposée dans cet exercice.
Seule une partie de l’image Landsat-5 est montrée.
Les caractéristiques techniques les plus importantes des 3 images utilisées sont reprises
dans les Figure 104, Figure 105 et Figure 106.
Numéro de la bande Longueur d’onde (µm) Nom de la bande Résolution spatiale (m)
Band 1 0.45-0.52 Blue 30
Band 2 0.52-0.60 Green 30
Band 3 0.63-0.69 Red 30
Band 4 0.76-0.90 Near Infrared (NIR) 30
Band 5 1.55-1.75 Shortwave Infrared (SWIR) 1 30
Band 6 10.40-12.50 Thermal 120 (30)*
Band 7 2.08-2.35 Shortwave Infrared (SWIR) 2 30
Figure 104 : Caractéristiques des bandes spectrales des images Landsat-5 Thematic Mapper
™. (Source : https ://www.usgs.gov/faqs/what-are-band-designations-landsat-satellites ?qt-news_science_products=0#qt-
news_science_products)
* La bande TM 6 a été acquise avec une résolution spatiale de 120 mètres, mais elle a été
rééchantillonnée à 30 mètres de résolution spatiale.
Numéro de la bande Longueur d’onde (µm) Nom de la bande Résolution spatiale (m)
Band 1 443 Coastal aerosol 60
Band 2 490 Blue 10
Band 3 560 Green 10
Band 4 665 Red 10
Band 5 705 Vegetation red edge 20
Band 6 740 Vegetation red edge 20
Band 7 783 Vegetation red edge 20
Band 8 842 NIR 10
Band 8A 865 Narrow NIR 20
Band 9 940 Water vapour 60
Band 10 1375 SWIR – Cirrus 60
Band 11 1610 SWIR 20
Band 12 2190 SWIR 20
Figure 105 : Caractéristiques des bandes spectrales des images Sentinel-2 (Source :
https ://sentinel.esa.int/web/sentinel/user-guides/sentinel-2-msi/resolutions/spatial)
164
La Figure 106 ci-dessous reprend les caractéristiques principales de 3 images satellites
utilisées dans cet exercice.
S2B_MSIL1C_20180707T080609_N0206_
R078_T35NRA_20180707T120150
Nom du satellite LANDSAT-5 TM Sentinel-2
Lieu dans le système WRS WRS_PATH = 173 Tile T35NRA
Landsat (Path, Row) et WRS_ROW = 60
« Tile » de Sentinel-2
Système de coordonnées EPSG 32635 EPSG 32635
WGS 84 / UTM zone 35N WGS 84 / UTM zone 35N
Etendue spatiale 170 km north-south by 183 100 km * 100 km
approximative d’une km east-west
image/tile
Date d’acquisition 7 août 1987 12 juin 2018
7 juillet 2018
Résolution spatiale 30 m 10 m – 20 m – 60 m
Nombre de bandes 7 12
spectrales
Couverture nuageuse Partielle Partielle
Format original des images .tif .jp2
lors du téléchargement
Poids numérique du dossier 138 Mo / 460 Mo 782 Mo / 782 Mo
compressé/décompressé
contenant une image
Site web de référence pour https://earthexplorer.usgs.g https://scihub.copernicus.eu/d
le téléchargement ov/ hus/#/home
Page web reprenant les https://www.usgs.gov/land- https://sentinel.esa.int/web/sentinel/user-
resources/nli/landsat/landsat-5?qt- guides/sentinel-2-msi
caractéristiques techniques science_support_page_related_con=0#qt-
des images science_support_page_related_con
Figure 106 : Caractéristiques principales des 3 images utilisées dans cet exercice
165
8.5.4. Téléchargement des images avec prétraitements dans la foulée via
le plugin SCP dans QGIS
! Méthode proposée dans le cadre de ce manuel :
Landsat-5 : utilisez les indications reprises dans cette section pour télécharger et
prétraiter l’image Landsat-5 via le plugin SCP dans QGIS (relativement rapide).
Sentinel-2 : à ne pas faire pour les images Sentinel-2 car ces images prétraitées sont
mises à votre disposition dans le dossier de formation «
5_CLASSIFICATION_SUPERVISEE_ANALYSE_DIACHRONIQUE\DATA\IMAGES\SENTINE
L_2\ » en format « .tif » (confer la section 8.5.2.2 page 149) (méthode avancée
uniquement)
Ouvrez ou reprenez votre projet QGIS tel que préparé selon les indications de la section
« 8.3 Préparation du projet QGIS » page 129.
Menu QGIS: SCP > Télécharger de produits (3 onglets) (Figure 107 à Figure 110)
Figure 107 : Interface de l’outil « Menu SCP > Télécharger de Produits », onglet « Rechercher »
Ce menu permet de ne télécharger qu’une partie des bandes spectrales des images
satellites, ce qui est bien pratique afin de ne pas télécharger de bandes inutiles. Dans le
cadre de cet exercice, seules les bandes spectrales suivantes seront téléchargées :
166
Sentinel-2 : 2, 3, 4, 5, 8, 11 (pas les bandes 1, 6, 7, 8A, 9, 10, 12)
La procédure à suivre pour le téléchargement et le prétraitement des images via SCP dans
QGIS est reprise ci-dessous. Elle est à exécuter une fois pour l’image Landsat-5 et une fois
pour les 2 images Sentinel-2. Attention, il faut être connecté à internet !
8.5.4.1. Connexion aux sites web sources des images satellites via le plugin
SCP dans QGIS
167
Figure 108 : Onglet « Identifiant » du « Menu SCP > Télécharger de Produits > Identifiants »
168
ou moins important et de résulter en un bug de QGIS ou en une saturation de
votre disque dur, en fonction du nombre d’images concernées. Pour
supprimer de la liste les images que vous ne souhaitez pas télécharger,
sélectionnez-les dans la liste, puis cliquez sur le bouton .
o Spécificités aux images Sentinel-2
Le filtre « T35NRA » permet de ne rechercher que les images Sentinel-
2 correspondant à cette « tile »
Faites la recherche pour les 2 images Sentinel-2 (2 dates) l’une après
l’autre pour avoir en définitive les 2 images d’intérêt dans la liste des
produits
ATTENTION ! EN DATE DU 31 mars 2020, il n’est plus possible de
rechercher de cette manière des images Sentinel-2 remontant à avant
décembre 2018 (archives de plus de 1 an et 4 mois) (confer
https://scihub.copernicus.eu/userguide/LongTermArchive). Vous
pouvez cependant essayer. Peut-être que la situation a changé entre-
temps. Sinon préférez la solution proposée dans la section 8.5.2.4
page 158 avec téléchargement directement depuis le site web de
référence. Vous pourrez alors faire les prétraitements comme indiqué
dans la section 8.6 page 176. Une autre option est d’utiliser des
images Sentinel-2 plus récentes.
169
Figure 109 : Menu « Rechercher » du plugin SCP dans QGIS paramétré pour la recherche des
images Landsat-5 (haut) et Sentinel-2 (bas, avec l’ancienne interface de SCP) utiles à cet exercice
avec affichage des images résultant des recherches
8.5.4.3. Options de téléchargement des images via le plugin SCP dans QGIS
Menu SCP > Télécharger de Produits > Ligne « Télécharger » en bas de l’interface
(Figure 109)
o Décochez toujours « Seulement si les aperçus sont chargés en couche » car
cette option semble empêcher tout téléchargement,
o Cochez « Prétraitement des images » pour appliquer les prétraitements qui
seront sélectionnés dans la section suivante
170
o Cochez « Charger les bandes dans QGIS » pour charger automatiquement les
bandes dans QGIS après téléchargement et prétraitements
! NE PAS ENCORE LANCER LE TELECHARGEMENT ! Il faut d’abord paramétrer les
prétraitements comme indiqué dans la section suivante.
171
8.5.4.4. Prétraitement des images dans la foulée du téléchargement via le
plugin SCP dans QGIS
Figure 111 : Paramétrisation des prétraitements à appliquer aux images Landsat-5 et Sentinel-2
172
8.5.4.5. Téléchargement des images avec prétraitement dans la foulée via le
plugin SCP dans QGIS
Figure 112 : Interface de visualisation QGIS et panneau « Couches » contenant les 6 bandes
Landsat-5 prétraitées (préfixe « RT_ ») et réorganisées dans le « groupe » « LANDSAT5_19870807 »
Les données sont également enregistrées dans le répertoire choisi ci-dessus, sous la forme
de 2 dossiers par image (exemple à la Figure 113 pour Landsat-5):
173
Un dossier pour l’image brute et les données auxiliaires
o Landsat-5 : 120 Mo
o Sentinel-2 : ~500 Mo
Un dossier pour l’image prétraitée (dossier avec suffixe « _con » ou « _converted »)
o Landsat-5 : 289 Mo
o Sentinel-2 : 1.24 Go par image (bandes 2, 3, 4, 5, 8, 11, toutes les bandes
considérées ont été ré échantillonnées à 10m)
Figure 113 : Répertoire « LANDSAT_5 » dans l’explorateur Windows après téléchargement avec
prétraitement dans la foulée via le plugin SCP dans QGIS, contenant les bandes originales
décompressées (en haut) et prétraitées (en bas, dossier « converted » avec préfixe « « RT_ »)
Enfin, comme demandé, des « Jeu de bandes » ont été créés pour chaque image prétraitée
(Figure 114). Pour le vérifier :
Menu SCP > Jeu de bandes (« Band set 1, Band set 2, etc)
174
Figure 114 : Divers « Jeux de bandes » (« Band set ») visibles dans le plugin SCP > menu « Jeu de
bandes » après téléchargement et prétraitements des images avec création d’un jeu de bandes,
pour Landsat-5 (en haut) et Sentinel-2 (avec ici plus de bandes que recommandé dans ce manuel)
(en bas, avec une ancienne interface du plugin SCP).
L’utilité et la gestion d’un « Jeu de bandes » sont décrites dans la section 8.6.3 page 178.
175
8.6. Prétraitements des images satellites non couplés au
téléchargement via le plugin SCP dans QGIS
Les prétraitements ont pour objet de préparer les images avant la réalisation des
traitements thématiques (classification supervisée et analyse diachronique dans ce cas-ci).
! Méthode proposée dans le cadre de ce manuel : cette section ne doit pas être réalisée !
(confer la section 8.5.2.2 page 149)
Ni pour Landsat-5 (car la décompression a déjà été réalisé via la section 8.5.4 page
166 et suivantes),
Ni pour Sentinel-2 (car les images déjà décompressées et prétraitées sont mises à
disposition dans le dossier « …\DATA\IMAGES\SENTINEL_2\)
Les dossiers contenant les images Landsat-5 ou Sentinel-2 sont téléchargés à partir des sites
web de référence (confer section 8.5.2 page 143) sur le disque dur de votre ordinateur sous
format compressé.
176
8.6.2. Conversion de radiance vers réflectance / correction atmosphérique
DOS1 et création d’un « Jeu de bandes »
! Méthode proposée dans le cadre de ce manuel : cette section ne doit pas être réalisée !
(confer la section 8.5.2.2 page 149)
Ni pour Landsat-5 (car ce prétraitement a déjà été réalisé via la section 8.5.4 page
166 et suivantes),
Ni pour Sentinel-2 (car cela prend beaucoup de temps et les images prétraitées sont
mises à disposition dans le dossier « …\DATA\IMAGES\SENTINEL_2\)
Cette section présente la procédure à suivre dans le plugin SCP pour réaliser la conversion
des valeurs de pixels depuis des valeurs de radiance vers des valeurs de réflectance /
correction atmosphérique DOS1 et la création d’un « Jeu de bandes » sur les images
Landsat-5 et Sentinel-2 dans le cas où les images ont été directement téléchargées via les
sites web source présentés dans la section 8.5.2 et non via le plugin SCP.
177
Figure 115 : Menu « SCP > Pré-traitement > Sentinel-2 » paramétré pour un prétraitement à
appliquer sur toutes les bandes spectrales de Sentinel-2 (dans le cadre de cet exercice, il est
conseillé de ne pas travailler avec toutes les bandes) (cette illustration correspond à l’ancienne
interface de SCP).
8.6.3. Création d’un « Jeu de bandes » et affichage des images dans QGIS
! Méthode proposée dans le cadre de ce manuel (confer la section 8.5.2.2 page 149):
Landsat-5 : un « Jeu de bandes » a en principe déjà été créé pour l’image Landsat-5
via la section 8.5.4 page 166. Cependant si vous deviez refaire un jeu de bande, vous
trouverez les indications nécessaires ci-dessous.
Sentinel-2 : utilisez les indications de cette section pour créer un « Jeu de bandes »
pour chaque image Sentinel-2 fournies dans le dossier
« …\DATA\IMAGES\SENTINEL_2 ».
Visualisation des images d’un « Jeu de bandes »: une fois les « Jeux de bandes »
créés, affichez-les dans QGIS en suivant les indications ci-dessous. Par soucis de
rapidité, inutile de
o « Créer un raster virtuel de band »
o « Créer un raster assemblant les bandes »
Et optez plutôt pour la méthode de création d’un
o « Virtual band set »
Cette étape ne constitue pas à proprement parler un « prétraitement » mais est présentée
ici pour respecter la logique chronologique du chemin méthodologique proposé.
178
Une fois les images décompressées et éventuellement transformées en réflectance (sections
8.6.1 et 8.6.2) elles peuvent être affichées dans QGIS.
Cette section explique comment créer un « Jeu de bandes » et comment l’utiliser pour la
visualisation des images dans QGIS.
179
bandes 2, 3, 4, 5, 6, 7, 8, 8A, 11, 12. Sinon vous pouvez entrer ces
informations manuellement via un double-clic sur la case à éditer.
ATTENTION ! Cet outil n’ordonne pas toujours parfaitement les
bandes! Vérifiez toujours le résultat et corrigez l’ordre des bandes le
cas échéant avec les flèches ou le bouton qui sera utile dans
certains cas mais pas dans tous
Cette même interface (menu SCP > Jeu de bandes) permet également de créer un raster
virtuel et/ou un raster assemblant les bandes comme détaillé ci-dessous. Ces opérations ne
sont cependant pas nécessaires et ne seront pas commentées dans le cadre de ce manuel.
Si vous avez créé un « jeu de bandes » comme indiqué à la section précédente ou qu’un
« jeu de bandes » est disponible dans le plugin SCP, vous pourrez l’utiliser pour créer un
« Jeu de bandes virtuel » ou « Virtual band set » via le plugin SCP. Ce « Virtual band set »
permettra d’alterner facilement/rapidement des compositions colorées d’un même raster et
de réaliser des améliorations de contraste. Pour ce faire :
Menu SCP > Jeu de bandes : sélectionner l’onglet (menu « détails du jeu de bandes »
ou « Band set definition ») du « Jeu de bandes » pour lequel créer un « Virtual band
set »
Dans l’interface principale de QGIS > Barre d’outils SCP > boutons
(ou via le menu SCP > Outils basic > Liste RGB) > inscrire
180
manuellement (ou sélectionner) la/les compositions colorées souhaitées (ex : 3-2-1
ou 5-3-2) dans le menu déroulant
ATTENTION ! Les numéros font référence à l’ordre des bandes dans la table
du « Jeu de bandes » actif (numéros dans la marge gauche de cette table) ET
PAS AUX NUMEROS OFFICIELS DES BANDES SPECTRALES ! En conséquence, si
vous voulez visualiser une composition « vraies couleurs » ou « couleurs
naturelles » (confer section 8.7.4.1 page 213), vous devrez indiquer dans le
menu déroulant présenté ci-avant les numéros de lignes du jeu de bandes qui
correspondent aux bandes spectrales des longueurs d’onde Rouge, Verte et
Bleue, dans cet ordre.
Après un petit délai, le « Virtual band set X » vient s’ajouter dans le panneau
« Couches » de QGIS (« X » fait référence au numéro du « Jeu de bandes »
dans le menu « Jeu de bandes »).
Attention, si vous avez plusieurs « Virtual band set X » la modification de la
composition colorée via le menu déroulant ne s’applique QUE sur le jeu de
bandes actif/sélectionné dans le menu « Jeu de bandes » !
permettront d’appliquer des améliorations de contrastes qui seront
calculées sur la partie de l’image actuellement visible dans l’interface de
visualisation spatiale de QGIS.
181
8.6.4. Découpage spatial de l’image sur la zone d’intérêt (en anglais : «
Spatial subset », « clip », « crop »)
! Méthode proposée dans le cadre de ce manuel : (confer la section 8.5.2.2 page 149) :
Landsat-5: effectuez un découpage spatial sur l’image Landsat-5 à l’aide de cette
section
Sentinel-2 : à ne pas faire car les images mises à disposition dans le dossier
« …\DATA\IMAGES\SENTINEL_2\... » sont déjà découpées spatialement (car sinon
elles sont trop volumineuses)
Les images satellites téléchargées couvrent parfois une zone plus grande que la zone
d’intérêt. Dans ce cas, il convient, afin de réduire la taille des données manipulées et par là
d’accélérer la vitesse des traitements et de réduire la place que prendront les données
générées sur votre disque dur, de découper spatialement l’image pour ne conserver que la
zone d’intérêt.
Définition d’un jeu de bandes à découper (si le jeu de bandes existe déjà, pas
besoin de le recréer):
o Confer section 8.6.3.1 ci-dessus
Découpage spatial d’un « Jeu de bandes » créé préalablement (Figure 116):
o Menu : SCP > Prétraitement > Découper plusieurs rasters
Sélectionnez le numéro du « Jeu de bandes » à découper (le numéro
du jeu est disponible via le menu « Jeu de bandes »)
Préfixe de sortie : « clip » : le préfixe qui sera utilisé pour le nom des
bandes découpées.
Pour définir la zone à découper:
Utilisez, dans le cadre de cet exercice, le shapefile
« ZONE_ETUDE_KASESE.shp ». Attention le shapefile doit être
dans le même système de coordonnées que les bandes raster à
découper !
Une autre possibilité serait de délimiter la zone de manière
interactive en cliquant sur puis, en cliquant, dans la fenêtre
de visualisation spatiale de QGIS, avec un clic-gauche sur le
coin supérieur gauche de la zone à délimiter, et avec un clic-
droit sur le coin inférieur droit de la zone à délimiter. Ou entrer
les coordonnées manuellement.
182
~ 10 secondes pour 6 bandes Landsat-5 (37 Mo)
~ 2min30 pour 6 bandes 2-3-4-5-8-11 Sentinel-2 en résolution
originale (10 et 20 m) (430 Mo)
(~ 4min30 pour 10 bandes Sentinel-2 à 10 m (650 Mo))
o Les bandes converties sont automatiquement
enregistrées dans le répertoire de sortie
chargées dans QGIS (nom de fichier commençant par le préfixe
« clip_ »).
Attention ! Veillez à organiser correctement vos données dans
le panneau « Couches » de QGIS, éventuellement à l’aide du
bouton « Ajouter un groupe »
Visualisation du « Jeu de bandes » découpé spatialement
o Confer les indications de la section 8.6.3 page 178.
Figure 116 : Menu de SCP permettant de réaliser un découpage spatial sur un « Jeu de bandes »
Lorsque les différentes bandes spectrales d’une image satellite sont disponibles séparément,
sous la forme de fichiers indépendants il peut être pratique, avant de faire d’autres
traitements, de combiner les différentes bandes spectrales en une seule image unique multi-
bandes qui sera de la sorte plus facilement manipulable car l’ensemble des traitements
pourront s’appliquer simultanément sur l’ensemble des bandes spectrales regroupées.
183
Cette combinaison de bandes spectrales peut se faire via une opération
d’assemblage/regroupement/collage des bandes (anglais : « Layer stacking »).
Fonction SCP : SCP > Pré-traitement > Assembler les bandes raster
Fonction QGIS : Raster > Divers > Fusionner… > et cocher « Placer chaque fichier en
entrée dans une bande séparée » ;
Si l’outil le permet, classez les bandes spectrales dans le bon ordre avant l’assemblage et
choisissez la résolution spatiale de l’image résultante. Il est conseillé, dans le cas des images
LANDSAT, d’utiliser 30 m de résolution spatiale, même si vous considérez dans l’assemblage
une bande de meilleure résolution spatiale (exemple : la bande 8 des images LANDSAT8 qui
est à 15m). Attribuer une plus haute résolution spatiale ne fera qu’alourdir fortement votre
image et sera probablement inutile. Vous pouvez aussi exclure la bande à 15 m.
8.6.6. Pan-sharpening
Ce prétraitement ne sera pas utilisé dans le cadre de cet exercice. Il est mentionné à titre
d’information.
La technique de pan-sharpening est une méthode qui permet de fusionner une bande
spectrale de plus haute résolution spatiale (par exemple, la bande panchromatique des
images LANDSAT à 15m) avec une ou plusieurs bandes spectrales de moins haute résolution
spatiale (par exemple, les autres bandes des images LANDSAT à 30m) afin d’obtenir une
image qui présente à la fois la meilleure résolution spatiale (15m) et les qualités spectrales
des bandes à 30m. L’image résultante facilite la photo-interprétation de l’image.
Menu QGIS : Traitement > Boite à outils > GDAL > Divers raster > Pansharpening
184
Le rééchantillonnage consiste notamment au changement de la taille des pixels d’un raster
pour par exemple en diminuer la résolution spatiale ou pour harmoniser la taille des pixels
des différentes bandes spectrales constituant un raster unique multi-bandes.
Dans le cadre de cet exercice, il s’est avéré que les images Sentinel-2 à 10m de résolution
spatiale sont assez lourdes à traiter et que certains traitements, en particulier la
classification, prennent beaucoup trop de temps à cette résolution spatiale.
La fonction de rééchantillonnage n’est pas directement disponible dans SCP. L’on pourra
donc utiliser une fonction de rééchantillonnage du menu principale de QGIS.
Pour rééchantillonner toutes les bandes d’une image Sentinel-2 à 30 m dans QGIS (Figure
117):
185
Figure 117 : Rééchantillonnage de 6 bandes d’une image Sentinel-2 à 30 mètres (-tr 30 30) via le
traitement par lot de la fonction « Convertir… » de QGIS
Lorsque la zone d’étude est couverte par plusieurs images voisines, il peut être intéressant
de les fusionner/mosaïquer afin d’obtenir une seule grande image couvrant toute la zone et
plus facilement manipulable.
Fonction QGIS : Raster > Divers > Fusionner… >
186
8.6.9. Mosaïquage de 2 images Sentinel-2 superposées pour
réduire/supprimer la couverture nuageuse
! Méthode proposée dans le cadre de ce manuel : (confer la section 8.5.2.2 page 149) :
Landsat-5: cette image n’est pas concernée par cette section
Sentinel-2 : suivez les indications de cette section pour mosaïquer les 2 images
Sentinel-2
Dans le cas de cet exercice, les 2 images Sentinel-2 les moins ennuagées identifiées pour
l’étude présentent tout de même une couverture nuageuse relativement importante et
gênante pour l’analyse diachronique. Il convient dès lors de les mosaïquer afin de produire
une image présentant une couverture nuageuse réduite.
Le mosaïquage permet de remplacer les parties d’une image sur lesquelles un masque (de
nuages par exemple) a été appliqué, par les pixels d’une autre image (idéalement, sans
nuage aux endroits masqués) pour former une image avec une couverture nuageuse
amoindrie (Figure 118).
187
Figure 118 : Principe de mosaïquage d’une image masquée (image A) avec une autre image (image
B)
188
Sentinel-2 20180612 Sentinel-2 20180707
Figure 119 : Composition vraies couleurs des 2 images Sentinel-2 (au-dessus) avec affichage de la
couche vectorielle de masque de nuages « MSK_CLOUDS_B00.gml » fournie avec les images avec
les nuages denses (opaques) en bleu et les cirrus (voile léger) en vert (en-dessous).
Il est clairement visible que ce désennuagement ne sera pas total (partie nord-ouest de
l’image). Il peut alors être utile d’utiliser des masques de nuages améliorés (confer plus loin)!
https://fromgistors.blogspot.com/p/user-manual.html?spref=scp
189
devra très vite fastidieuse avec l’augmentation du nombre et la complexité de la
forme des nuages à masquer
Une combinaison de ces techniques
8.6.9.3.1. Création d’un masque de nuages à partir d’un fichier vectoriel de masque
existant
Les images Sentinel-2 sont toujours accompagnées d’une bande d’information appelée
« MSK_CLOUDS_B00.gml » qui est un fichier vectoriel identifiant très grossièrement les
endroits ennuagés sur l’image (Figure 119). Ce fichier est décrit ici
https://earth.esa.int/web/sentinel/technical-guides/sentinel-2-msi/level-1c/cloud-masks.
Nous pouvons y lire que le fichier vectoriel peut prendre 2 valeurs : dense (opaque) ou cirrus
(léger voile) (Figure 119 et Figure 120), valeurs reprises dans la table d’attributs de la couche
vectorielle.
Figure 120 : Aperçu de la table d’attributs du fichier vectoriel de masque de nuage d’une image
Sentinel-2 avec présence des 2 types de nuages : CIRRUS et OPAQUE
Remarquez à la Figure 119 que pour l’image Sentinel-2 20180707 ce masque est
particulièrement imprécis ! Cela constituerait une faiblesse/limitation de l’approche si ce
masque était utilisé. Dans une situation réelle, il est indiqué d’identifier la couverture
nuageuse de manière plus précise. Dans le cadre de cet exercice le masque de nuages
réellement utilisé sera celui créé selon la méthode de la section 8.6.9.3.2 page 192. La suite
de cette section-ci est donnée pour info.
Pour créer un masque de nuages pour l’image 20180612 à partir de la couche vectorielle de
masque fournie avec l’image:
190
l’image. Remarquez que sur la zone d’étude, il n’y a que des nuages
« opaques/denses » et pas de « cirrus ».
L’exporter au format shapefile : clic-droit sur le fichier « MSK_CLOUDS_B00
MaskFeature » dans le panneau Couches de QGIS > Exporter > Sauvegarder les
entités sous… > Choisir le format shapefile, le répertoire de sortie et le nommer par
exemple « MASQUE_NUAGE_20180612.shp » et l’ajouter dans le projet QGIS.
Menu SCP > Pré-traitement > Vecteur vers raster >
o (attention éviter répertoire long/ avec espace)
o Sélectionner un vecteur: cliquez sur pour rafraichir la liste de couches
disponibles et sélectionnez le masque créé à l’étape précédente «
MASQUE_NUAGE_20180612.shp »
o Cochez « Utilisez une valeur constante » et indiquez « 1 » pour que le raster
de masque produit prenne une valeur de 1 pour les nuages
o Sélectionnez le raster de référence : cliquez sur , puis sélectionnez une
bande de l’image Sentinel-2 découpée spatialement correspondant au
masque utilisé, par exemple
« clip_RT_L1C_T35NRA_A015519_20180612T082541_B02.tif » pour créer un
raster de la même taille et aligné sur le raster Sentinel-2.
Le raster de masque de nuages est produit et prend une valeur de « 1 » pour les nuages et
« NoData » pour les autres zones (consultable via l’outil de QGIS).
191
8.6.9.3.2. Création d’un masque de nuages par classification des nuages dans l’image à
masquer
Un masque de nuage peut être réalisé par classification des nuages présents sur l’image à
masquer.
Cette étape de classification est décrite à la section 8.8.8 page 240 vu qu’il s’agit
essentiellement d’un processus de classification puis raffinement de la classification.
Dans le cadre du TP, vous utiliserez simplement la couche de masque résultant de cette
étape qui a déjà été faite pour vous et qui est disponible dans le répertoire
« …\DATA\IMAGES\SENTINEL_2\SENTINEL_2_MASQUE_NUAGES_20180612.tif ».
Une méthode de seuil aurait également pu être utilisée. Elle consiste à simplement
reclassifier les pixels d’une bande en-dessous/au-dessus d’une valeur seuil vers une
constante. Dans ce cas-ci, on pourrait par exemple reclassifier la bande bleu, celle qui
présente le plus de contraste entre sol et nuages, en identifiant la meilleure valeur seuil
permettant de faire cette distinction entre sol et nuage. Ceci peut se faire facilement via la
calculatrice raster :
192
8.6.9.4. Appliquer le masque de nuages sur une image
! Méthode proposée dans le cadre de ce manuel : (confer la section 8.5.2.2 page 149):
Landsat-5: pas concernée par cette section
Sentinel-2 : image 20180612 seulement
Pour appliquer un masque de nuages, existant sous la forme d’un raster uni-bande dont
la/les valeurs de pixels correspondant aux zones à masquer sont clairement définies, sur
toutes les bandes de l’image 20180612 présentes dans le Jeu de bandes de cette image:
Figure 121 : Application d’un masque de nuages sur une image dans SCP. Remarque : le nom du
fichier correspondant à ce masque et fourni dans le dossier de cet exercice est en fait
« SENTINEL_2_MASQUE_NUAGES_20180612.tif » (différent de ce qui est indiqué dans cette
figure). (Ancienne interface de SCP)
193
8.6.9.5. Mosaïquer l’image masquée et la deuxième image
Pour remplacer les pixels « NoData » de l’image 20180612 masquée par les pixels de l’image
20180707 :
Menu SCP > Jeu de bandes > Ajouter un nouveau Jeu de bandes via le bouton >
Y ajouter les bandes masquées « masque_clip_RT_... 20180612…»
Menu SCP > Jeu de bandes > Assurez-vous d’avoir un Jeu de bandes qui correspond à
l’image qui sera utilisée pour remplacer les pixels masqués de l’image 20180612, soit
un jeu de bandes pour l’image 20180707 découpée spatialement, sinon créez-le
Menu SCP > Prétraitement > Mosaïque de Jeu de bandes
o Entrez le numéro du Jeu de bandes contenant les bandes masquées
20180612, suivi d’une virgule, suivi par le numéro du Jeu de bandes
contenant les bandes de remplacement (20180707)
o Préfixe « mosaique » (possibilité d’utiliser plusieurs rasters).
194
Image Sentinel-2 20180612
A
Image Sentinel-2 20180707
C
Image Sentinel-2 20180612 après
application du masque de nuages amélioré
(masque obtenu par classification des
nuages et raffinement) et remplacement des
pixels de ce masque par les pixels de l’image
Sentinel-2 20180707 via la technique du
mosaïquage
D
Figure 123 : Comparaison des 2 images Sentinel-2 initiales avec les images résultant du
mosaïquage. Les images mosaïquées présentent une couverture nuageuse amoindrie.
195
o Le masque de nuage créé à partir du masque vectoriel fourni avec l’image 20180612
ne masque pas tous les nuages, ce qui a pour conséquence qu’une partie d’entre eux
se retrouvent dans l’image résultant du mosaïquage (C).
o Ceci constitue une grosse limitation de la méthode car ces pixels ne seront pas
masqués et produiront un biais dans la matrice de transition !
Un néocanal est une nouvelle couche d’information (une nouvelle bande raster) calculée à
partir des bandes spectrales disponibles dans une image et qui permet, idéalement, de faire
ressortir certaines informations non directement ou moins facilement accessibles à partir
des bandes spectrales de base.
196
8.6.11. Géoréférencement
Cette opération ne sera pas utilisée dans le cadre de cet exercice, car pas nécessaire. Elle est
mentionnée à titre d’information.
Lors d’une analyse diachronique consistant en la comparaison, pixel par pixel, des
changements d’occupation du sol s’étant produit entre 2 images acquises à 2 dates
différentes, il est évidemment très important que les images à comparer soient parfaitement
géoréférencées l’une par rapport à l’autre, c’est-à-dire que les pixels de l’image initiale
superposant les pixels de l’image finale correspondent à une même zone, pixel par pixel.
Il est nécessaire de vérifier si le géoréférencement relatif des 2 images est suffisant ou pas.
Ceci peut simplement se faire via analyse visuelle dans QGIS via :
Il se peut que le géoréférencement initial des images téléchargées soit suffisant que pour
réaliser correctement une analyse diachronique.
Géoréférencement avec points de calage (en anglais : Ground Control Point (GCP))
Géoréférencement « image-to-image » : géoréférencement d’une image
relativement à une autre.
197
o https://docs.qgis.org/3.4/fr/docs/training_manual/forestry/map_georeferenc
ing.html?highlight=g%C3%A9or%C3%A9f%C3%A9rencement
des nuages, et parfois leurs ombres projetées au sol, sont présents sur l’image
une partie de l’image ne doit pas être analysée (par exemple : à l’extérieur des
limites d’un parc naturel pour lequel une analyse spécifique au parc est souhaitée)
Lors d’une analyse diachronique de l’évolution de l’occupation du sol, afin que les 2 images
classifiées restent comparables, un même masque devra être appliqué sur les 2 images. Ce
masque correspondra donc à la somme des masques créés pour chacune des deux images.
Les étapes à suivre seront donc :
198
Figure 124 : Principe de la création d’un masque de nuages multi-dates
199
https://www.usgs.gov/land-resources/nli/landsat/landsat-collection-1-level-1-
quality-assessment-band?qt-science_support_page_related_con=0#qt-
science_support_page_related_con
“The Collection 1 Level-1 QA band can be stretched to emphasize the light and dark pixels
for an immediate survey of a scene's general quality conditions.”
Section “Landsat 4-5 TM, Landsat 7 ETM+ Level-1 Possible Attributes, Pixel Values, and Pixel
Value Interpretations”, 2ème table:
On peut notamment résumer très grossièrement une partie de l’information contenue dans
cette table comme ceci:
[704-764] = cloud
[928-972] = cloud shadow
[1696-1740] = snow
Figure 125 : Visualisation dans QGIS de la bande BQA de l’image Landsat-5 renseignant sur la
couverture nuageuse
200
Une analyse attentive devrait être opérée afin de déterminer les pixels à conserver / à ne
pas conserver.
Dans ce cas-ci nous allons considérer que les pixels de la bande BQA ayant une valeur entre
[704 - 972] sont impactés par les nuages et sont donc à masquer.
Pour créer le masque de nuage correspondant, il suffit de reclassifier les valeurs du raster
avec:
Utilisez la fonction « r.reclass » qui dépend de GRASS et disponible dans Traitement >
Boite à outils ». Attention ! Il faut donc que votre projet QGIS ait été ouvert avec «
QGIS avec GRASS » pour que cette fonction soit disponible.
o Documentation ici https://grass.osgeo.org/grass70/manuals/r.reclass.html
Paramétrage comme à la Figure 126.
Règles de reclassification à indiquer :
o 704 thru 972 = 1
o *= 2
Nommez le masque raster reclassifié « MASQUE_NUAGE_LANDSAT5_BQA.tif » par
exemple
Figure 126 : Fonction r.reclass de GRASS dans QGIS pour créer le masque de nuages de l’image
Landsat-5 à partir de son fichier BQA
201
8.6.12.3. Création d’un masque de nuages pour l’image Sentinel-2 mosaïquée
Le masque de nuage à considérer pour l’image Sentinel-2 mosaïquée correspond à
l’intersection des masques de nuages améliorés des 2 images mosaïquées (20180612 et
20180707) eux-mêmes créés selon la méthode présentée à la section 8.8.8 page 240.
Figure 127 : Principe de construction d’un masque de nuages commun à 2 masques de nuages (noir
= NoData, blanc = valeur de 1)
202
Ouvrir dans QGIS les 2 rasters correspondant aux masques de nuages améliorés des 2
images Sentinel-2, dont les valeurs de pixels sont de 1 pour les nuages et 2 ailleurs.
Ces 2 masques ont déjà été produits pour vous et sont disponibles dans le répertoire
« …\ DATA\IMAGES\SENTINEL_2\ » sous les noms de :
o SENTINEL_2_MASQUE_NUAGES_20180612.tif (repris sous le nom
« MASQUE_NUAGE_20180612.tif » dans l’illustration ci-dessous)
o SENTINEL_2_MASQUE_NUAGES_20180707.tif (repris sous le nom
« MASQUE_NUAGE_20180707.tif dans l’illustration ci-dessous)
Utiliser la calculatrice raster de QGIS pour faire le produit des 2 masques :
o Raster > Raster Calculator… >
Appelez le raster résultant « MASQUE_NUAGE_2018_0612_0707_PRODUIT.tif » par
exemple
Les 4 valeurs possibles sont :
o 1 * 1 = 1 Nuage * Nuage = Nuage commun aux 2 images
o 1 * 2 = 2 Nuage * Autre = Autre
o 2 * 1 = 2 Autre * Nuage = Autre
o 2 * 2 = 4 Autre * Autre = Autre
Le raster résultant prendra donc une valeur de 1 dans ce cas-ci pour les nuages commun aux
2 images Sentinel-2, et 2 ou 4 dans les autres situations.
Figure 128 : Produit des 2 masques de nuages améliorés des 2 images Sentinel-2, dans la
calculatrice raster de QGIS, pour produire le masque de nuages commun aux 2 images Sentinel-2
(remarque : le nom des masques dans cette illustration ne correspond pas aux noms des masques
mis à disposition dans le dossier de cet exercice comme indiqué plus haut)
203
Pour faciliter la suite des opérations, ce raster peut être reclassifié selon les règles
reprises dans la figure ci-dessous afin de garder une valeur de 1 pour les nuages
communs aux 2 images et une valeur de 2 ailleurs
o Traitement > Boite à outils > r.reclass (dépend de GRASS)
o Raster en sortie =
« MASQUE_NUAGE_2018_0612_0707_PRODUIT_reclass.tif »
Figure 129 : Utilisation de la fonction GRASS « r.reclass » pour reclassifier le raster de masque de
nuages commun aux 2 images Sentinel-2
Réaliser une somme des 2 masques (valeurs : 1 = nuage ; 2 = autre) produits aux
sections précédente via la calculatrice raster de QGIS (Raster > Raster Calculator…)
o « MASQUE_NUAGE_2018_0612_0707_PRODUIT_reclass.tif »
o « MASQUE_NUAGE_LANDSAT5_BQA_CLIP.tif »
Le raster résultant prend les valeurs suivantes :
o 2 : nuages communs aux 2 images
o 3 : nuages sur une seule des 2 images
o 4 : absence de nuage sur les 2 images
Nom du raster résultant :
o « MASQUE_NUAGE_SOMME_2_L5_S2.tif »
204
8.6.12.5. Application du « masque fusionné » ou « masque somme » à
chacune des 2 images
Le « masque somme » doit ensuite être appliqué sur chacune des 2 images à considérer pour
l’analyse diachronique, pour les rendre comparables, soit les images :
Landsat-5 : clip_RT_LT05_L1TP_173060_19870807_20170211_01_T1_1987-08-07_...
Sentinel-2 : mosaique_masque_30m_clip_RT_T35NRA_20180612T080611_...
Menu SCP > Jeu de bandes
o Assurez-vous d’avoir un Jeu de bandes qui correspond à chacune des 2 images
sur lesquelles appliquer le « masque somme », sinon créez-les
Menu SCP > Prétraitement > Masque nuage
o Jeu de bandes : mettre le numéro du Jeu de bandes sur lequel appliquer le
masque
o Sélectionner la classification: sélectionner le «
MASQUE_NUAGE_SOMME_2_L5_S2.tif »
o Valeurs de classe du masque: Attention ! Mettre 2,3 (confer section
précédente)
o Décochez buffer
o Préfixe : masque_som
o Cliquez sur , sélectionnez le répertoire
Les images masquées avec le masque somme prenne une valeur de NoData aux endroits
masqués.
Figure 130 : Application du « masque somme » sur un jeu de bandes (ancienne interface de SCP)
205
Landsat-5 avec masque somme Sentinel-2 mosaïquée avec masque somme
206
8.7. Observation des images et identification des classes
d’occupation du sol
Les sections suivantes discutent de ces aspects et présentent diverses méthodes facilitant
l’analyse visuelle de l’image.
Les classes objectifs ou MACRO-CLASSES : ce sont les classes d’intérêt que vous
désirez analyser et/ou faire apparaitre dans la représentation cartographique finale
et qui seront donc reprises dans la légende de la carte d’occupation du sol produite.
Ces macro-classes peuvent englober plusieurs sous-classes spectrales qui
correspondent à des signatures spectrales différentes. Ce sont ces macro-classes qui
feront l’objet de la validation par matrice de confusion et qui seront étudiées dans la
matrice de transition.
Les sous-classes ou CLASSES SPECTRALES : ce sont les classes d’occupation du sol qui
présentent des signatures spectrales « différentes ». Ce sont les classes spectrales
qui seront utilisées pour le processus de classification par les algorithmes de
classification se basant sur la valeur des pixels. Il s’agit donc de définir autant de
classes spectrales qu’il existe de classes d’occupation du sol présentant des
signatures spectrales bien différentes. Ces « sous-classes spectrales » pourront être
fusionnées en « macro-classes » après la classification pour une meilleure
représentation cartographique (confer ci-dessous la section 8.10 page 247).
207
Figure 131 : Notions de « macro-classe » et « classes spectrale ». Source :
https://semiautomaticclassificationmanual.readthedocs.io/en/latest/remote_sensing.html#
classes-definition
Lors d’une analyse diachronique, les classes d’occupation du sol rencontrées dans les
différentes images à comparer (images acquises à différentes dates), peuvent différer en
partie (voire totalement dans une situation extrême). Par exemple, la 1ère image pourrait
présenter une zone forestière unique qui aurait été transformée totalement en zone agricole
et bâtit dans la 2ème image.
208
Lors de l’identification/sélection des classes spectrales à analyser, veillez à ne pas être trop
ambitieux et à ne pas surestimer la capacité de classification/distinction précise des
différentes classes d’occupation du sol par la méthode de classification que vous allez
utiliser… et donc à ne pas vouloir distinguer l’indistinguable ! Une classification d’image a
ses limites, qui sont parfois (souvent) bien vite atteintes… Il vous faudra donc trouver le juste
compromis entre le nombre de classes à identifier et la précision de classification
atteignable pour chacune de ces classes.
En outre, si, pour une raison ou l’autre, vous ne pouvez pas attribuer avec certitude une
classe d’occupation du sol à certains pixels, cela constituera une limitation à la
méthodologie de classification. Une classification de qualité est parfois impossible et il faut
l’accepter!
En particulier, lors de cette étape d’identification des classes d’occupation du sol, en plus
d’identifier les « macro-classes » et les « classes spectrales », veillez à créer des classes
spectrales pour :
Dans le cas de cet exercice, les figures ci-dessous reprennent les macro-classes et sous-
classes spectrales proposées pour chacune des 2 types d’images (Landsat-5 et Sentinel-2).
Notez que par soucis de concision, de simplicité et de rapidité ce nombre de classes a été
volontairement réduit.
209
Cette description des classes utilisée peut se faire via un tableau comme celui repris à la
Figure 133 (bien que les définitions soient ici assez sommaires). Remarquez cependant que
les couleurs des différentes classes apparaissant à l’écran peuvent fortement changer en
fonction de l’amélioration de contraste appliquée. Idéalement, il faudrait donc renseigner les
paramètres de l’amélioration de contraste utilisée, ce qui ne sera pas fait ici.
Couvert
Couvert Couvert végétal
végétal
végétal faible faible, vert clair
faible
Eau
Roche, présente
Roche Roche
en montagne
Sol nu
210
Sol nu brûlé
Sol nu brûlé
mauve
mauve foncé
foncé
Figure 133 : Exemple de description des « macro-classes » et « classes spectrales » des classes
d’occupation du sol pouvant être utilisées lors de la classification supervisée de l’image Landsat-5
utilisée dans le cadre de cet exercice.
Une autre manière de présenter les classes d’occupation du sol étudiée, peut-être plus
représentative et informative, est de présenter quelques extraits/morceaux de l’image
satellite à classifier, en composition vraies couleurs, à grande échelle (zoom in) et annotée
de numéros renvoyant à une légende identifiant les classes (Figure 134).
211
Figure 134 : Exemple de présentation des macro-classes d’occupation du sol utilisées lors d’une
classification supervisée d’une partie de l’île de Negros aux philippines. Source :
https://orbi.uliege.be/handle/2268/33068, page 61.
212
8.7.4. Observation des images en 2D
Voici ci-dessous une liste non exhaustive de techniques de visualisation d’image permettant
d’en améliorer l’analyse visuelle.
213
Vraies couleurs RGB = B3-B2-B1 Fausses couleurs RGB = B4-B7-B3
Figure 135 : Image Landsat-5 TM en composition « vraies couleurs » (à gauche) et en composition
« fausses couleurs » (un exemple parmi les nombreuses possibles) (à droite), toutes deux avec
amélioration de contraste.
Figure 136 : Exemple de visualisation d’une zone forestière présentant 2 types de forêts, forêt
feuillue au Nord-Ouest, et forêt résineuse (conifères) au Sud-Est, pendant la saison de végétation,
en composition vraies couleurs (gauche) et fausses couleurs (Infrarouge-Rouge-Vert) (droite).
Source : Image sentinel-2 d’une région au Sud de la ville de Namur en Belgique en juillet 2019
affichée dans le viewer « TERRASCOPE » https://viewer.terrascope.be/terrascope/.
214
8.7.4.2. Amélioration de contraste
L’amélioration de contraste consiste en une transformation virtuelle des données afin
d’utiliser toute la gamme d’intensité de couleur disponible pour visualiser les données à
l’écran. Cette amélioration ne modifie donc pas vos données sources (la valeur des pixels).
Elle ne fait qu’attribuer différemment la palette de couleurs disponibles à votre image afin
de faire ressortir certains éléments plus clairement.
215
8.7.5. Observation multi-temporelle
L’analyse de la dynamique/variation saisonnière (ou simplement temporelle) des classes
d’occupation du sol à l’aide de plusieurs images acquises à différentes période de l’année
(ou à différents moments) peut également aider à l’identification précise de certaines classes
d’occupation du sol, comme illustré à la Figure 138 qui montre :
Eté, forêt à feuilles persistantes verte et forêt Hiver, forêt à feuilles persistantes verte et forêt
décidue verte décidue brune
Figure 138 : Exemple de visualisation en composition vraies couleurs d’une zone forestière
présentant 2 types de forêts, forêt feuillue décidue (feuilles caduques) au Nord-Ouest, et forêt
résineuse à feuilles persistantes (conifères) au Sud-Est, pendant la saison de végétation en été (à
gauche) et en hiver (à droite). Source : Images sentinel-2 d’une région au Sud de la ville de Namur
en Belgique affichées dans le viewer « TERRASCOPE » https://viewer.terrascope.be/terrascope/.
216
8.7.6. Observation des images en 3D via une représentation du relief
naturel
Une représentation 3D du relief naturel du paysage à classifier (Figure 140) peut aider
l’opérateur à mieux comprendre l’organisation de ce paysage et par là l’aider dans sa tâche
de photo-interprétation ayant pour but final l’identification des différentes classes
d’occupation du sol.
Pour visualiser une couche d’information sous la forme d’une représentation 3D du relief
naturel, il est nécessaire d’avoir à disposition une couche d’information correspondant à ce
relief naturel, souvent appelée Modèle Numérique de Terrain (MNT) (en anglais : Digital
Elevation Model (DEM)).
Dans le cadre de cet exercice le MNT utilisé est celui présenté à la section 8.4.2.1 page 135
ci-dessus.
217
Figure 139 : Interface de construction d’une vue 3D avec la « Vue Cartographique 3D » de QGIS
L’extension (à installer) QGIS2threejs, qui après installation (via le menu Extension >
Installer/Gérer les Extensions), est disponible dans le menu « Internet > QGIS2threejs
> QGIS2threejs Exporter ». Cette extension, en plus de permettre une représentation
3D, permet également d’exporter cette représentation sous la forme de webmap
(carte au format web, interactive (zoomable et orientable en 3D))
218
Figure 140 : Représentation 3D en vraies couleurs de la zone d’étude pour l’image Landsat-5 (au-
dessus) et Sentinel-2 mosaïquée (en-dessous) produites avec la Vue Cartographique 3D de QGIS
219
8.7.7. Utilisation de données complémentaires pour la photo-
interprétation d’une image
La photo-interprétation (reconnaissance des différentes occupations du sol par analyse
visuelle) d’une image est parfois difficile :
Soit parce que l’image ne présente pas une qualité suffisante (résolution spatiale et
spectrale adéquate),
Soit à cause d’un manque d’expérience du photo-interpréteur ou d’un manque de
connaissance de la zone.
Images satellites de « Google Earth Pro » (exemple sur la zone d’étude à la Figure 77
page 127). « Google Earth Pro » dispose d’images satellites à très haute résolution
spatiale pour une partie sans cesse grandissante de la planète et permet d’accéder
assez facilement à des images satellites d’archive en composition vraies couleurs via
l’outil « Images d’archive » (confer section 3 page 31 pour une présentation détaillée
de « Google Earth Pro »). La très haute résolution spatiale permettra parfois de
mettre un nom précis sur des classes d’occupation du sol qui, à plus faible résolution
spatiale (exemple : sur une image Landsat à 30 m ou une image Sentinel-2 à 10 m), se
distinguent spectralement (couleurs différentes) mais sont difficilement
reconnaissables (Figure 141 ci-dessous).
Cartes topographiques
Une enquête de terrain avec relevé GPS
Divers documents donnant des informations sur l’occupation du sol de la zone.
220
Très haute résolution spatiale : ~20 cm Haute résolution spatiale : 10 m
Forêts feuillue et résineuse reconnaissables Forêts feuillue et résineuse non reconnaissables
Figure 141 : Exemple de visualisation d’une zone forestière présentant 2 types de forêts, forêt
feuillue au Sud-Est (cimes rondes), et forêt résineuse (conifères) au Nord-Ouest (cimes pointues),
pendant la saison de végétation, en composition vraies couleurs, en très haute résolutions spatiale
(~20 cm) (à gauche, dans Google Earth Pro) et à 10 m de résolution spatiale (à droite, Sentinel-2
dans le viewer TERRASCOPE) (Sud de la ville de Namur en Belgique).
221
8.8. Classification supervisée par maximum de
vraisemblance
222
Figure 142 : image aérienne d’un verger (arbres = boules vertes, sol = pixels gris) : image originale
(à gauche) et image après segmentation (à droite).
Figure 143 : Principe de l’arbre de décision pour classifier une image satellite multispectrale
Notez également qu’une classification d’image peut utiliser des indices dérivés des bandes
spectrales de l’image à classifier (NDVI, fAPAR, etc) et faire intervenir une série de données
auxiliaires autres que l’image à classifier comme l’altitude, l’exposition, la pente et toute
autre information spatiale (Figure 143).
223
8.8.2. La classification supervisée : un processus itératif
La mise au point d’une méthode de classification supervisée pour un cas particulier
nécessitera en général de réaliser une série de tests et de corriger itérativement la méthode
en fonction de la qualité des résultats obtenus. Cette méthode par essais-erreurs peut être
dirigée essentiellement par :
224
choisi) caractérisant chaque classe spectrale, typiquement la moyenne, l’écart-type, les
minimum et maximum et de la sorte obtenir la « signature spectrale » de chaque classe
spectrale à partir des différentes bandes spectrales de l’image considérée.
Les ROI décrivant une classe spectrale donnée doivent être spectralement
(relativement) homogènes : ils ne peuvent pas s’étendre sur plusieurs classes
différentes (Figure 144-gauche et centre) et ne peuvent pas inclure de
« mixels » (pixel mixte couvrant plusieurs classes d’occupation du sol). Typiquement
les pixels se trouvant à l’interface de 2 classes d’occupation du sol sont à éviter. Par
exemple, lors de la sélection des ROI sur un fleuve, ne pas considérer des pixels à
cheval ou se trouvant top proche de la berge car leur réflectance sera/risque d’être
influencée par la berge (Figure 144-droite).
Figure 144 : Exemples de mauvaises zones d’entrainement ou ROI couvrant plusieurs classes
d’occupation du sol différentes ou incluant des « mixels » ou « pixels mixtes ».
Il est souvent préférable que ces aires d’entrainement soient bien réparties sur toute
l’image. En effet, ceci permettra de prendre en compte :
o Une légère différence spectrale entre une même classe d’occupation du sol
dans des régions différentes (gradient ou classes pédologiques, variation des
espèces végétales ou de leur état de développement, etc).
o …
225
Il faut sélectionner un nombre suffisant de pixels dans les ROI d’une classe spectrale
donnée pour qu’ils soient REPRESENTATIFS de la classe spectrale considérée. Le
nombre de pixels nécessaire dépendra donc de la complexité spectrale de chaque
classe et de la complexité de l’image. Par exemple :
226
général préférable de créer une classe spectrale par mode/bosse/groupe de la
courbe (mais ceci pourrait être nuancé en fonction des « classes objectif » de la
classification).
Figure 145 : Courbes de la distribution de la fréquence des valeurs de pixels unimodale (vert) et tri-
modale (bleue). Axe X = valeurs des pixels, par exemple proportionnelles à de la réflectance. Axe y
= Nombre de pixels.
S’il y a des nuages et qu’il n’ont pas fait l’objet d’un masquage préalable, faites une
classe spectrale « Nuages ».
ATTENTION ! Avant de créer des ROI, il faut créer un « Jeu de bandes » et que ce jeu soit
repris dans l’onglet actif du menu « SCP > Jeu de bandes » : de cette manière, les
« signatures spectrales » seront calculées sur ce « Jeu de bandes ».
Si ce n’est pas encore fait, créez le « Jeu de bandes » se rapportant à l’image à classifier
via le menu QGIS :
Ce « Jeu de bandes » doit inclure toutes les bandes spectrales utiles à la classification de
l’image.
La création d’un « Jeu de bandes » a déjà été expliquée en détails à la section 8.6.3 page
178.
Cliquez sur l’onglet de ce « Jeu de bandes » pour l’activer via le menu « SCP > Jeu de
bandes »
227
Facultatif :
Classez les bandes du « Jeu de bandes » par ordre croissant via le bouton
Une fois le « jeu de bandes » créé, affichez une composition colorée de l’image satellite à
classifier dans QGIS comme indiqué à la section 8.6.3.2 page 180.
ATTENTION ! Avant de créer des ROI, il faut paramétrer les « Préférences » du projet QGIS
pour définir le Système de Coordonnées de Référence (SCR) qui sera attribué par défaut
aux fichiers vectoriels nouvellement créés.
Lors de la création d’un fichier de ROI dans SCP, qui est un fichier vectoriel, le système de
coordonnées de ce fichier sera automatiquement défini sur base des « Préférences »
paramétrées de votre projet QGIS.
Pour éviter un bug lors de la création des ROI dans SCP, il faut que le fichier vectoriel de ROI
automatiquement créé par SCP soit exprimé dans le même SCR que l’image à classifier et
que le projet QGIS lui-même (confer à ce propos la remarque importante de la section 8.3.2
page 131).
Paramétrez les préférences de votre projet QGIS comme souhaité via le menu QGIS
« Préférences » (Figure 146) :
Préférences > Options > SCR > SCR pour les couches > SCR par défaut pour les
couches : indiquer : « SCR du projet : EPSG : 32635… » (l’EPSG 32635 est valable
pour l’exemple du manuel. A adapter si vous travaillez sur une autre zone d’étude
utilisant un autre SCR).
Préférences > Options > SCR > SCR pour les couches > Quand une nouvelle couche
est créée ou… : indiquer « Utiliser le SCR du projet ».
228
Figure 146 : Paramétrage du menu QGIS « Préférences > Options > SCR > SCR pour les couches »
avant la création de ROI avec le plugin SCP.
Pour rappel, en cas de nécessité consultez l’aide en ligne, notamment la page dédié à la
description du panneau « Menu CP » disponible ici :
https://fromgistors.blogspot.com/p/user-manual.html?spref=scp
Allez dans le panneau « Menu SCP > Entrée données entrainement > Créez une
nouvelle donnée d’entrainement via le bouton > Enregistrez dans votre
répertoire de travail et donnez un nom, par exemple
« ROI_LANDSAT_5_CALIBRATION.scp » (Figure 147).
Attention, pour éviter toute perte de données, il ne faut pas éditer cette couche avec les
fonctions traditionnelles d’édition de QGIS mais uniquement via les outils du plugin SCP.
229
Figure 147 : Panneau « Menu SCP » du plugin SCP avec une série de ROI créés en vue de la
classification de l’image Landsat-5 (remarque : une classe spectrale « Nuages » aurait dû être
ajoutée à la liste des classes spectrales)
Ces boutons créent une « ROI temporaire » qu’il faudra ensuite sauvegarder.
230
Fenêtre « Menu SCP » > Entrée données entrainement > ROI & Liste des signatures >
et compléter « MC ID » et « Nom de MC », « C ID » et « Nom de C ».
MC = « Macro Classe »
o MC ID : un numéro unique pour chaque macro classe. Plusieurs ROI peuvent
avoir le même MC ID.
o Nom de MC : nom de la MC
C = « Classe spectrale »
o C ID : un numéro unique pour chaque classe spectrale. Plusieurs ROI peuvent
avoir le même C ID (auquel cas ils seront tous considérés pour définir la
signature spectrale de la classe « C ID »)
o Nom de C : nom de la classe spectrale
Cette attribution peut se faire également après la sauvegarde de la ROI temporaire
en cliquant-gauche 2 fois espacées sur la case à modifier dans le tableau des ROI.
Remarque : les lettres reprises dans la colonne « Type » dans le panneau des ROI signifie :
Clic sur (à cliquer après chaque dessin de ROI pour la sauvegarder et pouvoir
passer à la ROI suivante).
Les sélectionner dans la table des ROI puis clic sur . Cette fonction n’est
probablement pas utile pour la classification en tant que telle mais se révèle
intéressante si vous désirez réaliser un graphique des signatures spectrales avec une
seule signature par classe spectrale (confer section 8.8.7 page 237).
232
sous un autre nom
o ATTENTION ! L’utilisation de cet outil semble parfois déboucher sur un
message d’erreur « Warning [22]: The coordinate system of training input is
different from the input image. Please create a new training input”. Ce
message apparait lorsque, après avoir créé un nouveau ROI1.scp via le bouton
ET avoir créé un ou plusieurs polygones ROI dans ce fichier ROI1.scp, l’on
tente de soit créer un autre ROI2.scp avec ou d’ouvrir un ROI3.scp
préexistant avec le bouton . Une solution est de fermer le projet QGIS et de
le ré-ouvrir : vous pourrez alors ré-ouvrir des ROI.scp préexistant sans
problème ou en créer un nouveau.
Pour rappeler des ROI créés et sauvegardés précédemment ou changer de jeu de ROI
lorsque plusieurs jeux de ROI sont disponibles :
Remarque : les signatures spectrales des classes spectrales doivent avoir été calculées
(confer section 8.8.4.2 page 229).
Menu SCP > Traitement de bande > Classification > (Figure 148)
o Sélectionner le numéro du jeu de bandes sur lequel s’applique la classification
o Choisir MC ID ou C ID en fonction du type de classe à utiliser pour la
prévisualisation
o Choisir l’algorithme de classification
Dans le menu « Classification preview » de la barre d’outils SCP
o Cliquez sur puis sur l’endroit de l’image pour lequel vous voulez pré-
visualiser la classification et la prévisualisation apparait (Figure 148)
o Découvrez éventuellement l’utilité des autres boutons de prévisualisation en
les testant
o A tester avec MC ID puis C ID coché
233
Figure 148 : Aperçu de la classification pour une partie de l’image grâce à l‘outil de prévisualisation
de la classification du plugin SCP.
Menu SCP > Traitement de bande > Classification > interface à paramétrer comme
désiré, et dans le cadre de cet exercice (Figure 149) :
o Sélectionner le numéro du jeu de bandes sur lequel s’applique la classification
o Cochez « MC ID » pour que la classification représente les MC et pas les
classes spectrales.
o Choisir la méthode par « Maximum de vraisemblance » (à vous de tester les
autres méthodes et de comparer la qualité des résultats : il semblerait que la
méthode « Plus proche voisin » donne souvent de meilleurs résultats)
o Laissez les autres paramètres par défaut.
o Cliquez sur
o Choisir le répertoire et donner un nom à l’image en sortie
234
Figure 149 : Menu SCP > Traitement de bande > Classification
La classification se réalise.
Landsat-5 : ~4 minutes
235
Figure 150 : Premier essai de classification (de mauvaise qualité) de l’image Landsat-5
Figure 151 : Premier essai de classification (de mauvaise qualité) de l’image Sentinel-2 mosaïquée
236
8.8.7. Etude de la séparabilité spectrale des classes spectrales
La séparabilité spectrale des différentes classes spectrales peut être analysée afin de
quantifier/d’objectiver le caractère séparable ou non séparable de ces classes. 2 méthodes
sont disponibles dans QGIS :
https://fromgistors.blogspot.com/p/user-manual.html?spref=scp
237
8.8.7.1. Graphique des signatures spectrales
Un graphique des signatures spectrales peut se réaliser en
Sélectionnant les ROI d’intérêt dans l’onglet « Entrée données d’entrainement »du
panneau « Menu SCP »
Puis en cliquant sur le bouton « Ajouter les signatures spectrales surlignées au
graphique » pour créer les graphiques des signatures spectrales
Une fenêtre « SCP : Graphique signature spectrale » s’ouvre et son onglet « Graphique »
présente les signatures spectrales (Figure 152). Le bouton permet d’ajuster le graphique
correctement.
Figure 152 : Signatures spectrales de différentes classes d’occupation du sol dans le plugin SCP de
QGIS (après fusion des ROI correspondant à une même classe spectrale)
238
8.8.7.2. Indices de séparabilité spectrale
4 méthodes de mesures de la séparabilité spectrale sont proposées dans le plugin SCP :
239
8.8.8. Classification rapide des nuages en vue de la réalisation d’un
masque de nuages amélioré
! Méthode proposée dans le cadre de ce manuel :
Landsat-5 : à ne pas faire
Sentinel-2 : la méthode décrite dans cette section a été utilisée pour créer les 2
masques de nuages améliorés pour les 2 images Sentinel-2. Ces masques sont
disponibles dans le dossier « …\DATA\IMAGES\SENTINEL_2\ » et s’appellent
SENTINEL_2_MASQUE_NUAGES_20180612.tif et
SENTINEL_2_MASQUE_NUAGES_20180707.tif et prennent pour valeurs 1 pour les
nuages et 2 pour le reste.
Comme mentionné à la section 8.6.9.3 page 189, un masque de nuages peut être réalisé par
classification de l’image présentant les nuages et devant être masquée.
Voici, pour information, la procédure qui a été suivie pour réaliser les masques de nuages
améliorés des images Sentinel-2.
240
Un aperçu de la classification peut être utilisé pour vérifier la qualité obtenue. L’outil
« Map swipe tool » (plugin de QGIS) peut être utile pour comparer la
classification avec l’image de départ
La classification est ensuite produite et sa qualité vérifiée visuellement (~1 minute)
Figure 154 : Création des ROI de calibration pour la classification des nuages (ancienne interface de
SCP)
La classification, si réussie, permet d’identifier les nuages beaucoup mieux que le masque de
nuages fourni avec les images Sentinel-2 1C. Cette classification ne sera cependant
certainement pas parfaite. Il convient dès lors de la raffiner, en fonction des défauts qu’elle
présente (trop ou trop peu de pixels classifiés comme nuage).
Pour raffiner cette classification, les techniques suivantes ont été utilisées (à adapter selon
la qualité de classification obtenue):
241
Figure 155 : Outils de tamisage de QGIS (pour fusionner les petits groupes de pixels isolés avec la
classe voisine majoritaire)
Figure 156 : Outil de dilation de classification utilisé pour dilater la classe nuage (ancienne interface
de SCP)
242
L’image résultante est inspectée visuellement pour estimer la qualité d’identification des
nuages.
Faire des tests en modifiant le paramétrage si le résultat n’est pas satisfaisant.
L’identification des nuages est de loin meilleure à celle fournie dans le fichier vectoriel initial
(Figure 157).
Image Sentinel-2 20180707 bande bleue Image Sentinel-2 20180707 bande bleue
avec masque de nuages vectoriel fourni
avec l’image
Figure 158 : Masques de nuages améliorés résultant de la classification et raffinement pour les 2
images Sentinel-2
243
8.9. Amélioration de la classification par reclassification
conditionnelle et données auxiliaires
8.9.1. Introduction
Une classification existante peut être modifiée/améliorée par reclassification de certains
pixels sur base d’une comparaison avec les valeurs de pixels d’autres rasters.
Il peut être intéressant dans certains cas de reclassifier une classification sur base d’une
autre classification de la même zone mais d’une autre date. Par exemple la classe « forêt »
de la classification correspondant au temps 2 pourrait être reclassifiée en :
« forêt ancienne » si une forêt existait déjà au temps 1 (40 ans plus tôt par exemple)
« forêt jeune » s’il n’y avait pas de forêt au temps 1
244
8.9.3. Reclassifications conditionnelles en pratique
Pour reclassifier des pixels d’une classification en fonction des valeurs que prennent un ou
plusieurs autres rasters, utilisez le :
Menu QGIS SCP > Calcul de bande > Règles de décision (Figure 159)
Il faut que les rasters contenant les paramètres utiles à la reclassification conditionnelle
soient présents dans le projet QGIS et donc visible dans la « Liste de bandes » de cette
fenêtre.
(Figure 160).
Figure 159 : Reclassification conditionnelle d’une classification avec l’outil « Règles de décision » de
SCP (raster19 = la classification, raster59 = MNT et raster29 = raster de proximité aux rivières (zone
tampon de 60 m) (ancienne interface de SCP)
245
Figure 160 : Classification avant (gauche) et après (droite) reclassification conditionnelle faisant
apparaitre les forêts d’altitude (vert clair en montagne) et les forêts riveraines (encadré)
Figure 161 : 1ère classification de l’image Sentinel-2 mosaïquée avec reclassification conditionnelle
pour les forêts d’altitude et riveraines
246
8.10. Post-traitements d’une image classifiée pour en
améliorer le rendu visuel ou la raffiner
Plusieurs opérations de post-traitement peuvent être appliquées sur une image classifiée en
vue d’en améliorer le rendu visuel, de corriger certaines erreurs ou de parfaire la
classification, par exemple :
247
Figure 162 : Exemple de l’effet de l’application d’un filtre majoritaire ou tamisage sur une
classification: à gauche : classification originale avec petits groupes de pixels isolés (=effet
ressemblant à un mélange de grains de poivre et de sel); à droite : classification filtrée (filtre
majoritaire) ou tamisée avec les petits groupes de pixels isolés qui ont été fusionnés dans la classe
majoritaire voisine (outils utilisé : « Tamiser » de QGIS avec un seuil de 4).
248
8.11. Evaluation de la précision de la classification par
« matrice de confusion » ou « table de contingence »
Une comparaison visuelle entre la classification et une donnée de référence qui peut
être l’image ayant servi à la classification, une autre classification, une autre données
qui contient une information sur l’occupation du sol.
o Avantage : rapide, peut donner un rapide aperçu général de la qualité de la
classification
o Désavantage : peut être assez subjective et difficile à mener en détails si la
zone à observer est grande et les classes nombreuses et présentant une
configuration spatiale complexe, par exemple via l’enchevêtrement de
nombreuses petites classes sous forme de zones irrégulières (versus
répartition claire de peu de classes sous forme de grandes zones régulières)
L’utilisation d’une « matrice de confusion » ou « table de contingence »
o Avantage : est plus objective car elle permet de quantifier la précision
obtenue à travers une série d’indicateurs chiffrés (précision globale, précision
classe par classe, etc) (confer plus loin). A la condition que l’échantillonnage
soit correctement réalisé, la matrice de confusion permettra une estimation
de la qualité réelle de la classification (pas de sous-/sur-estimation de la
qualité)
o Désavantage : prend plus de temps à réaliser car il faut créer/échantillonner
les zones de validation et leur attribuer une classe
Les ROI de CALIBRATION : elles sont utilisées pour CALIBRER la classification, c’est-à-
dire pour calculer les signatures spectrales de chaque classe spectrale qui seront
ensuite utilisées pour attribuer une classe à chacun des pixels de l’image au moment
de la classification
249
Les ROI de VALIDATION : elles sont utilisées pour VALIDER/ESTIMER LA QUALITE de
la classification, via la production d’une matrice de confusion, et par COMPARAISON
(au sein des ROI de validation) entre:
o l’occupation du sol « classifiée » ou « prédite » de la CLASSIFICATION et
o l’occupation du sol associée aux zones de validation lors de la création des
ROI de VALIDATION
La matrice de confusion est donc calculée UNIQUEMENT sur base des zones de l’image
classifiée identifiées dans les ROI de VALIDATION, d’où l’importance que celles-ci soient
suffisamment nombreuses que pour être représentatives de la qualité de l’ensemble de la
classification.
250
Figure 163 : Principes de calibration et validation d’une classification supervisée
251
8.11.1.2. Stratégies d’échantillonnage des zones de validation (ROI de
validation)
Stratégies d’échantillonnage
Echantillonnage dirigé : les points sont positionnés à des endroits choisis par
l’opérateur. Cette méthode n’est pas idéale car l’opérateur aura naturellement
tendance à choisir pour la validation des zones qu’il peut facilement photo-
interpréter et qui sont probablement plus facilement classifiables que d’autres. Ceci
aura pour conséquence de fausser les résultats de la validation en présentant une
précision de classification surestimée. Il est préférable que les points soient
positionnés aléatoirement.
Echantillonnage aléatoire simple : les points sont positionnés aléatoirement dans
l’espace
Echantillonnage aléatoire stratifié (méthode utilisée dans le cadre de cet exercice) :
un certain nombre de points, en partie proportionnel à la superficie de chaque classe,
est positionné aléatoirement au sein de chaque classe
Echantillonnage régulier ou systématique sous forme de grille régulière ou de
transects : les points sont positionnés selon une grille ou un transect (ligne)
Echantillonnage spatialement conditionné par l’accessibilité des zones sur le
terrain : par exemple, positionnement des points le long de/à proximité des voies de
communication (route, chemin, rivière) ou en évitant des zones d’insécurité, ou trop
difficilement accessibles.
Il faut veiller à ce que les zones de validation soient suffisamment nombreuses et réparties
de manière homogène sur la carte AFIN QUE CELLES-CI PERMETTENT DE CALCULER UNE
MATRICE DE CONFUSION QUI SOIT BIEN REPRESENTATIVE DE LA PRECISION DE
CLASSIFICATION DE L’ENSEMBLE DE L’IMAGE. Par exemple, si vous utilisez une seule petite
zone de validation pour une classe donnée et que par chance tous les pixels de cette zone
sont bien classés, la précision de classification obtenue pour cette classe par la matrice de
confusion sera de 100% alors qu’il est probable que d’autres pixels appartenant à cette
classe à d’autres endroits de l’image soient eux mal classifiés. Dans ce cas votre zone de
252
validation ne sera pas représentative de l’ensemble de la classe d’occupation du sol et
l’estimation de la précision de classification de cette classe sera donc faussée.
Classe inconnue ?
S’il existe des zones de validation pour lesquelles la classe d’occupation du sol réelle n’est
pas connue (photo-interprétation trop difficile, pas d’information venant du terrain, etc) il
devient impossible d’en réaliser la validation. A l’opérateur alors d’estimer si ces zones sont
négligeables (par exemple parce qu’elles couvrent un très faible pourcentage de la surface
de l’image), ou pas, et de documenter les limitations de la méthode de validation utilisée.
253
classification simpliste ne considérant que quelques classes fortement
contrastées (par exemple, les classes sol nu, végétation, eau)
Du nombre de classes considérées
o Une précision globale de 50% pour une classification ne considérant que 2
classes, équivaut à une classification aléatoire des pixels (1 chance sur 2
d’être correctement classifié)
o Une précision globale de 50% pour une classification considérant 10 classes,
représente déjà un certain niveau de qualité de classification
De la surface relative des différentes classes
o Par exemple, une petite imprécision sur une classe couvrant 90% de la zone
étudiée aura peut-être plus d’importance qu’une grande imprécision sur une
classe ne couvrant que 5% de la zone.
https://fromgistors.blogspot.com/2019/09/Accuracy-Assessment-of-Land-Cover-
Classification.html
https://fromgistors.blogspot.com/2019/09/Accuracy-Assessment-of-Land-Cover-
Classification.html
254
Dans le cadre de cet exercice, par soucis de simplification et de rapidité, une approche
simplifiée est proposée.
Le nombre de point total à utiliser pour la validation (N) sera arbitrairement fixé à 100 (ce
qui peut être considéré comme très peu). Dans une situation réelle il faudra suivre les
recommandations du tutoriel de SCP pour calculer le nombre de points total (N) !
La répartition de ces 100 points au sein des différentes macro-classes que l’on veut valider
se fera selon la règle proposée dans le tutoriel de SCP, à savoir :
𝑁
+ 𝑁 ∗ 𝑊𝑖
𝑁𝑖 = ( 𝑐 )
2
Où :
Ni correspond à la valeur moyenne entre une distribution égale du nombre de points entre
chaque classe (N/c) et une distribution directement proportionnelle à la superficie de
chaque classe (N*Wi).
N = 100
c = nombre de MC = 8 (ne pas considérer la classe « 0 - Unclassified »)
Wi peut être obtenu via le
o Menu SCP > Post-traitement > Rapport de classification
Sélectionner la classification à valider
Lancer le calcul et choisir le répertoire de sortie + donner un nom au
fichier texte csv produit, par exemple
« Rapport_de_classification_Landsat5.csv »
Wi correspond à la colonne « Percentage % » du fichier csv produit.
Le nombre de pixels par classe (Ni) peut ensuite être facilement calculé selon la
formule reprise ci-dessus et implémentée dans la feuille Excel « Calcul
échantillonnage stratifié pour validation.xlsx » disponible dans le dossier « DATA »
de cet exercice (Figure 164). Vous adapterez cette feuille Excel en fonction du
nombre de classes considérées, de leur Wi et éventuellement du N que vous avez
choisi ou calculé.
255
Figure 164 : Calcul du nombre de points de validation par classe directement en partie
proportionnel à la superficie de chaque classe (échantillonnage stratifié) dans la feuille Excel «
Calcul échantillonnage stratifié pour validation.xlsx »
o Fenêtre « Menu SCP » > Entrée données d’entrainement > ROI & Liste des
signatures : renseigner “MC ID” = 0, “Nom de MC” vide, “C ID” = 0, et “Nom
256
de C” vide. Ces valeurs seront automatiquement reprises pour tous les ROI de
validation générés dans l’étape suivante. Ces valeurs devront être ensuite
éditées via photo-interprétation. Indiquer une valeur de 0 ou vide permettra
de distinguer plus facilement les ROI déjà édités de ceux non encore édités.
o SCP > Outils basic > Création de ROI multiples (Figure 165)
Paramétrez la fenêtre comme désiré pour la création des ROI de
validation d’une classe donnée. Un exemple est donné dans la Figure
165 et sa légende.
« Nombre de points » correspond au Ni calculé précédemment.
Cliquez sur ce qui créé les points et les affiche dans la
table
Répéter l’opération pour toutes les macro-classes en adaptant à
chaque fois le nombre de points à créer et la valeur du raster de
classification qui désigne la classe concernée (Figure 164)
Décocher « Calculer les signatures » (car pas besoin)
Lancer la création
(Attention ! Il faut que la classification et le projet QGIS soient dans le
même SCR !, Confer section 8.3.2 page 131)
Les points créés viennent s’ajouter dans la Fenêtre « Menu SCP » >
Entrée données d’entrainement > ROI & Liste des signatures avec dans
la colonne « Type » : ROI ( = sans signature spectrale calculée) (Figure
166) mais avec un « MC ID » et « C ID » = 0 et sans « Nom ».
257
Figure 165 : Fenêtre « SCP > Outils basic > Création de ROI multiples ». La paramétrisation de cet
exemple a créé 14 points dans la classe correspondant à une valeur de pixel de 6 dans le raster de
la classification qui est repris dans le « Jeu de bandes » numéro 3.
Affichez dans QGIS une composition colorée de l’image utilisée pour produire la
classification à valider ou une autre source d’information qui pourrait être utilisée
pour attribuer les classes aux ROI de validation
Fenêtre « Menu SCP » > Entrée données d’entrainement > ROI & Liste des
signatures (Figure 166)
o Double-cliquez sur une ligne de la table des ROI de validation pour zoomer
automatiquement sur ce ROI
o Dé-zoomez ensuite jusqu’à être capable d’identifier par photo-interprétation
la classe correspondante à ce pixel
o Entrez les « MC ID » (et « C ID » éventuellement, ce qui n’est pas nécessaire si
vous faites une validation sur les MC uniquement) dans la table des ROI pour
le ROI considéré en cliquant UNE FOIS de plus sur la cellule à éditer
o Répétez l’opération pour tous les ROI de validation
o Si un pixel est un mixel (pixel composé de plusieurs occupations du sol) une
solution approximative est de considérer la classe la plus importante en
surface dans le pixel
o Si la photo-interprétation n’est pas possible, soit attribuez la classe la plus
probable, soit attribuez un nouveau code MC ID qui identifiera les pixels non
photo-interprétables (méthode non encore testée).
258
o Le fichier de ROI est sauvegardé automatiquement si la case correspondante
est cochée.
Figure 166 : Attribution des identifiants des macro-classes (« MC ID ») aux points de validation
(générés de manière aléatoire stratifiée) par photo-interprétation. Les points de validation
correspondent aux pixels noirs dans l’interface de visualisation spatiale de QGIS. Dans cette
illustration, seuls 2 ROI de validation ont déjà reçu un MC ID (de 3 et 5).
259
Figure 167 : Paramétrage de la création de la matrice de confusion
Figure 168 : Matrice de confusion produite et les statistiques dérivées. (Remarque : cette matrice
de confusion n’a pas été créée avec un N (nombre de points de validation théorique total (confer
section 8.11.2.3)) = 100 mais avec un N = 50 lors d’un essai préalable. Normalement un N = 100
aurait dû être utilisé pour cette illustration pour être raccord avec le reste des explications).
(Ancienne interface de SCP)
La matrice exprimée en « pixel count », c’est-à-dire en nombre pixels, est plus facilement
interprétable pour une première approche.
La matrice « area based error matrix » est exprimée en pourcentage de la surface totale de
l’image classifiée et c’est cette matrice qui est utilisée pour dériver les paramètres
d’évaluation de la précision de la classification dans SCP. Quelques informations sont
disponibles à ce propos dans la documentation officielle en ligne.
260
La manière d’interpréter une matrice de confusion est présentée en détails dans la section
suivante 8.11.2.7 page 262.
261
8.11.2.7. Matrice de confusion : construction et interprétation, un exemple
illustré
Considérons une image satellite couvrant une zone d’étude caractérisée par 3 types d’occupation
du sol à identifier par classification:
Eau (rivière, étang)
Prairie
Matrice de confusion : construction et interprétation (1/4)
Forêt
L’image est classifiée et des zones de validation (ou « Ground truth Region Of Interest ») sont
créées pour chacune des classes d’occupation du sol (Figure 169). Ces zones correspondent à des
régions pour lesquelles l’occupation du sol est connue avec certitude, soit via photo-
interprétation, soit via enquête de terrain avec relevés GPS par exemple. La précision de la
classification de l’ensemble de l’image est ensuite évaluée en comparant, dans ces zones de
validation, l’occupation du sol réelle (définie par l’opérateur) avec celle identifiée par le processus
de classification.
Figure 169 : Image classifiée de la zone d’étude. Les encadrés noirs correspondent aux zones de
validation ou « Ground truth Region Of Interest ». Les classes auxquelles sont assignées chacune de ces
zones de validation sont mentionnées au-dessus des encadrés.
262
La matrice de confusion (ou « table de contingence ») (Figure 170) permet d’évaluer la précision
de la classification (de manière globale et pour chacune des classes). Elle est calculée à partir des
pixels correctement/incorrectement classifiés dans les zones de validation mais est sensée être
représentative de la précision de la classification de l’ensemble de l’image. D’où l’importance de
créer suffisamment de zones de validation, correctement réparties sur l’ensemble de l’image, et
reprenant l’ensemble des réalités spectrales de chacune des classes de l’image (exemple : créer
Matrice de confusion : construction et interprétation (2/4)
des zones de validation de la classe forêt dans les différents massifs forestiers de l’image qui
présentent probablement de petites nuances spectrales (variation de l’exposition, des essences,
de l’ensoleillement, etc).
Zones de validation
Erreurs de Précision
Eau Prairie Forêt Total
commission d’utilisation
Eau 17 0 1 18 1/18 (6%) 17/18 (94%)
Classification
263
Les erreurs d’omission :
L’erreur d’omission d’une classe donnée est calculée en divisant le nombre total de pixels
incorrectement classifiés dans les zones de validation de cette classe (« omis » de la classe
étudiée) par le nombre total de pixels des zones de validation de cette classe. Dans la matrice de
confusion, ce calcul s’effectue sur base de la colonne correspondant à la classe étudiée.
« Eau » : 0 pixel n’ont été mal classifiés sur les 17 pixels des zones de validation de la classe
« Eau », soit 0/17 = 0% d’erreur d’omission (0% des pixels des zones de validation de la
classe « Eau » n’ont été omis de la classe « Eau »)
Matrice de confusion : construction et interprétation (3/4)
« Prairie » : 6 pixel ont été mal classifiés (dans la classe « Forêt ») sur les 33 pixels des
zones de validation de la classe « Prairie », soit 6/33 = 18% d’erreur d’omission (18% des
pixels des zones de validation de la classe « Prairie » ont été omis de la classe « Prairie »
(et classifié dans la classe « Forêt ») alors qu’ils auraient du être classifiés dans la classe
« Prairie »)
« Forêt » : 4 pixel ont été mal classifiés (dans les classes « Eau » et « Prairie») sur les 30
pixels des zones de validation de la classe « Forêt », soit 4/30 = 13% d’erreur d’omission
(13% des pixels des zones de validation de la classe « Forêt » ont été omis de la classe
« Forêt » (et classifié dans les classes « Eau » et « Prairie») alors qu’ils auraient du être
classifiés dans la classe « Forêt »)
Les précisions de production :
La précision de production d’une classe donnée est calculée en divisant le nombre total de pixels
correctement classifiés dans les zones de validation de cette classe par le nombre total de pixels
des zones de validation de cette classe. Elle est égale à (1-l’erreur d’omission).
264
« Prairie » : 3 pixels ont été mal classifiés (en tant que « Prairie ») alors qu’ils
appartiennent (d’après les zones de validation) à d’autres classes (ici, la classe « Forêt »),
sur un total de 30 pixels classifiés dans la classe « Prairie » dans l’ensemble des zones de
validation de l’ensemble des classes, soit 3/30 = 10% d’erreur de commission (10% des
pixels classifiés dans la classe « Prairie » appartiennent en réalité à d’autres classes (ici la
classe « Forêt»)
« Forêt » : 6 pixels ont été mal classifiés (en tant que « Forêt ») alors qu’il appartiennent
(d’après les zones de validation) à d’autres classes (ici, la classe « Prairie »), sur un total de
32 pixels classifiés dans la classe « Forêt » dans l’ensemble des zones de validation de
l’ensemble des classes, soit 6/32 = 19% d’erreur de commission (19% des pixels classifiés
Matrice de confusion : construction et interprétation (4/4)
dans la classe « Forêt » appartiennent en réalité à d’autres classes (ici la classe « Prairie»)
Les précisions d’utilisation :
La précision d’utilisation d’une classe donnée est calculée en divisant le nombre total de pixels
correctement classifiés dans cette classe dans l’ensemble des zones de validation de cette classe
par le nombre total de pixels classifiés dans cette classe dans l’ensemble des zones de validation
de l’ensemble des classes. Elle est égale à (1-l’erreur de commission).
La précision d’utilisation mesure la probabilité qu’un pixel appartienne réellement à une certaine
classe lorsque le processus de classification a classifié ce pixel dans cette classe.
Autrement dit :
Probabilité qu’un pixel qui a été classifié dans une classe « a » appartienne réellement (dans le
monde réel) à cette classe « a ».
Autrement dit :
« En tant qu’utilisateur d’une classification (celui qui lit/utilise la classification), je me demande,
pour un pixel classifié dans une classe « a », quelle est la probabilité qu‘il appartienne réellement
(dans le monde réel) à la classe « a ». »
Erreur globale de la classification :
L’erreur globale de la classification est calculée en divisant le nombre de pixels mal classifiés dans
l’ensemble des zones de validations de l’ensemble des classes par le nombre de pixels de
l’ensemble des zones de validation de l’ensemble des classes. L’erreur globale = (1-la précision
globale).
Erreur globale : 10 pixels mal classifiés dans l’ensemble des zones de validation de
l’ensemble des classes divisé par les 80 pixels de l’ensemble des zones de validation de
l’ensemble des classes = 12.5 % d’erreur globale.
Précision globale de la classification :
La précision globale de la classification est calculée en divisant le nombre de pixels correctement
classifiés dans l’ensemble des zones de validation de l’ensemble des classes par l’ensemble des
pixels de l’ensemble des zones de validation de l’ensemble des classes.
Précision globale : 70 pixels correctement classifiés dans l’ensemble des zones de
validations de l’ensemble des classes divisé par les 80 pixels de l’ensemble des zones de
validation de l’ensemble des classes = 87.5 % de précision globale. La précision globale =
(1- l’erreur globale)
Plus d’infos sur les matrices de confusion : Centre Canadien de télédétection
http://cct.rncan.gc.ca/glossary/index_f.php?id=3124
265
8.12. Matrice de transition
(Remarque : le menu SCP > Post-traitement > Changement d’occupation des sols produit les
mêmes résultats mais moins bien présentés, sans matrice)
266
Figure 171 : Matrice de transition entre les images Landsat-5 et Sentinel-2 (ancienne interface de
SCP)
267
8.13. Graphiques exprimant l’évolution de l’occupation du
sol
2 types de graphiques sont typiquement utilisés pour représenter les résultats issus d’une
analyse diachronique de l’occupation du sol d’une région.
Ces graphiques peuvent facilement être réalisés à partir des résultats du menu :
Les fichiers CSV pourront être éventuellement édité/remanié et les graphiques réalisés :
La feuille Excel utilisée pour produire les graphiques ci-dessous est fournie dans le répertoire
« …\5_CLASSIFICATION_SUPERVISEE_ANALYSE_DIACHRONIQUE\DATA\GRAPHIQUES DE
TRANSITION.xlsx »
268
8.13.1. Histogramme par classe
Double bâtonnet pour chaque classe.
Figure 172 : Graphique de type histogramme qui permet de comparer l’évolution des pourcentages
de la superficie totale de la zone d’étude des différentes classes d’occupation du sol entre 2 dates.
Figure 173 : Graphique de type camembert qui permet de comparer l’évolution des pourcentages
de la superficie totale de la zone d’étude des différentes classes d’occupation du sol entre 2 dates.
269
8.14. Réalisation de cartes d’occupation du sol - Mise en
page cartographique
Les images classifiées, ayant subi éventuellement l’un ou l’autre post-traitement, pourront
ensuite faire l’objet d’une mise en page cartographique si celles-ci doivent être présentées
dans un rapport par exemple.
La mise en page cartographique consiste en l’habillage des images classifiées à l’aide d’une
série d’éléments tels légende, barre d’échelle, flèche du Nord, graticule de
géoréférencement, carte de localisation ou rectangle d’emprise et éléments textuels tels
titre, auteur, date de création de la carte, système de coordonnées, source et date des
données, etc
L’interface de mise en page cartographique d’un projet QGIS est accessible via le menu :
Projet > Nouvelle mise en page… > Créer un titre de mise en page
Pour plus d’information et indications sur la mise en page cartographique, le lecteur est
invité à consulter
Lors de l’insertion des cartes issues de l’analyse diachronique dans un rapport, veiller à ce
que:
270
8.15. Discussion et limitations de l’approche
Quelques éléments importants sont brièvement repris ci-dessous :
271
9. Analyse de données SAR SENTINEL-1
(télédétection active) pour la
cartographie d’inondation dans SNAP (en
anglais)
SNAP ( https://step.esa.int/main/download/snap-download/)
ATTENTION ! Le logiciel SNAP de l’ESA, gratuit, nécessite pour certaines
Logiciel ou
opérations une certaine « puissance » de l’ordinateur que nombres «
application
d’ordinateurs étudiants » n’ont pas, ce qui résulte parfois en des temps
de traitements à rallonge, voire des bugs.
Temps
~ 1h30 à 2h00
approximatif
272
Aucun. Cependant, certaines connaissances théoriques préalables sur la
Prérequis télédétection RADAR permettront à l’apprenant de mieux comprendre
les manipulations réalisées.
https://www.youtube.com/watch?v=derOXkPCH80
273
o Mode VV
Polarization mode of the SAR signal: Vertical emitted by the sensor
and Vertically received by the sensor
o Sensor mode IW
Interferometric Wide acquisition mode (Figure 174): confer here for
more information: https://sentinel.esa.int/web/sentinel/user-
guides/sentinel-1-sar/acquisition-modes
o Product type GRD
Ground Range Detected
More information here:
https://sentinel.esa.int/web/sentinel/missions/sentinel-1/data-
products
Remarque: in GRDH (confer image filename), the H stands for the
resolution class, which is High (H) in this case)
274
9.4. Download of the 2 Sentinel-1 images (not in the video)
Go to the “Copernicus Open Access Hub” https://scihub.copernicus.eu/dhus/#/home
Create an account and then log in
Search and download the 2 images
o You can use the image name (confer below) without “.zip” extension in the
search menu to easily find the images, one by one (Figure 175)
o The 2 images used are (table below):
Figure 175 : Copernicus Open Access Hub with one of the 2 Sentinel-1 images identified over
Myanmar
275
9.5. SNAP software - Download and installation (not in the
video)
You can download the free SNAP software here:
o https://step.esa.int/main/download/snap-download/
o Choose the « Sentinel Toolboxes » or “All Toolboxes” in the table and the
download adapted to your operating system
Install SNAP on your computer and make the updates after installation
Figure 176 : Outlook of the SNAP (STEP) software download web page
276
9.6. Opening of the 2 Sentinel-1 images in SNAP (video
timing: 0 min 30 sec)
Simply drag and drop the zip images to the “Product Explorer” tab of SNAP (there is
no need to unzip the images!)
In the “World View” tab, you can see that the images are over a part of Myanmar
In the “Product Explorer” tab, expand for each image the sub-folder “Bands >
Amplitude_VV” and double-clic on “Amplitude_VV” which will open the images in
the viewer
o Open first the Archive image in viewer 1
o Open then the Crisis image in viewer 2
Remarque: notice that there are 2 bands per image, Amplitude and Intensity.
Definition:
o Amplitude_VV image
“A SAR signal contains amplitude and phase information. Amplitude is
the strength of the radar response” (source:
https://sentinel.esa.int/web/sentinel/user-guides/sentinel-1-
sar/product-overview/interferometry).
“Measure of the strength of a signal, and in particular the strength or
height of an electromagnetic wave (units of voltage).” (source:
https://earth.esa.int/handbooks/asar/CNTR5-2.html)
o Intensity_VV image
Intensity = Amplitude_VV * Amplitude_VV
Observe the images (Figure 178):
o You can zoom or move into the images by using the mouse wheel or the
panning tool
o In the crisis image:
Flooded area: appears darker because of a low backscatter due to
specular reflection over the smooth water surfaces: the signal get
reflected away from the sensor.
The surrounding areas: are much rougher and look brighter.
Images distortion is visible in the mountain areas of the image due to
the layover and foreshortening effects (Figure 177)
Foreshortening Layover
277
Figure 177 : Illustration of the Foreshortening and Layover effects: geometric distortions in
mountainous regions (Source: https://earth.esa.int/web/guest/missions/esa-operational-eo-
missions/ers/instruments/sar/applications/radar-courses/content-2/-
/asset_publisher/qIBc6NYRXfnG/content/radar-course-2-slant-range-ground-range#sixteen)
278
9.8. Visual analysis of the 2 images side by side (video
timing: 3 min 09 sec)
Menu: Window > Tile evenly
You should synchronize the 2 windows (and optionally the mouse pointer) by using,
in the “Navigation” tab of the bottom left panel, the 2 related tools button
Figure 178 : The archive/normal situation (left) and crisis/flood (right) Sentinel-1 images over
Myanmar displayed evenly in SNAP interface
279
9.9. Analysis of images Metadata (video timing: 3 min 15
sec)
In the “Product Explorer” tab, open the subfolders “Metadata > Abstracted
Metadata” and double-click on it to open it (Figure 179). You can see that:
o Both images were acquired in the same geometry which is essential for flood
mapping:
PASS = ASCENDING: satellite ascending
incidence_near (30.9°) & incidence _far (46.1°): same incidence angle
o PRODUCT_TYPE = GRD
o ACQUISITION_MODE = IW
o mds1_tx_rx_polar = VV (and no info for mds2…, mds3…, mds4…)
o etc
There is a shift in Azimuth (in the North-South direction, the 2 images do not
perfectly cover the same area) but this is not a problem
To know what image is open in what tab: look at the number between [] both in the
“Product Explorer” tab and in the visualization windows
Figure 179 : Outlook of some Abstracted Metadata of one Sentinel-1 image. The […] means that
some metadata are skipped from this illustration
280
9.10. Spatial subset of the images, on an area common to
the 2 images, where the main floods occur (video
timing: 3 min 53 sec)
Identify in the navigation tab an area which is common to both images: zoom and
pan in order to display that common area in the viewers (Figure 178)
Click, in the geographical interface, on the image you want to process in order to
select/active it, in order to use this image current extent for the spatial subset
Menu: Raster > Subset (video timing: 4 min 03 sec)
o The subset corresponds to the extent of the viewer. You can visually check
that in the tool interface
o Click OK
o Repeat the operation for both images
281
o Keep the most important information: date product, processing level etc
o Remove the automatically created prefix (“subset_0_of_”)
o And give also a suffix “crop” (or “subset”) in order to easily identify the
output image as the subset image
o Select your working directory where to store the output images
o The output filenames of the 2 images are, in the video:
S1A_IW_GRDH_1SSV_20150320T114745_crop.zip
S1A_IW_GRDH_1SSV_20150904T114747_crop.zip
o Click on Save
The 2 subseted images are now present in you working directory
282
9.12. Multilooking (video timing: 6 min 30 sec)
Purpose: to reduce the speckle and reduce the dimension of the image to speed up
the processing time. This will also result in a loss of accuracy (spatial resolution) but
this is not a problem in this particular case because the flooded area is large and we
are not particularly interested in having a high resolution flood map.
Menu: Radar > SAR Utilities > Multilooking (this directory is a bit different than in
the tutorial video due to different SNAP software versions) (Figure 181)
o In the “I/O Parameters” window
Keep default added suffix “ML” in the output filename
Select your working directory as output directory
o In the “Processing Parameters” window
Set the Number of Range/Azimuth Looks to 3 by 3
o Click on Run
o Repeat the process for the other image
(Remarque: after this operation, the subseted images cannot be synchronized well
with the original images because a size reduction has been applied, reason why all
products were closed in the previous “intermediate steps”)
283
9.13. Radiometric Calibration (video timing: 7 min 35 sec)
Calibration is essential to enable the comparison of the 2 images
Calibration transforms Digital Numbers (DN) to a physical quantity which is in this
case “sigma0 backscatter”
Definition:
o “Backscatter: […] The normalised measure of the radar return from a
distributed target is called the backscatter coefficient, or sigma nought , and
is defined as per unit area on the ground.”
o “Sigma Nought: Scattering coefficient, or the conventional measure of the
strength of radar signals reflected by a distributed scatterer, usually
expressed in dB. It is a normalised dimensionless number, comparing the
strength observed to that expected from an area of one square meter.”
o (Source: https://sentinel.esa.int/web/sentinel/user-guides/sentinel-1-
sar/definitions)
Menu: Radar > Radiometric > Calibrate (Figure 182)
o In the “I/O Parameters” window
Select the “ML” images as input
Keep default added suffix “Cal” in the output filename
Select your working directory as output directory
o In the “Processing Parameters” window
Check “Output sigma0 band” (= ratio instant/receive backscatter per
unit area in ground range)
o Click on Run
o Repeat the process for the other image
284
9.14. Intermediate steps: visualization (video timing: 8 min
22 sec)
Visualize the 2 sigma0 images in the viewer
o Double-click on each calibrated images in the “Product Explorer” window to
open them in the SNAP viewer (“Sigma0_VV” is the name of the calibrated
images)
o Menu: Window > Tile Evenly to display the images side by side
Figure 183 : Subseted, Multilooked and Radiometrically Calibrated (Sigma0_VV) Sentinel-1 images
displayed side by side
285
o Right click on the band to convert (Sigma0_VV) > Linear to/from dB > Yes to
create a new virtual band
o Repeat that for the other image
Open the 2 converted images in dB by double-clicking on them
Menu: Window > Tile evenly (you can rearrange the windows if you want) (Figure
184)
In dB, there is a much better distinction between land and water pixels
We can see 2 peaks in the dB images histograms:
o 1 small low values peak for pixels over water
o 1 big high values peak for land pixel
Close all images from the viewer but keep the images in the “Product Explorer” tab
Figure 184 : Comparison of the Sigma0_VV Sentinel-1 images before (above) and after (below) the
“Linear to/from dB” image histogram conversion for the archive (left) and flood (right) images
286
9.15. Terrain Correction (video timing: 10 min 22 sec)
Purpose:
o To project the image into a coordinates system, and
o To correct for distortion due to the mountainous terrain
Menu: Radar > Geometric > Terrain correction > Range-Doppler Terrain Correction
> all parameters by default > Run (default suffix = TC) (Figure 185)
o In the “I/O Parameters” window
Select the “ML_cal” images as input
Keep default added suffix “TC” in the output filename
Select your working directory as output directory
o In the “Processing Parameters” window
Keep all default parameters (WGS84 map projection)
o Click on Run
o Repeat the process for the other image
Figure 185 : Parametrization of the Range Doppler Terrain Correction tool of SNAP
287
9.16. Visualization of the terrain corrected images (video
timing: 11 min 32 sec)
Convert the “TC” bands to dB (confer method explained above)
Convert the virtual created dB bands to real file
o Right click on the dB band > convert band
Save to dB band to file
o Select the dB band > Main menu File > Save Product: this saves the dB band
to the image
Visualize the dB terrain corrected bands
o The image has been projected into a coordinate system: the TC image
orientation has changed! (Figure 186)
o In order to be able to compare the non TC an TC images visually, you have to
unsynchronize the views with the dedicated button in the “Navigation” tab
(bottom left panel), and then open each image in a window that you can
display side by side. You will maybe have to unzoom a lot from one of the two
images to adjust the visualization.
o There is no more distortion in the mountain area (Figure 187)
o (You could also do a contrast stretch with the “Color Manipulation” tab
(bottom left panel) to highlight only the pixels over the land)
288
Figure 186 : Comparison of the archive Sentinel-1 images before (left) and after (right) the Terrain
Correction: change of the image orientation due to the projection into a coordinate system
Figure 187 : Comparison of the archive Sentinel-1 images before (left) and after (right) the Terrain
Correction for a mountainous area: the geometric distortion have been corrected
289
9.17. Combination of the Archive and Crisis images to
produce a RGB composite of the flooded area,
preliminary step: images stack (video timing: 13 min
25 sec)
Purpose of the RGB composite: to enable the distinction of permanent water bodies
and flooded areas.
Remarque: The RGB composite is effectively produced in the section 9.19 page 293.
Application of a stack by simply using the products geolocation
o Menu: Radar > Coregistraction > Stack Tools > Crate Stack tabs
o In window “1-ProductSet-Reader”
Click on the “Add opened” button > and keep only the TC file by
removing unwanted files or open only the 2 images of interest with
the Add button
o In window “2-CreateStack”
Choose “Product Geolocation” for “Initial offset method”
([…Remarque: for interferometry, a much more precise geolocalization
method should be used…])
o In window “3-Write”
Remove part of the filename which is not common to both images: the
acquisition dates
o Select Run
o Close window
290
Figure 188 : Parametrization of the « Create Stack » tool of SNAP
291
9.18. Overlay of the stacked images for visual comparison of
the 2 images (optional)
(This step is not really necessary in the general workflow) (video timing: 15 min 33
sec)
o We get the stack image file which is one single image file that contains 4
layers coming from the 2 archive and crisis images in the “Bands” subfolder
o Open the 2 images in dB
o Apply a contrast stretch with the histogram focused on the higher pixel
values peak
o To overlay the 2 images in the same viewer:
Click on the Archive image in the viewer to select it
Go to “Layer manager” (right panel) > + icon >Image of Band / Tie-
point Grid > Next > overlay the 2nd image in dB (Crisis image) on the
first one > Finish
o You can then compare both images by checking unchecking one image in the
right panel “Layer manager” or by playing with the transparency slider at the
bottom of the “Layer manager”
Figure 189 : Overlay of the 2 TC dB images in 1 viewer and use of the « Layer Manager » panel
allowing to check/uncheck one of the 2 overlaid layers as well as applying a transparency in order
to compare the 2 images
292
9.19. RGB composite image (video timing: 16 min 50 sec)
Purpose: to distinguish between flooded area and permanent water bodies
Select the name of the stack file in the “Product Explorer” window
Menu: Window > Open RGB window
o Red band: select Archive image
o Green and Blue bands: select Crisis image
o Click OK
o View the flood map (Figure 190)
Figure 190 : RGB composite using, in the Red band, the Archive image, and in the Green and Blue
bands, the Crisis image
293
9.20. Interpretation of the RGB composite image (video
timing: 17 min 28 sec)
The RGB composite image is presented in Figure 190.
Flooded areas
o Appears in red because, given the selected RGB composite above, where
floods occur,
the Archive image (in the red channel) has higher backscatter (no
flood) than
the Crisis image (in the green and blue channel) (low backscatter for
flooded areas)
o So for flooded areas, there is a high value in the red channel and low value in
the green and blue channels.
Non-flooded land
o Appears in tones of gray as there is approximately the same backscatter for
the 2 images (no difference), and that that “same” information is associated
to all 3 channels
Permanent water bodies
o Appears as uniform dark areas as there is a low backscatter values for both
Archive and Crisis images, both associated to the 3 RGB channels
Some part of the land are cyan
o This translates a higher response in the green and blue channels
corresponding to the crisis image, than in the Archive image. This may be due
to particular ground cover which is not related to flood.
294
9.21. Exportation of the RGB composite of the flood map in
another format and visualization of the result in
GOOGLE EARTH PRO (video timing: 19 min 22 sec)
(Menu: File > Export > Geotiff). This step is not done in the video. This is just an
example.
The visualization of the RGB composite of the flood map in GOOGLE EARTH PRO is
presented in (Figure 191).
Right click on the view > Export View as Google Earth KMZ
o Record the output in your working directory and call it for example “Flood”
o Browse to you working directory where the file “Flood.kmz” has been saved
o If you have Google Earth Pro installed on your computer, you can simply
double-click on this file and it will automatically open in Google Earth at the
right position.
o (If you do not have Google Earth Pro installed on your computer, you can get
it here for example: https://www.google.com/intl/fr_ALL/earth/versions/
(bottom of page))
o You can then compare the flood map with Google Earth Pro imagery and
also check the registration (georeferencing) of the flood map, which is good
in this case (by deselecting the layer or playing with the layer transparency).
In Google Earth Pro you can also apply an oblique view and ask for a 3D
rendering of the natural relief (that can be exaggerated: menu Tools >
Options > 3D View > Relief > and set the elevation factor to 3).
295
Figure 191 : Outlook of the RGB composite flood map overlaid in Google Earth Pro with an oblique
3D view
296
9.22. Autre tutoriel: Satellite Remote Sensing of Flood
Monitoring and Management (NASA ARSET)
Ce tutoriel en anglais présente différentes méthodes d’identification des zones inondées à
parti d’images SAR Sentinel-1 dans le logiciel SNAP et Google Earth Engine. Les méthodes
proposées sont :
Composition RGB multi-temporelle dans SNAP (même méthode que dans la section 9
de ce manuel).
Classification par seuillage simple dans SNAP
Classification supervisée par « random forest » dans SNAP
Classification par seuillage simple dans Google Earth Engine
https://arset.gsfc.nasa.gov/water/workshops/india2018
SAR Applications for Flood Mapping / SAR Application for Mapping the Kerala Floods,
tutoriel SNAP + Sentinel-1, PDF :
https://arset.gsfc.nasa.gov/sites/default/files/disasters/IndiaFlood18/sar_hands_on_
india-final.pdf
Rappel théorique utile sur le SAR :
https://arset.gsfc.nasa.gov/sites/default/files/disasters/IndiaFlood18/D2P2-
Overview-SAR-final.pdf
297
10. SPIRITS : Analyse d’une série temporelle
d’images NDVI basse résolution dans
SPIRITS
SPIRITS (https://mars.jrc.ec.europa.eu/asap/download.php)
ATTENTION ! Ce logiciel est déjà inclus dans le dossier mis à disposition
Logiciel ou « TIME_SERIES_ANALYSIS\SPIRITS\SpiritsExtract_152\ ». Pas besoin de
application le re-télécharger donc ! De plus il n’y a pas de procédure d’installation. Il
doit simplement être disponible dans le répertoire dans lequel il se
trouve actuellement.
298
Temps
~ 4 à 6 heures
approximatif
Prérequis Aucun
299
10.2. Introduction et motivation
Contextualisation :
Pour en savoir plus, jetez un œil rapide au site web suivant (Figure 193):
o https://mars.jrc.ec.europa.eu/asap/
o https://mars.jrc.ec.europa.eu/asap/about.php
o ASAP is an online decision support system for early warning about hotspots
of agricultural production anomaly (crop and rangeland), developed by the
JRC for food security crises prevention and response planning anticipation.
o https://mars.jrc.ec.europa.eu/asap/country.php?cntry=217
o https://mars.jrc.ec.europa.eu/asap/wexplorer/
300
o A explorer en quelques minutes
A travers cet exercice vous mènerez une analyse qui vous permettra notamment de :
Produire des cartes d’anomalie des conditions de végétation sur le Sénégal pour
l’année 2002 : ce type d’analyse, si menée en temps réel lors d’une période de
sécheresse par exemple, pourrait permettre d’objectiver, de quantifier et de
cartographier le mauvais état relatif des conditions de végétation en tout lieu du
Sénégal et par là contribuer au déclanchement d’un éventuel processus d’aide
alimentaire d’urgence et/ou permettre l’activation d’indemnisation des agriculteurs
impactés, via un système d’assurances agricoles.
Produire des images de paramètres phénologiques de la végétation (dates de
début/fin de la saison de végétation, valeur maximum du NDVI, etc). Ce type
d’images permet de caractériser les conditions de végétation au cours d’une saison
et peuvent être utilisées par exemple dans un modèle de prévision des rendements
agricoles.
https://land.copernicus.eu/global/
301
Choisir « Vegetation » : https://land.copernicus.eu/global/themes/vegetation
Jeter un œil à :
Définition du NDVI
Section « NDVI characteristics » pour le NDVI 1km
o Access
o Algorithm
302
o Technical
Arriver sur
https://land.copernicus.vgt.vito.be/PDF/portal/Application.html#Browse;Root=513186;Colle
ction=1000085;Time=NORMAL,NORMAL,-1,,,-1,,
303
S’enregistrer via « Register »
o Bien noter son « user name » et « password »
o Finaliser la procédure d’enregistrement
304
Une fois inscrit, se loguer : bouton « Login » en haut à droite
Définir les critères de recherche d’images :
o Taper « Senegal » + Go dans le menu de recherche géographique pour définir
le ROI du Sénégal
(ou une zone un peu plus large pour avoir une marge autour du
Sénégal, pour un meilleur style des cartes ?)
o Collection : « NDVI 1km V2.2 Global »
o Date : depuis le 01/04/1998 jusque maintenant
Cliquer sur « Search »
Arriver sur la page présentant les images correspondant aux critères de recherche
o 774 produits disponibles le 5/10/2019
o Nombre d’images en détail :
1998 : à partir de 19980401 = 9 * 3 = 27 images
1999-2018 = 20 * 36 = 720 images
2019 : jusque 20190921 = 9*3 = 27 images
TOTAL : 774 images
Remarquez que le nom des fichiers d’images NDVI contient l’information sur la période
correspondante à l’enregistrement des valeurs présentes dans l’image (une décade, 10
305
jours). Cette information est encodée comme suit :
YYYYMMDD,
par exemple 19980401,
où les 4 premiers chiffres correspondent à l’année, les 2 suivants au mois et les 2
derniers à la décade 1, 2 ou 3 du mois sous le format 01, 11, 21 respectivement.
Cette convention reste valide pour l’ensemble de l’exercice, y compris dans le logiciel
SPIRITS.
Cliquer sur « Select all xxx products » pour sélectionner toutes les images (774
produits disponibles le 5/10/2019) (si tel est votre but)
306
Cliquer sur “Order now”
Arriver à
307
Une fois votre commande prête, vous recevrez un email avec les instructions à suivre pour
le téléchargement.
308
Les images téléchargées sont présentes dans le dossier
« C:\TIME_SERIES_ANALYSIS\DOWNLOADED_NDVI_IMAGES\ZIPPED » :
309
Pour pouvoir utiliser ces images, il faut les décompresser ou dézipper. Cela pourra se faire
manuellement pour chaque image, mais cela prendrait beaucoup de temps.
Un fichier de commande DOS a été créé pour dé-zipper toutes les images en un seul clic. Ce
fichier est dans le répertoire « C:\TIME_SERIES_ANALYSIS\DOWNLOADED_NDVI_IMAGES\ »
et s’appelle « UNZIP_NDVI_IMAGES.bat ». Un fichier de commande est un simple fichier
texte mais avec l’extension « .bat » ou « .cmd ». Il permet notamment d’exécuter des
logiciels par ligne de commande.
310
description des commandes (pour information) ici :
https://sevenzip.osdn.jp/chm/cmdline/commands/extract.htm
une fois téléchargé, le décompresser et le placer dans le répertoire
voulu pour obtenir ceci :
« C:\TIME_SERIES_ANALYSIS\DOWNLOADED_NDVI_IMAGES\7zip_co
mmand_line\7za920\7za.exe » comme indiqué dans le fichier de
commande
Une fois votre fichier de commande correctement paramétré, vos images et le logiciel 7-zip
correctement placé dans vos répertoires,
ATTENTION ! 2 images ont une structure de nom de fichier différente des autres, avec la
mention « V2.2.2 » alors que tous les autres sont en « V2.2.1 », confer image ci-dessous.
311
Comme ceci peut provoquer un bug lors de l’importation dans SPIRITS (étape ultérieure), il
faut :
• Renommer ces 2 images en remplaçant « V2.2.2 » par « V2.2.1 » dans les images
copiées.
Vos images décompressées sont prêtes pour être utilisées dans SPIRITS !
https://www.naturalearthdata.com/downloads/
o Choisir « Cultural 10m »)
312
10.5. Aperçu des données dans QGIS (facultatif)
Installer QGIS : confer https://qgis.org/fr/site/
! SPIRITS (logiciel, manuel et tutorial officiels) a déjà été téléchargé et décompressé pour
vous et se trouvent dans le dossier “C:\TIME_SERIES_ANALYSIS\SPIRITS\ ».
https://mars.jrc.ec.europa.eu/asap/download.php
Version 1.5.2, February 2018.
Est disponible également à cette adresse un jeu de données exemple qui est utilisé
dans le tutoriel officiel (non fourni dans le cadre de ce cours étant donné le volume
de ce jeu de données)
“Important note: SPIRITS requires a Java virtual machine (JVM) to run. If not already
installed, users can download Java version 1.7 or higher.” Java peut être téléchargé
gratuitement via le lien ci-dessous. Après téléchargement, installez JVM.
https://www.java.com/fr/download/
313
Double cliquer sur
« C:\TIME_SERIES_ANALYSIS\SPIRITS\SpiritsExtract_152\Spirits.jar » et attendre
quelques secondes
(Remarque : si « Spirits.jar » n’était pas associé au bon programme de lancement par
défaut, et ne s’ouvrait donc pas correctement, il faut associer « spitis.jar » à
« javaw.exe » qui peut probablement être trouvé dans un répertoire ressemblant à
ceci : « C:\Program Files\Java\jre1.8.0_171\bin\javaw.exe ».)
314
Le répertoire à utiliser pour cet exercice dans
SPIRITS est le répertoire :
“C:\TIME_SERIES_ANALYSIS\SPIRITS\SENEGAL\”
315
Abréviation Nom des Fonction
fichiers
TNT Task files Permet d’exécuter un scénario (.SNS) sur une période donnée
SNS Scenario files Contient la paramétrisation nécessaire pour effectuer une action
(ex : importation des fichiers avec certaines valeurs des
paramètres de cette fonction)
UNT User Tool files /
LNL List of files Dossier contenant des fichiers texte contenant des listes de
fichiers…
Liste de quoi… ?
Utile pour… ?
CNC Chart files Dossier contenant des templates de graphiques
QNQ Quick Look files Dossier contenant les « map template » ou mise en page de carte
prédéfinie
DATA Image diles Dossier contenant les images NDVI avec 1 sous-dossier par
traitement (importation, lissage temporel,…)
VAR VAR/MTA files /
RUM RUM files RUM veut dire: “Regional and Unmixed Mean values” et désigne
les valeurs qui sont extraites des images par région, par exemple
316
la valeur moyenne du NDVI par zone d’intérêt (ex : par
département du Sénégal).
SPX Specification Paramétrge aditionel de certaines fonctions
files (*.SPX)
PNG PNG files Dossier contenant des images n’étant pas les images satellites
analysées:
cartes produites
graphique
logo et autres illustrations à faire apparaitre sur une carte
…
TXT TXT files Dossier contenant des fichiers texte. Par exemple les fichiers
utilisés pour définir la structure de la base de données qui
accueillera les statistiques calculées.
VEC Vector files Dossier contenant les fichiers vectoriels (shapefile) utilisés pour la
mise en page cartographique ou fichiers devant être rastérisé
pour l’extraction de statistiques.
RST Spatial Dossier contenant les régions d’intérêts rastérisées (en format
Reference Set raster) car SPIRITS a besoin de ce format raster pour pouvoir faire
files l’extraction.
• Ensuite :
• Soit, créer 2 scénarios comme dans l’image ci-dessous. La partie de l’interface
non montrée (…) peut rester vide.
• Soit, pour aller plus vite ouvrir le scénario pré-enregistré pour vous via la
commande « File > Open > » de la fenêtre « Generic Import Scenario » et
choisir le scénario suivant :
• Noms des scénarios :
• Nom du scenario pour importer les images SPOT VEGETATION:
• “IMPORT_and_CONVERT_VGT_tiff_to_img_ENVI.sns”
317
• Nom du scenario pour importer les images PROBA-V:
• “IMPORT_and_CONVERT_PROBAV_tiff_to_img_ENVI.sns”
• La seule différence entre ces 2 scénarios se situe au niveau du suffixe des
« input files ».
• Par convention, pour l’output :
• Préfixe = cgls (= Copernicus Global Land Service)
• Suffixe = i (=NDVI)
• Une partie des infos de la rubrique « Spectral » peuvent être trouvée sur la
page présentée précédemment :
https://land.copernicus.eu/global/products/ndvi
319
Figure 194 : Principe de la conversion des DN (Digital Number des pixels) en valeurs de NDVI
320
• Après exécution de l’importation, vérifier que l’exécution s’est opérée
correctement :
• Via l’onglet « Results »
• point noir : OK
• point rouge : ERROR
• Via le nombre de fichiers présents dans le dossier dans lequel les images ont
été importées (dossier
« C:\TIME_SERIES_ANALYSIS\SPIRITS\SENEGAL\DATA\SEN\NDVI_1km_V2.2.1
»)
• 774 images * 3 fichiers/image = 2322 fichiers doivent être présents
La série complète d’images NDVI a été importée dans SPIRITS dans le format ENVI. Elle
se trouve dans le répertoire
« C:\TIME_SERIES_ANALYSIS\SPIRITS\SENEGAL\DATA\SEN\NDVI_1km_V2.2.1 ».
321
La description détaillée du « header file » est faite dans le manuel officiel de SPIRITS à la
page 21.
Dans un premier temps, vous allez visualiser vos images via la création en chaine
(automatique) de carte pour chaque image.
322
10.6.5.1. Création d’un « map template » ou template de carte (QNQ)
Pour créer des cartes en chaine, il faut d’abord créer un « map template », c’est-à-dire une
mise en page cartographique. Elle pourra être ensuite utilisée pour l’ensemble des images.
Pour gagner du temps, ouvrez un map template préalablement préparé pour vous :
Soit paramétrez la fenêtre comme vous voulez en réutilisant le map template (QNQ)
créé dans la section précédente et sauvegarder cette tâche comme un nouveau
« task files » (TNT)
Soit, pour gagner du temps, ouvrez le TASK FILES (TNT)
« CREATE_MAPS_SEN_NDVI_2013_2015.tnt » préalablement créé pour vous.
Modifiez éventuellement les dates de début et de fin selon la période pour laquelle
vous voulez produire des cartes (et dans ce cas, sauvegardez votre nouveau TNT (File
> Save as)).
323
Cliquer sur « Execute » pour produire la série de cartes.
Après exécution de la création de cartes en série, vérifier que l’exécution s’est opérée
correctement :
324
325
10.6.6. Extraction de statistiques à partir des images NDVI
Pour « SENSORS » :
326
Pour « VARIABLES » : définir 2 variables :
327
o Les autres variables déjà présentes dans ce fichier texte correspondent à des
variables phénologiques et seront utilisée ultérieurement (confer section
10.6.9 page 359 et suivantes)
Pour « REGIONS » :
Pas possible d’importer de « Regions sets » faire un ajout manuel (confer image ci-
dessous). Ceci correspond à la région générale : on lui donne un nom.
Importer ensuite les « Regions in Regions set » à partir du shapefile des
départements du Sénégal qui est dans le dossier
« C:\TIME_SERIES_ANALYSIS\SPIRITS\SENEGAL\VEC\Senegal_deps.shp »
o Confer le paramétrage dans l’image ci-dessous
o Ceci correspond aux sous-régions dans la région générale : on indique à quoi
correspondent les régions qui seront utilisées pour l’extraction de
statistiques : on leurs donne des noms.
328
o Notez que la colonne du shapefile appelée « SNNIV3_ID », qui varie de 1 à 35
et qui identifie chacun des départements du Sénégal par un identifiant
unique, est utilisée comme identifiant des régions dans la base de données.
Figure 195 : Aperçu du shapefile sur le Sénégal et sa table d’attributs avec 35 départements avec
les identifiants de 1 à 35 dans la colonne « SNNIV3_ID »
329
Figure 196 : Importation des « Regions » dans SPIRITS à partir du shapefile « Senegal_deps.shp »
Pour « Classes » :
330
En résumé, les 4 éléments de la base de données sont créés :
331
332
10.6.6.2. Création de l’image raster des zones d’intérêt pour l’extraction des
statistiques
Pour pouvoir faire l’extraction des statistiques des images NDVI par zones d’intérêt (les
départements du Sénégal), SPIRITS a besoin que ces zones d’intérêt lui soit fournies en
format raster (image, matrice). SPIRITS ne peut utiliser directement un fichier vectoriel
shapefile.
Il faut donc créer une image raster à partir du fichier vectoriel des départements du Sénégal.
Cette image doit être compatible en termes de système de coordonnées, résolution spatiale
et étendue avec les images sur lesquelles l’utiliser (images NDVI dans notre cas).
Pour plus d’information sur cette étape se référer au tutorial officiel de SPIRITS, page 61
« Exercise 4-6 Rasterize Shapefiles ».
Pour ce faire:
333
confer image ci-dessous
Cliquer sur « Execute » pour lancer la rastérisation
Vérifier le résultat :
o Dans le panneau « Results » de SPIRITS (point noir)
o Dans le dossier
« C:\TIME_SERIES_ANALYSIS\SPIRITS\SENEGAL\DATA\SEN\RASTERIZED\
SEN_DEPATEMENTS_RASTERIZED.img » (présence de la couche)
o Dans QGIS ou un autre logiciel : afficher l’image rastérisée (facultatif)
334
Figure 197 : Visualisation du fichier « SEN_DEPATEMENTS_RASTERIZED.img » dans le logiciel ENVI
avec palette de couleurs « EOS A ».
10.6.6.3. Extraction des statistiques sur les valeurs de NDVI à partir des
images
Spécifications de l’extraction de statistiques : choix de la méthode (SPU)
Plusieurs méthodes d’extraction des valeurs des images par zone sont proposées dans
SPIRITS. Un fichier de « spécification » doit être établi afin de spécifier le type de méthode
qui sera utilisé.
« no land use classification », soit une valeur moyenne par zone (par département)
(= regional means ("Overall Mean")) sans tenir compte de l’occupation du sol.
Pour plus d’information sur les autres méthodes, voir les pages 249 et 251 du manuel officiel
de SPIRITS :
« Hard classification method : one mean value per class per region” (“class” = classe
d’occupation du sol)
“Area fractions method: one mean value per class per region, weighted by the AFI of
that class.” (AFI = Area Fraction Images)
335
Analysis > Database > Extract > Specification
Paramétrer le SPU :
Avant de procéder à l’extraction, un scenario d’extraction doit être créé. Celui-ci permet de
préciser d’autres paramètres de l’extraction.
Pour ce faire :
Paramétrer le scénario :
336
ATTENTION ! Remarque importante !
Les ID du « Region Sets » de la base de données (« Region Sets » appelé « DPTS_SEN » dans
ce cas-ci et préalablement construit à partir du shapefile des départements du Sénégal avec
des ID de 1 à 35) doivent correspondre aux ID de l’image utilisée pour l’extraction de
statistiques (l’image rastérisée à partir du shapefile des départements du Sénégal dans ce
cas-ci « SEN_DEPATEMENTS_RASTERIZED.img »). C’est cette correspondance qui permet de
placer les valeurs extraites dans les bonnes « cases » correspondantes de la base de
données.
“In case the extracted data will be uploaded into the projects RUM database, the ID's in this
Regions IMG must agree with the Id's of the Regions contained in the selected Regions Set
from the database: this Id is one of the links between the data extracted by the Extract RUM
tool and the data imported into the RUM database.” (page 246, SPIRITS official manual)
337
Analysis > Database > Extract > Time series
338
10.6.6.4. Visualisation des profils temporels de NDVI par région
Maintenant que les statistiques ont été calculées, extraites, et stockées dans la base de
données du projet SPIRITS, il est possible de la consulter, notamment via la création de
graphiques.
Pour visualiser un profil temporel de l’évolution du NDVI entre 1998 et 2019 pour une
région donnée :
339
Dagana, Podor, Louga : forte variabilité inter-annuelle : Nord du Sénégal, faibles
précipitations !
Remarquer la variabilité du pic de NDVI d’une année à l’autre, en particulier pour 2002 où le
pic est faible.
340
Dakar et Pikine : décroissant. Urbanisation ?
341
342
Source: Analysis of mean climate conditions in Senegal (1971–98), S Fall, D Niyogi, FHM
Semazzi - Earth Interactions, 2006, journals.ametsoc.org,
http://landsurface.itap.purdue.edu/publications-protected/J41.pdf
343
Pour visualiser un profil temporel de l’évolution du NDVI entre 1998 et 2019 pour une
région donnée avec superposition des années:
344
Bruit du signal NDVI et lissage temporel
On peut remarquer que le signal NDVI « brut », celui utilisé pour l’extraction des statistiques,
présente une certaine irrégularité dans son évolution intra-annuelle : il y a un certain bruit
du aux conditions atmosphériques changeantes d’une décade à l’autre.
Pour diminuer l’impact des conditions atmosphériques sur le signal et de la sorte produire
un signal qui traduit mieux les conditions de végétation réelles, il est possible d’effectuer
un « lissage temporel » du signal NDVI. Ceci sera fait dans la section suivante.
345
10.6.7. Lissage temporel des images NDVI et extraction des statistiques
La courbe d’évolution temporelle du NDVI dérivée de chacun des pixels de l’image (courbe
noire « Original » sur le graphique) présente un certain bruit. Cette courbe correspond à
l’évolution du NDVI pour un pixel donné sur la série d’images NDVI utilisée.
Cette courbe peut être lissée via de nombreuses méthodes (courbe rouge « S-G filter » sur le
graphique (algorithme de Savitzky-Golay)).
Les valeurs de NDVI lissées (courbe rouge) sont ensuite utilisées pour produire de nouvelles
images de NDVI lissées temporellement.
Source: Construction of smooth daily remote sensing time series data: a higher spatio
temporal resolution perspective, Zhuokun Pan, Yueming Hu, Bin Cao, 2017, DOI
10.1186/s40965-017-0038-z, https://link.springer.com/content/pdf/10.1186%2Fs40965-
017-0038-z.pdf
Soit paramétrer la tâche de lissage temporel comme indiqué dans l’image ci-dessous
et sauvegarder éventuellement votre « Task file » (TNT)
346
Soit, pour gagner du temps, ouvrir (File > Open) le « Task file » (TNT) préalablement
établi pour vous
« C:\TIME_SERIES_ANALYSIS\SPIRITS\SENEGAL\TNT\SMOOTH_NDVI_TEMPORALY_SE
N_1998_2019.tnt »
o Suffixe k : par convention, qui signifie qu’il s’agit d’un NDVI lissé
temporellement
o Output flags : choisir dans ce cas-ci « Copy UNI-flags » car nos images
respectent ce système. Confer manuel officiel de Spirits page 25 section
« 2.3.2. The UNIflags system ».
347
Cliquer sur « Execute » pour lancer le lissage temporel (temps d’exécution : 1 minute)
Vérifier le résultat :
o Dans le panneau « Results » de SPIRITS (point noir)
o Dans le dossier
« C:\TIME_SERIES_ANALYSIS\SPIRITS\SENEGAL\DATA\SEN\NDVI_1km_V2.2.1
_SMOOTHED\ » (présence des images)
Le même type d’opération a déjà été effectué précédemment sur les images NDVI brutes
dans la section 10.6.6 page 326 et suivantes.
348
Création de la structure de la base de données :
o créer une variable sNDVI (smoothed NDVI) : déjà fait à la section 10.6.6.1.
Création du scenario d’extraction à partir du scénario précédant:
o Analysis > Database > Extract > Scenario
o Soit
File > Open > RUM_Extract_scenario_NDVI_SEN.sns
Modifier le Scénario comme indiqué dans l’image ci-dessous et
sauvegarder éventuellement votre scénario
o Soit, pour gagner du temps,
Ouvrir (File > Open) le « scénario » (SNS) préalablement établi pour
vous «
C:\TIME_SERIES_ANALYSIS\SPIRITS\SENEGAL\SNS\RUM_Extract_scen
ario_NDVI_SMOOTHED_SEN.sns »
o Attention à :
« input path »: smoothed NDVI
« output path »: NDVI_SMOOTHED
« suffix » : k (k = par convention, qui signifie qu’il s’agit d’un NDVI lissé
temporellement)
« Variable » : sNDVI
349
Extraction des statistiques des images NDVI lissées temporellement par zone
d’intérêt
o Analysis > Database > Extract > Time series
o Paramétrer la tâche d’extraction :
Soit paramétrer la tâche d’extraction comme indiqué dans l’image ci-
dessous et sauvegarder éventuellement votre « Task file » (TNT)
Soit, pour gagner du temps, ouvrir (File > Open) le « Task file » (TNT)
préalablement établi pour vous «
C:\TIME_SERIES_ANALYSIS\SPIRITS\SENEGAL\SNS\EXTRACT_RUM_SE
N_NDVI_SMOOTHED_1998_2019.tnt »
o Cliquer sur “Execute” (temps d’exécution: 7 minutes)
o Vérifier le résultat :
Dans le panneau « Results » de SPIRITS (point noir)
Dans le dossier
“C:\TIME_SERIES_ANALYSIS\SPIRITS\SENEGAL\RUM\SEN\NDVI_SMOO
THED\ » (présence des fichiers de résultat)
Une fois les statistiques extraites à partir des images lissées temporellement, elles peuvent
être explorées, via des graphiques notamment, comme similairement fait pour les images
NDVI brutes à la section 10.6.6.4 page 339 mais en choisissant cette fois la variable
« sNDVI » !
Exemple pour Dagana dans l’image ci-dessous (images NDVI lissées temporellement mais
l’option de lissage du graphique n’a pas été utilisée ici).
350
Ces images temporellement lissées sont de meilleure qualité que les images brute dans le
sens où elles expriment mieux le comportement réel de la végétation à travers le temps car
elles sont moins impactées par les variations des conditions atmosphériques.
Dans le cadre de cet exercice nous allons créer une carte d’anomalie du NDVI pour chaque
décade de l’année 2002 pour le Sénégal en comparant le NDVI de chaque décade de l’année
2002 à la valeur moyenne observée sur toute la série d’images disponibles (1998-2019) pour
chaque décade. Cette valeur moyenne sera appelée « Long Term Average » (LTA). Les étapes
à suivre sont :
Calcul du LTA
Calcul des images d’anomalie
Visualisation des images d’anomalie
351
Pour calculer les images NDVI décadaires « Long Term Average » pour le Sénégal pour la
période 1999-2018 (1998 et 2019 ne sont pas considérées vu que ces séries ne sont pas
complètes) :
Vérifier le résultat :
o Dans le panneau « Results » de SPIRITS (point noir)
o Dans le dossier
“C:\TIME_SERIES_ANALYSIS\SPIRITS\SENEGAL\DATA\SEN\NDVI_1km_V2.2.1_
352
SMOOTHED_LTA\ » (présence de 36 images (72 fichiers (img+hdr)), soit 1 par
décade, pour chaque type de résultat demandé (mean, median, etc))
ATTENTION ! Les résultats de l’analyse « Long Term Average » utilisent
un code année allant de 1950 à 1964 qui a une signification précise :
confer illustration ci-dessous.
Les images de la valeur moyenne auront comme code année « 1962 »
Si pour un endroit et une décade donnés, on observe un NDVI de 0.2 et que le LTA
NDVI correspondant est de 0.8, alors on obtient une différence relative de (0.2-
0.8)/0.8 = -0.6/0.8 = -75%, c’est-à-dire que le NDVI est de 75% inférieur à la LTA.
353
Si NDVI = 0.8 et LTA = 0.2, alors la RD = (0.8-0.2)/0.2 = 300%, soit un NDVI supérieur
de 300% à la LTA.
Processing > Temporal > Anomalies > Time series > et création d’un “New”
scenario (bouton “New”)
Ensuite,
o Soit paramétrer le « Difference scenario » comme indiqué dans la fenêtre ci-
dessous, sauvegarder éventuellement le « scenario file » (SNS)
o Soit ouvrir (File > Open) le « scenario file » (SNS) préalablement établi pour
vous disponible dans le répertoire
« C:\TIME_SERIES_ANALYSIS\SPIRITS\SENEGAL\SNS\
RELATIVE_DIFFERENCE_SNDVI_SEN_2002_LTA_1999_2018.sns »
Ensuite :
354
o Dans « Time Series », indiquer l’année (ou les années) pour lesquelles calculer
des images d’anomalie.
Soit ouvrir (File > Open) le “task file” (TNT)
“RELATIVE_DIFFERENCE_SNDVI_SEN_2002_LTA_1999_2018 »
Exécuter la tâche (temps d’exécution : 10 secondes)
Vérifier le résultat :
o Dans le panneau « Results » de SPIRITS (point noir)
o Dans le dossier
“C:\TIME_SERIES_ANALYSIS\SPIRITS\SENEGAL\DATA\SEN\NDVI_1km_V2.2.1_
RELATIVE_DIFFERENCE_2002\ » (présence des 72 fichiers (36 images (1 par
décade)))
Pour gagner du temps, utiliser le map template préalablement créé pour vous et disponible
dans le dossier « C:\TIME_SERIES_ANALYSIS\SPIRITS\SENEGAL\QNQ\
Map_SEN_TEMPLATE_RELATIVE_DIFFERENCE_2002.qnq ». Posez-vous cependant la
question du meilleur choix des couleurs de la légende et des classes de valeurs à utiliser :
Pour le choix des couleurs, s’intéresser à la gamme de variation des valeurs des
images d’anomalie, accessible via le header file de ces images :
o
o Minimum = (0.01*25) - 1.25 = 0.25 - 1.25 = -1, soit - 100%, ce qui correspond
bien à la limite théorique inférieur de la différence relative.
o Maximum = (0.01*250) - 1.25 = 2.5 - 1.25 = 1.25, soit + 125%. !! Cette valeur
maximale constitue une limitation dans la représentation des anomalies
355
positives (confer l’exemple (extrême) chiffré ci-dessus où on arrive à une DR
de +300%). (confer à ce propos le tutoriel officiel de SPIRTIS page 85 et le
manuel de SPIRITS page 170 : il est probable que les valeurs dépassant cette
limite supérieure soient automatiquement converties vers cette limite
supérieure).
Cette carte d’anomalie (ci-dessus) de fin août 2002 (20020821) met bien en évidence les
valeurs de NDVI très inférieure à la moyenne LTA (-50% à -70%), en particulier dans
l’extrême nord du Sénégal (département de Dagana).
Une fois le « map template » réalisé, vous pouvez l’enregistrer et l’utiliser pour produire
toutes les cartes d’anomalie pour toutes les 36 décades de 2002. Pour ce faire :
356
Ensuite:
Soit paramétrer la fenêtre “Create Quick Look” comme dans l’image ci-dessous,
sauvegarder éventuellement le task file et exécuter la tâche
Soit, pour gagner du temps, ouvrez (File > Open) le « task file » (TNT) créé
préalablement pour vous et disponible dans le répertoire
« C:\TIME_SERIES_ANALYSIS\SPIRITS\SENEGAL\TNT\
RELATIVE_DIFFERENCE_SNDVI_SEN_2002_LTA_1999_2018.tnt », et exécutez-le
(temps d’exécution : 30 secondes).
Vérifier le résultat :
o Dans le panneau « Results » de SPIRITS (point noir)
o Dans le dossier “ C:\TIME_SERIES_ANALYSIS\SPIRITS\SENEGAL\PNG\SENEGAL
NDVI RELATIVE DIFFERENCE 2002\ » (présence des 36 cartes pour les 36
décades de 2002) (confer image ci-dessous (cartes rognées ici)).
357
L’ensemble de ces cartes d’anomalie (ci-dessus) pour 2002 mettent bien en évidence les
valeurs de NDVI qui sont en majorité inférieures à la moyenne LTA (couleur rouge), en
particulier dans l’extrême nord du Sénégal en juillet, août et septembre.
Notez la possibilité de créer des cartes d’anomalie sur des périodes d’un mois via
l’agrégation temporelle préalable des images NDVI et le calcul de LTA mensuel. Une carte
mensuelle a l’avantage (par rapport à une carte décadaire) de représenter via une seule
image les conditions de végétation pour un mois complet (il faudrait 3 images décadaires
aussi non).
358
10.6.9. Analyse des variables phénologiques
Les dates de début/de fin de la saison de végétation et les valeurs de NDVI associées
(4 paramètres)
La longueur de la saison (période), l’amplitude du NDVI entre le début et la fin de la
saison, différente surfaces (intégrales) sous la courbe de NDVI, etc
Ces paramètres ont une connotation « phénologique », c’est-à-dire qu’ils peuvent être
utilisés pour caractériser le développement « phénologique » de la végétation.
Ces paramètres peuvent être utilisés dans des modèles qui tentent de caractériser l’état de
la végétation ou encore de prédire les rendements agricoles par exemple.
359
10.6.9.2. Calcul des variables phénologiques à partir des images NDVI
temporellement lissées (résultat = des images correspondant à des
paramètres phénologiques)
Dans SPIRITS, la fonction « Pheno » accessible via le menu
permet de calculer 12 « paramètres phénologiques ». Ils sont décrits dans la figure ci-
dessous. Remarque : le « Pok » n’est pas à proprement parlé un « paramètre
phénologique » !
Figure 198 : Paramètres phénologiques calculables dans SPIRITS avec la fonction « Pheno »
accessible via le menu « Processing >Temporal > Times Statistics ». Source : page 146 du manuel
officiel de SPIRITS.
Attention, malheureusement, cet outil doit être exécuté autant de fois qu’il n’y a d’années
à traiter et le paramétrage adapté à chaque année (confer Figure ci-dessous).
360
Afin que le calcul des paramètres phénologiques soit le moins possible impacté par la longue
saison sèche qui caractérise la plupart des départements du Sénégal, il est possible de définir
une fenêtre temporelle dans laquelle calculer ces paramètres phénologiques. Cette fenêtre
temporelle devrait idéalement être adaptée en fonction des différences de dynamique des
saisons de végétation dans les différents départements, si ces différences sont significatives.
Ces considérations dépassent cependant le cadre de cet exercice de découverte de SPIRITS.
Dans le cadre de cet exercice nous nous contenterons de définir une fenêtre temporelle
identique pour tous les départements, entre avril (0401, décade 10) et décembre (1221,
décade 36) de chaque année, période qui semble d’ailleurs inclure complétement la saison
de végétation dans le centre et le nord du Sénégal sur base des graphiques présentés à la
section 10.6.6.4 page 339.
361
Cliquer sur « Execute » pour lancer le calcul des paramètres phénologiques (temps
d’exécution : 4 secondes par année)
Vérifier le résultat :
o Dans le panneau « Results » de SPIRITS (point noir)
o Dans le dossier
« C:\TIME_SERIES_ANALYSIS\SPIRITS\SENEGAL\DATA\SEN\NDVI_1km_V2.2.1
_PHENO\ » (présence de 528 fichiers = 22 années (1998-2019) * 12
paramètres (= 12 images par saison) * 2 fichiers (header + img).
Remarque : pour 2019, étant donné la série incomplète (qui s’arrête à 20190911), les
paramètres phénologiques ne sont pas comparables à ceux des autres années !
Le résultat de cette analyse correspond donc à 12 images par saison, avec pour une saison,
1 image par paramètre phénologique.
362
10.6.9.3. Visualisation des images correspondant aux paramètres
phénologiques
Une fois les images des paramètres phénologiques produites, il est opportun de les visualiser
pour :
Pour visualiser les images de paramètres phénologiques, des « map template » peuvent être
créé selon la même procédure que détaillée dans la section 10.6.5.1 page 323.
363
date exprimée en pourcentage de 0 à 100 (200*0.5=100). Les valeurs
de DN réellement observées varient, dans le cas de l’exemple ci-
dessous, entre 4 et 196, soit entre 2% et 98%.
Sur cette base une symbologie allant de 0 à 100% par pas de 10% semble
adaptée (confer image ci-dessus). Cependant cette symbologie ne met pas
clairement en évidence la variabilité spatiale de la variable représentée.
On voit dans l’exemple ci-dessous que très peu de pixels sont concernés
par les faibles valeurs. Il convient dès lors d’adapter la symbologie pour
avoir une plus grande variation de couleurs pour les valeurs plus élevées.
Afin d’adapter au mieux la symbologie, il faut se faire une meilleure idée
de la distribution des valeurs des pixels. Cela peut se faire via la création
364
d’un histogramme (graphique de la fréquence de distribution des valeurs
de pixels) via l’outil :
o Analysis > Images > Histogram > Tool
365
Sur lequel on peut voir que il n’y a que ~ 2.5% de pixels (colonne CUM%)
qui présentent des valeurs de DN ≤ 89 (colonne V = valeurs de DN) (valeur
de % de date associée = 44.5, colonne Y)
la symbologie peut être adaptée pour appliquer la plus grande
variabilité de couleurs aux DN compris entre 90 et 190 (valeurs de % de
dates associées = 45-95, colonne Y) : confer image ci-dessous :
l’amélioration est légère !
366
Le rétrécissement de la gamme de valeurs à la gamme 50% à 75%
augmente le contraste mais les pixels hors de cette gamme ne sont plus
colorés (confer image ci-dessous).
367
Remarque : la symbologie à 3 couleurs présente un bug : elle n’attribue
pas la « From colour » à la « From value » et remplace la « From colour »
par une couleur intermédiaire, ce qui annihile la valeur ajoutée de cette
option.
Les 11 « paramètres phénologiques » sont présentés pour l’année 2002 dans la figure ci-
dessous.
368
369
Remarque: les images Aup et Adn sont quasiment binaires:
Aup : valeurs de 90 et 91
Adn : valeurs de 89 et 90
Pour ce faire,
370
10.6.9.5. Extraction des statistiques par zone (RUM) pour les variables
phénologiques
Création de « fichiers liste » ou « list file » des images des paramètres phénologiques
Pour extraire les statistiques à partir des images correspondant aux paramètres
phénologiques, il faut préalablement créer des « fichiers listes » (« list file »).
Chaque fichier liste listera l’ensemble des 22 images (car 22 années étudiées, entre 1998 et
2019) correspondant à un paramètres phénologique donné.
(Confer page 31 et suivantes du manuel officiel de SPIRITS pour plus d’information sur la
création de « list file »).
371
Pour créer ces fichiers liste :
Analysis > Database > Extract > Tool > Onglet List of files > New
Sélectionner les 22 images correspondant à un paramètre phénologique
Sauvegarder la liste (File > Save as) dans le répertoire
« C:\TIME_SERIES_ANALYSIS\SPIRITS\SENEGAL\LNL\ » avec comme nom, par
exemple pour le paramètre Adn, « Adn_1998_2019 ».
Ensuite,
o Soit répéter l’opération pour les autres paramètres,
o Soit faites de simples copier/coller du fichier texte de liste et adaptez-les via
remplacement de texte.
o Soit se contenter des fichiers listes préalablement créés pour vous dans le
répertoire « C:\TIME_SERIES_ANALYSIS\SPIRITS\SENEGAL\LNL\ »
Rester dans le menu « Analysis > Database > Extract > Tool > Onglet List of files”
Ensuite :
o Soit paramétrer la fenêtre comme indiqué dans l’image ci-dessous, exécuter,
et éventuellement sauvegarder le « task file » correspondant dans le dossier
« TNT ». Puis répéter l’opération pour chaque variable phénologique.
372
!! Ne pas mettre de préfixe et suffixe ici !
o Soit ouvrir (File > Open) le « task file » préalablement créé pour vous pour le
paramètre « Adn » et disponible dans le répertoire
« C:\TIME_SERIES_ANALYSIS\SPIRITS\SENEGAL\TNT\ ». L’exécuter. Puis
exécuter l’outil pour les 10 autres paramètres phéno en modifiant à chaque
fois dans l’interface de l’outil, le fichier liste input et la variable (section
« Parameters »)!
Vérifier le résultat :
o Dans le panneau « Results » de SPIRITS (point noir)
o Dans le dossier
“C:\TIME_SERIES_ANALYSIS\SPIRITS\SENEGAL\RUM\SEN\PHENO\ » (présence
des 242 fichiers de résultat (11 paramètres * 22 années = 242)
10.6.9.6. Visualisation des statistiques par zone (RUM) pour les variables
phénologiques
Une fois extraites, les statistiques peuvent être explorées via consultation de la base de
données et la réalisation de graphiques, comme indiqué préalablement à la section 10.6.6.4
page 339.
373
Analysis > Database > Browse
374
On retrouve pour le paramètre Vmx (maximum NDVI value) le même comportement que
celui observé à la section 10.6.6.4 page 339 pour les départements de Dagana (Vmx faible et
irrégulier), Pikine (Vmx décroissant avec le temps) et Ziguinchor (Vmx haut et régulier).
Cet outil permet de détecter le nombre de saisons par année, les dates de début, maximum
et fin de saison, de même que les valeurs de NDVI (par exemple) associées et la longueur
de la saison. Confer le manuel officiel de Spirits, page 148 et suivante.
375
11. BFAST : Breaks For Additive Season and
Trend (dans R et RSTUDIO)
376
Il est fortement préférable que l’apprenant ait déjà quelques notions de
Prérequis programmation, idéalement dans R. Le non initié risque fort de ne pas
s’y retrouver.
11.1. Introduction
Les sections suivantes présentent d’abord le contenu du dossier de l’exercice, le package R
« BFAST » et ses concepts principaux, le tutoriel source sur lequel est basé cet exercice, la
manière d’installer les logiciels R et RSTUDIO, la présentation des images MODIS «
MOD13Q1 » qui seront utilisées pour l’exercice, et l’outil web MODIS « Fixed Sites Subsets
Tool » qui permet d’accéder facilement à certaines données MODIS. Il n’est pas nécessaire
d’explorer en détails tous les liens qui sont donnés dans ces sections.
Ce n’est qu’à la section « 11.8 Utilisation de BFAST dans R et RSTUDIO » page 387 que
l’exercice pratique commence réellement via l’utilisation du code R « R_CODE_BFAST.R ».
Pour bien comprendre ce que le code R fait, les indications des sections précédentes doivent
avoir été parcourues.
377
11.3. Présentation de BFAST (package de R)
« BFAST, Breaks For Additive Season and Trend » permet la décomposition de séries
temporelles en différentes composantes : une tendance (anglais : « trend »), un aspect
saisonnier (anglais : « season »), et reste (anglais : « remainder ») et propose une séries de
méthodes pour détecter et caractériser les changements (dont des cassures (anglais :
« breaks » qui sont des changements brutaux et importants) au sein d’une série temporelle.
Source: http://bfast.r-forge.r-project.org/
Figure 199 : « Fitted seasonal, trend and remainder (i.e. estimated noise) components for a 16-day
MODIS NDVI time series (data series) of a pine plantation”
378
Source : J. Verbesselt et al. / Remote Sensing of Environment 114 (2010) 106–115
Figure ci-dessus: “Simulated 16-day MODIS NDVI time series. The period from 2004 until
mid-2010 (i.e., the time step just before the simulated break), is considered the history
period and the period after the simulated break is the monitoring period (grey background).
The monitoring period contains 6 observations. The result of the monitoring approach is
shown: A stable history period is identified within the history period (i.e., 2007 until mid-
379
2010) and used to model and predict the normal data variation (blue dashed line) to enable
disturbance detection. Here, a disturbance is detected after 3 observations in the
monitoring period (green vertical line). Applications: Deforestation, forest health
monitoring and phenological change detection within time series of spatio-temporal data
sets (satellite images)."
380
11.4. Tutoriel source
L’exercice qui est présenté dans cette section « BFAST » de ce manuel de travaux pratique se
base en très grande majorité sur le tutoriel de Jan Verbesselt disponible ici :
Ce tutoriel ne doit pas être parcouru vu que toutes les informations nécessaires à la
réalisation de l’exercice sont reprises dans ce manuel et dans le code R mis à disposition. Il
pourra cependant peut être s’avérer une ressource utile en cas de problème.
381
11.5. Installation des 2 logiciels R et RSTUDIO
BFAST est un package du logiciel R. Il s’utilise donc dans l’environnement de R et il vous faut
donc installer R sur votre ordinateur.
RSTUDIO est un logiciel qui sert d’interface pour travailler avec R. Il faut également installer
RSTUDIO sur votre ordinateur.
382
11.6. Présentation des images MODIS « MOD13Q1 »
Le produit NDVI de MODIS est repris dans le produit MODIS appelé « MOD13Q1 ». Ce
produit sera utilisé dans ce tutoriel comme source de donnée du NDVI pour construire la
série temporelle. Cette section présente brièvement les images MODIS et ce produit en
particulier.
https://modis.gsfc.nasa.gov/
https://modis.gsfc.nasa.gov/about/
o MODIS (or Moderate Resolution Imaging Spectroradiometer) is a key
instrument aboard the Terra and Aqua satellites.
o Terra MODIS and Aqua MODIS are viewing the entire Earth's surface every 1
to 2 days, acquiring data in 36 spectral bands, or groups of wavelengths
Technical specifications: https://modis.gsfc.nasa.gov/about/specifications.php
o Spatial Resolution: 250 m (bands 1-2, red & NIR), 500 m (bands 3-7), 1000 m
(bands 8-36)
List of MODIS products:
o https://modis.gsfc.nasa.gov/data/dataprod/
MODIS Vegetation Index Products (NDVI and EVI): MOD13!
o https://modis.gsfc.nasa.gov/data/dataprod/mod13.php
MOD13Q1 v006
o https://lpdaac.usgs.gov/products/mod13q1v006/
o MODIS/Terra Vegetation Indices 16-Day L3 Global 250 m
o NDVI + EVI (EVI: improved sensitivity over high biomass regions)
o 2 vegetation layers + 1 quality layer and 1 quality reliability layer for the
indices NDVI and EVI + 4 reflectance bands (1 (red), 2 (near-infrared), 3 (blue),
and 7 (mid-infrared)), as well as 4 other layers 12 bands
383
Source de l’illustration : https://sites.google.com/site/chempendix/em-spectrum
384
11.7. Présentation du MODIS « Fixed Sites Subsets Tool »
Le MODIS « Fixed Sites Subsets Tool » est une plateforme web qui donne un accès facile à
certaines données MODIS pour des sites prédéfinis. Cet exercice utilisera des données
MODIS provenant de ces sites prédéfinis.
Site web :
https://modis.ornl.gov/sites/
Dans le tutoriel web de Jan VERBESSELT, c’est le site 'Gelderland Loobos' aux Pays-
Bas qui est utilisé.
385
Dans ce manuel, c’est le site « ke_laikipia_mpala » dans le centre du Kenya qui sera
utilisé. Trouvez et sélectionnez ce site via l’interface spatiale ou le menu de
recherche.
Mais libre à vous d’utiliser un autre site. Choisissez cependant un site encore actif
(certains ne le sont plus). Idéalement trouvez un site présentant une évolution
temporelle du NDVI non régulière (présence d’une variation saisonnière, d’une ou
plusieurs tendances et cassures) afin que l’analyse de cette série avec BFAST (section
suivante) produise des résultats « intéressants ». ATTENTION, il vous faudra alors
adapter le nom du site lors de l’utilisation de R (confer section suivante).
Remarquez également :
“Under 'Download data' tab and table you can see the Pixel Numbering Scheme
where the blue pixel (i.e. 545) is the "site pixel" (i.e. the middle of the spatial
subset).”
ATTENTION ! Dans le cadre de cet exercice, la série temporelle d’images NDVI MODIS ne
sera pas téléchargée manuellement à partir de ce site web « Fixed Sites Subsets Tool » de
MODIS MAIS plutôt via l’utilisation du code R mis à disposition (confer section suivante) qui,
lui, téléchargera automatiquement les données à partir de ce site web. Certaines des
informations identifiées ci-dessus (« Site ID », etc) sont reprises dans le code R.
386
11.8. Utilisation de BFAST dans R et RSTUDIO
Dans cette section, vous allez utiliser l’outil BFAST pour réaliser une analyse de série
temporelle d’images NDVI MODIS dans R et RSTUDIO pour un site MODIS prédéfinis (confer
ci-dessus).
Suivez les indications reprises sous forme de commentaires (#) dans ce fichier. Tentez de
comprendre comment fonctionne ce code R en l’exécutant morceau par morceau et
visualisez les résultats intermédiaires et finaux. Faites éventuellement varier certains
paramètres afin de comprendre leur impact sur les résultats de l’analyse.
Quelques résultats sont repris ci-dessous pour le pixel central (N° 545) du site
« ke_laikipia_mpala » dans le centre du Kenya.
387
Résultats de la fonction bfastmonitor() en fonction de différents paramétrages
388
12. Autres outils d’analyses de séries
temporelles d’images satellites
Section donnée à titre d’information
Example: https://lpdaacsvc.cr.usgs.gov/appeears/help
389
TIMESAT software
o http://web.nateko.lu.se/timesat/timesat.asp
PolyTrend
o Web-app : http://polytrend.gis.lu.se/ web application for time series analysis
o R-package: https://cran.r-project.org/web/packages/PolyTrend/index.html
DBEST R package
o https://cran.r-project.org/web/packages/DBEST/index.html
390
13. Introduction à GOOGLE EARTH ENGINE
(GEE)
Logiciel ou
Google Earth Engine (GEE) (https://earthengine.google.com/)
application
Temps
~ 3-4 heures
approximatif
391
13.1. Préalable !!
ATTENTION! Un délai de quelques jours peut être nécessaire pour obtenir un accès à GEE.
Faites donc la demande d’accès plusieurs jours (1 semaine) avant la séance de TP (une
demande faite le 21/11/2019 a pris une dizaine d’heures pour être acceptée).
To get access to Earth Engine, please fill out the form at signup.earthengine.google.com. You
will receive an email titled "Welcome to Google Earth Engine" with instructions for getting
started. » (FAQ GEE).
Pas besoin d’installer de logiciel : les traitements se font dans le cloud GEE.
https://earthengine.google.com/
392
13.2.2. Vidéo d’introduction (2 minutes)
13.2.3. Timelapse
GEE Timelapse est une « vidéo mondiale zoomable » qui permet de visualiser l’évolution de
l’occupation du sol de la planète pendant ces 35 dernières années (1984 - 2018)
o https://earthengine.google.com/timelapse/
393
13.2.4. Exemples de cas d’étude avec GEE
https://earthengine.google.com/case_studies/, notamment
o Global Forest Cover Change
o Global Forest Watch
o Global Surface Water
o …
https://developers.google.com/earth-engine/datasets/
394
13.2.6. 2 Interfaces: EXPLORER et CODE EDITOR
GEE propose 2 interfaces, présentées ici https://earthengine.google.com/platform/
13.2.6.1. EXPLORER
o https://explorer.earthengine.google.com/#workspace
o “The Explorer is a lightweight web app for exploring our data catalog and
running simple analyses.”
https://code.earthengine.google.com/
“The Code Editor is a web-based IDE (Integrated Development Environment) for
writing and running scripts.”
Documentation sur le code editor:
o https://developers.google.com/earth-engine/playground
395
Découvrez les différentes parties de l’interface du CODE EDITOR grâce au “Feature tour”.
396
13.2.7. Documentations
La documentation sur GEE est disponible via :
https://developers.google.com/earth-engine/
o GUIDE: intro, info générale, tutoriels, etc
o REFERENCE : définition précise des algorithmes disponibles
o TUTORIALS : tutoriels
o EDU : autres tutoriels
o DATA CATALOG : catalogue de données
397
13.3. Exercices
o https://drive.google.com/file/d/1NRREEjp_KFaJcmQtqYcof2rGvHcUBBZz/view
o aussi disponible en PDF dans le dossier de formation « GEE_0 GEE Explorer and
Supervised Classification.pdf »
Suivez ce tutoriel
o https://explorer.earthengine.google.com/#workspace
Adaptation du tutoriel :
o Zoomez sur un pays/région qui vous intéresse, par exemple la région de Kasese en
Ouganda (confer exercice de classification SCP dans QGIS)
o Choisissez dans le « data catalog » les 2 jeux de données suivants :
o « Landsat 7 Collection 1 Tier 1 32-Day TOA Reflectance Composite »
o « Landsat 8 Collection 1 Tier 1 32-Day TOA Reflectance Composite »
o Et constatez :
o La présence importante de nuage d’où l’utilité de travailler avec des
composites désennuagés (comme ce sera le cas plus loin dans l’exercice).
o Le problème des stries de Landsat 7 (https://www.usgs.gov/faqs/what-
landsat-7-etm-slc-data?qt-news_science_products=0#qt-
news_science_products)
o Comparez la disponibilité des images pour :
o La période de 32 jours en cours de réalisation (une partie du territoire
seulement est couvert : variable en fonction de la date actuelle par rapport à
la période de 32 jours)
o Une période de 32 jours précédente (100% du territoire est couvert)
o Informations détaillées sur les bandes de Landsat-7 (utiles pour la réalisation des
compositions fausses et vraies couleurs) :
o https://developers.google.com/earth-
engine/datasets/catalog/LANDSAT_LE07_C01_T1_TOA
o Attention ! Le téléchargement du produit finit ne semble pas fonctionner !
L’interface reste à l’étape de « Preparing download… » indéfiniment…
398
Pour en savoir plus sur ce qu’est le « Landsat TOA Percentile Composite » :
Définitions :
Tutoriel de référence :
o https://drive.google.com/file/d/1oD1jbWfjxXFWFXbBs87QoXi1ReO_cvLZ/view
o page 19
o Aussi disponible en PDF dans le dossier de formation sous le nom de :
o « Google Earth Outreach Introductory Tutorial.pdf »
399
13.3.1.3. Raster de classification (2 minutes)
Objectif :
Tutoriel de référence :
o https://drive.google.com/file/d/1oD1jbWfjxXFWFXbBs87QoXi1ReO_cvLZ/view
o page 14
o Aussi disponible en PDF dans le dossier de formation sous le nom de :
o « Google Earth Outreach Introductory Tutorial.pdf »
13.3.2.1. Découvrir et visualiser une collection via le code editor (45 minutes)
Les collections de données sont disponibles ici :
https://developers.google.com/earth-engine/datasets
400
Figure 202 : Aperçu de la collection Landsat-8 Surface reflectance
401
« DOC » et descendez jusqu’à trouver celle qui se rapporte à
« ee.ImageCollection »
402
Code Signification et commentaires
dataset.select('NDVI') Cette commande semble inutile car dans ce cas-ci la
collection ne contient qu’une seule bande, celle du NDVI.
Pour le vérifier, modifiez cette ligne de commande comme
ceci :
o var colorized = dataset;//.select('NDVI');
et exécutez le script. Y a-t-il une difference?
Le signe « // » permet de mettre en commentaire ce qui se
trouve à sa droite , commentaire qui ne sera pas pris en
considération au moment de l’exécution du code.
Multi-line comments can start with /* and end with */
Remarque : il ne semble pas possible d’exécuter une partie
de code seulement (une ligne par exemple). Pour ne pas
exécuter certaines lignes d’un code, il faut les mettre sous
forme de commentaires.
palette: [ Pour la palette de couleur, les « Hex color codes » sont
'FFFFFF', 'CE7E45', 'DF923D', 'F1B555',
'FCD163', '99B718', '74A901', utilisés. Ces codes sont documentés ici par exemple :
'66A000', '529400', '3E8601', '207401',
'056201', '004C00', '023B01',
https://htmlcolorcodes.com/
'012E01', '011D01', '011301'
]
403
13.3.2.2. Tutoriel officiel « The Earth Engine API »
Suivez le tutoriel officiel « The Earth Engine API » à partir de la section « Visualizing Images
and Image Bands » disponible ici :
o https://developers.google.com/earth-engine/tutorial_api_02
Copiez-Collez les lignes de codes dans le CODE EXPLORER, lisez les instructions et
descriptions et faites l’effort de comprendre ce que vous faites. N’hésitez pas à :
o utiliser l’onglet REFERENCE de ce site ou l’onglet DOC du CODE EDITOR pour trouver
plus d’information sur les fonctions utilisées
o Consulter le DATA CATALOG pour obtenir plus d’information sur les données
utilisées.
Objectifs : afficher correctement une image à bande unique, SRTM dans ce cas-ci, via le
contrôle des paramètres de visualisation.
404
13.3.2.2.2. Computations using Images (20 min)
Slope Sin(Aspect)
405
13.3.2.2.3. Image Collections (20 minutes)
Filtres spatial et temporel sur une collection Composition RGB par défaut, de l’image
Landsat-8 dans la « Console » et (entre présentant le moins de nuages pour la zone
autres) metadata « CLOUD_COVER » d’une et la période d’intérêt identifiée par les 2
image filtres
406
Collection Landsat-8, de l’image…(idem) Graphique des valeurs de la bande BQA
pour un pixel pour les 38 images de la
collection Landsat-8 2016 (2720 = pas de
nuage*)
* La signification des valeurs de la bande BQA de Landsat-8 est reprise sur la page web
suivante :
https://www.usgs.gov/land-resources/nli/landsat/landsat-collection-1-level-1-
quality-assessment-band?qt-science_support_page_related_con=0#qt-
science_support_page_related_con
dans la table « Landsat 8 OLI/ OLI-TIRS Level-1 Possible Attributes, Pixel Values, and
Pixel Value Interpretations »
407
Composite « pixel le plus récent » de la collection d’images Landsat-8 2016
Composite « pixel median » de la collection d’images Landsat-8 2016 masquée par les pixels
classifiés comme « water » et « no data » (according to the Hansen et al. dataset) (pixels
masqués = pixels transparents : on voit Google maps à la place du composite median).
408
Mosaïque de :
composite « pixel median » de la collection d’images Landsat-8 2016 masquée par
les pixels classifiés comme « water » et « no data » (according to the Hansen et al.
dataset)
une image des zones classifiés comme « water » et « no data » (according to the
Hansen et al. dataset) colorées en bleu foncé
https://developers.google.cn/earth-engine/landsat
409
14. Visionnage de la vidéo « HOLOGLOBE »
Cette vidéo est réalisée en grande partie à partir de données issues de satellites
d’observation de la terre. Entre autres, les phénomènes suivants sont illustrés:
Cette vidéo, et une série d’informations complémentaires, sont disponibles sur le site web
http://svs.gsfc.nasa.gov/vis/a000000/a000100/a000155/.
410
14.2. Traduction des commentaires de la vidéo
HOLOGLOBE_v3
Source : http://svs.gsfc.nasa.gov/vis/a000000/a000100/a000155/
Durée : 6 minutes 19 secondes
Nous voyons maintenant les océans, les surfaces terrestres, l’atmosphère, la vie, le tout
interconnecté dans un système mondial. La terre est baignée dans l’énergie solaire. Les
océans, surfaces terrestres et atmosphères absorbent et sont réchauffés par l’énergie
solaire. La chaleur absorbée par les océans et transportée par les courants marins est
continuellement libérée vers l’atmosphère. La chaleur et l’humidité libérées par les océans
dirigent la circulation atmosphérique et le climat. L’humidité dans l’atmosphère forme les
nuages qui couvrent en moyenne 40% de la surface terrestre à tout moment donné.
“The pattern in cloud motion in this time lapse sequence shows how earth's winds, moving
bands or zones, which define regional winds directions” (difficile à comprendre et traduire).
411
au large des côtes africaines (2m32) et se déplacent à l’ouest à travers l’Océan Atlantique
vers les Caraïbes et la côte est des USA.
Le climat mondial est également influencé par la manière dont la chaleur est répartie par les
océans. Sont montés ici les changements de températures de surface des océans sur une
période de 9 ans. Les surfaces d’eau les plus chaudes sont indiquées en rouge, les plus
froides en bleu. Avec ce genre de données, il est possible de détecter l’augmentation
anormale de température qui arrive lorsque les eaux froides en temps normal de l’est de
l’Océan Pacifique sont remplacées par des eaux plus chaudes. Ce phénomène est connu sous
le nom d’El Nino. El Nino pourrait perturber le climat à une échelle mondiale causant de
vastes inondations et sécheresses.
Drainant les océans […], il y a la surface solide de la terre avec la croute terrestre qui est
divisée en « Hautes terres » et « Basses terres »
Les Hautes terres sont les terres émergées qui forment les continents
Les Basses terres forment les bassins océaniques.
La croute terrestre n’est pas une coquille fixe et continue. Elle est brisée en une mosaïque
de plaques mouvantes. Au fur et à mesure que ces plaques bougent et se frottent les unes
aux autres, elles libèrent d’énormes quantités d’énergie sous la forme de tremblements de
terre. Les points jaunes montrent la localisation des tremblements de terre qui ont eu lieu
entre 1980 et 1995 avec une magnitude supérieure à 4.5 sur l’échelle de Richter. Ces
tremblements de terre indiquent clairement les limites des plaques de la croute terrestre.
Lorsque les plaques océaniques s’entrechoquent avec les plaques terrestres […] cela forme
les volcans. Les triangles rouges indiquent les éruptions volcaniques enregistrées qui se sont
passées entre 1960 et 1995. Comme pour les tremblements de terre, la plupart des volcans
sont localisés le long des limites des plaques terrestres.
Les océans, l’atmosphère et les continents jouent tous un rôle critique pour le maintient de
la vie durable sur terre. Les caractéristiques des courants marins, des vents, et de la
topographie peuvent faire d’une région un désert et d’une autre une région verte […]. Les
changements de saisons affectent grandement la végétation terrestre. Les couleurs vertes
et jaunes montrent comment les plantes grandissent et dépérissent au cours des saisons
d’une année. Les vastes déserts du nord de l’Afrique, de la péninsule arabique et de l’Asie
centrale sont facilement reconnaissables à cause de l’absence de végétation.
L’influence des hommes sur la terre peut également être vue depuis l’espace. Une vision de
la terre pendant la nuit permet de voir les lumières des campements, villages et villes qui
sont éclairés. Les endroits fortement peuplés apparaissent en blanc. Les incendies et
l’agriculture sur brulis (slash and burn) et la combustion de gaz dans les champs de pétroles
sont montrés en rouge et jaune. Depuis ce point de vue, l’importance de l’activité humaine
sur la planète peut être évaluée. L’entièreté du globe est concernée.
Un suivi global de la terre nous aide à comprendre les interactions complexes entre les
atmosphères, les océans et les terres ainsi que leurs impacts sur la vie. Seule une vue
mondiale nous permet d’apprécier notre terre en tant que planète.
412
15. Sites web intéressants
Quelques exemples de détection de changement avec des images Landsat, vidéo, anglais.
« A Planetary Perspective: With Landsat and Google Earth Engine ».
http://www.youtube.com/watch?v=Ezn1ne2Fj6Y
Quelques exemples d’applications environnementales à partir des images Landsat, site web,
anglais.
« Environmental Watch with Landsat satellites ».
http://www.nasa.gov/content/goddard/environmental-watch-with-nasa-s-landsat-satellites/
15.2. Divers
Programme belge “STEREO”
Belgium has its very own national program supporting research in Earth observation
called STEREO (Support To Exploitation and Research in Earth Observation).
Documentation sur tous les projets STEREO disponible ici:
https://eo.belspo.be/en/stereo-in-action
http://www.youtube.com/watch?v=NNGDj9IeAuI&feature=youtu.be&hd=1.
Site web du film : http://www.homethemovie.org/. « Ce film développe le lien qui
unit l'homme à la Terre. Conçu comme un carnet de voyages, il est constitué
uniquement d'images aériennes et d'une voix off » (Source : Wikipedia).
413
Simulation 3D de la couverture nuageuse mondiale durant 7 jours sur base de données
réelles, vidéo.
https://www.youtube.com/watch?annotation_id=annotation_397244329&feature=i
v&src_vid=JZXErLns1mM&v=zlqjz9OEhk0
414
16. Exemples d’applications de la
télédétection spatiale pour la gestion des
risques et des catastrophes
http://www.esa.int/spaceinvidéos/Vidéos/2016/02/Sentinels_for_Copernicus
Satellites and examples of application
Les données satellites, aériennes, et stations au sol, intégrées dans un même
système : COPERNICUS
Cœur du système = famille de satellites Sentinelle
The role of the different 6 sentinel satellites
https://emergency.copernicus.eu
Copernicus Emergency Management Service (Copernicus EMS) provides information
for emergency response in relation to different types of disasters, including
meteorological hazards, geophysical hazards, deliberate and accidental man-made
disasters and other humanitarian disasters as well as prevention, preparedness,
response and recovery activities. The Copernicus EMS is composed of:
o an on-demand mapping component providing rapid maps for emergency
response and
o risk & recovery maps for prevention and planning and of the early warning and
monitoring component which includes systems for floods, droughts and forest
fires:
415
o Rapid Mapping consists of the provision of geospatial information within
hours or days from the activation in support of emergency management
activities immediately following a disaster. Standardised mapping
products are provided:
to ascertain the situation before the event (reference product),
to roughly identify and assess the most affected locations (first
estimate product),
assess the geographical extent of the event (fdelineation product)
or
to evaluate the intensity and scope of the damage resulting from
the event (grading product).
o Examples here:
https://emergency.copernicus.eu/mapping/list-of-
components/EMSR388
o Map of active risk:
https://emergency.copernicus.eu/mapping/map-of-activations-
risk-and-recovery
FLOODS: The European and Global Flood Awareness Systems (EFAS & GloFAS)
o https://emergency.copernicus.eu/
o Short introductory vidéo explaining the structure of EFAS and the main
objectives using floods in Eastern Europe in May 2010
o Passer le début: combine précipitations + conditions hydrologiques et fait
un warning !
FOREST FIRE: The European Forest Fire Information System (EFFIS)
o https://effis.jrc.ec.europa.eu/
o Visualisation des risques d’incendies:
o https://effis.jrc.ec.europa.eu/static/effis_current_situation/public/index.h
tml
The EMS Drought Observatory (DO)
o https://emergency.copernicus.eu/
Plusieurs exemples de l’utilisation des satellites dans le cadre de la gestion des crises :
https://www.esa.int/Applications/Observing_the_Earth/Copernicus/Emergency_man
agement
416
Services for emergency management response will help mitigating the effects of
natural and manmade disasters such as floods, forest fires and earthquakes and
contribute to humanitarian aid exercises.
http://sertit.unistra.fr/
o Vidéo, 2 min, français
http://fr.euronews.com/2016/06/03/copernicus-ou-quand-l-imagerie-satellite-
permet-d-aider-a-la-gestion-des
o « Copernicus ou quand l'imagerie satellite permet d'aider à la gestion des
catastrophes naturelles, 8 juin 2016 »
o Le SERTIT, impliqué dans le service européen de gestion des situations d’urgence
Copernicus, explique en vidéo sur la chaîne Euronews comment l’imagerie
satellite permet de fournir une aide à la gestion des catastrophes naturelles,
notamment dans le contexte d’inondations en France, Allemagne et en Belgique.
417
16.2. Exemples supplémentaires par outils ou thématique
Google Earth Engine (GEE)
o https://www.youtube.com/watch?v=kIYHGkSb-fU
Cas d’études
o https://earthengine.google.com/case_studies/
o Global Forest Cover Change
o Malaria Risk Mapping
Certains sites web mettent à disposition des fichiers dont l’extension est « .KMZ » ou
« .KML », formats qui permettent la visualisation de ces données dans GOOGLE EARTH
facilement et rapidement sans devoir télécharger l’image brute. Voici quelques exemples :
Inondation
o Australian city of Rockhampton, 2011
o Capteur : Advanced Spaceborne Thermal Emission and Reflection Radiometer
(ASTER) on NASA’s Terra
o http://earthobservatory.nasa.gov/IOTD/view.php?id=48456&src=eorss-iotd
o Cliquer sur “Google Earth file (KML)”
Typhon
o http://www.gearthblog.com/blog/archives/2016/10/acquisition-imagery-natural-
disasters-improving.html
o Hurricane Matthews, 2016, widespread devastation, crossing over Haiti, the
Bahamas and then up the east coast of the United States.
418
o Cliquer sur: “You can see the satellite imagery in Google Earth with this KML file
from Google Crisis Response”
o Dans Google Earth : cocher uniquement : « Hurrican Matthews 2016 > Haiti >
Digital Globa imagery > Les Irois, Haiti »
Tremblement de terre
o http://www.gearthblog.com/blog/archives/2016/09/post-earthquake-
kumamoto-google-earth-3d.html
o Large earthquakes in the City of Kumamoto, Japan, 2016
o Images acquises après le séisme: visualisation des destructions
o Google Earth tour of the area showing all the light blue roofs, which you can view
in Google Earth with this KML file or see in the YouTube vidéo below.
Déforestation
Déforestation au Brésil, Etat du Rondonia, en détails, vidéo haute qualité (30 ‘’)
o https://www.youtube.com/watch?v=JsIB81sLe2w
Autre courtes vidéos sur différents endroits déforestés entre 2000 et 2012, avec
LANDSAT
o https://www.youtube.com/channel/UCGnX8ABNZYuZrhzT-sJVaeg
Monitoring Forests From the Ground to the Cloud , vidéo 5 min, Système d’alerte de
déforestation via Google Earth Engine, exemples internationaux,
o https://www.youtube.com/watch?v=ymKHb3WJLz4
419
Mine
Appalachian Mountaintop Removal in Google Earth & Maps (vidéo en anglais, 5 min)
o https://www.youtube.com/watch?v=aiSzOiGFa-0
Visible également dans GOOGLE EARTH PRO Desktop > Menu sensibilisation
mondiale > extraction minière à ciel ouvert dans les Appalaches
Masse d’eau
Lac Chad: UNEP & Google Earth highlights environmental change, vidéo 1 minute
o https://www.youtube.com/watch?v=JXW29zsr6xg
Inondation
Annual Arctic Sea Ice Minimum 1979-2015 with Area Graph, vidéo
o http://svs.gsfc.nasa.gov/cgi-bin/details.cgi?aid=4435&button=popular
Weekly Animation of Arctic Sea Ice Age with Graph of Ice Age by Percent of Total:
1984 – 2016, vidéo + textes:
o http://svs.gsfc.nasa.gov/cgi-bin/details.cgi?aid=4509&button=popular
420
Feux
https://earthdata.nasa.gov/earth-observation-data/near-real-time/firms
Température du globe
Urbanisation
Tsunami
Littoral
Marée noire
Volcan
Sécurité alimentaire
http://www.fews.net/fr
Avertissements précoces
421
Suivi des conditions de végétation, précipitations, climat, niveau des étendues d’eau,
etc dans l’optique de déclencher des avertissements vers les autorités d’un pays,
d’une région en cas de problème/crise
o https://earlywarning.usgs.gov/fews
Refugiés
Story Map about Rohingya refugees (by the The UN Refugee Agency (UNHCR))
o https://www.esri.com/en-us/arcgis/products/esri-story-maps/contest/winners-
gallery/2018-winners
Divers
ENVI, Les images satellites pour la gestion des catastrophes naturelles, vidéo
français,
o https://www.youtube.com/watch?v=Tp2hIeeSVXE
MapGive
https://mapgive.state.gov/
o Volunteers like you can help trace roads, buildings, and houses using
OpenStreetMap that will aid in humanitarian missions.
422
17. Recherche d’images satellites sur le web
Une des premières étapes de tout travail en télédétection spatiale est de se procurer les
images satellites à partir desquelles votre étude pourra être menée. Cette section à pour
objectif de vous introduire à la recherche d’images sur le web.
Avant de commencer votre recherche d’images, il est important que vous cerniez clairement
quel sont vos objectifs (en termes de résultats) et quels sont vos moyens (connaissance,
argent, temps).
Au cours de vos recherches, vous devrez trouver un juste compromis entre ce que le web
vous « offre » et ce que vous recherchez. En particulier, vous ferez très attention aux critères
suivants:
423
o https://asterweb.jpl.nasa.gov/
o https://asterweb.jpl.nasa.gov/characteristics.asp
o Données d’ALTITUDE : DEM (Digital Elevation Model) ou MNT (Modèle
Numérique de Terrain) pour l’ensemble de la planète à ~30 mètres de
résolution spatiale via l menu « Data Sets > Digital Elevation > SRTM »
o Toute une série de cartes thématiques (occupation du sol,…)
o RADAR
o …
424
les données dans un format plus facile à utiliser (ENVI, géotiff, ERMapper ...) Logiciel
« Crop VGT » disponible pour les pays non africains.
425
AUTRES (gratuit):
Africa data dissemination service
o http://earlywarning.usgs.gov/adds/imgbulks3.php?imgtype=nd&spextent=a
Site de BELSPO: liste de satellites et liens
o http://eo.belspo.be/Directory/Satellites.aspx
o exemples : CORONA, IKONOS2, NOAA, Airborne platform-Apex, MSG3,
PROBA-CHRIS, ENMAP
LISTE DE SITES WEB listant des sites web fournissant des images satellites
15 Free Satellite Imagery Data Sources
o https://gisgeography.com/free-satellite-imagery-data-list/
50 Satellites You Need To Know: Earth Satellite List
o https://gisgeography.com/earth-satellite-list/
Resources for Finding and Using Satellite Images
o https://gijn.org/resources-for-finding-and-using-satellite-images/
426
18. Ressources complémentaires en ligne
Sont reprises ci-dessous une sélection de ressources en ligne pour l’apprentissage des
techniques d’analyses des données d’observation de la terre
Gamme très large et très variée de tutoriaux prêts à l’emploi pour apprendre les techniques
d’analyses des données d’observation de la terre.
Up to date.
https://arset.gsfc.nasa.gov/
http://step.esa.int/main/doc/tutorials/
Le site web présente cependant quelques ratés pour le moment (octobre 2019).
Page d’accueil
https://eo-college.org/
Introduction au cours
https://eo-college.org/courses/echoes-in-space/lessons/history/topic/what-is-radar/
https://www.youtube.com/channel/UCka-UP2X3y3JpJRRvARzEhQ
427
Liens vers une séries de cours RADAR supplémentaires
https://eo-college.org/resources/#results
https://earth.esa.int/web/guest/missions/esa-operational-eo-
missions/ers/instruments/sar/applications/radar-courses
https://www.servirglobal.net/
“SERVIR aims to strengthen the ability of governments and other development stakeholders
to incorporate Earth observations and geospatial technology into decision-making through
workshops, training events, and other capacity building activities.”
https://www.servirglobal.net/training-
materials/search?country=NaN&themes=&pageIndex=1&pageSize=6
18.5. EO LEARN
En anglais.
Série de tutoriaux pour apprendre les techniques d’analyses des données d’observation de la
terre. Cependant certaines leçons semblent dater, elles se focalisent surtout sur l’utilisation
du logiciel BILKO et elles se limitent parfois à retaper un code préétabli.
https://www.learn-eo.org/
428
www.bilko.org
The aim of this guide is to help non-experts in providing a starting point in the decision
process for selecting an appropriate Earth Observation (EO) solution.
https://fromgistors.blogspot.com/
429
19. Newsletters auxquelles s’abonner
Pour rester informé des dernières nouveautés dans le domaine de la télédétection, vous
pouvez vous inscrire aux newsletters suivantes :
430
20. Quelques logiciels de télédétection
spatiale
Nom Prix Site web
Télédétection
SNAP Gratuit https://step.esa.int/main/toolboxes/snap/
Logiciel de l’ESA pour les images « Sentinel » et +
ENVI Payant http://www.exelisvis.com/ProductsServices/ENVI/ENVI.aspx
IDRISI Payant http://www.clarklabs.org/products/idrisi.cfm
ERDAS Payant http://geospatial.intergraph.com/products/ERDAS-
IMAGINE/Details.aspx
QGIS Gratuit http://www.qgis.org/fr/site/
Via des extensions de Quantum GIS : GDAL/OGR, GRASS, SAGA, etc.
R Gratuit http://www.r-project.org/
Via des packages dédiés aux rasters : raster, rgdal, maptools, etc
431
21. Liste de logiciels utilisés dans ce manuel
QGIS 3.10.0 ou version ultérieure
o https://qgis.org/en/site/
o Dont les plugins/extensions suivant :
« Semi-Automatic Classification Plugin » (SCP)
SPIRITS + JAVA:
o https://mars.jrc.ec.europa.eu/asap/download.php#box5
o SPIRITS ne doit pas être installé en tant que telle mais "déposé" et
« décompressé » sur l’ordinateur
o SPIRITS nécessite que JAVA soit installé : « Spirits requires a Microsoft
Windows system with a Java 7 Runtime Environment (JRE 1.7) or higher
properly installed. » (source : manuel officiel de SPIRITS)
https://www.java.com/fr/download/
Remarque : SPIRITS démarre via « javaw.exe « (localisé par exemple
dans « C:\Program Files\Java\jre1.8.0_171\bin\javaw.exe »)
R + Rstudio
o https://www.rstudio.com/
o https://www.r-project.org/
SNAP 7 (ESA)
o https://step.esa.int/main/download/snap-download/
o --> voir le tableau et la ligne "All Toolboxes"
o Faire les mises à jour des plugins après installation
432