Voisin naturel (3D Analyst)
Récapitulatif
Interpole une surface raster à partir de points à l'aide de la méthode par voisins naturels.
Utilisation
-
Si le centre des cellules du périmètre du raster en sortie figure en dehors de l'enveloppe convexe (définie par les points en entrée), ces cellules ont la valeur NoData. Si un point en entrée est compris dans l'une de ces cellules du périmètre et que le centre de la cellule figure en dehors de l'enveloppe convexe, la cellule a la valeur NoData.
-
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.
-
Cet outil a une limite d'environ 15 millions de points en entrée. Si la classe d'entités en entrée contient un très important nombre de points (environ 15 millions ou plus), l'outil peut ne pas parvenir à générer un résultat.
Vous pouvez contourner cette limite en traitant votre zone d'étude en plusieurs sections et en réalisant un mosaïquage des résultats dans un seul grand jeu de données raster. Assurez-vous que les deux sections se superposent. Vous pouvez également utiliser un jeu de données de MNT pour stocker et visualiser des points et des surfaces contenus dans des milliards de points de mesure.
Il est préférable que les données en entrée se trouvent dans un système de coordonnées projetées plutôt que dans un système de coordonnées géographiques.
-
Vous pouvez également utiliser un jeu de données du TIN. Créez tout d'abord un TIN à partir de vos données source. Convertissez ensuite le TIN résultant en raster à l'aide de l'outil TIN vers raster, à l'aide de l'option Voisins naturels. Ceci s'avère particulièrement utile si vous avez des lignes de fracture ou une zone de données irrégulièrement façonnée.
Syntaxe
| 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_raster |
Raster de surface interpolé en sortie. | Raster Layer |
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 |
Exemple de code
Cet exemple entre un fichier de formes ponctuel et interpole la surface en sortie en tant que raster TIFF.
import arcpy
from arcpy import env
env.workspace = "C:/data"
arcpy.NaturalNeighbor_3d("ca_ozone_pts.shp", "ozone",
"C:/output/nnout.tif", 2000)
Cet exemple entre un fichier de formes ponctuelles et interpole la surface en sortie en tant que raster Grid.
# Name: NaturalNeighbor_3d_Ex_02.py
# Description: Interpolate a series of point features onto
# a rectangular raster using Natural Neighbor interpolation.
# Requirements: 3D Analyst Extension
# Import system modules
import arcpy
from arcpy import env
# Set environment settings
env.workspace = "C:/data"
# Set local variables
inPntFeat = "ca_ozone_pts.shp"
zField = "ozone"
outRaster = "C:/output/nnout"
cellSize = 40000
# Check out the ArcGIS 3D Analyst extension license
arcpy.CheckOutExtension("3D")
# Execute NaturalNeighbor
arcpy.NaturalNeighbor_3d(inPntFeat, zField, outRaster, cellSize)