3D-GA-Layer in NetCDF (Geostatistical Analyst)

Mit der Geostatistical Analyst-Lizenz verfügbar.

Zusammenfassung

Exportiert geostatistische 3D-Layer, die mit dem Werkzeug Empirical Bayesian Kriging 3D erstellt wurden, in das netCDF-Format (*.nc-Datei). Die Ausgabedatei wird als Voxel-Layer in einer lokalen Szene angezeigt.

Weitere Informationen zu Voxel-Layern

Abbildung

Der Voxel-Layer zeigt die Ergebnisse einer 3D-Interpolation an.

Verwendung

• 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.

• Alle geostatistischen Eingabe-Layer müssen geostatistische Layer in 3D sein, die mit Empirical Bayesian Kriging 3D erstellt wurden.

• Zum Auswählen der Positionen für Vorhersage und Export haben Sie die folgenden Möglichkeiten:

  • Vorhersage in 3D-Gitterpunkte: Die Vorhersagen werden in Gitterpunkte in 3D exportiert. Damit diese Option verwendet wird, geben Sie 3D-Gitterpunkte im Parameter Export-Positionen an. Die Parameter X-Abstand, Y-Abstand und Höhenabstand geben die Entfernung zwischen den einzelnen Gitterpunkten in den Dimensionen X, Y und Z (Höhe) an.
  • Vorhersage in benutzerdefinierte 3D-Positionen: Die Vorhersagen werden in benutzerdefinierte 3D-Positionen, die durch 3D-Punkt-Features definiert werden, exportiert. Damit diese Option verwendet wird, geben Sie Benutzerdefinierte 3D-Punkte im Parameter Export-Positionen an und fügen dem Parameter 3D-Punktpositionen die 3D-Punkt-Features, die die vorherzusagenden Positionen angeben, hinzu.

    Tipp:

    Sie können das Werkzeug Feature zu 3D nach Attribut verwenden, um aus 2D-Punkten mit Höhenfeld 3D-Punkte zu erstellen.

• Als Standard-Ausgabeausdehnung der Ausgabe-netCDF-Datei wird die Vereinigung der Ausdehnungen aller geostatistischen Eingabe-Layer in 3D verwendet.

• 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.

• 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.

• 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.

• 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
Geostatistische Eingabe-Layer in 3D	Die geostatistischen 3D-Layer, die in eine Ausgabe-netCDF-Datei exportiert werden sollen.	Geostatistical Layer
Ausgabe-netCDF-Datei	Die Ausgabe-netCDF-Datei, die die exportierten Werte aus den geostatistischen Eingabe-Layern enthält. Die Ergebnisse der einzelnen geostatistischen Layer werden als unterschiedliche Variablen in der netCDF-Datei gespeichert.	File
Export-Positionen (optional)	Gibt die Positionen an, die aus dem Parameterwert Geostatistische Eingabe-Layer in 3D exportiert werden sollen. Wenn Sie die Option Gitterpunkte in 3D auswählen, müssen Sie für die Parameter X-Abstand, Y-Abstand und Höhenabstand die Werte angeben, die die Entfernung zwischen den einzelnen Gitterpunkten in allen Dimensionen darstellen. Wenn Sie die Option Benutzerdefinierte 3D-Punkte auswählen, müssen Sie im Parameter 3D-Punktpositionen die 3D-Punkt-Features angeben, die die zu exportierenden Positionen darstellen. 3D-Gitterpunkte—Die Vorhersage-Positionen sind 3D-Gitterpunkte. Dies ist die Standardeinstellung. Benutzerdefinierte 3D-Punkte—Die Vorhersage-Positionen werden durch benutzerdefinierte 3D-Punkt-Features definiert.	String
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
3D-Punktpositionen (optional)	Die 3D-Punkt-Features, die die zu exportierenden Positionen darstellen. Die Höhen der Punkt-Features müssen im Geometrieattribut in "Shape.X" gespeichert sein.	Feature Layer
Ausgabevariablen (optional)	Gibt die Ausgabetypen für die Werte von Geostatistische Eingabe-Layer in 3D an. Sie können für jeden dieser Layer einen oder mehrere Ausgabetypen angeben oder einen Ausgabetyp für alle geostatistischen Eingabe-Layer verwenden. Standardmäßig werden die Vorhersagen für alle Layer exportiert. Wenn andere Ausgabetypen exportiert werden sollen, dann müssen Sie im ersten Eintrag der Wertetabelle den zu exportierenden Layer angeben (oder Alle auswählen, um alle Layer anzugeben). Geben Sie im zweiten Eintrag der Wertetabelle den Ausgabetyp an. Wenn Sie Wahrscheinlichkeit oder Quantil als Ausgabetyp auswählen, dann müssen Sie im dritten Eintrag der Wertetabelle den Schwellenwert (für Wahrscheinlichkeit) bzw. den Quantilwert (für Quantil) angeben. Wenn Sie Vorhersage oder Standardfehler der Vorhersage als Ausgabetyp auswählen, dann können Sie den dritten Eintrag in der Wertetabelle leer lassen. Weitere Informationen zu den geostatistischen Ausgabetypen	Value Table
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

