IDW 3D (Geostatistical Analyst)

Zusammenfassung

Interpoliert die Werte von 3D-Punkten mithilfe von IDW (Inverse Distance Weighting, inverse Distanzgewichtung) und erstellt einen Voxel-Layer sowie eine Quelldatei (.nc) der vorhergesagten Werte.

Abbildung

Verwendung

• Im Vergleich zum Werkzeug Empirical Bayesian Kriging 3D (EBK 3D), das ebenfalls eine 3D-Interpolation durchführt, ist IDW 3D ein schnelleres und einfacheres Werkzeug, bei dem keine Annahmen zur Verteilung oder zu Trends der Datenwerte getroffen werden. IDW 3D ist eine exakte Interpolationsmethode, d. h. die vorhergesagte 3D-Oberfläche verläuft genau durch die gemessenen Werte von Eingabepunkten. Daher ist dies ein nützliches Visualisierungswerkzeug für unregelmäßige 3D-Punkte.

  IDW 3D erzeugt in der Regel weniger genaue Vorhersagen als EBK 3D und reagiert besonders empfindlich auf geclusterte Eingabepunkte. IDW 3D kann keine Standardfehler vorhergesagter Werte erzeugen. Daher wird eine Schätzung der Unsicherheit der Vorhersagen nicht unterstützt.

• Das Werkzeug sagt Werte an jeder neuen Position anhand eines gewichteten Durchschnitts der Werte der Eingabepunkte, die sich innerhalb der 3D-Suchnachbarschaft der vorhergesagten Position befinden, in 3D vorher. Die Gewichtung für die einzelnen benachbarten Punkte ist die inverse Entfernung (1 geteilt durch die Entfernung) zur vorhergesagten Position mit einer Potenz (Exponent). Die Gewichtungen werden so normalisiert, dass sie im gewichteten Durchschnitt die Summe 1 ergeben.

  Weitere Informationen zu IDW

  • , wobei k die Anzahl der Nachbarn, ω_i die Gewichtung von Nachbar i und z_i der gemessene Wert von Nachbar i ist.

  • , wobei d_i die euklidische 3D-Entfernung zur vorhergesagten Position für Nachbar i und p der Potenzwert ist.

• Wenn das Werkzeug in einer lokalen Szene mit demselben horizontalen und vertikalen Koordinatensystem wie in den Eingabe-Features verwendet wird, wird der Szene ein Voxel-Layer hinzugefügt, in dem Sie die Ergebnisse interaktiv auswerten können. Sie können auch die Ausgabe-netCDF-Datei als Voxel-Layer hinzufügen. Verwenden Sie dazu das Werkzeug Multidimensionalen Voxel-Layer erstellen oder das Dialogfeld Multidimensionalen Voxel-Layer hinzufügen.

  Sie können die Ausgabe-netCDF-Datei mit dem Werkzeug Raster kopieren in ein multidimensionales Raster konvertieren. Sie können sie auch mit dem Werkzeug NetCDF-Feature-Layer erstellen oder mit dem Werkzeug NetCDF-Raster-Layer erstellen einer Karte oder Szene als Feature- bzw. Raster-Layer hinzufügen.

• Summenstatistiken auf Grundlage der Leave-One-Out-Kreuzvalidierung werden als Geoverarbeitungsmeldungen angezeigt, um die Genauigkeit und Zuverlässigkeit der Vorhersagen zu bewerten. Die folgenden Summenstatistiken werden angezeigt:

  • Anzahl: Die Anzahl der Features mit Kreuzvalidierungsergebnissen. Dieser Wert kann sich von der Anzahl der Eingabe-Features unterscheiden, wenn einige Features NULL-Werte oder lagegleichen Positionen aufweisen oder wenn keine benachbarten Features gefunden werden können.
  • Mean Error: Der Mittelwert der Kreuzvalidierungsfehler. Mit dieser Statistik wird die Modellverzerrung gemessen. Dieser Wert sollte so nah an Null sein wie möglich. Positive Werte weisen auf eine Tendenz zur Überschätzung (vorhergesagte Werte sind größer als die gemessenen Werte) und negative Werte auf eine Tendenz zur Unterschätzung hin.
  • Root Mean Square Error: Die Quadratwurzel der durchschnittlichen quadrierten Kreuzvalidierungsfehler. Mit dieser Statistik wird die Vorhersagegenauigkeit gemessen. Der Wert sollte möglichst klein sein. Der Wert stellt eine Schätzung der durchschnittlichen Differenz zwischen den vorhergesagten und den gemessenen Werten dar. Für die Interpolation der Temperatur in Grad Celsius bedeutet ein RMS-Fehler von 1,5 zum Beispiel, dass die Vorhersagen voraussichtlich um durchschnittlich etwa 1,5 Grad von den tatsächlichen Werten abweichen.

  Weitere Informationen zur Kreuzvalidierung

