如何一次读取一个月的OMI Grid HDF 文件，并求均值

一股夏天的风

初次接触meteoinfo和meteoinfolab，学习了两个教程，编写代码实现了OMI臭氧柱总量的中国地区的图像显示，感觉软件上手较快，教程也很详尽，受益匪浅。
浏览了论坛的所有文章，但是还有一个问题没有解决，请大家帮忙。
1、我仅处理了一个HDF-EOS5文件，见下图和代码。如何写代码一次读一个月的数据，然后计算月平均值，并图像输出呢？

folder = 'F:/OMI/'
fns = 'OMI-Aura_L3-OMDOAO3e_2005m0101_v003-2011m1109t085643.he5'
fn = folder + fns
f = addfile(fn)
axesm()
m_china = shaperead('D:/BaiduNetdiskDownload/MeteoInfo/MeteoInfo_Java_1.4.2_Files/MeteoInfo/map/huabei.shp')
geoshow(m_china, linecolor=(0,0,255))
vname = 'ColumnAmountO3'
v = f[vname]
data = v[:,:]
layer = imshowm(data, 20)
masklayer(m_china, [layer])
title('OMI - ' + vname)
xlim(95,128k)
ylim(32,55)
ylabel('Latitude')
xlabel('Longitude')
colorbar(layer)
show()

︶ㄣ安定■丶

可以看看addfiles函数

︶ㄣ安定■丶

http://bbs.06climate.com/forum.p ... &extra=page%3D1

fabula

谢谢！很有用

MeteoInfo

show() 是多余的，可以去掉。

可以用经纬度限定数组的范围会算短运算时间，比如：data = v['32:55', '95:128']

可以用循环读取多个文件并取出数组，对数组进行累加，然后除以累加的次数就是平均。数组可以像一个标量一样进行各种运算。

一股夏天的风

datadir = 'F:/OMI'
st = datetime.datetime(2017,1,1)
et = datetime.datetime(2017,1,2)
fns = []
while st <= et:
    fn = os.path.join(datadir, ' OMI-Aura_L3-OMTO3e_' + st.strftime('%Y%m%d').he5')
    fns.append(fn)
    st = st + datetime.timedelta(hours=24)
fs = addfiles(fns)
vname = 'ColumnAmountO3'
v = f[vname]
data = v['32:55', '95:128']
ndata = mean(data, axis=0)
#Plot
axesm()
lhuabei=shaperead('D:/BaiduNetdiskDownload/MeteoInfo/MeteoInfo_Java_1.4.2_Files/MeteoInfo/map/huabei.shp ')
geoshow(lhuabei, edgecolor='k')
layer = imshowm(ndata)
colorbar(layer)
运行出错：
SyntaxError: no viable alternative at input 'fns'
完全是照着教程修改的，不知道应该怎么改

MeteoInfo

一股夏天的风 发表于 2017-3-29 21:06
datadir = 'F:/OMI'
st = datetime.datetime(2017,1,1)
et = datetime.datetime(2017,1,2)

可能是代码书写格式的问题吧，你重新将代码些以便试试

一股夏天的风

MeteoInfo 发表于 2017-3-29 22:30
可能是代码书写格式的问题吧，你重新将代码些以便试试

重新写，还是遇到同样的问题，SyntaxError: no viable alternative at input 'fns'

MeteoInfo

一股夏天的风 发表于 2017-3-30 11:24
重新写，还是遇到同样的问题，SyntaxError: no viable alternative at input 'fns'

fn = os.path.join(datadir, ' OMI-Aura_L3-OMTO3e_' + st.strftime('%Y%m%d').he5')
这个语句有问题，你自己调试一下吧

一股夏天的风

MeteoInfo 发表于 2017-3-30 12:07
fn = os.path.join(datadir, ' OMI-Aura_L3-OMTO3e_' + st.strftime('%Y%m%d').he5')
这个语句有问题， ...

谢谢老师指点，查看学习了Python的语言，修改了语句，再次运行时错误如下
Traceback (most recent call last):
  File "<iostream>", line 10, in <module>
  File "D:\BaiduNetdiskDownload\MeteoInfo\MeteoInfo_Java_1.4.2_Files\MeteoInfo\pylib\mipylib\dataset\midata.py", line 68, in addfiles
    return DimDataFiles(dfs)
  File "D:\BaiduNetdiskDownload\MeteoInfo\MeteoInfo_Java_1.4.2_Files\MeteoInfo\pylib\mipylib\dataset\dimdatafile.py", line 352, in __init__
    ftimes.append(ds.gettime(0))
  File "D:\BaiduNetdiskDownload\MeteoInfo\MeteoInfo_Java_1.4.2_Files\MeteoInfo\pylib\mipylib\dataset\dimdatafile.py", line 171, in gettime
    t = self.dataset.getDataInfo().getTimes().get(idx)     
        at org.meteoinfo.data.meteodata.DataInfo.getTimes(DataInfo.java:154)
        at sun.reflect.NativeMethodAccessorImpl.invoke0(Native Method)
        at sun.reflect.NativeMethodAccessorImpl.invoke(Unknown Source)
        at sun.reflect.DelegatingMethodAccessorImpl.invoke(Unknown Source)
        at java.lang.reflect.Method.invoke(Unknown Source)
请老师指点，应该是我addfiles用的不对