ArcGIS教程:按圆形区域提取 (空间分析)
摘要
基于圆提取栅格像元。
用法
· 可通过像元的中心来确定该像元是位于圆的内部还是位于圆的外部。如果中心位于圆弧的内部,则即使部分像元落在圆外,也会将此像元视为完全处于圆内。
· 未选择的像元位置被赋予 NoData 值。
· 当输入为多波段栅格时,将输出一个新的多波段栅格。对输入多波段栅格中的每一单个波段都会进行相应分析。
默认输出格式是 Esri Grid 堆栈。请注意,Esri Grid 堆栈的名称不能以数字开头、不能使用空格,并且长度不能超过 9个字符。
· 如果输入内容为根据多波段栅格(超过三个波段)创建的图层,提取操作将只考虑图层加载(符号化)的波段。因此,输出的多波段栅格只会具有三个波段,分别对应于输入图层显示中使用的多波段栅格。
· 如果输入栅格为整型,那么输出栅格也为整型。如果输入栅格为浮点型,则输出栅格也为浮点型。
代码实例
按圆形区域提取 (ExtractByCircle) 示例 1(Python 窗口)
本示例将在某一点位置方圆 500 米的半径范围内提取像元。
import arcpy
from arcpy import env
from arcpy.sa import *
env.workspace = “C:/sapyexamples/data”
outExtCircle = ExtractByCircle(“elevation”, arcpy.Point(482838.823, 222128.982),
500, “INSIDE”)
outExtCircle.save(“c:/sapyexamples/output/extcircle”)
按圆形区域提取 (ExtractByCircle) 示例 2(独立脚本)
本示例将在某一点位置方圆 1,000 米的半径范围内提取像元。
# Name: ExtractByCircle_Ex_02.py
# Description: Extracts the cells of a raster based on a circle.
# 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 = (“elevation”)
centerPoint = arcpy.Point(482838.823, 222128.982)
circRadius = 1000
extractType = “INSIDE”
# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension(“Spatial”)
# Execute ExtractByCircle
outExtCircle = ExtractByCircle(inRaster, centerPoint, circRadius,
extractType)
# Save the output
outExtCircle.save(“c:/sapyexamples/output/extcircle02”)
转载自:https://blog.csdn.net/dsac1/article/details/49276329