• Die Eingabe-Features müssen 3D-Punkte mit Höhen sein, die im Geometrieattribut Shape.Z gespeichert sind. Mit dem Werkzeug Feature in 3D nach Attribut können Sie 2D-Punkt-Features mit einem Höhenfeld in 3D-Punkt-Features konvertieren.

  Es wird empfohlen, dass die Eingabe-Features ein vertikales Koordinatensystem verwenden, das die Z-Koordinaten genau definiert. Sie können den Punkten mit dem Werkzeug Projektion definieren ein vertikales Koordinatensystem zuweisen.

• Verwenden Sie den Parameter Ausgabe-Feature-Class der Kreuzvalidierung, um die Kreuzvalidierungsfehler für jeden Eingabepunkt zu untersuchen. Die gemessenen Werte und die Vorhersagen der Kreuzvalidierung werden als Felder in der Feature-Class gespeichert.

  Die Feature-Class enthält zwei Scatterplots zur Untersuchung von Trends in den Ergebnissen der Kreuzvalidierung:

  • Kreuzvalidierung: "Vorhergesagt" gegenüber "Gemessen": Zeigt die Vorhersagen der Kreuzvalidierung im Vergleich zu den gemessenen Werten an. Wenn die vorhergesagten Werte etwa den gemessenen Werten entsprechen (was auf genaue Interpolationsergebnisse schließen lässt), sollten die Punkte im Scatterplot eine Linie mit der Neigung 1 bilden.
  • Kreuzvalidierung: "Gemessen" gegenüber "Fehler": Zeigt die gemessenen Werte im Vergleich zu den Kreuzvalidierungsfehlern an. Wenn die Fehler von den gemessenen Werten unabhängig sind, weisen die Punkte im Scatterplot keine Muster oder Trends auf, und die Trendlinie ist flach (Neigung etwa gleich 0). Trendlinien mit negativer (abnehmender) Neigung deuten auf einen Glättungseffekt im Interpolationsmodell hin, d. h. das Modell weist eine Tendenz zur Unterschätzung großer und Überschätzung kleiner Werte auf.

  Weitere Informationen zur Interpretation von Kreuzvalidierungsdiagrammen

• Die Eingabe-Features und die Ausschneideraster für minimale und maximale Höhe müssen sich in einem projizierten Koordinatensystem befinden. Wenn die Punkte oder Raster ein geographisches Koordinatensystem mit Koordinaten für Breitengrad und Längengrad haben, müssen sie mit dem Werkzeug Projizieren oder Raster projizieren projiziert werden.

• Beim Anlegen des 3D-Punktgitters, das die Voxel repräsentiert, wird der erste Punkt an der minimalen X-, minimalen Y- und minimalen Z-Koordinate der Ausgabeausdehnung erstellt (standardmäßig die Ausdehnung der Eingabe-Features). Die verbleibenden Punkte werden erstellt, indem unter Verwendung der Parameter X-Abstand, Y-Abstand, und Höhenabstand die Entfernungen in den einzelnen Dimensionen der Ausgabeausdehnung iterativ durchlaufen werden. Wenn sich mit einem dieser Abstände die entsprechende Dimension der Ausgabeausdehnung nicht gleichmäßig aufteilen lässt, wird jenseits der Ausgabeausdehnung eine Zeile oder Spalte mit Punkten erstellt. Wenn zum Beispiel für X die Ausgabeausdehnung 0 bis 10 angegeben wurde und der Parameter X-Abstand auf 3 festgelegt wurde, dann enthält die Ausgabe fünf Zeilen in der X-Ausdehnung: 0, 3, 6, 9 und 12. Wenn sich mit den Abständen die X- oder die Y-Ausdehnung nicht gleichmäßig aufteilen lässt, wird eine zusätzliche Zeile bzw. Spalte mit Punkten erstellt.

