Tutorial: Delimitação de Rios e Bacias Hidrográficas

Agora vamos usar essas coordenadas para calcular a área a montante (ou seja, bacia) que produz descarga neste ponto. 

1. Na caixa de ferramentas de processamento pesquise por Upslope Area e escolha SAGA | Terrain Analysis - Hydrology | Upslope Area.

2. Preencha a caixa de diálogo como na figura abaixo. Use as coordenas que você capturou próximo ao ponto de saída (exutório). Para Elevação use a camada dem_fill. Use o método padrão Deterministic D8. Nomeie a saída Upslope Area catchment.sdat. Clique em Executar. Clique fechar após concluir o algorítimo.

O resultado deve se parecer com a captura de tela na figura abaixo se você aplicar zoom na extensão da camada.

Para sobrepor o limite da bacia hidrográfica com outros dados, é melhor convertê-lo de raster em vetor (polígono).

3. Para converter a camada raster em vetor, vá ao menu principal  escolha Raster | Converter | Raster para Polígono (Poligonizar).

4. Certifique-se de escolher a entrada e chamar a saída Rur_catchment.shp. Clique Executar. Clique Fechar para voltar à tela principal.

5. Observe o resultado. Verifique também a tabela de atributo (clique com o botão direito do mouse no nome da camada e escolha Abrir tabela de atributo).

Qual valor é atribuído às células dentro do polígono de captação? Qual é o valor para quem está de fora?

No cálculo com o algorítimo UPSLOPE, as células pertencentes à bacia hidrográfica obtêm o valor 100, enquanto as outras células obtêm o valor 0. Durante a conversão para polígonos, pode acontecer que sejam introduzidos erros de geometria. Se você encontrar mais de um recurso com um valor de 100, isso indica um erro de geometria (topologia incorreta), porque o limite do polígono faz um laço. Veja a figura abaixo. Isso pode gerar erros quando usamos o polígono para geoprocessamento.

Claro que estamos interessados apenas na área de captação, então temos que remover o polígono externo.

6. Na tabela de atributos, alterne para o modo de edição usando o botão , então selecione a feição que deseja remover clicando na linha. a seleção será destacada em amarelo no mapa. Clique no botão  para deletar a feição selecionada, desative a edição clicando no botão   novamente, e salve as alterações.

7. Agora remova todas as camadas desnecessárias da lista de camadas para que tenhamos apenas channels, Rur_catchment, dem_fill, e OSM Standard (nessa ordem).

Para visualização, é melhor cortar as camadas até o limite da área de captação.

8. Vamos primeiro cortar a camada vetorial channels para ver apenas os trechos que estão dentro da área de captação. Acesse o menu principal e selecione Vetor | Geoprocessamento | Recortar.

9. Preencha a caixa de diálogo como a figura abaixo para usar a camada de captação como um "cortador de biscoitos" para recortar a camada de entrada de canais a partir dos limites da camada Rur_catchment. Chame a camada Recortada de channels_clip.shp. Clique Executar para executar a ferramenta. Clique em fechar para retornar à tela principal.

Se a receber um erro relacionado à geometria, poderá corrigi-lo. Vá para a caixa de ferramentas de processamento e escolha Geometria do Vetor | Corrigir Geometrias. Outras solução é tentar ferramentas da Caixa de Ferramentas de Processamento que seja menos sensíveis a erros geométricos.

10. Remova a camada channels da lista camadas e arraste a camada channels_clip para o topo da lista.

11. Corte o MDE da mesma maneira que fizemos na  seção 6 (Criar um Subconjunto do MDE). A única diferença é que temos que fornecer um valor sem dados que esteja fora do intervalo de elevações. Você pode usar -9999.