python_进行栅格计算_加法_批量处理

import arcpy  
from arcpy.sa import *  
arcpy.CheckOutExtension("spatial")  
arcpy.gp.overwriteOutput=1   
#定义输出路径 
outpath="e:/ccc/"   
outfilename="mean"
#更改工作空间将，将"e:/bbb"更改为处理cdf文件时的输出路径
arcpy.env.workspace = "e:/bbb"
files=arcpy.ListFiles()  
for file in files:
#将"e:/bbb"更改为处理cdf文件时的输出路径
    wk="e:/bbb/"+file
    arcpy.env.workspace = wk
    rasters = arcpy.ListRasters()
    Sum1=0 
    for raster in rasters:  
       Sum1=Sum1+Raster(raster)  
#也可以用outCellStatistics = CellStatistics(list_raster, "SUM", "DATA") 来计算平均或者求和
       Sum1.save(outpath+outfilename+file)
    print file

##print("Done,please close")  

按年为单位的需要计算每年的平均值，也就是说要计算每一年每个格网12个月的平均值。ArcGIS通过Cell Statistics可以实现，这里使用python进行批量处理