import arcpy
from arcpy.sa import *
arcpy.CheckOutExtension("spatial")
arcpy.gp.overwriteOutput=1
#custom
arcpy.env.workspace="G:\\Phenology of 30 Years\\GIMMS 3g\\15Length\\1Length\\"
#custom
outpath="G:\\Phenology of 30 Years\\GIMMS 3g\\15Length\\2mean_len\\"
#custom
outfilename="mean";n=30;
Sum=0
files=arcpy.ListRasters()
for file in files:
Sum=Sum+Raster(file)
(Sum/n).save(outpath+outfilename)
print("Done,please close")
摘抄至互联网。
n指栅格个数。要是n可以通过遍历自动获取就更好了,有时间改进一下