MeteoInfoLab脚本示例：沿特定维平均

MeteoInfo

MeteoInfoLab的minum模块中有mean()函数来计算多维数组的平均值：
mean(x, axis=None)

参数x是欲进行平均计算的数组，axis参数指定沿某个特定坐标维进行平均。
>>>  a = array([[1, 2], [3, 4]])
>>> mean(a)    #axis缺省为None，结果为数组所有元素的平均值
2.5
>>> mean(a, axis=0)    #数组沿第一维进行平均
array([2.0, 3.0])
>>> mean(a, axis=1)    #数组沿第二维进行平均
array([1.5, 3.5])


再看看实际多维气象数据平均的例子，打开一个GrADS格式数据文件，从3维变量PS中获取包含4个时次（0:3）的3维数组，然后沿时间维进行平均，最后绘制平均结果图形：

1. fn = 'D:/Temp/GrADS/model.ctl'f = addfile(fn)
   
2. data = f['PS'][0:3,:,:]
   
3. ave = mean(data, axis=0)
   
4. #Plot
   
5. axesm()
   
6. lworld = shaperead('D:/Temp/map/country1.shp')
   
7. geoshow(lworld)
   
8. layer = contourfm(ave)
   
9. colorbar(layer)
   
10. st = f.gettime(0)
    
11. et = f.gettime(3)
    
12. title('Mean surface pressure (' + st.strftime('%Y%m%d%H') + ' - ' + \
    
13.     et.strftime('%Y%m%d%H') + ')')

fn = 'D:/Temp/GrADS/model.ctl'f = addfile(fn)<br /> data = f['PS'][0:3,:,:]<br /> ave = mean(data, axis=0)<br /> #Plot<br /> axesm()<br /> lworld = shaperead('D:/Temp/map/country1.shp')<br /> geoshow(lworld)<br /> layer = contourfm(ave)<br /> colorbar(layer)<br /> st = f.gettime(0)<br /> et = f.gettime(3)<br /> title('Mean surface pressure (' + st.strftime('%Y%m%d%H') + ' - ' + \<br />     et.strftime('%Y%m%d%H') + ')')

Linda.W

感谢王老师

huizhch

我怎么觉得mean函数当axis=0时，是对各列求均值。