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

You may also like...