Méthodes de Machine Learning Pour Le Suivi de L'occupation Du Sol Des Deltas Du Viêt-Nam
Méthodes de Machine Learning Pour Le Suivi de L'occupation Du Sol Des Deltas Du Viêt-Nam
L'UNIVERSITE
DE BRETAGNE OCCIDENTALE
ECOLE DOCTORALE N° 598
Sciences de la Mer et du littoral
Spécialité : Géomatique
Par
Chi-Nguyen LAM
Méthodes de Machine Learning pour le suivi de l’occupation du sol des
deltas du Viêt-Nam
Composition du Jury :
Rapporteurs avant soutenance :
Directrice de thèse
Dr Thomas CORPETTI, CNRS, LETG-Rennes UMR 6554, Dr. HDR Simona NICULESCU, LETG-Brest UMR 6554, UBO.
Université de Rennes 2
Co-directeur de thèse
Prof. Sébastien GADAL, CNRS ESPACE UMR 7300, Aix- Prof. Bernard POTTIER, Lab-STICC, Département
Marseille Université Informatique, UBO
Enfin et surtout, je tiens à remercier tous les membres de ma famille, qui m'ont
donné la motivation pour continuer sur la voie de la recherche scientifique. Merci à
mon épouse qui s'est occupée de la famille en mon absence qui a sympathisé et m'a
motivé à continuer ma passion scientifique. Je tiens à remercier tous Yao et mes amis
qui m'ont soutenu pendant mon séjour en France.
Chi-Nguyen LAM
2
Table de matière
Introduction générale .................................................................................................... 9
Contexte général et Problématique .......................................................................... 9
Objectifs ............................................................................................................ 10
Structure de la thèse .............................................................................................. 11
Partie 1 : Le système deltaïque vietnamien ................................................................ 12
Chapitre 1. Le système du Delta du Mékong ......................................................... 14
1.1. Type d’occupation du sol du delta du Mékong ................................................ 15
1.1.1. La guerre du Viêt-Nam ............................................................................ 15
1.1.2. Les activités aquacoles et les mangroves dans le delta du Mékong ........ 16
1.1.3. Activités agricoles dans le delta du Mékong ............................................ 17
1.1.4. Impact des conditions climatiques sur le delta du Mékong ...................... 18
1.2. Risque d’inondation dans le delta du Mékong ................................................. 18
1.2.1. La sécheresse ......................................................................................... 20
1.2.2. L’intrusion des eaux salines .................................................................... 21
Chapitre 2. Le système du fleuve Rouge ............................................................... 22
2.1. Type d’occupation du sol du delta du fleuve Rouge ........................................ 23
2.1.1. Les activités aquacoles ........................................................................... 23
2.1.2. Les activités agricoles ............................................................................. 23
2.1.3. La législation et Gestion spatiale du delta du fleuve Rouge ..................... 24
2.2. Impact des conditions climatiques sur le delta du fleuve Rouge ...................... 25
2.2.1. Risque d’inondation dans le delta du fleuve rouge................................... 25
2.2.2. Intrusion des eaux salines ....................................................................... 25
Chapitre 3. L’occupation du sol dans les deltas du Viêt-Nam ................................ 27
3.1. L’étude de changement de l’occupation/couverture du sol .............................. 27
3.2. Les données satellitaires et la détection de l’occupation du sol ....................... 28
3.3. L’étude du risque d’inondation ........................................................................ 28
3.3.1. Étude de l’occupation du sol pour l’évaluation du risque d’inondation...... 29
3.3.2. Les données satellitaires et l’évaluation du risque d’inondation ............... 30
Chapitre 4. Satellites Landsat et Sentinel .............................................................. 35
4.1. Le satellite Landsat ......................................................................................... 35
4.2. La constellation Sentinel ................................................................................. 36
Partie 2 : Démarches Méthodologiques ...................................................................... 39
Chapitre 5. Structuration des données et Prétraitement......................................... 40
5.1. Les images SAR : Sentinel-1 .......................................................................... 40
5.1.1. Les équations de radar ............................................................................ 40
5.1.2. Les images SAR...................................................................................... 43
5.1.3. Les images Sentinel-1 ............................................................................. 43
5.2. Les images optiques : Sentinel-2 .................................................................... 47
3
5.3. Les données de série temporelle .................................................................... 48
5.3.1. Les données optiques de série temporelle .............................................. 48
5.3.2. Les données radars de série temporelle .................................................. 48
5.3.3. Synergie des données Optiques et RADARS .......................................... 49
Chapitre 6. Machine Learning (Apprentissage automatique) ................................. 51
6.1. Définition Apprentissage Automatique ............................................................ 52
6.2. Définition Apprentissage Profond .................................................................... 53
6.3. La segmentation des images satellites............................................................ 54
6.3.1. L’algorithme K-means ............................................................................. 54
6.3.2. Meanshift................................................................................................. 57
6.3.3. Multi-résolution Segmentation (MRS) ...................................................... 60
6.3.4. Superpixel (SLIC) .................................................................................... 62
6.4. Les méthodes de classification ....................................................................... 64
6.4.1. Support Vector Machine : SVM ............................................................... 65
6.4.2. Random forest ......................................................................................... 68
6.5. Méthodes de Deep Learning ........................................................................... 70
6.5.1. Multilayer Perceptron – MLP ................................................................... 70
6.5.2. Convolutional neural network : CNN ........................................................ 76
6.6. Evaluation de la précision des classifications .................................................. 78
6.6.1. La matrice de confusion .......................................................................... 79
6.6.2. L’estimation des précisions ..................................................................... 81
6.6.3. L’estimation des superficies. ................................................................... 82
6.6.4. Les exemples .......................................................................................... 82
Chapitre 7. L’Analyse orientée objet d'images satellite .......................................... 87
7.1. Émergence de l'approche « orientée objet » ................................................... 87
7.2. Définition du concept de GEOBIA ................................................................... 87
7.3. Evaluation de GEOBIA.................................................................................... 88
7.3.1. Avantages de GEOBIA ............................................................................ 88
7.3.2. Inconvénient de GEOBIA ........................................................................ 88
7.4. Paradigme de GEOBIA ................................................................................... 88
7.4.1. Etape de travail dans le cadre de GEOBIA .............................................. 89
7.4.2. Méthodes de machine learning................................................................ 89
7.4.3. Méthodes de deep learning ..................................................................... 90
Partie 3 : De l’acquisition des données aux Applications ............................................ 91
Chapitre 8. Détection des changements de l’occupation du sol dans le delta du fleuve
Rouge (Viêt-Nam) ...................................................................................................... 92
8.1. Zone d'étude ................................................................................................... 92
8.2. Données mises en œuvre ............................................................................... 94
8.3. Méthodologie .................................................................................................. 96
8.3.1. Calcul du NDVI ........................................................................................ 96
4
8.3.2. Segmentation des images ....................................................................... 98
8.3.3. Classification des images ...................................................................... 100
8.4. Résultats et discussion ................................................................................. 102
8.4.1. Résultat des Random Forest ................................................................. 102
8.5. Conclusions .................................................................................................. 105
Chapitre 9. Étude et cartographie de l’inondation au niveau du delta du Mékong .106
9.1. Méthodologie ................................................................................................ 106
9.2. Les données SAR et la détection de l’extension de l’inondation : .................. 107
9.2.1. Le suivi des inondations des rizières : ................................................... 108
9.3. Zone d’étude ................................................................................................. 109
9.3.1. La culture du Riz au delta du Mékong ................................................... 110
9.3.2. Régime hydro-météorologique du Mékong ............................................ 110
9.4. Données utilisées.......................................................................................... 111
9.4.1. Les données satellitaires ....................................................................... 111
9.4.2. Données environnementales ................................................................. 112
9.5. Méthode d'apprentissage automatique et classification ................................. 113
9.5.1. Segmentation d'images et extraction de caractéristiques ...................... 113
9.5.2. Construction du modèle du CNN ........................................................... 114
9.5.3. Random Forest (RF) ............................................................................. 115
9.5.4. Évaluation de la précision ...................................................................... 115
9.6. Résultats ....................................................................................................... 116
9.6.1. Profils de rétrodiffusion de rizières ........................................................ 116
9.7. Validation et comparaison des modèles ........................................................ 120
9.7.1. L'évaluation de la précision du modèle CNN ......................................... 120
9.7.2. L'évaluation de la précision du classificateur RF ................................... 122
9.8. Résultats des modèles : cartographie inondée et inondable ......................... 125
9.8.1. Cartographie de l'étendue des inondations ............................................ 125
9.8.2. Cartographie des zones inondables ...................................................... 128
9.8.3. Discussion ............................................................................................. 131
Conclusion générale ..................................................................................................132
Bibliographie ...........................................................................................................135
ANNEXE ...........................................................................................................152
5
Liste des figures
Figure 1 : Système des deltas au Viêt-Nam..........................................................................13
Figure 2 : Photo de trois sorties d'avions C-123 ...................................................................16
Figure 3 : Le delta du fleuve Rouge au Vietnam ...................................................................22
Figure 4 : Polarisations horizontale (noire) et verticale (rouge) .............................................41
Figure 5 : Les modes du système Sentinel 1 ........................................................................44
Figure 6 : Graphe de prétraitement des images Sentinel-1 ..................................................45
Figure 7 : Illustration des relations entre apprentissage automatique ...................................51
Figure 8 : Le processus général de l'algorithme d'apprentissage supervisé..........................53
Figure 9 : L'architecture globale du Convolutional neural network (CNN) .............................54
Figure 10 : Le courbe des clusters k de 1 à 9, le coude est trouvé par K = 3 ........................56
Figure 11 : Image optique Sentinel-2 sur TranDe au 20 février 2020 ....................................58
Figure 12 : Kmeans (A) K=5 , (B) K=10 , (C) K=15 ...............................................................58
Figure 13 : La courbe des clusters k de 1 à 25 .....................................................................58
Figure 14 : Segmentation de MeanShift avec la bande passante .........................................60
Figure 15 : Segmentation par multi résolution ......................................................................62
Figure 16 : Segmentations par SLIC.....................................................................................64
Figure 17 : Hyperplan de SVM (a) données séparables, (b) données non-séparables ........66
Figure 18 : Hyperplan de SVM (a) données non-linéaires, (b) données transformées .........67
Figure 19 : Exemple de forêt d’ arbres de décision ...............................................................70
Figure 20 : Un neurone dans le cerveau humain ..................................................................71
Figure 21: Un neurone artificiel.............................................................................................72
Figure 22 : Graphiques des fonctions d'activation ................................................................72
Figure 23 : Modèle d'un réseau de neurones........................................................................73
Figure 24 : Illustration de l'étape Feedforward ......................................................................74
Figure 25 : Illustration de l'étape backpropagation. ...............................................................75
Figure 26 : Modèle CNN avec 02 couches convolutives .......................................................77
Figure 27 : Exemple de l'opération convolutive. ....................................................................78
Figure 28 : Fonctions Max-Pool et Avg-Pool.........................................................................78
Figure 29 : Carte de l'occupation du sol de province An Giang et province Dong Thap ........83
Figure 30 : L'occupation du sol en 2019 à gauche et en 2020 à droite. ................................85
Figure 31 : Relation entre les objets et la résolution spatiale. ...............................................87
Figure 32 : Processus de GEOBIA .......................................................................................89
Figure 33 : Modèle Générale de CNN. .................................................................................90
Figure 34 : Delta du fleuve Rouge au Viêt-Nam....................................................................93
Figure 35 : Chaîne de prétraitement et de traitement. ..........................................................98
Figure 36 :Segmentation ESP2. ...........................................................................................99
Figure 37 : Modèle MLP .....................................................................................................101
Figure 38 : Occupation du sol en 2016, 2017, 2018 à TienHai, ThaiBinh au Viet-Nam. ......104
6
Figure 39 : Prétraitement et traitement d’image SAR Sentinel-1. ........................................107
Figure 40 : Représentation schématique des mécanismes de rétrodiffusion ......................108
Figure 41 : Zone d'etudes - Province AnGiang et Province DongThap ...............................110
Figure 42 : Classificateur CNN généralisé avec l'architecture 2D-CNN. .............................114
Figure 43 : Profil temporel VH et VV. ..................................................................................118
Figure 44 : a) Pluies quotidiennes et niveaux d'eau quotidiens du Mékong ........................119
Figure 45 : Cartographie des inondations par classification CNN par mois ........................126
Figure 46 : schéma synthétique a) Précipitations mensuelles moyennes ..........................127
Figure 47 : Cartographie des inondations. ..........................................................................130
Figure 48 : Points d’observation et de validation dans le delta du Mékong. ........................152
Figure 49 : Profils VH temporel des classes de l’occupation du sol en 2019. ......................153
Figure 50 : Profils VV temporel des classes de l’occupation du sol en 2019. ......................153
Figure 51 : Profils VH temporel des classes de l’occupation du sol en 2020. ......................154
Figure 52 : Profils VV temporel des classes de l’occupation du sol en 2020 .......................154
7
Liste des tableaux
Tableau 1 : Polarisations utilisées ........................................................................................42
Tableau 2 : Modes de système Sentinel-1 ............................................................................44
Tableau 3 : Récapitulatif des résolutions spectrale, spatiale et temporelle de Sentinel 2 ......47
Tableau 4 : Séries temporelles 2016, 2017,2018 au delta du fleuve Rouge. .......................48
Tableau 5 : Organisation des jeux de données SAR ............................................................49
Tableau 6 : Noyaux de SVM non linéaire..............................................................................67
Tableau 7 : Matrice de confusion de k classes. ....................................................................80
Tableau 8 : Etendue de la matrice de confusion ...................................................................80
Tableau 9 : Matrice de confusion estimée. ...........................................................................81
Tableau 10 : Matrice de confusion de la carte d'occupation du sol. ......................................83
Tableau 11 : Précisions des classifications de l'occupation du sol. .......................................83
Tableau 12 : Matrice de confusion estimée ..........................................................................84
Tableau 13 : Estimations des précisions de la carte de l'occupation du sol. .........................84
Tableau 14 : Matrice de changements de l'occupation du sol de 2019 à 2020. ....................85
Tableau 15 : Précisions estimées de la classification par CNN.............................................86
Tableau 16 : Ensemble de données satellitaires disponibles. ...............................................95
Tableau 17 : Résolutions spatiales et spectrales du satellite Sentinel-2. ..............................95
Tableau 18 : Résolutions spatiales et spectrales du satellite Landsat-8. ..............................96
Tableau 19 : Calendrier du riz dans la zone de la côtière du delta du fleuve Rouge .............97
Tableau 20 : Calendrier des cultures de riz et passages des Satellites ................................98
Tableau 21 : Précision globale et indice de kappa par classe .............................................103
Tableau 22 : Comparaison des précisions globales par le modèle CNN. ...........................120
Tableau 23 : Précisions de l'utilisateur par le modèle CNN ................................................121
Tableau 24 : Précisions du producteur par le modèle CNN ...............................................121
Tableau 25 : Comparaison des précisions globales par le modèle RF................................123
Tableau 26 : Précisions de l'utilisateur par le modèle RF....................................................123
Tableau 27 : Précisions du producteur par le modèle RF ...................................................124
Tableau 28 : Analyse spatiales des zones impactées par les inondations. .........................127
Tableau 29 : Superficie de LULC avant la saison des pluies en 2019 et 2020. ...................128
Tableau 30 : Superficie des LULC non inondés en 2019 et 2020 .......................................130
Tableau 31: Nombre de points d’observation. ....................................................................152
Tableau 32 : Code Python du modèle Multilayer Perceptron ..............................................155
Tableau 33 : Code Python du modèle Convolutional Neural Network .................................155
8
Introduction générale
Les deltas sont parmi les zones les plus importantes au monde, fournissant une
gamme variée de services écosystémiques essentiels au bien-être socio-économique.
Ces deltas sont, en effet, les plus propices à l'habitation et l’occupation humaine sur la
Terre. Elles connaissent d’importants changements de l’occupation et de couverture
du sol en raison de la concentration de la population et des activités humaines (Bianchi
et Allison, 2009). Ces occupations peuvent aussi bien avoir un impact positif et
bénéfique pour le bien-être social et économique, comme elles peuvent avoir un
impact négatif sur l’écosystème naturel et les questions environnementales (Lambin
et Geist, 2006). Par ailleurs ces zones, qui sont le reflet des besoins humains, sont
sévèrement affectées par les changements environnementaux et les activités
anthropiques (Tilman, 2001). En effet, dans ce contexte de changements
environnementaux et socio-économiques, l'étude de l’occupation du sol au niveau des
zones deltaïques a dûment attiré l'attention des scientifiques (Minderhoud et al., 2018).
Bien que ces données présentent un grand potentiel pour la prise de décisions
dans des domaines critiques tels que la gestion des risques (inondation, sécheresse,
érosion, etc.), elles ne sont pas exploitées à leur plein potentiel en fonction des besoins
spécifiques (Kansakar et Hossain, 2016). En effet, le développement des techniques
de traitement automatique et reproductible des données de télédétection, est
nécessaire pour rendre le suivi des changements au niveau des zones deltaïques plus
pertinent.
Objectifs
10
celles de la segmentation d’image et de l’apprentissage automatique peu
profond et profond.
Structure de la thèse
Pour mener à terme notre travail et atteindre ces objectifs, cette thèse se
présente en trois parties.
11
Partie 1 : Le système deltaïque vietnamien
12
Les deltas et les fleuves sont parmi les zones les plus attractives pour
l’occupation et la population humaine. Par conséquent, elles connaissent une forte
concentration des activités socio-économiques à fort potentiel, qui se développent et
évoluent en fonction des besoins humains et des changements environnementaux
globaux (Tilman, 2001). Ces occupations peuvent aussi bien avoir un impact positif et
bénéfique pour le bien-être social et économique comme elles peuvent avoir un impact
négatif sur l’écosystème naturel et les questions environnementales (Lambin et Geist,
2006)
13
Chapitre 1. Le système du Delta du Mékong
Le fleuve Mékong est parmi les plus grands à l’échelle mondiale, et le troisième
fleuve le plus long d’Asie. Il prend sa source dans la chaîne de l’Himalaya et se jette
dans la mer de Chine au Viêt-Nam, en parcourant plus de 4 000 km et en irriguant 6
pays. En Chine, où il est appelé Lancang Jiang, il traverse la région du Kham entre le
Tibet et le Sichuan, puis il pénètre en Birmanie, et au Laos où il prend le nom de
Mékong (le fleuve mère), pour tracer enfin vers la Thaïlande.
Avec son climat tropical humide, le delta du Mékong est caractérisé par des
températures mensuelles moyennes constamment élevées (25-29°C) et de fortes
précipitations saisonnières (1 200-2 300 mm) (Hung, 2017), par conséquent le débit
du fleuve Mékong présente de fortes variations saisonnières.
Par ailleurs, le Viêt-Nam est un pays agricole d’excellence, selon les données
de l'Organisation des Nations Unies pour l'alimentation et l'agriculture (FAO) il est
classé cinquième producteur de riz et troisième exportateur mondial de produits
aquatiques (FAO, 2016). Le delta de Mékong est l’une des régions les plus fertiles, et
l'une des zones agricoles les plus productives au Viêt-Nam (Wilder et Phuong, 2002)
14
et du monde. Ses terres fertiles ainsi que le climat chaud et humide du sud-est
asiatique ont largement favorisé la riziculture, d’où le surnom de « Bol de riz du Viêt-
Nam ». Par ailleurs, le delta du Mékong contient les plus grandes zones aquacoles,
qui représentent 70 % de sa superficie, en fournissant plus de 52 % de la production
totale du delta (Truong et Do, 2018). Elle contribue à environ 55 % de la production
aquatique nationale et à environ 51 % des exportations vietnamiennes des produits
aquatiques (Quang, 2002).
Dans les années 1954 à 1975, pendant que le Viêt-Nam était en guerre, l’armée
américaine organisait des attaques chimiques et pulvérisait des défoliants sur de
grandes surfaces terrestres (Figure 2) au sud du pays, y compris les forêts, les
espaces boisés, les marécages et les mangroves (Dorian, 2015). Le but était
d'éradiquer les abris militaires et les sites de stockage de nourriture vietnamiens mais
ceci n’a pas été sans conséquence sur les écosystèmes et l’environnement. Le delta
du Mékong n’a pas échappé à ces attaques chimiques, car il y abritait de nombreuses
stations d'approvisionnement en nourriture ainsi que des bases militaires pour le
déploiement stratégique.
15
En 1968, environ 600 000 ha de terres ont été pulvérisés avec des produits
chimiques et plus de 20 % des zones forestières du sud du Viêt-Nam ont été défoliées,
dont environ 104 939 ha (36 %) de mangroves (NAS, 1974).
Figure 2 : Photo de trois sorties d'avions C-123 dispersant des herbicides lors d'une mission Ranch Hand
au-dessus du Sud-Vietnam
Le delta du Mékong possède un littoral de plus de 600 km, avec des conditions
naturelles favorables au développement de l'aquaculture. En effet, à partir de 1980
des bassins aquacoles sont apparus dans certaines zones côtières du delta du
Mékong, et l'aquaculture a commencé à se développer à un rythme important, en dépit
des mangroves.
16
Pour encadrer cette expansion et restaurer les systèmes de mangrove, une
réglementation a été instaurée en 1991 (64/[Link] (28/3/1991), autorisant les
agriculteurs à consacrer 30 % des terres allouées à l'aquaculture ou à l'agriculture, et
à préserver ou restaurer plus de 70 % des forêts de mangrove. Cela a conduit à
l'émergence d'un nouveau type d’occupation du sol, qui est le système complexe
mangrove-aquaculture (Ha et al., 2014). Cependant, ce système de gestion n'a pas
permis de réduire la surface d'aquaculture, où il a été constaté un coût écologique
important dans le delta du Mékong. Le gouvernement vietnamien a augmenté les
contrôles.
En effet, dans les années 2000 et 2001, les résolutions 09/NQ-CP et 1116,
respectivement, ont encouragé les agriculteurs à convertir plutôt les terres salines
côtières en fermes aquatiques en développant les activités de pisciculture et l'élevage
de crevettes. Avec la demande croissante en produits aquacoles, les agriculteurs sont
passés de l’agriculture à l’aquaculture (Trung et al., 2007).
17
Après 1998, une transformation massive des terres de culture dans le delta du
Mékong en aquaculture a été remarquée. En conséquence, les superficies agricoles
ont rapidement diminué.
Dans le delta du Mékong, les inondations sont causées par les débits élevés du
fleuve résultant de fortes pluies sur le bassin versant. Ces inondations sont aggravées
par l'élévation du niveau de la mer (Gupta, 2009). Le changement climatique entraîne
un déséquilibre du régime des précipitations, ce qui provoque des changements dans
le cycle hydrologique. Combiné au réchauffement qui a des impacts directs sur
l'évapotranspiration, la génération de ruissellement peut être modifiée (Vörösmarty et
al., 2000). D'autre part, les activités anthropiques telles que les barrages et l'irrigation
pourraient perturber le débit naturel des rivières et modifier les caractéristiques des
régimes hydrologiques (Hall et al., 2014).
1 Accord intervenu entre les agences spatiales participantes visant à fournir de l'information et des données
obtenues de l'espace pour appuyer les efforts humanitaires lors de situations d'urgence causées par des
catastrophes majeures.
19
1.2.1 La sécheresse
Durant l’année 2020, plus de 362 000 ha de rizières auraient été détruits, selon
les estimations du gouvernement vietnamien, à cause du phénomène de la
sécheresse qui affecte le delta du Mékong.
20
de manque d'eau deviendront graves. Ce sont les conditions initiales de la sécheresse.
Et si la sécheresse dure, ses impacts sur les secteurs socio-économiques seront plus
grands et resteront plus longtemps (IMHEN, 2011).
Les eaux salines pénètrent à l'intérieur des terres par les différents bras du
Mékong et les canaux sur 20 à 65 km de la côte (Trung, 2006). La mer de Chine
méridionale a une marée semi-diurne, c.à.d. deux pleines mers et deux basses mers
par jour avec un marnage quotidien moyen entre 3,5 m et 4,5 m dans la mer de l'Est
(NEDECO, 1993). Pendant la saison sèche, le débit principal du fleuve Mékong est
assez faible et son niveau d'eau est bas. En outre, durant la pleine mer, l'eau saline
s'infiltre profondément à l'intérieur des terres, sur une superficie de 500 000 ha sur
toute la zone du delta du Mékong ce qui a de graves conséquences sur les rizicultures
et les autres cultures (White, 2002).
21
Chapitre 2. Le système du fleuve Rouge
Cette région est caractérisée par un climat tropical de mousson avec une saison
humide, de mai à octobre, et une saison sèche, de novembre à avril. La saison estivale
est chaude et très humide, avec des températures moyennes allant de 27°C à 29°C,
tandis que la saison hivernale est fraîche et sèche avec des températures mensuelles
moyennes allant de 16°C à 21°C (Li et al., 2006).
22
La majorité de la population et de l'activité économique de la moitié nord du Viêt-
Nam est concentrée dans le delta du fleuve Rouge et près de la moitié de la production
de riz du pays y est cultivée. En effet, la production alimentaire dans le delta est
principalement basée sur la récolte de riz durant la seconde moitié de l'année
(Niculescu et Lam, 2019).
L'utilisation actuelle des terres dans la région côtière est orientée vers
l'agriculture, avec le riz comme principale culture. Deux cultures de riz et une culture
sèche sont produites par an (Do et Egashira, 2002) (patates douces, arachides, tabac,
maïs, etc.). Dans certaines zones très basses, salées et récemment défrichées, on
23
cultive temporairement du scirpe pour améliorer la qualité du sol avant de le mettre en
culture.
Au contraire, les zones résidentielles sont situées sur des terrains plus élevés,
hors de portée des inondations. On y trouve des villages où ils alternent la culture des
arbres à fruits et arbres d'ombrage comme le bananier, l'eucalyptus et le filao. Dans
les zones côtières salines et les estuaires, la végétation se développe naturellement
pour former des forêts. Cependant, la formation végétale prédominante est la forêt de
palétuviers, principalement du type Sonneratia, avec une hauteur allant de 1 à 3,5 m.
Dans certaines zones de haute altitude, le filao est cultivé pour la production de bois
et pour la protection des villages de l'intérieur contre les tempêtes.
24
récolte. Ce décret est appliqué dans le sud du Viêt-Nam depuis 1988. Il n'est entré en
vigueur que deux ans plus tard dans le nord du pays, où les dirigeants conservateurs,
cherchant à conserver leur pouvoir et leurs privilèges, se sont opposés à toute forme
de réforme de la gestion des terres.
Les intrusions des eaux salines sont principalement dues aux marées hautes.
Dans le Golfe du Tonkin la marée est de type diurne, avec un seul cycle (une marée
descendante et une marée montante) par jour et une amplitude qui suit un cycle lunaire
de 14 jours et diminue progressivement de 4 à 2 m du nord au sud (Fang et al., 1999).
25
Le niveau de la marée à basse altitude, à l'embouchure de la branche principale du
fleuve Rouge, varie de 0,5 m pendant les marées de mortes eaux à 2,5 m pendant les
marées de vives eaux.
26
Chapitre 3. L’occupation du sol dans les deltas du Viêt-Nam
Les zones deltaïques sont des écosystèmes très productifs qui fournissent un
moyen de subsistance aux populations locales, néanmoins elles sont considérées
comme étant biologiquement et écologiquement complexes (Cosslett et Cosslett,
2014). Au Viêt-Nam, les évolutions actuelles du contexte, à la fois climatique,
économique et sociétal, font de la gouvernance des changements de l’occupation des
sols au niveau des systèmes deltaïques une question fondamentale. Le modèle de
couverture des sols dans ces zones côtières/deltaïques, a changé de manière
significative au cours des décennies (Sakamoto et al., 2009) (voir titres 2.1 et 3.1).
L’analyse se base sur des données multi-date qui sont des sources
d'information précieuses pour connaître la trajectoire du changement du type
d'utilisation et de couverture des sols, tout en décelant les facteurs entraînant la
destruction des écosystèmes, tels que l’expansion de l’aquaculture et l’agriculture en
dépit du système naturel des deltas du Viêt-Nam.
Par ailleurs, plusieurs études ont été menées pour détecter l'expansion de
l'élevage de crevettes à l'aide de la télédétection optique (Tong et al., 2004 ;
Sakamoto et al., 2009 ; Vo et al., 2013). Dans ce double sens agriculture-aquaculture,
des études concernant les changements spatio-temporels des systèmes d'aquaculture
intensive à l'aide de techniques de télédétection, telle que celle de Sakamoto et al.
(2009), ont montré que l'aquaculture s'est rapidement développée, et que la proportion
de riz à triple culture a diminué. Ce constat a été révélé grâce aux données satellitaires
disponibles depuis 1972.
Nous détaillons les propriétés de ces capteurs satellites, dont les données sont
largement utilisées pour les études de l’occupation du sol dans les deltas du Viêt-Nam,
dans le chapitre qui suit.
Les inondations à travers le monde, sont parmi les catastrophes naturelles les
plus destructives affectant, plus particulièrement, les terres basses, à l’instar des
zones côtières. Les pays d'Asie du Sud-Est sont vulnérables aux inondations,
notamment pendant la saison des pluies. Ainsi, ce risque constitue une des plus
28
grandes menaces du changement climatique dans les deltas du Viêt-Nam (Voir titres
2.2 et 3.2). Les impacts et la fréquence des inondations devraient augmenter à l'avenir
en raison du développement urbain et du changement climatique (Kundzewicz et al.,
2014). En effet, des organisations telles que la Commission du Mékong et le Centre
asiatique de préparation aux catastrophes mettent en œuvre des systèmes régionaux
de prévision des inondations en combinant de manière synergique des données et des
modèles hydrologiques (Ahmed et al., 2017).
29
d’inondation. Dans le delta du Mékong, l’estimation des étendues d’eau (zonage) peut
affecter les rizières qui constituent un enjeu économique incontournable pour le pays,
touchées par les inondations, ces enjeux socio-économiques sont devenus
vulnérables.
Le but est de développer des prévisions fiables des inondations. Pendant les
dernières décennies, il y a eu des avancées majeures dans le domaine de la
30
télédétection, en particulier de la télédétection de micro-ondes, et de la modélisation
des inondations.
Une deuxième catégorie de données radar est représentée par les données
InSAR interférométrie qui part de deux acquisitions d’images complexes (SLC) avec
des angles d’enregistrement légèrement différents. Après l’enregistrement des deux
images initiales, on calcule la corrélation croisée complexe normalisée. Sa magnitude,
31
appelée cohérence et sa phase sont les critères les plus connus en interféromètrie
SAR (Maître, 2001). La cohérence est utilisée dans la plupart des cas pour la
description de la stabilité temporelle du signal acquis. Alors que la phase inclut une
composante géométrique liée à la distance entre la cible et le capteur en donnant des
informations sur l’altitude et le déplacement (éventuel) de la cible (Massonet et
Rabaute, 1993, Massonet et Feigle, 1998). Afin d’éviter la phase de l’atmosphère,
Ferreti et al. (2001) propose l’analyse de toute une série d’interférogrammes et la
phase des pixels stables. Cette technique porte le nom de Permanent Scatterer
Technique. A la base de cette analyse, c’est le choix des pixels en fonction de leur
amplitude pour la mesure des déplacements. Dans le domaine de surveillance des
inondations, le potentiel de l’interférométrie SAR a été démontré depuis les années 90
(Geutner et al., 1996). Les cartes de cohérence interférométrique InSAR ont donné
des informations utiles sur l’ampleur des inondations du fleuve Rhin.
Selon certains chercheurs (Vidal, 2005), pour être fiables, les modèles
hydrauliques doivent être contraints, par le biais d’un calage, grâce à des
connaissances concernant la crue (niveaux d’eau, volume transitant de la plaine au
moment du passage de la crue, etc.). Selon d’autres auteurs (Aronica et al., 1998 ;
32
Pappenberger et al., 2005 ; Hostache et al., 2006 ; 2010) ces données sont la plupart
du temps difficiles à acquérir sur le terrain. Une solution peut être envisagée :
l’exploitation des niveaux d’eau en période d’inondation estimés à partir d’images
satellites radar. Des méthodes basées sur la mesure du coefficient de rétrodiffusion
radar sur des données à haute résolution ont fait l’objet de développement en bande
L (Ulaby et al., 1983 ; Kim et al., 2014 ; Chapman et al., 2015) et en bande C (Fellah
et al., 1994, 2003 ; Brown et al., 1997 ; Matgen et al., 2007 ; Ramsey et al., 2013) ;
des niveaux d’eau ont été estimés par croisement entre des limites de zones inondées
extraites d’images satellitaires radar toujours en bande C et des cartes topographiques
(Brakenridge et al., 1998 ; Hostache, 2009 ; Mason et al., 2009, 2012). Les incertitudes
obtenues (entre 1 à 3 m) ont été considérées comme incompatibles avec la précision
requise par les modélisations hydrauliques. Pour améliorer l’extraction du rivage dans
des secteurs où la végétation est inondée, Horrit et al. (2002) utilisent la taille de la
végétation dérivée des données du LiDAR. Une approche similaire est développée par
Puech et Raclot (2002) à partir de photographies aériennes, Hostache et al. (2006,
2009, 2010) appliquent une méthode de cartographie sur une image SAR en utilisant
des données photogrammétriques de haute précision (incertitude de 30 cm en
moyenne).
33
Ces techniques incluent : simple interprétation visuelle (Oberstadler et al., 1997),
seuillage d’histogrammes d’images (Brivio et al., 2002 ; Matgen et al., 2004),
algorithme de classification automatique (Bonn et Dixon, 2005 ; Pulvirenti et al., 2011 ;
Kuenzer et al., 2013 ; Niculescu et al., 2009, 2015), algorithme de texture d’image
(Schumann et al., 2009), méthodes de détection multitemporelle de changements
(Laugier et al., 1997 ; Niculescu et al., 2009, 2010 ; 2015 ; Martinis et Rieke, 2015).
De nombreuses études (Marinelli et al., 1997; Nico et al., 2000 ; Yésou et al., 2000 ;
Alsdorf et al., 2001) montrent l’apport important de la variation de cohérence de phase,
obtenue à partir d’une paire interférométrique d’images radar, en complétant l’étude
de l’amplitude du signal.
Les images optiques disponibles gratuitement, telles que Landsat (TM) peuvent
être utiles pour l’analyse du risque grâce à leur résolution spatiale et leur intervalle
d’acquisition, mais aussi grâce à l’archive et la capture continue sur les mêmes zones
qui rendent ces données très pratiques pour la cartographie des masses d'eau. En
effet, les données de la bande 7 de Landsat sont utiles pour distinguer les plans d'eau
de la masse terrestre environnante, avec un taux d'erreur d'environ 5 % (Dhara et al.,
2020). Bien que les données Landsat aient une résolution spatiale relativement plus
élevée que les données MODIS, leur temps de revisite est plus long que celui du
MODIS. Ces derniers comprennent deux instruments clés, les satellites Terra et Aqua.
Terra MODIS capture des images de la surface de la terre tous les jours ou tous les
deux jours, et acquiert des données dans 36 bandes spectrales (Ticehurst et al., 2014;
Revilla-Romero et al., 2015) ce qui reste très pratique pour la cartographie des zones
inondables. Par ailleurs, des images optiques Sentinel-2 A et B, ont également été
utilisées dans plusieurs études des risques d’inondation.
34
Chapitre 4. Satellites Landsat et Sentinel
Il faut noter aussi que parmi les facteurs qui ont contribué au succès technique
des données Landsat c’est la combinaison de capteurs avec des domaines spectraux
façonnés pour l'observation de la Terre, ce qui est efficace dans une grande variété
de contextes d'application. Par ailleurs, le programme a permis d'accumuler des
archives de données volumineuses et de maintenir des données historiques. Ces
données historiques sont une exigence fondamentale pour la surveillance de
l'environnement afin de prévenir les impacts environnementaux néfastes avant qu'ils
ne deviennent irréparables (Manandhar et al., 2009). Entre 1972 et 2013, huit satellites
Landsat sont lancés, permettant de capturer plusieurs millions d'images. L’archive de
ces données est fournie par le Centre national d'observation et de science des
ressources terrestres (USGS EROS) du Service géologique des États-Unis.
35
Le satellite Landsat-8, une collaboration entre la NASA et l’United States
Geological Survey (USGS), a été lancé en février 2013. Il couvre la Terre tous les
16 jours avec des images de 185 km x 185 km, en 16 bits, comptant 11 bandes
spectrales : 9 dans le visible (8 multispectrales de résolution 30 m ; 1 panchromatique
à 15 m) et 2 thermiques (60 m). L'équipement de télédétection installé sur le LCDM
(Landsat Data Continuity Mission) comprend un radiomètre à balayage multicanaux
OLI (Operational Land Imager) et un radiomètre IR à deux canaux TIRS (Thermal
Infrared Sensor). L'imageur terrestre opérationnel (OLI) offre deux nouvelles bandes
spectrales par rapport à l'instrument Landsat-7 ETM+, l'une spécialement conçue pour
la détection des cirrus (bande 9, nouvelle bande proche infrarouge (NIR)) et l'autre
pour l'observation des zones côtières (bande 1, nouveau canal bleu dans le visible). Il
mesure dans les parties visibles, NIR et infrarouge à ondes courtes (SWIR) du spectre
électromagnétique et offre une résolution spatiale panchromatique de 15 mètres et
multi spectrale de 30 mètres. Le capteur infrarouge thermique (TIRS) a été ajouté à la
charge utile de Landsat 8 afin de poursuivre l'imagerie thermique et de prendre en
charge des applications émergentes telles que les mesures du taux
d'évapotranspiration pour la gestion de l'eau. Le TIRS a été construit par le Goddard
Space Flight Center de la NASA, et sa durée de vie est de trois ans. Les données
TIRS de 100 m sont enregistrées avec les données OLI pour créer des produits de
données Landsat-8 de 12 bits corrigés sur le plan radiométrique, géométrique et du
terrain. Les images des satellites LANDSAT sont libres de droit et diffusées par le US
Geological Survey (USGS). Dans nos applications les images Lansat 8 (OLI) sont
utilisées pour l’analyse et la cartographie des changements de l’occupation du sol dans
le delta du fleuve Rouge.
36
l'atmosphère, des océans et des terres, et les données recueillies sont utiles dans de
nombreuses applications. Les missions Sentinel embarquent une gamme de
technologies telles que des instruments radar et d'imagerie multi-spectrale. En effet,
la mission Sentinel-1 (Sentinel-1A et Sentinel-1B) embarque une technologie radar en
orbite polaire, jour et nuit, destinée aux services terrestres et océaniques. Sentinel-1A
a été lancé le 3 avril 2014 et Sentinel-1B le 25 avril 2016. Sentinel-2 (Sentinel-2A et
Sentinel-2B) est une mission d'imagerie multi-spectrale haute résolution en orbite
polaire pour la surveillance des terres. Elle fournit, par exemple, des images de la
végétation, de la couverture du sol et de l'eau, des voies navigables intérieures et des
zones côtières. Sentinel-2 peut également fournir des informations aux services
d'urgence. Sentinel-2A a été lancé le 23 juin 2015 et Sentinel-2B a suivi le 7 mars
2017. Sentinel-3 (Sentinel-3A et Sentinel-3B) est une mission multi-instruments qui
soutient les systèmes de prévision océanique, ainsi que la surveillance de
l'environnement et du climat. Elle vise à mesurer la topographie de la surface de la
mer, la température de la surface de la mer et de la terre, la couleur de l'océan et de
la terre avec une précision et une fiabilité de haut niveau. Il existe d’autres missions à
divers objectifs opérationnels comme le Sentinel-4 et le Sentinel-5 qui sont une charge
consacrée à la surveillance de l'atmosphère, ainsi que le Sentinel-6 qui a été mis en
orbite le 21 novembre 2020, il transporte un altimètre radar pour mesurer la hauteur
de la surface de la mer à l'échelle mondiale, principalement pour l'océanographie
opérationnelle et les études climatiques.
Par ailleurs, des études ont testé la capacité du Sentinel-1 (SAR) dans le delta
du Mékong ce qui a été particulièrement pertinent pour cette région tropicale, car les
capacités tout temps du SAR permettent de faire face à la fréquente couverture
nuageuse qui rend la télédétection optique difficile. En effet, le SAR joue un rôle crucial
dans la télédétection en raison de sa capacité de pénétration des nuages et de
détection (Bouvet et Le Toan, 2011).
37
En outre, Sentinel-1 est le seul satellite SAR disponible gratuitement et constitue
donc une ressource précieuse pour les pays équatoriaux à revenus faibles et moyens
dont les budgets pour l'achat de données satellitaires sont limités.
38
Partie 2 : Démarches Méthodologiques
39
Chapitre 5. Structuration des données et Prétraitement
Dans ce travail, nous avons utilisé des données SAR de Sentinel-1. Avant la
présentation de ces données, une explication du système radar s’avère nécessaire.
𝐺 2 𝜆2 𝜎
𝑃𝑟 = 𝑃𝑒 (1)
(4𝜋)3 𝑅 4
Avec :
𝜆 : la Longueur d’onde
40
𝑟 2 𝑒 2
⃗⃗⃗ | et |𝐸
|𝐸 ⃗⃗⃗ | représentent de l’intensité de l’onde reçue et de l’onde émise.
Figure 4 : Polarisations horizontale (noire) et verticale (rouge) d'une onde électromagnétique plane
L'onde électromagnétique (𝐸⃗ 𝑒 ) est émise par une source en deux polarisations
V et H puis l’onde reçue (𝐸⃗ 𝑟 ) est diffusée par une cible. Les deux ondes sont
décomposées comme suit :
𝑖2𝜋𝑟
𝑒𝑥𝑝(− ) (4)
𝐸⃗ 𝑟 = 𝜆 . |𝑆|.⃗⃗⃗⃗𝐸 𝑒
𝑟
41
En remplaçant des vecteurs de formule (3) par la formule (4) et une
décomposition|𝑆|, la matrice de rétrodiffusion est donnée pour définir le processus de
transformation de l’onde (𝐸⃗ 𝑒 ) en une onde réfléchie (𝐸⃗ 𝑟 ). Elle est définie comme suit :
(Cloude et Zebker, 2010 ; Lee et Pottier, 2009).
𝑖2𝜋𝑟
𝑒𝑥𝑝(− ) (5)
[𝐸𝐻𝑟 𝐸𝑉𝑟 ] = 𝜆 . [𝑆 𝑆𝐻𝑉 𝑆𝑉𝐻 𝑆𝑉𝑉 ][𝐸𝐻𝑒 𝐸𝑉𝑒 ]
𝐻𝐻
𝑟
Les ondes radars sont caractérisées par deux paramètres : la fréquence de
fonctionnement, souvent exprimée en GigaHertz (GHz), et l’état de polarisation. Par
convention, le champ électrique est pris comme référence pour déterminer l’état de
polarisation défini par la forme de la courbe décrite, en fonction du temps, par la
projection de l’extrémité du champ électrique dans le plan transverse de propagation
(Frédéric Baup, 2007).
Polarisation simple HH ou VV ou HV ou VH
42
a) Rétrodiffusion directe de surface provenant du sol ; b) rétrodiffusion
volumique provenant du feuillage des arbres, arbustes et de la litière végétale si
présente ; c) rétrodiffusion directe de surface provenant de tiges, de branches
principales, de troncs ; d) rétrodiffusion double-bond provenant de l’interaction entre
les troncs et les autres composantes ; e) rétrodiffusion multiple provenant de
l’interaction du signal incident avec plusieurs composantes.
Dans ce travail, les données Sentinel-1 GDR sont utilisées pour l’étude et la
détection du type de l’occupation du sol. La série temporelle de Sentinel-1 offre une
opportunité unique de surveiller systématiquement les occupations du sol à un cycle
de répétition hebdomadaire. En outre, la continuité des données Sentinel-1 est
garantie jusqu'en 2030. Plusieurs études (ex. Abdikan et al., 2016; Balzter et al., 2015;
Lee et Pottier, 2009; Longepe et al., 2011) ont pu démontrer les différents avantages
des données SAR.
Les données Sentinel-1 sont distribuées par le Copernicus Open Access Hub
et comprennent :
43
Niveau-2 ou Level-2 : Ocean (OCN) : ce sont des données pour la
détection des paramètres géophysiques de l'océan.
44
Dans un graphe, le début du processus consiste à utiliser des fonctions qui
permettent de lire les images et d’enregistrer des résultats. Ainsi, les fonctions qui
permettent de traiter les signaux sont mis en ordre sériel ou parallèle. Dans nos études,
les fonctions utilisées sont :
Figure 6 : Graphe de prétraitement des images Sentinel-1 à servir de détecter des occupations du sol
dans le delta du Mékong
45
Correction de terrain: les corrections géométriques.
Les données SAR sont généralement captées avec un angle de vue variable
supérieur à 0 degré, ce qui donne des images présentant une certaine distorsion liée
à la géométrie de l'orientation latérale. Les corrections de terrain visent à compenser
ces distorsions afin que la représentation géométrique de l'image soit aussi proche
que possible du monde réel. La correction de terrain Doppler à distance est une
correction des distorsions géométriques causées par la topographie, telles que le
raccourcissement et les ombres, en utilisant un modèle numérique d'élévation pour
corriger l'emplacement de chaque pixel. L'opérateur de correction de terrain par effet
Range Doppler met en œuvre la méthode d'ortho rectification (Schubert et al., 2015)
pour le géocodage de scènes SAR à partir d'images en géométrie radar.
Le Speckle, qui apparaît dans les images SAR comme un bruit granulaire, est
dû à l'interférence des ondes réfléchies par de nombreux diffuseurs élémentaires
(Lee et al., 1994a). Le Speckle filter est une procédure qui permet d'améliorer la
qualité de l'image en réduisant le Speckle (bruit). Lorsqu'une telle procédure est
effectuée à un stade précoce du traitement des données SAR, le Speckle n'est pas
propagé dans les processus en cours (c'est-à-dire, Correction de terrain ou
conversion en dB). Le Speckle filter n'est pas conseillé lorsqu'il s'agit d'identifier de
petites structures spatiales ou la texture de l'image, car il risque de supprimer ces
informations. Le filtre de Lee raffiné s'est avéré plus efficace, comparant aux autres
Speckle filter à produit unique, pour l'interprétation visuelle, en raison de sa capacité
à préserver les bords, les caractéristiques linéaires, les cibles ponctuelles et les
informations de texture (Lee et al., 1994b). Plus récemment, des Speckle filter de
séries temporelles ont été développés pour réduire le bruit (Speckle), en tirant parti
des multi-observations SAR dans le temps. L'outil SNAP met à disposition plusieurs
filtres de Speckle : 'Boxcar', 'Median', 'Frost', 'Gamma Map', 'Lee', 'Refined Lee', 'Lee
Sigma', 'IDAN'.
46
Conversion to décibel (dB) : les transformations logarithmiques dB.
De plus, les produits Sentinel-1 GRD peuvent être de données MSI de Sentinel-
2, afin de promouvoir l'utilisation de constellations virtuelles de satellites par le biais de
la fusion de données. Dans le processus de prétraitement on peut inclure la fusion des
données SAR avec des données optiques.
Sentinel-2
Largeur λ
20 65 35 30 15 15 20 115 20 20 30 90 180
(nm)
Résolution
spatiale 60 10 10 10 20 20 20 10 20 60 60 20 20
(m)
Résolution
5 jours
spatiale
47
principaux niveaux disponibles pour notre zone d'étude sont de niveau 1C (le niveau
1C fournit des images ortho rectifiées) et 2A.
Les images de séries temporelles fournies par les satellites Sentinel-2 offrent
une opportunité unique pour la cartographie de la végétation (Mercier et al., 2019). Le
satellite Sentinel-2, a un temps de revisite de trois jours et une résolution spatiale de
10 mètres. Cependant, l'acquisition d'images optiques au cours des périodes clés de
surveillance peut être limitée en raison de leur vulnérabilité au temps pluvieux ou
nuageux. Dans ce contexte, le radar à synthèse d'ouverture, indépendant de
l'illumination solaire et de la couverture nuageuse, est une forme de télédétection
active. Le radar à synthèse d'ouverture (SAR) peut être utilisé comme une importante
source de données alternative ou complémentaire (Bovenga, 2020; Moreira et al.,
2013; Schmullius et Evans, 1997; Tomiyasu, 1978). Les systèmes SAR enregistrent
l'amplitude et la phase du signal rétrodiffusé, qui dépend des propriétés physiques et
électriques de l'objet imagé (par exemple, la rugosité du terrain, la permittivité).
Récemment, l'imagerie SAR série temporelle en bande C a été étudiée pour la
surveillance de la végétation (Dobrinić et al., 2020; Phan, 2018). Dans cette thèse, on
48
a étudié les inondations en 2019, 2020 dans le delta du Mékong à l'aide des données
Sentinel-1 - SAR de série temporelle. Des jeux de données sont construits pour deux
saisons (inondation et non-inondation) chaque année. Une carte de l’occupation du
sol va être construite avec le jeu de données Non-inondation (de décembre à mai). Et
les cartes de l’inondation vont être extraites des jeux de données de chaque mois qui
sont dans la période d’inondation de juin à novembre (Tableau 5).
Tableau 4 : Séries temporelles 2016, 2017,2018 au delta du fleuve Rouge.
49
Tableau 5 : Organisation des jeux de données SAR de série temporelle au delta du Mékong en 2019-2020.
Juin
Juillet
Chaque mois 5 images – 20
Août
inondation bandes de double polarisation VV-
Novembre
VH
Octobre
Novembre
50
Chapitre 6. Machine Learning (Apprentissage automatique)
51
6.1. Définition Apprentissage Automatique
Sur la base des mappages sous-jacents entre les données d'entrée et les
résultats attendus présentés pendant la phase d'apprentissage de l'apprentissage
automatique, les algorithmes d'apprentissage automatique peuvent être classés dans
six catégories (Awad et Khanna, 2015) : l’apprentissage supervisé, l’apprentissage
non-supervisé, l’apprentissage semi-supervisé, la méthodologie d'apprentissage par
renforcement, l’apprentissage traductif, l’inférence inductive. Les algorithmes
d'apprentissage automatique fréquemment appliqués dans les recherches de
télédétection appartiennent à deux catégories d'algorithmes : l'apprentissage
supervisé et l'apprentissage non-supervisé. L'efficacité de l'application de ces
algorithmes a été largement prouvée dans les recherches (Blaschke et al., 2014; Hay
et Castilla, 2008; Lary et al., 2016).
52
Figure 8 : Le processus général de l'algorithme d'apprentissage supervisé
53
Multilayer perceptron, Convolution Neural Network, Auto-Encoder (Zhu et al., 2017).
Dans cette thèse, le Convolutional neural network est appliqué pour détecter de
l’occupation du sol au delta du Mékong au Viêt-Nam. On a utilisé le modèle AlexNet –
Réseau de neurones convolutifs entièrement connecté (FCN - Figure 9)
Figure 9 : L'architecture globale du Convolutional neural network (CNN) comprend une couche d'entrée,
plusieurs couches de convolution et de max-pooling alternées, une couche entièrement connectée et une couche
de classification.
En 1976, Kettig et Landgrebe ont appliqué des segmentations sur des images
satellites par des extractions et des classifications des objets homogènes (Kettig et
Landgrebe, 1976). Par la suite, des techniques de traitement des images satellites ont
été développées. Il y a plusieurs travaux qui ont tenté d’appliquer des techniques de
segmentation d’images pour des applications en télédétection. Après l’arrivée des
images de haute résolution, des algorithmes de segmentation ont été développés de
manière significative (Blaschke, 2010; Schiewe, 2012). La plupart des travaux
antérieurs ont catégorisé la segmentation d'image en 3 approches : basée contour,
basée région, basée clustering.
54
La fonction de similarité, c’est un facteur important dans le processus de
l’algorithme K-means. Elle est représentée par le calcul de distance entre les objets
avec N-dimension. La distance euclidienne est souvent utilisée dans ces calculs.
La distance de Minkowski :
𝑛
𝑝
𝐷𝑖𝑠𝑡(𝑥, 𝑦) = √∑ (𝑥𝑖 − 𝑦𝑖 )𝑝 (6)
𝑖=1
La distance Euclidienne :
La distance Manhattan :
𝑛
∑ ∑ 𝐷𝑖𝑠𝑡(𝑥𝑖 , 𝜇𝑖 )2 (9)
𝑥𝑖 ∈ 𝐶𝑖
Les résultats ne sont pas uniques après chaque exécution. Cela dépend des
centroïdes initiaux. Il existe plusieurs approches pour collecter les centroïdes, à savoir
des approches de points aléatoires, des points stratifiés, des points d’incertitude, la
densité de probabilité (Catlett, 1991 ; Khan et Ahmad, 2004 ; Lewis et Catlett, 1994).
55
La méthode consiste à tracer la variation expliquée en fonction du nombre de clusters,
et à choisir « le coude » de la courbe comme le nombre de clusters à utiliser.
Avec :
56
Les avantages de K-means :
Pour que la classification par K-means soit efficace, il est nécessaire de spécifier
le nombre de clusters (au début de l’exécution de l’algorithme). Il est également difficile
de comparer la qualité des clusters produits.
6.3.2 Meanshift
57
Figure 11 : Image optique Sentinel-2 sur TranDe au 20 février 2020
La recherche d’un ensemble des noyaux multi-variés par des fonctions des
densités de probabilité Gaussien sur des informations spatiales et des signaux
(Demirović, 2019; Michel et al., 2015).
58
𝐶 𝑥𝑠 2 𝑥𝑟
𝐾(𝑥) = 𝑝 𝑘 (‖ ‖ ) 𝑘 (‖ ‖2 ) (11)
ℎ𝑠 2 ℎ𝑟 ℎ𝑠 ℎ𝑟
Avec :
Informations spatiales :
Le regroupement des pixels autour des noyaux les plus proches (maximum
probabilité).
59
constatons qu’il y a des objets qui ne changent pas, c’est-à-dire qu’ils sont
convergents. En outre, une segmentation sous forme de grands objets a été constatée
dans la Figure 14(C) et en revanche, des objets de taille plus fine sont révélés dans
les figures Figure 14(A) et Figure 14(B). A ce niveau on se pose la question : « quelle
est la meilleure bande pour la segmentation pour chaque image satellite ? »
Figure 14 : Segmentation de MeanShift avec la bande passante (A) K = (16,32), (B) K = (05,15), (C) K =
(03,05)
L’objectif de l’algorithme MRS est de fusionner des objets d’une image, de sorte
que l’objet soit précédé d'un pixel. À chaque étape, une paire d’objets est fusionnée
en un plus grand objet. La décision de la fusion dépend de la condition d’homogénéité
qui représente la similarité des objets adjacents. A son tour, la possibilité d’une fusion
des paires d’objets adjacents dépend du coût de la fusion qui doit être minimum. La
procédure est arrêtée quand il n’y a plus de possibilité de fusion dans une itération.
Le coût de fusion est composé par des valeurs homogénéités de couleur 𝐻𝑐𝑜𝑙𝑜𝑟
et de forme 𝐻𝑠ℎ𝑎𝑝𝑒 .
60
𝐻𝑑𝑖𝑓𝑓 = 𝑤𝑐𝑜𝑙𝑜𝑟 ∗ 𝐻𝐶𝑜𝑙𝑜𝑟 + (1 − 𝑤𝐶𝑜𝑙𝑜𝑟 )𝐻𝑠ℎ𝑎𝑝𝑒 (12)
L'hétérogénéité spectrale (ℎ𝑜𝑏𝑗 ) d'un objet est donnée par la somme des
écarts-types (𝜎𝑖 ) des valeurs spectrales de chaque bande i multipliée par leur poids
(𝑤𝑖 ):
𝑛𝑏𝑎𝑛𝑑𝑒
ℎ𝑜𝑏𝑗 = ∑ 𝑤𝑖 𝜎𝑖 (13)
𝑖=1
𝑙
ℎ𝑐𝑝𝑡 =
√𝑛𝑜𝑏𝑗
L'hétérogénéité douceur (ℎ𝑠𝑚𝑜𝑜𝑡ℎ ) d’un objet est composée par le nombre
des pixels (𝑛𝑜𝑏𝑗 ) dans l’objet et le périmètre (𝑏) de la boîte délimitant l’objet.
𝑙
ℎ𝑠𝑚𝑜𝑜𝑡ℎ =
𝑏
Le coût de fusion de forme (𝐻𝑠ℎ𝑎𝑝𝑒 ) est basé sur la somme des coûts de fusion
de compacité et du lissage :
61
Les inconvénients de MRS :
Les objets ne sont pas réels, il s’agit de rassemblements des pixels homogènes,
et les objets réels sont fusionnés à cause de la confusion spectrale des bandes des
images.
Figure 15 : Segmentation par multi résolution (A) Scale = 10, (B) Scale = 20 (C) Scale = 40
Parmi les différents algorithmes proposés pour générer des superpixels, SLIC
se trouve sur la liste des algorithmes montrant des performances supérieures (Stutz et
al., 2018). Il présente des avantages en termes de simplicité, de vitesse de calcul et
d'efficacité de la mémoire. C’est un algorithme basé sur le clustering inspiré d’autres
méthodes telles que k-means initialisé par des centroïdes en utilisant des informations
spectrales et spatiales. Les paramètres qui doivent être fixés pour l’exécution de SLIC
sont le paramètre K, qui indique le nombre de superpixels à générer, et le paramètre
m, qui contrôle la compacité des superpixels.
62
sur une grille carrée de taille S × S sont générés. Cette étape est construite selon
l’algorithme K-Means. Une fois que tous les emplacements des centroïdes de cluster
initiaux sont confirmés, une procédure itérative est appliquée pour mettre à jour les
centres de cluster. La procédure itérative peut être répartie en trois étapes :
Dans SLIC, les informations de l’espace couleurs RGB sont remplacées par
des informations couleurs CIELAB. La distance D est calculée par la composition de
la distance de couleur avec la distance du spatial :
La distance couleur :
La distance spatiale :
63
Plus la valeur de m est élevée, plus la proximité spatiale est accentuée et plus
la grappe est compacte.
(A) (B)
Figure 16 : Segmentations par SLIC avec la taille minimum (A) S = 16 pixels (B) S=8 pixels
2 : Les Itérations
Pour chaque cluster 𝐶𝑘 : les pixels les plus proches des centroïdes
sont attribués en fonction des distances et à partir d'un voisinage
carré 2S × 2S autour du centre du cluster
Le calcul des nouveaux centres de cluster et l'erreur résiduelle E
64
image Landsat au début des années 1970. Parmi eux, on trouve les algorithmes
d’apprentissage automatique peu profond (Machine Learning) comme la méthode de
maximum de vraisemblance (Maximum Likelihood ML), les algorithmes de réseaux
neuronaux (ex. Multilayer Perceptron MLP), des arbres de décision (ex. Random
forest-RF) et l’algorithme de Support Vector Machine (SVM). Par ailleurs, suivant la
tendance du développement des méthodes d’intelligence artificielle, plusieurs
algorithmes d’apprentissage profond sont appliqués pour la classification des images
(ex. Convolutional neuron network CNN).
Le principe de (SVM) est basé sur la division, par un hyperplan séparateur, d’un
ensemble en deux sous-ensembles marqués par des signes +1 ou -1. Il s’agit de
déterminer un hyperplan optimal de façon à ce que la distance aux exemples
d’apprentissage les plus proches (vecteurs de support) soit maximale. Cette distance
est appelée “marge”.
𝑤∙𝑥+𝑏 =0
Deux hyperplans de support :
65
𝑤 ∙ 𝑥𝑖 + 𝑏 ≥ +1 pour 𝑦𝑖 = +1
𝑤 ∙ 𝑥𝑖 + 𝑏 ≤ −1 pour 𝑦𝑖 = −1
2 2
𝑚𝑎𝑟𝑔𝑒 = =
‖𝑤‖
√∑𝑘𝑖=1 𝑤𝑖2
𝑤 ∙ 𝑥𝑖 + 𝑏 ≥ +1 − 𝜉𝑖 pour 𝑦𝑖 = +1
𝑤 ∙ 𝑥𝑖 + 𝑏 ≤ −1 + 𝜉𝑖 pour 𝑦𝑖 = −1
Avec ∀ 𝜉𝑖 ≥ 0
(a) (b)
66
Support Vector Machine non-linéaire
Généralement, les données sont des espaces non linéaires (Figure 18). Des
vecteurs supporteurs sont implémentés de façon à ce que les entrées x i soient
transformées par une fonction noyau 𝛷: 𝑅 𝑘 → 𝐻 qui projette les données dans un
espace dans lequel elles seront considérées comme linéairement séparables (Figure
18). Ainsi, un nouveau SVM linéaires est appliqué sur le nouvel espace H.
(a) (b)
Figure 18 : Hyperplan de SVM (a) données non-linéaires, (b) données transformées par noyaux
polynomiales avec 2-degrés
vote majoritaire « max-wins voting » soit en post traitement des résultats grâce à
l’estimation de probabilités a posteriori.
67
Tableau 6 : Noyaux de SVM non linéaire.
Linéaire 𝐾(𝑢, 𝑣) = 𝑢. 𝑣
Le degré d et la gamma
Polynomial 𝐾(𝑢, 𝑣) = (𝛾. 𝑢. 𝑣 + 𝑐)𝑑
𝛾 sont à préciser
2
Radial basis fonction (RBF) Gaussien 𝐾(𝑢, 𝑣) = 𝑒 −𝛾‖𝑢−𝑣‖ Gamma 𝛾 à préciser
Cette méthode est considérée plutôt pratique au-delà des arguments théoriques
ou empiriques avancés pour sa bonne performance.
Depuis les années 1980, les algorithmes d'arbre de décision ID3, CART, C45
ont prouvé leur efficacité en utilisant des données numériques et symboliques
(Quinlan, 1986; Salzberg, 1994). Ces algorithmes donnent des résultats de bonne
précision pour les données combinées multi-sources et hétérogènes. Afin de
performer ces algorithmes d’arbre de décision et d’éliminer certains de leurs
inconvénients (Breiman, 2001) a développé un algorithme appelé Random forest
68
construit à partir d'arbres de décision générés de sous-échantillons aléatoires en
utilisant l’approche du Bagging.
L'arbre de décision est une structure hiérarchique générée par des fonctions
statistiques et probabilistes. Le but du processus de construction d'arbres de décision
est d'utiliser des règles de division et d'arrêt de calcul qui peuvent être partitionnées
en ensembles d'éléments homogènes. L’algorithme de référence est celui de
(Breiman, 1998) décrivant la méthode CART (Classification et Regression Tree).
Néanmoins, la paternité de cette méthode est attribuée à (Morgan et Sonquist, 1963)
Le principe des arbres de décision repose sur la prédiction d’une classe en sortie
après un processus de séparation successive des données. L’objectif est de prédire
des groupes les plus homogènes possibles par rapport à la variable à prédire. La mise
en place de cette méthode est considérée comme étant facile et rapide, mais elle peut
être sensible à l’apprentissage exagéré.
69
Echantillons
E1 E2 ……… Ek
Y1 Y2 ……… Yk=
Un élément xi est appliqué sur tous les arbres dans la forêt, et la prévision
de chaque arbre est sauvegardée.
70
automatique supervisé qui aide à résoudre les problèmes de la classification binaire.
Il est proposé par (Rosenblatt, 1957).
Dans le neurone de computation (Figure 21), le neurone reçoit des signaux qui
sont transmis par des dendrites (des entrées). Ce neurone intègre des signaux par
une opération de multiplication des valeurs d’entrées X avec une matrice de poids W.
Ensuite, une fonction d’activation f est appliquée sur la valeur d’intégration avant de
l’envoyer au neurone suivant.
71
Figure 21: Un neurone artificiel
RELU :
𝑓(𝑥) = (0, 𝑥)
Sigmoïde :
1
𝑓(𝑥) =
1 + 𝑒 −𝑥
Tanh :
𝑒 𝑥 − 𝑒 −𝑥
𝑓(𝑥) =
𝑒 𝑥 + 𝑒 −𝑥
Les valeurs sont transférées de couche à couche (Figure 24). Si 𝑋 (𝑙−1) sont des
sorties de couche (𝑙 − 1), les sorties 𝑋 (𝑙) de couche (𝑙) sont calculées par l’application
d’une fonction d’activité f sur la somme des valeurs de sorties de la couche (𝑙 − 1):
73
Dans la dernière couche, il existe une fonction qui permet d’évaluer la qualité
du modèle. La fonction est souvent appelée loss-function ou cost-function, et la
valeur calculée par cette fonction est appelée "loss" ou "cost". L’objectif du loss-
function est de minimiser l'erreur de prévision : Mean Square Error, Cross Entropy
Loss.
Le Mean Square Error, MSE, est calculé par la moyenne des différences au
carré des valeurs prévues et réelles. Sa valeur est importante quand le modèle
présente plus d’erreurs. Son meilleur résultat est égal à zéro (0). Le Mean Square
Error (MSE) est donné par :
𝑛
1
𝐶𝑜𝑠𝑡 = ∑ ‖𝑦𝑖 − 𝑦̂𝑖 ‖2
𝑛
𝑖=1
74
parmi les algorithmes d’apprentissage d’un réseau de neurones. Cet algorithme se
base sur les optimisations des valeurs de loss-function (Cost) en utilisant des
𝜕𝐶𝑜𝑠𝑡
dérivées partielles
𝜕𝑤
de la fonction Cost par rapport à chaque poids 𝑤. Les
(𝐿)
𝛻𝑋 𝐶𝑜𝑠𝑡 = {𝛿𝑖 }
Les erreurs de couche de sortie (𝐿) sont calculées avec l’opération Hadamard
product ⨀ données par :
75
𝛿 (𝑙) = (𝑊 (𝑙+1) 𝛿 (𝑙+1) )⨀ 𝑓 ′ (𝑋 (𝑙) )
La mise à jour des poids w de chaque couche (𝑙) dans l’itération 𝑡 est donnée
par :
(𝑙) (𝑙) (𝑙)
𝑊𝑡 = 𝑊𝑡−1 − 𝜂 . 𝛿𝑡
En résumé, l’algorithme backpropagation se présente comme suit :
76
(pool-layer) est la fusion des caractéristiques dans une fenêtre MxM en un pixel par
des opérations maximales ou moyennes. Après une couche de regroupement, les
feature-maps sont réduites. Ses bases caractéristiques sont conservées pour l’étape
suivante.
Figure 26 : Modèle CNN avec 02 couches convolutives, 02 couches de regroupement, une couche
perceptron. Les images de l’entrée ont la taille 8x8, et la sortie a 05 classes de classification.
Dans une couche convolutive (Figure 27, Figure 28), les images ou les feature-
maps sont convolutées avec un ensemble K noyaux 𝑊 = {𝑊1 , 𝑊2 , … , 𝑊𝐾 } et 𝛾 =
{𝑏1 , 𝑏2 , … , 𝑏𝑘 }. La dimension de chaque noyau est k. Le feature-map de couche (𝑙) est
calculé par une opération convolutive ⨂ avec le noyau k qui est donnée par :
(𝑙) (𝑙−1) (𝑙−1)
𝑋𝑘 = 𝑓(𝑊𝑘 ⨂ 𝑋 (𝑙−1) + 𝑏𝑘 )
Dans la Figure 26, l’image d’origine est une image de Sentinel-2, la mission de
classification est la détection des occupations du sol de la province Trande au delta
du Mekong. Le modèle CNN est installé par plusieurs couches convolutives et couches
de regroupement. Dans la première couche convolutive, il y a 64 noyaux de 3x3
dimensions. Six feature-maps sont sélectionnées avec des caractéristiques très
particulières. Les caractéristiques sont utilisées pour classifier des classes : eau (a),
bâtiment (c), des végétations (d), des non végétations (f), etc.
77
quartier rectangulaire, la norme 𝐿2 d'un quartier rectangulaire ou une moyenne
pondérée basée sur la distance du pixel central.
A partir des travaux de (Olofsson et al., 2014, 2013) l’objectif de ce chapitre est
de présenter une méthode simple de l’utilisation des informations extraites d’un
processus d’apprentissage automatique. Ceci dans une optique d’évaluation de la
qualité de la carte de l’occupation du sol et d’estimation des superficies des classes
ou du changement de l’occupation du sol pour construire des intervalles de confiance
qui reflètent l’incertitude de cette estimation de superficies obtenues.
𝑛𝑖𝑖
𝑃𝑖 = (16)
𝑛.𝑖
La précision globale – OAA (OverAll Accuracy)
∑𝑘𝑖=1 𝑛𝑖𝑖
𝑂= (17)
𝑁
Tableau 7 : Matrice de confusion de k classes.
Classes 1 2 3 ⋯ 𝑘 total
⋮ ⋮ ⋮ ⋮ ⋱ 𝑛1𝑘 𝑛1.
𝑛𝑖𝑗
𝑝𝑖𝑗 = 𝑊𝑖 (18)
𝑛𝑖.
⋮ ⋮ ⋮ ⋮ ⋱ ⋮ ⋮ ⋮ ⋮
80
Le Tableau 9 est construit avec des proportions 𝑝𝑖𝑗 . Les précisions et ses
estimations sont calculées directement par ce tableau.
Tableau 9 : Matrice de confusion estimée.
⋮ ⋮ ⋮ ⋮ ⋱ ⋮ ⋮
𝑝𝑖𝑖
̂𝑖 =
𝑈 (19)
𝑝𝑖.
Estimation des précisions du producteur (Producer Accuracy) :
𝑝𝑖𝑖
𝑃̂𝑖 = (20)
𝑝.𝑖
𝑂̂ = ∑ 𝑝𝑖𝑖 (21)
𝑖=1
̂:
Estimation des variances pour les précisions de l’utilisateur 𝑈
̂𝑖 (1 − 𝑈
𝑈 ̂𝑖 )
𝑉̂ (𝑈
̂𝑖 ) = (22)
(𝑛𝑖. − 1)
Estimation des variances pour les précisions du producteur 𝑃̂:
81
1 𝑛.𝑗 (1 − 𝑃̂𝑗 )𝑈
̂𝑗 (1 − 𝑈
̂𝑗 )
𝑉̂ (𝑃̂𝑗 ) = 2[
̂.𝑗
𝑁 (𝑛𝑗. − 1)
𝑛𝑖𝑗 𝑛𝑖𝑗
𝑘 𝑛.𝑖 2 𝑛 (1 − ) (23)
𝑛𝑖.
+ 𝑃̂𝑗 ∑ 𝑖.
]
𝑛𝑖. − 1
𝑖≠𝑗
𝑛.𝑖
̂.𝑗 = ∑𝑘𝑖=1
Avec 𝑁 𝑛𝑖𝑗
𝑛𝑖.
𝑘
𝑝𝑖𝑗 (1 − 𝑝𝑖𝑗 )
𝑆̂(𝑝.𝑗 ) = √∑ 𝑊𝑖 2 = (27)
(𝑛.𝑗 − 1)
𝑖=1
Dans cet exemple, une carte d’occupation du sol (Figure 29) est générée par
une application de la méthode CNN – Convolutional Neuron Network pour détecter
cinq classes d’occupation du sol qui sont : l’eau permanente, les rizières, les
bâtiments, les vergers, les forêts. Des données de validation ont été générées sous
forme de polygones, ou région d’intérêt (ROI), depuis une carte de référence. Le
processus de l’estimation des précisions commence en construisant le tableau de
confusion étendue (Tableau 10):
82
Tableau 10 : Matrice de confusion de la carte d'occupation du sol au delta du Mékong en 2019.
Superficie
Classes 1 2 3 4 5 Pixels Poids (Wi)
(km2)
Précision de Précision du
Nom de Classe Code de classe
l'utilisateur producteur
Figure 29 : Carte de l'occupation du sol de province An Giang et province Dong Thap au delta du
Mékong en 2019.
83
Ainsi à l’étape suivante, on calcule la matrice de confusion par les formules (18),
(27) (Tableau 12)
Classes 1 2 3 4 5 𝑝𝑖.
Dans la dernière étape, des précisions estimées sont calculées par des formules
(19)(20)(21)(22)(23)(24)(26). Ainsi, on obtient le tableau d’évaluation des précisions
de la carte de l’occupation du sol dans la zone d’étude (deux provinces AnGiang et
DongThap au delta du Mékong) Tableau 13.
84
Tableau 13 : Estimations des précisions de la carte de l'occupation du sol.
Des comparaisons pixel par pixel entre les deux cartes sont présentées dans le
Tableau 14. Des valeurs de (𝑛11 , 𝑛22 , 𝑛33 , 𝑛44 , 𝑛55 ) sont des superficies sans
changement. Celle 𝑛𝑖𝑗 représente le nombre des pixels qui change de 𝑐𝑙𝑎𝑠𝑠𝑒 𝑖 en 2019
à 𝑐𝑙𝑎𝑠𝑠𝑒 𝑗 en 2020. Ici dans cet exemple, des superficies qui changent sont évaluées
par les valeurs estimées dans le Tableau 15 qui donne les proportions du
changements de chaque classe.
Classe 1 2 3 4 5 Total
85
Dans le Tableau 15, la précision globale 𝑂𝐴𝐴 = 69 ± 0.01 signifie que le
pourcentage des superficies qui ne changent pas est de l’ordre de 69% avec une
variance égale à 0.01 %. A la ligne 3, la précision de l’utilisateur de la classe bâtiments,
UA = 45.07 ± 0.04 % montre la proportion des changements (1 − 𝑈𝐴) qui est égale à
55 %. La classe bâtiments est la classe qui a connu plus de changements par rapport
aux autres classes.
Tableau 15 : Précisions estimées de la classification par CNN.
86
Chapitre 7. L’Analyse orientée objet d'images satellite
Figure 31 : Relation entre les objets et la résolution spatiale. (a) basse résolution : pixels nettement plus
grands que les objets (b) résolution moyenne : les tailles des pixels et des objets sont du même ordre, les
techniques pixel par pixel sont appropriées.
Dans leur travaux Hay et Castilla (2008) ont défini le concept comme suit :
« L'analyse d'images basée sur les objets géographiques (GEOBIA) est une sous-
discipline de la science de l'information géographique (GEOScience) consacrée au
développement de méthodes automatisées pour partitionner l'imagerie de
télédétection en objets d'image significatifs, et évaluer leurs caractéristiques à travers
les échelles spatiales, spectrales et temporelles, afin de générer de nouvelles
informations géographiques dans un format prêt pour le SIG. »
87
7.3. Evaluation de GEOBIA
Les évaluations de GEOBIA sont données par Hay et Castilla (2008), elles sont
poursuivies ensuite (Blaschke, 2010; Hay et al., 1996). Grâce à ces travaux, on
connaît les avantages et les inconvénients de GEOBIA.
Certains des algorithmes d’extraction des objets offrent des options trop
compliquées, ce qui oblige l'analyste à prendre le temps de les régler.
Les facteurs d’un paradigme de GEOBIA sont : des données, des méthodes de
segmentations, des méthodes de classifications (Baatz et al., 2008; Blaschke et al.,
2014; Castilla et Hay, 2008).
88
7.4.1 Etape de travail dans le cadre de GEOBIA
Les données utilisées dans GEOBIA comprennent des images satellites qui sont
des données entrées. Et elles sont enregistrées sous forme de fichiers RASTER. Les
pixels dans une image sont organisés en 2-dimensions. Chaque pixel est un vecteur
K-valeurs où K est le nombre de bande dans l’image optique. Dans le cas des images
SAR il y a un nombre doubles de polarisations. Et dans la troisième figure de cas, ce
sont les images fusionnées où K est le total des bandes. La structure des données
satellites est présentée dans les sections Chapitre 5, Chapitre 6.
Les objets sont des résultats de l’étape de segmentation (Figure 32) et sont
sauvegardés dans un fichier GIS sous forme des polygones et ses coordonnées. On
a besoin d’une sous-étape pour extraire les caractéristiques (les moyennes, les écarts-
types, la superficie, ...). Les échantillons entraînés et les échantillons validés sont
collectés à partir des objets extraits dans l'étape de segmentation avec des
informations d’observation.
Méthodes non-supervisées
89
fenêtre m x m ou des algorithmes d’apprentissage automatique non-supervisés. Dans
cette thèse, les méthodes de segmentations sont présentées dans la section 6.3. Les
algorithmes de segmentation donnent plusieurs paramètres. Et les optimisations des
paramètres améliorent l'homogénéité intra-segment, l'hétérogénéité inter-segment et
la précision de la classification (Gao et al., 2011; Yang et al., 2015)
Méthodes supervisées
90
Partie 3 : De l’acquisition des données aux
Applications
91
Chapitre 8. Détection des changements de l’occupation du sol
dans le delta du fleuve Rouge (Viêt-Nam)
92
Figure 34 : Delta du fleuve Rouge au Viêt-Nam.
Dans les zones résidentielles, situées sur des terrains plus élevés et hors de
portée des inondations, l’occupation des terres est caractérisée par la présence des
villages ou l’on trouve des vergers saisonniers et des arbres d'ombrage comme le
bananier, l'eucalyptus et le filao, ainsi que quelques arbres fruitiers. Dans certaines
zones de haute altitude, le filao est cultivé pour la production de bois et pour la
protection des villages de l'intérieur contre les tempêtes. Dans la région côtière,
l’occupation des terres est principalement de type agricole, avec le riz comme
principale culture. Dans certaines zones très basses, salées et récemment défrichées,
le scirpe est temporairement cultivé pour améliorer la qualité du sol avant de laisser
place aux cultures de riz. Au contraire, dans les zones côtières salines et les estuaires,
la végétation se développe naturellement pour former des forêts. Cependant, la
formation végétale prédominante est la forêt de palétuviers, principalement du type
Sonneratia, avec une hauteur allant de 1 à 3,5 m partie importante de cette zone est
actuellement couverte par des étangs d'aquaculture pour la culture de crevettes et de
crabes, utilisant des méthodes industrielles aussi bien que traditionnelles (Niculescu
et Lam, 2019).
93
8.2. Données mises en œuvre
Pour cette application, nous avons utilisé les données du satellite optique
Sentinel-2 et les données optiques de Landsat-8 (Tableau 16). La principale source
de données provient des satellites Sentinel-1. Toutes les images Sentinel-2 (niveau
1C) couvrant l'ensemble de la zone d'étude ont été téléchargées depuis le Sentinel
Scientific Hub de l'ESA (ESA, 2018). Ce satellite est capable de fournir des images
avec une largeur au sol de 290 km, et une périodicité de 5 jours. Les acquisitions sont
réalisées selon trois résolutions en fonction de la bande spectrale : 10, 20 et 60 m
(Tableau 17).
Par ailleurs des images Landsat ont été téléchargées depuis le site officiel
Landsat [Link] Le satellite Landsat-8 transporte une charge
utile à deux capteurs, l'imageur terrestre opérationnel (OLI) et le capteur infrarouge
thermique (TIRS). L'OLI dispose de deux bandes supplémentaires de longueurs
d'onde réfléchissantes : une nouvelle bande bleue à courte longueur d'onde (0,43 -
0,45 μm) destinée à améliorer la sensibilité à la chlorophylle et aux autres matières en
suspension dans les eaux côtières et à récupérer les propriétés des aérosols
atmosphériques, et une nouvelle bande infrarouge à ondes courtes (1,36 - 1,39 μm)
pour la détection des cirrus. Les autres bandes OLI sont spectralement plus étroites
dans la plupart des cas que les bandes ETM+ correspondantes (Roy et al., 2014)
(Tableau 18).
Les dates de prise de vue des images satellites disponibles pour la zone
d’étude, ont été choisies en fonction des principales étapes de la culture du riz dans le
delta du fleuve Rouge.
94
Tableau 16 : Ensemble de données satellitaires disponibles.
1 05/26/2018 Sentinel-2
2 06/20/2018 Sentinel-2
3 07/05/2018 Sentinel-2
4 11/02/2018 Sentinel-2
5 12/17/2017 Sentinel-2
6 09/17/2017 Landsat-8
7 07/30/2017 Sentinel-2
8 04/21/2017 Sentinel-2
9 06/15/2016 Sentinel-2
10 04/23/2016 Landsat-8
11 06/15/2016 Sentinel-2
12 04/23/2016 Landsat-8
Bande passante
Sentinel-2 Bands longueur d’onde (µm) Résolution (m)
(nm)
95
Tableau 18 : Résolutions spatiales et spectrales du satellite Landsat-8.
(micromètres) (mètres)
8.3. Méthodologie
GEOBIA est un paradigme évolutif avec des outils, des logiciels, des méthodes,
des règles et un langage spécifique et il est de plus en plus utilisé pour les études qui
ont besoin de conceptualiser et de formaliser les connaissances représentant la réalité
basée sur l’emplacement (Blaschke et al., 2014). Le processus de classification de
cette méthode utilise un processus de segmentation hiérarchique au stade du
prétraitement afin que les caractéristiques de l'objet puissent être complétées par des
informations supplémentaires telles que la texture, le contexte et d'autres informations
liées à l'objet.
96
Les indices de végétation ont été calculés pour l'image de réflectance Sentinel-
2A, il est défini par la réflectance de la bande rouge (RED) et de la bande proche
infrarouge (NIR) car elles détectent des profondeurs très différentes à travers les
couverts végétaux. La bande Rouge se situe dans la région de forte absorption de la
chlorophylle, tandis que la bande NIR présente une réflectance élevée du couvert
végétal dans cette zone.
Saison de
Période culture croissance récolte
culture
Une fusion ou empilement des couches a été effectué en ajoutant toutes les
couches d'images optiques traitées, et ce après une normalisation des couches à la
même plage de valeurs. Les couches d'images suivantes ont été empilées :
Tableau 20 : Calendrier des cultures de riz et passages des Satellites (bleu Landsat, jaune Sentinel -2).
L'ESP permet une estimation rapide des paramètres d'échelle pour une
segmentation multirésolution dans l'environnement logiciel Definiens. ESP, analyse la
variance locale (LV) de l'hétérogénéité des objets des images satellite. En nous
98
appuyant sur les résultats de Kim et al. (2008), nous avons utilisé le concept de LV
(Woodcock et Strahler, 1987) pour créer un outil qui informe sur la sélection des
paramètres d'échelle appropriés pour la segmentation dans Definiens Developer.
99
8.3.3 Classification des images
Random Forest
100
Multi-Layer Perceptron (MLP)
L’ensemble de neurones de la même couche est sans connexion entre eux. Les
neurones sont reliés entre eux, de couche en couche, par des connexions pondérées.
Ce sont les poids de ces connexions qui régissent le fonctionnement du réseau et qui
« programment » une application de l'espace des entrées à l'espace des sorties au
moyen d'une transformation non linéaire.
Pour notre application, nous avons opté pour un modèle (Figure 37) avec 20
entrées de bandes, 32 neurones sur la première couche cachée, 64 neurones sur la
deuxième, 32 neurones sur la troisième et 5 classes sur la couche de sortie. Le
processus d'apprentissage passe par deux étapes : l'anticipation et la rétro-
propagation.
101
8.4. Résultats et discussion
Pour les résultats de l’algorithme Random Forest, nous constatons qu’en ce qui
concerne la distribution des classes (type) d'utilisation des terres dans le district, les
champs agricoles (rizières) situés à l'intérieur des digues couvrent plus de 40% de la
surface totale pendant les trois années d'observation. L'eau couvre environ de 23 % à
28 % de la surface et est principalement représentée par la mer, le fleuve Rouge et
les rivières Lân et Tra Ly. Les villages, y compris les jardins, couvrent environ 14 % à
15 % des espaces et sont situés au milieu des rizières. Les sols humides ne se
trouvent pas vraiment à l'intérieur des digues, sauf au bord des canaux d'irrigation,
mais plutôt à l'extérieur des digues, aux embouchures de la rivière Tra Ly au nord et
du fleuve Rouge au sud du district. Les mangroves sont identifiées sur les îles dans
les estuaires de ces mêmes rivières au cours de 2016, mais elles sont devenues très
clairsemées en 2018. En 2018, de nouvelles mangroves ont été plantées autour des
bassins d'aquaculture. Les arbustes immergés sont dispersés à l'extérieur des digues,
parfois à l'intérieur, couvrent environ 1 à 3 % de la surface. En outre, les marais salants
couvrent très peu de territoire, principalement au centre du district et à côté des digues.
D'autre part, les étangs d'aquaculture sont souvent situés à l'extérieur des
digues, mais parfois à l'intérieur, et le pourcentage de surface qu'ils occupent n'a pas
changé en 2016 mais a augmenté d'environ 5 % en 2017 et 8 % en 2018. Les légumes
en général se trouvent souvent dans les basses terres près des rives du fleuve Rouge.
Les résultats montrent une très bonne performance de classification pour toutes
les combinaisons (indice de kappa 0,90 et 0,97). En suivant la procédure standard de
calcul de la précision (Congalton et Green, 2019), nous avons estimé une précision
globale de 93 % (2016), 97 % (2017) et 98 % (2018) et un coefficient de Kappa de
0.90 (2016) et 0,97 (2017 et 2018). Les erreurs de classification par omission et
commission ont été constatées principalement dans la classe « légumes » pour 2016
et dans la classe « zones humides » pour les trois années (Tableau 21).
102
Tableau 21 : Précision globale et indice de kappa par classe pour les classifications par forêts aléatoires.
Rizière 99 100 97
Aquaculture 95 95 99
Végétation basse 93 91 47
Zone humide 76 81 71
Par ailleurs, nous constatons que les zones humides sont confondues avec
l'aquaculture. En Asie notamment, les bassins d'aquaculture sont dominants à
proximité des deltas et des estuaires des fleuves où la culture du riz est également
une pratique importante. L'inondation temporaire des champs pendant la saison de
culture des rizières et les caractéristiques de forme similaires aux étangs d'aquaculture
peuvent créer des erreurs de classification importantes entre ces types lorsque des
approches uni-temporelles sont appliquées. Une série chronologique dense
d'acquisitions Sentinel-2 et Landsat-8 peut aider à distinguer les plans d'eau
temporairement inondés des étangs d'aquaculture inondés relativement stables.
En examinant les changements par les images satellites (Figure 38) dans les
municipalités du district de Tien Hai nous constatons que ces municipalités présentant
les mêmes conditions biophysiques (telles que Nam Phu et Dong Long) ont affiché
plusieurs tendances similaires en matière de changement d'utilisation des terres, du
fait qu'elles possèdent les mêmes écosystèmes côtiers. De plus, leurs différences sont
souvent causées par les activités humaines, à savoir, la fondation récente de Nam
Phu contribuant à l'expansion de ses champs agricoles, et la création hâtive de
l'aquaculture de Nam Phu encourageant le développement d'étangs et la coupe
importante des mangroves. Par ailleurs, les communes ayant des conditions
biophysiques différentes (comme Tây Tiên et les communes de Dong Long) n'ont pas
les mêmes types d'occupation des sols. En effet, les facteurs physiques et
anthropiques n'ont pas le même impact sur les changements d'utilisation des terres.
103
Nous constatons que, malgré l'enchevêtrement des traitements d'images qui
doivent être améliorés, la détection des changements d'utilisation des terres à l'aide
des images satellites est très utile. Nous avons pu voir comment chaque type
d'utilisation des terres a changé dans le district au cours des années observées afin
d'expliquer la relation entre la libéralisation économique et les changements
d'utilisation des terres. Ces changements reflètent l'impact des facteurs physiques (tels
que les caractéristiques des écosystèmes côtiers qui favorisent les mangroves, les
arbustes immergés, les marais salants et la production aquacole) et des facteurs
anthropiques (la libéralisation économique, la nouvelle politique de stimulation de
l'aquaculture et les plans de développement économique du gouvernement) sur
l'utilisation des terres. Les facteurs anthropiques ont un impact plus important sur les
changements que les facteurs physiques, notamment dans l'augmentation rapide des
étangs d'aquaculture, qui a entraîné à son tour des changements dans les autres types
d'utilisation des terres.
104
8.5. Conclusions
105
Chapitre 9. Étude et cartographie de l’inondation au niveau du
delta du Mékong
Dans le delta du Mékong, les terres touchées par les inondations sont les
rizières. Des études antérieures ont montré que l’identification des rizières touchées
par l’inondation n'est pas une tâche facile (Phan et al., 2018), car il faut distinguer les
eaux permanentes des eaux de crue. Le delta du Mékong est si important sur le plan
écologique, économique et social que de nombreuses études ont été menées au fil
des ans à l'aide de données d'observation de la terre pour caractériser la culture du
riz. La majorité de ces études se sont concentrées sur la délimitation de la distribution
des cultures de riz en utilisant des images satellites optiques (Kontgis et al., 2015;
Morinaka et al., 2006) et des images satellites radar (Liew et al., 1998 ; Karila et al.,
2014 ; Nguyen et al., 2016 ; Kontgis et al., 2017). Ces dernières années, différentes
approches ont été adoptées, notamment le seuillage basé sur l'analyse de la
phénologie du riz et l'apprentissage automatique.
9.1. Méthodologie
106
examinés en matière de leurs performances sur deux groupes d'occupation du sol :
les zones inondées et les zones inondables (rizières).
Dans l'imagerie SAR, l'eau apparaît généralement comme une surface lisse
avec une intensité de rétrodiffusion bien définie. La signature de rétrodiffusion permet
de détecter les plans d'eau. Ces derniers constituent un ‘’réflecteur‘’ spéculaire de
l'impulsion radar, qui entraîne un signal minimal renvoyé au satellite (Jung et al., 2010
; Schlaffer et al., 2015). La rugosité de la surface de l'eau, créée par de fortes pluies
ou par le vent, peut provoquer une rétrodiffusion du signal radar, ce qui augmente la
possibilité que les zones inondées ne soient pas mises en évidence (Alsdorf et al.,
2007 ; Jung et al., 2010). La rétrodiffusion du signal SAR varie également en fonction
de l'angle d'incidence (IA) (donnée en général par les fournisseurs de production) et
des variations de l'angle d'incidence local (LIA) dues à la topographie de la cible et à
l'angle d'incidence (Wilusz et al., 2017). L'intensité de la rétrodiffusion peut être
influencée par les conditions environnementales telles que la topographie du paysage
et les ombres. L'identification de l’extension de l’eau peut être difficile dans les zones
107
où d'autres objets dépassent et interagissent avec le signal radar. Ainsi, il est difficile
de déterminer un seuil général de rétrodiffusion.
En effet, l'environnement peut jouer un rôle très important. Par exemple, l'eau
peut être masquée par la végétation. La couverture végétale, lotus et herbes
aquatiques, peut être une source d'incertitude dans la cartographie de l'étendue de
l'inondation. Selon certains auteurs (Bouvet & Le Toan, 2011), il y a normalement de
grandes zones d'herbes aquatiques et de lacs de lotus répandus pendant la saison
des inondations. Les types de cultures dominants dans les zones d'étude de la
province d'An Giang et de Dong Thap sont le riz et les arbres fruitiers à 2 et 3 saisons.
Les données SAR ont la réponse spécifique de rétrodiffusion radar de la végétation
structurée verticalement sur un sol inondé ou humide et donc la capacité de distinguer
le riz des autres classes de couverture du sol. L'interaction entre une onde
électromagnétique radar et la végétation implique principalement trois mécanismes :
la réflexion directe (spéculaire), le double rebond, et la diffusion volumique (Figure 40).
Le riz est la principale culture dans notre zone d'étude et les rizières peuvent
être affectées par des événements d'inondation. La cartographie de l'étendue des
inondations est un processus qui permet d'identifier les zones terrestres impactées par
la flooding. Dans le même temps, la production de cartes d'utilisation des terres /
couverture des terres (LULC) est devenue un moyen fréquent et répandu et apparaît
comme un indicateur clé pour le suivi des risques d'inondation (An Thi Ngoc Dang &
Lalit Kumar, 2017 ; Mojaddadi et al., 2017). De ce fait, elles peuvent être intégrées
dans les bases de données sur les inondations pour réaliser un zonage des lieux à
risque et déterminer la vulnérabilité au risque d’inondation.
Figure 40 : Représentation schématique des mécanismes de rétrodiffusion dans la bande C pour le riz
(Oryza sativa) avec trois phases principales dans le delta du Mékong (phase végétative, phase reproductive et
phase de maturation) : (a) réflexion directe (spéculaire) ; (b) double rebond ; et (c) diffusion volumique.
108
Au cours de la première période (phase végétative), les rizières sont
recouvertes d'eau pendant la majeure partie de leur cycle de croissance. Les champs
ne sont inondés que pendant certaines périodes (émergence et tallage de la phase
végétative et booting, épiaison pour la phase reproductive). Le sol reste humide (non
inondé) pendant le début du tallage et l'élongation et l'initiation de la panicule pour la
phase reproductive et pendant la phase de maturation (maturation). La rétrodiffusion
radar des champs inondés est faible en raison de la réflectance spéculaire de la
surface de l'eau.
Le Mékong est le septième plus long fleuve d'Asie et le douzième plus long du
monde. Son delta couvre une superficie de près de 40 000 Km2 dans 12 provinces et
municipalités habitées par près de 18 millions de personnes, soit 20% de la population
du pays.
109
Figure 41 : Zone d'etudes - Province AnGiang et Province DongThap au delta du Mékong.
Le delta du Mékong est une vaste plaine inondable et l'une des zones agricoles
les plus fertiles d'Asie. En effet, le sud du Viêt-Nam est une grande région agricole, qui
est spécialisée dans la culture du riz. En 2020, malgré la sécheresse et la salinisation,
le delta du Mékong, grenier du Viêt-Nam, a produit 24 millions de tonnes de riz à partir
de ses 1,5 million d'hectares de terres cultivées, soit un rendement moyen par hectare
de 6 tonnes. En 2020, les exportations de riz du delta du Mékong s'élevaient à 6
millions de tonnes. (Le Courrier du Viêt-Nam, 27/05/2021). Le delta du Mékong
possède des écosystèmes rizicoles principalement irrigués avec trois saisons de
culture principales : saison de culture hiver-printemps, été-automne et automne-hiver
(Thi Hoa Phan, 2018). Par rapport à la saison des pluies et des inondations, nous nous
intéressons davantage à la saison culturale automne-hiver. Selon le calendrier de la
riziculture de cette saison, le riz est planté en juillet-août, et récolté en décembre-
janvier. Selon Thi Hoa Phan, 2018, la culture d'automne-hiver pendant la saison des
pluies a la productivité la plus faible.
110
temps et dans l'espace. D'un mois à l'autre, la quantité de pluie enregistrée peut
doubler.
Le bas du delta du Mékong, avec ses rivières, est influencé par la mousson
d'Asie du Sud-Est. Le régime hydrologique présente de très grandes variations
annuelles et interannuelles. La mousson, avec théoriquement une saison sèche et une
saison des pluies (maximum en juillet) donne lieu à des précipitations à la mi-avril, s'il
n'y a pas de dégel, le tout avec de fortes variations d'une année sur l'autre en fonction
de la force de la mousson humide. Ensuite, les précipitations de convergence
prolongent la saison des pluies ; plus constantes, elles ajoutent des inondations
d'octobre à novembre. Et enfin, les typhons, qui génèrent le potentiel hydrique maximal
en novembre-décembre. Les inondations les plus fréquentes dans le delta sont
principalement induites par le régime des crues du Mékong ainsi que par l'écoulement
de surface. La construction de systèmes de digues à partir du début des années 1990
pour protéger les champs des inondations et permettre une troisième culture dans
certaines parties du delta, contribue à des perturbations du système hydrologique. Ce
n'est pas le cas partout dans notre zone d'étude.
111
Les deux satellites Sentinel-1 sont placés sur la même orbite à environ 900 km
de la surface de la Terre. Leur déphasage permet une meilleure résolution temporelle.
Ces images en bande C et de longueur d'onde d'environ 4,5 cm sont acquises en deux
polarisations VH et VV. La polarisation est l'orientation des champs
électromagnétiques de la longueur d'onde.
Sentinel-1 génère pour chaque acquisition deux bandes avec une polarisation
verticale en émission et horizontale en réception, et l'autre verticale en émission et en
réception. La polarisation de Sentinel-1 est également appelée polarisation
monochromatique, ce qui signifie que nous aurons toujours les mêmes amplitudes et
les mêmes longueurs d'onde. Les données sont acquises sur l'orbite N° 18 qui passe
au-dessus du delta du Mékong. Le numéro de l'orbite est un paramètre important car
il permet d'avoir des images sur la même zone d'étude et au même moment.
Cependant, les orbites ascendantes ou descendantes permettent l'acquisition. La
direction de l'orbite est un paramètre important à prendre en compte pour une
meilleure interprétation des images.
Afin de comparer les valeurs des images radar avec les précipitations, les
données quotidiennes de quatre stations météorologiques de 2019 à 2020, notamment
celles relatives à la précipitation, ont été téléchargées depuis la plateforme NOAA :
[Link] et le niveau d'eau du Mékong depuis la Mekong River
Commision [Link]
112
9.5. Méthode d'apprentissage automatique et classification
Les CCNs sont devenus populaires en raison de leur capacité à résoudre des
problèmes de classification tels que la reconnaissance d'images et la classification de
séries temporelles. LeCunet al., 1989 ont obtenu de très bons résultats en proposant
un réseau de neurones convolutifs à partir d'un modèle appliquant la rétropropagation.
L'idée de développer des réseaux de neurones convolutifs était initialement basée sur
la connectivité locale, dans laquelle chaque nœud est connecté uniquement à une
région locale de l'entrée (Borovykh et al., 2017). Le réseau résultant possède de
nombreuses connexions mais relativement peu de paramètres libres.
113
9.5.2 Construction du modèle du CNN
L'architecture des CNNs est une structure comportant une série de couches :
couche convolutive, couche de regroupement, couche de perceptron. Les réseaux de
neurones convolutifs (CNNs) que nous proposons sont constitués de deux couches
convolutives et deux couches de regroupement (Max pooling layer) (Figure 42).
Le processus de convolution des matrices est répété sur les cartes neuronales
(neural maps) construites pour créer de nouvelles couches cachées profondes. La
couche de regroupement (pool-layer) a pour rôle de fusionner les caractéristiques de
la fenêtre en un pixel par des opérations de maximum ou de moyenne.
D'après Kim, 2017, les réseaux neuronaux convolutifs peuvent être considérés
comme des réseaux neuronaux artificiels multicouches feedforward entraînables, qui
comprennent plusieurs étapes d'extraction de caractéristiques. Les couches
convolutives caractérisent chaque étape d'extraction de caractéristiques avec des
114
filtres d'apprentissage, des couches de mise en commun et des fonctions d'activation
ou des couches de non-linéarité (Gebrehiwot, 2019). Un autre élément important dans
l'ensemble du processus est l'unité linéaire rectifiée ou ReLU, qui remplace les valeurs
négatives du réseau en cours de formation par un zéro.
L'algorithme Random Forest (RF) a été appliqué pour la classification des zones
d'inondation et des zones inondables, en utilisant le logiciel libre Orfeo ToolBox (OTB),
dans sa version intégrée à QGIS. Ce classificateur construit par Breiman (2001) est
un algorithme non paramétrique d'apprentissage automatique basé sur les arbres de
décision. Il est fréquemment utilisé pour la classification d'images SAR, que ce soit en
zone forestière (Balzter et al. 2015 ; Rüetschi et al., 2018) ou en zone de culture
(Ghazaryan et al., 2018). Il a l'avantage de ne pas être influencé par le problème du
sur-apprentissage qui se produit lorsqu'un modèle est trop basé sur les échantillons
d'entraînement.
115
Plusieurs études fournissent également un intervalle de confiance pour une
estimation ponctuelle de la précision (Niculescu et al, 2020). Un choix courant est un
intervalle de confiance de 95 % (Olofsson et al, 2014 ; Stehman & Foody, 2019). Dans
les évaluations de la précision de la couverture terrestre, la méthode standard pour
construire un intervalle de confiance consiste à supposer une distribution normale ou
une estimation ponctuelle avec un écart-type égal à l'estimation résultante de la marge
d'erreur, c'est-à-dire la racine carrée de la variance estimée divisée par la taille de
l'échantillon (Newcombe, 1998). Dans cette étude, les calculs des intervalles de
confiance ont été automatisés dans le script python inclus dans QGIS.
9.6. Résultats
Les valeurs de rétrodiffusion sont les plus faibles pour VH (b), légèrement plus
élevées pour VV (c) en 2020. VH a la plus grande amplitude et présente des pics de
rétrodiffusion importants, notamment pour la rizière inondée. La variation des valeurs
du coefficient de rétrodiffusion montre l'influence de l'humidité de surface sur le signal
radar, notamment pour la saison humide et pour les rizières inondées. Par exemple, à
partir de juin, les valeurs du coefficient de rétrodiffusion sont très faibles pour la rizière
inondée avec des valeurs comprises entre -25 et -30db en VH et entre -20 et -25db en
116
VV (Figure 43). Le coefficient a montré des valeurs faibles jusqu'à la fin de la saison
humide (jusqu'en novembre) pour la rizière inondée. La saison rizicole automne-hiver
est dominée dans certaines zones par les inondations et cette saison peut donc être
considérée comme moins productive. La période d'inondation qui coïncide avec la
saison du riz d'automne-hiver doit être considérée comme une période de "repos" dans
le cycle de croissance du riz, comme une période avec une production de riz réduite.
En saison sèche pour les deux rizières différentes et en saison humide pour la
rizière non inondée, la variation du coefficient en VH et VV des deux rizières différentes
s'explique par la phénologie différente du riz sur les deux rizières.
117
les rizières non inondées une augmentation des valeurs des coefficients de
rétrodiffusion pour la période végétative et reproductive des trois saisons rizicoles, et
deux saisons rizicoles pour la rizière inondée.
118
Selon les Figures 43 et 42, la saison des pluies en 2019 commence plus tard
(fin juin) et se termine plus tôt (en novembre). Par ailleurs, les précipitations maximales
ne dépassent pas 45 mm et les quantités les plus importantes ont été enregistrées à
partir du mois d'août. Le niveau maximum d'environ 3,5 m a été enregistré à la mi-
septembre. Les valeurs de rétrodiffusion sont les plus faibles pour VH (b), légèrement
plus élevées pour VV (c) en 2019, tout comme en 2020. Les valeurs élevées du
coefficient de rétrodiffusion correspondent à une rizière au développement végétatif
119
pendant la phase de croissance du riz au cours des trois saisons pour les rizières non
inondées, et des deux saisons pour les rizières inondées. Les deux polarisations
illustrent cette sensibilité du signal RADAR à la croissance du riz. Par ailleurs, nous
constatons également que le signal RADAR est sensible au moment de la récolte. En
effet, une diminution du signal a été observée après la récolte (Figure 44).
Précision globale
Mois
2019 2020
120
Tableau 23 : Précisions de l'utilisateur par le modèle CNN chaque mois pendant la saison humide en en
2019 et 2020
Précision de l’utilisateur
Tableau 24 : Précisions du producteur par le modèle CNN chaque mois pendant la saison humide en 2019
et 2020.
Précision du producteur
En général, les PA par mois de la saison des pluies pour les zones inondées
sont inférieurs aux PA des zones non inondées. Cette tendance se vérifie pour les
deux années. En 2019, les valeurs de PA des zones inondées pour les 7 mois
observés varient de 79,8% (en juillet) à 97,73 (septembre et octobre). En effet, durant
cette année 2019, le mois de juillet correspond à l'extension minimale des zones
121
inondées tandis que les mois de septembre et octobre correspondent à l'extension
maximale. Pour les zones non inondées, les erreurs d'Omission en 2019 vont de
99,09% pour l'un des mois les plus inondés (novembre) à 100% pour le mois le moins
inondé (juillet).
Ainsi, on remarque qu’en 2020, les valeurs de l'AP pour les zones inondées par
mois suivent la même logique, avec les valeurs les plus faibles pour les mois
d'extension moindre des zones inondées (juillet 88,34%) et les valeurs les plus élevées
pour les mois d'extension maximale des zones inondées (octobre avec 99,59% puis
novembre avec 98,92%). Durant la même année 2020, la tendance des valeurs de PA
est la même, en effet, des valeurs élevées, 99,99% pour le mois de juin et 99,48%
pour le mois d’octobre qui marque plus d’inondation.
Selon les résultats du tableau 22, les valeurs de l'UA sont plus cohérentes et
plus élevées pour les zones non inondées que pour les zones inondées pour les deux
saisons des deux années. En effet, on note qu’en 2019, des valeurs faibles (90,20 %)
ont été enregistrées pour la période de décembre à la fin de la saison humide, c’est-
à-dire les mois qui ont marqué moins d’inondation, tandis que ces valeurs de UA sont
plus élevées (99,89 %) pour le mois d’août, considéré comme un mois intermédiaire
pour le phénomène de l’inondation. Des UA supérieures à 99% ont été constatés dans
le cas des zones non inondées en 2019, et ce, pour tous les mois sauf ceux avec une
extension maximale des zones inondées (octobre et novembre avec des valeurs
supérieures à 98%).
En 2020, dans le cas des zones inondées, les valeurs les plus faibles (88,66%)
ont été enregistrées durant le mois d’extension maximale des inondations (en
novembre). La valeur la plus élevée (99,70%) de l'UA est atteinte au mois de juin avec
très peu d'inondations. En ce qui concerne les zones non inondées, les valeurs de l'UA
sont sans exception supérieures à 99% avec un maximum en septembre et un
minimum en novembre.
122
Les résultats suivant les calculs de l'AO et leurs intervalles de confiance
appliquant les classifications du modèle RF pour la série temporelle S1 montrent la
différence de la précision du classificateur RF pour les mois de la saison humide
(Tableau 25).
Tableau 25 : Comparaison des précisions globales par le modèle RF chaque mois en de la saison
humide en 2019 et 2020.
Précision globale
Mois
2019 2020
Tableau 26 : Précisions de l'utilisateur par le modèle RF pour les mois de la saison humide en 2019 et
2020.
Précision de l’utilisateur
2019 2020
Mois
Non-
inondée Non-inondée inondée
inondée
123
Tableau 27 : Précisions du producteur par le modèle RF pour les mois de la saison humide en 2019 et
2020.
Précision du producteur
En 2019, une précision globale supérieure à 99% a été obtenue pour le mois de
juin tandis qu’en 2020 cette précision a été notée pour les mois de juin et juillet. Par
ailleurs, on note qu’avec le classificateur CNN, la précision la plus faible était de 96%
(novembre 2020), alors qu’avec le classificateur RF, les précisions les plus faibles sont
entre 92% (décembre 2019) et 93% (novembre 2020). En revanche, les deux
algorithmes de classification ont obtenu une précision globale supérieure à 92%, sauf
que la précision globale obtenue par le modèle CNN est supérieure à celle obtenue
par les classificateurs RF. Ainsi, la performance de classification globale montre que
le modèle CNN obtient une meilleure classification que le RF, pour tous les mois de la
période d'inondation.
Les PA (Tableau 27) pour les deux années des zones inondées, sont inférieurs
aux PA des zones non inondées. On constate la même chose pour l'application du
CNN. Il y a deux exceptions, les PA de novembre et décembre pour les zones
inondées en 2019 sont plus élevés par rapport aux PA des zones non inondées :
96,88% et 94,29% en novembre et 96,72% et 91,96% en décembre. Cette valeur de
96,72% au mois de décembre 2019 obtenue suite à l'application de la RF est même
supérieure au PA du CNN pour le même mois et la même année (89,99%). Une
situation similaire se présente pour le mois de juillet 2019, durant lequel la valeur du
PA est plus faible avec le CNN qu'avec l'application de la RF (79,80% et 91,18%).
124
Avec ces exceptions, en général, les PA obtenus suite à l'application du CNN sont plus
élevés pour les deux années et pour les deux classes majeures que les PA obtenus
suite à l'application du RF.
Les valeurs de l'UA (Tableau 26) sont significativement plus faibles pour les
zones inondées par rapport aux zones non inondées pour les deux années. Au mois
de décembre 2019, nous avons obtenu la valeur minimale de 62,05% pour les zones
inondées alors qu'avec le CNN nous avons obtenu 90,20%. Des valeurs inférieures
ont également été enregistrées en 2020 pour les mois d'octobre et de novembre
(84,15% et 81,45%). L'UA maximale était obtenu au d’août 2019, qui atteignait 99,75%
pour les zones inondées et 99,72% en octobre 2020 pour les zones non inondées,
alors qu’avec le CNN la valeur maximale de l'UA pour les zones inondées est toujours
notée pour le mois d'août (99,89%).
Dans cette étude, les CNN sont directement utilisés comme classificateurs pour
la cartographie des inondations régionales dans le delta du Mékong. L'application de
l'algorithme CNN sur les cinq images par mois pendant la saison des pluies permet de
délimiter correctement la classe d'inondation sur l'ensemble de la zone inondée.
125
En comparant l'étendue des inondations, nous constatons que les zones
inondées en 2019 sont considérablement plus grandes que les zones inondées en
2020. Ainsi, en 2019, les mois marqués par les périodes de grandes inondations
étaient septembre, octobre et novembre, tandis qu'en 2020, il y a un décalage
temporaire d'un mois et les mois les plus inondés sont donc octobre, novembre et
décembre. Le mois de juillet, qui est un mois de la saison humide, présente la plus
faible extension par les zones d'inondation pour les deux années. Les mois de
novembre et octobre de la saison humide présentent une extension maximale des
zones inondées en 2019 et en 2020 respectivement.
(A (B
) )
Figure 45 : Cartographie des inondations par classification CNN par mois dans la saison
humide de juin à décembre : (A) Zones inondées en 2019 et (B) Zones inondées en 2020
Par ailleurs, nous avons effectué une analyse spatiale basée sur les pixels pour
déterminer la surface en Km2 qui sont correctement délimités avec les cartes
résultantes du modèle CNN (Figure 46).
126
implique une organisation particulière de la culture du riz pour cette troisième saison
(Figure 46).
(A) (B)
Figure 46 : schéma synthétique a) Précipitations mensuelles moyennes ; niveaux d'eau mensuels
maximums en 2019 (A) et 2020 (B) ; b) Zone inondée en km2 pour 2019 (A) et 2020 (B).
Il faut noter que le riz cultivé dans ces régions est un riz profond ou flottant, dont
l'eau est fournie par les crues des rivières. Ainsi, les rendements restent faibles car ils
sont largement dépendants des conditions climatiques. Avec 8% des surfaces
cultivées, ce type de riz ne fournit que 3% de la production mondiale, par conséquent
il est progressivement remplacé par le riz irrigué au fur et à mesure des aménagements
hydrauliques.
Ainsi, le tableau 28 présente les résultats de cette analyse spatiale sous forme
de pourcentages relatifs au nombre total de pixels.
Tableau 28 : Analyse spatiales des zones impactées par les inondations.
2019 2020
classe Superficie(km2) Pourcentage Superficie(km2) Pourcentage
Eau 25.4227 1% 30.0524 1%
Rizière 2615.265 58% 2075.6704 53%
Village 409.666 9% 556.608 14%
Jardin 1208.079 27% 1133.9639 29%
Forêt 220.5616 5% 127.0772 3%
Inondable 4478.9943 100% 3923.3719 100%
A travers l’analyse spatiale (Tableau 28), nous constatons que les zones
inondées en 2019 (4478,9 km2) sont plus importantes que les zones inondées en 2020
127
(3923,3 km2). Ces différences dans l'étendue des inondations pourraient indiquer que
la crue de 2019 doit être analysée en corrélation avec le niveau maximal du Mékong
à la station Tan Chau en 2019 de près de 4 m par rapport au niveau maximal enregistré
en 2020 de moins de 3 m à la même station.
Dans cette section l’analyse se concentre, dans un premier temps, sur les
classes d’occupation/couverture du sol (LULC) avant la saison humide (Figure 47b)
pour déterminer quelles classes d'utilisation des terres ont été inondées et, dans un
second temps, sur les LULC non inondées (Figure 47c). Cette analyse s'est effectuée
à partir des informations sur l'extension maximale des zones inondées (partie des
zones inondables) pour deux mois, à savoir d’octobre 2019 et de novembre 2020
(Figure 47a).
2019 2020
classe Superficie(km2) Pourcentage Superficie(km2) Pourcentage
Eau 260.3604 2% 318.9245 3%
Rizière 5941.3851 53% 6060.6319 55%
Village 690.8062 6% 426.1447 4%
Jardin 1615.0679 15% 1272.2426 11%
Forêt 2605.3455 24% 3035.0214 27%
L'analyse spatiale de la carte d'occupation du sol montre que la zone bâtie, les
rivières / zones humides et les rizières sont plus importantes et plus claires en sorties
des classificateurs CNN que RF. On note, ainsi, que les forêts ont tendance à se
concentrer dans les zones humides de la zone d'étude.
Durant les deux périodes d’octobre 2019 et de novembre 2020, les rizières
occupent plus de la moitié de la zone d'étude (53% en 2019 et 55% en 2020), suivies
128
par les forêts (24% en 2019 et 27% en 2020) et les jardins avec des superficies, 15%
et 11% du paysage, respectivement pour 2019 et 2020. Le pourcentage du couvert
urbain était de 6% en 2019 et 4% en 2020, tandis que l'eau représente 2 à 3 % du
paysage.
Plus de la moitié des rizières ont été inondées durant les deux années (58% en
2019 et 53% en 2020). En 2020, la superficie de rizières est plus importante qu'en
2019, cependant pour cette dernière la superficie des rizières impactée par les
inondations était plus grande. Ces inondations ont eu un impact significatif sur une
autre classe d'utilisation des terres ayant une valeur économique locale importante. A
l'exception des forêts qui ont été moins impactées en 2020, la situation des autres
types d’utilisation du sol est plus impactée en 2020 ; on note ainsi que l'urbain a été
beaucoup plus touché par l'inondation en 2020 (14%) par rapport à la situation de 2019
(9%), et de même pour les jardins où 29% et 27% de leur superficie est impactée en
2020 et 2019 respectivement.
Les superficies des LULC non-inondées en 2019 et 2020 sont présentées dans
le tableau 30.
129
Tableau 30 : Superficie des LULC non inondés en 2019 et 2020
2019 2020
classe Superficie(km2) Pourcentage Superficie(km2) Pourcentage
Eau 228.4982 2% 204.5155 2%
Rizière 2554.6812 23% 2866.4995 26%
Village 409.666 4% 556.608 5%
Jardin 1208.079 11% 1522.9639 14%
Forêt 2220.561 20% 2027.0772 18%
Inondable 4478.9943 40% 3923.3719 35%
Zones
inondées (a)
Occupation
du sol avant
les
inondations
(b)
Zones
inondées et
non
inondées ©
2019 2020
130
9.8.3 Discussion
131
Conclusion générale
Le présent travail a comme objectif principal d'étudier en détail les méthodes
d'apprentissage automatique supervisé et non supervisé appliquées dans le domaine
de l'environnement. C’est une contribution dans l’étude de l’occupation du sol des
zones fragiles telles que les systèmes deltaïques. Il s’agit d’un développement
méthodologique basé sur l’apport des nouvelles méthodes de télédétection pour le
suivi de l’occupation du sol au niveau de deux grands Delta au Vietnam, le delta du
fleuve Rouge et le delta du Mékong. Ainsi, la thèse a permis de mettre en évidence
les progrès qui ont été réalisés dans l'utilisation de données satellitaires de haute
résolution spatiale, ainsi que des méthodes de l’Intelligence Artificielle (Machine et
Deep Learning).
Pour aboutir à nos objectifs nous avons utilisé des données multi-capteurs et
des séries temporelles. Une modélisation de l'organisation des données satellites
optiques et radar a été effectuée afin d’améliorer les résultats d'analyse d'images. Ces
données ont été utilisées pour deux sortes d’application, la cartographie du type de
végétation et le suivi de l'extension de l'inondation et des zones inondables. En effet,
par rapport à des études d'analyse d'images uniques, les données optiques de séries
temporelles et les méthodes de traitement ici développées ont démontré leur efficacité
et faisabilité dans le processus de cartographie de la végétation du delta du fleuve
Rouge. Par ailleurs, les séries temporelles des données SAR se sont avérées
efficaces pour surveiller et cartographier les zones inondables au niveau du delta du
Mékong. Pour cette dernière application une approche systémique regroupant les
données environnementales et météorologiques ainsi que des données de
télédétection a permis d'étudier les zones inondables et de cartographier l'extension
de l’inondation. Cette approche met en avant la contribution des données satellites et
les technologies de télédétection; et aide à tracer les plans de prévention du risque et
orienter les politiques de gestion spatiales des zones deltaïques.
Par ailleurs, les méthodes d'apprentissage non supervisées ont été appliquées
pour effectuer des segmentations et extraire des objets géographiques d’images
satellitaires. Des algorithmes de segmentation ont été étudiés et appliqués lors du
développement méthodique afin d’identifier le meilleur choix pour notre domaine
d’application. Ces algorithmes ont été vérifiés et validés via une approche comparative
et à travers plusieurs essais.
132
Ainsi, nous avons constaté que le choix de la méthode de segmentation affecte
grandement les résultats de l'analyse d'images. Chaque méthode de segmentation est
optimisée et adaptée aux données en ajustant ces paramètres. Grâce à l'évaluation et
à la comparaison des résultats, la méthode de segmentation multi-résolution
(Mutiresolution Segmentation) a été choisie pour être utilisée pour le processus de
classification. Ce choix repose sur des mesures de similarité qui comprennent la
mesure de similarité spectrales et des couleurs ainsi que sur des coefficients de
rétrodiffusion et la mesure de similarité de la forme et la taille de chaque objet
géographique.
133
d'images. A partir de là, il est possible de construire un modèle GEOBIA avec le
changement de l'étape de prétraitement des données pour utiliser plus efficacement
les applications de la méthode d'apprentissage profond de CNN.
134
Bibliographie
Abdikan, S., Balik Sanli, F., Üstüner, M., Calò, F., 2016. LAND COVER MAPPING
USING SENTINEL-1 SAR DATA, ISPRS - International Archives of the
Photogrammetry, Remote Sensing et Spatial Information Sciences.
[Link]
Achanta, R., Shaji, A., Smith, K., Lucchi, A., Fua, P., Süsstrunk, S., 2010. SLIC
superpixels. Technical report, EPFL.
[Link], M., Shafri, H., Hamedianfar, A., 2017. New semi-automated mapping of
asbestos cement roofs using rule-based object-based image analysis et
Taguchi optimization technique from WorldView-2 images. International
Journal of Remote Sensing 38, 467–491.
[Link]
Awad, M., Khanna, R., 2015. Machine Learning, in: Awad, M., Khanna, R. (Eds.),
Efficient Learning Machines: Theories, Concepts, et Applications for Engineers
et System Designers. Apress, Berkeley, CA, pp. 1–18.
[Link]
Baatz, M., Hoffmann, C., Willhauck, G., 2008. Progressing from object-based to object-
oriented image analysis, in: Blaschke, T., Lang, S., Hay, G.J. (Eds.), Object-
Based Image Analysis: Spatial Concepts for Knowledge-Driven Remote
Sensing Applications, Lecture Notes in Geoinformation et Cartography.
Springer, Berlin, Heidelberg, pp. 29–42. [Link]
77058-9_2
Balzter, H., Cole, B., Thiel, C., Schmullius, C., 2015. Mapping CORINE Land Cover
from Sentinel-1A SAR et SRTM Digital Elevation Model Data using Random
Forests. Remote Sensing 7, 14876–14898.
[Link]
Bé, N., Van, T., Hằng, T., Ân, N., 2017. Ảnh hưởng của xâm nhập mặn đến công tác
quản lý nguồn tài nguyên nước trong sản xuất nông nghiệp tại huyện Long
Phú, tỉnh Sóc Trăng. Can Tho University, Journal of Science 52, 104.
[Link]
135
Bengio, Y., Courville, A., Vincent, P., 2014. Representation Learning: A Review et New
Perspectives. arXiv:1206.5538 [cs].
Benz, U.C., Hofmann, P., Willhauck, G., Lingenfelder, I., Heynen, M., 2004. Multi-
resolution, object-oriented fuzzy analysis of remote sensing data for GIS-ready
information. ISPRS Journal of Photogrammetry et Remote Sensing,
Integration of Geodata et Imagery for Automated Refinement et Update of
Spatial Databases 58, 239–258.
[Link]
Bianchi, T., Allison, M., 2009. Large-river delta-front estuaries as natural “recorders” of
global environmental change. Proceedings of the National Academy of
Sciences. [Link]
Blaschke, T., 2010. Object based image analysis for remote sensing. ISPRS Journal
of Photogrammetry et Remote Sensing 65, 2–16.
[Link]
Blaschke, T., Hay, G.J., Kelly, M., Lang, S., Hofmann, P., Addink, E., Queiroz Feitosa,
R., van der Meer, F., van der Werff, H., van Coillie, F., Tiede, D., 2014.
Geographic Object-Based Image Analysis - Towards a new paradigm. ISPRS
Journal of Photogrammetry et Remote Sensing 87, 180–191.
[Link]
Boualleg, Y., Farah, M., Farah, I.R., 2019. Remote Sensing Scene Classification Using
Convolutional Features et Deep Forest Classifier. IEEE Geoscience et Remote
Sensing Letters 16, 1944–1948. [Link]
Bovenga, F., 2020. Special Issue “Synthetic Aperture Radar (SAR) Techniques et
Applications.” Sensors (Basel) 20, 1851. [Link]
Breiman, L. (Ed.), 1998. Classification et regression trees, 1. CRC Press repr. ed.
Chapman & Hall/CRC, Boca Raton, Fla.
Brown, P.R., Tuan, V.V., Nhan, D.K., Dung, L.C., Ward, J., 2018. Influence of
livelihoods on climate change adaptation for smallholder farmers in the
136
Mekong Delta Vietnam. null 16, 255–271.
[Link]
Ca, V.T., Vongvisessomjai, S., Asaeda, T., 1994. Study on Salinity Intrusion in the Red
River Delta. Environmental Systems Research 22, 213–218.
[Link]
Castilla, G., Hay, G. J., 2008. Image objects et geographic objects, in: Blaschke, T.,
Lang, S., Hay, Geoffrey J. (Eds.), Object-Based Image Analysis: Spatial
Concepts for Knowledge-Driven Remote Sensing Applications, Lecture Notes
in Geoinformation et Cartography. Springer, Berlin, Heidelberg, pp. 91–110.
[Link]
Cian, F., Marconcini, M., Ceccato, P., 2018. Normalized Difference Flood Index for
rapid flood mapping: Taking advantage of EO big data. Remote Sensing of
Environment 209, 712–730. [Link]
Cloude, S., Zebker, H., 2010. Polarisation: Applications in Remote Sensing. Physics
Today - PHYS TODAY 63. [Link]
Comaniciu, D., Ramesh, V., Meer, P., 2001. The variable bandwidth mean shift et data-
driven scale selection, in: Proceedings Eighth IEEE International Conference
on Computer Vision. ICCV 2001. Presented at the Proceedings Eighth IEEE
International Conference on Computer Vision. ICCV 2001, pp. 438–445 vol.1.
[Link]
Congalton, R., Green, K., 2019. Assessing the Accuracy of Remotely Sensed Data:
Principles et Practices, Third Edition. [Link]
Cortes, C., Vapnik, V., 1995. Support-vector networks. Mach Learn 20, 273–297.
[Link]
Cosslett, T.L., Cosslett, P.D., 2014. Major Threats to Mekong Delta: Climate Change
et Mainstream Dams, in: Cosslett, T.L., Cosslett, P.D. (Eds.), Water Resources
et Food Security in the Vietnam Mekong Delta, Natural Resource Management
et Policy. Springer International Publishing, Cham, pp. 75–96.
[Link]
Dai, Q., Wang, L., Xu, Q., Zhang, Y., 2014. An mean shift algorithm with adaptive
bandwidth et weight selection for high spatial remotely sensed imagery
segmentation, in: 2014 IEEE Geoscience et Remote Sensing Symposium.
137
Presented at the 2014 IEEE Geoscience et Remote Sensing Symposium, pp.
1592–1595. [Link]
Dang, T.H., Coynel, A., Orange, D., Blanc, G., Etcheber, H., Le, L.A., 2010. Long-term
monitoring (1960–2008) of the river-sediment transport in the Red River
Watershed (Vietnam): Temporal variability et dam-reservoir impact. Science
of The Total Environment 408, 4654–4664.
[Link]
Đạt T.Q., Trung N.H., Likitdecharote K., 2012. MÔ PHỎNG XÂM NHẬP MẶN ĐỒNG
BẰNG SÔNG CỬU LONG DƯỚI TÁC ĐỘNG MỰC NƯỚC BIỂN DÂNG VÀ
SỰ SUY GIẢM LƯU LƯỢNG TỪ THƯỢNG NGUỒN. Tạp chí Khoa học
Trường Đại học Cần Thơ Số 21b, 141–150.
Delgado, J.M., Apel, H., Merz, B., 2010. Flood trends et variability in the Mekong river.
Hydrology et Earth System Sciences 14, 407–418.
[Link]
Dhara, S., Dang, T., Parial, K., Lu, X.X., 2020. Accounting for Uncertainty et
Reconstruction of Flooding Patterns Based on Multi-Satellite Imagery et
Support Vector Machine Technique: A Case Study of Can Tho City, Vietnam.
Water 12, 1543. [Link]
Do, N.H., Egashira, K., 2002. Irrigation et drainage systems et their impact on land use
et rice production after renovation (Doi Moi) in the Red River delta of Vietnam
- A case study in Tien Son District, Bac Ninh Province. Journal- Faculty of
Agriculture Kyushu University 46, 411–422.
Dobrinić, D., Gašparović, M., Medak, D., 2021. Sentinel-1 et 2 Time-Series for
Vegetation Mapping Using Random Forest Classification: A Case Study of
Northern Croatia. Remote Sensing 13, 2321.
[Link]
Dorian, B., 2015. L’utilisation de l’agent orange durant la guerre du Vietnam et ses
conséquences sur la santé.
Drǎguţ, L., Tiede, D., Levick, S.R., 2010. ESP: a tool to estimate scale parameter for
multiresolution image segmentation of remotely sensed data. International
138
Journal of Geographical Information Science 24, 859–871.
[Link]
Dumanski, J., Pieri, C., 2000. Land quality indicators: research plan. Agriculture,
Ecosystems & Environment 81, 93–102. [Link]
8809(00)00183-3
Fang, G., Kwok, Y.-K., Yu, K., Zhu, Y., 1999. Numerical simulation of principal tidal
constituents in the South China Sea, Gulf of Tonkin et Gulf of Thailand.
Continental Shelf Research 19, 845–869. [Link]
4343(99)00002-3
FAO, 2016. FAO Yearbook. Fishery et Aquaculture Statistics 2014. FAO annuaire.
Statistiques des pêches et de l’aquaculture 2014. FAO anuario. Estadísticas
de pesca y acuicultura 2014, FAO Yearbook of Fishery et Aquaculture
Statistics. FAO, Rome, Italy.
Filipponi, F., 2019. Sentinel-1 GRD Preprocessing Workflow. Proceedings 18, 11.
[Link]
Foley, J., Defries, R., Asner, G., Barford, C., Bonan, G., Carpenter, S., Chapin III, F.S.,
Coe, M., Daily, G., Gibbs, H., Helkowski, J., Holloway, T., Howard, E.,
Kucharik, C., Monfreda, C., Patz, J., Prentice, I., Ramankutty, N., Snyder, P.,
2005. Global Consequences of Land Use. Science (New York, N.Y.) 309, 570–
4. [Link]
Foley, J.A., DeFries, R., Asner, G.P., Barford, C., Bonan, G., Carpenter, S.R., Chapin,
F.S., Coe, M.T., Daily, G.C., Gibbs, H.K., Helkowski, J.H., Holloway, T.,
Howard, E.A., Kucharik, C.J., Monfreda, C., Patz, J.A., Prentice, I.C.,
Ramankutty, N., Snyder, P.K., 2005. Global Consequences of Land Use.
Science 309, 570–574. [Link]
Foody, G.M., Mathur, A., 2004. Toward intelligent training of supervised image
classifications: directing training data acquisition for SVM classification.
Remote Sensing of Environment 93, 107–117.
[Link]
Fourie, C., Schoepfer, E., 2014. Data Transformation Functions for Expanded Search
Spaces in Geographic Sample Supervised Segment Generation. Remote
Sensing 6, 3791–3821. [Link]
Frenken, K., 2012. Mekong River Basin.” In Irrigation in Southern et Eastern Asia in
Figures: Aquastat Survey. FAO Water Reports.
Fu, T., Ma, L., Li, M., Johnson, B.A., 2018. Using convolutional neural network to
identify irregular segmentation objects from very high-resolution remote
sensing imagery. JARS 12, 025010. [Link]
139
Fukunaga, K., Hostetler, L., 1975. The estimation of the gradient of a density function,
with applications in pattern recognition. IEEE Transactions on Information
Theory 21, 32–40. [Link]
Gao, Y., Mas, J., Kerle, N., Pacheco, J., 2011. Optimal region growing segmentation
et its effect on classification accuracy. International Journal of Remote Sensing
32, 3747–3763. [Link]
Gilard, O., 2006. Risques d’inondation dans le delta du fleuve Rouge. Herodote no
121, 42–54.
Gislason, P.O., Benediktsson, J.A., Sveinsson, J.R., 2006. Random Forests for land
cover classification. Pattern Recognition Letters, Pattern Recognition in
Remote Sensing (PRRS 2004) 27, 294–300.
[Link]
Gong, P., Wang, J., Yu, Le, Zhao, Yongchao, Zhao, Yuanyuan, Liang, L., Niu, Z.,
Huang, X., Fu, H., Liu, S., Li, C., Li, X., Fu, W., Liu, C., Xu, Y., Wang, X.,
Cheng, Q., Hu, L., Yao, W., Zhang, Han, Zhu, P., Zhao, Z., Zhang, Haiying,
Zheng, Y., Ji, L., Zhang, Y., Chen, H., Yan, A., Guo, J., Yu, Liang, Wang, L.,
Liu, X., Shi, T., Zhu, M., Chen, Y., Yang, G., Tang, P., Xu, B., Giri, C., Clinton,
N., Zhu, Z., Chen, Jin, Chen, Jun, 2013. Finer resolution observation et
monitoring of global land cover: first mapping results with Landsat TM et ETM+
data. International Journal of Remote Sensing 34, 2607–2654.
[Link]
Gu, H., Li, H., Yan, L., Lu, X., 2015. A Framework for Geographic Object-Based Image
Analysis (GEOBIA) based on geographic ontology. ISPRS - International
Archives of the Photogrammetry, Remote Sensing et Spatial Information
Sciences XL-7/W4, 27–33. [Link]
2015
Gupta, A., 2009. Chapter 3 - Geology et Landforms of the Mekong Basin, in: Campbell,
I.C. (Ed.), The Mekong, Aquatic Ecology. Academic Press, San Diego, pp. 29–
51. [Link]
Ha, T.T.P., van Dijk, H., Visser, L., 2014. Impacts of changes in mangrove forest
management practices on forest accessibility et livelihood: A case study in
mangrove-shrimp farming system in Ca Mau Province, Mekong Delta,
Vietnam. Land Use Policy 36, 89–101.
[Link]
Hall, J., Arheimer, B., Borga, M., Brázdil, R., Claps, P., Kiss, A., Kjeldsen, T.R.,
Kriaučiūnienė, J., Kundzewicz, Z.W., Lang, M., Llasat, M.C., Macdonald, N.,
McIntyre, N., Mediero, L., Merz, B., Merz, R., Molnar, P., Montanari, A.,
Neuhold, C., Parajka, J., Perdigão, R. a. P., Plavcová, L., Rogger, M., Salinas,
J.L., Sauquet, E., Schär, C., Szolgay, J., Viglione, A., Blöschl, G., 2014.
Understanding flood regime changes in Europe: a state-of-the-art assessment.
140
Hydrology et Earth System Sciences 18, 2735–2772.
[Link]
Hamada, M.A., Kanat, Y., Abiche, A.E., 2019. Multi-Spectral Image Segmentation
Based on the K-means Clustering 9, 4.
Haruyama, S., 1995. Geomorphologic features of the Song Hong delta. J Int Relat Stud
21, 1–13.
Hay, G., Castilla, G., 2008. Geographic Object-Based Image Analysis (GEOBIA): A
new name for a new discipline, in: Object-Based Image Analysis - Spatial
Concepts for Knowledge-Driven Remote Sensing Applications. pp. 75–89.
[Link]
Hay, G.J., 1994. Visualizing 3-D texture: a three dimensional structural approach to
model forest texture. Can. J. Remote Sensing 90–101.
Hay, G.J., Niernann, K.O., Goodenough, D.G., 1997. Spatial thresholds, image-
objects, et upscaling: A multiscale evaluation. Remote Sensing of Environment
62, 1–19. [Link]
Hinton, G.E., Osindero, S., Teh, Y.-W., 2006. A Fast Learning Algorithm for Deep Belief
Nets. Neural Computation 18, 1527–1554.
[Link]
Hirabayashi, Y., Mahendran, R., Koirala, S., Konoshima, L., Yamazaki, D., Watanabe,
S., Kim, H., Kanae, S., 2013. Global flood risk under climate change. Nature
Climate Change 3, 816–821. [Link]
Hoanh, C.T., , N.D. Phong, T.P. Tuong, J.W. Gowing, N.V. Ngoc, L.T. An, et N.V.
MinhPhong, T.P. Tuong, J.W. Gowing, N.V. Ngoc, L.T. An,, Minh, et N.V.,
2002. Water Management For Sustainable Agricultural And Aquacultural
Development In Bac Lieu Province. Paper presented at the Mid-term Review
IRRI, University of Newcastle, ICLARM, University of Can Tho, Sub-Institute
of Water Resources Planning, Department of Agriculture & Rural
Development, Bac Lieu Province, Vietnam. Presented at the Workshop in Bac
Lieu June 2002.
141
Hong Quang, N., Tuan, V.A., Thi Thu Hang, L., Manh Hung, N., Thi The, D., Thi Dieu,
D., Duc Anh, N., Hackney, C.R., 2020. Hydrological/Hydraulic Modeling-Based
Thresholding of Multi SAR Remote Sensing Data for Flood Monitoring in
Regions of the Vietnamese Lower Mekong River Basin. Water 12, 71.
[Link]
Huang, C., Davis, L., Townshend, J., 2002. An assessment of support vector machines
for land cover classification. International Journal of Remote Sensing 23.
[Link]
Hung, B.V., 2017. Identify the Major Reasons to Cause Vulnerability to Mekong Delta
Under the Impacts of Drought et Climate Change, in: Kurisu, F., Ramanathan,
AL., Kazmi, A.A., Kumar, M. (Eds.), Trends in Asian Water Environmental
Science et Technology. Springer International Publishing, Cham, pp. 211–222.
[Link]
IMHEN, 2011. Scenarios of climate change et sea level rises for Viet Nam. Ministry of
Natural Resources et Environments, Hanoi.
Ienco, D., Interdonato, R., Gaetano, R., Ho Tong Minh, D., 2019. Combining Sentinel-
1 et Sentinel-2 Satellite Image Time Series for land cover mapping via a multi-
source deep learning architecture. ISPRS Journal of Photogrammetry et
Remote Sensing 158, 11–22. [Link]
Jozdani, S., Johnson, B., Chen, D., 2019. Comparing Deep Neural Networks,
Ensemble Classifiers, et Support Vector Machine Algorithms for Object-Based
Urban Land Use/Land Cover Classification. Remote Sensing 11.
[Link]
Kansakar, P., Hossain, F., 2016. A review of applications of satellite earth observation
data for global societal benefit et stewardship of planet earth. Space Policy 36,
46–54. [Link]
Kavitha, A.V., Srikrishna, A., Satyanarayana, C., 2021. A Review on Detection of Land
Use et Land Cover from an Optical Remote Sensing Image. IOP Conf. Ser.:
Mater. Sci. Eng. 1074, 012002. [Link]
899X/1074/1/012002
Khan, S.S., Ahmad, A., 2004. Cluster center initialization algorithm for K-means
clustering. Pattern Recognition Letters 25, 1293–1302.
[Link]
Kondolf, G.M., Rubin, Z.K., Minear, J.T., 2014. Dams on the Mekong: Cumulative
sediment starvation. Water Resources Research 50, 5158–5169.
[Link]
142
Kontgis, C., Schneider, A., Ozdogan, M., 2015. Mapping rice paddy extent et
intensification in the Viêt-Namese Mekong River Delta with dense time stacks
of Landsat data. Remote Sensing of Environment 169, 255–269.
[Link]
Kuemmerle, T., Chaskovskyy, O., Knorn, J., Radeloff, V.C., Kruhlov, I., Keeton, W.S.,
Hostert, P., 2009. Forest cover change et illegal logging in the Ukrainian
Carpathians in the transition period from 1988 to 2007. Remote Sensing of
Environment 113, 1194–1207. [Link]
Lambin, E.F., Geist, H., n.d. Land-use et Land-cover Change: Local Processes et
Global Impacts. Springer London NetLibrary, Inc. [distributor], Guildford
Boulder.
Lambin, E.F., Geist, H.J., 2006. Land-Use et Land-Cover Change: Local Processes et
Global Impacts. Springer Science & Business Media.
Lam-Dao, N., Pham-Bach, V., Nguyen-Thanh, M., Pham-Thi, M.-T., Hoang-Phi, P.,
2011. Change Detection of Land use et Riverbank in Mekong Delta, Vietnam
using Time Series Remotely Sensed Data. jore 2, 370–374.
[Link]
Lary, D.J., Alavi, A.H., Gandomi, A.H., Walker, A.L., 2016. Machine learning in
geosciences et remote sensing. Geoscience Frontiers, Special Issue:
Progress of Machine Learning in Geosciences 7, 3–10.
[Link]
Le, T.P.Q., Garnier, J., Gilles, B., Sylvain, T., Van Minh, C., 2007. The changing flow
regime et sediment load of the Red River, Viet Nam. Journal of Hydrology 334,
199–214. [Link]
LeCun, Y., Bengio, Y., Hinton, G., 2015. Deep Learning. Nature 521, 436–44.
[Link]
Lee, J., Jurkevich, L., Dewaele, P., Wambacq, P., Oosterlinck, A., 1994a. Speckle
filtering of synthetic aperture radar images: A Review. Remote Sensing
Reviews 8. [Link]
Lee, J., Jurkevich, L., Dewaele, P., Wambacq, P., Oosterlinck, A., 1994b. Speckle
filtering of synthetic aperture radar images: A Review. Remote Sensing
Reviews 8. [Link]
Lee, J.S., Pottier, E., 2009. Polarimetric Radar Imaging : From basics to applications.,
Optical Science et Engineering. CRC Press, Taylor & Francis Editor,.
Lewis, D.D., Catlett, J., 1994. Heterogeneous Uncertainty Sampling for Supervised
Learning, in: In Proceedings of the Eleventh International Conference on
Machine Learning. Morgan Kaufmann, pp. 148–156.
143
Li, Z., Saito, Y., Matsumoto, E., Wang, Y., Tanabe, S., Lan Vu, Q., 2006. Climate
change et human impact on the Song Hong (Red River) Delta, Vietnam, during
the Holocene. Quaternary International, Precipitation: it’s influence on
Quaternary events in eastern Asia 144, 4–28.
[Link]
Liu, D., Xia, F., 2010. Assessing object-based classification: Advantages et limitations.
Remote Sensing Letters 1, 187–194.
[Link]
Liu, S., Li, X., Chen, D., Duan, Y., Ji, H., Zhang, L., Chai, Q., Hu, X., 2020.
Understanding Land use/Land cover dynamics et impacts of human activities
in the Mekong Delta over the last 40 years. Global Ecology et Conservation
22, e00991. [Link]
Longepe, N., Rakwatin, P., Isoguchi, O., Shimada, M., Uryu, Y., Yulianto, K., 2011.
Assessment of ALOS PALSAR 50 m Orthorectified FBD Data for Regional
Land Cover Classification by Support Vector Machines. IEEE Transactions on
Geoscience et Remote Sensing 49, 2135–2150.
[Link]
Lowe, B., Kulkarni, A., 2015. Multispectral Image Analysis Using Random Forest.
International Journal on Soft Computing 6, 1–14.
[Link]
Luu, T.N.M., Garnier, J., Billen, G., Orange, D., Némery, J., Le, T.P.Q., Tran, H.T., Le,
L.A., 2010. Hydrological regime et water budget of the Red River Delta
(Northern Vietnam). Journal of Asian Earth Sciences 37, 219–228.
[Link]
Lv, X., Ming, D., Chen, Y., Wang, M., 2018. Very high resolution remote sensing image
classification with SEEDS-CNN et scale effect analysis for superpixel CNN
classification. International Journal of Remote Sensing 40, 1–26.
[Link]
Ma, L., Du, B., Chen, H., Soomro, N.Q., 2016. Region-of-Interest Detection via
Superpixel-to-Pixel Saliency Analysis for Remote Sensing Image. IEEE
Geoscience et Remote Sensing Letters 13, 1752–1756.
[Link]
Ma, L., Li, M., Ma, X., Cheng, L., Du, P., Liu, Y., 2017. A review of supervised object-
based land-cover image classification. ISPRS Journal of Photogrammetry et
Remote Sensing 130, 277–293.
[Link]
Ma, L., Liu, Y., Zhang, X., Ye, Y., Yin, G., Johnson, B.A., 2019. Deep learning in remote
sensing applications: A meta-analysis et review. ISPRS Journal of
Photogrammetry et Remote Sensing 152, 166–177.
[Link]
144
MacQueen, J., 1967. Some methods for classification et analysis of multivariate
observations. Presented at the Proceedings of the Fifth Berkeley Symposium
on Mathematical Statistics et Probability, Volume 1: Statistics, The Regents of
the University of California.
Manandhar, R., Odeh, I.O.A., Ancev, T., 2009. Improving the Accuracy of Land Use et
Land Cover Classification of Landsat Data Using Post-Classification
Enhancement. Remote Sensing 1, 330–344.
[Link]
Marceau, D.J., Howarth, P.J., Dubois, J.M., Gratton, D.J., 1990. Evaluation Of The
Grey-level Co-occurrence Matrix Method For Land-cover Classification Using
Spot Imagery. IEEE Transactions on Geoscience et Remote Sensing 28, 513–
519. [Link]
Mercier, A., Betbeder, J., Rumiano, F., Baudry, J., Gond, V., Blanc, L., Bourgoin, C.,
Cornu, G., Ciudad, C., Marchamalo, M., Poccard-Chapuis, R., Hubert-Moy, L.,
2019. Evaluation of Sentinel-1 et 2 Time Series for Land Cover Classification
of Forest–Agriculture Mosaics in Temperate et Tropical Landscapes. Remote
Sensing 11, 979. [Link]
Michel, J., Youssefi, D., Grizonnet, M., 2015. Stable Mean-Shift Algorithm et Its
Application to the Segmentation of Arbitrarily Large Remote Sensing Images.
IEEE Transactions on Geoscience et Remote Sensing 53, 952–964.
[Link]
Minderhoud, P.S.J., Coumou, L., Erban, L.E., Middelkoop, H., Stouthamer, E., Addink,
E.A., 2018. The relation between land use et subsidence in the Vietnamese
Mekong delta. Science of The Total Environment 634, 715–726.
[Link]
Michel, J., Youssefi, D., Grizonnet, M., 2015. Stable Mean-Shift Algorithm et Its
Application to the Segmentation of Arbitrarily Large Remote Sensing Images.
IEEE Transactions on Geoscience et Remote Sensing 53, 952–964.
[Link]
Moreira, A., Prats-Iraola, P., Younis, M., Krieger, G., Hajnsek, I., Papathanassiou, K.P.,
2013. A tutorial on synthetic aperture radar. IEEE Geoscience et Remote
Sensing Magazine 1, 6–43. [Link]
Morgan, J.N., Sonquist, J.A., 1963. Problems in the Analysis of Survey Data, et a
Proposal. Journal of the American Statistical Association 58, 415–434.
[Link]
Morinaka, Y., Sakamoto, T., Inukai, Y., Agetsuma, M., Kitano, H., Ashikari, M.,
Matsuoka, M., 2006. Morphological Alteration Caused by Brassinosteroid
145
Insensitivity Increases the Biomass et Grain Production of Rice1. Plant
Physiology. [Link]
Mujabar, S., Chandrasekar, N., 2011. A Shoreline Change analysis along the Coast
between Kanyakumari et Tuticorin, India using Remote sensing et GIS.
Arabian Journal of Geosciences 6. [Link]
4
NAS, 1974. NAS The Effects of Herbicides in South Vietnam National Academy of
Sciences, Washington (1974) Part a - Summary et Conclusions.
NEDECO, 1993. Mekong Delta Master Plan Study. A Perspective for Suitable
Development of Land et Water Resources, Volume I. World Bank.
Nguyen, T.T.H., Bie, C.A.J.M.D., Ali, A., Smaling, E.M.A., Chu, T.H., 2012. Mapping
the irrigated rice cropping patterns of the Mekong delta, Vietnam, through
hyper-temporal SPOT NDVI image analysis. International Journal of Remote
Sensing 33, 415–434. [Link]
NIAPP, 2001. Adjustment plan for agriculture, aquaculture, salt et forest production in
South of National Road 1A, Bac Lieu province. [in Vietnamese], Ho Chi Minh,
Vietnam.
Niculescu, S., Billey, A., Talab Ou Ali, H., 2019. Random Forest Classification using
Sentinel-1 et Sentinel-2 series for vegetation monitoring in the Pays de Brest
(France). SPIE j 10783.
Niculescu, S., Lam, C.N., 2019. Geographic Object-Based Image Analysis of Changes
in Land Cover in the Coastal Zones of the Red River Delta (Vietnam). Journal
of Environmental Protection 10, 413–430.
[Link]
Olofsson, P., Foody, G.M., Herold, M., Stehman, S.V., Woodcock, C.E., Wulder, M.A.,
2014. Good practices for estimating area et assessing accuracy of land
change. Remote Sensing of Environment 148, 42–57.
[Link]
Olofsson, P., Foody, G.M., Stehman, S.V., Woodcock, C.E., 2013. Making better use
of accuracy data in land change studies: Estimating accuracy et area et
quantifying uncertainty using stratified estimation. Remote Sensing of
Environment 129, 122–131. [Link]
Olofsson, P., Kuemmerle, T., Griffiths, P., Knorn, J., Baccini, A., Gancz, V., Viorel, B.,
Houghton, R., Abrudan, I., Woodcock, C., 2011. Carbon implications of forest
restitution in post-socialist Romania. Environmental Research Letters 6,
045202. [Link]
Pham, L.H., Pham, L.T.H., Dang, T.D., Tran, D.D., Dinh, T.Q., 2021. Application of
Sentinel-1 data in mapping land-use et land cover in a complex seasonal
146
landscape: a case study in the coastal area of Vietnamese Mekong Delta.
Geocarto International 0, 1–18.
[Link]
Pham, Q.S., 2004. Study of development of Red River – Thai Binh River estuaries on
the basis of Remote Sensing’s Information et GIS for rational exploitation of
use of territory. PhD thesis on Protection, rational use et generation of natural
resources. National University.
Phan, T.H., 2018. Rice monitoring using radar remote sensing (phdthesis). Université
Paul Sabatier - Toulouse III.
Phan, H., Le Toan, T., Bouvet, A., Lam, N., Duy, T., Zribi, M., 2018. Mapping of Rice
Varieties et Sowing Date Using X-Band SAR Data. Sensors 18, 316.
[Link]
Phung, D., Nguyen-Huy, T., Tran, N.N., Tran, D.N., Doan, V.Q., Nghiem, S., Nguyen,
N.H., Nguyen, T.H., Bennett, T., 2021. Hydropower dams, river drought et
health effects: A detection et attribution study in the lower Mekong Delta
Region. Climate Risk Management 32, 100280.
[Link]
Pokhrel, Y., Burbano, M., Roush, J., Kang, H., Sridhar, V., Hyndman, D.W., 2018. A
Review of the Integrated Effects of Changing Climate, Land Use, et Dams on
Mekong River Hydrology. Water 10, 266. [Link]
Potere, D., Schneider, A., 2007. A critical look at representations of urban areas in
global maps. GeoJournal 69, 55–80. [Link]
9102-z
Quang, N.N., 2002. Vietnam et the sustainable development of the Mekong river basis.
Water Science et Technology 45, 261–266.
[Link]
Rekik, A., Zribi, M., Benjelloun, M., Hamida, A. ben, 2006. A k-Means Clustering
Algorithm Initialization for Unsupervised Statistical Satellite Image
Segmentation, in: 2006 1ST IEEE International Conference on E-Learning in
Industrial Electronics. Presented at the 2006 1ST IEEE International
Conference on E-Learning in Industrial Electronics, pp. 11–16.
[Link]
Revilla-Romero, B., Hirpa, F.A., Pozo, J.T., Salamon, P., Brakenridge, R.,
Pappenberger, F., De Groeve, T., 2015. On the Use of Global Flood Forecasts
et Satellite-Derived Inundation Maps for Flood Monitoring in Data-Sparse
Regions. Remote Sensing 7, 15702–15728.
[Link]
147
Rodriguez-Galiano, V.F., Ghimire, B., Rogan, J., Chica-Olmo, M., Rigol-Sanchez, J.P.,
2012. An assessment of the effectiveness of a random forest classifier for land-
cover classification. ISPRS Journal of Photogrammetry et Remote Sensing 67,
93–104. [Link]
Rouse, J., Haas, R.H., Deering, D., Schell, J.A., Harlan, J., 1973. Monitoring the Vernal
Advancement et Retrogradation (Green Wave Effect) of Natural Vegetation.
[Great Plains Corridor]. undefined.
Roy, D.P., Wulder, M.A., Loveland, T.R., C.e., W., Allen, R.G., Anderson, M.C., Helder,
D., Irons, J.R., Johnson, D.M., Kennedy, R., Scambos, T.A., Schaaf, C.B.,
Schott, J.R., Sheng, Y., Vermote, E.F., Belward, A.S., Bindschadler, R.,
Cohen, W.B., Gao, F., Hipple, J.D., Hostert, P., Huntington, J., Justice, C.O.,
Kilic, A., Kovalskyy, V., Lee, Z.P., Lymburner, L., Masek, J.G., McCorkel, J.,
Shuai, Y., Trezza, R., Vogelmann, J., Wynne, R.H., Zhu, Z., 2014. Landsat-8:
Science et product vision for terrestrial global change research. Remote
Sensing of Environment 145, 154–172.
[Link]
Rumelhart, D., Hinton, G.E., Williams, R.J., 1986. Learning representations by back-
propagating errors. Nature. [Link]
Sakamoto, T., Van Phung, C., Kotera, A., Nguyen, K.D., Yokozawa, M., 2009. Analysis
of rapid expansion of inland aquaculture et triple rice-cropping areas in a
coastal area of the Vietnamese Mekong Delta using MODIS time-series
imagery. Landscape et Urban Planning 92, 34–46.
[Link]
Salzberg, S.L., 1994. C4.5: Programs for Machine Learning by J. Ross Quinlan.
Morgan Kaufmann Publishers, Inc., 1993. Mach Learn 16, 235–240.
[Link]
Schmullius, C.C., Evans, D.L., 1997. Review article Synthetic aperture radar (SAR)
frequency et polarization requirements for applications in ecology, geology,
hydrology, et oceanography: A tabular status quo after SIR-C/X-SAR.
International Journal of Remote Sensing 18, 2713–2722.
[Link]
Schubert, A., Small, D., Miranda, N., Geudtner, D., Meier, E., 2015. Sentinel-1A
Product Geolocation Accuracy: Commissioning Phase Results. Remote
Sensing 7, 9431–9449. [Link]
148
Singh, M., Ahuja, N., 2002. Mean-Shift Segmentation with Wavelet-based Bandwidth
Selection. [Link]
Smajgl, A., Toan, T.Q., Nhan, D.K., Ward, J., Trung, N.H., Tri, L.Q., Tri, V.P.D., Vu,
P.T., 2015. Responding to rising sea levels in the Mekong Delta. Nature
Climate Change 5, 167–174. [Link]
Spruce, J., Bolten, J., Mohammed, I.N., Srinivasan, R., Lakshmi, V., 2020. Mapping
Land Use Land Cover Change in the Lower Mekong Basin From 1997 to 2010.
Front. Environ. Sci. 8. [Link]
Stibig, H.-J., Achard, F., Fritz, S., 2004. A new forest cover map of continental
southeast Asia derived from SPOT-VEGETATION satellite imagery. Applied
Vegetation Science 7, 153–162. [Link]
109X.2004.tb00606.x
Strahler, A.H., Boschetti, L., Foody, G.M., Friedl, M.A., Hansen, M.C., Herold, M.,
Morisette, J.T., Stehman, S.V., 2006. Global Land Cover Validation:
Recommendations for Evaluation et Accuracy Assessment of Global et Cover
Maps.
Stutz, D., Hermans, A., Leibe, B., 2018. Superpixels: An Evaluation of the State-of-the-
Art. Computer Vision et Image Understanding 166, 1–27.
[Link]
Taylor, P., 2014. The khmer lands of Vietnam: Environment, cosmology et sovereignty.
NUS Press, Singapore.
Thanh Noi, P., Kappas, M., 2018. Comparison of Random Forest, k-Nearest Neighbor,
et Support Vector Machine Classifiers for Land Cover Classification Using
Sentinel-2 Imagery. Sensors 18, 18. [Link]
Thompson, J.R., Laizé, C.L.R., Green, A.J., Acreman, M.C., Kingston, D.G., 2014.
Climate change uncertainty in environmental flows for the Mekong River.
Hydrological Sciences Journal 59, 935–954.
[Link]
Thuan, N.T., 2006. The dry season flow in the Mekong Delta: The report collection of
the Science Technology conference. The Institute of Hydrology Technology et
Environments.
Ticehurst, C., Guerschman, J.P., Chen, Y., 2014. The Strengths et Limitations in Using
the Daily MODIS Open Water Likelihood Algorithm for Identifying Flood
Events. Remote Sensing 6, 11791–11809. [Link]
Tiede, D., Lang, S., Hoffmann, C., 2008. Domain-specific class modelling for one-level
representation of single trees, in: Blaschke, T., Lang, Stefan, Hay, G.J. (Eds.),
Object-Based Image Analysis: Spatial Concepts for Knowledge-Driven
Remote Sensing Applications, Lecture Notes in Geoinformation et
149
Cartography. Springer, Berlin, Heidelberg, pp. 133–151.
[Link]
Tilman, D., Fargione, J., Wolff, B., D’Antonio, C., Dobson, A., Howarth, R., Schindler,
D., Schlesinger, W.H., Simberloff, D., Swackhamer, D., 2001. Forecasting
Agriculturally Driven Global Environmental Change. Science 292, 281–284.
[Link]
Tomiyasu, K., 1978. Tutorial review of synthetic-aperture radar (SAR) with applications
to imaging of the ocean surface. Proceedings of the IEEE 66, 563–583.
[Link]
Tran, H., Tran, T., Kervyn, M., 2015. Dynamics of Land Cover/Land Use Changes in
the Mekong Delta, 1973–2011: A Remote Sensing Analysis of the Tran Van
Thoi District, Ca Mau Province, Vietnam. Remote Sensing 7, 2899–2925.
[Link]
Trung, N., 2006. Comparing land use planning approaches in the Mekong Delta,
Vietnam /.
Trung, N.H., Tri, L.Q., Mensvoort, M.E.F. van, Bregt, A.K., 2007. Application of GIS in
land-use planning: a case study in the Coastal Mekong Delta of Vietnam.
International Journal of Geoinformatics 3, 1–8.
Truong, T.D., Do, L.H., 2018. Mangrove forests et aquaculture in the Mekong river
delta. Land Use Policy 73, 20–28.
[Link]
van Maren, D.S., 2007. Water et sediment dynamics in the Red River mouth et
adjacent coastal zone. Journal of Asian Earth Sciences, Morphodynamics of
the Red River Delta, Vietnam 29, 508–522.
[Link]
Västilä, K., Kummu, M., Sangmanee, C., Chinvanno, S., 2010. Modelling climate
change impacts on the flood pulse in the Lower Mekong floodplains.
JOURNAL OF WATER AND CLIMATE CHANGE 1, 67–86.
[Link]
Vörösmarty, C.J., Green, P., Salisbury, J., Lammers, R.B., 2000. Global Water
Resources: Vulnerability from Climate Change et Population Growth. Science
289, 284–288. [Link]
Wackernagel, M., Schulz, N.B., Deumling, D., Linares, A.C., Jenkins, M., Kapos, V.,
Monfreda, C., Loh, J., Myers, N., Norgaard, R., Randers, J., 2002. Tracking
the ecological overshoot of the human economy. PNAS 99, 9266–9271.
[Link]
150
Wang, K., Franklin, S.E., Guo, X., He, Y., McDermid, G.J., 2009. Problems in remote
sensing of landscapes et habitats. Progress in Physical Geography: Earth et
Environment 33, 747–768. [Link]
Wang, W., Lu, H., Leung, L.R., Li, H.-Y., Zhao, J., Tian, F., Yang, K., Sothea, K., 2017.
Dam Construction in Lancang-Mekong River Basin Could Mitigate Future
Flood Risk From Warming-Induced Intensified Rainfall. Geophysical Research
Letters 44, 10,378-10,386. [Link]
Weng, Q., 2009. Remote Sensing et GIS Integration: Theories, Methods, et
Applications.
White, I., 2002. Water management in the Mekong delta: changes, conflicts et
opportunities, in: UNESCO’s International Hydrological Programme. Technical
Documents in Hydrology, No 61.
Wilder, M., Phuong, N.T., 2002. The status of aquaculture in the Mekong Delta region
of Vietnam: Sustainable production et combined farming systems. Fisheries
science 68, 847–850. [Link]
Woodcock, C.E., Strahler, A.H., 1987. The factor of scale in remote sensing. Remote
Sensing of Environment 21, 311–332. [Link]
4257(87)90015-0
Xian, G., Crane, M., 2005. Assessments of urban growth in the Tampa Bay watershed
using remote sensing data. Remote Sensing of Environment 97, 203–215.
[Link]
Yang, G., Pu, R., Zhang, J., Zhao, C., Feng, H., Wang, J., 2013. Remote sensing of
seasonal variability of fractional vegetation cover et its object-based spatial
pattern analysis over mountain areas. ISPRS Journal of Photogrammetry et
Remote Sensing 77, 79–93. [Link]
Yang, J., He, Y., Weng, Q., 2015. An Automated Method to Parameterize
Segmentation Scale by Enhancing Intrasegment Homogeneity et
Intersegment Heterogeneity. IEEE Geoscience et Remote Sensing Letters 12,
1282–1286. [Link]
Zhou, Chellappa, 1988. Computation of optical flow using a neural network, in: IEEE
1988 International Conference on Neural Networks. Presented at the IEEE
1988 International Conference on Neural Networks, pp. 71–78 vol.2.
[Link]
Zhou, H., Li, X., Schaefer, G., Celebi, M.E., Miller, P., 2013. Mean shift based gradient
vector flow for image segmentation. Computer Vision et Image Understanding
117, 1004–1016. [Link]
Zhu, X.X., Tuia, D., Mou, L., Xia, G.-S., Zhang, L., Xu, F., Fraundorfer, F., 2017. Deep
learning in remote sensing: a review. IEEE Geosci. Remote Sens. Mag. 5, 8–
36. [Link]
151
ANNEXE
152
2. Calcul du coefficient de rétrodiffusion
153
Figure 51 : Profils VH temporel des classes de l’occupation du sol en 2020.
154
3. Paramétrage des modèles d’apprentissage automatique
155
Titre : Méthodes de Machine Learning pour le suivi de l’occupation du sol des deltas du Viêt-Nam.
Mots clés : méthodes d’apprentissage automatique, delta du Vietnam, suivi de l’occupation du sol,
Title : Machine Learning methods for monitoring land use in deltas of Vietnam
Abstract : Socio-economic development in Vietnam Delta. Indeed, experimental work has been carried
is associated with the existence of large fluvial out by verifying and evaluating the contribution of
deltas. Furthermore, environmental factors such as multi-sensor image processing through various
dryness and flooding have an important role in the image segmentation approaches and machine/deep
change of land use/land cover within these deltas. learning algorithms. Thus, a Convolutional Neural
These changes have an impact on the natural and Network (CNN) model adapted to the context of the
economic balance of the country. In this perspective, study demonstrated its robustness for the detection
the objectives of the present thesis are to suggest and mapping of land use in order to characterise the
processing methods of satellite data for an efficient flood hazard and analyse the issues at risk.
mapping and monitoring of land use in the two main
deltas of Vietnam, the Red River and the Mekong
L'adoption d'un système intégré agriculture-aquaculture est cruciale pour maximiser l'utilisation des ressources limitées face à l'expansion de l'aquaculture et à la réduction des terres agricoles. Cela permet d'optimiser les rendements agricoles tout en soutenant les activités économiques essentielles dans le delta .
Le changement climatique provoque un déséquilibre dans le régime des précipitations et dans le cycle hydrologique du delta du Mékong, influençant l'évapotranspiration et le ruissellement. L'élévation du niveau de la mer et l'exploitation des barrages modifient la disponibilité en eau douce, entraînant des altérations dans l'utilisation des terres .
Le delta du Mékong fait face à de nombreux problèmes environnementaux, dont le manque d'eau douce, l'intrusion d'eau salée, l'acidification des sols, et la pollution des eaux estuariennes. Ces problèmes, aggravés par les changements climatiques et le développement socio-économique, menacent l'équilibre écologique de la région .
L'intelligence artificielle est intégrée dans le traitement des données de télédétection principalement à travers des algorithmes de machine learning et de deep learning, qui permettent d'améliorer l'analyse et l'interprétation des images satellites. Par exemple, des méthodes comme les Support Vector Machines, les réseaux neuronaux et les Random Forests sont utilisées pour des applications comme la classification et la segmentation des images . Les réseaux convolutifs neuronaux (CNN) sont particulièrement efficaces pour extraire automatiquement des caractéristiques importantes des données brutes, facilitant ainsi la détection d'occupation du sol et la cartographie d'inondations régionales . De plus, les techniques GEOBIA, qui incluent l'utilisation de méthodes d'apprentissage automatique pour l'analyse orientée objet, sont appliquées pour étudier les changements de la couverture terrestre, particulièrement dans les zones deltaïques du Vietnam . Enfin, l'intégration de ces technologies avancées permet une gestion et une planification plus efficaces des catastrophes naturelles comme les inondations, grâce à une cartographie précise et à jour .
Les algorithmes de Deep Learning, notamment CNN, ont amélioré la classification des images de télédétection en offrant une meilleure capacité d'analyse et de reconnaissance de formes complexes. Ils surpassent souvent les méthodes d'apprentissage automatique traditionnelles en termes de précision et d'efficacité, comme discuté dans diverses études comparatives .
Les changements dans les pratiques agricoles ont transformé l'occupation du sol dans le delta du Mékong principalement par le développement massif de l'aquaculture aux dépens des terres agricoles traditionnelles, notamment les rizières. L'industrie aquacole, encouragée par le gouvernement vietnamien, a poussé à la conversion des terres agricoles en zones d'aquaculture, comme l'élevage de crevettes, ce qui a entraîné une diminution significative des zones de riziculture . En plus, des conditions climatiques, telles que l'intrusion salines, ont également affecté la riziculture et d'autres cultures, ce qui a incité à adapter l'utilisation des terres . La politique et des réglementations ont également joué un rôle, cherchant à gérer l'équilibre entre l'aquaculture et la conservation des mangroves . Par ailleurs, l'urbanisation et le développement d'infrastructures ont contribué à la transformation spatiale et économique de la région, aggravant les pressions sur l'environnement et l'occupation du sol .
Les techniques de télédétection, notamment grâce aux images satellites multi spectrales, SAR et hyperspectrales, contribuent significativement à la gestion des ressources dans les deltas du Viêt-Nam en permettant de suivre et cartographier l'occupation du sol . Elles aident à évaluer les impacts des risques d'inondation sur les pratiques agricoles, particulièrement la culture du riz , et permettent de distinguer les surfaces inondées des terres agricoles, ce qui est crucial pour la gestion des crises comme les inondations . De plus, la télédétection soutient la gestion des ressources aquacoles en documentant les changements spatio-temporels des systèmes d'aquaculture intensive . Ces outils facilitent également l'adaptation au changement climatique en suivant les intrusions salines, lesquelles affectent considérablement l'environnement et l'utilisation des terres . Enfin, les avancées dans l'apprentissage automatique permettent une meilleure classification et segmentation des images satellites, offrant des cartes précises pour la gestion intégrée et durable des zones sensibles .
Les facteurs influençant l'occupation du sol dans le delta du Mékong incluent des éléments socio-économiques comme le développement de l'aquaculture et l'agriculture, principalement la riziculture, encouragés par les politiques nationales . Les transformations des terres agricoles en zones d'aquaculture ont été notables depuis 1998 . Les conditions environnementales jouent aussi un rôle, avec des phénomènes comme l'intrusion des eaux salines pendant la saison sèche, et les inondations dues aux fortes pluies du bassin versant et à l'élévation du niveau de la mer . Le changement climatique exacerbe ces effets en modifiant le régime des précipitations et la disponibilité en eau . Ces facteurs naturels et anthropiques combinés ont un impact significatif sur l'équilibre écologique et économique de la région .
Les défis liés à l'utilisation de la télédétection pour la gestion des risques naturels incluent la nécessité de développer des techniques de traitement automatique des données pour une exploitation optimale dans des contextes spécifiques comme les inondations, la sécheresse ou l'érosion . Malgré le potentiel élevé pour la prise de décisions, les données de télédétection ne sont pas toujours utilisées au maximum de leurs capacités, en partie à cause des contraintes techniques liées à l'évaluation et à la validation d'images satellites multi-capteurs , ainsi que des défis dans la modélisation et la structuration des modèles hydrodynamiques pour les inondations . La sélection précise des algorithmes d'apprentissage automatique pour différents types de données reste également complexe et requiert un ajustement minutieux des paramètres pour obtenir des résultats optimaux .
La guerre du Viêt-Nam a eu des impacts significatifs sur l'environnement du delta du Mékong, en particulier en raison des attaques chimiques massives menées par l'armée américaine. Environ 600 000 hectares de terres ont été pulvérisés avec des défoliants, détruisant plus de 20 % des zones forestières du sud du Viêt-Nam, y compris une part importante des mangroves. Plus de 2 millions d'hectares de forêts et de mangroves côtières ont été détruits, ce qui a eu un effet dévastateur sur les écosystèmes . Environ 14 % des missions chimiques visaient même les rizières et autres cultures, entraînant la perte de 300 000 tonnes de récoltes et la destruction de 236 000 hectares de terres agricoles . Ces actions ont considérablement altéré la structure écologique et la fertilité des sols dans le delta, exacerbant les défis environnementaux et économiques de la région. Enfin, l'intrusion des eaux salines, exacerbée par les altérations de l'environnement naturel, continue de poser un problème pour l'agriculture, en particulier pendant la saison sèche, affectant ainsi les rizières et réduisant par conséquent la production agricole .