Trend (Spatial Analyst)
Resumen
Interpola una superficie de ráster a partir de puntos utilizando una técnica de tendencia.
Uso
-
As the order of the polynomial is increased, the surface being fitted becomes progressively more complex. A higher-order polynomial will not always generate the most accurate surface; it is dependent on the data.
-
The optional RMS file output contains information on the RMS (root mean square) error of the interpolation. This information can be used to determine the best value to use for the polynomial order, by changing the order value until you get the lowest RMS error. See How Trend works for information on the RMS file.
-
For the LOGISTIC option of Type of regression, the z-value field of input point features should have codes of zero (0) and one (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 Adquirir eventos de la caja de herramientas de Estadística espacial es útil para identificar cualquier punto coincidente en los datos.
Sintaxis
| Parámetro | Explicación | Tipo de datos |
in_point_features |
Entidades de punto 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. If the regression type is Logistic, the values in the field can only be 0 or 1. | Field |
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) |
The order of the polynomial. This must be an integer between 1 and 12. A value of 1 will fit a flat plane to the points, and a higher value will fit a more complex surface. The default is 1. | Long |
regression_type (Opcional) |
The type of regression to be performed.
| String |
out_rms_file (Opcional) |
File name for the output text file that contains information about the RMS error and the Chi-Square of the interpolation. The extension must be ".txt". | File |
Valor de retorno
| Nombre | Explicación | Tipo de datos |
| out_raster |
Ráster de superficie interpolado de salida. | Raster |
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
from arcpy.sa import *
env.workspace = "C:/sapyexamples/data"
outTrend = Trend("ozone_pts.shp", "ozone", 2000, 2, "LINEAR")
outTrend.save("C:/sapyexamples/output/trendout.tif")
En este ejemplo se introduce un shapefile de punto y se interpola la superficie de salida como ráster de cuadrícula.
# Name: Trend_Ex_02.py
# Description: Interpolate a series of point features
# onto a rectangular raster using a trend technique.
# 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
PolynomialOrder = 2
regressionType = "LINEAR"
# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")
# Execute Trend
outTrend = Trend(inPointFeatures, zField, cellSize,
PolynomialOrder, regressionType)
# Save the output
outTrend.save("C:/sapyexamples/output/trendout02")