Disponible avec une licence Spatial Analyst.
Disponible avec une licence 3D Analyst.
Résumé
Interpole une surface raster à partir de points à l'aide du 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 paramètres avancés 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.
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.
Syntaxe
Kriging_3d (in_point_features, z_field, out_surface_raster, semiVariogram_props, {cell_size}, {search_radius}, {out_variance_prediction_raster})| Paramètre | Explication | Type 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.
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 Fixed. Variable utilise un rayon de recherche variable afin de rechercher un 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 :
| 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 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
Informations de licence
- ArcGIS Desktop Basic: Requiert 3D Analyst ou Spatial Analyst
- ArcGIS Desktop Standard: Requiert 3D Analyst ou Spatial Analyst
- ArcGIS Desktop Advanced: Requiert 3D Analyst ou Spatial Analyst