Couche GA 3D vers NetCDF (Geostatistical Analyst)

Disponible avec une licence Geostatistical Analyst.

Synthèse

Exporte une ou plusieurs couches géostatistiques 3D créées à l’aide de l’outil Krigeage bayésien empirique 3D au format netCDF (fichier *.nc). Le principal objectif de cet outil est de préparer les couches géostatistiques 3D pour les visualiser en tant que couche voxel dans une scène locale.

En savoir plus sur les couches voxel

Illustration

Couche géostatistique 3D exportée vers une couche voxel
La couche voxel affiche les résultats d’une interpolation 3D.

Utilisation

  • Le paramètre Fichier NetCDF en sortie n’ajoute pas de couches automatiquement sur la carte ou scène. Pour afficher la sortie sous forme de couche voxel, elle doit être ajoutée à une scène locale avec le même système de coordonnées horizontales et verticales que le fichier netCDF. Elle peut également être ajoutée à une carte ou scène en tant que couche d’entités ou raster à l’aide de l’outil Générer une couche d’entités NetCDF ou Générer une couche raster NetCDF, respectivement.

  • Pour afficher le fichier netCDF en sortie en tant que couche voxel, définissez le paramètre Exporter les localisations sur Maillages de points en 3D. La couche voxel requiert des localisations maillées pour s’afficher.

  • Toutes les valeurs Couches géostatistiques 3D en entrée doivent être des couches géostatistiques 3D créées avec l’outil Krigeage bayésien empirique 3D.

  • Vous pouvez choisir les localisations à prédire et exporter de l’une des façons suivantes :

    • Prédire des maillages de points en 3D : les prédictions seront exportées vers des maillages de points en 3D. Pour utiliser cette option, spécifiez Maillages de points en 3D pour le paramètre Exporter les localisations. Les paramètres Espacement X, Espacement Y et Espacement d’élévation indiquent la distance entre chaque point maillé dans les dimensions x, y et z (élévation).
    • Prédire les localisations 3D personnalisées : les prédictions seront exportées vers des localisations 3D personnalisées définies par des entités ponctuelles 3D. Pour utiliser cette option, spécifiez Points 3D personnalisés pour le paramètre Exporter les localisations et ajoutez des entités ponctuelles 3D au paramètre Localisations ponctuelles 3D pour spécifier les localisations à exporter. Les élévations des localisations ponctuelles 3D doivent être stockées dans l’attribut géométrique Shape.Z. Si les élévations sont stockées sous forme de champ attributaire et non d’attribut Shape.Z, vous pouvez les convertir en points 3D à l’aide de l’outil Entité en 3D par attribut.

  • L’étendue en sortie par défaut du paramètre Fichier netCDF en sortie est l’union des étendues de toutes les Couches géostatistiques 3D en entrée.

  • Si le paramètre Maillages de points en 3D est spécifié pour le paramètre Exporter les localisations, le premier point est créé à la coordonnée x minimale, y minimale et z minimale de l’environnement Étendue en sortie. Les points restants sont créés en itérant les distances des paramètres Espacement X, Espacement Y et Espacement d’élévation à travers les dimensions de l’étendue en sortie. Si certaines des distances d’espacement ne divisent pas équitablement la dimension correspondante de l’étendue en sortie, une ligne ou colonne de points est créée au-delà de l’étendue en sortie. Par exemple, si l’étendue en sortie de x est spécifiée comme comprise entre 0 et 10 et que le paramètre Espacement X a la valeur 3, la sortie aura 4 lignes dans l’étendue x : 0, 3, 6, 9 et 12. De même, une ligne ou colonne supplémentaire est créée si les distances d’espacement ne divisent pas équitablement les étendues y ou z.

  • Le paramètre Polygones de zone d’étude en entrée permet de limiter les résultats d’interpolation à une zone d’étude spécifique au lieu d’afficher un cube entier de résultats. Des points seront créés dans le fichier netCDF en sortie uniquement pour les localisations situées dans la zone d’étude.

  • Si des polygones de zone d’étude en entrée sont fournis, l’étendue de la zone d’étude est utilisée comme étendue en sortie par défaut, et les valeurs des paramètres Espacement X, Espacement Y et Espacement d’élévation sont recalculées en fonction de cette étendue. Cela permet de s’assurer que la sortie occupera entièrement la zone d’étude par défaut.

