Sensibilité du semi-variogramme (Geostatistical Analyst)

Disponible avec une licence Geostatistical Analyst.

Synthèse

Cet outil réalise une analyse de la sensibilité sur les valeurs prévues et les erreurs standards associées en changeant les paramètres de semi-variogramme du modèle (la pépite, le seuil partiel et les portées majeure/mineure) selon un pourcentage des valeurs d’origine.

L’outil fait appel à une source de modèle géostatistique pour renseigner les valeurs initiales de la pépite, du seuil partiel et des portées majeure/mineure. La sortie de l’outil est une table indiquant les valeurs de paramètre qui ont été utilisées ainsi que les valeurs prévues et les valeurs d’erreur standard obtenues. Si la sortie affiche de grandes fluctuations alors que de petites modifications ont été apportées aux valeurs des paramètres du modèle, la sortie n’est pas vraiment fiable. En revanche, si les modifications de la sortie sont limitées, vous pouvez faire confiance aux prévisions du modèle et prendre des décisions en vous appuyant sur sa sortie.

Utilisation

  • La source de modèle géostatistique est soit une couche géostatistique, soit un modèle géostatistique (XML).

  • Dans la table en sortie, PtID correspond à la valeur OID dans la classe d’entités Localisations d’observation des points en entrée et le champ RndParam identifie les paramètres du semi-variogramme qui ont été modifiés de façon aléatoire.

  • Cet outil utilise un générateur de nombres aléatoires dans son fonctionnement. La valeur initiale utilisée peut être contrôlée dans l’environnement Générateur de nombres aléatoires.

    • Si une valeur initiale égale à 0 est utilisée (valeur par défaut), chaque fois que l’outil est exécuté, un jeu différent de nombres aléatoires est utilisé et une sortie différente est générée.
    • Si la valeur initiale aléatoire est définie sur un nombre fixe supérieur à 0, l’outil produira la même sortie à chaque exécution tant que la valeur initiale n’est pas modifiée.

    Remarque :

    Seul le type de générateur de nombres aléatoires Mersenne Twister est pris en charge. Si Algorithme de collecte ACM 599 ou Standard C Rand est choisi, Mersenne Twister est utilisé à la place.

  • Dans les scripts Python, la classe ArcPy GeostatisticalDatasets est utile pour renseigner le paramètre Jeu(x) de données en entrée.

  • Pour les formats de données prenant en charge les valeurs nulles (par exemple, les classes d’entités de géodatabase fichier), une valeur nulle est ignorée pour indiquer qu’une prédiction ne peut pas être établie pour cette localisation et que la valeur doit être ignorée lorsqu’elle est utilisée comme entrée. Pour les formats de données qui ne prennent pas en charge les valeurs nulles, par exemple les shapefiles, la valeur -1.7976931348623158e+308 est utilisée (c’est la valeur négative de la constante définie C++ DBL_MAX) pour indiquer qu’une prévision ne peut pas être établie pour cette localisation.

Paramètres

ÉtiquetteExplicationType de données
Source de modèle géostatistique en entrée

Source de modèle géostatistique à analyser.

File; Geostatistical Layer
Jeu(x) de données en entrée

Nom des jeux de données en entrée et noms de champ utilisés lors de la création de la couche en sortie.

Geostatistical Value Table
Localisations d’observation des points en entrée

Localisations ponctuelles auxquelles l’analyse de la sensibilité est effectuée.

Feature Layer
Plage de la pépite (% de la valeur du modèle)
(Facultatif)

Pourcentage soustrait et ajouté au paramètre Pépite pour créer une plage en vue de la sélection ultérieure d’un paramètre Pépite aléatoire.

Double
Nombre de calculs pour la pépite
(Facultatif)

Nombre de valeurs de pépite échantillonnées de façon aléatoire à partir de la plage de la pépite.

Long
Plage du seuil partiel (% de la valeur du modèle)
(Facultatif)

Pourcentage soustrait et ajouté au paramètre Seuil partiel pour créer une plage en vue de la sélection d’un seuil partiel aléatoire.

Double
Nombre de calculs pour le seuil partiel
(Facultatif)

Nombre de valeurs de seuil partiel échantillonnées de façon aléatoire à partir de la plage du seuil partiel.

Long
Plage de la portée majeure (% de la valeur du modèle)
(Facultatif)

Pourcentage soustrait et ajouté au paramètre Portée majeure pour créer une plage en vue de la sélection d’une portée majeure aléatoire.

Double
Nombre de calculs pour la portée majeure
(Facultatif)

Nombre de valeurs de portée majeure échantillonnées de façon aléatoire à partir de la plage de la portée majeure.

Long
Plage de la portée mineure (% de la valeur du modèle)
(Facultatif)

Pourcentage soustrait et ajouté au paramètre Portée mineure pour créer une plage en vue de la sélection d’une portée mineure aléatoire.

Double
Nombre de calculs pour la portée mineure
(Facultatif)

