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

3. Calculer la pente en degré et en radian

Pour calculer TWI nous avons besoin de la pente dans les radians. Avec les outils d’analyse raster, nous pouvons calculer la pente en degrés et en pourcentages. Donc, nous allons d’abord calculer la pente en degrés et ensuite convertir les valeurs en radians.

1. Dans le menu principal, allez à Raster | Analysis | Slope...


2. Dans la boîte de dialogue Pente, choisi dtm_voidfilled comme couche d’entrée ‘Input layer’. Cocher la case d’utilisation de la formule ZevenbergenThorne au lieu de celui de Horn. La formule ZevenbergenThorne marche mieux avec moins de relief, lors que nous utilisons celle de Horn pour les zones raides. Choisir slopedegrees.tif comme résultat de la couche de Pente.



3. Cliquez sur Exécuter. Fermez le dialogue après traitement.

Parce que plus tard, nous allons calculer le logarithme naturel (ln) nous devons vérifier s’il ya des pixels avec une pente de 0 degrés.

4. Cliquez à droite sur slopedegrees et choisissez Propriétés.

5. Dans la boîte de dialogue Propriétés de la couche, accédez à l’onglet Informations et vérifiez la valeur minimum (STATISTICS_MINIMUM).



Dans ce cas, la pente minimale est supérieure à 0, donc nous pouvons procéder. Toutefois, si la pente du minumum est de 0, vous devez effectuer les étapes suivantes, que vous pouvez également appliquer à ces données, car il n’aura aucun effet.

6. Cliquez sur Annuler pour fermer la boîte de dialogue Propriétés de la couche.

7. Dans le menu principal, aller à ‘Raster | Raster Calculator’.



8. Dans la calculatrice Raster, créez l’équation suivante:

( "slopedegrees@1" <= 0 ) * 1 + ( "slopedegrees@1" > 0 ) * "slopedegrees@1"

Cela signifie: si les cellules ont une pente inférieure ou égale à 0, alors donner la valeur 1 (True) et multiplier par 1. Pour les cellules qui sont plus grandes que 0, assigner la valeur 1 (True) et multiplier avec sa valeur de pente. Le résultat sera 1 pour les cellules <= 0 et la pente d’origine pour les cellules 0.



9. Enregistrez le résultat sur slopedegrees_modified.tif et cliquez sur OK.

Pour passer des degrés aux radians, nous devons multiplier les valeurs de pixels par 0.01745

 

10. Retournez à la calculatrice Raster et calculez l’équation suivante :

"slopedegrees_modified@1" * 0.01745


11. Enregistrez le résultat sur sloperadians.tif et cliquez sur OK.

12. Styler le résultat du jaune au rouge.

13. Enlever les autres couches de pente.


Slope in radians

L’étape suivante consiste à calculer la zone de pente amont contributive pour chaque pixel du MDT.