17.15. Recortar e mesclar camadas raster¶

Para esta lição, vamos calcular uma camada de inclinação para uma área em torno de uma área da cidade, que é dado em uma camada vetorial com um único polígono. O DEM base é dividido em duas camadas raster que, em conjunto, cobrem uma área muito maior do que o entorno da cidade que queremos trabalhar. Se você abrir o projeto correspondente a esta lição, você vai ver algo como isto

Essas camadas têm dois problemas:

• Eles cobrem uma área que é muito grande para o que queremos (nós estamos interessados ​​em uma região menor ao redor do centro da cidade)

• They are in two different files (the city limits fall into just one single raster layer, but, as it’s been said, we want some extra area around it).

Ambos são facilmente solucionáveis com os geoalgoritmos adequados.

Primeiro. Nós criamos um retângulo definindo a área que nós queremos. Para fazer isto, nós criamos uma camada contendo a caixa delimitadora da camada com os limites da área da cidade, e então lhe aplicamos um buffer, de modo a ter uma camada raster que cobre um pouco mais do que o estritamente necessário.

Para calcular a caixa delimitadora, nós podemos usar o algoritmo Polígono a partir da extensão da camada

Para aplicar o buffer, nós usamos o Buffer de distância fixa, com os seguintes valores de parâmetros.

Aqui está a caixa delimitadora obtida usando os parâmetros acima

É uma caixa arredondada, mas podemos facilmente chegar a caixa equivalente com ângulos retos, executando o algoritmo Polígono a partir da extensão da camada nele. Poderíamos ter aplicado um buffer dos limites da cidade em primeiro lugar, e depois calcular a extensão retângulo, economizando uma etapa.

Você vai notar que os rasters tem uma projeção diferente do vetor. Devemos, portanto, reprojetá-los antes de prosseguir, usando a ferramenta Deformar (reprojeção).

With this layer that contains the bounding box of the raster layer that we want to obtain, we can crop both of the raster layers, using the Clip grid with polygon algorithm.

Uma vez que as camadas foram cortadas, eles podem ser mescladas usando o algoritmo GDAL Merge.

O tamanho da célula é necessário para a camada mesclada. Vamos usar o mesmo da entrada. Você não precisa saber com antecedência antes de chamar o algoritmo. Basta clicar no botão direito - tamanho do campo de texto e você terá uma janela para introduzir pequenas fórmulas matemáticas, e uma lista de valores freqüentemente usados, entre eles os tamanhos de células e coordenadas delimitadoras, de todas as camadas disponíveis.

Com isso, temos o DEM final que queremos.

Agora é hora de calcular a camada de declividade.

A camada de inclinação pode ser calculado com o algoritmo Slope, Aspect, Curvature, mas o DEM obtido na última etapa não é adequado como entrada, uma vez que valores de elevação estão em metros, mas o tamanho da célula não é expresso em metros (a camada usa um CRS com coordenadas geográficas). É necessária uma reprojeção. Para reprojetar uma camada raster, o algoritmo *Deformar (reprojetar) * pode ser usado novamente. Nós reprojetamos em um CRS com medidores como unidades (por exemplo, 3857), para que possamos, então, calcular correctamente a declividade, tanto com SAGA ou GDAL.

Aqui está o DEM reprojetado.

Com o novo DEM, a declividade pode ser agora calculada.

E aqui está a camada resultante de declividade.

A declividade produzida pelo algoritmo Slope, Aspect, Curvature é expresso em radianos, mas graus são uma unidade mais prático e comum. O algoritmo Conversões métricas nos ajudará a fazer a conversão (mas caso você não sabia que existia o algoritmo, você pode usar a calculadora raster que nós já usamos).

Reprojetando a camada de declividade convertida novamente com o Reprojetar camada raster, temos a camada final que queríamos.

Os processos de reprojeção podem fazer com que a camada final possa conter dados fora da caixa delimitadora que foi calculada em um dos primeiros passos. Isso pode ser resolvido por recorte-o novamente, como fizemos para obter o DEM base.