import arcpy
from arcpy import env
from arcpy.sa import *
import os
env.workspace = "E:\ArcGIS\Projects\ET0\ET0.gdb" #设置环境
arcpy.env.cellSize = "MAXOF"#设置像元大小
env.extent = E:\ArcGIS\Projects\cailiao\AICGIS\zhongyao\mian.shp" #设置范围
env.mask = "E:\ArcGIS\Projects\cailiao\AICGIS\zhongyao\mian.shp" #设置掩膜
inPointFeatures = "E:\ArcGIS\Projects\ET0\ET0.gdb\ET0_year" #点要素
zField = "1965" #z字段
for i in range(50):
j = zField + i
outIDW = Idw(inPointFeatures,str(j))
outIDW.save(r"C:\Users\WH\Desktop\ET0\year\ET0\ET0"+str(j)+".tif")