Interpolation polynomiale locale (Outils Geostatistical Analyst)

Résumé

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 Local Polynomial Interpolation (Interpolation polynomiale locale)

Utilisation

Utilisez l’outil Local Polynomial Interpolation (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. Dans les sciences de la Terre toutefois, la variable d’intérêt présente généralement une variation à court terme en plus d’une tendance à long terme. Lorsque l’ensemble de données présente une variation à court terme, les cartes Local Polynomial Interpolation (Interpolation polynomiale locale) peuvent capturer la variation à court terme.

Etiqueter	Explication	Type 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 Environments (Environnements) par le paramètre Cell Size (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
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. Standard Major semiaxis (Demi-grand axe) : valeur du demi-grand axe du voisinage de recherche. Minor semiaxis (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. Max neighbors (Nbre max. de voisins) : nombre maximal de voisins qui sont utilisés pour estimer la valeur à la localisation inconnue. Min neighbors (Nbre min. de voisins) : nombre minimal de voisins qui sont utilisés pour estimer la valeur à la localisation inconnue. Sector Type (Type de secteur) : géométrie du voisinage. One sector (Un secteur) : ellipse unique. Four sectors (Quatre secteurs) : ellipse divisée en quatre secteurs. Four sectors shifted (Quatre secteurs orientés) : ellipse divisée en quatre secteurs et orientée à 45 degrés Eight sectors (Huit secteurs) : ellipse divisée en huit secteurs. Lisser Major semiaxis (Demi-grand axe) : valeur du demi-grand axe du voisinage de recherche. Minor semiaxis (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. Smoothing factor (Facteur de lissage) : l’option Smooth Interpolation (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 (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. Max neighbors (Nbre max. de voisins) : nombre maximal de voisins qui sont utilisés pour estimer la valeur à la localisation inconnue. Min neighbors (Nbre min. de voisins) : nombre minimal de voisins qui sont utilisés pour estimer la valeur à la localisation inconnue. Sector Type (Type de secteur) : géométrie du voisinage. One sector (Un secteur) : ellipse unique. Four sectors (Quatre secteurs) : ellipse divisée en quatre secteurs. Four sectors shifted (Quatre secteurs orientés) : ellipse divisée en quatre secteurs et orientée à 45 degrés Eight sectors (Huit secteurs) : ellipse divisée en huit secteurs. Circulaire lissé Radius (Rayon) : longueur du rayon du cercle de recherche. Smoothing factor (Facteur de lissage) : l’option Smooth Interpolation (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. Exponentiel—La fonction croît ou décroît proportionnellement. Gaussien—Fonction de courbe en cloche qui décroît rapidement vers moins l’infini ou plus l’infini. Quartique—Fonction polynomiale de quatrième ordre. Epanechnikov—Fonction parabolique discontinue. Polynôme de cinquième ordre—Fonction polynomiale de cinquième ordre. Constante—Fonction d’indicatrice.	String
Bandwidth (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. Activé—Les prédictions et les erreurs standards de prédiction ne sont pas créées là où les prédictions sont instables. Décoché—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.	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évision—Les surfaces de prédiction sont générées à partir des valeurs interpolées. Erreur standard de prévision—Les surfaces d’erreur standard sont générées à partir des erreurs standard des valeurs interpolées. Nombre de conditions—La 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})

Nom	Explication	Type 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 Environments (Environnements) par le paramètre Cell Size (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 Search Neighborhood (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 sont utilisés pour estimer la valeur à la localisation inconnue. `nbrMin`—Nombre minimal de voisins qui sont 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 Smooth Interpolation (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 sont utilisés pour estimer la valeur à la localisation inconnue. `nbrMin`—Nombre minimal de voisins qui sont 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 Smooth Interpolation (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. `EXPONENTIAL`—La fonction croît ou décroît proportionnellement. `GAUSSIAN`—Fonction de courbe en cloche qui décroît rapidement vers moins l’infini ou plus l’infini. `QUARTIC`—Fonction polynomiale de quatrième ordre. `EPANECHNIKOV`—Fonction parabolique discontinue. `POLYNOMIAL5`—Fonction polynomiale de cinquième ordre. `CONSTANT`—Fonction 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. `USE_CONDITION_NUMBER`—Les prédictions et les erreurs standards de prédiction ne sont pas créées là où les prédictions sont instables. `NO_USE_CONDITION_NUMBER`—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.	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 ? `PREDICTION`—Les surfaces de prédiction sont générées à partir des valeurs interpolées. `PREDICTION_STANDARD_ERROR`—Les surfaces d’erreur standard sont générées à partir des erreurs standard des valeurs interpolées. `CONDITION_NUMBER`—La 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.ga.LocalPolynomialInterpolation(
    "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)

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

Environnements

Taille de cellule, Points coïncidents, Espace de travail courant, Etendue, Transformations géographiques, Masque, Système de coordonnées en sortie, Facteur de traitement parallèle, Raster de capture

Informations de licence