Nombre de valeurs de portée mineure échantillonnées de façon aléatoire à partir de la plage de la portée mineure.

Si le paramètre Anisotropie a été défini dans la source de modèle géostatistique en entrée, une valeur est obligatoire.

Long
Table en sortie

Table stockant les résultats de la sensibilité.

Table

arcpy.ga.GASemivariogramSensitivity(in_ga_model_source, in_datasets, in_locations, {nugget_span_percents}, {nugget_calc_times}, {partialsill_span_percents}, {partialsill_calc_times}, {range_span_percents}, {range_calc_times}, {minrange_span_percents}, {minrange_calc_times}, out_table)
NomExplicationType de données
in_ga_model_source

Source de modèle géostatistique à analyser.

File; Geostatistical Layer
in_datasets

Un objet GeostatisticalDatasets.

Par ailleurs, il peut s’agir d’une chaîne d’éléments séparés par des points virgules. Chaque élément est composé des composants suivants :

  • Chemin d’accès au catalogue et nom attribués au jeu de données, ou nom d’une couche dans la table des matières actuelle, suivi d’un espace.
  • Suite de noms de champ séparés par un espace. Dans le cas d’un raster, les valeurs de cellule sont utilisées.
Geostatistical Value Table
in_locations

Localisations ponctuelles auxquelles l’analyse de la sensibilité est effectuée.

Feature Layer
nugget_span_percents
(Facultatif)

Pourcentage soustrait et ajouté au paramètre Pépite pour créer une plage en vue de la sélection ultérieure d’un paramètre Pépite aléatoire.

Double
nugget_calc_times
(Facultatif)

Nombre de valeurs de pépite échantillonnées de façon aléatoire à partir de la plage de la pépite.

Long
partialsill_span_percents
(Facultatif)

Pourcentage soustrait et ajouté au paramètre Seuil partiel pour créer une plage en vue de la sélection d’un seuil partiel aléatoire.

Double
partialsill_calc_times
(Facultatif)

Nombre de valeurs de seuil partiel échantillonnées de façon aléatoire à partir de la plage du seuil partiel.

Long
range_span_percents
(Facultatif)

Pourcentage soustrait et ajouté au paramètre Portée majeure pour créer une plage en vue de la sélection d’une portée majeure aléatoire.

Double
range_calc_times
(Facultatif)

Nombre de valeurs de portée majeure échantillonnées de façon aléatoire à partir de la plage de la portée majeure.

Long
minrange_span_percents
(Facultatif)

Pourcentage soustrait et ajouté au paramètre Portée mineure pour créer une plage en vue de la sélection d’une portée mineure aléatoire.

Double
minrange_calc_times
(Facultatif)

Nombre de valeurs de portée mineure échantillonnées de façon aléatoire à partir de la plage de la portée mineure.

Si le paramètre Anisotropie a été défini dans la source de modèle géostatistique en entrée, une valeur est obligatoire.

Long
out_table

Table stockant les résultats de la sensibilité.

Table

Exemple de code

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

Effectuer une analyse de la sensibilité en faisant varier les valeurs Pépite, Seuil partiel et Portée.

import arcpy
arcpy.env.workspace = "C:/gapyexamples/data"
arcpy.GASemivariogramSensitivity_ga(
    "C:/gapyexamples/data/kriging.lyr", "C:/gapyexamples/data/ca_ozone_pts.shp OZONE",
    "C:/gapyexamples/data/obs_pts.shp", "", "", "", "", "", "", "", "", 
    "C:/gapyexamples/output/outtabSS")
Exemple 2 d’utilisation de la fonction SemivariogramSensitivity (script autonome)

Effectuer une analyse de la sensibilité en faisant varier les valeurs Pépite, Seuil partiel et Portée.

# Name: SemivariogramSensitivity_Example_02.py
# Description: The semivariogram parameters Nugget, Partial Sill and Range can
#   be varied to perform a sensitivity analysis.
# Requirements: Geostatistical Analyst Extension

# Import system modules
import arcpy

# Set environment settings
arcpy.env.workspace = "C:/gapyexamples/data"

# Set local variables
inLayer = "C:/gapyexamples/data/kriging.lyr"
inData = "C:/gapyexamples/data/ca_ozone_pts.shp OZONE"
inObs = "C:/gapyexamples/data/obs_pts.shp"
nugPercents = ""
nugCalc = ""
sillPercents = ""
sillCalc = ""
rangePercents = ""
rangeClac = ""
minrangePercent = ""
midrangeCalc = ""
outTable = "C:/gapyexamples/output/outtabSS"

# Execute SemivariogramSensitivity
arcpy.GASemivariogramSensitivity_ga(inLayer, inData, inObs, nugPercents, nugCalc,
                                    sillPercents, sillCalc, rangePercents, rangeClac,
                                    minrangePercent, midrangeCalc, outTable)

Informations de licence

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

Rubriques connexes