Bug principali nelle statistiche di ArcGIS Zonal?

Aggiornamento: il bug è stato corretto nella versione ArcGIS 10.4

Sto usando ArcGIS 10.2.2 per determinare le statistiche zonali per un numero di zone. Se ci sono dei NoData nel valore raster, voglio che i risultati della zona siano "NoData", esattamente come pubblicizzato dalla descrizione degli strumenti. La descrizione di questo strumento indica:

DATI - All'interno di una determinata zona, solo le celle che hanno un valore nel valore raster di input verranno utilizzate per determinare il valore di output per quella zona. Le celle NoData nel raster Value verranno ignorate nel calcolo statistico.

NODATA - All'interno di una determinata zona, se nel raster Value sono presenti celle NoData, si ritiene che non vi siano informazioni sufficienti per eseguire calcoli statistici per tutte le celle in quella zona; pertanto, l'intera zona riceverà il valore NoData sul raster di output.

Dai un'occhiata al mio setup in questa immagine:

Sto usando l'opzione NODATA con un valore raster che ha un pixel NoData e quindi mi aspetto che il valore della zona risultante (zona 61154) sia 'NoData'. Invece, ottengo un valore di 12,74 (arrotondato a 13 nell'immagine), il che mi confonde su due livelli: in primo luogo, mi aspettavo "NoData", e in secondo luogo, il valore risultante di 12,74 è matematicamente impossibile, perché la media non può essere maggiore rispetto al valore massimo nel valore raster, che in questo caso è 10.

Se sto usando l'opzione DATA, ottengo un valore di circa 9.1, il che ha senso. Lo abbiamo testato su diversi set di dati, computer e versioni ArcGIS.

Cosa mi sto perdendo qui?

Modifica / Commento aggiuntivo: ho appena notato che l'attributo 'Count' è sbagliato anche per quella particolare zona. Esistono effettivamente 421 celle in quella zona, ma lo strumento ha contato solo 297. Calcolando 421 meno 297 si ottiene 124 - stranamente, questa è la "posizione" in cui si trova il pixel NoData, se si contano i pixel dalla parte superiore sinistra a quella inferiore proprio nella zona. Lo strumento potrebbe sbagliare il conteggio delle celle (troppo basso), il che potrebbe spiegare l'aumento della media.

Modifica: ecco un link ai dati che sto usando.

Modifica: Dan Patterson e io abbiamo fatto qualche altro debug qui sul forum ESRI.

Esiste un bug che sembra corrispondere a quello che stai riscontrando - è registrato come BUG-000084883 - L'opzione 'Ignora NoData nei calcoli' in Zonal Statistics come strumento Table {e Zonal Statistics} non viene onorata quando spuntata, producendo risultati errati.

Si verifica con 10.3 e 10.2.2 ma non 10.1. Hai provato lo strumento con questa versione?

È un bug. Qualcosa di terribilmente sbagliato nel conteggio delle cellule.

La media corretta (9.0452380952381) volte il numero corretto di celle non vuote (420) diviso per 297 (ovvero un conteggio di celle riportato dallo strumento ) risulta in 12.7912457912458. Questa è una media errata riportata dallo strumento.

Risultati del mio test sulle griglie delle dimensioni dei giocattoli:

Simile a un'altra risposta , sposta i dati raster in array mascherati NumPy per calcolare le tue statistiche. Supponendo due raster sovrapposti con la stessa forma, questo è semplice:

import numpy as np
zones = arcpy.RasterToNumPyArray("zones")
value = np.ma.masked_equal(arcpy.RasterToNumPyArray("value"),
                           arcpy.Raster("value").noDataValue)
print("Zone\tCount\tNoData\tMean")
for z in np.unique(zones):
    sel = (zones == z)
    print z, sel.sum(), value.mask[sel].sum(), value[sel].mean()

Spettacoli:

Zone    Count   NoData  Mean
61131   53   0   8.92452830189
61154   421   1   9.04523809524
61207   1   0   8.0
61317   35   0   7.2
61644   644   0   7.90838509317
61677   12   0   7.41666666667
61789   7   0   9.0
61871   193   0   7.98445595855
187472   349   0   8.5787965616