Come funzionano esattamente le statistiche zonali?


10

Sono un po 'confuso su come funzionano gli strumenti di statistica zonale in QGIS (o lo strumento "Statistiche zonali" sotto gli strumenti Raster dallo strumento QGIS geoalgorithms o lo strumento "Statistiche raster da poligoni" sotto gli strumenti SAGA Vector <-> Raster ).

Fondamentalmente ho un livello di mappa raster e alcuni piccoli poligoni in un livello separato. Voglio assegnare a ciascun poligono il valore del raster riportato di seguito. Nel caso in cui un poligono si sovrapponga a più celle raster, voglio solo avere il massimo.

Tuttavia, l'output non ha alcun senso e non riesco a superare il ragionamento di quei risultati (che ovviamente non sono il massimo della cella raster sottoposta a denuncia):

inserisci qui la descrizione dell'immagine

Qualche idea su cosa sto facendo di sbagliato o cosa non ho ancora preso in considerazione? I risultati con l'API python in qgis.analysis.QgsZonalStatistics()( secondo questa domanda qui ) sono molto migliori, ma offre solo count, meane sum... ma per qualche motivo no maxo min.


1
Ho appena provato questo in ArcMap convertendo il mio Raster in Punti e sembra che lo strumento QGIS Zonal Stats stia facendo lo stesso internamente e usi solo i casi in cui un punto interseca un poligono. Tuttavia, sembra non esserci un'opzione di corrispondenza appropriata in QGIS (qualcosa di paragonabile a quelli offerti nello strumento di unione spaziale di ArcMap, ad esempio INTERSECT, WITHIN_A_DISTANCE, CLOSEST, ...).
GeoEki,

Come indizio per capirlo, ho calcolato il "conteggio" per un livello raster approssimativo (all'incirca un pixel è nel poligono) con le statistiche zonali di QGIS e ho ottenuto un valore di 0,45 - sorprendentemente non un numero intero. Ciò suggerisce che QGIS potrebbe in qualche modo ponderare i valori dei pixel. Questo potrebbe essere il motivo per cui il max / min non corrisponde a quello che ti aspetti, probabilmente una cattiva implementazione di questa ponderazione durante il calcolo di max / min.
Jon

Vedi gis.stackexchange.com/a/281753/18189 per una descrizione di come funziona l'algoritmo delle statistiche zonali QGIS.
dbaston,

Risposte:


4

Aggiornamento come aprile 2018 e QGIS 2.18

La versione corrente del plugin Zonal Statistics può fornire:

Count: to count the number of pixels
Sum: to sum the pixel values
Mean: to get the mean of pixel values
Median: to get the median of pixel values
StDev: to get the standard deviation of pixel values
Min: to get the minimum of pixel values
Max: to get the maximum of pixel values
Range: to get the range (max - min) of pixel values
Minority: to get the less represented pixel value
Majority: to get the most represented pixel value
Variety: to count the number of distinct pixel values

Risposta originale

In breve, questo approccio con QGIS, per quanto ne so, non è possibile visualizzare il massimo dei pixel raster di conseguenza. Solo il conteggio, la media e la somma. Ho fornito alcune risorse che trattano questo.

Ecco la pagina doc per il plugin Zonal Statistics Plugin

Con il plug-in delle statistiche zonali icon_zonal_statistics, puoi analizzare i risultati di una classificazione tematica. Consente di calcolare diversi valori dei pixel di un livello raster con l'aiuto di un livello vettoriale poligonale (vedere figure_zonal_statistics). È possibile calcolare la somma, il valore medio e il conteggio totale dei pixel all'interno di un poligono. Il plug-in genera colonne di output nel livello vettoriale con un prefisso definito dall'utente.

Inoltre, questo tutorial QGIS che campiona i dati raster usando punti o poligoni fornisce un campione di dati per ottenere il ritorno delle statistiche zonali. Ho incluso screenshot che forniscono come esempio del flusso di lavoro e dei risultati.

Abilita i plugin delle statistiche zonali. Questo è un plug-in principale, quindi è già installato. Consulta Utilizzo dei plug-in per sapere come abilitare il plug-in principale

inserisci qui la descrizione dell'immagine

inserisci qui la descrizione dell'immagine

inserisci qui la descrizione dell'immagine

Al termine dell'elaborazione, seleziona il livello delle contee. Utilizza lo strumento Identifica e fai clic su qualsiasi poligono della contea. Vedrai tre nuovi attributi aggiunti al layer: ZS_count, ZS_mean e ZS_sum. Questi attributi contengono rispettivamente il conteggio dei pixel raster, la media dei valori dei pixel raster e la somma dei valori dei pixel raster. Poiché siamo interessati alla temperatura media, il campo ZS_mean sarà quello da utilizzare.

inserisci qui la descrizione dell'immagine


2

A partire da QGIS 3.0 (possibilmente in precedenza, non sono sicuro), sono disponibili le seguenti statistiche:

  • Contare
  • Somma
  • Significare
  • Mediano
  • Deviazione standard
  • Minimo
  • Massimo
  • Gamma
  • Minoranza
  • Maggioranza (modalità)
  • Varietà
  • Varianza
Utilizzando il nostro sito, riconosci di aver letto e compreso le nostre Informativa sui cookie e Informativa sulla privacy.
Licensed under cc by-sa 3.0 with attribution required.