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). Hauptsächlich dient dieses Werkzeug zum Vorbereiten der geostatistischen 3D-Layer für die Visualisierung als Voxel-Layer in einer lokalen Szene.

Weitere Informationen zu Voxel-Layern

Abbildung

Geostatistischer 3D-Layer nach dem Export in einen Voxel-Layer
Der Voxel-Layer zeigt die Ergebnisse einer 3D-Interpolation an.

Verwendung

  • Der Parameter Ausgabe-netCDF-Datei fügt der Karte oder Szene nicht automatisch Layer hinzu. Wenn die Ausgabe als Voxel-Layer dargestellt werden soll, dann muss er mit dem horizontalen und vertikalen Koordinatensystem der netCDF-Datei einer lokalen Szene hinzugefügt werden. Er kann auch mit dem Werkzeug NetCDF-Feature-Layer erstellen oder mit dem Werkzeug NetCDF-Raster-Layer erstellen als Feature- bzw. Raster-Layer in einer Karte oder Szene hinzugefügt werden.

  • Wenn die Ausgabe-netCDF-Datei als Voxel-Layer ausgegeben werden soll, dann legen Sie den Parameter Export-Positionen auf 3D-Gitterpunkte fest. Der Voxel-Layer benötigt gerasterte Positionen, damit er dargestellt werden kann.

  • Alle geostatistischen 3D-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 zu exportierenden Positionen angeben, hinzu. Die Höhen der 3D-Punktpositionen müssen im Geometrieattribut in "Shape.X" gespeichert sein. Wenn die Höhen nicht in "Shape.X", sondern als Attributfeld gespeichert sind, dann können Sie sie mit dem Werkzeug Feature zu 3D nach Attribut in 3D-Punkte konvertieren.

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

  • Wenn 3D-Gitterpunkte für den Parameter Export-Positionen angegeben wurde, dann wird der erste Punkt an der minimalen X-, minimalen Y- und minimalen Z-Koordinate der Umgebung für die Ausgabeausdehnung erstellt. 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 4 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.

  • Der Parameter Eingabe-Untersuchungsgebiet-Polygone kann verwendet werden, um die Interpolationsergebnisse auf ein bestimmtes Untersuchungsgebiet zu beschränken, anstatt einen vollständigen Würfel mit allen Ergebnissen anzuzeigen. In der Ausgabe-netCDF-Datei werden Punkte nur für die Positionen innerhalb des Untersuchungsgebietes erstellt.

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

Parameter

BeschriftungErläuterungDatentyp
Geostatistische Eingabe-Layer in 3D

Die geostatistischen 3D-Layer, die in die Ausgabe-netCDF-Datei exportiert werden sollen. Wenn mehr als ein Layer angegeben wurde, ist die Ausgabe eine multivariate NetCDF-Datei.

Geostatistical Layer
Ausgabe-netCDF-Datei

Die Ausgabe-netCDF-Datei, die die exportierten Werte aus den geostatistischen Eingabe-Layern in 3D enthält.

File
Export-Positionen
(optional)

Gibt die Positionen an, die aus den geostatistischen Eingabe-Layern in 3D exportiert werden sollen. Sie können in 3D-Gitterpunkte exportieren oder benutzerdefinierte 3D-Punkt-Features zur Darstellung der Export-Positionen angeben. Wenn Sie Gitterpunkte in 3D auswählen, dann 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 Benutzerdefinierte 3D-Punkte auswählen, dann müssen Sie im Parameter 3D-Punktpositionen die 3D-Punkt-Features angeben, die die zu exportierenden Positionen darstellen.

  • 3D-GitterpunkteDie Vorhersage-Positionen sind 3D-Gitterpunkte. Dies ist die Standardeinstellung.
  • Benutzerdefinierte 3D-PunkteDie 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 geostatistischen 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

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})
NameErläuterungDatentyp
in_3d_geostat_layers
[in_3d_geostat_layers,...]

Die geostatistischen 3D-Layer, die in die Ausgabe-netCDF-Datei exportiert werden sollen. Wenn mehr als ein Layer angegeben wurde, ist die Ausgabe eine multivariate NetCDF-Datei.

Geostatistical Layer
out_netcdf_file

Die Ausgabe-netCDF-Datei, die die exportierten Werte aus den geostatistischen Eingabe-Layern in 3D enthält.

File
export_locations
(optional)

Gibt die Positionen an, die aus den geostatistischen Eingabe-Layern in 3D exportiert werden sollen. Sie können in 3D-Gitterpunkte exportieren oder benutzerdefinierte 3D-Punkt-Features zur Darstellung der Export-Positionen angeben. Wenn Sie Gitterpunkte in 3D auswählen, dann 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 Benutzerdefinierte 3D-Punkte auswählen, dann müssen Sie im Parameter 3D-Punktpositionen die 3D-Punkt-Features angeben, die die zu exportierenden Positionen darstellen.

  • 3D_GRIDDED_POINTSDie Vorhersage-Positionen sind 3D-Gitterpunkte. Dies ist die Standardeinstellung.
  • CUSTOM_3D_POINTSDie 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 geostatistischen 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

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)

Lizenzinformationen

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

Verwandte Themen