用Python实现.shp边界批量剪裁TIF影像，处理出来有一定的误差

因为.nc原影像转为tiff格式也是代码处理的，所以剪裁后里面的数据有一定的误差。目前没有找到相关原因

以下代码可以用ARCGIS自带的GUI编译器编译，也可以用ARCGIS PRO自带的python编译，也是在网上的找的相关代码进行修改的，源码现在是找不到了，网上有很多，你们再查找的过程中应该都已经查看过。

以下是可实现代码案例：

import arcpy
import os

arcpy.CheckOutExtension("spatial")

arcpy.env.workspace = "H:/LAB/LABDATA/ECMWF/3/TIF/si10" # 工作目录

rasters = arcpy.ListRasters("*", "TIF") # 获取所有的TIF影像

mask= "H:/LAB/LABDATA/landsat/shp/caohai2013.shp" # 获取掩膜边界

for raster in rasters:

    print(raster)

    out= "H:/LAB/LABDATA/ECMWF/3/TIF/caohai/si10/"+"clip_"+raster # 输出路径及文件名

    arcpy.Clip_management(raster,"#", out, mask, "0", "ClippingGeometry") # arcpy调用ARCMAP里的Clip工具

    print("clip_"+raster+"  has done")

print("All done")

运行结果：

结果仅供参考，如果中间不报错，有上面这些信息，应该就是正确的！