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

ArcGIS教程:用矩形提取 (空间分析)

ArcGIS教程:用矩形提取 (空间分析)
可通过像元的中心来确定该像元是位于矩形的内部还是矩形的外部。如果中心位于矩形轮廓的内部,则即使部分像元落在矩形之外,也会将此像元视为完全处于矩形内部。

摘要

    基于矩形提取栅格像元。

用法

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

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

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

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

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

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

语法

    ExtractByRectangle (in_raster, rectangle, {extraction_area})

参数

说明

数据类型

in_raster

提取像元的输入栅格。

Raster Layer

rectangle

extent

用于定义待提取区域的矩形。可以使用范围对象来指定坐标。

对象形式为:

·         Extent(XMinYMinXMaxYMax)

其中 XMin  YMin 定义待提取区域左下方的坐标,XMax  YMax 定义右上方的坐标。

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

Extent

extraction_area

(可选)

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

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

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

String

返回值

名称

说明

数据类型

out_raster

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

Raster

代码实例

    用矩形提取 (ExtractByRectangle) 示例 1Python 窗口)

    本示例会将矩形范围以外的像元提取到新栅格。

import arcpy

from arcpy import env

from arcpy.sa import *

env.workspace = "C:/sapyexamples/data"

rectExtract = ExtractByRectangle("elevation",

                                 Extent(477625, 213900, 486400, 224200),

                                 "OUTSIDE")

rectExtract.save("c:/sapyexamples/output/extrect")

    用矩形提取 (ExtractByRectangle) 示例 2(独立脚本)

    本示例会将矩形范围以外的像元提取到新栅格。

# Name: ExtractByRectangle_Ex_02.py

# Description:

# 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"

inRectangle = Extent(477625, 213900, 486400, 224200)

# Check out the ArcGIS Spatial Analyst extension license

arcpy.CheckOutExtension("Spatial")

# Execute ExtractByRectangle

rectExtract = ExtractByRectangle(inRaster, inRectangle, "INSIDE")

# Save the output

rectExtract.save("c:/sapyexamples/output/extrect02")

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