Extraer por puntos (Spatial Analyst)
Resumen
Extrae las celdas de un ráster según un conjunto de puntos de coordenadas.
Uso
-
A las ubicaciones de las celdas que no están seleccionadas se les asigna el valor NoData.
-
Cuando un ráster multibanda se especifica como entrada, se creará un ráster multibanda nuevo como salida. Cada banda individual del ráster multibanda de entrada se analizará como corresponda.
El formato de salida predeterminado es una Pila de Esri Grid. Note que el nombre de una pila de Esri Grid no puede comenzar con un número, utilizar espacios o tener más de 9 caracteres de largo.
-
Si la entrada es una capa creada desde un ráster multibanda con más de tres bandas, la operación de extracción considerará únicamente las bandas que fueron cargadas (simbolizadas) por la capa. Como resultado, el ráster multibanda de salida sólo puede tener tres bandas, que corresponden a las utilizadas en la visualización de la capa de entrada.
-
Si el ráster de entrada es un entero, el ráster de salida será un entero. Si la entrada es de punto flotante, la salida será de punto flotante.
Sintaxis
| Parámetro | Explicación | Tipo de datos |
in_raster |
El ráster de entrada desde el cual se extraerán las celdas. | Raster Layer |
points [point,...] |
Una lista de Python de objetos de clase de Punto que denota las ubicaciones en donde se extraerán los valores del ráster. Los objetos de punto se especifican en una lista de pares de coordenadas x,y. El formato del objeto es:
Los puntos se encuentran en las mismas unidades de mapa que in_raster. | Point |
extraction_area (Opcional) |
Identifica si debe extraer celdas según los puntos de ubicación especificados (dentro) o fuera de los puntos de ubicación (fuera).
| String |
Valor de retorno
| Nombre | Explicación | Tipo de datos |
| out_raster |
Ráster de salida que contiene los valores de las celdas extraídos del ráster de entrada. | Raster |
Ejemplo de código
En este ejemplo se extraen celdas de un ráster en función de las coordenadas de puntos especificadas.
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")
En este ejemplo se extraen celdas de un ráster en función de las coordenadas de puntos especificadas.
# 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")