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

本文介绍了如何使用Python和GDAL库处理WRF输出的兰伯特双标准纬线投影数据。通过设置椭球体参数、投影中心参数,实现了从WRF数据的GrADS二进制码到地理坐标的投影转换。通过逐步调整和验证,最终降低误差,实现准确的坐标转换。
摘要由CSDN通过智能技术生成

d14db6971e8b39d78d747cb356f377e1.png

有一种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的小伙

  • 0
    点赞
  • 8
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值