16.2. Análise raster

16.2.1. Calculadora raster

The Raster Calculator in the Raster menu allows you to perform calculations on the basis of existing raster pixel values (see Fig. 16.20). The results are written to a new raster layer in a GDAL-supported format.

../../../_images/raster_calculator1.png

Fig. 16.20 Calculadora raster

The Raster bands list contains all loaded raster layers that can be used. To add a raster to the raster calculator expression field, double click its name in the Fields list. You can then use the operators to construct calculation expressions, or you can just type them into the box.

In the Result layer section, you will need to define an output layer. You can:

  • checkbox Create on-the-fly raster instead of writing layer to disk:

    • If unchecked, the output is stored on the disk as a new plain file. An Output layer path and an Output format are required.

    • If checked, a virtual raster layer, i.e. a raster layer defined by its URI and whose pixels are calculated on-the-fly, is created. It’s not a new file on disk; the virtual layer is still connected to the rasters used in the calculation meaning that deleting or moving these rasters would break it. A Layer name can be provided, otherwise the calculation expression is used as such. Removing the virtual layer from the project deletes it, and it can be made persistent in file using the layer Export ► Save as… contextual menu.

  • Define the Spatial extent of the calculation based on an input raster layer extent, or on custom X,Y coordinates

  • Set the Resolution of the layer using columns and rows number. If the input layer has a different resolution, the values will be resampled with the nearest neighbor algorithm.

  • With the checkbox Add result to project checkbox, the result layer will automatically be added to the legend area and can be visualized. Checked by default for virtual rasters.

The Operators section contains all available operators. To add an operator to the raster calculator expression box, click the appropriate button. Mathematical calculations (+, -, *, … ) and trigonometric functions (sin, cos, tan, … ) are available. Conditional expressions (=, !=, <, >=, … ) return either 0 for false or 1 for true, and therefore can be used with other operators and functions.

Dica

Veja também o algoritmo Calculadora Raster.

16.2.1.1. Exemplos

Converter valores de elevação de metros para pés

Para criar um raster de elevação em pés a partir de um raster em metros, você precisa usar o fator de conversão de metros para pés: 3,28. A expressão é:

Usando uma máscara

Se você deseja mascarar partes de um raster – digamos, por exemplo, porque você só está interessado em elevações acima de 0 metros – você pode usar a seguinte expressão para criar uma máscara e aplicar o resultado a um raster em um passo.

Em outras palavras, para cada célula maior ou igual a 0, a expressão condicional é avaliada como 1, o que mantém o valor original multiplicando-o por 1. Caso contrário, a expressão condicional é avaliada como 0, que define o valor raster como 0. Isso cria a máscara em vôo.

Classify a Raster

Se você deseja classificar um raster – por exemplo, em duas classes de elevação, é possível usar a expressão a seguir para criar um raster com dois valores 1 e 2 em um passo.

("[email protected]" < 50) * 1 + ("[email protected]" >= 50) * 2

Em outras palavras, para cada célula menor que 50, defina seu valor como 1. Para cada célula maior ou igual a 50, defina seu valor como 2.

Or you can use the IF operator.

if ( elevation@1 < 50 , 1 , 2 )

16.2.2. Alinhamento de raster

Essa ferramenta é capaz de receber vários rasters como entrada e alinhá-los perfeitamente, ou seja:

  • reprojetar para o mesmo SRC,

  • reamostrar para o mesmo tamanho de célula e deslocamento na grade,

  • recortar a uma região de interesse,

  • redimensionar valores quando requerido.

Todos os rasters serão salvos em outros arquivos.

Primeiro, abra as ferramentas de Raster ► Alinhar Raster… e pressione no botão symbologyAdd Adicionar novo raster para escolher um raster existente no QGIS. Selecione um arquivo de saída para salvar o raster após o alinhamento, o método de reamostragem e se as ferramentas devem Redimensionar valores de acordo com o tamanho da célula. O método de reamostragem pode ser (veja Fig. 16.21):

  • Vizinho mais próximo

  • Bilinear (2x2 kernel)

  • Cúbica (núcleo 4x4): Aproximação da Convolução Cúbica

  • B-Spline Cúbica (kernel 4x4): Aproximação da B-Spline Cúbica

  • Lanczos (6x6 kernel): Interpolação Lanczos com janela sinc

  • Média: computa a média de todos os pixels contribuintes exceto aqueles SEM DADO

  • Moda: seleciona o valor que aparece com mais frequência de todos os pontos amostrados

  • Máximo, Mínimo, Mediana, Primeiro quartil (Q1) ou Terceiro Quartil (Q3) de todos os pixels que contribuintes exceto aqueles SEM DADO

../../../_images/raster_align_edit.png

Fig. 16.21 Selecione as opções de reamostragem raster

Na caixa de diálogo principal Alinhar raster, você pode ainda symbologyEdit Editar configurações do arquivo ou symbologyRemove Remover um arquivo existente da lista de camadas raster. Você também pode escolher uma ou mais opções adicionais (veja Fig. 16.22):

  • Selecionar a Camada de refêrencia,

  • Transforme para um novo SRC,

  • Configure um diferente Tamanho de célula,

  • Configure um diferente Deslocamento de grade,

  • Recortar à extensão: pode ser definida pelo usuário, ligada a uma camada ou à tela do mapa

  • Tamanho de saída,

  • Adicionar raster alinhado à tela do mapa.

../../../_images/raster_align.png

Fig. 16.22 Alinhamento de raster