Trend (Spatial Analyst)
Zusammenfassung
Interpoliert eine Raster-Oberfläche anhand von Punkten mithilfe einer Trendmethode.
Weitere Informationen zur Funktionsweise des Werkzeugs "Trend"
Verwendung
-
Wird die Polynom-Ordnung erhöht, wird auch die angepasste Oberfläche zunehmend komplexer. Ein Polynom höherer Ordnung liefert nicht immer die präziseste Oberfläche. Dies hängt von den Daten ab.
-
Die optionale RMS-Dateiausgabe enthält Informationen zum RMS-Fehler der Interpolation. Diese Informationen können verwendet werden, um den optimalen Wert für die Polynom-Ordnung zu ermitteln, indem der Ordnungswert so lange geändert wird, bis Sie den niedrigsten RMS-Fehler erhalten. Informationen zur RMS-Datei finden Sie unter Funktionsweise des Werkzeugs "Trend".
-
Für die Option LOGISTIC des Regressionstyps muss das Z-Wertfeld der Eingabe-Punkt-Features die Codes Null (0) und Eins (1) aufweisen.
-
Einige Eingabe-Datasets weisen mehrere Punkte mit denselben XY-Koordinaten auf. Wenn die Werte der Punkte an der gemeinsamen Position identisch sind, werden sie als Duplikate betrachtet und haben keinerlei Auswirkung auf die Ausgabe. Falls die Werte nicht identisch sind, werden sie als lagegleiche Punkte betrachtet.
Die verschiedenen Interpolationswerkzeuge verarbeiten diese Datenbedingung möglicherweise unterschiedlich. In einigen Fällen wird beispielsweise der erste lagegleiche Punkt für die Berechnung verwendet, während in anderen Fällen der letzte Punkt verwendet wird. Dies kann dazu führen, dass einige Positionen im Ausgabe-Raster andere Werte enthalten als Sie erwarten. Die Lösung besteht darin, Ihre Daten vorzubereiten, indem Sie diese lagegleichen Punkte entfernen. Das Werkzeug Ereignisse erfassen in der Toolbox "Spatial Statistics" unterstützt Sie bei der Ermittlung aller lagegleichen Punkte in Ihren Daten.
Syntax
| Parameter | Erläuterung | Datentyp |
in_point_features |
Die Eingabe-Punkt-Features mit den Z-Werten, die in ein Oberflächen-Raster interpoliert werden. | Feature Layer |
z_field |
Ein Feld mit einem Höhen- oder Größenwert für jeden Punkt. Dies kann ein Zahlenfeld oder das Shape-Feld sein, wenn die Eingabe-Punkt-Features Z-Werte enthalten. Wenn der Regressionstyp LOGISTIC ist, können die Werte im Feld nur 0 oder 1 sein. | Field |
cell_size (optional) |
Die Zellengröße im zu erstellenden Ausgabe-Raster. Dies ist der Wert in der Umgebung, wenn dieser explizit festgelegt wurde. Andernfalls wird die Breite oder Höhe (der kleinere Wert von beiden) der Ausdehnung der Eingabe-Punkt-Features im Eingaberaumbezug dividiert durch 250 verwendet. | Analysis Cell Size |
order (optional) |
Die Polynom-Ordnung. Dies muss eine Ganzzahl von 1 bis 12 sein. Mit dem Wert 1 wird eine flache Ebene an den Punkten angebracht. Ein höherer Wert ergibt eine komplexere Oberfläche. Die Standardeinstellung ist 1. | Long |
regression_type (optional) |
Der Typ der durchzuführenden Regression.
| String |
out_rms_file (optional) |
Der Dateiname für die Ausgabetextdatei, die Informationen zum RMS-Fehler und dem Chi-Quadrat der Interpolation enthält. Die Erweiterung muss TXT sein. | File |
Rückgabewert
| Name | Erläuterung | Datentyp |
| out_raster |
Das Ausgabe-Raster für die interpolierte Oberfläche. | Raster |
Codebeispiel
In diesem Beispiel wird ein Punkt-Shapefile eingegeben und die Ausgabeoberfläche als TIFF-Raster interpoliert.
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")
In diesem Beispiel wird ein Punkt-Shapefile eingegeben und die Ausgabeoberfläche als Grid-Raster interpoliert.
# 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")