[2018年最新整理]ArcGIS教程：提取栅格像元几种方式.docxVIP

ArcGIS教程：栅格像元的提取一、按属性提取摘要　基于逻辑查询提取栅格像元。插图用法　· 如果 Where 子句的求值结果是 True，则将为该像元位置返回初始输入值。如果其求值结果是 False，则将为像元位置指定 NoData。　· 对于输出栅格，将舍弃输入栅格中的所有额外项(值和计数除外)。　· 如果在 Where 子句中指定了除输入栅格的值以外的某一项，则将为此像元位置返回初始输入值。　· 当输入为多波段栅格时，将输出一个新的多波段栅格。对输入多波段栅格中的每一单个波段都会进行相应分析。　默认输出格式是 Esri Grid 堆栈。请注意，Esri Grid 堆栈的名称不能以数字开头、不能使用空格，并且长度不能超过 9个字符。　· 如果输入内容为根据多波段栅格(超过三个波段)创建的图层，提取操作将只考虑图层加载(符号化)的波段。因此，输出的多波段栅格只会具有三个波段，分别对应于输入图层显示中使用的多波段栅格。　· 如果输入栅格为整型，那么输出栅格也为整型。如果输入栅格为浮点型，则输出栅格也为浮点型。语法　ExtractByAttributes (in_raster, where_clause)参数说明数据类型in_raster提取像元的输入栅格。Raster Layerwhere_clause用于选择栅格像元子集的逻辑表达式。表达式遵循?SQL?表达式的一般格式。SQL Expression返回值名称说明数据类型out_raster包含提取像元值的输出栅格。Raster代码实例　按属性提取 (ExtractByAttributes) 示例 1(Python 窗口)　本示例将基于逻辑查询从栅格中提取像元，其中高程大于 1,000 米。import?arcpyfrom?arcpy?import?envfrom?arcpy.sa?import?*env.workspace =?C:/sapyexamples/dataattExtract = ExtractByAttributes(elevation,?VALUE 1000)attExtract.save(c:/sapyexamples/output/attextract) 按属性提取 (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 modulesimport?arcpyfrom?arcpy?import?envfrom?arcpy.sa?import?*# Set environment settingsenv.workspace =?C:/sapyexamples/data# Set local variablesinRaster =?elevationinSQLClause =?VALUE 1000# Check out the ArcGIS Spatial Analyst extension licensearcpy.CheckOutExtension(Spatial)# Execute ExtractByAttributesattExtract = ExtractByAttributes(inRaster, inSQLClause)# Save the outputattExtract.save(c:/sapyexamples/output/attextract02)二、按圆形区域提取摘要　基于圆提取栅格像元。用法　· 可通过像元的中心来确定该像元是位于圆的内部还是位于圆的外部。如果中心位于圆弧的内部，则即使部分像元落在圆外，也会将此像元视为完全处于圆内。　· 未选择的像元位置被赋予 NoData 值。　· 当输入为多波段栅格时，将输出一个新的多波段栅格。对输入多波段栅格中的每一单个波段都会进行相应分析。　默认输出格式是 Esri Grid 堆栈。请注意，Esri Grid 堆栈的名称不能以数字开头、不能使用空格，并且长度不能超过 9个字符。　· 如果输入内容为根据多波段栅格(超过三个波段)创建的图层，提取操作将只考虑图层加载(符号化)的波段。因此，输出的多波段栅格只会具有三个波段，分别对应于输入图层显示中使用的多波段栅格。　· 如果输入栅格为整型，那么输出栅格也为整型。如果输入栅格为浮点型，则输出栅格也为浮点型。语法　ExtractByCircle (i