Statistiques par bloc (ArcGIS Spatial Analyst)
Récapitulatif
Cet outil permet de partitionner une entrée en blocs non superposés et de calculer les statistiques relatives aux valeurs dans chaque bloc. La valeur est attribuée à toutes les cellules dans chaque bloc dans la sortie.
Pour en savoir plus sur le fonctionnement de l'outil Statistiques par bloc
Illustration
 |
OutRas = BlockStatistics(InRas1, NbrRectangle(3,3,MAP), "MAXIMUM", "") |
Utilisation
-
Si le raster en entrée est de type virgule flottante, seules les statistiques Moyenne, Maximum, Minimum, Plage, ECARTYPE et Somme sont disponibles. Les statistiques Majorité, Minorité, Médiane et Variété ne sont pas autorisées. Si le raster en entrée est de type entier, tous les types de statistiques sont disponibles.
-
Lorsque vous définissez un voisinage en anneau, circulaire ou en forme de secteur, et compte tenu de la taille de ce voisinage, les cellules qui ne sont pas perpendiculaires à l'axe x ou y ne sont pas prises en compte dans les calculs. Toutefois, les emplacements de ces cellules sont associés à la valeur obtenue suite aux calculs du voisinage, car elles sont comprises dans le rectangle d'emprise minimale (ou le bloc en sortie) de ces types de voisinages circulaires.
-
Si le raster en entrée est un nombre entier, le raster en sortie sera un nombre entier. Il existe une exception dans le cas des types de statistiques Moyenne ou ECARTYPE. Dans ce cas, le raster en sortie est toujours à virgule flottante. Si le type d'entrée est à virgule flottante, la sortie le sera également pour tous les types de statistiques disponibles.
-
Les types de voisinage Irrégulier et Pondération nécessitent qu'un fichier de noyau soit spécifié. Les fichiers de noyau doivent avoir une extension de fichier .txt.
Pour plus d'informations sur la création et l'utilisation d'un fichier noyau, consultez les sections Irrégulier et Pondération dans la rubrique Fonctionnement des statistiques par bloc.
-
Dans le cas d'une statistique de type Médiane, si le nombre de cellules dans le bloc est impair, les valeurs sont classées et la valeur centrale représente la médiane sous la forme d'un entier. Si le nombre de cellules dans le bloc est équivalent, les valeurs sont classées. Les deux valeurs centrales représentent une moyenne, soit l'entier le plus proche.
-
Dans le cas d'une statistique de type Majorité, les cellules dans lesquelles il n'existe aucune valeur (deux valeurs ou plus dans un bloc sont associées, puisqu'elles ont le plus grand nombre de cellules dans la valeur) sont associées à la valeur NoData. Dans le cas d'une statistique de type Minorité, les cellules dans lesquelles il n'existe pas de valeur minoritaire sont également associées à la valeur NoData.
-
Lorsque le type de statistiques est Moyenne, Minorité, Ecart type ou Somme, le type de voisinage peut être associé au paramètre Poids.
Syntaxe
| Paramètre | Explication | Type de données |
in_raster |
Raster faisant l'objet de calculs avec l'outil Statistiques par bloc. | Raster Layer |
neighborhood (Facultatif) |
La classe Voisinage dicte la forme de la zone autour de chaque cellule utilisée pour le calcul des statistiques. Les différents types de voisinage disponibles sont NbrAnnulus, NbrCircle, NbrRectangle, NbrWedge, NbrIrregular, et NbrWeight. Formes des voisinages :
Le voisinage par défaut est un voisinage NbrRectangle carré avec une largeur et une hauteur de 3 cellules. | Neighborhood |
statistics_type (Facultatif) |
Type de statistique à calculer.
Le type de statistiques par défaut est MEAN. | String |
ignore_nodata (Facultatif) |
Indique si les valeurs NoData sont ignorées par le calcul de statistiques.
| Boolean |
Valeur renvoyée
| Nom | Explication | Type de données |
| out_raster |
Raster des statistiques du bloc en sortie. | Raster |
Exemple de code
Cet exemple décrit le calcul de la valeur d'une cellule minimale dans chaque voisinage d'anneau non superposé dans le raster GRID en entrée.
import arcpy
from arcpy import env
from arcpy.sa import *
env.workspace = "C:/sapyexamples/data"
nbr = NbrAnnulus(1, 3, "MAP")
outBlockStat = BlockStatistics("block", nbr, "MINIMUM", "")
outBlockStat.save("C:/sapyexamples/output/blockstat")
Cet exemple décrit le calcul de la valeur d'une cellule minimale dans chaque voisinage d'anneau non superposé dans le raster GRID en entrée.
# Name: BlockStatistics_Ex_02.py
# Description: Calculates statistics for a nonoverlapping
# neighborhood.
# Requirements: Spatial Analyst Extension
# Import system modules
import arcpy
from arcpy import env
from arcpy.sa import *
# Set environment settings
env.workspace = "C:/sapyexamples/data"
# Set local variables
inRaster = "block"
nbr = NbrAnnulus(1, 3, "MAP")
# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")
# Execute BlockStatistics
outBlockStat = BlockStatistics(inRaster, nbr, "MINIMUM", "NODATA")
# Save the output
outBlockStat.save("C:/sapyexamples/output/blockstat")