ArcGIS 10 栅格数据批量值提取到点

ArcGIS 10 栅格数据批量值提取到点

在网上找了好久都没有找到批量值提取到点的代码，自己试着改了改别人的代码，也问了问大神，终于弄出来了，现在把代码分享给大家：

import arcpy
from arcpy import env
from arcpy.sa import *

arcpy.CheckOutExtension("spatial")

arcpy.gp.overwriteOutput=1 #检查许可

arcpy.env.workspace = "G:\\YR_R\\pre_month"#输入文件夹

rasters = arcpy.ListRasters("*", "tif")#提取列表中的.tif文件

mask= "G:\\YR_R\\verify\\ALT_ver.shp"#输入裁剪区域

out= "G:\\YR_R\\pre_point\\"#输出名称和路径

for raster in rasters:#循环栅格

    print(raster)#输出栅格
    out_put_shp = out + raster[:-4]+".shp"
    

    ExtractValuesToPoints(mask, raster, out_put_shp,"INTERPOLATE","VALUE_ONLY")#值提取到点工具

    print("has done")#输出提醒

print("!")#文件中数据完成提醒