• Mit den Parametern Eingabe-Untersuchungsgebiet-Polygone, Ausschneide-Raster für minimale Höhe und Ausschneide-Raster für maximale Höhe kann die Analyse auf ein bestimmtes Untersuchungsgebiet und zwischen zwei Höhenoberflächen begrenzt werden. Alle Voxel außerhalb dieser Grenzen haben keinen Wert und werden nicht angezeigt. Wenn die Punkte beispielsweise in einem Meeresschutzgebiet liegen, können Sie einen Voxel-Layer erstellen, der nur innerhalb eines Polygons des Schutzgebietes (Untersuchungsgebiet), über dem Meeresboden (Raster für minimale Höhen) und unter der Thermokline (Raster für maximale Höhe) angezeigt wird.

  Es gibt verschiedene Überlegungen zur Verwendung von Höhenoberflächen als Raster für minimale oder maximale Höhe. Image-Services, Web-Höhen-Layer und Web-Bilddaten-Layer haben die langsamste Performance. Bei einer großen Anzahl von Abfragen können Fehler auftreten. Raster, die als lokale Dateien auf der Festplatte gespeichert werden, haben die schnellste Performance. Sie werden empfohlen, wenn Voxel-Layer mit hoher Auflösung über große räumliche Ausdehnungen erstellt werden sollen.

• Wenn die Eingabe-Features eine Auswahl haben, werden die Parameter X-Abstand, Y-Abstand und Höhenabstand neu berechnet, während das Werkzeug ausgeführt wird, basierend auf der Ausdehnung der selektierten Features. Die neu berechneten Werte werden als Warnmeldungen ausgedruckt, wenn das Werkzeug beendet wird. Wenn Sie einen Wert für einen Abstandsparameter manuell angeben (oder eine Ausgabeausdehnung angeben), wird der Wert nicht neu berechnet.

• Wenn Eingabe-Untersuchungsgebiet-Polygone angegeben wurden, dann wird die Ausdehnung des Untersuchungsgebiets als Standard-Ausgabeausdehnung verwendet. Anschließend werden basierend auf dieser Ausdehnung die Werte der Parameter X-Abstand und Y-Abstand neu berechnet. Dadurch wird sichergestellt, dass die Ausgabe standardmäßig das gesamte Untersuchungsgebiet ausfüllt.

Parameter

Dialogfeld

