Como encontrar a área da bacia hidrográfica para cada ponto da rede fluvial usando DEM？

Atualização # 1 : Obrigado, Micha, segui suas instruções e utilizei com sucesso os módulos no GRASS.

Eu uso r.watershed para obter a varredura da bacia e a linha de drenagem:

e eu uso r.water.outlet para obter o divisor de águas de uma saída:

Eu tenho toda a área da bacia do ponto de saída verde, mas tudo o que eu quero é a pequena parte da área original da bacia, como esta: parte B é o que eu quero

Esta foto esquemática mostra a ideia do que eu quero:

Posso atingir esse objetivo usando r.water.outlet ou outros módulos no GRASS (ou tudo o que posso fazer é usar a ferramenta de interseção nas ferramentas de geoprocessamento no QGIS)

Como título, quero encontrar a área da bacia hidrográfica, mas usando o QGIS.

Tenho os seguintes arquivos para trabalhar:

1. shapefile de rede de fluxo (polilinha)
2. shapefile de ponto recortado do shapefile de polilinha
3. uma ilha inteira (o fluxo está dentro da ilha) DEM (eu tenho 2 tipos: .lan & .tiff)

Quero encontrar a área total da bacia hidrográfica localizada a montante de todos os pontos da rede, aqui está o diagrama esquemático:

as linhas verde, vermelha e amarela delimitam a bacia hidrográfica de cada ponto (aqui eu só uso 3 pontos como exemplo) e quero calcular sua área.

Já obtive o arquivo DEM e devo cortar a região que contém o fluxo, devo poligonizá- lo primeiro？ Posso usar diretamente este formulário (.lan ou .tiff) para trabalhar com GRASS？

Definitivamente, você precisará de um raster de elevação (DEM) para fazer o que quiser. Com o DEM, você pode usar dois módulos GRASS para obter áreas de drenagem individuais para pontos ao longo da rede fluvial, como segue: Primeiro calcule uma grade de direção do fluxo do DEM com r.watershed:

r.watershed elev=DEM thresh=<your threshold> drain=flow_dir_grid

Veja o manual r.watershed para obter mais detalhes Próximo loop através dos valores XY dos pontos ao longo da rede fluvial para calcular cada bacia hidrográfica com o módulo r.water.outlet.

r.water.outlet drain=flow_dir_grid basin=basin1 easting=X1 northing=Y1

Detalhes na página de manual do r.water.outlet

Para fazer o loop, você poderá obter algumas idéias nesta postagem