ArcGIS教程：提取栅格像元几种方式(4)

地理国情监测云平台

arcpy.Point(743500, 4322000),arcpy.Point(742900, 4321800)] # Check out the ArcGIS Spatial Analyst extension license arcpy.CheckOutExtension(\) # Execute ExtractByPolygon

extPolygonOut = ExtractByPolygon(inRaster, polyPoints, \) # Save the output

extPolygonOut.save(\)

六、用矩形提取

摘要

基于矩形提取栅格像元。 用法

· 可通过像元的中心来确定该像元是位于矩形的内部还是矩形的外部。如果中心位于矩形轮廓的内部，则即使部分像元落在矩形之外，也会将此像元视为完全处于矩形内部。 · 未选择的像元位置被赋予 NoData 值。

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

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

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

ExtractByRectangle (in_raster, rectangle, {extraction_area}) 参数 in_raster rectangle extent 说明 提取像元的输入栅格。 用于定义待提取区域的矩形。可以使用范围对象来指定坐标。 对象形式为： 数据类型 Raster Layer Extent ? Extent(XMin，YMin，XMax，YMax) 其中 XMin 和 YMin 定义待提取区域左下方的坐标，XMax 和 YMax 定义右上方的坐标。 将指定坐标使用与输入栅格相同的地图单位。 extraction_area 标识要提取输入矩形内部还是输入矩形外部的像元。 北京数字空间科技有限公司

String 地理国情监测云平台

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

说明 包含提取像元值的输出栅格。 数据类型 Raster 名称 out_raster 代码实例

用矩形提取 (ExtractByRectangle) 示例 1(Python 窗口) 本示例会将矩形范围以外的像元提取到新栅格。

import arcpy

from arcpy import env from arcpy.sa import *

env.workspace = \ rectExtract = ExtractByRectangle(\,

Extent(477625, 213900, 486400, 224200), \) rectExtract.save(\) 用矩形提取 (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 = \ # Set local variables inRaster = \

inRectangle = Extent(477625, 213900, 486400, 224200) # Check out the ArcGIS Spatial Analyst extension license

北京数字空间科技有限公司

地理国情监测云平台

arcpy.CheckOutExtension(\) # Execute ExtractByRectangle

rectExtract = ExtractByRectangle(inRaster, inRectangle, \) # Save the output

rectExtract.save(\)

北京数字空间科技有限公司