栅格数据季节求和（需根据数据数量进行调整）

#coding:utf-8
import arcpy
from arcpy.sa import *
arcpy.CheckOutExtension("spatial")
arcpy.gp.overwriteOutput=1

arcpy.env.workspace=r"I:\硕士论文数据\GPP_MOD17A2HGF裁剪_去掉Nodata"
outpath1=unicode("I:/硕士论文数据/GPP_MOD17A2HGF季节合成_1/spring/",'utf-8')
outpath2=unicode("I:/硕士论文数据/GPP_MOD17A2HGF季节合成_1/summer/",'utf-8')
outpath3=unicode("I:/硕士论文数据/GPP_MOD17A2HGF季节合成_1/fall/",'utf-8')

for y in range(2000, 2020):
    y=str(y)
    SpringFiles = arcpy.ListRasters("*A"+y+"*Gpp_500m.tif")
    outfilename1 = y+"_spring"+"_MOD17A2GF_GPP.tif"
    Sum1 = 0
    for i in range(7,19):
        file=SpringFiles[i]
        Sum1 = Sum1 + Raster(file)
    (Sum1).save(outpath1 + outfilename1)

    outfilename2 = y+"_summer"+"_MOD17A2GF_GPP.tif"
    Sum2 = 0
    for j in range(19,30):
        file=SpringFiles[j]
        Sum2 = Sum2 + Raster(file)
    (Sum2).save(outpath2 + outfilename2)

    outfilename3 = y+"_fall"+"_MOD17A2GF_GPP.tif"
    Sum3 = 0
    for k in range(30,42):
        file=SpringFiles[k]
        Sum3 = Sum3 + Raster(file)
    (Sum3).save(outpath3 + outfilename3)

print("Done,please close")