arcpy 实现将值提取到点（批量）

奇丶迹

已于 2023-03-20 09:30:44 修改

阅读量2.4k

点赞数

分类专栏： arcpy小工具文章标签： arcgis

于 2022-10-06 22:05:47 首次发布

本文链接：https://blog.csdn.net/qq_42940285/article/details/127188637

版权

现有一dem 以及多个shp点图层，需要将DEM的高程值添加到点图层中。若单个作业，则需要用到arcgis中的值提取至点工具在空间分析模块-提取分析-值提取至点，分别输入点要素与栅格图层即可。

因要提取的图层较多，故需进行批量处理。

代码如下：

# -*- coding: utf-8 -*-
# @Time  :2022/10/6 20:18
# @Author:奇迹

import arcpy
from arcpy import env
from arcpy.sa import *           #取出空间分析工具中的所有工具
from glob import glob            #方便遍历shp文件
import sys
reload(sys)
sys.setdefaultencoding('utf-8')


path = r"E:\shp"                
list = glob(path+r'\xian\*.shp')    #遍历该目录下的所有shp文件
env.workspace = path
input_tif = path + r'\dem\eastsouth.tif'   #输入的栅格图层
output = path + r'\xianres\{name}'           #输出的路径 

for lis in list:
    output = path + r'\xianres\{name}'         #提取输出文件的文件名
    prefix =