Étiquette | Explication | Type de données |
Entités en entrée | Entités ponctuelles en entrée contenant le champ à interpoler. | Feature Layer |
Champ d’élévation | Champ des entités en entrée contenant la valeur d’élévation de chaque point en entrée. Si les valeurs d’élévation sont stockées sous forme d’attributs géométriques dans Shape.Z, il est recommandé d’utiliser ce champ. Si les valeurs d’élévation sont stockées dans un champ attributaire, elles doivent indiquer la distance par rapport au niveau moyen de la mer. Les valeurs positives indiquent une distance au-dessus du niveau de la mer et les valeurs négatives une distance en dessous du niveau de la mer. | Field |
Champ de valeur | Champ des entités en entrée contenant les valeurs mesurées qui seront interpolées. | Field |
Couche géostatistique en sortie | Couche géostatistique en sortie qui affichera le résultat de l’interpolation. | Geostatistical Layer |
Unités du champ d’élévation (Facultatif) | Unités du champ d’élévation. Si Shape.Z est fourni comme champ d’élévation, les unités concordent automatiquement avec les unités z du système de coordonnées verticales.
| String |
Champ d’erreur de mesure (Facultatif) | Champ des entités en entrée contenant les valeurs d’erreur de mesure de chaque point. La valeur doit correspondre à un écart type de la valeur mesurée de chaque point. Utilisez ce champ si les valeurs d’erreur de mesure ne sont pas identiques pour chaque point. Des données mesurées à l’aide de différents appareils sont fréquemment à l’origine d’erreurs de mesure non constantes. Un appareil peut être plus précis qu’un autre, ce qui implique qu’il génère des erreurs de mesure plus faibles. Par exemple, un thermomètre arrondit au degré le plus proche tandis qu’un autre arrondit au dixième de degré le plus proche. La variabilité des mesures est généralement fournie par le fabricant de l’appareil de mesure ou déterminée de manière empirique. Laissez ce paramètre vide en l’absence de valeurs d’erreur de mesure ou si ces dernières sont inconnues. | Field |
Type de modèle de semi-variogramme (Facultatif) | Modèle de semi-variogramme qui sera utilisé pour l’interpolation.
| String |
Type de transformation (Facultatif) | Type de transformation qui sera appliqué aux entités en entrée.
| String |
Taille des sous-ensembles (Facultatif) | Taille des sous-ensembles. Les données en entrée seront automatiquement divisées en sous-ensembles avant le traitement. Ce paramètre contrôle le nombre de points que contiendra chaque sous-ensemble. | Long |
Facteur de superposition de surface du modèle local (Facultatif) | Facteur représentant le degré de superposition entre les modèles locaux (également appelés sous-ensembles). Chaque point en entrée peut être placé dans plusieurs sous-ensembles et le facteur de superposition spécifie le nombre moyen de sous-ensembles dont fera partie chaque point. Une valeur élevée du facteur de superposition produit une surface en sortie plus lisse, mais augmente également le temps de traitement. Les valeurs doivent être comprises entre 1 et 5. La superposition réelle qui sera utilisée sera généralement supérieure à cette valeur et chaque sous-ensemble contiendra par conséquent le même nombre de points. | Double |
Nombre de semi-variogrammes simulés (Facultatif) | Nombre de semi-variogrammes simulés qui seront utilisés pour chaque modèle local. En augmentant le nombre de simulations, les calculs du modèle seront plus stables, mais le temps de calcul sera plus long. | Long |
Ordre de suppression des tendances (Facultatif) | Spécifie l’ordre de suppression des tendances dans la direction verticale. Pour la plupart des données 3D, les valeurs des points changent plus rapidement verticalement qu’horizontalement. La suppression des tendances dans la direction verticale permet de compenser ce comportement et de stabiliser les calculs.
| String |
Facteur d’inflation de l’élévation (Facultatif) | Valeur constante qui est multipliée par la valeur Champ d’élévation avant de procéder à l’établissement de sous-groupes et à l’estimation du modèle. Pour la plupart des données 3D, les valeurs des points changent plus rapidement verticalement qu’horizontalement. Ce facteur étire les localisations des points de telle sorte qu’une unité de distance verticale soit statistiquement équivalente à une unité de distance horizontale. Les localisations des points reviendront à leur position d’origine avant le renvoi du résultat de l’interpolation. Cette correction est nécessaire pour estimer avec précision le modèle de semi-variogramme ainsi que les voisins appropriés pour le paramètre Voisinage de recherche. Le facteur d’inflation de l’élévation est sans unité et fournira les mêmes résultats quelles que soient les unités de la coordonnée x, y ou z des points en entrée. Si aucune valeur n’est fournie pour ce paramètre, une valeur est calculée au moment de l’exécution à l’aide d’une estimation de probabilité maximale. La valeur apparaîtra sous forme de message de géotraitement. La valeur calculée au moment de l’exécution sera comprise entre 1 et 1 000. Vous pouvez toutefois fournir des valeurs comprises entre 0,01 et 1 000 000. Si la valeur calculée est égale à 1 ou 1 000, vous pouvez indiquer des valeurs hors de cette plage et choisir une valeur basée sur la validation croisée. | Double |
Voisinage de recherche (Facultatif) | Indique le nombre et l’orientation des voisins qui seront utilisés pour prévoir des valeurs aux nouvelles localisations. Standard 3D
| Geostatistical Search Neighborhood |
Élévation en sortie par défaut (Facultatif) | Élévation par défaut de la valeur du paramètre Couche géostatistique en sortie. La couche géostatistique s’affichera sous forme d’une surface horizontale à une élévation donnée et ce paramètre indique cette élévation. Lors de sa création, l’élévation de la couche géostatistique peut être modifiée à l’aide du curseur de plage. | Double |
Type de surface en sortie (Facultatif) | Type de surface pour stocker les résultats d’interpolation.
| String |
Valeur de quantile (Facultatif) | Valeur de quantile pour laquelle la couche en sortie sera générée. | Double |
Type de seuil de probabilité (Facultatif) | Indique si la probabilité de dépassement ou de non-dépassement du seuil spécifié par une valeur doit être calculée.
| String |
Seuil de probabilité (Facultatif) | Valeur du seuil de probabilité. Si aucune valeur n’est fournie, la médiane (50e quantile) des données en entrée sera utilisée. | Double |
Disponible avec une licence Geostatistical Analyst.
Synthèse
Interpole des points 3D avec la méthodologie de krigeage bayésien empirique. Tous les points doivent avoir des coordonnées x, y et z et une valeur mesurée à interpoler. La sortie est une couche géostatistique 3D qui se calcule et s’affiche sous forme de transect 2D à une élévation donnée. Vous pouvez modifier l’élévation de la couche à l’aide du curseur de plage, ce qui met à jour la couche afin de représenter les prévisions interpolées pour la nouvelle élévation.
Les applications possibles de l’interpolation 3D sont les suivantes :
- Les océanographes créent des cartes de l’oxygène dissous et de la salinité à différentes profondeurs océaniques.
- Les chercheurs en sciences atmosphériques modélisent la pollution et les gaz à effet de serre dans l’atmosphère.
- Les géologues prédisent des propriétés géologiques souterraines telles que des concentrations en minéraux et la porosité.
Illustration
Utilisation
Les entités en entrée peuvent être fournies comme suit :
- Entités ponctuelles 3D dont les élévations sont stockées sous forme d’attribut géométrique dans Shape.Z
- Entités ponctuelles 2D dont les élévations sont stockées dans un champ attributaire
Il est recommandé de fournir des entités ponctuelles 3D car toutes les unités et conversions d’unités peuvent s’effectuer automatiquement. Vous pouvez convertir des entités ponctuelles 2D avec un champ d’élévation en entités ponctuelles 3D à l’aide de l’outil Entité en 3D par attribut.
Les couches géostatistiques en 3D peuvent être affichées en tant que couches voxel à l’aide de l’outil Couche GA 3D vers NetCDF. Elles peuvent également prévoir les points cibles en 3D et être exportées sous forme de rasters et isolignes d’entité à n’importe quelle élévation. Plusieurs rasters à différentes élévations peuvent également être simultanément exportés et enregistrés sous forme de jeu de données raster multidimensionnel.
Toutes les entités en entrée doivent se trouver dans un système de coordonnées projetées. Si les points sont stockés dans un système de coordonnées géographiques avec des coordonnées de latitude et de longitude, ils doivent être projetés à l’aide de l’outil Projeter avant d’utiliser cet outil.
Un voisinage de recherche Standard 3D est utilisé pour calculer les prévisions. Toutes les distances qui sont utilisées pour trouver les voisins seront calculées dans le système de coordonnées étirées après l’application de la valeur du paramètre Facteur d’inflation de l’élévation. Pour plus d’informations, consultez Modifications horizontales et verticales des valeurs de données.
Paramètres
arcpy.ga.EmpiricalBayesianKriging3D(in_features, elevation_field, value_field, out_ga_layer, {elevation_units}, {measurement_error_field}, {semivariogram_model_type}, {transformation_type}, {subset_size}, {overlap_factor}, {number_simulations}, {trend_removal}, {elev_inflation_factor}, {search_neighborhood}, {output_elevation}, {output_type}, {quantile_value}, {threshold_type}, {probability_threshold})
Nom | Explication | Type de données |
in_features | Entités ponctuelles en entrée contenant le champ à interpoler. | Feature Layer |
elevation_field | Champ des entités en entrée contenant la valeur d’élévation de chaque point en entrée. Si les valeurs d’élévation sont stockées sous forme d’attributs géométriques dans Shape.Z, il est recommandé d’utiliser ce champ. Si les valeurs d’élévation sont stockées dans un champ attributaire, elles doivent indiquer la distance par rapport au niveau moyen de la mer. Les valeurs positives indiquent une distance au-dessus du niveau de la mer et les valeurs négatives une distance en dessous du niveau de la mer. | Field |
value_field | Champ des entités en entrée contenant les valeurs mesurées qui seront interpolées. | Field |
out_ga_layer | Couche géostatistique en sortie qui affichera le résultat de l’interpolation. | Geostatistical Layer |
elevation_units (Facultatif) | Unités du champ d’élévation. Si Shape.Z est fourni comme champ d’élévation, les unités concordent automatiquement avec les unités z du système de coordonnées verticales.
| String |
measurement_error_field (Facultatif) | Champ des entités en entrée contenant les valeurs d’erreur de mesure de chaque point. La valeur doit correspondre à un écart type de la valeur mesurée de chaque point. Utilisez ce champ si les valeurs d’erreur de mesure ne sont pas identiques pour chaque point. Des données mesurées à l’aide de différents appareils sont fréquemment à l’origine d’erreurs de mesure non constantes. Un appareil peut être plus précis qu’un autre, ce qui implique qu’il génère des erreurs de mesure plus faibles. Par exemple, un thermomètre arrondit au degré le plus proche tandis qu’un autre arrondit au dixième de degré le plus proche. La variabilité des mesures est généralement fournie par le fabricant de l’appareil de mesure ou déterminée de manière empirique. Laissez ce paramètre vide en l’absence de valeurs d’erreur de mesure ou si ces dernières sont inconnues. | Field |
semivariogram_model_type (Facultatif) | Modèle de semi-variogramme qui sera utilisé pour l’interpolation.
| String |
transformation_type (Facultatif) | Type de transformation qui sera appliqué aux entités en entrée.
| String |
subset_size (Facultatif) | Taille des sous-ensembles. Les données en entrée seront automatiquement divisées en sous-ensembles avant le traitement. Ce paramètre contrôle le nombre de points que contiendra chaque sous-ensemble. | Long |
overlap_factor (Facultatif) | Facteur représentant le degré de superposition entre les modèles locaux (également appelés sous-ensembles). Chaque point en entrée peut être placé dans plusieurs sous-ensembles et le facteur de superposition spécifie le nombre moyen de sous-ensembles dont fera partie chaque point. Une valeur élevée du facteur de superposition produit une surface en sortie plus lisse, mais augmente également le temps de traitement. Les valeurs doivent être comprises entre 1 et 5. La superposition réelle qui sera utilisée sera généralement supérieure à cette valeur et chaque sous-ensemble contiendra par conséquent le même nombre de points. | Double |
number_simulations (Facultatif) | Nombre de semi-variogrammes simulés qui seront utilisés pour chaque modèle local. En augmentant le nombre de simulations, les calculs du modèle seront plus stables, mais le temps de calcul sera plus long. | Long |
trend_removal (Facultatif) | Spécifie l’ordre de suppression des tendances dans la direction verticale. Pour la plupart des données 3D, les valeurs des points changent plus rapidement verticalement qu’horizontalement. La suppression des tendances dans la direction verticale permet de compenser ce comportement et de stabiliser les calculs.
| String |
elev_inflation_factor (Facultatif) | Valeur constante qui est multipliée par la valeur Champ d’élévation avant de procéder à l’établissement de sous-groupes et à l’estimation du modèle. Pour la plupart des données 3D, les valeurs des points changent plus rapidement verticalement qu’horizontalement. Ce facteur étire les localisations des points de telle sorte qu’une unité de distance verticale soit statistiquement équivalente à une unité de distance horizontale. Les localisations des points reviendront à leur position d’origine avant le renvoi du résultat de l’interpolation. Cette correction est nécessaire pour estimer avec précision le modèle de semi-variogramme ainsi que les voisins appropriés pour le paramètre Voisinage de recherche. Le facteur d’inflation de l’élévation est sans unité et fournira les mêmes résultats quelles que soient les unités de la coordonnée x, y ou z des points en entrée. Si aucune valeur n’est fournie pour ce paramètre, une valeur est calculée au moment de l’exécution à l’aide d’une estimation de probabilité maximale. La valeur apparaîtra sous forme de message de géotraitement. La valeur calculée au moment de l’exécution sera comprise entre 1 et 1 000. Vous pouvez toutefois fournir des valeurs comprises entre 0,01 et 1 000 000. Si la valeur calculée est égale à 1 ou 1 000, vous pouvez indiquer des valeurs hors de cette plage et choisir une valeur basée sur la validation croisée. | Double |
search_neighborhood (Facultatif) | Indique le nombre et l’orientation des voisins à l’aide de la classe SearchNeighborhoodStandard3D. Standard 3D
| Geostatistical Search Neighborhood |
output_elevation (Facultatif) | Élévation par défaut de la valeur du paramètre out_ga_layer. La couche géostatistique s’affichera sous forme d’une surface horizontale à une élévation donnée et ce paramètre indique cette élévation. Lors de sa création, l’élévation de la couche géostatistique peut être modifiée à l’aide du curseur de plage. | Double |
output_type (Facultatif) | Type de surface pour stocker les résultats d’interpolation. Pour en savoir plus sur les types de surface en sortie, reportez-vous à Quels types de surface en sortie peuvent générer les modèles d’interpolation ?
| String |
quantile_value (Facultatif) | Valeur de quantile pour laquelle la couche en sortie sera générée. | Double |
threshold_type (Facultatif) | Indique si la probabilité de dépassement ou de non-dépassement du seuil spécifié par une valeur doit être calculée.
| String |
probability_threshold (Facultatif) | Valeur du seuil de probabilité. Si aucune valeur n’est fournie, la médiane (50e quantile) des données en entrée sera utilisée. | Double |
Exemple de code
Interpoler une classe d’entités ponctuelles 3D à l’aide de la fonction EmpiricalBayesianKriging3D.
import arcpy
arcpy.ga.EmpiricalBayesianKriging3D("my3DLayer", "Shape.Z", "myValueField", "myGALayer", "METER", "",
"POWER", "NONE", 100, 1, 100, "NONE", "",
"NBRTYPE=Standard3D RADIUS=10000 NBR_MAX=15 NBR_MIN=10 SECTOR_TYPE=ONE_SECTOR",
"", "PREDICTION", 0.5, "EXCEED", None)
Interpoler une classe d’entités ponctuelles 3D à l’aide de la fonction EmpiricalBayesianKriging3D.
# Name: EBK3D_Example_02.py
# Description: Interpolates 3D points.
# Requirements: Geostatistical Analyst Extension
# Author: Esri
# Import system modules
import arcpy
# Set local variables
in3DPoints = "C:/gapyexamples/input/my3DPoints.shp"
elevationField = "Shape.Z"
valueField = "myValueField"
outGALayer = "myGALayer"
elevationUnit = "METER"
measurementErrorField = "myMEField"
semivariogramModel = "LINEAR"
transformationType = "NONE"
subsetSize = 80
overlapFactor = 1.5
numSimulations = 200
trendRemoval = "FIRST"
elevInflationFactor = 20
radius = 10000
maxNeighbors = 15
minNeighbors = 10
sectorType = "FOUR_SECTORS"
searchNeighborhood = arcpy.SearchNeighborhoodStandard3D(radius, maxNeighbors, minNeighbors, sectorType)
outputElev = 1000
outputType = "PREDICTION"
# Check out the ArcGIS Geostatistical Analyst extension license
arcpy.CheckOutExtension("GeoStats")
# Execute Empirical Bayesian Kriging 3D
arcpy.ga.EmpiricalBayesianKriging3D(in3DPoints, elevationField, valueField, outGALayer, elevationUnit, myMEField,
semivariogramModel, transformationType, subsetSize, overlapFactor, numSimulations,
trendRemoval, elevInflationFactor, searchNeighborhood, outputElev, outputType)
Environnements
Informations de licence
- Basic: Nécessite Geostatistical Analyst
- Standard: Nécessite Geostatistical Analyst
- Advanced: Nécessite Geostatistical Analyst
Rubriques connexes
Vous avez un commentaire à formuler concernant cette rubrique ?