17.30. Interpolazione e curve di livello

Modulo contribuito da Paolo Cavallini - Faunalia

Nota

Questo capitolo mostra come utilizzare diversi backend per calcolare interpolazioni differenti.

17.30.1. Interpolazione

Il progetto mostra un gradiente di precipitazioni, da sud a nord. Usiamo diversi metodi di interpolazione, tutti basati sul vettore points.shp, parametro RAIN:

Avvertimento

Imposta la dimensione della cella a 500 per tutte le analisi.

  • GRASS ► v.surf.rst

  • SAGA ► Multilevel B-Spline Interpolation

  • SAGA ► Inverse Distance Weighted [Distanza inversa a una potenza; Potenza: 4; Raggio di ricerca: Globale; Raggio di ricerca: tutti i punti]

  • GDAL ► Grid (Inverse Distance to a power) [Power:4]

  • GDAL ► Grid (Moving average) [Radius1&2: 50000]

Quindi misura la variazione tra i metodi e correlala alla distanza dai punti:

  • GRASS ► r.series [Unselect Propagate NULLs, Aggregate operation: stddev]

  • GRASS ► v.to.rast.value on points.shp

  • GDAL ► Proximity

  • GRASS ► r.covar per mostrare la matrice di correlazione; verifica la significatività della correlazione, ad esempio con http://vassarstats.net/rsig.html.

Pertanto, le aree lontane dai punti avranno un’interpolazione meno accurata.

17.30.2. Curve di livello

Vari metodi per disegnare le curve di livello [sempre passo=10] sul raster stddev:

  • GRASS ► r.contour.step

  • GDAL ► Contour

  • SAGA ►Contour lines from grid [NB: in alcune vecchie versioni di SAGA, l’output shp non è valido, bug noto].