Beschriftung	Erläuterung	Datentyp
Eingabe-Features	Die 3D-Punkt-Features, die das zu interpolierende Feld enthalten. Die Punkte müssen sich in einem projizierten Koordinatensystem befinden.	Feature Layer
Wertfeld	Das Feld aus den Eingabe-Features mit den gemessenen Werten, die interpoliert werden sollen.	Field
Ausgabe-netCDF-Datei	Die Ausgabe-netCDF-Datei, die die vorhergesagten Werte in einem 3D-Gitternetz enthält. Diese Datei kann als Datenquelle eines Voxel-Layers verwendet werden.	File
Potenz (optional)	Der Potenzwert, der beim Berechnen von Vorhersagen zur Gewichtung der Werte von benachbarten Features verwendet wird. Bei Wahl einer höheren Potenz haben näher gelegene Punkte mehr Einfluss. Der Wert muss zwischen 1 und 100 liegen. Die Standardeinstellung ist 2.	Double
Höheninflationsfaktor (optional)	Ein konstanter Wert, der vor dem Suchen von Nachbarn und Berechnen von Entfernungen mit den Z-Koordinaten der Eingabe-Features multipliziert wird. Bei den meisten 3D-Daten ändern sich die Werte der Punkte schneller vertikal als horizontal. Dieser Faktor streckt die Positionen der Punkte, sodass eine vertikale Entfernungseinheit einer horizontalen Entfernungseinheit entspricht. Die Positionen der Punkte werden zurück an ihre ursprünglichen Positionen verschoben, bevor das Ergebnis der Interpolation zurückgegeben wird. Wird kein Wert angegeben, wird er während der Ausführung des Werkzeugs geschätzt und als Geoverarbeitungsmeldung angezeigt. Der geschätzte Wert wird durch Minimierung des RMS-Kreuzvalidierungsfehlers ermittelt. Der Wert muss zwischen 1 und 1.000 liegen.	Double
Ausgabe-Feature-Class der Kreuzvalidierung (optional)	Eine Feature-Class der Kreuzvalidierungsstatistik für jeden Eingabepunkt. Die Feature-Class enthält zwei Scatterplots.	Feature Class
X-Abstand (optional)	Der Abstand zwischen den einzelnen Gitterpunkten in der X-Dimension. Mit dem Standardwert werden 40 Punkte entlang der X-Ausdehnung der Ausgabe erstellt.	Linear Unit
Y-Abstand (optional)	Der Abstand zwischen den einzelnen Gitterpunkten in der Y-Dimension. Mit dem Standardwert werden 40 Punkte entlang der Y-Ausdehnung der Ausgabe erstellt.	Linear Unit
Höhenabstand (optional)	Der Abstand zwischen den einzelnen Gitterpunkten in der Höhendimension (Z). Mit dem Standardwert werden 40 Punkte entlang der Z-Ausdehnung der Ausgabe erstellt.	Linear Unit
Eingabe-Untersuchungsgebiet-Polygone (optional)	Die Polygon-Features, die das Untersuchungsgebiet darstellen. In der Ausgabe-netCDF-Datei werden nur die Punkte gespeichert, die sich innerhalb des Untersuchungsgebietes befinden. Bei Visualisierung als Voxel-Layer werden nur die Voxel innerhalb des Untersuchungsgebietes in der Szene dargestellt. Ob sich Punkte innerhalb oder außerhalb des Untersuchungsgebietes befinden, wird anhand ihrer X- und Y-Koordinaten bestimmt.	Feature Layer
Ausschneide-Raster für minimale Höhe (optional)	Das zum Ausschneiden des unteren Teils des Voxel-Layers verwendete Höhenraster. Nur Voxeln über diesem Höhenraster werden Vorhersagen zugewiesen. Wenn Sie z. B. ein Geländehöhenraster verwenden, wird der Voxel-Layer nur oberirdisch angezeigt. Es kann auch für die Oberfläche von Grundgestein oder den Boden einer Schiefer-Lagerstätte verwendet werden. Das Raster muss sich in einem projizierten Koordinatensystem befinden, und die Höhenwerte müssen in derselben Einheit vorliegen wie die vertikale Einheit des Rasters.	Raster Layer
Ausschneide-Raster für maximale Höhe (optional)	Das zum Ausschneiden des oberen Teils des Voxel-Layers verwendete Höhenraster. Nur Voxeln unter diesem Höhenraster werden Vorhersagen zugewiesen. Wenn Sie z. B. ein Geländehöhenraster verwenden, wird der Voxel-Layer nur unterirdisch angezeigt. Es kann auch zum Ausschneiden von Voxeln oberhalb eines eingeschränkten Luftraums verwendet werden. Das Raster muss sich in einem projizierten Koordinatensystem befinden, und die Höhenwerte müssen in derselben Einheit vorliegen wie die vertikale Einheit des Rasters.	Raster Layer
Suchnachbarschaft (optional)	Gibt die Anzahl und Ausrichtung der Nachbarn an, die zum Vorhersagen von Werten an neuen Positionen verwendet werden sollen. Standard 3D Max. Anzahl Nachbarn: Die maximale Anzahl der Nachbarn pro Sektor, die zum Schätzen des Wertes an der unbekannten Position verwendet werden sollen. Min. Anzahl Nachbarn: die minimale Anzahl der Nachbarn pro Sektor, die zum Schätzen des Wertes an der unbekannten Position verwendet werden sollen. Sektortyp: Die Geometrie der 3D-Nachbarschaft. Sektoren werden verwendet, um sicherzustellen, dass Nachbarn in jeder Richtung um die Vorhersageposition verwendet werden. Sektoren aller Typen werden aus den platonischen Körpern gebildet. 1 Sektor (Kugel): Verwendet werden die nächstgelegenen Nachbarn in jeder Richtung. 4 Sektoren (Tetraeder): Der Raum wird in 4 Regionen geteilt. Verwendet werden dann Nachbarn in jeder dieser 4 Regionen. 6 Sektoren (Würfel): Der Raum wird in 6 Regionen geteilt. Verwendet werden dann Nachbarn in jeder dieser 6 Regionen. 8 Sektoren (Oktaeder): Der Raum wird in 8 Regionen geteilt. Verwendet werden dann Nachbarn in jeder dieser 8 Regionen. 12 Sektoren (Dodekaeder): Der Raum wird in 12 Regionen geteilt. Verwendet werden dann Nachbarn in jeder dieser 12 Regionen. 20 Sektoren (Ikosaeder): Der Raum wird in 20 Regionen geteilt. Verwendet werden dann Nachbarn in jeder dieser 20 Regionen. Radius: Die Länge des Radius der Suchnachbarschaft. Wird kein Wert angegeben, wird während der Ausführung des Werkzeugs ein Wert geschätzt und als Geoverarbeitungsmeldung angezeigt.	Geostatistical Search Neighborhood

