如何利用python将wrfout的数据从兰勃特投影处理成等经纬网格数据

august2015

使用python，可以读取wrfout的nc格式数据，为兰勃特投影。为了和格点降水实况数据进行统计分析，如何将其处理成等经纬网格的数据？请大家指教。谢谢！

eeeee

可以用cartopy
https://scitools.org.uk/cartopy/ ... RS.transform_points

august2015

eeeee 发表于 2019-9-4 12:00
可以用cartopy
https://scitools.org.uk/cartopy/docs/v0.14/crs/index.html#cartopy.crs.CRS.transform_p ...

能具体指点一下吗？我是python新手，在网上查到的cartopy多用于绘图，处理数据插值如何使用呢？谢谢。

eeeee

1. import cartopy.crs as ccrs
   
2. import numpy as np
   
3. #构造原始投影，投影参数自己改
   
4. src = ccrs.LambertConformal(central_longitude=0, central_latitude=0)
   
5. #构造目标投影
   
6. target = ccrs.PlateCarree()
   
7. #测试数据
   
8. x = np.arange(10)
   
9. y = np.arange(10)
   
10. output = target.transform_points(src, x, y)
    
11. lon = output.T[0]
    
12. lat = output.T[1]

复制代码

august2015

eeeee 发表于 2019-9-4 15:16

非常感谢。如果是降水值，也能对应的插值到等经纬网格上吗？

eeeee

august2015 发表于 2019-9-4 16:03
非常感谢。如果是降水值，也能对应的插值到等经纬网格上吗？

这个不是插值，是坐标系转换，你要插值的话需要先转换成经纬度再插值

august2015

eeeee 发表于 2019-9-4 15:16

非常感谢。如果是降水值，也能对应的插值到等经纬网格上吗？

言墨

楼主 求问转换成功了吗

言墨

eeeee 发表于 2019-9-4 15:16

求问 那个原始投影的参数  根据什么改？

kisscate

感谢