Vecino natural (3D Analyst)
Resumen
Interpola una superficie de ráster a partir de puntos utilizando una técnica de vecinos naturales.
Uso
-
Si el centro de celda de las celdas del perímetro del ráster de salida cae fuera de la envoltura convexa (definida por los puntos de entrada), entonces a esas celdas se le asignarán valores NoData. Si un punto de entrada cae dentro de una de estas celdas del perímetro y el centro de la celda cae fuera de la envoltura convexa, aún se le asignará un valor de NoData a la celda.
-
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.
-
Esta herramienta tiene un límite de aproximadamente 15 millones de puntos de entrada. Si la clase de entidad de entrada contiene una cantidad extremadamente grande de puntos (alrededor de 15 millones o más), la herramienta puede fallar al crear un resultado.
Puede evitar este límite procesando el área de estudio en varias secciones y creando un mosaico de los resultados para generar un sólo dataset ráster grande. Asegúrese de que haya alguna superposición entre las secciones. Alternativamente, puede utilizar un Dataset de terreno para almacenar y visualizar puntos y superficies comprimidos de mil millones de puntos de medición.
Se recomienda que los datos de entrada estén en un sistema de coordenadas proyectadas en vez de en un sistema de coordenadas geográficas.
-
Un enfoque alternativo es utilizar un dataset de TIN. Primero, cree un TIN a partir de los datos de origen. Luego, convierta el TIN resultante a un ráster con la herramienta De TIN a ráster mediante la opción Vecinos naturales. Esto es particularmente útil si tiene líneas de corte o un área de datos con formas irregulares.
Sintaxis
| 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 |
out_raster |
Ráster de superficie interpolado de salida. | Raster Layer |
cell_size (Opcional) |
El tamaño de celda con el que se creará el ráster de salida. Este será el valor del entorno si se establece explícitamente; de lo contrario, será el valor más bajo del ancho o de la altura de la extensión de las entidades de punto de entrada, en la referencia espacial de entrada, dividido por 250. | Analysis Cell Size |
Ejemplo 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
env.workspace = "C:/data"
arcpy.NaturalNeighbor_3d("ca_ozone_pts.shp", "ozone",
"C:/output/nnout.tif", 2000)
En este ejemplo se introduce un shapefile de punto y se interpola la superficie de salida como ráster de cuadrícula.
# Name: NaturalNeighbor_3d_Ex_02.py
# Description: Interpolate a series of point features onto
# a rectangular raster using Natural Neighbor interpolation.
# Requirements: 3D Analyst Extension
# Import system modules
import arcpy
from arcpy import env
# Set environment settings
env.workspace = "C:/data"
# Set local variables
inPntFeat = "ca_ozone_pts.shp"
zField = "ozone"
outRaster = "C:/output/nnout"
cellSize = 40000
# Check out the ArcGIS 3D Analyst extension license
arcpy.CheckOutExtension("3D")
# Execute NaturalNeighbor
arcpy.NaturalNeighbor_3d(inPntFeat, zField, outRaster, cellSize)