Interpolation par noyaux avec interruptions (Geostatistical Analyst)

Disponible avec une licence Geostatistical Analyst.

Synthèse

Fenêtre mobile prédictive qui utilise la distance la plus courte entre les points de sorte que les points situés de part et d’autre des interruptions linéaires soient connectées.

Fonctionnement de l’outil Interpolation par noyau avec interruptions

Utilisation

  • L’entité d’interruption absolue emploie une approche basée sur une distance non euclidienne au lieu d’une ligne de visée. L’approche basée sur une ligne de visée nécessite qu’une ligne droite entre la localisation mesurée et la localisation où la prédiction est requise n’intersecte pas l’entité interruption. Si la distance autour de l’interruption se trouve dans les spécifications du voisinage de recherche, elle est considérée comme faisant partie de cette approche basée sur une distance non euclidienne.

  • Le temps de traitement dépend de la complexité de la géométrie des classes d’entités d’interruption. Les outils qui figurent dans le jeu d’outils Généralisation permettent de créer une nouvelle classe d’entités en lissant ou en supprimant certaines de ces entités.

  • Pour les fonctions de noyau Exponentiel, Gaussien et Constant, un facteur de lissage est appliqué de telle sorte que les noyaux possèdent un rayon fini égal à la bande passante spécifiée.

Paramètres

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

Entités ponctuelles en entrée contenant les valeurs z à interpoler.

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 en entrée contiennent des valeurs z ou des valeurs m.

Field
Couche géostatistique en sortie
(Facultatif)

Couche géostatistique produite. Cette couche est une sortie obligatoire uniquement si aucun raster en sortie n’est demandé.

Geostatistical Layer
Raster en sortie
(Facultatif)

Raster en sortie. Ce raster est une sortie obligatoire uniquement si aucune couche géostatistique en sortie n’est demandée.

Raster Dataset
Taille de cellule en sortie
(Facultatif)

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

Cette valeur peut être définie de manière explicite dans Environnements par le paramètre Taille de cellule.

Si elle n’est pas définie, il s’agit de la valeur la plus petite entre la largeur et la hauteur de l’étendue de l’entité ponctuelle en entrée, dans la référence spatiale en entrée, divisée par 250.

Analysis Cell Size
Entités d’interruption absolue en entrée
(Facultatif)

Entités d’interruption absolue utilisant des distances non euclidiennes plutôt que des distances de ligne de visée.

Feature Layer
Fonction de noyau
(Facultatif)

Fonction de noyau utilisée dans la simulation.

  • ExponentielLa fonction croît ou décroît proportionnellement.
  • GaussienFonction de courbe en cloche qui décroît rapidement vers moins l’infini/plus l’infini.
  • QuartiqueFonction polynomiale de quatrième ordre.
  • EpanechnikovFonction parabolique discontinue.
  • Polynôme de cinquième ordreFonction polynomiale de cinquième ordre.
  • ConstantFonction d’indicatrice.
String
Bande passante
(Facultatif)

Utilisée pour spécifier la distance maximale à laquelle les points de données sont utilisés pour la prédiction. Lorsque la bande passante augmente, le biais de prédiction augmente et la variance de prédiction diminue.

Double
Ordre du polynôme
(Facultatif)

Définit le degré du polynôme.

Long
Paramètre de régularisation
(Facultatif)

Utilisé pour la stabilisation numérique de la solution du système d’équations linéaires. Il n’a pas d’effet sur les prédictions en cas de données réparties de manière régulière sans interruptions. Les prédictions concernant des zones dans lesquelles les données se trouvent à proximité d’une entité interruption ou isolées par des interruptions peuvent être instables et ont tendance à nécessiter des valeurs de paramètre de régularisation relativement élevées.

Double
Type de surface en sortie
(Facultatif)

Type de surface pour stocker les résultats d’interpolation.

  • PrévisionLes surfaces de prédiction sont générées à partir des valeurs interpolées.
  • Erreur standard de prédictionLes surfaces d’erreur standard sont générées à partir des erreurs standard des valeurs interpolées.
String

arcpy.ga.KernelInterpolationWithBarriers(in_features, z_field, {out_ga_layer}, {out_raster}, {cell_size}, {in_barrier_features}, {kernel_function}, {bandwidth}, {power}, {ridge}, {output_type})
NomExplicationType de données
in_features

