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 fichier en sortie s’affiche en tant que couche voxel dans une scène locale.

En savoir plus sur les couches voxel

Illustration

La couche voxel affiche les résultats d’une interpolation 3D.

Utilisation

• Si l’outil s’exécute dans une scène locale avec les mêmes systèmes de coordonnées horizontales et verticales que les entités en entrée, une couche voxel est ajoutée à la scène pour vous permettre d’explorer interactivement les résultats. Vous pouvez également ajouter le fichier netCDF en sortie en tant que couche voxel à l’aide de l’outil Générer une couche voxel multidimensionnelle ou de la boîte de dialogue Ajouter une couche voxel multidimensionnelle.

  Vous pouvez convertir le fichier netCDF en sortie en raster multidimensionnel à l’aide de l’outil Copier un raster. Vous pouvez également l’ajouter à une carte 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.

• Toutes les couches géostatistiques 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 fournissez des entités ponctuelles 3D dans le paramètre Localisations ponctuelles 3D pour spécifier les localisations à prévoir.

    Conseil :

    Vous pouvez utiliser l’outil Entité en 3D par attribut pour créer des points 3D à partir de points 2D avec un champ d’élévation.

• L’étendue en sortie par défaut du fichier netCDF en sortie est l’union des étendues de toutes les couches géostatistiques 3D en entrée.

• Vous pouvez utiliser les paramètres Polygones de zone d’étude en entrée, Raster de découpage de l’élévation minimale et Raster de découpage de l’élévation maximale pour limiter l’analyse à une zone d’étude spécifique et entre deux surfaces d’élévation. Les voxels hors de ces limites n’auront pas de valeur et ne seront pas affichés. Par exemple, si les points se trouvent dans une réserve marine, vous pouvez créer une couche voxel qui ne s’affiche que dans un polygone de la réserve (zone d’étude), au-dessus du plancher océanique (raster d’élévation minimale) et sous la thermocline (raster d’élévation maximale).

  Divers aspects doivent être pris en compte lors de l’utilisation de surfaces d’élévation comme rasters d’élévation minimale et maximale. Les services d’imagerie, les couches d’élévation Web et les couches d’imagerie Web présentent les performances les plus faibles et des erreurs peuvent se produire lorsque le nombre de requêtes est élevé. Les rasters enregistrés en tant que fichiers locaux sur le disque offrent de meilleures performances et sont recommandés lors de la création de couches voxel haute résolution couvrant de grandes étendues spatiales.

• Lors de la mise en page de la grille 3D de points qui représentera les voxels, le premier point est créé avec les coordonnées x, y et z minimales de l’étendue en sortie (par défaut, l’étendue des entités en entrée). Les points restants sont créés en itérant les distances des paramètres Espacement sur l’axe x, Espacement sur l’axe y et Espacement sur l’axe d’élévation sur les dimensions de l’étendue en sortie. Si certaines 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 pour x est spécifiée comme comprise entre 0 et 10 et que le paramètre Espacement sur l’axe x a la valeur 3, la sortie comportera 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.

• Les rasters de découpage de l’élévation minimale et maximale doivent se trouver dans un système de coordonnées projetées. Si les points ou les rasters possèdent un système de coordonnées géographiques avec des coordonnées de latitude et de longitude, ils doivent être projetés dans un système de coordonnées projetées à l’aide de l’outil Projeter ou Projeter un raster.

• 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 sur l’axe x et Espacement sur l’axe y 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

Boîte de dialogue

