Disponible con una licencia de Spatial Analyst.
Disponible con una licencia de 3D Analyst.
Resumen
Interpola una superficie de ráster a partir de puntos utilizando una técnica de distancia inversa ponderada (IDW).
Uso
El valor de salida para una celda que utiliza la distancia inversa ponderada (IDW) se limita al rango de valores utilizados para la interpolación. Dado que la IDW es un promedio de distancia ponderada, el promedio no puede ser mayor que la entrada máxima o inferior que la entrada mínima. Por lo tanto, no puede crear crestas o valles si estos extremos aún no se han muestreado (Watson y Philip 1985).
Los mejores resultados de la IDW se obtienen cuando la muestra es lo suficientemente densa respecto de la variación local que intenta simular. Si la muestra de los puntos de entrada es dispersa o despareja, los resultados pueden no representar en forma suficiente la superficie deseada (Watson y Philip 1985).
La influencia de un punto de entrada sobre un valor interpolado es isotrópica. Dado que la influencia de un punto de entrada sobre un valor interpolado se relaciona con la distancia, la IDW no preserva la cresta (Philip y Watson 1982).
El Tamaño de celda de salida se puede definir mediante un valor numérico u obtenerse desde un dataset ráster existente. Si el tamaño de celda no se ha especificado explícitamente como el valor del parámetro, se deriva del entorno Tamaño de celda si se ha especificado. Si no se ha especificado el tamaño de celda del parámetro o el tamaño de celda del entorno, pero se ha configurado el entorno Ráster de alineación, se usa el tamaño de celda del ráster de alineación. Si no se especifica nada, el tamaño de celda se calcula a partir de la anchura o la altura de la extensión dividido por 250, donde la extensión se encuentra en el Sistema de coordenadas de salida especificado en el entorno.
Si el tamaño de celda se especifica usando un valor numérico, la herramienta lo usa directamente para el ráster de salida.
Si el tamaño de celda se especifica mediante un dataset ráster, el parámetro muestra la ruta del dataset ráster en lugar del valor del tamaño de celda. El tamaño de celda de ese dataset ráster se usa directamente en el análisis, siempre que la referencia espacial del dataset sea la misma que la referencia espacial de salida. Si la referencia espacial del dataset es diferente a la referencia espacial de salida, se proyecta en función del Método de proyección de tamaño de celda seleccionado.
Algunos datasets de entrada pueden tener algunos puntos con las mismas coordenadas x,y. Si los valores de los puntos de una ubicación común son los mismos, se consideran duplicados y no afectan a la salida. Si los valores son diferentes, se consideran puntos coincidentes.
Las distintas herramientas de interpolación pueden manejar esta condición de datos de maneras distintas. Por ejemplo, en algunos casos el primer punto coincidente encontrado se utiliza para el cálculo; en otros casos, se utiliza el último punto encontrado. Esto puede causar que algunas ubicaciones del ráster de entrada tengan valores distintos a los que puede esperar. La solución es preparar los datos quitando estos puntos coincidentes. La herramienta Capturar eventos de la caja de herramientas de Estadística espacial es útil para identificar cualquier punto coincidente en los datos.
La opción de barreras se utiliza para especificar la ubicación de entidades lineales conocidas por interrumpir la continuidad de la superficie. Estas entidades no tienen valores z. Acantilados, fallas y terraplenes son ejemplos típicos de barreras. Las barreras limitan el conjunto seleccionado de puntos de muestra de entrada utilizado para interpolar valores z de salida a esas muestras sobre el mismo lado de la barrera como la celda de procesamiento actual. La separación por una barrera se determina mediante un análisis de línea de visión entre cada par de puntos. Esto significa que no se requiere separación topológica para que dos puntos se excluyan mutuamente de la región de influencia del otro. Los puntos de muestra de entrada que yacen exactamente sobre la línea de la barrera se incluyen en el conjunto seleccionado de muestras para ambos lados de la barrera.
Las entidades de barrera se introducen como entidades de polilínea. IDW sólo utiliza las coordenadas x,y para la entidad lineal; por lo tanto, no es necesario proporcionar valores z para los lados izquierdo y derecho de la barrera. Se ignorarán todos los valores z.
El uso de barreras extenderá el tiempo de procesamiento de manera significativa.
Esta herramienta tiene un límite de aproximadamente 45 millones de puntos de entrada. Si la clase de entidad de entrada contiene más de 45 millones de puntos, la herramienta puede fallar al crear un resultado. Puede evitar este límite interpolando el área de estudio en varias partes, asegurándose de que haya alguna superposición en los bordes, y creando luego un mosaico de los resultados para generar un sólo dataset ráster grande. Alternativamente, puede utilizar un dataset de terreno para almacenar y visualizar puntos y superficies comprimidos de mil millones de puntos de medición.
Si tiene la extensión de Geostatistical Analyst, podrá procesar datasets más grandes con la versión de la herramienta IDW que está disponible ahí.
Los datos de la entidad de entrada deben contener por lo menos un campo válido.
Consulte Entornos de análisis y Spatial Analyst para obtener detalles adicionales sobre los entornos de geoprocesamiento que se aplican a esta herramienta.
Referencias:
Philip, G. M. y D. F. Watson. "A Precise Method for Determining Contoured Surfaces." Australian Petroleum Exploration Association Journal 22: 205–212. 1982.
Watson, D. F., and G. M. Philip. "A Refinement of Inverse Distance Weighted Interpolation." Herramienta de geoprocesamiento 2:315–327. 1985.
Sintaxis
Idw(in_point_features, z_field, {cell_size}, {power}, {search_radius}, {in_barrier_polyline_features})| Parámetro | Explicación | Tipo de datos |
in_point_features | Entidades de puntos de entrada que contienen los valores z que se interpolarán en un ráster de superficie. | Feature Layer |
z_field | Campo que contiene un valor de altura o magnitud para cada punto. Puede ser un campo numérico o el campo Forma si las entidades de punto de entrada contienen valores z. | Field |
cell_size (Opcional) | El tamaño de celda del ráster de salida que se va a crear. Este parámetro se puede definir mediante un valor numérico u obtenerse desde un dataset ráster existente. Si el tamaño de celda no se ha especificado explícitamente como valor de parámetro, se usa el valor del tamaño de celda del entorno si se ha especificado; de lo contrario, se utilizan reglas adicionales para calcularlo a partir de las otras entradas. Consulte la sección Uso para obtener más detalles. | Analysis Cell Size |
power (Opcional) | El exponente de distancia. Controla la significancia de los puntos circundantes sobre el valor interpolado. Una potencia mayor tiene como resultado una influencia menor de los puntos lejanos. Puede ser cualquier número real mayor que 0, pero los resultados más razonables se obtendrán con valores de 0,5 a 3. El valor predeterminado es 2. | Double |
search_radius (Opcional) | La clase Radius define cuáles de los puntos de entrada se utilizarán para interpolar el valor para cada celda en el ráster de salida. Hay dos tipos de clases de radios: RadiusVariable y RadiusFixed. Un radio de búsqueda Variable se utiliza para encontrar una cantidad específica de puntos de muestra de entrada para la interpolación. El tipo Fijo utiliza una distancia fija especificada dentro de la cual todos los puntos de entrada se utilizarán para la interpolación. El tipo Variable es la opción predeterminada.
| Radius |
in_barrier_polyline_features (Opcional) | Entidades de polilínea que se utilizarán como corte o límite en la búsqueda de puntos de muestra de entrada. | Feature Layer |
Valor de retorno
| Nombre | Explicación | Tipo de datos |
| out_raster | Ráster de superficie interpolado de salida. Siempre es un ráster de punto flotante. | Raster |
Muestra de código
En este ejemplo se introduce un shapefile de punto y se interpola la superficie de salida como ráster TIFF.
import arcpy
from arcpy import env
from arcpy.sa import *
env.workspace = "C:/sapyexamples/data"
outIDW = Idw("ozone_pts.shp", "ozone", 2000, 2, RadiusVariable(10, 150000))
outIDW.save("C:/sapyexamples/output/idwout.tif")En este ejemplo se introduce un shapefile de punto y se interpola la superficie de salida como ráster de cuadrícula.
# Name: IDW_Ex_02.py
# Description: Interpolate a series of point features onto a rectangular
# raster using Inverse Distance Weighting (IDW).
# 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
inPointFeatures = "ca_ozone_pts.shp"
zField = "ozone"
cellSize = 2000.0
power = 2
searchRadius = RadiusVariable(10, 150000)
# Execute IDW
outIDW = Idw(inPointFeatures, zField, cellSize, power, searchRadius)
# Save the output
outIDW.save("C:/sapyexamples/output/idwout02")Entornos
Información de licenciamiento
- Basic: Requiere Spatial Analyst o 3D Analyst
- Standard: Requiere Spatial Analyst o 3D Analyst
- Advanced: Requiere Spatial Analyst o 3D Analyst