『公告』 预祝您龙年大吉,万事如意, 过节期间, 大家如需数据服务,请拨打400 或直接添加客服微信,再祝大家龙年,心想事成。
关注我们 新浪 腾讯

ArcGIS教程:按圆形区域提取 (空间分析)

ArcGIS教程:按圆形区域提取 (空间分析)
基于圆提取栅格像元。

摘要

    基于圆提取栅格像元。

用法

·      可通过像元的中心来确定该像元是位于圆的内部还是位于圆的外部。如果中心位于圆弧的内部,则即使部分像元落在圆外,也会将此像元视为完全处于圆内。

·      未选择的像元位置被赋予 NoData 值。

·      当输入为多波段栅格时,将输出一个新的多波段栅格。对输入多波段栅格中的每一单个波段都会进行相应分析。

默认输出格式是 Esri Grid 堆栈。请注意,Esri Grid 堆栈的名称不能以数字开头、不能使用空格,并且长度不能超过 9 个字符。

·      如果输入内容为根据多波段栅格(超过三个波段)创建的图层,提取操作将只考虑图层加载(符号化)的波段。因此,输出的多波段栅格只会具有三个波段,分别对应于输入图层显示中使用的多波段栅格。

·      如果输入栅格为整型,那么输出栅格也为整型。如果输入栅格为浮点型,则输出栅格也为浮点型。

语法

    ExtractByCircle (in_raster, center_point, radius, {extraction_area})

参数

说明

数据类型

in_raster

提取像元的输入栅格。

Raster Layer

center_point

指示用于定义提取区域的圆的中心坐标 (x,y) 的点类。

该类的形式为:

·         Point (x, y)

将指定坐标使用与输入栅格相同的地图单位。

Point

radius

用于定义提取区域的圆半径。

将以地图单位指定半径,并且与输入栅格的单位相同。

Double

extraction_area

(可选)

标识要提取输入圆内部还是输入圆外部的像元。

·         INSIDE — 指定应选择输入圆内部的像元并将其写入输出栅格的关键字。圆形区域外部的所有像元都将在输出栅格中获得 NoData 值。

·         OUTSIDE — 指定应选择输入圆外部的像元并将其写入输出栅格的关键字。圆形区域内部的所有像元都将在输出栅格中获得 NoData 值。

String

返回值

名称

说明

数据类型

out_raster

包含提取像元值的输出栅格。

Raster

代码实例

    按圆形区域提取 (ExtractByCircle) 示例 1Python 窗口)

    本示例将在某一点位置方圆 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")

      京ICP备2025132830号-1 京公网安备 号