17.16. Análise hidrológica

Nota

En esta lección vamos a realizar algunos análisis hidrológicos. Este análisis será utilizado en algunas de las siguientes lecciones, como se constituye un muy buen ejemplo de un flujo de trabajo de análisis, y lo utilizaremos para mostrar algunas características avanzadas.

Objetivos: Começando com um MDE, iremos extrair uma rede de canais, delinear bacias hidrográficas e calcular algumas estatísticas.

1. Lo primero es cargar el proyecto con los datos de la lección, que solo contiene un MDT.

   

2. O primeiro módulo a ser executado é Área de captação (em algumas versões do SAGA é chamado de Acumulação de vazão (Top Down)). Você pode utilizar qualquer um dos outros nomeados Área de captação. Eles têm algoritmos diferentes por baixo, mas os resultados são basicamente os mesmos.

3. Selecione o MDE no campo Elevação, e deixe os valores padrão para o resto dos parâmetros.

   

   Alguns algoritmos calculam muitas camadas, mas a camada Área de captação é a única que usaremos. Se quiser, você pode se livrar dos outros.

   El renderizado de la capa no es muy informativa.

   

   Para saber o motivo, você pode dar uma olhada no histograma e verá que os valores não estão distribuídos uniformemente (existem algumas células com valor muito alto, aquelas correspondentes à rede de canais). Use o algoritmo Raster calculator para calcular o logaritmo da área de valor da bacia hidrográfica e você obterá uma camada com muito mais informações

   

4. A área de captação (também conhecida como acúmulo de fluxo) pode ser usada para definir um limite para a iniciação do canal. Isso pode ser feito usando o algoritmo Channel network

   • Initiation grid: use the catchment area layer and not the logarithm one.

   • Limite inicial: 10.000.000

   • Guilabel: “Tipo de Inicialização”: Maior do que

   

   If you increase the Initiation threshold value, you will get a more sparse channel network. If you decrease it, you will get a denser one. With the proposed value, this is what you get.

   

   The image above shows just the resulting vector layer and the DEM, but there should be also a raster layer with the same channel network. That raster will be, in fact, the layer we will be using.

5. Now, we will use the Watersheds basins algorithm to delineate the subbasins corresponding to that channel network, using as outlet points all the junctions in it. Here is how you have to set the corresponding parameters dialog.

   

   Y esto es lo que obtendrá.

   

6. This is a raster result. You can vectorise it using the Vectorising grid classes algorithm.

   

   

Ahora, vamos a tratar de calcular estadísticas sobre los valores de elevación en una de las subcuencas. La idea es tener una capa que simplemente represente la elevación dentro de esa subcuenca y luego pasarla al módulo que calcula estas estadísticas.

1. First, let’s clip the original DEM with the polygon representing a subbasin. We will use the Clip raster 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.

   1. Selecionar um polígono

      

   2. Chame o algoritmo de recorte com os seguintes parâmetros:

      

      The element selected in the input field is, of course, the DEM we want to clip.

      Obteremos algo como isto.

      

2. This layer is ready to be used in the Raster layer statistics algorithm.

   

   As estatísticas resultantes são as seguintes.

   

Vamos a utilizar tanto el procedimiento de cálculo de cuenca y el cálculo de las estadísticas en otras lecciones, para averiguar cómo otros elementos pueden ayudar a automatizar ambos y trabajar más eficazmente.