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é comme entrée, un nouveau raster multicanal est créé comme sortie. Chaque canal individuel dans le raster multicanal en entrée sera 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 des espaces ni comporter plus de 9 caractères.
-
Si l'entrée est une couche créée à partir d'un raster multicanal comportant plus de trois canaux, l'opération d'extraction considérera uniquement les canaux qui ont été chargés (symbolisés) par la couche. En conséquence, le raster multicanal en sortie peut avoir uniquement 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 à virgule flottante, la sortie sera également à 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")