Disponible con una licencia de 3D Analyst.
Resumen
Interpola una superficie de ráster a partir de puntos utilizando una técnica de tendencia.
Uso
A medida que aumenta el orden de la función polinómica, la superficie que se ajusta se vuelve cada vez más compleja. Una función polinómica de mayor orden no siempre generará la superficie más precisa, esto depende de los datos.
La salida del archivo RMS opcional contiene información sobre el error RMS (error cuadrático medio) de la interpolación. Esta información se puede utilizar para determinar el mejor valor para el orden de la función polinómica, cambiando el valor del orden hasta obtener el error RMS más bajo. Consulte Cómo funciona Tendencia para obtener información sobre el archivo RMS.
Para la opción LOGISTIC de Tipo de regresión, el campo valor z de las entidades de punto de entrada debe tener códigos de cero (0) y uno (1).
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.
Sintaxis
Trend_3d (in_point_features, z_field, out_raster, {cell_size}, {order}, {regression_type}, {out_rms_file})| 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. Si el tipo de regresión es logística, los valores en el campo sólo pueden ser 0 o 1. | Field |
out_raster | Ráster de superficie interpolado de salida. Siempre es un ráster de punto flotante. | 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 |
order (Opcional) | El orden de la función polinómica. Debe ser un entero entre 1 y 12. Un valor de 1 ajustará un plano chato a los puntos y un valor superior ajustará una superficie más compleja. El valor predeterminado es 1. | Long |
regression_type (Opcional) | El tipo de regresión que se va a realizar.
| String |
out_rms_file (Opcional) | El nombre del archivo para el archivo de texto de salida que contiene información sobre el error RMS y el chi-cuadrado de la interpolación. La extensión debe ser .txt. | File |
Muestra de código
Ejemplo 1 de Tendencia (ventana de Python)
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.Trend_3d("ca_ozone_pts.shp", "ozone",
"C:/output/trendout.tif", 2000, 2, "LINEAR")
Ejemplo 2 de Tendencia (secuencia de comandos independiente)
En este ejemplo se introduce un shapefile de punto y se interpola la superficie de salida como ráster de cuadrícula.
# Name: Trend_3d_Ex_02.py
# Description: Interpolate a series of point features onto a
# rectangular raster using a trend technique.
# Requirements: 3D Analyst Extension
# Import system modules
import arcpy
from arcpy import env
# Set environment settings
env.workspace = "C:/data"
# Set local variables
inPointFeatures = "ca_ozone_pts.shp"
zField = "ozone"
outRaster = "C:/sapyexamples/output/trendout02"
cellSize = 2000.0
PolynomialOrder = 2
regressionType = "LINEAR"
# Check out the ArcGIS 3D Analyst extension license
arcpy.CheckOutExtension("3D")
# Execute Trend
arcpy.Trend_3d(inPointFeatures, zField, outRaster, cellSize,
PolynomialOrder, regressionType)
Entornos
Información sobre licencias
- ArcGIS Desktop Basic: Requiere 3D Analyst o Análisis espacial
- ArcGIS Desktop Standard: Requiere 3D Analyst o Análisis espacial
- ArcGIS Desktop Advanced: Requiere 3D Analyst o Análisis espacial