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

地理国情监测云平台

ArcGIS教程：栅格像元的提取

一、按属性提取

摘要

基于逻辑查询提取栅格像元。 插图

用法

· 如果 Where 子句的求值结果是 True，则将为该像元位置返回初始输入值。如果其求值结果是 False，则将为像元位置指定 NoData。

· 对于输出栅格，将舍弃输入栅格中的所有额外项(值和计数除外)。

· 如果在 Where 子句中指定了除输入栅格的值以外的某一项，则将为此像元位置返回初始输入值。 · 当输入为多波段栅格时，将输出一个新的多波段栅格。对输入多波段栅格中的每一单个波段都会进行相应分析。

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

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

ExtractByAttributes (in_raster, where_clause)

参数 说明 数据类型 北京数字空间科技有限公司

地理国情监测云平台

in_raster 提取像元的输入栅格。 表达式遵循 SQL 表达式的一般格式。 返回值

Raster Layer SQL Expression where_clause 用于选择栅格像元子集的逻辑表达式。

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

按属性提取 (ExtractByAttributes) 示例 1(Python 窗口)

本示例将基于逻辑查询从栅格中提取像元，其中高程大于 1,000 米。

import arcpy

from arcpy import env from arcpy.sa import *

env.workspace = \

attExtract = ExtractByAttributes(\, \) attExtract.save(\)

按属性提取 (ExtractByAttributes) 示例 2(独立脚本)

本示例将基于逻辑查询从栅格中提取像元，其中高程大于 1,000 米。

# Name: ExtractByAttributes_Ex_02.py

# Description: Extracts the cells of a raster based on a logical query. # 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 = \ inSQLClause = \

# Check out the ArcGIS Spatial Analyst extension license arcpy.CheckOutExtension(\) # Execute ExtractByAttributes

attExtract = ExtractByAttributes(inRaster, inSQLClause) # Save the output

attExtract.save(\)

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

地理国情监测云平台

二、按圆形区域提取

摘要

基于圆提取栅格像元。 用法

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

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

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

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

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

参数 in_raster center_point 说明 提取像元的输入栅格。 指示用于定义提取区域的圆的中心坐标 (x,y) 的点类。 该类的形式为： 数据类型 Raster Layer Point ? Point (x, y) 将指定坐标使用与输入栅格相同的地图单位。 radius 用于定义提取区域的圆半径。 将以地图单位指定半径，并且与输入栅格的单位相同。 extraction_area 标识要提取输入圆内部还是输入圆外部的像元。 (可选) Double String ? INSIDE — 指定应选择输入圆内部的像元并将其写入输出栅格的关键字。圆形区域外部的所有像元都将在输出栅格中获得 NoData 值。 ? OUTSIDE — 指定应选择输入圆外部的像元并将其写入输出栅格的关键字。圆形区域内部的所有像元都将在输出栅格中获得 NoData 值。 返回值

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

地理国情监测云平台

名称 out_raster 代码实例

按圆形区域提取 (ExtractByCircle) 示例 1(Python 窗口) 本示例将在某一点位置方圆 500 米的半径范围内提取像元。

说明 包含提取像元值的输出栅格。 数据类型 Raster import arcpy

from arcpy import env from arcpy.sa import *

env.workspace = \

outExtCircle = ExtractByCircle(\, arcpy.Point(482838.823, 222128.982), 500, \) outExtCircle.save(\) 按圆形区域提取 (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 = \ # Set local variables inRaster = (\)

centerPoint = arcpy.Point(482838.823, 222128.982) circRadius = 1000 extractType = \

# Check out the ArcGIS Spatial Analyst extension license arcpy.CheckOutExtension(\) # Execute ExtractByCircle

outExtCircle = ExtractByCircle(inRaster, centerPoint, circRadius, extractType) # Save the output

outExtCircle.save(\)

三、按掩膜提取

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

地理国情监测云平台

摘要

提取所定义掩膜区域内的栅格像元。 插图

用法

· 按掩膜提取工具与设置掩膜环境大体相同，但二者的区别在于输入掩膜只适用于即时实例，而在环境中设置的掩膜可在发生更改或被禁用之前应用于所有工具。

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

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

· 如果输入内容为根据多波段栅格(超过三个波段)创建的图层，提取操作将只考虑图层加载(符号化)的波段。因此，输出的多波段栅格只会具有三个波段，分别对应于输入图层显示中使用的多波 …… 此处隐藏：2441字，全部文档内容请下载后查看。喜欢就下载吧 ……