MeteoInfoLab脚本示例：GRIB2雷达反射率数据

MeteoInfo

示例读取GRIB2格式的雷达反射率数据并绘图。

脚本程序：

1. f = addfile('D:/Temp/grib/1kmradar_201512241600.grib2')
   
2. data = f['Base_reflectivity_surface_layer'][0,:,:]
   
3. t = f.gettime(0)
   
4. axesm(projinfo=f.proj)
   
5. lstates = shaperead('D:/Temp/Map/states.shp')
   
6. geoshow(lstates, facecolor=[220,220,220], edgecolor='k')
   
7. layer = imshowm(data, 20, proj=f.proj, order=1)
   
8. colorbar(layer)
   
9. title('Base reflectivity (' + t.strftime('%Y-%m-%d %H:%M') + ')')

f = addfile('D:/Temp/grib/1kmradar_201512241600.grib2')<br /> data = f['Base_reflectivity_surface_layer'][0,:,:]<br /> t = f.gettime(0)<br /> axesm(projinfo=f.proj)<br /> lstates = shaperead('D:/Temp/Map/states.shp')<br /> geoshow(lstates, facecolor=[220,220,220], edgecolor='k')<br /> layer = imshowm(data, 20, proj=f.proj, order=1)<br /> colorbar(layer)<br /> title('Base reflectivity (' + t.strftime('%Y-%m-%d %H:%M') + ')')

风子

看着就跟Python代码一样一样的，包括最后的图

兰溪之水

王老师，WSR98D雷达才是刚需啊

MeteoInfo

兰溪之水 发表于 2016-1-16 19:45
王老师，WSR98D雷达才是刚需啊

我对雷达数据没什么经验，只是在GrADS的email list里看到有人问帖子里的数据为什么用GrADS打开会很慢，顺便试了一下。对于这种高分辨率且有投影（比如Lambert）的数据，GrADS在读取数据数组后要将数据投影到等经纬度投影，会比较慢。

MeteoInfo

风子 发表于 2016-1-16 19:07
看着就跟Python代码一样一样的，包括最后的图

就是Python呀

schliezer

老师你好，想问下这个addfile是那个模块里的？应该import 哪个模块呢？

MeteoInfo

schliezer 发表于 2016-7-7 17:01
老师你好，想问下这个addfile是那个模块里的？应该import 哪个模块呢？

addfile在minum模块中，在启动MeteoInfoLab的时候已经import过了。

schliezer

MeteoInfo 发表于 2016-7-7 21:59
addfile在minum模块中，在启动MeteoInfoLab的时候已经import过了。

好吧，看来还得自己在ipython里装遍pygrib 或者minum

schliezer

schliezer 发表于 2016-7-7 22:23
好吧，看来还得自己在ipython里装遍pygrib 或者minum

我试了下，好像并不存在minum这个包，无论pip install 还是conda install 都显示找不到此包

MeteoInfo

schliezer 发表于 2016-7-7 23:12
我试了下，好像并不存在minum这个包，无论pip install 还是conda install 都显示找不到此包

这个是MeteoInfoLab，是用Java和Jython开发的，不是传统的CPython。虽然语法都一样，但模块是互不相通的。minum是用Jython开发的模块，不能在CPython中运行。