Interpolation polynomiale locale (Geostatistical Analyst)

Disponible avec une licence Geostatistical Analyst.

Synthèse

Ajuste le polynôme d’ordre spécifié (zéro, premier, second, troisième et ainsi de suite), chacun figurant dans les voisinages superposés spécifiés, pour produire une surface en sortie.

Fonctionnement de l’outil Interpolation polynomiale locale

Utilisation

  • Utilisez l’outil Interpolation polynomiale locale lorsque votre jeu de données montre une variation à court terme.

  • L’interpolation polynomiale globale est utile pour créer des surfaces lisses et identifier des tendances à long terme dans le jeu de données. Cependant, dans le domaine des sciences de la Terre, la variable d’intérêt a généralement une variation à court terme en plus d’une tendance à long terme. Lorsque le jeu de données présente une variation à court terme, les cartes d’interpolation polynomiale locale peuvent capturer la variation à court terme.

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
Ordre du polynôme
(Facultatif)

Ordre de la tendance polynomiale.

Long
Voisinage de recherche
(Facultatif)

Définit les points de voisinage qui seront utilisés pour contrôler la sortie. Standard est la valeur par défaut.

Standard

  • Demi-grand axe - Valeur du demi-grand axe du voisinage de recherche.
  • Demi-petit axe - Valeur de la taille du demi-petit axe du voisinage de recherche.
  • Angle — Angle de rotation de l’axe (cercle) ou du demi-grand axe (ellipse) de la fenêtre mobile.
  • Nbre max. de voisins - Nombre maximal de voisins qui seront utilisés pour estimer la valeur à la localisation inconnue.
  • Nbre min. de voisins - Nombre minimal de voisins qui seront utilisés pour estimer la valeur à la localisation inconnue.
  • Type de secteur — Géométrie du voisinage.
    • Un secteur — Ellipse unique.
    • Quatre secteurs — Ellipse divisée en quatre secteurs.
    • Quatre secteurs orientés - Ellipse divisée en quatre secteurs et orientée à 45 degrés
    • Huit secteurs — Ellipse divisée en huit secteurs.

Lisser

  • Demi-grand axe - Valeur du demi-grand axe du voisinage de recherche.
  • Demi-petit axe - Valeur de la taille du demi-petit axe du voisinage de recherche.
  • Angle — Angle de rotation de l’axe (cercle) ou du demi-grand axe (ellipse) de la fenêtre mobile.
  • Facteur de lissage — L’option Interpolation lisse crée une ellipse extérieure et une ellipse intérieure à une distance égale au semi-grand axe multiplié par le facteur de lissage. Les points se trouvant à l’extérieur de la plus petite ellipse, mais à l’intérieur de la plus grande ellipse sont pondérés à l’aide d’une fonction sigmoïde avec une valeur comprise entre zéro et un.

Circulaire standard

  • Rayon — Longueur du rayon du cercle de recherche.
  • Angle — Angle de rotation de l’axe (cercle) ou du demi-grand axe (ellipse) de la fenêtre mobile.
  • Nbre max. de voisins - Nombre maximal de voisins qui seront utilisés pour estimer la valeur à la localisation inconnue.
  • Nbre min. de voisins - Nombre minimal de voisins qui seront utilisés pour estimer la valeur à la localisation inconnue.
  • Type de secteur — Géométrie du voisinage.
    • Un secteur — Ellipse unique.
    • Quatre secteurs — Ellipse divisée en quatre secteurs.
    • Quatre secteurs orientés - Ellipse divisée en quatre secteurs et orientée à 45 degrés
    • Huit secteurs — Ellipse divisée en huit secteurs.

Circulaire lissé

  • Rayon — Longueur du rayon du cercle de recherche.
  • Facteur de lissage — L’option Interpolation lisse crée une ellipse extérieure et une ellipse intérieure à une distance égale au semi-grand axe multiplié par le facteur de lissage. Les points se trouvant à l’extérieur de la plus petite ellipse, mais à l’intérieur de la plus grande ellipse sont pondérés à l’aide d’une fonction sigmoïde avec une valeur comprise entre zéro et un.
Geostatistical Search Neighborhood
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 ou 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
Utiliser le seuil de l’indice de conditionnement spatial
(Facultatif)

