用点提取 (空间分析)
摘要
基于一组坐标点提取栅格像元。
用法
-
未选择的像元位置被赋予 NoData 值。
-
当指定多波段栅格为输入时,创建一个新的多波段栅格作为输出。对输入多波段栅格中的每一单个波段都会进行相应的分析。
默认输出格式是 Esri Grid 堆栈。请注意,Esri Grid 堆栈的名称不能以数字开头、不能使用空格,并且长度不能超过 9 个字符。
-
如果输入内容为根据多波段栅格(超过三个波段)创建的图层,提取操作将只考虑图层加载(符号化)的波段。因此,输出的多波段栅格只能具有三个波段,分别对应于输入图层显示中使用的多波段栅格。
-
如果输入栅格为整型,那么输出栅格也为整型。如果输入栅格为浮点型,则输出栅格也为浮点型。
语法
ExtractByPoints (in_raster, points, {extraction_area})
参数 | 说明 | 数据类型 |
in_raster |
从其中提取像元的输入栅格。 | Raster Layer |
points [point,...] |
点类对象的 Python 列表用于指示要提取栅格值的位置。 点对象均在 x,y 坐标对列表中指定。对象形式为:
点所使用的地图单位与输入栅格相同。 | Point |
extraction_area (可选) |
标识是基于指定点位置(内部)提取像元还是基于点位置外部(外部)提取像元。
| String |
返回值
名称 | 说明 | 数据类型 |
out_raster |
包含从输入栅格中提取的像元值的输出栅格。 | Raster |
代码示例
用点提取 (ExtractByPoints) 示例 1(Python 窗口)
本示例将基于指定点坐标提取栅格的像元。
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")
用点提取 (ExtractByPoints) 示例 2(独立脚本)
本示例将基于指定点坐标提取栅格的像元。
# 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")
相关主题
许可信息
ArcView: 需要 Spatial Analyst
ArcEditor: 需要 Spatial Analyst
ArcInfo: 需要 Spatial Analyst
7/10/2012