Disponible avec une licence Spatial Analyst.
Résumé
Extrait les cellules d'un raster selon un polygone.
Illustration
Utilisation
Pour effectuer une extraction selon un polygone dans une classe d'entités au lieu de fournir une série de paires x,y, vous pouvez utiliser l'outil Extraction par masque.
Le centre de la cellule permet de déterminer si une cellule est à l'intérieur ou à l'extérieur d'un polygone. Si le centre se trouve dans l'arc du polygone, on considère que la cellule entière se trouve à l'intérieur, même si ses proportions débordent de ce polygone.
Le polygone est limité à 1 000 sommets. Les sommets du polygone doivent être saisis dans le sens des aiguilles d'une montre. Le premier et le dernier sommet doivent être identiques. Cela permet de fermer chaque polygone lorsque vous en utilisez plusieurs. Si les derniers points ne sont pas identiques, le polygone sera fermé automatiquement. Les arcs du polygone peuvent s'entrecroiser, mais les polygones déformés ne sont pas recommandés.
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 un raster de géodatabase. Si le format en sortie est une pile Grid Esri, notez que le nom ne peut pas commencer par un nombre, utiliser des espaces ni comporter plus de neuf 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.
Pour plus d'informations sur les environnements de géotraitement qui s'appliquent à cet outil, reportez-vous à la rubrique Environnements d'analyse et Spatial Analyst.
Syntaxe
ExtractByPolygon (in_raster, polygon, {extraction_area})
Paramètre | Explication | Type de données |
in_raster | Raster en entrée à partir duquel les cellules seront extraites. | Raster Layer |
polygon [point,...] | Un polygone (ou plusieurs) qui définit la surface du raster en entrée à extraire. Chaque partie d'un polygone est une liste de sommets définie par les classes Point. Le cas échéant, une classe Polygon peut permettre de définir une liste de parties de polygone. Les points sont spécifiées en paires de coordonnées x,y dans les mêmes unités de carte que le raster en entrée. L'objet a la forme suivante :
Notez que la dernière coordonnée doit être la même que la première pour fermer un polygone. | Point |
extraction_area (Facultatif) | Détermine l'extraction des cellules à l'intérieur ou à l'extérieur du polygone en entrée.
| 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
Premier exemple d'utilisation de l'outil ExtractByPolygon (fenêtre Python)
Cet exemple décrit l'extraction des cellules d'un raster selon les coordonnées du polygone que vous avez indiquées.
import arcpy
from arcpy import env
from arcpy.sa import *
polyPoints = [arcpy.Point(743050, 4321275), arcpy.Point(743100, 4321200),
arcpy.Point(743500, 4322000),arcpy.Point(742900, 4321800)]
env.workspace = "C:/sapyexamples/data"
extPolygonOut = ExtractByPolygon("soil", polyPoints, "INSIDE")
extPolygonOut.save("c:/sapyexamples/output/extpoly")
Deuxième exemple d'utilisation de l'outil ExtractByPolygon (script autonome)
Cet exemple décrit l'extraction des cellules d'un raster selon les coordonnées du polygone que vous avez indiquées.
# Name: ExtractByPolgyon_Ex_02.py
# Description: Extracts the cells of a raster based on a polygon.
# 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"
polyPoints = [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 ExtractByPolygon
extPolygonOut = ExtractByPolygon(inRaster, polyPoints, "INSIDE")
# Save the output
extPolygonOut.save("c:/sapyexamples/output/extpoly02")
Environnements
Informations de licence
- ArcGIS Desktop Basic: Requiert Spatial Analyst
- ArcGIS Desktop Standard: Requiert Spatial Analyst
- ArcGIS Desktop Advanced: Requiert Spatial Analyst