Abgeleitete Ausgabe

Beschriftung	Erläuterung	Datentyp
Anzahl	Die Gesamtzahl der verwendeten Stichproben.	Long
Mean Error	Die gemittelte Differenz zwischen den gemessenen und den vorhergesagten Werten.	Double
Root Mean Square	Gibt an, wie gut das Modell die gemessenen Werte vorhersagt.	Double
Ausgabe-Voxel-Layer	Ein Voxel-Layer der vorhergesagten Werte.	Voxel Layer

Python

arcpy.ga.IDW3D(in_features, value_field, out_netcdf_file, {power}, {elev_inflation_factor}, {out_cv_features}, {x_spacing}, {y_spacing}, {elevation_spacing}, {in_study_area}, {min_elev_raster}, {max_elev_raster}, {search_neighborhood})

Name	Erläuterung	Datentyp
in_features	Die 3D-Punkt-Features, die das zu interpolierende Feld enthalten. Die Punkte müssen sich in einem projizierten Koordinatensystem befinden.	Feature Layer
value_field	Das Feld aus den Eingabe-Features mit den gemessenen Werten, die interpoliert werden sollen.	Field
out_netcdf_file	Die Ausgabe-netCDF-Datei, die die vorhergesagten Werte in einem 3D-Gitternetz enthält. Diese Datei kann als Datenquelle eines Voxel-Layers verwendet werden.	File
power (optional)	Der Potenzwert, der beim Berechnen von Vorhersagen zur Gewichtung der Werte von benachbarten Features verwendet wird. Bei Wahl einer höheren Potenz haben näher gelegene Punkte mehr Einfluss. Der Wert muss zwischen 1 und 100 liegen. Die Standardeinstellung ist 2.	Double
elev_inflation_factor (optional)	Ein konstanter Wert, der vor dem Suchen von Nachbarn und Berechnen von Entfernungen mit den Z-Koordinaten der Eingabe-Features multipliziert wird. Bei den meisten 3D-Daten ändern sich die Werte der Punkte schneller vertikal als horizontal. Dieser Faktor streckt die Positionen der Punkte, sodass eine vertikale Entfernungseinheit einer horizontalen Entfernungseinheit entspricht. Die Positionen der Punkte werden zurück an ihre ursprünglichen Positionen verschoben, bevor das Ergebnis der Interpolation zurückgegeben wird. Wird kein Wert angegeben, wird er während der Ausführung des Werkzeugs geschätzt und als Geoverarbeitungsmeldung angezeigt. Der geschätzte Wert wird durch Minimierung des RMS-Kreuzvalidierungsfehlers ermittelt. Der Wert muss zwischen 1 und 1.000 liegen.	Double
out_cv_features (optional)	Eine Feature-Class der Kreuzvalidierungsstatistik für jeden Eingabepunkt. Die Feature-Class enthält zwei Scatterplots.	Feature Class
x_spacing (optional)	Der Abstand zwischen den einzelnen Gitterpunkten in der X-Dimension. Mit dem Standardwert werden 40 Punkte entlang der X-Ausdehnung der Ausgabe erstellt.	Linear Unit
y_spacing (optional)	Der Abstand zwischen den einzelnen Gitterpunkten in der Y-Dimension. Mit dem Standardwert werden 40 Punkte entlang der Y-Ausdehnung der Ausgabe erstellt.	Linear Unit
elevation_spacing (optional)	Der Abstand zwischen den einzelnen Gitterpunkten in der Höhendimension (Z). Mit dem Standardwert werden 40 Punkte entlang der Z-Ausdehnung der Ausgabe erstellt.	Linear Unit
in_study_area (optional)	Die Polygon-Features, die das Untersuchungsgebiet darstellen. In der Ausgabe-netCDF-Datei werden nur die Punkte gespeichert, die sich innerhalb des Untersuchungsgebietes befinden. Bei Visualisierung als Voxel-Layer werden nur die Voxel innerhalb des Untersuchungsgebietes in der Szene dargestellt. Ob sich Punkte innerhalb oder außerhalb des Untersuchungsgebietes befinden, wird anhand ihrer X- und Y-Koordinaten bestimmt.	Feature Layer
min_elev_raster (optional)	Das zum Ausschneiden des unteren Teils des Voxel-Layers verwendete Höhenraster. Nur Voxeln über diesem Höhenraster werden Vorhersagen zugewiesen. Wenn Sie z. B. ein Geländehöhenraster verwenden, wird der Voxel-Layer nur oberirdisch angezeigt. Es kann auch für die Oberfläche von Grundgestein oder den Boden einer Schiefer-Lagerstätte verwendet werden. Das Raster muss sich in einem projizierten Koordinatensystem befinden, und die Höhenwerte müssen in derselben Einheit vorliegen wie die vertikale Einheit des Rasters.	Raster Layer
max_elev_raster (optional)	Das zum Ausschneiden des oberen Teils des Voxel-Layers verwendete Höhenraster. Nur Voxeln unter diesem Höhenraster werden Vorhersagen zugewiesen. Wenn Sie z. B. ein Geländehöhenraster verwenden, wird der Voxel-Layer nur unterirdisch angezeigt. Es kann auch zum Ausschneiden von Voxeln oberhalb eines eingeschränkten Luftraums verwendet werden. Das Raster muss sich in einem projizierten Koordinatensystem befinden, und die Höhenwerte müssen in derselben Einheit vorliegen wie die vertikale Einheit des Rasters.	Raster Layer
search_neighborhood (optional)	Gibt die Anzahl und Ausrichtung der Nachbarn bei Verwendung der Klasse SearchNeighborhoodStandard3D an. Standard 3D radius: Die Länge des Radius der Suchnachbarschaft. Wird kein Wert angegeben, wird während der Ausführung des Werkzeugs ein Wert geschätzt und als Geoverarbeitungsmeldung angezeigt. nbrMax: Die maximale Anzahl der Nachbarn pro Sektor, die zum Schätzen des Wertes an der unbekannten Position verwendet werden sollen. nbrMin: Die minimale Anzahl der Nachbarn pro Sektor, die zum Schätzen des Wertes an der unbekannten Position verwendet werden sollen. sectorType: Die Geometrie der 3D-Nachbarschaft. Sektoren werden verwendet, um sicherzustellen, dass Nachbarn in verschiedenen Richtungen um die Vorhersageposition verwendet werden. Sektoren aller Typen werden aus den platonischen Körpern gebildet. ONE_SECTOR: Verwendet werden die nächstgelegenen Nachbarn in jeder Richtung. FOUR_SECTORS: Der Raum wird in 4 Regionen geteilt. Verwendet werden dann Nachbarn in jeder dieser 4 Regionen. SIX_SECTORS: Der Raum wird in 6 Regionen geteilt. Verwendet werden dann Nachbarn in jeder dieser 6 Regionen. EIGHT_SECTORS: Der Raum wird in 8 Regionen geteilt. Verwendet werden dann Nachbarn in jeder dieser 8 Regionen. TWELVE_SECTORS: Der Raum wird in 12 Regionen geteilt. Verwendet werden dann Nachbarn in jeder dieser 12 Regionen. TWENTY_SECTORS: Der Raum wird in 20 Regionen geteilt. Verwendet werden dann Nachbarn in jeder dieser 20 Regionen.	Geostatistical Search Neighborhood