Abgeleitete Ausgabe

Beschriftung	Erläuterung	Datentyp
Ausgabe-Voxel-Layer	Ein Voxel-Layer der vorhergesagten Werte.	Voxel Layer

Python

arcpy.ga.GALayer3DToNetCDF(in_3d_geostat_layers, out_netcdf_file, {export_locations}, {x_spacing}, {y_spacing}, {elevation_spacing}, {in_points_3d}, {output_variables}, {in_study_area}, {min_elev_raster}, {max_elev_raster})

Name	Erläuterung	Datentyp
in_3d_geostat_layers [in_3d_geostat_layers,...]	Die geostatistischen 3D-Layer, die in eine Ausgabe-netCDF-Datei exportiert werden sollen.	Geostatistical Layer
out_netcdf_file	Die Ausgabe-netCDF-Datei, die die exportierten Werte aus den geostatistischen Eingabe-Layern enthält. Die Ergebnisse der einzelnen geostatistischen Layer werden als unterschiedliche Variablen in der netCDF-Datei gespeichert.	File
export_locations (optional)	Gibt die Positionen an, die aus dem Parameterwert in_3d_geostat_layers exportiert werden sollen. Wenn Sie die Option 3D_GRIDDED_POINTS auswählen, müssen Sie für die Parameter x_spacing, u_spacing und elevation_spacing die Werte angeben, die die Entfernung zwischen den einzelnen Gitterpunkten in allen Dimensionen darstellen. Wenn Sie die Option CUSTOM_3D_POINTS auswählen, müssen Sie im Parameter in_points_3d die 3D-Punkt-Features angeben, die die zu exportierenden Positionen darstellen. 3D_GRIDDED_POINTS—Die Vorhersage-Positionen sind 3D-Gitterpunkte. Dies ist die Standardeinstellung. CUSTOM_3D_POINTS—Die Vorhersage-Positionen werden durch benutzerdefinierte 3D-Punkt-Features definiert.	String
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_points_3d (optional)	Die 3D-Punkt-Features, die die zu exportierenden Positionen darstellen. Die Höhen der Punkt-Features müssen im Geometrieattribut in "Shape.X" gespeichert sein.	Feature Layer
output_variables [[layer_name, output_type, quantile_probability_value],...] (optional)	Gibt die Ausgabetypen für die Werte von Geostatistische Eingabe-Layer in 3D an. Sie können für jeden dieser Layer einen oder mehrere Ausgabetypen angeben oder einen Ausgabetyp für alle geostatistischen Eingabe-Layer verwenden. Standardmäßig werden die Vorhersagen für alle Layer exportiert. Wenn andere Ausgabetypen exportiert werden sollen, dann müssen Sie im ersten Eintrag der Wertetabelle den zu exportierenden Layer angeben (oder Alle auswählen, um alle Layer anzugeben). Geben Sie im zweiten Eintrag der Wertetabelle den Ausgabetyp an. Wenn Sie Wahrscheinlichkeit oder Quantil als Ausgabetyp auswählen, dann müssen Sie im dritten Eintrag der Wertetabelle den Schwellenwert (für Wahrscheinlichkeit) bzw. den Quantilwert (für Quantil) angeben. Wenn Sie Vorhersage oder Standardfehler der Vorhersage als Ausgabetyp auswählen, dann können Sie den dritten Eintrag in der Wertetabelle leer lassen. Weitere Informationen zu den geostatistischen Ausgabetypen	Value Table
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

