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.

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.

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.