17.16. Analyse hydrologique¶

Dans cette leçon, nous allons faire une analyse hydrographique. En commençant avec un MNE, nous allons extraire un réseau de canaux, délimiter des bassins versants et calculer des statistiques.

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

The first module to execute is Catchment area. You can use anyone of the others named Catchment area. They have different algorithms underneath, but the results are basically the same.

Sélectionnez le MNE dans le champ Élévation, et laissez les valeurs par défaut pour le reste des paramètres.

Certains algorithmes calculent plusieurs couchent, mais celui de la Surface de bassin versant est le seul que nous utiliserons.

Vous pouvez vous débarrasser des autres si vous voulez.

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

Pour savoir pourquoi, vous pouvez regarder l’histogramme et vous verrez que des valeurs ne sont pas distribuées uniformément (il y a quelques cellules avec une très grande valeur, celles qui correspondent au réseau de canaux). Le calcul du logarithme de la valeur du bassin versant produit une couche qui transmet beaucoup plus d’information (vous pouvez le faire en utilisant la calculatrice raster).

Le basin versant (aussi connu comme l’accumulation de flux) peut être utilisé pour mettre une limite pour l’initiation de canaux. Voici comment vous devez le mettre en place (notez le Seuil d’initiation Plus grand que 10.000.000).

Utilisez la couche originale bassin versant, pas celle du logarithme. Celle-ci était juste à des fins de rendu.

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

L’image ci-dessus montre simplement la couche vectorielle de résultat et le MNE, mais il devrait aussi y avoir un raster avec le même réseau de canaux. Ce raster-là sera en fait celui que nous utiliserons.

Maintenant, nous utiliserons l’algorithme Basins versants pour délimiter les sous-bassins correspondant à ce réseau de canaux, en utilisant comme points de sortie toutes les jonctions. Voici comment vous devez configurer la boîte de dialogue des paramètres correspondants.

Et voici ce que vous obtiendrez.

Il en résulte un raster. Vous pouvez le vectoriser en utilisant l’algorithme Vectoriser les classes de grille.

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.

First, let’s clip the original DEM with the polygon representing a subbasin. We will use the Clip grid with polygon algorithm. If we select a single subbasin polygon and then call the clipping algorithm, we can clip the DEM to the area covered by that polygon, since the algorithm is aware of the selection.

Sélectionnez un polygone,

et appelez l’algorithme de découpage avec les paramètres suivants :

L’élément sélectionné dans le champ d’entrée est, bien sûr, le MNE que nous voulons découper.

Vous obtiendrez quelque chose comme ça.

Cette couche est prête à être utilisée dans l’algorithme Statistiques de la couche raster.

Les résultats des statistiques sont les suivants.

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.