Paramètres

ÉtiquetteExplicationType de données
Couches géostatistiques 3D en entrée

Couches géostatistiques 3D qui seront exportées vers le Fichier netCDF en sortie. Si plusieurs couches sont fournies, la sortie est un fichier netCDF multivarié.

Geostatistical Layer
Fichier NetCDF en sortie

Fichier netCDF en sortie contenant les valeurs exportées du paramètre Couches géostatistiques 3D en entrée.

File
Exporter les localisations
(Facultatif)

Indique les localisations à exporter du paramètre Couches géostatistiques 3D en entrée. Vous pouvez exporter des maillages de points en 3D ou fournir des entités ponctuelles 3D personnalisées pour représenter les localisations exportées. Si vous choisissez Maillages de points en 3D, vous devez fournir des valeurs pour les paramètres Espacement X, Espacement Y et Espacement d’élévation qui représentent la distance entre chaque point maillé dans toutes les dimensions. Si vous choisissez Points 3D personnalisés, vous devez fournir des entités ponctuelles dans le paramètre Localisations ponctuelles 3D représentant les localisations à exporter.

  • Maillages de points en 3DLes localisations de prédiction sont des maillages de points en 3D. Il s’agit de l’option par défaut.
  • Points 3D personnalisésLes localisations de prédiction sont définies par des entités ponctuelles 3D personnalisées.
String
Espacement X
(Facultatif)

Espacement entre chaque point maillé dans la dimension x. La valeur par défaut crée 40 points le long de l’étendue x en sortie.

Linear Unit
Espacement Y
(Facultatif)

Espacement entre chaque point maillé dans la dimension y. La valeur par défaut crée 40 points le long de l’étendue y en sortie.

Linear Unit
Espacement d’élévation
(Facultatif)

Espacement entre chaque point maillé dans la dimension d’élévation (z). La valeur par défaut crée 40 points le long de l’étendue z en sortie.

Linear Unit
Localisations ponctuelles 3D
(Facultatif)

Entités ponctuelles 3D représentant les localisations à exporter. Les élévations des entités ponctuelles doivent être stockées dans l’attribut géométrique Shape.Z.

Feature Layer
Variables en sortie
(Facultatif)

Indique les types en sortie du paramètre Couches géostatistiques 3D en entrée. Vous pouvez spécifier un ou plusieurs types en sortie pour chacune des couches ou appliquer un type en sortie à toutes les couches géostatistiques en entrée. Par défaut, les prédictions de toutes les couches sont exportées

Pour exporter d’autres types de sortie, spécifiez la couche à exporter (ou choisissez Tous pour indiquer toutes les couches) dans la première entrée de la table de valeurs. Spécifiez le type de sortie dans la deuxième entrée de la table de valeurs. Si vous avez choisi Probabilité ou Quantile comme type de sortie, spécifiez la valeur de seuil (pour la probabilité) ou la valeur de quantile (pour le quantile) dans la troisième entrée de la table de valeurs. Si vous avez choisi Prédiction ou Erreur standard de prédiction comme type de sortie, vous pouvez laisser la troisième entrée vide dans la table de valeurs.

En savoir plus sur les types de sorties géostatistiques.

Value Table
Polygones de zone d’étude en entrée
(Facultatif)

Entités surfaciques qui représentent la zone d’étude. Seuls les points qui se trouvent dans la zone d’étude sont enregistrés dans le fichier netCDF en sortie. En cas d’affichage sous forme de couche voxel, seuls les voxels dans la zone d’étude s’affichent dans la scène. Les points sont déterminés comme étant à l’intérieur ou à l’extérieur de la zone d ’étude avec uniquement leurs coordonnées x et y.

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})
NomExplicationType de données
in_3d_geostat_layers
[in_3d_geostat_layers,...]

Couches géostatistiques 3D qui seront exportées vers le Fichier netCDF en sortie. Si plusieurs couches sont fournies, la sortie est un fichier netCDF multivarié.

