想到最近几年最红的程式大概就是Google earth了,而在大气界中尺度里最有名的大概就是WRF了,要如何把两者mapping在一起,其实国外的学者已经把他们统合了,目前已经可以一起Display出来了,但是目前只有二维场可以被显示出来(有被释放出来)。如果想看3D特效的同好们可能就要去下载Vis5D或Vapor,如果有点钱的朋友也可以买Avizo这套软体(功能较强大),上述这三套软体可以呈现出模式的3D图形,小弟虽然有用过Vis5D&Vapor但是仅仅试用而已,而Avizo也有去上过一堂课(不过早就忘光光了),其实功能都差不多,有兴趣的人可以往这方面多多学习,因为这毕竟是未来的趋势。

     我想应该已经有人可以利用3D的wrfout放到google earth上面了,只是如果我会的话我应该就是Google的员工了(我也希望我是),所以目前仍然只能转成2D的wrfout,如果真的要看3D,个人比较推荐应该就是要用vapor这套软体了,因为他免钱又可以在window上面使用,废话不多说,赶快来介绍如何把wrfout_2D直接转程wrf.kml吧。

1. 资料来源来自于这位丹佛大学的学者[https://github.com/jbeezley/wrf2kmz],真的要感谢有人这么不藏私(也难怪国外学术发展比较先进,这可是一个月左右才发表的新版本)。

2. 下载里面的资料(wrf2kmz&readme),因为我是用linux安装的,所以想用window的人不要问我

3.接下来就是重点了,需要安装python这套软体,这个跟perl、qt、cvs都是程式语言,千万不要问我说为啥他不用简单的bash安装就好,就跟你说中文还英文一样,没有道理的,你爸说国语你就说国语,天生如此。而python还需要下面三个套装的lib

    simplekml      : http://code.google.com/p/simplekml/
    matplotlib     : http://matplotlib.sourceforge.net/
    netcdf4-python : http://code.google.com/p/netcdf4-python/

ps : python的安装lib方式 python setup.py install,这三个安装过程就要完全看经验了,我也是搞了一天才完成(linux每台状况都不一样阿)

4. 第三步骤完成之后其实就简单啦,直接键入python wrf2kmz.py wrfout U10 V10 ....就会出现这样的讯息

could not compile reprojection module  (出现这个讯息主要是因为mapping问题,因为每个domain设定(参考点)不同,而作者也有一套reprojection的程式,但是需要下载f2py这个转档的程式,另外还需要一些.py档案(类似include),我有安装这些东西,但是可能是因为版本的问题,所以还是有一些错误讯息,但是mapping的误差有变好)

Creating grround overlay from U10

Creating grround overlay from V10

5.会出现wrf.kml,直接丢到google earth就可以看图啦 !! 

这边给大家看个范例

范例

我输出了地表两公尺温度还有潜热跟土地利用还有地表高度。

如果安装有问题可以问我阿 ~~ 可是只限于我认识的人 !!!

相关文章