Méthodes d'extraction
Quel que soit le logiciel SIG, le principe d'extraction est le même que celui détaillé dans la vidéo suivante.
Simulation : Etape 1
Sur la base des exemples fournis ci-dessous (vidéos et captures d'écran) sur le bassin du Giessen (Alsace, France), appliquer successivement les étapes clés sur le bassin versant de la Doller (en téléchargeant les bases de données à la fois sur la Doller et le Giessen).
Etape 1 : La première étape d'extraction automatique consiste à combler les zones endoréiques, c'est-à-dire les pixels dont l'altitude est plus basse que les 8 pixels voisins. Ces zones provoquent l'arrêt des calculs de propagation amont-aval et ne permettent pas d'extraire correctement les réseaux.
La vidéo suivante illustre l'outil de calcul de Grass dans QGIS r.fill.dir permettant de combler les zones endoréiques :
Conseil :
L'exécution de l'outil r.fill.dir peut ne pas aboutir au comblement de toutes les zones endoréiques comme le montre cette illustration : |  |
Il faut alors relancer l'outil à partir de la couche obtenue et ainsi de suite jusqu'à résolution (c'est-à-dire comblement) de toutes les zones endoréiques : |  |
Simulation : Etape 2
Etape 2 : Une fois le mnt traité (c'est-à-dire ses zones endoréiques comblées), il est nécessaire d'utiliser la fonction r.watershed illustrée ci-dessous :
La couche d'entrée (mntbis) a été comblée, le seuil de 400 correspond au nombre de pixels amont pour extraire le début du réseau (la notion de seuil sera rediscuté dans la section sur le réseau hydrographique).
Les couches de sortie sont :
|  |
|  |
|  |
|  |
Simulation : Etape 3
Etape 3 : Il est ensuite possible d'utiliser la fonction r.water.outlet pour extraire la zone amont d'un exutoire dont vous devez définir les coordonnées (x, y) comme illustrée dans la figure suivante sur le bassin du Giessen (Alsace, France). Il est nécessaire d'utiliser le raster dirflux400 créé avec le précédent outil (r.watershed)
Simulation : Etape 4
Etape 4 : La couche raster ainsi obtenue peut être transformée en vecteur pour des traitements spatiaux ultérieurs ou pour créer des cartes.
Cette opération s'appelle une vectorisation.
La vidéo suivante illustre la conversion du bassin versant du Giessen du format raster au format vecteur :
Il est possible au format vecteur de calculer des champs (ou colonnes) des propriétés géométriques clés comme la superficie et le périmètre comme illustré dans la vidéo suivante sur le bassin du Giessen.
