Mit der Spatial Analyst-Lizenz verfügbar.
Mit der 3D Analyst-Lizenz verfügbar.
Zusammenfassung
Interpoliert eine Raster-Oberfläche anhand von Punkten mittels IDW (Inverse Distance Weighting).
Verwendung
Der Ausgabewert einer Zelle bei Verwendung von IDW (Inverse Distance Weighting) ist auf den Bereich der Werte begrenzt, die zum Interpolieren verwendet werden. Da "IDW" ein gewichteter Entfernungsdurchschnitt ist, kann der Durchschnitt nicht größer als die höchste und nicht kleiner als die niedrigste Eingabe sein. Aus diesem Grund kann die Funktion keine Bergrücken oder Täler erstellen, wenn diese Extreme nicht bereits erfasst wurden (Watson und Philip, 1985).
Die besten Ergebnisse mit "IDW" werden erzielt, wenn die Referenzpunkterfassung hinsichtlich der lokalen Variation ausreichend dicht ist, die Sie zu simulieren versuchen. Wenn die Erfassung von Eingabepunkten verstreut oder ungleichmäßig ist, bilden die Ergebnisse unter Umständen die gewünschte Oberfläche nicht ausreichend ab (Watson und Philip, 1985).
Der Einfluss eines Eingabepunktes auf einen interpolierten Wert ist isotrop. Da der Einfluss eines Eingabepunktes auf einen interpolierten Wert entfernungsbezogen ist, ist IDW nicht werterhaltend (Philip und Watson, 1982).
Die Ausgabezellengröße kann über einen numerischen Wert definiert oder aus einem vorhandenen Raster-Dataset abgerufen werden. Wird die Zellengröße nicht explizit über den Parameterwert angegeben, wird sie aus der Umgebung "Zellengröße" abgeleitet, falls diese angegeben wurde. Wenn der Parameter und die Umgebung für die Zellengröße nicht angegeben wurden, aber die Umgebung Fang-Raster festgelegt wurde, wird die Zellengröße des Fang-Rasters verwendet. Wenn kein Wert angegeben wird, wird die Zellengröße aus der Breite oder Höhe der Ausdehnung (je nachdem was kürzer ist) berechnet, indem der Wert durch 250 dividiert wird. Dabei wird die Ausdehnung in der Umgebung in Ausgabekoordinatensystem angegeben.
Wenn die Zellengröße mit einem numerischen Wert angegeben wird, wird dieser vom Werkzeug direkt für das Ausgabe-Raster verwendet.
Wenn die Zellengröße mit einem Raster-Dataset angegeben wird, zeigt der Parameter anstelle des Zellengrößenwerts den Pfad des Raster-Datasets an. Die Zellengröße dieses Raster-Datasets wird direkt in der Analyse verwendet, vorausgesetzt der Raumbezug stimmt mit dem Ausgabe-Raumbezug überein. Wenn der Raumbezug des Datasets nicht mit dem Ausgabe-Raumbezug übereinstimmt, wird er basierend auf der ausgewählten Projektionsmethode für Zellengröße projiziert.
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.
Die Option "Eingabe-Barrieren-Features" dient der Angabe der Position linearer Features, die bekanntermaßen die Kontinuität einer Oberfläche unterbrechen. Diese Features haben keine Z-Werte. Klippen, Verwerfungen und Dämme sind typische Beispiele für Barrieren. Barrieren beschränken die ausgewählte Menge der Eingabereferenzpunkte zum Interpolieren von Z-Werten auf diejenigen Referenzpunkte, die sich auf derselben Seite der Barriere wie die aktuell bearbeitete Zelle befinden. Die Trennung durch eine Barriere wird mithilfe einer Sichtfeldanalyse zwischen jedem Punktepaar bestimmt. Dies bedeutet, dass die topologische Trennung für zwei Punkte nicht erforderlich ist, um diese aus dem gegenseitigen Einflussbereich auszuschließen. Eingabereferenzpunkte, die sich exakt auf der Trennlinie befinden, werden in die ausgewählte Referenzmenge für beide Seiten der Barriere einbezogen.
Barrieren-Features werden als Polylinien-Features eingegeben. Die Funktion IDW verwendet für das lineare Feature die XY-Koordinaten. Daher ist es nicht erforderlich, Z-Werte für die linke und rechte Seite der Barriere anzugeben. Angegebene Z-Werte werden ignoriert.
Durch das Verwenden von Barrieren erhöht sich die Verarbeitungsdauer beträchtlich.
Die Begrenzung für dieses Werkzeug liegt bei etwa 45 Millionen Eingabepunkten. Wenn die Eingabe-Feature-Class mehr als 45 Millionen Punkte enthält, erzielen Sie mit diesem Werkzeug kein Ergebnis. Um dies zu vermeiden, können Sie das Untersuchungsgebiet in mehreren Teilen interpolieren, wobei Sie sicherstellen, dass die Kanten bis zu einem gewissen Grad überlappen. Anschließend können Sie die Ergebnisse mosaikartig in einem großen Raster-Dataset zusammenfassen. Alternativ können Sie ein Terrain-Dataset zum Speichern und Visualisieren von Punkten und Oberflächen verwenden, die Milliarden von Messpunkten enthalten.
Wenn Sie über die Erweiterung "Geostatistical Analyst" verfügen, sind Sie möglicherweise in der Lage, größere Datasets mit dem dort verfügbaren Werkzeug IDW zu verarbeiten.
Die Eingabe-Feature-Daten müssen mindestens ein gültiges Feld enthalten.
Quellen:
Philip, G. M., and 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." Geoprocessing 2: 315–327. 1985.
Syntax
Idw(in_point_features, z_field, out_raster, {cell_size}, {power}, {search_radius}, {in_barrier_polyline_features})| Parameter | Erklärung | 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. | Field |
out_raster | Das Ausgabe-Raster für die interpolierte Oberfläche. Es handelt sich stets um ein Gleitkomma-Raster. | Raster Dataset |
cell_size (optional) | Die Zellengröße des Ausgabe-Rasters, das erstellt wird. Dieser Parameter kann über einen numerischen Wert definiert oder aus einem vorhandenen Raster-Dataset abgerufen werden. Wenn die Zellengröße nicht explizit als Parameterwert angegeben wurde, wird der Zellengrößenwert der Umgebung verwendet, sofern dieser angegeben wurde. Andernfalls werden zusätzliche Regeln verwendet, um ihn aus anderen Eingaben zu berechnen. Weitere Informationen finden Sie im Abschnitt zur Verwendung. | Analysis Cell Size |
power (optional) | Der Exponent der Entfernung. Bestimmt die Bedeutung der umgebenden Punkte für den interpolierten Wert. Bei Wahl einer höheren Potenz haben entfernte Punkte weniger Einfluss. Sie können jede reelle Zahl größer als 0 verwenden. Die besten Ergebnisse liefern jedoch Werte von 0,5 bis 3. Die Standardeinstellung ist 2. | Double |
search_radius (optional) | Definiert, welcher der Eingabepunkte zum Interpolieren des Wertes jeder Zelle im Ausgabe-Raster verwendet wird. Es gibt zwei Möglichkeiten, die Suchumgebung anzugeben: Variable und Fixed. Variable verwendet einen variablen Suchradius, um eine bestimmte Anzahl von Eingabereferenzpunkten für die Interpolation zu finden. Fixed verwendet eine bestimmte feste Entfernung, innerhalb der alle Eingabepunkte verwendet werden. Variable ist die Standardeinstellung. Die Syntax für diese Parameter ist wie folgt:
| Radius |
in_barrier_polyline_features (optional) | Polylinien-Features, die als Bruchkante oder Barriere verwendet werden sollen, um die Suche nach Eingabereferenzpunkten zu begrenzen. | Feature Layer |
Codebeispiel
In diesem Beispiel wird ein Punkt-Shapefile eingegeben und die Ausgabeoberfläche als TIFF-Raster interpoliert.
import arcpy
from arcpy import env
env.workspace = "C:/data"
arcpy.Idw_3d("ozone_pts.shp", "ozone", "C:/output/idwout.tif", 2000, 2, 10)In diesem Beispiel wird ein Punkt-Shapefile eingegeben und die Ausgabeoberfläche als Grid-Raster interpoliert.
# Name: IDW_3d_Ex_02.py
# Description: Interpolate a series of point features onto a
# rectangular raster using Inverse Distance Weighting (IDW).
# 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:/output/idwout01"
cellSize = 2000.0
power = 2
searchRadius = 150000
# Execute IDW
arcpy.Idw_3d(inPointFeatures, zField, outRaster, cellSize,
power, searchRadius)Umgebungen
Lizenzinformationen
- Basic: Erfordert 3D Analyst oder Spatial Analyst
- Standard: Erfordert 3D Analyst oder Spatial Analyst
- Advanced: Erfordert 3D Analyst oder Spatial Analyst