我可能在这里问错了问题，而且我远非一个问题专家，但是在创建横截面时，是否可能遵循主干线?文档中指出，它需要一个垂直的横截面切片沿测地线通过一个规则网格上的给定数据。这似乎表明它将遵循大圆。
或者，这是否都是在调用横截面方法之前将我的数据置于某个投影中的问题?如果有，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 ]
MetPy的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))
如果你的数据在其他投影上，你将需要获得正确转换的坐标并使用它们。