python--计算纬度/经度格式的网格点之间的实际距离

前言

• 一般在处理海洋、气象数据时，往往会通过一些收支方程来对某一气象要素进行诊断分析，而大部分的诊断方程往往都会包含对于变量的偏导。比如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

以上两种方法亲测可行，之间直接跑通