In questo articolo giocherò un po’ con i raster e il calcolatore raster per dimostrare quanto è semplice manipolare parti di un raster anche con leggi di variazioni non usuali.
Supponiamo di avere un raster e di voler riclassificare porzioni dello stesso seguendo una legge di variazione non lineare come quella descritta sotto:

i valori di Rmin e Rmax devono essere relativi alle porzioni di raster delimitati da poligoni (vedi figura sotto):

Per calcolare Rmin e Rmax del raster in ogni poligono, basti utilizzare il geoalgoritmo: statistiche zonali presenti nel processing di QGIS:

il risultato sarà un nuovo layer poligonale con geometria di partenza e con altri campi relativi alle statistiche:

creiamo i raster r_min (Rmin) r_max (Rmax) con Rasterizza dal progessing di QGIS, prendendo i valori dal layer statistiche zonali:

creiamo, ora, altri tre raster con il calcolatore raster di QGIS:

r_delta:

r_min+r_delta:

r_min+2r_delta:

a questo punto abbiamo tutti i raster per applicare la legge di variazione non lineare, l’equazione che useremo nel calcolatore raster è la seguente (che rappresenta la legge descritta all’inizio):

( ( “f_ndvi@1” > “r_min@1” ) AND ( “f_ndvi@1” <= “r_min_+_r_delta@1” ) ) * 1 +
( ( “f_ndvi@1” > “r_min_+_r_delta@1”) AND ( “f_ndvi@1” <= “r_min_+_r_2delta@1” ) ) * 2 +
( ( “f_ndvi@1” > “r_min_+_r_2delta@1”) AND ( “f_ndvi@1” <= “r_max@1” ) ) * 3

l’equazione è caratterizzata da tre addendi – (A)*1+(B)*2+(C)*3 – in questa equazione solo un addendo per volta sarà vero (1) , per esempio A, e gli altri due B e C saranno falsi (0) e cosi via per gli altri casi.

ecco di seguito il raster risultante:

Tutto il processo si potrebbe sintetizzare in uno modello.

ConQGIStando!!!