17.16. Analyse hydrologique

Note

Dans cette leçon nous allons réaliser une analyse hydrographique. Cette analyse sera utilisée dans certaines des leçons suivantes, car elle constitue un très bon exemple d’un flux de travail d’analyse, et nous l’utiliserons pour démontrer certaines fonctionnalités avancées.

Objectifs : En partant d’un MNT, nous allons extraire un réseau de canaux, délimiter les bassins versants et calculer quelques statistiques.

  1. La première chose à faire est de charger le projet avec les données de la leçon qui contient uniquement un MNE.

    ../../../_images/dem.png
  2. Le premier module à exécuter est Zone de captage (dans certaines versions de SAGA, il est appelé Accumulation de flux (Top Down)). Vous pouvez utiliser n’importe lequel des autres modules nommés Zone de captage. Ils ont des algorithmes différents, mais les résultats sont fondamentalement les mêmes.

  3. Sélectionnez le MNT dans le champ Elevation, et laissez les valeurs par défaut pour le reste des paramètres.

    ../../../_images/catchmentarea.png

    Certains algorithmes calculent plusieurs couches, mais la couche Bassin versant est la seule que nous utiliserons. Vous pouvez vous débarrasser des autres couches si vous le souhaitez.

    Le rendu de la couche n’est pas très instructif.

    ../../../_images/catchmentlayer.png

    Pour savoir pourquoi, vous pouvez jeter un coup d’œil à l’histogramme et vous verrez que les valeurs ne sont pas distribuées de manière uniforme (il y a quelques cellules avec une valeur très élevée, celles qui correspondent au réseau de canaux). Utilisez l’algorithme Raster calculator pour calculer le logarithme de la surface de la valeur du bassin versant et vous obtiendrez une couche avec beaucoup plus d’informations.

    ../../../_images/catchmentlayerlog.png
  4. La zone de captage (également connue sous le nom d’accumulation de flux) peut être utilisée pour fixer un seuil pour l’initialisation du canal. Ceci peut être fait en utilisant l’algorithme Réseau de canaux.

    • Grille d’initialisation : utilisez la couche du bassin versant et non celle du logarithme.

    • Seuil d’initiation : 10.000.000

    • Type d’initialisation : Supérieur à

    ../../../_images/channelnetwork.png

    Si vous augmentez la valeur Seuil d’initialisation, vous obtiendrez un réseau de canaux plus clairsemé. Si vous la diminuez, vous obtiendrez un réseau plus dense. Avec la valeur proposée, voici ce que vous obtenez.

    ../../../_images/channelnetworklayer.png

    L’image ci-dessus montre juste la couche vecteur résultante et le MNT, mais il devrait y avoir aussi une couche raster avec le même réseau de canaux. Cette couche raster sera, en fait, la couche que nous allons utiliser.

  5. Maintenant, nous allons utiliser l’algorithme bassins versants pour délimiter les sous-bassins correspondant à ce réseau de canaux, en utilisant comme points de sortie toutes les jonctions qu’il contient. Voici comment vous devez définir le dialogue des paramètres correspondants.

    ../../../_images/watersheds.png

    Et voici ce que vous obtiendrez.

    ../../../_images/watershedslayer.png
  6. Il s’agit d’un résultat de type raster. Vous pouvez le vectoriser en utilisant l’algorithme Vectorisation des classes de grille.

    ../../../_images/vectorising.png
    ../../../_images/watershedslayervector.png

Maintenant, essayons de calculer des statistiques sur les valeurs d’élévation dans un des sous-bassins. L’idée est d’avoir une couche qui représente uniquement l’élévation au sein de ce sous-bassin et ensuite de la passer au module qui calcule ces statistiques.

  1. Tout d’abord, nous allons découper le MNT original avec le polygone représentant un sous-bassin. Nous allons utiliser l’algorithme découper un raster avec un polygone. Si nous sélectionnons un seul polygone de sous-bassin et appelons ensuite l’algorithme de découpage, nous pouvons découper le MNT dans la zone couverte par ce polygone, puisque l’algorithme est conscient de la sélection.

    1. Sélectionnez un polygone

      ../../../_images/selectedpolygon.png
    2. Appelez l’algorithme d’écrêtage avec les paramètres suivants :

      ../../../_images/clipgrid.png

      L’élément sélectionné dans le champ de saisie est, bien entendu, le MNT que nous voulons découper.

      Vous obtiendrez quelque chose comme ça.

      ../../../_images/clippeddem.png
  2. Cette couche est prête à être utilisée dans l’algorithme Statistiques des couches raster.

    ../../../_images/rasterstats.png

    Les résultats des statistiques sont les suivants.

    ../../../_images/stats.png

Nous utiliserons tant la procédure de calculs de bassin que le calcul statistique dans d’autres leçons pour découvrir comment d’autres éléments peuvent nous aider à les automatiser les deux et à travailler plus efficacement.