15.2. Analisi raster

15.2.1. Calcolatore raster

Il Calcolatore Raster nel menu Raster ti permette di eseguire calcoli sulla base dei valori dei pixel raster esistenti (vedi Fig. 15.17). I risultati vengono scritti in un nuovo layer raster in un formato supportato da GDAL.

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

Fig. 15.17 Calcolatore raster

La lista Bande raster elenca i raster caricati che possono essere utilizzati. Per aggiungere un raster nella finestra Espressione del calcolatore di raster, fai doppio click sul suo nome nella finestra Bande raster. Puoi usare gli operatori per costruire un’espressione oppure puoi scriverla direttamente nella finestra delle Espressione del calcolatore di raster.

Nella sezione Risultato del layer devi definire il nome del raster in output. Hai diverse opzioni per scegliere l’estensione dell’area di calcolo: sulla base del raster in input, su coordinate X e Y oppure sulla base del numero di righe e colonne, per impostare la risoluzione finale del raster in output. Se il raster in input ha una risoluzione diversa, i valori verranno ricampionati con l’algoritmo del vicino più prossimo.

La sezione Operatori contiene tutti gli operatori disponibili. Per aggiungere un operatore alla casella di espressione del calcolatore raster, clicca sul pulsante appropriato. Sono disponibili calcoli matematici (+, -, *, … ) e funzioni trigonometriche (sin, cos, tan, … ). Le espressioni condizionali (=, !=, <, >=, … ) restituiscono 0 per falso o 1 per vero, e quindi possono essere usate con altri operatori e funzioni.

Selezionando la casella di controllo checkbox Aggiungi al progetto il raster finale verrà aggiunto automaticamente alla legenda e lo potrai visualizzare sulla mappa.

Suggerimento

Vedi anche l’algoritmo Calcolatore Raster.

15.2.1.1. Esempi

Convertire unità di misura dell’elevazione da metri a piedi

Per creare un raster con altimetria in piedi a partire da un raster con altimetria in metri devi utilizzare il fattore di conversione di 3,28 piedi per metro. L’espressione è:

"elevation@1" * 3.28

Utilizzare una maschera

Se vuoi usare una maschera su parti di un raster - perché, per esempio, sei interessato solo ai valori di altezza superiori allo 0 - puoi usare l’espressione seguente per creare in un solo passaggio una maschera e applicare il risultato al raster:

("elevation@1" >= 0) * "elevation@1"

In altre parole, per ogni cella maggiore o uguale a 0 imposta il valore uguale a 1, in questo modo si mantiene il valore originario moltiplicandolo per 1. Altrimenti imposta 0, il che imposta il valore del raster a 0. In questo modo creerai la maschera al volo.

Se vuoi classificare un raster, ad esempio, in due classi di elevazione, puoi utilizzare la seguente espressione per creare un raster con due valori 1 e 2 in un solo passo:

("elevation@1" < 50) * 1 + ("elevation@1" >= 50) * 2

In altre parole, per ogni cella minore di 50 imposta il valore su 1. Per ogni cella maggiore o uguale a 50 imposta il valore su 2.

15.2.2. Allineamento Raster

Questo strumento è in grado di prendere diversi raster come input e di allinearli perfettamente, il che significa:

  • riproiettare allo stesso SR,

  • ricampionare alla stessa dimensione di cella e offset nella griglia,

  • ritagliare a una regione di interesse,

  • ridimensionare i valori quando necessario.

Tutti i raster saranno salvati in altri file.

Per prima cosa, apri gli strumenti da Raster ► Allinea Raster… e clicca sul pulsante signPlus Aggiungi raster per scegliere un raster esistente in QGIS. Seleziona un file di output per salvare il raster dopo l’allineamento, il metodo di ricampionamento e se gli strumenti hanno bisogno di Riscala i valori in base alla dimensione della cella. Il metodo di ricampionamento può essere (vedi Fig. 15.18):

  • Vicino più prossimo (Nearest Neighbor)

  • Bilineare (2x2 kernel)

  • Cubica (4x4 kernel): Approssimazione Convoluzione Cubica

  • Cubica B-Spline (4x4 kernel): Approssimazione B-Spline Cubica

  • Lanczos (6x6 kernel): interpolazione limitata a bande di Lanczos

  • Media: calcola la media di tutti i pixel che non contengono NODATA

  • Moda: seleziona il valore che appare più spesso tra tutti i punti campionati

  • Massimo, Minimo, Mediana, Primo Quartile (Q1) o Terzo Quartile (Q3) di tutti i pixel diversi da NODATA

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

Fig. 15.18 Selezione Opzioni di ricampionamento raster

Nella finestra di dialogo principale Allinea raster, puoi scegliere symbologyEdit Modifica impostazioni file o signMinus Rimuovi un file esistente dalla lista dei layer raster. Puoi anche scegliere una o più opzioni (vedi Fig. 15.19):

  • Seleziona Layer di riferimento,

  • Trasforma in un nuovo SR,

  • Settaggio ad una diversa Dimensione cella,

  • Settaggio ad un differente Offset reticolo,

  • Ritaglia all’estensione: può essere definito dall’utente, ritagliato a un layer o all’estensione della mappa

  • Dimensioni in uscita,

  • Aggiungi i raster allineati alla mappa.

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

Fig. 15.19 Allineamento Raster