Entités ponctuelles en entrée contenant les valeurs z à interpoler.

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 en entrée contiennent des valeurs z ou des valeurs m.

Field
out_ga_layer
(Facultatif)

Couche géostatistique produite. Cette couche est une sortie obligatoire uniquement si aucun raster en sortie n’est demandé.

Geostatistical Layer
out_raster
(Facultatif)

Raster en sortie. Ce raster est une sortie obligatoire uniquement si aucune couche géostatistique en sortie n’est demandée.

Raster Dataset
cell_size
(Facultatif)

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

Cette valeur peut être définie de manière explicite dans Environnements par le paramètre Taille de cellule.

Si elle n’est pas définie, il s’agit de la valeur la plus petite entre la largeur et la hauteur de l’étendue de l’entité ponctuelle en entrée, dans la référence spatiale en entrée, divisée par 250.

Analysis Cell Size
in_barrier_features
(Facultatif)

Entités d’interruption absolue utilisant des distances non euclidiennes plutôt que des distances de ligne de visée.

Feature Layer
kernel_function
(Facultatif)

Fonction de noyau utilisée dans la simulation.

  • EXPONENTIALLa fonction croît ou décroît proportionnellement.
  • GAUSSIANFonction de courbe en cloche qui décroît rapidement vers moins l’infini/plus l’infini.
  • QUARTICFonction polynomiale de quatrième ordre.
  • EPANECHNIKOVFonction parabolique discontinue.
  • POLYNOMIAL5Fonction polynomiale de cinquième ordre.
  • CONSTANTFonction d’indicatrice.
String
bandwidth
(Facultatif)

Utilisée pour spécifier la distance maximale à laquelle les points de données sont utilisés pour la prédiction. Lorsque la bande passante augmente, le biais de prédiction augmente et la variance de prédiction diminue.

Double
power
(Facultatif)

Définit le degré du polynôme.

Long
ridge
(Facultatif)

Utilisé pour la stabilisation numérique de la solution du système d’équations linéaires. Il n’a pas d’effet sur les prédictions en cas de données réparties de manière régulière sans interruptions. Les prédictions concernant des zones dans lesquelles les données se trouvent à proximité d’une entité interruption ou isolées par des interruptions peuvent être instables et ont tendance à nécessiter des valeurs de paramètre de régularisation relativement élevées.

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 ?

  • PREDICTIONLes surfaces de prédiction sont générées à partir des valeurs interpolées.
  • PREDICTION_STANDARD_ERRORLes surfaces d’erreur standard sont générées à partir des erreurs standard des valeurs interpolées.
String

Exemple de code

Exemple 1 d’utilisation de l’outil KernelInterpolationWithBarriers (fenêtre Python)

Interpolez les entités ponctuelles sur un raster rectangulaire à l’aide d’une classe d’entités interruption.

import arcpy
arcpy.env.workspace = "C:/gapysamples/data"
arcpy.KernelInterpolationWithBarriers_ga("ca_ozone_pts", "OZONE", "outKIWB",
                                         "C:/gapyexamples/output/kiwbout", "2000",
                                         "ca_outline", "QUARTIC", "", "", "50", "PREDICTION")
Exemple 2 d’utilisation de l’outil KernelInterpolationWithBarriers (script Python autonome)

Interpolez les entités ponctuelles sur un raster rectangulaire à l’aide d’une classe d’entités interruption.

# Name: KernelInterpolationWithBarriers_Example_02.py
# Description: Kernel Interpolation with Barriers is a moving window predictor
#   that uses non-Euclidean distances.
# Requirements: Geostatistical Analyst Extension

# Import system modules
import arcpy

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

# Set local variables
inPointFeatures = "ca_ozone_pts.shp"
zField = "ozone"
outLayer = "outKIWB"
outRaster = "C:/gapyexamples/output/kiwbout"
cellSize = 2000.0
inBarrier = "ca_outline.shp"
kernelFunction = "QUARTIC"
bandwidth = ""
power = ""
ridgeParam = "50"
outputType = "PREDICTION"

# Execute KernelInterpolationWithBarriers
arcpy.KernelInterpolationWithBarriers_ga(inPointFeatures, zField, outLayer, outRaster,
                                         cellSize, inBarrier, kernelFunction, bandwidth,
                                         power, ridgeParam, outputType)

Informations de licence

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

Rubriques connexes