Tutoriel: Créer une carte du niveau d’eaux souterraines à partir de données de forage et d’un MNT
13. Interpolation des niveaux de la nappe en raster
Nous avons maintenant les niveaux de la nappe dans les forages. Nous pouvons interpoler les points en raster.
Nous allons utiliser deux méthodes : pondération de la distance inverse (IDW) et les Polygones de Thiessen.Commençons par l’interpolation par IDW
1. Dans la boite à outils de traitement, choisir Interpolation | IDW interpolation.
2. Dans la fenêtre de dialogue de IDW Interpolation
choisir la couche Boreholes_Z.shp couche vecteur. Utiliser GW_level_m
comme attribut
de l’interpolation. Cliquer pour ajouter ces réglages à la liste. Garder le
coefficient de distance P à 2 pour une fonction de décroissance exponentielle.
Sélectionner l’étendu (extent) à partir de la couche Stampriet study area.
Changer la taille des pixels à 30 m. Enregistrer le raster interpolé dans le
dossier du projet en tant que GW_level_IDW.tif.
3. Cliquer sur Exécuter (run). Cliquer sur Close à la fin du traitement.
Stylons le résultat avec un rampe de couleur.
4. Sélectionner la couche GW_level_IDW layer dans le panneau des couches et cliquer pour ouvrir le panneau pour styler les couches (Styling panel)5. Dans le panneau pour styler les couches changer le renderer à Singleband pseudocolor et choisir All Color Ramps | RdBu.
Le panneau de stylisation de la couche devrait ressembler à ceci.
Et la couche interpolée devrait ressembler à la figure
ci-dessous. Vous pour essayer plusieurs rampes de couleur et classes pour
améliorer la visualisation. Vous pouvez voir que l’interpolateur IDW créé une
sortie lisse, mais comme le nombre de forages est faible, il y a de nombreux biais
en direction des forage individuels et les aberrations.
Appliquons à présent les polygones de Thiessen (au connu comme voisin le plus proche ou tesselation de Voronoi).
6. Dans le menu principal choisir Raster | Analysis | Grid (Nearest Neighbor).
Dans la grille de dialogue (plus proche voisin)
choisir la couche Boreholes_Z.shp comme couche de points. Sous Advanced
(avancée) choisir GW_level_m comme champ de valeur Z. Enregistrer la sortie
interpolée en tant que GW_level_Thiessen.tif. Garder le reste par défaut.
8. Cliquer sur Exécuter (run). Cliquer sur close pour fermer la fenêtre de dialogue.
9. Styler le résultat.
Cet algorithme ne couvre pas toute la zone d’étude mais se limite à qu’on appelle enveloppe convexe de la couche d’entrée des points. Il donne un résultat moins lisse comparé à l’interpolation IDW. Cependant, avec des données rares, les plus proches valeurs peuvent être la meilleure approximation.
Nous pouvons terminer le projet en stylant les puits et exportant toutes les couches que nous voulons garder dans le GeoPackage.