平均10年的数据与Python

Question

日降水量pcraster地图栈我在为日降水量值pcraster地图堆型3652地图

pptn0000.001 
pptn0000.002 
... 
pptn0003.652 

据我了解，有一个pcrcalc命令，但我缺乏的例子在这种类型的地图上进行数学运算。 我如何在python中对这些地图进行平均并将一个栅格地图作为输出？ 谢谢

Answer 1

那么，我刚刚做了这个以获得TRMM使命（17年，1998 - 2014年）的'气候降水'。 Python可以直接使用pcraster库读取栅格地图。例如：

import os 
    import fnmatch 
    import subprocess 
    from os import listdir 
    from pcraster import * 

    PrecDaily = '/home/someData/precipitation/pptn' 
    OutSave = '/home/someData/precipitation/pptn/meanMap' 
    TypeFile = 'pptn????.???' 

    os.chdir(PrecDaily) 
    PrecDailyFiles = [] 
    for iListFile in sorted(os.listdir('.')): 
     if fnmatch.fnmatch(iListFile, TypeFile): 
      PrecDailyFiles.append(iListFile) 
    n=len(PrecDailyFiles) 

    FirstRasterMap = readmap(PrecDailyFiles[0]) 
    SUM = FirstRasterMap 
    for iMeanRaster in xrange(1, n): 
     PCRasterMap = readmap(PrecDailyFiles[iMeanRaster]) 
     SUM = SUM + PCRasterMap 
    MEAN = SUM/n 

    os.chdir(OutSave) 
    report(MEAN, 'pptn_mean.map')