Tutoriel: Calculer l’indice d’humidité topographique à partir d’un MNT

2. Interpoler les lacunes dans le MNT

Le MDT fourni est dérivé des données LIDAR des Pays-Bas (AHN3). Les données de la projection néerlandaise d’Amersfoort (EPSG : 28992). Le MDT, cependant, a des vides. Les vides sont des pixels sans données. Ces vides doivent être interpolés avant que nous puissions continuer.

1. Lance QGIS.

2. Ajoutez le MDT fourni (dtm.tif) à un projet vide.

Vous pouvez clairement voir les vides/lacunes.

voids

3. Accédez à la boîte à outils traitement : dans le menu principal, choisissez Processing | Toolbox.

4. Dans la boîte à outils de traitement, sélectionnez GDAL | Fill nodata.

5. Dans la boîte de dialogue ‘Fill nodata’ choisir dtm comme couche d’entré ‘Input layer’. Conservez les autres paramètres par défaut (vous pouvez jouer avec ceux-ci si les vides sont plus grands). Enregistrez le résultat en tant que ‘dtm_voidfilled.tif’.

6. Cliquez sur Exécuter ‘run’. Fermez le dialogue après traitement.

7. Retirez la couche dtm du panneau des couches. 

Dans le résultat, vous pouvez clairement voir que les vides ont été interpolés.

voids interpolated

8. Stylez le DTM avec un pseudocolore à bande unique. Vous pouvez utiliser une rampe de topographie du catalogue cpt-city. Mélangez le DTM avec le ‘hillshade’ (voir d’autres tutoriels sur la façon de le faire).

Styled DTM

Maintenant que le DTM est prêt, nous pouvons procéder au calcul de la pente.

Notez que dans QGIS 3.14 vous devez définir la projection de la couche raster après chaque calcul. Ce n’est pas le cas dans les versions LTR 3.10.