¿Cómo funcionan exactamente las estadísticas zonales?

10

Estoy un poco confundido sobre cómo funcionan las herramientas de estadísticas zonales en QGIS (ya sea la herramienta 'Estadísticas zonales' en las herramientas Ráster de la caja de herramientas de geoalgoritmos QGIS o la herramienta 'Estadísticas ráster de polígonos' en la herramienta Vector SAGA <-> Herramientas ráster )

Básicamente tengo una capa de mapa ráster y algunos polígonos pequeños en una capa separada. Quiero asignar a cada polígono el valor del ráster a continuación. En caso de que un polígono se superponga con varias celdas ráster, solo quiero tener el máximo.

Sin embargo, la salida no tiene ningún sentido y no puedo respaldar el razonamiento de esos resultados (que obviamente no son el máximo de la celda ráster subyacente):

ingrese la descripción de la imagen aquí

¿Alguna idea de lo que estoy haciendo mal o lo que aún no consideré? Los resultados con la API de Python en qgis.analysis.QgsZonalStatistics()(de acuerdo con esta pregunta aquí ) son mucho mejores, pero solo ofrece count, meany sum... pero por alguna razón no maxo min.

GeoEki
fuente
1
Acabo de intentar esto en ArcMap al convertir mi Ráster en Puntos y parece que la herramienta de Estadísticas zonales de QGIS está haciendo lo mismo internamente y usa solo los casos en los que un punto se cruza con un polígono. Sin embargo, parece que no hay una opción de coincidencia apropiada en QGIS (algo comparable a los ofrecidos en la herramienta de unión espacial de ArcMap, por ejemplo, INTERSECT, WITHIN_A_DISTANCE, CLOSEST, ...).
GeoEki
Como pista para resolver esto, calculé el "recuento" para una capa ráster gruesa (aproximadamente un píxel está en el polígono) con las Estadísticas zonales de QGIS y obtuve un valor de 0.45, sorprendentemente no es un número entero. Esto sugiere que QGIS puede de alguna manera ponderar los valores de píxeles. Esta podría ser la razón por la cual el máximo / mínimo no corresponde a lo que espera, posiblemente una implementación deficiente de esta ponderación al calcular el máximo / mínimo.
Jon
Consulte gis.stackexchange.com/a/281753/18189 para obtener una descripción de cómo funciona el algoritmo de estadísticas zonales QGIS.
dbaston

Respuestas:

4

Actualización a abril de 2018 y QGIS 2.18

La versión actual de Zonal Statistics Plugin puede proporcionar:

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

Respuesta original

En resumen, este enfoque con QGIS, por lo que sé, no es posible ver el máximo de píxeles de trama como resultado. Solo el recuento, la media y la suma. He proporcionado algunos recursos que cubren esto.

Aquí está la página de documentación del complemento Zonal Statistics Plugin

Con el complemento de estadísticas zonales icon_zonal_statistics, puede analizar los resultados de una clasificación temática. Le permite calcular varios valores de los píxeles de una capa ráster con la ayuda de una capa de vector poligonal (consulte la figura estadísticas_zonales). Puede calcular la suma, el valor medio y el recuento total de los píxeles que están dentro de un polígono. El complemento genera columnas de salida en la capa vectorial con un prefijo definido por el usuario.

Además, este tutorial QGIS Muestreo de datos ráster utilizando puntos o polígonos proporciona una muestra de datos para lograr el retorno de las estadísticas zonales. He incluido capturas de pantalla que proporcionan como ejemplo de su flujo de trabajo y resultados.

Habilite los complementos de estadísticas zonales. Este es un complemento central, por lo que ya está instalado. Consulte Uso de complementos para saber cómo habilitar el complemento principal

ingrese la descripción de la imagen aquí

ingrese la descripción de la imagen aquí

ingrese la descripción de la imagen aquí

Una vez que finaliza el procesamiento, seleccione la capa de los condados. Use la herramienta Identificar y haga clic en cualquier polígono del condado. Verá tres nuevos atributos agregados a la capa: ZS_count, ZS_mean y ZS_sum. Estos atributos contienen el recuento de píxeles de trama, la media de los valores de píxeles de trama y la suma de los valores de píxeles de trama, respectivamente. Como estamos interesados ​​en la temperatura promedio, el campo ZS_mean será el que se usará.

ingrese la descripción de la imagen aquí

whyzar
fuente
2

A partir de QGIS 3.0 (posiblemente antes, no estoy seguro), las siguientes estadísticas están disponibles:

  • Contar
  • Suma
  • Media
  • Mediana
  • Desviación Estándar
  • Mínimo
  • Máximo
  • Rango
  • Minoría
  • Mayoría (modo)
  • Variedad
  • Diferencia
dbaston
fuente