Krigeage (3D Analyst)

Niveau de licence :BasicStandardAdvanced

Récapitulatif

Interpole une surface raster à partir de points à l'aide du krigeage.

Pour en savoir plus sur la fonction Krigeage

Utilisation

Syntaxe

Kriging_3d (in_point_features, z_field, out_surface_raster, semiVariogram_props, {cell_size}, {search_radius}, {out_variance_prediction_raster})
ParamètreExplicationType de données
in_point_features

Les entités points en entrée contenant les valeurs z à interpoler dans un raster de surface.

Feature Layer
z_field

Champ contenant une valeur de hauteur ou de magnitude pour chaque point.

Il peut s'agir d'un champ numérique ou du champ Forme, si les entités ponctuelles en entrée contiennent des valeurs Z.

Field
out_surface_raster

Raster de surface interpolé en sortie.

Raster Dataset
semiVariogram_props
kriging_model

Modèle de semi-variogramme à utiliser.

Il existe deux modèles de krigeage : ordinaire et universel. Le modèle ordinaire présente cinq types de semi-variogramme disponibles. Le modèle universel présente deux types de semi-variogramme disponibles. Chaque semi-variogramme comporte plusieurs paramètres facultatifs qui peuvent également être définis.

  • Semi-variogramme ordinaire :
    • Spherical – Modèle de semi-variogramme sphérique. Il s'agit de l'option par défaut.
    • Circular – Modèle de semi-variogramme circulaire.
    • Exponential - Modèle de semi-variogramme exponentiel.
    • Gaussian - Modèle de semi-variogramme gaussien (ou de distribution normale).
    • Linear - Modèle de semi-variogramme linéaire avec un seuil.
  • Modèle de semi-variogramme universel :
    • LinearDrift - Krigeage universel avec dérive linéaire.
    • QuadraticDrift - Krigeage universel avec dérive quadratique.
  • Une fois le modèle de semi-variogramme défini, les paramètres restants sont les mêmes pour les krigeages ordinaire et universel. Il s'agit des suivants :
    • Lag size - La valeur par défaut est la taille de cellule du raster en sortie.
    • MajorRange - Représente une distance au-delà de laquelle il n'y a que peu ou pas de corrélation.
    • PartialSill - Différence entre la pépite et le seuil.
    • Nugget - Représente l'erreur et la variation à des échelles spatiales trop petites pour être détectées. L'effet pépite est considéré comme une discontinuité à l'origine.

Le semi-variogramme prend la forme d'une chaîne de texte :

"{semivariogramType},{lagSize},{majorRange},{partialSill},{nugget}"

Par exemple :

"Circular, 2000, 2,6, 542"
KrigingModel
cell_size
(Facultatif)

Taille des cellules qui sera utilisée pour la création du raster en sortie.

Il s'agit de la valeur dans l'environnement s'il est explicitement défini, ou de la plus petite largeur ou hauteur de l'étendue des entités ponctuelles en entrée, dans la référence spatiale en entrée, divisée par 250.

Analysis Cell Size
search_radius
(Facultatif)

Indique les points en entrée utilisés pour interpoler la valeur de chaque cellule dans le raster en sortie.

Il existe deux moyens de spécifier la recherche de voisinage : Variable et Constant.

Variable utilise un rayon de recherche variable afin de rechercher un nombre spécifié de points d'échantillonnage en entrée pour l'interpolation. Constant utilise une distance constante spécifiée dans laquelle tous les points en entrée sont utilisés. Variable est l'option par défaut.

La syntaxe de ces paramètres est la suivante :

  • Variable, number_of_points, maximum_distance, où :
    • number_of_points - Valeur entière spécifiant le nombre de points d'échantillonnage en entrée les plus proches à utiliser pour effectuer l'interpolation. La valeur par défaut est 12.
    • maximum_distance - Spécifie la distance, en unités de carte, de limite de la recherche des points d'échantillonnage en entrée les plus proches. La valeur par défaut est la longueur de la diagonale de l'étendue.
  • Fixed, distance, minimum_number_of_points, où :
    • distance - Spécifie la distance en tant que rayon au sein duquel les points d'échantillonnage en entrée permettent d'effectuer l'interpolation. La valeur du rayon est exprimée en unités de carte. Le rayon par défaut est de 5 fois la taille de cellule du raster en sortie.
    • minimum_number_of_points – Entier définissant le nombre minimal de points à utiliser pour l'interpolation. La valeur par défaut est 0.

      Si le nombre de points requis n'est pas trouvé dans la distance spécifiée, la distance de recherche augmente jusqu'à ce que le nombre minimal de points spécifié soit trouvé.

      Le rayon de recherche est augmenté lorsque cela est nécessaire jusqu'à ce que la valeur de l'option minimum_number_of_points soit comprise dans ce rayon ou que l'étendue du rayon traverse l'étendue inférieure (sud) et/ou supérieure (nord) du raster en sortie. La valeur NoData est attribuée à tous les emplacements qui ne remplissent pas la condition susmentionnée.

Radius
out_variance_prediction_raster
(Facultatif)

Raster en sortie facultatif où chaque cellule contient les valeurs de semi-variance prévues pour cette localisation.

Raster Dataset

Exemple de code

1er exemple d'utilisation de l'outil Krigeage (fenêtre Python)

Cet exemple entre un fichier de formes ponctuelles et interpole la surface en sortie en tant que raster Grid.

import arcpy
from arcpy import env
env.workspace = "C:/data"
arcpy.Kriging_3d("ca_ozone_pts.shp", "OZONE", "c:/output/krigout",
                 "Spherical", 2000, "Variable 12")
2e exemple d'utilisation de l'outil Krigeage (script autonome)

Cet exemple entre un fichier de formes ponctuelles et interpole la surface en sortie en tant que raster Grid.

# Name: Kriging_3d_Ex_02.py
# Description: Interpolates a surface from points using kriging.
# Requirements: 3D Analyst Extension
# Import system modules

import arcpy
from arcpy import env

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

# Set local variables
inFeatures = "ca_ozone_pts.shp"
field = "OZONE"
outRaster = "C:/output/krigoutput02"
cellSize = 2000
outVarRaster = "C:/output/outvariance"
kModel = "CIRCULAR"
kRadius = 20000

# Check out the ArcGIS 3D Analyst extension license
arcpy.CheckOutExtension("3D")

# Execute Kriging
arcpy.Kriging_3d(inFeatures, field, outRaster, kModel, 
                 cellSize, kRadius, outVarRaster)

Environnements

Thèmes connexes

Informations de licence

ArcGIS for Desktop Basic: Requis 3D Analyst ou ArcGIS Spatial Analyst
ArcGIS for Desktop Standard: Requis 3D Analyst ou ArcGIS Spatial Analyst
ArcGIS for Desktop Advanced: Requis 3D Analyst ou ArcGIS Spatial Analyst
6/5/2014