栅格批量求平均（栅格镶嵌，不受缺失像元影响）

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

arcpy.env.workspace=r"I:\硕士论文数据\GPP_MOD17A2HGF裁剪_去掉Nodata"
outpath=unicode("I:/硕士论文数据/GPP_MOD17A2HGF年际平均/",'utf-8')

for y in range(2000, 2020):
    y=str(y)
    files = arcpy.ListRasters("*A"+y+"*Gpp_500m.tif")
    outfilename = "MOD17A2HGF_year_"+y+".tif"
    rast = ""
    for file in files:
        rast =rast + ";" + file
    arcpy.MosaicToNewRaster_management(rast, outpath, outfilename, "", "32_BIT_UNSIGNED", "", "1", "MEAN", "FIRST")

print("Done,please close")