python实现栅格计算器功能

对于Python不能直接调用arcgis中的栅格计算器工具。但是可以通过使用Raster函数实现。

如下面代码实现了每年12个月做一个累加，然后输出一个年值的功能：

for year in range(1960,1965):
    print year
    temp=0     #必须提前定义要不报错
    for month in range(1,13):
        print month
        temp=temp+Raster(str(year)+"_"+str(month)+".TIF")   #实现累加功能，通过Raster函数直接接受一个栅格数据的地址即可。
    temp.save("E:/aa/"+str(year)+"code.TIF")
print "OK"