Un bug majeur dans les statistiques ArcGIS Zonal?

Mise à jour: le bogue a été corrigé dans la version d'ArcGIS 10.4

J'utilise ArcGIS 10.2.2 pour déterminer les statistiques zonales pour un certain nombre de zones. S'il y a des NoData dans le raster de valeurs, je veux que les résultats de la zone soient "NoData", exactement comme annoncé par la description des outils. Cette description de l'outil indique:

DONNÉES - Dans une zone particulière, seules les cellules qui ont une valeur dans le raster de valeur d'entrée seront utilisées pour déterminer la valeur de sortie pour cette zone. Les cellules NoData du raster Value seront ignorées dans le calcul des statistiques.

NODATA - Dans une zone particulière, s'il existe des cellules NoData dans le raster de valeur, il est considéré qu'il n'y a pas suffisamment d'informations pour effectuer des calculs statistiques pour toutes les cellules de cette zone; par conséquent, la zone entière recevra la valeur NoData sur le raster en sortie.

Veuillez jeter un œil à ma configuration sur cette image:

J'utilise l'option NODATA avec un raster de valeurs qui a un pixel NoData, et je m'attends donc à ce que la valeur de zone résultante (zone 61154) soit «NoData». Au lieu de cela, j'obtiens une valeur de 12,74 (arrondie à 13 dans l'image), ce qui me confond à deux niveaux: Premièrement, je m'attendais à 'NoData', et deuxièmement, la valeur résultante de 12,74 est mathématiquement impossible, car la moyenne ne peut pas être plus grande que la valeur maximale dans le raster de valeurs, qui est de 10 dans ce cas.

Si j'utilise l'option DATA, j'obtiens une valeur d'environ 9,1, ce qui est logique. Nous l'avons testé sur différents ensembles de données, ordinateurs et versions d'ArcGIS.

Qu'est-ce que j'oublie ici?

Edit / Commentaire supplémentaire: Je viens de remarquer que l'attribut 'Count' est également incorrect pour cette zone particulière. Il y a en effet 421 cellules dans cette zone, mais l'outil n'a compté que 297. Calculer 421 moins 297 donne 124 - curieusement, c'est la "position" où se trouve le pixel NoData, si l'on compte les pixels du coin supérieur gauche au coin inférieur juste dans la zone. L'outil peut obtenir un nombre de cellules erroné (trop bas), ce qui pourrait expliquer l'augmentation de la moyenne.

Modifier: Voici un lien vers les données que j'utilise.

Edit: Dan Patterson et moi avons fait un débogage supplémentaire ici au forum ESRI.

Il y a un bogue qui semble correspondre à ce que vous rencontrez - il est enregistré sous le nom BUG-000084883 - L'option 'Ignorer NoData dans les calculs' dans l'outil Statistiques zonales en tant qu'outil de table {et l'outil Statistiques zonales} n'est pas respectée lorsqu'elle est cochée, ce qui produit résultats incorrects.

Cela se produit avec 10.3 et 10.2.2 mais pas 10.1. Avez-vous essayé l'outil avec cette version?

C'est un bug. Quelque chose de terrible avec le nombre de cellules.

La moyenne correcte (9,0452380952381) fois le nombre correct de cellules non vides (420) divisé par 297 (c'est-à-dire le nombre de cellules rapporté par l' outil ) donne 12,7912457912458. C'est une mauvaise moyenne rapportée par l'outil.

Résultats de mon propre test de grilles de taille de jouet:

Semblable à une autre réponse , déplacez les données raster dans des tableaux masqués NumPy pour calculer vos statistiques. En supposant deux rasters superposés de même forme, c'est simple:

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()

Spectacles:

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