3. Méthode d'analyse
géostatistique :
L'analyse géostatistiqueoffre une gamme d'outils
puissants pour l'analyse exploratoire des données spatiales et la
création des surfaces interpolées. Il permet de créer des
modèles continus à partir des données
échantillonnées et de prédire ainsi les valeurs sur des
régions où le prélèvement n'a pas pu être
effectué. De plus, l'analyse géostatistique accorde la
possibilité d'analyser les caractéristiques qualitatives et
quantitatives des données échantillonnées. (ESRI,
2010)
Nous avons fait l'analyse géostatistique pour
cartographier les données calibrées. En effet, nous avons choisi
la méthode géostatistique Krigeage ordinaire
(GAASCUEL-ODOUX, 1984 ;
VOLTZ, 1986 ; WALTER,
1990etBOIVOIN, 1991)
pour :
· construire et interpréter les
semi-variogrammes
· la représentation de la répartition de la
salinité des sols
· construiredes cartesd'isovaleurs
Nous avons réalisé l'analyse
géostatistique de données à l'aide du logiciel ARCGIS 9.3
en appliquant les étapes suivantes :
3.1. Mise en oeuvre des points d'échantillonnage.
· Ouvrir la base de données qui est
enregistrée sur Excel2010 par ARCGIS 9.3.
· Sélectionner les coordonnées X et Y (1)
et sélectionner le système de projection (2)
3
1
2
· Une fois qu'on a affiché les points, les
échantillons sont localisés (3) dans un système de
projection bien déterminé.
3.2. La cartographie de la surface de concentration des sels
dans le sol :
Nous allons créer les valeurs de conductivité
apparentepar l'interpolation. Nous allons faire la cartographie de la surface
de concentration des sels dans le sol en utilisant les paramètres de
Krigeage ordinaire.
4
3
2
1
· Nous cliquons sur « Geostatistical
Wizard » (1) pour commencer l'analyse géostatistique.
· Une boite de dialogue (2) qui sera affichée et
qui demandera la méthode d'interpolation et l'emplacement de
données qu'elles vont être interpolées. Nous choisissons la
méthode Krigeage et nous importons les valeurs de conductivité
électrique apparente et nous mettons en place les coordonnées X
et Y des données.
· Nous sélectionnons dans la boite de dialogue
« méthode géostatistique » la carte de
prédiction et nous choisissons le second ordre de suppression de la
tendance (3). Un polynôme de second degré sera adapté pour
tracerla courbe qui va préciser lesdirections géographiques de la
carte.
· La boite de dialogue (4) montre l'élimination de
la tendance en option standard.
8
7
6
5
a
b
· Cette étape est considérée comme
la plus importante puisque nous allons faire la modélisation
dusemivariogramme (5). Le nuage des points présentés dans
lesemivariogramme doit être corrélé en faisant cocher
anisotropy (a) ensuite nous modifions la direction de l'angle (b) pour
atteindre un maximum de points proches de la courbe modèle
· Dans la boite de dialogue
« Searchingneighborhood » (6) nous constatonsque les points
d'échantillonnages ainsi qu'une ellipse dans un emplacement par
défaut. Il va falloir orienter l'angle de direction comme il est
corrélé dans (1, b) lesemivariogramme
9
· L'étape (7) de la validation croisée
(Cross-validation) sert à fournir les prévisions les plus
précises pour indiquer si le modèle est raisonnable pour la
production d'une carte. En cliquant sur QQPlot (8), nous pouvons voir que
certaines valeurs tombent légèrement au-dessus et au dessous de
la ligne, mais la plupart des points très proche de la ligne droite en
pointillé, ce qui indique que les erreurs de prédiction sont
prêtesà être normalement distribuées.
· Une fois le taux d'erreur est très proche de 0,
nous pouvons passer à la dernière étape celle d'observer
le sommaire de la méthode (9) et la création de la carte de
répartition de la salinité des sols.
| 
