如何用meteoinfolab读取一个grib2

一股夏天的风

如何用meteoinfolab读取一个grib2里面的气压、温度、风场等各气象要素？并且用图表示呢？

MeteoInfo

这里有很多例子可以参考：http://bbs.06climate.com/forum.p ... &extra=page%3D1

另外建议仔细看看MeteoInfoLab的User guide （http://www.meteothinker.com/docs/meteoinfolab/user_guide.html），能极大的帮助理解MeteoInfoLab。

一股夏天的风

MeteoInfo 发表于 2017-4-17 23:45
这里有很多例子可以参考：http://bbs.06climate.com/forum.php?mod=viewthread&tid=36151&extra=page%3D1
...

[/img]

1. <p>参考一些例子后做出了850hpa的等温线，想叠加风场和高度场，单独读写风场，如下</p><p>fn = 'F:/GRIB/fnl_20121005_12_00.grib2'
   
2. f = addfile(fn)
   
3. u = f['U'][0,'500','15:55','70:140']
   
4. v = f['V'][0,'500','15:55','70:140']
   
5. speed = sqrt(u*u+v*v)
   
6. #Plot
   
7. axesm()
   
8. mlayer  = shaperead('D:/BaiduNetdiskDownload/MeteoInfo/MeteoInfo_Java_1.4.2_Files/MeteoInfo/map/country1.shp')
   
9. geoshow(mlayer)
   
10. layer = quiverm(u, v, speed, 10, size=8)
    
11. t = f.gettime(0)
    
12. title(' Wind field 500hpa (' + t.strftime('%Y-%m-%d %H:00') + ')')
    
13. colorbar(layer) </p><p>运行出错</p><p>Traceback (most recent call last):
    
14.   File "<iostream>", line 3, in <module>
    
15.   File "D:\BaiduNetdiskDownload\MeteoInfo\MeteoInfo_Java_1.4.2_Files\MeteoInfo\pylib\mipylib\dataset\dimvariable.py", line 72, in __getitem__
    
16.     if len(indices) != self.ndim:
    
17. AttributeError: DimVariable instance has no attribute 'ndim'</p><p>王老师，这条语句不对，u = f['U'][0,'500','15:55','70:140']，应该怎么改呢？</p>

复制代码

MeteoInfo

一股夏天的风 发表于 2017-4-18 16:19
[/img]

你首先需要了解数据里面包含什么变量，每个变量的维设置是什么。仔细看看这里：http://www.meteothinker.com/docs ... /data_tutorial.html

浪子

谢谢楼主分享

一股夏天的风

MeteoInfo 发表于 2017-4-18 17:16
你首先需要了解数据里面包含什么变量，每个变量的维设置是什么。仔细看看这里：http://www.meteothinker. ...

嗯，帮助里写
float U(T,Z,Y,X);
        U: description = "U"
float V(T,Z,Y,X);
        V: description = "V"
意思是风的两个分量表示方法，所以我参照写了下面的句子
u = f['U'][0,'500','15:55','70:140']

v = f['V'][0,'500','15:55','70:140']
，还是运行出错，估计是我的数据里面不包含风场，又试了气压场也不对，只有温度可以，心想不应该啊

MeteoInfo

一股夏天的风 发表于 2017-4-19 09:34
嗯，帮助里写
float U(T,Z,Y,X);
        U: description = "U"

在Console里输入  

>>> f

然后回车，看看文件里都有哪些变量，变量名称是什么