python - 获取xarray中点到点对角线上的所有网格值

Question

我希望从xarray一个点和另一个点之间的数据集中选择所有数据。这当然很容易，sel例如使用基于纬度/经度网格的数据集从全球数据集到大西洋海岸线的子集我可以做到ds_geo_subset = ds.sel( lat=slice(42, 20), lon=slice(-85,-65))，但我想选择墨西哥湾流前沿的所有数据，基本上是从到点（- 78 lon, 32 lat) 到 (-72 lon, 36 lat) 但我不知道该怎么做。最好通过给出一个完整的点数据集并从每个点获取数据来完成，还是可以将这两个端点提供给 xarray 并返回插值数据？

Answer 1

正如上面的评论所建议的，Metpy有一个特定的功能可以方便地完成这项任务。

from metpy.interpolate import cross_section

data = xr.open_dataset('my_data')

# Define start and end points:

start = (-78, 32)
end = (-72, 36)
cross = cross_section(data, start, end).set_coords(('lat', 'lon'))

这是一个绘制横截面的好例子。