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()
    Sum=0 
    for raster in rasters:  
       Sum=Sum+Raster(raster)  
       (Sum/len(rasters)).save(outpath+outfilename+file)
    print file
##  
##print("Done,please close")  
##files=arcpy.ListRasters()  
##  
##for file in files:  
##    Sum=Sum+Raster(file)  
##(Sum/n).save(outpath+outfilename)  
##  
##print("Done,please close")  

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

转载自:https://blog.csdn.net/u012131430/article/details/81130121

You may also like...