Abgeleitete Ausgabe

Name	Erläuterung	Datentyp
count	Die Gesamtzahl der verwendeten Stichproben.	Long
mean_error	Die gemittelte Differenz zwischen den gemessenen und den vorhergesagten Werten.	Double
root_mean_square	Gibt an, wie gut das Modell die gemessenen Werte vorhersagt.	Double
out_voxel_layer	Ein Voxel-Layer der vorhergesagten Werte.	Voxel Layer

Codebeispiel

IDW3D –Beispiel 1 (Python-Fenster)

Das folgende Python-Skript veranschaulicht, wie die Funktion IDW3D verwendet wird.

# Interpolate 3d oxygen measurements using IDW3D

arcpy.ga.IDW3D("OxygenPoints3D", "OxygenValue","outputNCDF.nc", "2",
         "", "outputCV.fc", "50 Meters", "50 Meters", "5 Meters",
         "MyStudyArea", "MinElevationRaster", "MaxElevationRaster",
         "NBRTYPE=Standard3D RADIUS=nan NBR_MAX=2 NBR_MIN=1 SECTOR_TYPE=TWELVE_SECTORS")

IDW3D –Beispiel 2 (eigenständiges Skript)

Das folgende Python-Skript veranschaulicht, wie die Funktion IDW3D verwendet wird.

