前言
- 一般在处理海洋、气象数据时，往往会通过一些收支方程来对某一气象要素进行诊断分析，而大部分的诊断方程往往都会包含对于变量的偏导。比如dTdx，DTdy等。而一般得到的经纬度的单位为弧度，而dx的单位为m。所以需要将单位进行转换，以下简单介绍一下计算的方法
1、使用metpy的函数
metpy.calc.lat_lon_grid_deltas(longitude, latitude, x_dim=- 1, y_dim=- 2, geod=None)
只需要传入经度和纬度数据即可得到dx和dy
2、第二种方法，通过numpy计算
xlon,ylat=np.meshgrid(lon,lat) dlony,dlonx=np.gradient(xlon) dlaty,dlatx=np.gradient(ylat) pi=3.14159265 re=6.37e6 dx=re*np.cos(ylat*pi/180)*dlonx*pi/180 dy=re*dlaty*pi/180
以上两种方法亲测可行，之间直接跑通
