Faire un semi variograme et une carte krigée avec surfer
Jérôme Mathieu
http://www.jerome.mathieu.freesurf.fr
avril 2004
Fichier de données
Faire un fichier excel avec les données organisée en colonnes: une colonne pour les X, une
pour les Y, une par variable à étudier.
La première ligne contient le nom des variables. Les virgules n'ont pas besoin d'être
transformées en ".".
Une fois les données enregistrées, lancer surfer.
Réalisation du semi- variograme dans surfer
Ouvrir surfer. Une feuille "plot" est créée. Dans le cas contraire la créer.
-1-
Ensuite choisir l'option Grid Variogram New Variogram
Et selectionner le fichier de données créé en excel.
Une fois le fichier de données sélectionné, une boîte de dialogue s'ouvre, pour paramétrer le
variogramme:
Vérifier que les variables sélectionnées correspondent bien à celle désirées. Laisser les
paramètres duplicates tels quels.
-2-
Dans l'onglet option, ajuster Max Lag Distance à la distance de prospection maximale du
semivariogramme. (Se traduit par la valeur maximale en abscisse dans le semi variogramme).
En général on rentre la moitié de la taille du carré d'étude.
Il est important de faire ce réglage dès le départ, sinon le variograme sera mal conditionné par
la suite.
Après avoir sélectionné OK, le variograme est créé, ajuté par un modèle linéaire de base, qu'il
va falloir améliorer.
Column C: temp
Direction: 0.0 Tolerance: 90.0
0.6
0.5
0.4
Variogram
0.3
0.2
0.1
0
0 10 20 30 40 50 60 70 80 90 100
Lag Distance
Variogramme brut
-3-
Ajuster le variogramme
En premier lieu il faut ajuster le nombre d'intervales et/ou leur taille. Il est intuitif de prendre
comme intervale ("lag)" la distance réelle entre les points d'échantillonnage
Pour cela faire un double clic sur le variogramme et une fenêtre s'ouvre:
On peut ici ajuster ces paramètres:
- Max Lag distance: distance maximale de prospection du semi variogramme (X max)
- Number of lags : nombre d'intervales d'études
- Lag width : taille des intervales (se règle automatiquement en fonction du nombre
d'intervales)
- Vertical scale: Y max du semi variogramme.
Lag direction:
permet de trouver d'éventuelles anisotropies. A utiliser une fois que le modèle a été défini.
Une fois que ces paramètres sont modifiés, il faut améliorer le modèle d'ajustement du
variogramme.
-4-
Ajuster le modèle
Une fois ces paramètres ajustés, on peut ajuster le modèle (courbe bleue). Pour cela
sélectionner l'onglet Model:
La partie Variogram Components décrit le modèle utilisé.
Par défaut Surfer commence avec un effet pépite 'nugget effect "(C0), et un modèle linéaire.
Très souvent ce modèle ne correspond pas et il faudra donc changer ces composants.
Pour cela, séléctionner le composant à enlever et selectionner remove.
Pour ajouter un nouveau composant, appuyer sur add, une fenêtre s'ouvre, avec différents
composants ajoutables:
Choisir le type de modèle désiré. (On prend souvent sphérique ou gaussien).
Une fois OK sélectionné, l'onglet Model est modifié en fonction du choix.
-5-
Paramétrer le modèle du variogramme
Il convient alors de paramétrer le modèle pour qu'il épouse la forme du semi variograme.
Personnellement je recommanderai plutôt l’utilisation du logiciel gstat pour cette étape.
Dans surfer, l'option autofit est possible, mais les résultats ne sont pas toujours très
convaincants. Si les résultats ne sont pas satisfaisant, il vaut mieux ajuster manuellement les
paramètres, en jouant sur les réglages proposés dans la boite d'affichage du modèle.
Cette boîte dépend du modèle choisit.
En premier lieu, vérifier que le composant à paramétrer est sélectionné (affiché en bleu).
Ensuite jouer sur les paramètres suivants:
- Scale : correspond à la valeur de la semivarinace où le plateau est sensé être atteint (hauteur)
- Length (A) : correspond à la valeur de x où le plateau est atteint: c'est la "range", ou la
distance maximale sur laquelle les données présentes une dépendance spatiale.
- Ratio et Angle permettent d'introduire de l'anisotropie (souvent inutile)
A chaque modification on peut visualiser le résultat en appuyant sur appliquer. Petit à petit
on finit par avoir un modèle reflétant le variogramme. (pour le moment je ne connais pas de
fonction dans Surfer permettant de tester l'ajustement du modèle au semivariogramme.
Column C: temp
Direction: 0.0 Tolerance: 90.0
0.55
0.5
0.45
0.4
0.35
Variogram
0.3
0.25
scale
0.2
0.15
length
0.1
0.05
0
0 10 20 30 40 50 60 70 80 90 100
Lag Distance
A cette étape le variogramme est fini, et peut être utilisé pour kriger les données, c'est à dire
pour faire une carte des valeurs.
-6-
Faire une carte krigée des données
Une fois le variogramme établis, et son modèle ajusté, on peut faire une carte des valeurs par
kriegeage, c'est à dire par interpolation des valeurs entre les points de mesure.
Il faut pour cela récupérer les paramètres du semivariogramme, puis créer les données
interpolées. Ensuite on peut afficher les données sous forme de carte.
Récupérér les paramètre du variogramme pour kriger
Sélectionner le variogramme dans le plot, puis choisir l'option :
Grid Variogram export Variogram
L'enregistrer sous un nom explicite
Ensuite il faut créér les donnnées interpolées selon les paramètres du variogramme modélisé.
-7-
Créér les données interpolées
choisir l'option Grid data
sélectionner le ficher excel de données brutes
Une fenêtre de dialogue s'ouvre alors:
-8-
La partie Grid line geometry décrit le fichier de données qui va être créé (données
interpolées). Normallement X direction et Y direction sont correctement réglés, et
correspondent aux coordonnées des points extrêmes de la carte à Krieger. Le nombre de
lignes et colonnes modifie la résolution de la carte (à vérifier).
Une fois que tout cela est bien réglé, il faut vérifier que les données qui vont être interpolées
sont bien celles qui nous interessent. Pour cela dans l'onglet data, vérifier les "data columns"
Ensuite il faut injecter les paramètres du modèle du variogramme pour faire l'interpolation:
Injection des paramètres du modèle du variogramme pour l'interpolation
des données
Toujours dans la fenetre Scattered Data Interpolation , (option Grid datat)
sélectionner le bouton Options de l'onglet "General". Cette boîte s'affiche :
Sélectionner Get Variogram
Le logiciel récupère alors automatiquement les paramètres du variogramme actif (il faut
vérifier que le bon est sélectionné quand il y a plusieurs variogrammes déjà crées!!).
Les paramètres s'affichent automatiquement.
Sélectionner OK.
Les données sont automatiquement générées et enregistrées dans un fichier .grd du meme
nom que le fichier de données excel.
-9-
Faire la carte krigée
Une fois le fichier de données interpolées créé (.grd), il est possible de réaliser la carte kriegée
des variables.
Pour cela selectionner Map Contour Map New Contour Map
Sélectionner le fichier de données interpolées (.grd) et faire les réglages classiques pour tout
graph.
Presser OK, la carte est réalisée. Pour l'editer, faire un double click dessus.
90
80
30.6
70
30.4
30.2
60 30
29.8
29.6
50
29.4
29.2
40
29
28.8
30 28.6
28.4
28.2
20
10
10 20 30 40 50 60 70 80 90
Carte kriégée finie
- 10 -
