metpy中流动矢量的地球相对和栅格相对变换

Question

我已经试着浏览了源代码，但是MetPy似乎不支持这一点。例如，我有一组相对于地球坐标系（速度、方向）的风矢量，我想将其转换为（u、v）相对于Lambert共形投影上的栅格

Answer 1

MetPy不处理坐标系变换，而是让CartoPy处理。它确实支持将速度和方向转换为组件。因此，您希望执行以下操作：

import cartopy.crs as ccrs
import metpy.calc as mpcalc

proj = ccrs.LambertConformal()

u_earth, v_earth = mpcalc.wind_components(speed, direction)
u_grid, v_grid = proj.transform_vectors(ccrs.PlateCarree(), lon, lat, u_earth, v_earth)

使用PlateCarree作为组件的源投影告诉CartoPy坐标和组件是lon/lat和地球相对的