Outdated version of the documentation. Find the latest one here.

17.15. Recortar e mesclar camadas raster

Nota

Nesta lição nós veremos outro exemplo de preparação de dados espaciais, para continuar usando geoalgoritmos em cenários do mundo real

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

../../../_images/medfordarea.png

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)

  • Elas estão em dois arquivos diferentes. (Os limites da cidade caem em apenas um uma camada raster simples, mas, como foi dito, queremos alguma área extra em torno dele).

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

../../../_images/bbox.png

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

../../../_images/buffer_dialog.png

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

../../../_images/buffer1.png

É 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.

../../../_images/buffer_squared.png

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).

../../../_images/warp1.png

Com essa camada que contém a caixa delimitadora do layer que queremos obter, podemos cortar ambas as camadas raster, usando o algoritmo Cortar Grade com Polígonos.

../../../_images/clip1.png

Uma vez que as camadas foram cortadas, estas podem ser mescladas com o algoritmo Mesclar camadas raster

../../../_images/merge1.png

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.

Nota: Você pode economizar tempo mesclando primeiro e depois cortando, e você vai evitar chamar o algoritmo de recorte duas vezes. No entanto, se há várias camadas de fusão e eles têm um tamanho bastante grande, você vai acabar com uma camada muito grande e mais tarde pode ser difícil de processar. Nesse caso, você pode ter que chamar o algoritmo de recorte várias vezes, o que pode ser demorado, mas não se preocupe, logo veremos que existem algumas ferramentas adicionais para automatizar essa operação. Neste exemplo, só temos duas camadas, e você não deve se preocupar com isso agora.

Com isso, temos o DEM final que queremos.

../../../_images/finaldem.png

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.

../../../_images/warpeddem.png

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

../../../_images/slope.png

E aqui está a camada resultante de declividade.

../../../_images/slopereproj.png

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).

../../../_images/metricconversions.png

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

Aviso

todo: Adicionar imagem

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.