Metpy 横截面测地线

Question

我可能在这里问错了问题，而且我远不是问题专家，但是在创建横截面时，是否可以遵循恒向线？该文档指出，它通过常规网格上的给定数据沿测地线获取垂直横截面切片。这似乎表明它将遵循大循环。

或者，这只是在调用横截面方法之前将我的数据放在某个投影中的问题吗？如果是这样，Metpy 中是否有任何设施可以做到这一点？我在其他 SO 搜索中发现了一些想法。

这些是我的一个横截面的纬度坐标。 我真正想要的是背阔肌保持不变。

From: (40.05, -80) lat, lon
To: (40.05, 20.06)

[40.05       40.49129217 40.92753221 41.35857652 41.78427825 42.20448739
 42.6190507  43.02781183 43.43061131 43.8272866  44.21767222 44.60159979
 44.97889817 45.34939353 45.71290958 46.06926767 46.41828699 46.75978478
 47.09357658 47.41947646 47.73729732 48.04685118 48.34794952 48.64040361
 48.92402495 49.19862562 49.46401872 49.72001885 49.96644255 50.20310884
 50.42983969 50.64646057 50.85280099 51.04869508 51.23398211 51.40850708
 51.5721213  51.72468292 51.86605751 51.99611856 52.11474806 52.22183695
 52.31728562 52.40100434 52.47291364 52.53294472 52.58103974 52.61715208
 52.64124659 52.65329974 52.65329974 52.64124659 52.61715208 52.58103974
 52.53294472 52.47291364 52.40100434 52.31728562 52.22183695 52.11474806
 51.99611856 51.86605751 51.72468292 51.5721213  51.40850708 51.23398211
 51.04869508 50.85280099 50.64646057 50.42983969 50.20310884 49.96644255
 49.72001885 49.46401872 49.19862562 48.92402495 48.64040361 48.34794952
 48.04685118 47.73729732 47.41947646 47.09357658 46.75978478 46.41828699
 46.06926767 45.71290958 45.34939353 44.97889817 44.60159979 44.21767222
 43.8272866  43.43061131 43.02781183 42.6190507  42.20448739 41.78427825
 41.35857652 40.92753221 40.49129217 40.05      ]

Answer 1

MetPy's cross_section 目前仅计算起点和终点纬度/经度点之间的测地线/大圆路径以用于垂直横截面。要使用替代路径（例如恒向线），可以将底层 interpolate_to_slice 函数与数据集水平维度坐标中的点数组一起使用。例如，假设您的数据是一维纬度/经度坐标（即等角投影），您的示例变为

start_lat, start_lon = (40.05, -80)
end_lat, end_lon = (40.05, 20.06)
n_points = 100

lats = np.linspace(start_lat, end_lat, n_points)
lons = np.linspace(start_lon, end_lon, n_points)

cross = interpolate_to_slice(data, np.stack([lons, lats], axis=1))

如果您的数据位于其他投影上，则需要获取正确转换的坐标并改用这些坐标。