Krigeage (3D Analyst)

Disponible avec une licence Spatial Analyst.

Disponible avec une licence 3D Analyst.

Synthèse

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

L’outil Krigeage bayésien empirique fournit des fonctionnalités ou des performances améliorées.

Pour en savoir plus sur la fonction Krigeage

Utilisation

  • Le processus de krigeage nécessite d'importantes ressources du processeur. La vitesse d'exécution dépend du nombre de points dans le jeu de données en entrée et la taille de la fenêtre de recherche.

  • Les valeurs faibles comprises dans la variance en sortie facultative du raster de prévision indiquent un haut degré de confiance dans la valeur de prévision. Des valeurs élevées peuvent signifier que davantage de points de données sont nécessaires.

  • Les types de krigeage universel supposent qu'il existe un composant structurel et que la tendance locale varie d'un emplacement à un autre.

  • Les Propriétés du semi-variogramme permettent de contrôler le semi-variogramme utilisé pour le krigeage. A l'origine, une Valeur de décalage par défaut est définie sur la taille de cellule en sortie par défaut. Les valeurs d'Etendue principale, de Seuil partiel et de Pépite sont calculées par défaut en interne si rien n'est spécifié.

  • La variance en sortie facultative du raster de prévision contient la variance de krigeage au niveau de chaque cellule de raster en sortie. Si l'on part du principe que les erreurs de krigeage sont normalement distribuées, il existe 95,5 % de probabilité que la valeur z réelle au niveau de la cellule corresponde à la valeur du raster prévue, plus ou moins deux fois la racine carrée de la valeur dans le raster de variance.

  • Le paramètre Output cell size (Taille de cellule en sortie) peut être défini par une valeur numérique ou obtenu à partir d’un jeu de données raster existant. Si la taille de cellule n’est pas explicitement spécifiée en tant que valeur de paramètre, elle est dérivée de l’environnement Cell Size (Taille de cellule), s’il est spécifié. Si la taille de cellule du paramètre ou la taille de cellule de l’environnement n’ont pas été spécifiées, mais que l’environnement Raster de capture a été défini, la taille de cellule du raster de capture est utilisée. Si aucun élément n’est spécifié, la taille de cellule résulte de la division par 250 de la largeur ou de la hauteur la plus petite (parmi les deux) de l’étendue spécifiée dans le système de coordonnées en sortie de l’environnement.

  • Si la taille de cellule est spécifiée à l’aide d’une valeur numérique, l’outil l’utilise directement pour le raster en sortie.

    Si la taille de cellule est spécifiée à l’aide d’un jeu de données raster, le paramètre affiche le chemin du jeu de données raster au lieu de la valeur de la taille de cellule. La taille de cellule de ce jeu de données raster sera utilisée directement dans l’analyse, à condition que la référence spatiale du jeu de données soit la même que la référence spatiale en sortie. Si la référence spatiale du jeu de données est différente de la référence spatiale en sortie, elle sera projetée en fonction de la valeur Cell Size Projection Method (Méthode de projection de la taille de cellule) spécifiée.

  • Certains jeux de données en entrée peuvent avoir plusieurs points avec les mêmes coordonnées x,y. Si les valeurs des points à l'emplacement commun sont les mêmes, elles sont considérées comme étant en double et n'ont aucune incidence sur la sortie. Si les valeurs sont différentes, elles sont considérées comme étant des points coïncidents.

    Les différents outils d'interpolation peuvent gérer cette condition de données différemment. Par exemple dans certains cas, le premier point coïncident détecté est utilisé pour le calcul, tandis que dans d’autres cas, c’est le dernier point détecté qui est utilisé. Cela peut entraîner des valeurs inattendues pour les emplacements du raster en sortie. La solution consiste à préparer vos données en supprimant ces points coïncidents. L'outil Collect Events de la boîte à outils Outils de statistiques spatiales sert à identifier tous points coïncidents de vos données.

  • Pour les formats de données prenant en charge les valeurs nulles, par exemple les classes d’entités d’une géodatabase fichier, une valeur nulle est ignorée si elle est utilisée en entrée.

