IDW (Geostatistical Analyst)

Disponible avec une licence Geostatistical Analyst.

Synthèse

Utilise l’interpolation pondérée par l’inverse de la distance pour prévoir la valeur d’une localisation non échantillonnée à l’aide des valeurs mesurées autour de la localisation de prévision. L’outil part de l’hypothèse que les objets qui sont proches les uns des autres se ressemblent davantage que ceux qui sont plus éloignés.

En savoir plus sur la fonction IDW

Utilisation

  • La valeur prévue se limite à la plage des valeurs utilisées pour l’interpolation. L’interpolation pondérée par l’inverse de la distance étant une distance moyenne pondérée, celle-ci ne peut pas être supérieure à la valeur de données en entrée la plus élevée, ni inférieure à la valeur de données en entrée la plus faible. Par exemple, elle ne peut pas prévoir les crêtes ou les vallées si le haut de la crête ou le bas de la vallée n’ont pas encore été échantillonnés.

  • L’interpolation pondérée de l’inverse de la distance peut générer des contours à de faibles distances autour de la localisation des données (effet connu sous le nom d’effet œil-de-bœuf).

  • Contrairement aux autres méthodes d’interpolation telles que le krigeage, l’interpolation pondérée par l’inverse de la distance n’émet aucune hypothèse explicite concernant les propriétés statistiques des données en entrée. L’interpolation pondérée par l’inverse de la distance est couramment utilisée lorsque les données en entrée ne répondent pas aux hypothèses statistiques des méthodes d’interpolation plus avancées.

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
Puissance
(Facultatif)

Exposant de la distance qui définit la rapidité selon laquelle les pondérations diminuent en fonction de la distance. Une puissance plus élevée engendre une moindre influence des points voisins distants. La valeur doit être supérieure ou égale à 1. Si aucune valeur n’est fournie, une valeur est estimée pendant l’exécution de l’outil et affichée dans un message. La valeur estimée est la valeur qui génère la plus petite erreur quadratique moyenne de validation croisée.

Double
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
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

arcpy.ga.IDW(in_features, z_field, {out_ga_layer}, {out_raster}, {cell_size}, {power}, {search_neighborhood}, {weight_field})
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)

Exposant de la distance qui définit la rapidité selon laquelle les pondérations diminuent en fonction de la distance. Une puissance plus élevée engendre une moindre influence des points voisins distants. La valeur doit être supérieure ou égale à 1. Si aucune valeur n’est fournie, une valeur est estimée pendant l’exécution de l’outil et affichée dans un message. La valeur estimée est la valeur qui génère la plus petite erreur quadratique moyenne de validation croisée.

Double
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
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

Exemple de code

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

Interpoler une série d’entités ponctuelles sur un raster.

import arcpy
arcpy.env.workspace = "C:/gapyexamples/data"
arcpy.ga.IDW("ca_ozone_pts", "OZONE", "outIDW", "C:/gapyexamples/output/idwout", "2000", "2",
             arcpy.SearchNeighborhoodStandard(300000, 300000, 0, 15, 10, "ONE_SECTOR"), "")
2e exemple d'utilisation de l'outil IDW (script autonome)

Interpoler une série d’entités ponctuelles sur un raster.

# Name: InverseDistanceWeighting_Example_02.py
# Description: Interpolate a series of point features onto a rectangular raster
#              using Inverse Distance Weighting (IDW).
# 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 = "outIDW"
outRaster = "C:/gapyexamples/output/idwout"
cellSize = 2000.0
power = 2

# Set variables for search neighborhood
majSemiaxis = 300000
minSemiaxis = 300000
angle = 0
maxNeighbors = 15
minNeighbors = 10
sectorType = "ONE_SECTOR"
searchNeighbourhood = arcpy.SearchNeighborhoodStandard(majSemiaxis, minSemiaxis,
                                                       angle, maxNeighbors,
                                                       minNeighbors, sectorType)

# Execute IDW
arcpy.ga.IDW(inPointFeatures, zField, outLayer, outRaster, cellSize, 
             power, searchNeighbourhood)

Informations de licence

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

Rubriques connexes