Geostatistical Layer
out_netcdf_file

Fichier netCDF en sortie contenant les valeurs exportées du paramètre Couches géostatistiques 3D en entrée.

File
export_locations
(Facultatif)

Indique les localisations à exporter du paramètre Couches géostatistiques 3D en entrée. Vous pouvez exporter des maillages de points en 3D ou fournir des entités ponctuelles 3D personnalisées pour représenter les localisations exportées. Si vous choisissez Maillages de points en 3D, vous devez fournir des valeurs pour les paramètres Espacement X, Espacement Y et Espacement d’élévation qui représentent la distance entre chaque point maillé dans toutes les dimensions. Si vous choisissez Points 3D personnalisés, vous devez fournir des entités ponctuelles dans le paramètre Localisations ponctuelles 3D représentant les localisations à exporter.

  • 3D_GRIDDED_POINTSLes localisations de prédiction sont des maillages de points en 3D. Il s’agit de l’option par défaut.
  • CUSTOM_3D_POINTSLes localisations de prédiction sont définies par des entités ponctuelles 3D personnalisées.
String
x_spacing
(Facultatif)

Espacement entre chaque point maillé dans la dimension x. La valeur par défaut crée 40 points le long de l’étendue x en sortie.

Linear Unit
y_spacing
(Facultatif)

Espacement entre chaque point maillé dans la dimension y. La valeur par défaut crée 40 points le long de l’étendue y en sortie.

Linear Unit
elevation_spacing
(Facultatif)

Espacement entre chaque point maillé dans la dimension d’élévation (z). La valeur par défaut crée 40 points le long de l’étendue z en sortie.

Linear Unit
in_points_3d
(Facultatif)

Entités ponctuelles 3D représentant les localisations à exporter. Les élévations des entités ponctuelles doivent être stockées dans l’attribut géométrique Shape.Z.

Feature Layer
output_variables
[[layer_name, output_type, quantile_probability_value],...]
(Facultatif)

Indique les types en sortie du paramètre Couches géostatistiques 3D en entrée. Vous pouvez spécifier un ou plusieurs types en sortie pour chacune des couches ou appliquer un type en sortie à toutes les couches géostatistiques en entrée. Par défaut, les prédictions de toutes les couches sont exportées

Pour exporter d’autres types de sortie, spécifiez la couche à exporter (ou choisissez Tous pour indiquer toutes les couches) dans la première entrée de la table de valeurs. Spécifiez le type de sortie dans la deuxième entrée de la table de valeurs. Si vous avez choisi Probabilité ou Quantile comme type de sortie, spécifiez la valeur de seuil (pour la probabilité) ou la valeur de quantile (pour le quantile) dans la troisième entrée de la table de valeurs. Si vous avez choisi Prédiction ou Erreur standard de prédiction comme type de sortie, vous pouvez laisser la troisième entrée vide dans la table de valeurs.

En savoir plus sur les types de sorties géostatistiques.

Value Table
in_study_area
(Facultatif)

Entités surfaciques qui représentent la zone d’étude. Seuls les points qui se trouvent dans la zone d’étude sont enregistrés dans le fichier netCDF en sortie. En cas d’affichage sous forme de couche voxel, seuls les voxels dans la zone d’étude s’affichent dans la scène. Les points sont déterminés comme étant à l’intérieur ou à l’extérieur de la zone d ’étude avec uniquement leurs coordonnées x et y.

Feature Layer

Exemple de code

Exemple 1 d’utilisation de la fonction GALayer3DToNetCDF (fenêtre Python)

Interpolez deux fois des points 3D et convertissez les sorties en fichier netCDF multivarié.

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 #")
Exemple 2 d’utilisation de la fonction GALayer3DToNetCDF (script autonome)

Interpolez deux fois des points 3D et convertissez les sorties en fichier netCDF multivarié. Exportez vers des maillages de points 3D personnalisés.

# 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)

Informations de licence

  • Basic: Nécessite Geostatistical Analyst
  • Standard: Nécessite Geostatistical Analyst
  • Advanced: Nécessite Geostatistical Analyst

Rubriques connexes