Paramètres

ÉtiquetteExplicationType de données
Entités ponctuelles en entrée

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

Feature Layer
Champ de la valeur Z

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
Raster de surface en sortie

Raster de surface interpolé en sortie.

Il s'agit toujours d'un raster à virgule flottante.

Raster Dataset
Propriétés du semi-variogramme

Modèle de semi-variogramme à utiliser. Il existe deux méthodes de krigeage : ordinaire et universelle.

Le krigeage ordinaire peut utiliser les modèles de semi-variogramme suivants :

  • SphericalModèle de semi-variogramme sphérique. Il s’agit de l’option par défaut.
  • CircularModèle de semi-variogramme circulaire.
  • ExponentialModèle de semi-variogramme exponentiel.
  • GaussianModèle de semi-variogramme gaussien ou de distribution normale.
  • LinearModèle de semi-variogramme linéaire avec un seuil.

Le krigeage universel peut utiliser les modèles de semi-variogramme suivants :

  • Linear with Linear driftKrigeage universel avec dérive linéaire.
  • Linear with Quadratic driftKrigeage universel avec dérive quadratique.

Des options sont disponibles dans la boîte de dialogue Advanced Parameters (Paramètres avancés). Ces paramètres sont les suivants :

  • Lag sizeLa valeur par défaut est la taille de cellule du raster en sortie.
  • Major rangeReprésente une distance au-delà de laquelle il n'y a que peu ou pas de corrélation.
  • Partial sillDifférence entre la pépite et le seuil.
  • NuggetRepré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.
KrigingModel
Taille de cellule en sortie
(Facultatif)

Taille de cellule du raster en sortie à créer.

Ce paramètre peut être défini par une valeur numérique ou obtenu à partir d’un jeu de données raster existant. Si la taille de cellule n’est pas explicitement spécifiée en tant que valeur de paramètre, la valeur de taille de cellule de l’environnement est utilisée si elle est spécifiée. Sinon, des règles supplémentaires sont utilisées pour la calculer à partir des autres entrées. Consultez la section Utilisation pour en savoir plus.

Analysis Cell Size
Rayon de recherche
(Facultatif)

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

Deux options sont disponibles : Variable et Constant. Variable est l'option par défaut.

  • Variable

    Utilise un rayon de recherche variable afin de rechercher un nombre spécifié de points d'échantillonnage en entrée pour l'interpolation.

    • Nombre de 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.
    • Distance maximale - 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.
  • Constant

    Utilise une distance constante spécifiée dans laquelle tous les points en entrée sont utilisés pour l'interpolation.

    • 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.

    • Nombre minimal de 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 Nombre minimal de 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
Raster de variance de prédiction en sortie
(Facultatif)

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

Raster Dataset

arcpy.ddd.Kriging(in_point_features, z_field, out_surface_raster, semiVariogram_props, {cell_size}, {search_radius}, {out_variance_prediction_raster})
NomExplicationType 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.

Il s'agit toujours d'un raster à virgule flottante.

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. On compte les suivantes :
    • 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 de cellule du raster en sortie à créer.

Ce paramètre peut être défini par une valeur numérique ou obtenu à partir d’un jeu de données raster existant. Si la taille de cellule n’est pas explicitement spécifiée en tant que valeur de paramètre, la valeur de taille de cellule de l’environnement est utilisée si elle est spécifiée. Sinon, des règles supplémentaires sont utilisées pour la calculer à partir des autres entrées. Consultez la section Utilisation pour en savoir plus.

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 Fixed.

Variable utilise un rayon de recherche variable afin de rechercher le nombre spécifié de points d’échantillonnage en entrée pour l’interpolation. Fixed utilise une distance constante spécifiée dans laquelle tous les points en entrée sont utilisés. Variable est la valeur 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 variance prévues pour cette localisation.

Raster Dataset

Exemple de code

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

Cet exemple prend en entrée un shapefile ponctuel et interpole la surface en sortie sous forme de 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 prend en entrée un shapefile ponctuel et interpole la surface en sortie sous forme de 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

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

Rubriques connexes