Option destinée à contrôler la création de la prédiction et des erreurs standards de prédiction là où les prédictions sont instables. Cette option n’est disponible que pour les polynômes de 1er, 2e et 3e ordre.

  • Désactivé : Les prédictions sont créées partout, y compris dans les zones où elles sont instables. Il s’agit de l’option par défaut.
  • Activé : Les prédictions et les erreurs standards de prédiction ne sont pas créées là où les prédictions sont instables.
Boolean
Seuil de l’indice de conditionnement spatial
(Facultatif)

Chaque matrice carrée inversible présente un indice de conditionnement indiquant le degré d’inexactitude de la solution des équations linéaires avec une petite variation des coefficients de la matrice (cela peut être lié à des données plus ou moins précises). Si l’indice de conditionnement est élevé, une faible variation des coefficients de la matrice entraîne une modification importante dans le vecteur solution.

Double
Pondération de champ
(Facultatif)

Utilisée pour mettre une observation en évidence. Plus la pondération est importante, plus l’impact est important sur la prédiction. Pour les observations coïncidentes, attribuer la pondération la plus importante à la mesure la plus fiable.

Field
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.
  • Nombre de conditionsLa surface de l’indice de conditionnement spatial indique la stabilité des calculs en une localisation particulière. Plus l’indice de conditionnement est élevé, plus la prédiction est instable ; les localisations affichant des indices de conditionnement élevés sont susceptibles de contenir des artefacts et des valeurs prédites irrégulières.
String

