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

地理国情监测云平台

据中的 NoData 像元指定 NoData 值。 返回值名称

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

按掩膜提取 (ExtractByMask) 示例 1(Python 窗口)

本示例将从栅格内提取由输入面 shapefile 要素类所定义的掩膜内的像元。

import arcpy

from arcpy import env from arcpy.sa import *

env.workspace = \

outExtractByMask = ExtractByMask(\, \) outExtractByMask.save(\) 按掩膜提取 (ExtractByMask) 示例 2(独立脚本)

本示例将从栅格内提取由输入面 shapefile 要素类所定义的掩膜内的像元。

# Name: ExtractByMask_Ex_02.py

# Description: Extracts the cells of a raster that correspond with the areas # defined by a mask.

# 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 = \ inMaskData = \

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

outExtractByMask = ExtractByMask(inRaster, inMaskData) # Save the output

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

地理国情监测云平台

outExtractByMask.save(\)

四、用点提取

摘要

基于一组坐标点提取栅格像元。 用法

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

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

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

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

ExtractByPoints (in_raster, points, {extraction_area}) 参数 in_raster points [point,...] 说明 提取像元的输入栅格。 点类对象的 Python 列表用于指示要提取栅格值的位置。 点对象均在 x,y 坐标对列表中指定。对象形式为： 数据类型 Raster Layer Point ? [point(x1,y1), point(x2,y2),...] 点所使用的地图单位与输入栅格相同。 extraction_area 标识是基于指定点位置（内部）提取像元还是基于点位置外部（外部）String (可选) 提取像元。 ? INSIDE — 指定将所选点落入的像元写入输出栅格的关键字。方框区域外部的所有像元都将在输出栅格中获得 NoData 值。 ? OUTSIDE — 指定应选择输入点外部的像元并将其写入输出栅格的关键字。 返回值名称 out_raster

说明 包含提取像元值的输出栅格。 数据类型 Raster 北京数字空间科技有限公司

地理国情监测云平台

代码实例

用点提取 (ExtractByPoints) 示例 1(Python 窗口) 本示例将基于指定点坐标提取栅格的像元。

import arcpy

from arcpy import env from arcpy.sa import *

env.workspace = \

pointList = [arcpy.Point(743050, 4321275), arcpy.Point(743100, 4321200), arcpy.Point(734500,4322000)]

outPointExtract = ExtractByPoints(\, pointList,\) outPointExtract.save(\) 用点提取 (ExtractByPoints) 示例 2(独立脚本) 本示例将基于指定点坐标提取栅格的像元。

# Name: ExtractByPoints_Ex_02.py

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

pointList = [arcpy.Point(743050, 4321275), arcpy.Point(743100, 4321200), arcpy.Point(743500, 4322000), arcpy.Point(742900, 4321800)] …… 此处隐藏：704字，全部文档内容请下载后查看。喜欢就下载吧 ……