ArcGIS：Extract Values to Points的ArcPy实践

Extract Values to Points又叫“值提取至点”，是地理信息数据处理中的一种常用方法，本文主要介绍了Extract Values to Points如何通过调用ArcPy模块在Python中进行实现。

用途

• 基于一组点要素，提取点对应的栅格像元值，并将这些值记录到输出点要素类的属性表中。

语法

• ExtractValuesToPoints (in_point_features, in_raster, out_point_features, {interpolate_values}, {add_attributes})

代码

需安装有ArcGIS，且只能用Python2进行运行（安装ArcGIS时一般都会默认安装Python2.7）。

# 导入相关的库
import arcpy
from arcpy import env
from arcpy.sa import *

# 设置工作空间（也就是输入矢量点数据和栅格数据的文件位置）
env.workspace = "C:/sapyexamples/data"

# 设置局部变量
inPointFeatures = "rec_sites.shp"  # 输入矢量点数据
inRaster = "elevation"  # 待提取的栅格数据（也可以是elevation.tif等形式）
outPointFeatures = "C:/sapyexamples/output/extractvaluespts.shp"  # 输出矢量点数据路径

# 检查ArcGIS Spatial Analyst扩展许可证
arcpy.CheckOutExtension("Spatial")

# 执行Extract Values To Points
# 设置interpolate_values（指定是否使用插值）
# NONE — 不应用任何插值法；将使用像元中心值。这是默认设置。
# INTERPOLATE — 将使用双线性插值法根据相邻像元的有效值计算像元值。除非所有相邻像元都为NoData，否则会在插值时忽略NoData值。
# 设置add_attributes（确定是否将栅格属性写入输出点要素数据集）
# VALUE_ONLY — 仅将输入栅格的值添加到点属性。这是默认设置。
# ALL — 输入栅格的所有字段（“计数”除外）都将添加到点属性。
ExtractValuesToPoints(inPointFeatures, inRaster, outPointFeatures, "INTERPOLATE", "VALUE_ONLY")

注意事项

• 输入点要素类中的所有字段均将包括在输出点要素类中。
• 输出要素类将添加一个名为 RASTERVALU 的新字段。

该字段不能已存在于输入要素的属性表中。否则，将发生错误。如果要保留原始信息，可在执行值提取至点之前向属性表添加一个新字段（例如，RASVAL1），计算该字段的值，然后删除原始的 RASTERVALU 字段。

• 对多波段栅格使用值提取至点时，RASTERVALU 字段将包含来自输入栅格的最后一个波段的值。要从多个栅格或某个多波段栅格数据集提取值，请使用多值提取至点工具。
• 对于属性表的 RASTERVALU 字段，若赋值栅格中的像元值为 NoData 像元，则该字段的值为 -9999。
• 插值选项决定了从栅格中获取值的方式。默认选项为使用被采样像元的中心值。插值选项将使用双线性插值法在像元中心内插一个值。
• 如果输入栅格为浮点型，则生成的输出点数据集将仅包含输入要素数据的属性和像元值（由插值选项确定）。
• 不支持将多点数据集作为输入。

最后

• 内容为本人在学习ArcGIS时的个人总结，若有不足之处，敬请大家批评指正！

参考资料

• Extract Values to Points帮助文档：https://desktop.arcgis.com/en/arcmap/10.4/tools/spatial-analyst-toolbox/extract-values-to-points.htm

最后

以上就是从容手机为你收集整理的ArcGIS：Extract Values to Points的ArcPy实践的全部内容，希望文章能够帮你解决ArcGIS：Extract Values to Points的ArcPy实践所遇到的程序开发问题。

如果觉得靠谱客网站的内容还不错，欢迎将靠谱客网站推荐给程序员好友。