arcpy.ga.LocalPolynomialInterpolation(in_features, z_field, {out_ga_layer}, {out_raster}, {cell_size}, {power}, {search_neighborhood}, {kernel_function}, {bandwidth}, {use_condition_number}, {condition_number}, {weight_field}, {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
power
(Facultatif)

Ordre de la tendance polynomiale.

Long
search_neighborhood
(Facultatif)

Définit les points de voisinage qui seront utilisés pour contrôler la sortie. Standard est la valeur par défaut.

Les éléments suivants sont des classes Voisinage de recherche : SearchNeighborhoodStandard, SearchNeighborhoodSmooth, SearchNeighborhoodStandardCircular et SearchNeighborhoodSmoothCircular.

Standard

  • majorSemiaxis - Valeur de la taille du demi-grand axe du voisinage de recherche.
  • minorSemiaxis - Valeur du demi-petit axe du voisinage de recherche.
  • angle — Angle de rotation de l’axe (cercle) ou du demi-grand axe (ellipse) de la fenêtre mobile.
  • nbrMax - Nombre maximal de voisins qui seront utilisés pour estimer la valeur à la localisation inconnue.
  • nbrMin - Nombre minimal de voisins qui seront utilisés pour estimer la valeur à la localisation inconnue.
  • sectorType - Géométrie du voisinage.
    • ONE_SECTOR - Ellipse unique.
    • FOUR_SECTORS - Ellipse divisée en quatre secteurs.
    • FOUR_SECTORS_SHIFTED - Ellipse divisée en quatre secteurs et orientée à 45 degrés
    • EIGHT_SECTORS - Ellipse divisée en huit secteurs.

Lisser

  • majorSemiaxis - Valeur de la taille du demi-grand axe du voisinage de recherche.
  • minorSemiaxis - Valeur du demi-petit axe du voisinage de recherche.
  • angle — Angle de rotation de l’axe (cercle) ou du demi-grand axe (ellipse) de la fenêtre mobile.
  • smoothFactor — L’option Interpolation lisse crée une ellipse extérieure et une ellipse intérieure à une distance égale au semi-grand axe multiplié par le facteur de lissage. Les points se trouvant à l’extérieur de la plus petite ellipse, mais à l’intérieur de la plus grande ellipse sont pondérés à l’aide d’une fonction sigmoïde avec une valeur comprise entre zéro et un.

Circulaire standard

  • radius — Longueur du rayon du cercle de recherche.
  • angle — Angle de rotation de l’axe (cercle) ou du demi-grand axe (ellipse) de la fenêtre mobile.
  • nbrMax - Nombre maximal de voisins qui seront utilisés pour estimer la valeur à la localisation inconnue.
  • nbrMin - Nombre minimal de voisins qui seront utilisés pour estimer la valeur à la localisation inconnue.
  • sectorType - Géométrie du voisinage.
    • ONE_SECTOR - Ellipse unique.
    • FOUR_SECTORS - Ellipse divisée en quatre secteurs.
    • FOUR_SECTORS_SHIFTED - Ellipse divisée en quatre secteurs et orientée à 45 degrés
    • EIGHT_SECTORS - Ellipse divisée en huit secteurs.

Circulaire lissé

  • radius — Longueur du rayon du cercle de recherche.
  • smoothFactor — L’option Interpolation lisse crée une ellipse extérieure et une ellipse intérieure à une distance égale au semi-grand axe multiplié par le facteur de lissage. Les points se trouvant à l’extérieur de la plus petite ellipse, mais à l’intérieur de la plus grande ellipse sont pondérés à l’aide d’une fonction sigmoïde avec une valeur comprise entre zéro et un.
Geostatistical Search Neighborhood
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 ou 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
use_condition_number
(Facultatif)

Option destinée à contrôler la création de la prédiction et des erreurs standards de prédiction là où les prédictions sont instables. Cette option n’est disponible que pour les polynômes de 1er, 2e et 3e ordre.

  • NO_USE_CONDITION_NUMBERLes prédictions sont créées partout, y compris dans les zones où elles sont instables. Il s’agit de l’option par défaut.
  • USE_CONDITION_NUMBERLes prédictions et les erreurs standards de prédiction ne sont pas créées là où les prédictions sont instables.
Boolean
condition_number
(Facultatif)

Chaque matrice carrée inversible présente un indice de conditionnement indiquant le degré d’inexactitude de la solution des équations linéaires avec une petite variation des coefficients de la matrice (cela peut être lié à des données plus ou moins précises). Si l’indice de conditionnement est élevé, une faible variation des coefficients de la matrice entraîne une modification importante dans le vecteur solution.

Double
weight_field
(Facultatif)

Utilisée pour mettre une observation en évidence. Plus la pondération est importante, plus l’impact est important sur la prédiction. Pour les observations coïncidentes, attribuer la pondération la plus importante à la mesure la plus fiable.

Field
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.
  • CONDITION_NUMBERLa surface de l’indice de conditionnement spatial indique la stabilité des calculs en une localisation particulière. Plus l’indice de conditionnement est élevé, plus la prédiction est instable ; les localisations affichant des indices de conditionnement élevés sont susceptibles de contenir des artefacts et des valeurs prédites irrégulières.
String

Exemple de code

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

Interpolez les entités ponctuelles sur un raster rectangulaire.

import arcpy
arcpy.env.workspace = "C:/gapyexamples/data"
arcpy.LocalPolynomialInterpolation_ga(
    "ca_ozone_pts", "OZONE", "outLPI", "C:/gapyexamples/output/lpiout", "2000",
    "2", arcpy.SearchNeighborhoodSmooth(300000, 300000, 0, 0.5), "QUARTIC", 
    "", "", "", "", "PREDICTION")
Exemple 2 d’utilisation de la fonction LocalPolynomialInterpolation (script autonome)

Interpolez les entités ponctuelles sur un raster rectangulaire.

# Name: LocalPolynomialInterpolation_Example_02.py
# Description: Local Polynomial interpolation fits many polynomials, each 
#              within specified overlapping neighborhoods. 
# 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 = "outLPI"
outRaster = "C:/gapyexamples/output/lpiout"
cellSize = 2000.0
power = 2
kernelFunction = "QUARTIC"
bandwidth = ""
useConNumber = ""
conNumber = ""
weightField = ""
outSurface = "PREDICTION"

# Set variables for search neighborhood
majSemiaxis = 300000
minSemiaxis = 300000
angle = 0
smoothFactor = 0.5
searchNeighbourhood = arcpy.SearchNeighborhoodSmooth(majSemiaxis, minSemiaxis,
                                                     angle, smoothFactor)


# Execute LocalPolynomialInterpolation
arcpy.LocalPolynomialInterpolation_ga(inPointFeatures, zField, outLayer, outRaster,
                                      cellSize, power, searchNeighbourhood,
                                      kernelFunction, bandwidth, useConNumber,
                                      conNumber, weightField, outSurface)

Informations de licence

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

Rubriques connexes