Étiquette	Explication	Type de données
Couches géostatistiques 3D en entrée	Couches géostatistiques 3D qui sont exportées vers un fichier netCDF en sortie.	Geostatistical Layer
Fichier NetCDF en sortie	Fichier netCDF en sortie contenant les valeurs exportées des couches géostatistiques en entrée. Les résultats de chaque couche géostatistique sont enregistrés sous différentes variables dans le fichier netCDF.	File
Exporter les localisations (Facultatif)	Indique les localisations à exporter depuis la valeur du paramètre Couches géostatistiques 3D en entrée. Si vous choisissez l’option 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 l’option Points 3D personnalisés, vous devez fournir des entités ponctuelles 3D dans le paramètre Localisations ponctuelles 3D représentant les localisations à exporter. Maillages de points en 3D—Les localisations de prédiction sont des maillages de points en 3D. Il s’agit de l’option par défaut. Points 3D personnalisés—Les 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 des valeurs 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
Raster de découpage de l’élévation minimale (Facultatif)	Raster d’élévation qui sera utilisé pour découper le bas de la couche voxel. Seuls les voxels se trouvant au-dessus de ce raster d’élévation se verront affecter des prévisions. Par exemple, si vous utilisez un raster d’élévation du sol, la couche voxel ne sera affichée qu’au-dessus de la surface. Ce raster peut aussi être utilisé pour les surfaces de substrat rocheux ou le fond d’un dépôt de schiste. Le raster doit se trouver dans un système de coordonnées projetées et l’unité des valeurs d’élévation doit être identique à l’unité verticale du raster.	Raster Layer
Raster de découpage de l’élévation maximale (Facultatif)	Raster d’élévation qui sera utilisé pour découper le haut de la couche voxel. Seuls les voxels se trouvant sous ce raster d’élévation se verront affecter des prévisions. Par exemple, si vous utilisez un raster d’élévation du sol, la couche voxel ne sera affichée que sous la surface. Ce raster peut aussi être utilisé pour découper des voxels sur la partie supérieure d’un espace aérien restreint. Le raster doit se trouver dans un système de coordonnées projetées et l’unité des valeurs d’élévation doit être identique à l’unité verticale du raster.	Raster Layer

Sortie obtenue

Étiquette	Explication	Type de données
Couche voxel en sortie	Couche voxel des valeurs prévues.	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})

Nom	Explication	Type de données
in_3d_geostat_layers [in_3d_geostat_layers,...]	Couches géostatistiques 3D qui sont exportées vers un fichier netCDF en sortie.	Geostatistical Layer
out_netcdf_file	Fichier netCDF en sortie contenant les valeurs exportées des couches géostatistiques en entrée. Les résultats de chaque couche géostatistique sont enregistrés sous différentes variables dans le fichier netCDF.	File
export_locations (Facultatif)	Indique les localisations à exporter depuis la valeur du paramètre in_3d_geostat_layers. Si vous choisissez l’option 3D_GRIDDED_POINTS, vous devez fournir des valeurs pour les paramètres x_spacing, u_spacing et elevation_spacing qui représentent la distance entre chaque point maillé dans toutes les dimensions. Si vous choisissez l’option CUSTOM_3D_POINTS, vous devez fournir des entités ponctuelles 3D dans le paramètre in_points_3d représentant les localisations à exporter. 3D_GRIDDED_POINTS—Les localisations de prédiction sont des maillages de points en 3D. Il s’agit de l’option par défaut. CUSTOM_3D_POINTS—Les 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 des valeurs 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
min_elev_raster (Facultatif)	Raster d’élévation qui sera utilisé pour découper le bas de la couche voxel. Seuls les voxels se trouvant au-dessus de ce raster d’élévation se verront affecter des prévisions. Par exemple, si vous utilisez un raster d’élévation du sol, la couche voxel ne sera affichée qu’au-dessus de la surface. Ce raster peut aussi être utilisé pour les surfaces de substrat rocheux ou le fond d’un dépôt de schiste. Le raster doit se trouver dans un système de coordonnées projetées et l’unité des valeurs d’élévation doit être identique à l’unité verticale du raster.	Raster Layer
max_elev_raster (Facultatif)	Raster d’élévation qui sera utilisé pour découper le haut de la couche voxel. Seuls les voxels se trouvant sous ce raster d’élévation se verront affecter des prévisions. Par exemple, si vous utilisez un raster d’élévation du sol, la couche voxel ne sera affichée que sous la surface. Ce raster peut aussi être utilisé pour découper des voxels sur la partie supérieure d’un espace aérien restreint. Le raster doit se trouver dans un système de coordonnées projetées et l’unité des valeurs d’élévation doit être identique à l’unité verticale du raster.	Raster Layer

Sortie obtenue

Nom	Explication	Type de données
out_voxel_layer	Couche voxel des valeurs prévues.	Voxel 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)

Environnements

Étendue, Système de coordonnées en sortie, Facteur de traitement parallèle, Transformations géographiques, Espace de travail courant

Informations de licence

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