3.4.3.5 Calcul de l'Indice de Salinité par Difference
Normalisée (NDSI)
Pour détecter la salinité à la surface du
sol, nous allons calculer un certain nombre d'indices mathématiques sur
les images satellitaires. Ces indices résultent de la combinaison de
plusieurs bandes de même spectre et qui ont été
développés par Abbas et Khan en 2007. Ils sont fondés sur
l'hypothèse que cette transformation rendrait l'information
demandée de salinité plus proéminente tout en supprimant
les effets d'autres utilisation et occupation des sols (Abbas et al.,
2013).L'indice utilisé dans le cadre de cette étude
est celle de salinité par différence normalisée
(NDSI). Ilest appliqué sur les images de Landsat 8.Il exploite
essentiellement la différence de réponse spectrale de la
végétation et des sols dans la bande rouge ® qui est
liée à l'absorption de la lumière par la chlorophylle et
la bande proche infrarouge (PIR) qui est liée à la densité
de la végétation verte. Cet indice permet de différencier
les zones de salinité, il est calculé par le PIR
et le Rougesous la formuleNDSI : (RED-PIR) /
RED+PIR).
Landsat 8 : NDSI = (B4 - B5) / (B4 +
B5)
Avec B4 la bande Rouge et B5 la bande Proche Infrarouge.
Les valeurs obtenus sont des nombres réels allant de -1
à 1.
3.4.3.6 Cartographie des zonesaffectées par la
dynamique saline
Il s'agit de l'identification et de l'estimation des zones
vulnérables face aux terres salines.
? Méthodologie de classification des zones
vulnérables par analyse de proximité
La répartition de la population en 2023, la
distribution des terres salées constante et apparue en 2024 sont
combinées afin de cartographier la population touchée par la
dynamique des terres salées. A partir de la sélection par
emplacement `select by location' par rapport aux terres salées, nous
avons ressortir cinq (5) catégories de populations de la plus
affectée à la moins touchée.
Tableau 10:
Classification de la population selon le degré d'exposition aux
terres salines
|
Catégories populations
|
Emplacement/Terres salées
|
Durée séjour dans terres
salées
|
|
Population extrêmement affectée
|
Intérieur des tannes constantes
|
52 ans
|
|
Population affectée
|
Intérieur des tannes apparues
|
-30 ans
|
|
Population à risque forte
|
A 50 m des tannes récentes
|
0
|
|
Population à risque faible
|
A 100 m des tannes récentes
|
0
|
|
Population sans risque
|
A plus 100 m des tannes récentes
|
0
|
Source : Diouf S. Aziz, Juin
2025
? Méthodologie de classification des zones
vulnérables par analyse multicritère
La cartographie des zones exposées à la
dynamique des terres salines a été obtenue par le processus
d'analyse hiérarchique (AHP). En effet, il s'agit d'une analyse
multicritère basée sur la pondération des facteurs
contribuant aux risques de la salinisation (pluviométrie, topographie,
le sol, texture, etc.). La méthode d'analyse hiérarchique
multicritère (AHP) développée par le mathématicien
Thomas Saaty (1977) est utilisée pour cartographier la
vulnérabilité liée aux terres salées dans la
Commune de Gandon. Sa méthode repose sur certains principes dont la
structuration hiérarchique selon les classes, les critères et
poids définis, la structuration des priorités selon les sous
critères, les rangs et la cohérence logique selon l'indice de
cohérence (IC) et le ratio de cohérence (RC) de la matrice de
comparaison par paire. Elle consiste à comparer les critères deux
à deux en termes d'importance relative par rapport à l'objectif
défini sur la base d'une échelle de pondération. Cette
échelle d'évaluation mesurant l'importance d'un
élément sur un autre.
Tableau 11 :
Echelle de pondération des couches
|
Comparaison d'un critère par rapport à un autre
|
Intensité de l'importance
|
|
Même importance que
|
1
|
|
Modérément plus important que
|
3
|
|
Fortement plus important que
|
5
|
|
Très fortement plus important que
|
7
|
|
Extrêmement plus important que
|
9
|
|
Modérément moins important que
|
1/3
|
|
Fortement moins important que
|
1/5
|
|
Très fortement moins important que
|
1/7
|
|
Extrêmement moins important que
|
1/9
|
Il convient d'abord de définir les facteurs relatifs
à l'extension des terres salées qui constituent en effet les
critères à prendre en considération. Les critères
choisis sont : la pluviométrie, la topographie, le pente, les types de
sol, la texture, les tannes, la densité de la population et le
réseau hydrographique. En fonction de l'échelle numérique,
les facteurs cités ci-dessus sont pondérés par comparaison
binaire.
Tableau 12 :
Matrice de pondération des facteurs selon l'échelle de
Saaty (1977)
|
Facteurs
|
P(mm)
|
Réseau Hydrographie
|
Texture
|
Sol
|
Topographie
|
Pente
|
Densité population
|
|
Pmm
|
1
|
1
|
3
|
3
|
5
|
7
|
9
|
|
Réseau Hydrographie
|
1
|
1
|
3
|
3
|
5
|
5
|
7
|
|
Texture
|
1/3
|
1/3
|
1
|
3
|
5
|
7
|
7
|
|
Sol
|
1/3
|
1/3
|
1/3
|
1
|
3
|
3
|
5
|
|
Topographie
|
1/5
|
1/5
|
1/5
|
1/3
|
1
|
3
|
5
|
|
Pente
|
1/7
|
1/5
|
1/7
|
1/3
|
1/3
|
1
|
3
|
|
Densité population
|
1/9
|
1/7
|
1/7
|
1/5
|
1/5
|
1/3
|
1
|
Le tableau ci-dessus nous a permis de définir les
vecteurs propres pour chaque facteur. Par la suite, le total des lignes pour
chaque facteur est calculé ainsi que la somme totale des lignes avant de
calculer les poids selon la formule: Poids du critère = Total ligne
facteur / Somme des lignes.
Tableau 13:
Calcul du poids relatifs aux différents facteurs après
comparaison deux à deux.
|
Facteurs
|
P(mm)
|
Réseau Hydrographie
|
Texture
|
Sol
|
Topographie
|
Pente
|
Densité population
|
Total Ligne
|
Poids
|
Poids %
|
|
Pmm
|
1
|
1
|
3
|
3
|
5
|
7
|
9
|
29,00
|
0,27
|
27
|
|
Réseau Hydrographie
|
1
|
1
|
3
|
3
|
5
|
5
|
7
|
25,00
|
0,23
|
23
|
|
Texture
|
1/3
|
1/3
|
1
|
3
|
5
|
7
|
7
|
23,67
|
0,22
|
22
|
|
Sol
|
1/3
|
1/3
|
1/3
|
1
|
3
|
3
|
5
|
13,00
|
0,12
|
12
|
|
Topographie
|
1/5
|
1/5
|
1/5
|
1/3
|
1
|
3
|
5
|
9,93
|
0,09
|
9
|
|
Pente
|
1/7
|
1/5
|
1/7
|
1/3
|
1/3
|
1
|
3
|
5,15
|
0,05
|
5
|
|
Densité population
|
1/9
|
1/7
|
1/7
|
1/5
|
1/5
|
1/3
|
1
|
2,13
|
0,02
|
2
|
|
Total
|
|
107,88
|
1,00
|
100
|
L'utilisation d'un SIG permet de spatialiser la
pondération de ces facteurs afin de cartographier les zones
vulnérables. Les facteurs convertis en format raster sont reclassifier.
Ensuite des indices statistiques ont été obtenue suite à
la combinaison de ces différents facteurs effectués avec l'outil
"Overlay_Weighted Overlay" du logiciel ArcGIS en attribuant chaque facteur son
poids en %.
Figure 18 :
Traitement par AHP dans Arcgis
Ce traitement multiple permet d'obtenir des intervalles
d'indices auxquels des classes de vulnérabilité sont
attribuées. Les résultats de l'analyse d'exposition aux terres
salées obtenue par l'AHP mettent en avant 4 classes hiérarchiques
de susceptibilité (très faible, faible, fort, très fort)
sur les facteurs contribuant aux risques d'atteintes par la salinité.Le
résultat est enfin converti en vecteur afin de calculer la part de
chaque classe en ha et sa représentation sous forme graphique.
DonnéesDémographiquesClimatiques
Images Satellitaire Landsat
Base de Données Spatiale
Localisation/Pédologie/Occupation sol
Correction, Calcul d'indices, classification
BV, Relief, Pente
Reclassification
QGIS
ArcGIS
Validation
ArcGIS
CE/pH/Texture (ISRIC)
Données vectorielles
Excel, ArcGIS
ArcGIS
ENVI
Image SRTM MNT
Calcul statistique
QGIS, Excel
ENVI
Classes thématiques
MCD, MLD, MPD
Analyse et visualisation
Cartes
Graphiques
Figure
19 :Schéma méthodologique
| 
