最近几天,我的朋友小宋试用了国家气象局开源的预报检验库meteva。
啊~~都是业内大佬
排名第一的刘凑华,是我和小宋的学长,所以我们经常亲切地管meteva叫:“老刘的库”
为了试用meteva的一些具体功能,小宋还特意准备了一套数据,包括micaps第3类格式的站点数据,还有一些卫星和雷达资料。我问她为什么要用卫星和雷达资料,直接搞一套EC资料不就可以测试格点数据的功能了吗?她说,不,我要用老刘的库试试QPE。
WOW~果然是进击的小宋。让我们祝她好运
值得一提的是,meteva程序库中除了检验部分,还有其他一些我们气象行业常用的算法,比如计算平均值、最大值、最小值、RMSE,根据风速风向计算UV分量,根据温度露点计算相对湿度,统计回归,还有站点与格点之间的相互插值,等等。
这就意味着,哪怕你并不需要做检验,也可以安装meteva来计算和插值哦!而且meteva插值速度还不慢呢!一举解决了python在气象方面的插值死慢的问题
此外,meteva还自带绘制填色图、曲线图、柱状图等命令,一句话就能画图,感觉棒棒哒!
有一点需要提醒一下大家:meteva的绘图模块都是基于matplotlib库的,在windows下绘图没什么问题,但是如果用shell远程linux来绘图的话,需要安装配套的xmanager,否则会因为无法调用xmanager而绘图失败。
还有要注意的是经纬度问题。meteva是将数据视作dataframe来处理的,而且由于python本身对于数据的存储顺序与ncl相反,用惯了ncl的人,容易在使用meteva时出现将经纬度写反了的情况。不放心的话,就及时将数据print出来看看吧。比如这里,小宋就把经纬度写反了:
(为了防止站点信息泄露,小宋还贴心地打上了马赛克发给我。。)