Abgeleitete Ausgabe

Name	Erläuterung	Datentyp
out_voxel_layer	Ein Voxel-Layer der vorhergesagten Werte.	Voxel Layer

Codebeispiel

GALayer3DToNetCDF: Beispiel 1 (Python-Fenster)

Zweimaliges Interpolieren von 3D-Punkten und Konvertieren der Ausgabe in eine multivariate netCDF-Datei.

import arcpy
arcpy.ga.EmpiricalBayesianKriging3D("my3DPoints1", "Shape.Z",
                                    "myValueField1", "my3DGALayer1")
arcpy.ga.EmpiricalBayesianKriging3D("my3DPoints2", "Shape.Z",
                                    "myValueField2", "my3DGALayer2")
arcpy.ga.GALayer3DToNetCDF("my3DGALayer1;my3DGALayer2", "outputNCDF.nc",
                           "3D_GRIDDED_POINTS", "50 Meters", "50 Meters", "5 Meters",
                           "", "<ALL> PREDICTION #")

GALayer3DToNetCDF: Beispiel 2 (eigenständiges Skript)

Zweimaliges Interpolieren von 3D-Punkten und Konvertieren der Ausgabe in eine multivariate netCDF-Datei. Exportieren in Gitterpunkte und benutzerdefinierte 3D-Punkte

# Name: GALayer3DToNetCDF_Example_02.py
# Description: Interpolates 3D points and exports to a netCDF file.
# Requirements: Geostatistical Analyst Extension
# Author: Esri

# Import system modules
import arcpy

# Allow overwriting output
arcpy.env.overwriteOutput = True

# Set up variables
in3DPoints1 = "C:/gapydata/inputs.gdb/my3DPoints1"
in3DPoints2 = "C:/gapydata/inputs.gdb/my3DPoints2"
elevationField1 = "Shape.Z"
elevationField2 = "Shape.Z"
valueField1 = "myValueField1"
valueField2 = "myValueField2"
outGALayer1 = "myGALayer1"
outGALayer2 = "myGALayer2"


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

# Execute Empirical Bayesian Kriging 3D twice
arcpy.ga.EmpiricalBayesianKriging3D(in3DPoints1, elevationField1, valueField1, outGALayer1)
arcpy.ga.EmpiricalBayesianKriging3D(in3DPoints2, elevationField2, valueField2, outGALayer2)


# Export predictions for first model and probability that second model exceeds 10
# Export to gridded 3D points

# Set up variables
in_3d_ga_layers = outGALayer1+";"+outGALayer2
out_ncdf = "C:/gapydata/outputs/outputNetCDF1.nc"
export_locations = "3D_GRIDDED_POINTS"
x_spacing = "50 Meters"
y_spacing = "50 Meters"
elev_spacing = "5 Meters"
custom_points = ""
out_vars = "myGALayer1 PREDICTION #;myGALayer2 PROBABILITY 10"

# Run tool.
arcpy.ga.GALayer3DToNetCDF(in_3d_ga_layers, out_ncdf, export_locations,
                           x_spacing, y_spacing, elev_spacing, custom_points, out_vars)


# Export standard errors for first model and 75th quantile for second model
# Export to custom 3D points

# Set up variables
in_3d_ga_layers = outGALayer1+";"+outGALayer2
out_ncdf = "C:/gapydata/outputs/outputNetCDF2.nc"
export_locations = "CUSTOM_3D_POINTS"
x_spacing = ""
y_spacing = ""
elev_spacing = ""
custom_points = "C:/gapydata/inputs.gdb/myCustom3DPoints"
out_vars = "myGALayer1 PREDICTION_STANDARD_ERROR #;myGALayer2 QUANTILE 0.75"

# Run tool.
arcpy.ga.GALayer3DToNetCDF(in_3d_ga_layers, out_ncdf, export_locations,
                           x_spacing, y_spacing, elev_spacing, custom_points, out_vars)

Umgebungen

Ausdehnung, Ausgabe-Koordinatensystem, Faktor für parallele Verarbeitung, Geographische Transformationen, Aktueller Workspace

Lizenzinformationen

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