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

MeteoInfo

一股夏天的风 发表于 2017-3-30 19:46
谢谢老师指点，查看学习了Python的语言，修改了语句，再次运行时错误如下
Traceback (most recent call  ...

addfiles函数要求文件中必须有时间维，也许你的数据文件缺少时间维。

MeteoInfo

一股夏天的风 发表于 2017-3-30 19:46
谢谢老师指点，查看学习了Python的语言，修改了语句，再次运行时错误如下
Traceback (most recent call  ...

可以不用addfiles，利用循环来做平均：

1. datadir = 'F:/OMI'
   
2. st = datetime.datetime(2017,1,1)
   
3. et = datetime.datetime(2017,1,2)
   
4. vname = 'ColumnAmountO3'
   
5. ave = None
   
6. n = 0
   
7. while st <= et:
   
8.     fn = os.path.join(datadir, ' OMI-Aura_L3-OMTO3e_' + st.strftime('%Y%m%d') + '.he5')
   
9.     f = addfile(fn)
   
10.     data = f[vname]['32:55', '95:128']
    
11.     if ave is None:
    
12.         ave = data
    
13.     else:
    
14.         ave = ave + data
    
15.     n += 1
    
16.     st = st + datetime.timedelta(days=1)
    
17. 
    
18. data = ave / n
    
19. 
    
20. #Plot
    
21. axesm()
    
22. lhuabei=shaperead('D:/BaiduNetdiskDownload/MeteoInfo/MeteoInfo_Java_1.4.2_Files/MeteoInfo/map/huabei.shp ')
    
23. geoshow(lhuabei, edgecolor='k')
    
24. layer = imshowm(data)
    
25. colorbar(layer)

datadir = 'F:/OMI'<br /> st = datetime.datetime(2017,1,1)<br /> et = datetime.datetime(2017,1,2)<br /> vname = 'ColumnAmountO3'<br /> ave = None<br /> n = 0<br /> while st <= et:<br />     fn = os.path.join(datadir, ' OMI-Aura_L3-OMTO3e_' + st.strftime('%Y%m%d') + '.he5')<br />     f = addfile(fn)<br />     data = f[vname]['32:55', '95:128']<br />     if ave is None:<br />         ave = data<br />     else:<br />         ave = ave + data<br />     n += 1<br />     st = st + datetime.timedelta(days=1)<br /> <br /> data = ave / n<br /> <br /> #Plot<br /> axesm()<br /> lhuabei=shaperead('D:/BaiduNetdiskDownload/MeteoInfo/MeteoInfo_Java_1.4.2_Files/MeteoInfo/map/huabei.shp ')<br /> geoshow(lhuabei, edgecolor='k')<br /> layer = imshowm(data)<br /> colorbar(layer)

一股夏天的风

MeteoInfo 发表于 2017-3-30 22:51
可以不用addfiles，利用循环来做平均：

王老师，非常感谢您的帮助，我试了您给的程序，运行错误如下；
> run script...
  File "<iostream>", line 8
    &Acirc;&nbsp; &Acirc;&nbsp; fn = os.path.join(datadir, 'OMI-Aura_L3-OMTO3e_' + st.strftime('%Y%m%d')+ '.he5')
    ^
SyntaxError: no viable alternative at character '&nbsp;'
我想大概是符号格式的问题，重写引号还是同样问题，是不是he5数据本身的问题？百思不得其解，请您给看看

MeteoInfo

一股夏天的风 发表于 2017-3-31 09:39
王老师，非常感谢您的帮助，我试了您给的程序，运行错误如下；
> run script...
  File "", line 8

你不要把格式符号也复制成代码了呀！

一股夏天的风

MeteoInfo 发表于 2017-3-31 09:43
你不要把格式符号也复制成代码了呀！

已经运行成功了，非常感谢您的耐心指导，自己的编程能力需要提高

一股夏天的风

MeteoInfo 发表于 2017-3-31 09:43
你不要把格式符号也复制成代码了呀！

王老师，如果我想统计数值中的最大和最小值，均值，应该用什么函数？在论坛里好像没有找到相关帖子，谢谢您

MeteoInfo

一股夏天的风 发表于 2017-3-31 21:08
王老师，如果我想统计数值中的最大和最小值，均值，应该用什么函数？在论坛里好像没有找到相关帖子，谢谢 ...

获取的数组有 max() 和 min() 函数来获取数组的最大最小值。比如：

max = data.max()
min = data.min()

MeteoInfo

一股夏天的风 发表于 2017-3-31 21:08
王老师，如果我想统计数值中的最大和最小值，均值，应该用什么函数？在论坛里好像没有找到相关帖子，谢谢 ...

12楼之前的脚本有点问题：

18行代码： data = data / n    改为：    data = ave / n

一股夏天的风

MeteoInfo 发表于 2017-3-31 22:31
12楼之前的脚本有点问题：

18行代码： data = data / n    改为：    data = ave / n

王老师，如何输出某一确定经纬度的数值？应该用什么语句？

MeteoInfo

一股夏天的风 发表于 2017-4-12 21:56
王老师，如何输出某一确定经纬度的数值？应该用什么语句？

类似这样，获取经度100，纬度40的数据：
v =data['40', '100']