Extraction par points (Spatial Analyst)
Récapitulatif
Extrait les cellules d'un raster selon un ensemble de points de coordonnées.
Utilisation
-
La valeur NoData est attribuée aux emplacements de cellule qui ne sont pas sélectionnés.
-
Lorsqu'un raster multicanal est spécifié en entrée, un nouveau raster multicanal est créé en sortie. Chaque canal individuel au sein du raster multicanal en entrée est analysé en conséquence.
Le format en sortie par défaut est une pile de grilles Esri. Notez que le nom d'une pile de grilles Esri ne peut pas commencer par un nombre, utiliser d'espaces ou comporter plus de 9 caractères.
-
Si l'entrée est une couche créée à partir d'un raster multicanal de plus de trois canaux, l'opération d'extraction ne prend en compte que les canaux chargés (symbolisés) par la couche. En conséquence, le raster multicanal en sortie ne peut comporter que trois canaux, qui correspondent à ceux utilisés dans l'affichage de la couche en entrée.
-
Si le raster en entrée est un nombre entier, le raster en sortie sera un nombre entier. Si l'entrée est en virgule flottante, la sortie est également en virgule flottante.
Syntaxe
| Paramètre | Explication | Type de données |
in_raster |
Raster en entrée à partir duquel les cellules seront extraites. | Raster Layer |
points [point,...] |
Une liste Python de classes d'objets Point représente les emplacements où les valeurs sont extraites du raster. Les points sont indiqués dans une liste de paires de coordonnées x,y. L'objet a la forme suivante :
Les points sont spécifiés dans les mêmes unités de carte que le paramètre in_raster. | Point |
extraction_area (Facultatif) |
Indique s'il faut extraire des cellules selon les emplacements des points spécifiés (à l'intérieur) ou à l'extérieur des emplacements des points (à l'extérieur).
| String |
Valeur renvoyée
| Nom | Explication | Type de données |
| out_raster |
Raster en sortie contenant les valeurs de cellule extraites du raster en entrée. | Raster |
Exemple de code
Cet exemple décrit l'extraction des cellules d'un raster selon les coordonnées de points spécifiées.
import arcpy
from arcpy import env
from arcpy.sa import *
env.workspace = "C:/sapyexamples/data"
pointList = [arcpy.Point(743050, 4321275), arcpy.Point(743100, 4321200),
arcpy.Point(734500,4322000)]
outPointExtract = ExtractByPoints("soil", pointList,"INSIDE")
outPointExtract.save("c:/sapyexamples/output/pntextract")
Cet exemple décrit l'extraction des cellules d'un raster selon les coordonnées de points spécifiées.
# Name: ExtractByPoints_Ex_02.py
# Description: Extracts the cells of a raster based on a set of points.
# Requirements: Spatial Analyst Extension
# Import system modules
import arcpy
from arcpy import env
from arcpy.sa import *
# Set environment settings
env.workspace = "C:/sapyexamples/data"
# Set local variables
inRaster = "soil"
pointList = [arcpy.Point(743050, 4321275),
arcpy.Point(743100, 4321200),
arcpy.Point(743500, 4322000),
arcpy.Point(742900, 4321800)]
# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")
# Execute ExtractByPoints
outPointExtract = ExtractByPoints("soil", pointList,"INSIDE")
# Save the output
outPointExtract.save("c:/sapyexamples/output/pntext")