# Name: IDW3D_Example_02.py
# Description: Creates a voxel layer source file from interpolated 3D points.
# Requirements: Geostatistical Analyst Extension
# Author: Esri



# Import system modules
import arcpy

# Allow overwriting output
arcpy.env.overwriteOutput = True

# Define 3D input points and value field to be interpolated
in3DPoints = "C:/gapydata/inputs.gdb/myOxygenPoints3D"
valueField = "OxygenValue"
outNetCDF = "C:/gapydata/outputs/OxygenMeasurementsVoxel.nc"
outCVFeatureClass = "C:/gapydata/outputs/outputCrossValidationErr.shp"

# Define power of IDW and elevation inflation factor
powerValue = "2"
elevInflation = ""

# Define voxel dimensions
xSpacing = "50 Meters"
ySpacing = "50 Meters"
elevSpacing = "5 Meters"


# Define study area, minimum clipping raster layer, and maximum clipping elevation layer
studyArea = "C:/gapydata/inputs.gdb/StudyAreaPolygon"
minElevRaster = "C:/gapydata/inputs.gdb/MinElevationClippingRaster"
maxElevRaster = "C:/gapydata/inputs.gdb/MaxElevationClippingRaster"

# Define the neighborhood
radius = ""
maxNeighbors = 2
minNeighbors = 1
sectorType = "TWELVE_SECTORS"
searchNeighborhood = arcpy.SearchNeighborhoodStandard3D(radius, maxNeighbors,
                     minNeighbors, sectorType)



# Check out the ArcGIS Geostatistical Analyst extension license
arcpy.CheckOutExtension("GeoStats")

# Execute IDW 3D
arcpy.ga.IDW3D(in3DPoints, valueField,outNetCDF,
                           powerValue, elevInflation, outCVFeatureClass,
                           xSpacing, ySpacing, elevSpacing,
                           studyArea, minElevRaster,
                           maxElevRaster, searchNeighborhood)

# Print messages
print(arcpy.GetMessages())

Umgebungen

Lagegleiche Punkte, Ausdehnung, Faktor für parallele Verarbeitung

Lizenzinformationen

• Basic: Erfordert Geostatistical Analyst
• Standard: Erfordert Geostatistical Analyst
• Advanced: Erfordert Geostatistical Analyst