Como as Estatísticas Zonais funcionam exatamente?

10

Estou um pouco confuso sobre como as ferramentas de estatística zonal funcionam no QGIS (a ferramenta 'Zonal Statistics' sob as ferramentas Raster da caixa de ferramentas QGIS geoalgorithms ou a ferramenta 'Raster statistics from polygons' na ferramenta SAGA Vector <-> Raster tools )

Basicamente, tenho uma camada de mapa raster e alguns polígonos pequenos em uma camada separada. Quero atribuir a cada polígono o valor da varredura abaixo. Caso um polígono se sobreponha a várias células raster, só quero ter o máximo.

No entanto, a saída não faz sentido e não posso ficar por trás do raciocínio desses resultados (que obviamente não são o máximo da célula de varredura subjacente):

insira a descrição da imagem aqui

Alguma idéia do que estou fazendo de errado ou do que ainda não considerei? Os resultados com a API python em qgis.analysis.QgsZonalStatistics()( de acordo com esta pergunta aqui ) são muito melhores, mas faz única oferta count, meane sum... mas por alguma razão não maxou min.

GeoEki
fonte
1
Apenas tentei isso no ArcMap convertendo meu Raster em pontos e parece que a ferramenta Estatísticas Zonais do QGIS está fazendo o mesmo internamente e usa apenas os casos em que um ponto cruza um polígono. No entanto, parece não haver uma opção de correspondência apropriada no QGIS (algo comparável aos oferecidos na ferramenta Spatial Join do ArcMap, por exemplo, INTERSECT, WITHIN_A_DISTANCE, CLOSEST, ...).
GeoEki
Como uma pista para descobrir isso, calculei a "contagem" de uma camada rasterizada grossa (aproximadamente um pixel está no polígono) com as Estatísticas Zonais do QGIS e obtive um valor de 0,45 - surpreendentemente, não um número inteiro. Isso sugere que o QGIS pode, de alguma forma, estar ponderando os valores de pixel. Talvez seja por isso que max / min não corresponda ao que você espera - possivelmente uma implementação ruim dessa ponderação ao calcular max / min.
22418 Jon
Consulte gis.stackexchange.com/a/281753/18189 para obter uma descrição de como o algoritmo de estatísticas zonais do QGIS funciona.
dbaston

Respostas:

4

Atualização em abril de 2018 e QGIS 2.18

A versão atual do Zonal Statistics Plugin pode fornecer:

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

Resposta original

Em suma, essa abordagem com o QGIS, até onde sei, não é possível visualizar o máximo de pixels rasterizados como resultado. Somente a contagem, média e soma. Eu forneci alguns recursos que cobrem isso.

Aqui está a página de documento do plug-in Zonal Statistics Plugin

Com o plugin de estatísticas zonais icon_zonal_statistics, você pode analisar os resultados de uma classificação temática. Permite calcular vários valores dos pixels de uma camada raster com a ajuda de uma camada vetorial poligonal (consulte figure_zonal_statistics). Você pode calcular a soma, o valor médio e a contagem total de pixels que estão dentro de um polígono. O plug-in gera colunas de saída na camada vetorial com um prefixo definido pelo usuário.

Além disso, este Tutorial de amostragem de dados raster do QGIS usando pontos ou polígonos fornece uma amostra de dados para obter retorno estatístico zonal. Eu incluí a captura de tela que eles fornecem como um exemplo do fluxo de trabalho e dos resultados.

Habilite os plug-ins de estatísticas de zonas. Este é um plug-in principal, por isso já está instalado. Consulte Usando plug-ins para saber como ativar o plug-in principal

insira a descrição da imagem aqui

insira a descrição da imagem aqui

insira a descrição da imagem aqui

Quando o processamento terminar, selecione a camada de condados. Use a ferramenta Identificar e clique em qualquer polígono do município. Você verá três novos atributos adicionados à camada: ZS_count, ZS_mean e ZS_sum. Esses atributos contêm a contagem de pixels de varredura, a média dos valores de pixels de varredura e a soma dos valores de pixels de varredura, respectivamente. Como estamos interessados ​​em temperatura média, será o campo ZS_mean a ser usado.

insira a descrição da imagem aqui

whyzar
fonte
2

A partir do QGIS 3.0 (possivelmente mais cedo, não tenho certeza), as seguintes estatísticas estão disponíveis:

  • Contagem
  • Soma
  • Significar
  • Mediana
  • Desvio padrão
  • Mínimo
  • Máximo
  • Alcance
  • Minoria
  • Maioria (modo)
  • Variedade
  • Variação
dbaston
fonte