python处理wrf气象数据_气象编程 | Python3之WRF的投影转换

python处理wrf气象数据_气象编程 | Python3之WRF的投影转换_第1张图片

有一种WRF输出的数据采用兰伯特双标准纬线投影,那么除非刚好需要同样的投影,想对这种数据进行处理的话往往要进行投影转换,WRF应该是有一套工具可以进行相关的处理,比如wrf-python,但是作为并不熟悉wrf、仅仅是使用WRF输出数据的小白,去使用WRF系工具的话学习成本就比较高了,如何用熟悉、更通用的工具实现这一投影转换呢?

难道不是设置几个投影参数,用常见的投影相关的包就可以实现了吗?

对,问题在于这个参数怎么设置?这个坑还是很坑的,好在最终找到一篇2018年的英文博客https://fabienmaussion.info/2018/01/06/wrf-projection/,加上我自己的尝试和理解,梳理出本文

采用WRF输出的数据格式为GrADS二进制码

关键

从我的角度来看,WRF的兰伯特双标准纬线投影有两个坑:椭球体参数

投影中心参数

完全不了解WRF输出的情况下

“啪”的一下,很快啊,观察到WRF输出的GrADS二进制数据对应的ctl文件中pdef如下所示1pdef  288 288 lcc  32.318  117.203  144.500  144.500  60.00000  30.00000  117.30000   3000.000   3000.000

然后查到pdef的lcc语法如下:

那么熟悉GDAL的小伙

你可能感兴趣的:(python处理wrf气象数据)