Estadísticas de bloque (ArcGIS 10.1 Spatial Analyst)
Resumen
Divide la entrada en bloques no superpuestos y calcula las estadísticas de los valores dentro de cada bloque. El valor se asigna a todas las celdas de cada bloque en la salida.
Ilustración
 |
OutRas = BlockStatistics(InRas1, NbrRectangle(3,3,MAP), "MAXIMUM", "") |
Uso
-
Si el ráster de entrada es de tipo de punto flotante, solo se encontrarán disponibles las estadísticas Media, Máximo, Mínimo, Rango, DTP y Suma, mientras que las estadísticas Mayoría, Minoría, Mediana y Diversidad no están permitidas. Si el ráster de entrada es de tipo entero, todos los tipos de estadísticas se encontrarán disponibles.
-
Cuando se especifica una vecindad circular, con forma de anillo o con forma de porción de círculo, según el tamaño de la vecindad, es probable que para el cálculo no se tengan en cuenta las celdas que no son perpendiculares al eje x o y. Sin embargo, estas ubicaciones de celda recibirán el valor resultante de los cálculos de la vecindad porque caen dentro del rectángulo mínimo de delimitación (o el bloque de salida) de estos tipos de vecindad circular.
-
Si el ráster de entrada es un entero, el ráster de salida será un entero. Existe una excepción para los tipos de estadística de valor medio o STD, para los cuales el ráster de salida siempre será un punto flotante. Si el tipo de entrada es flotante, la salida será de tipo flotante para todos los tipos de estadísticas disponibles.
-
Los tipos de Vecindad Irregular y Peso requieren que se especifique un Archivo kernel. Los archivos kernel deben tener una extensión de archivo .txt.
Consulte las secciones Irregular y Peso de Cómo funciona Estadísticas de bloques para obtener información sobre la creación y el uso de archivos kernel.
-
Para el tipo de estadística mediana, si la cantidad de celdas en el bloque es impar, se clasifican los valores y se toma el valor medio como la mediana, y es un entero. Si la cantidad de celdas en el bloque es par, se clasifican los valores y se promedian los dos valores medios al entero más próximo.
-
Para el tipo de estadística de mayoría, a las celdas que no poseen un valor de mayoría simple (es decir, dos o más valores dentro de un bloque están vinculados por tener la cantidad máxima de celdas con ese valor) se les asigna NoData. Para el tipo de estadística de minoría, a las celdas que no poseen un valor de minoría simple también se les asigna NoData.
-
Para el tipo de estadística de valor medio, minoría, desviación estándar o suma, el tipo Vecindad se puede establecer en Peso.
Sintaxis
| Parámetro | Explicación | Tipo de datos |
in_raster |
El ráster donde se realizan los cálculos de estadísticas en bloque. | Raster Layer |
neighborhood (Opcional) |
La clase de Vecindad indica la forma del área alrededor de cada celda utilizada para calcular la estadística. Los diferentes tipos de vecindad disponible son NbrAnnulus, NbrCircle, NbrRectangle, NbrWedge, NbrIrregular, y NbrWeight. A continuación, se presentan las formas de vecindad:
La vecindad predeterminada es un NbrRectangle cuadrado de ancho y altura de 3 celdas. | Neighborhood |
statistics_type (Opcional) |
El tipo de estadística que se calculará.
El tipo de estadísticas predeterminado es MEAN. | String |
ignore_nodata (Opcional) |
Indica si los valores NoData serán ignorados por el cálculo de estadísticas.
| Boolean |
Valor de retorno
| Nombre | Explicación | Tipo de datos |
| out_raster |
El ráster de estadísticas de bloques de salida. | Raster |
Ejemplo de código
En este ejemplo se calcula el valor de celda mínimo dentro de cada vecindad de anillo no superpuesto en el ráster GRID de salida.
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")
En este ejemplo se calcula el valor de celda mínimo dentro de cada vecindad de anillo no superpuesto en el ráster GRID de salida.
# 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")