获取X阵列中点到点的对角线上的所有栅格值

2024-03-28 20:20:09 发布

您现在位置:Python中文网/ 问答频道 /正文

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


Tags: 数据网格dsslice端点子集全球geo
1条回答
网友
1楼 · 发布于 2024-03-28 20:20:09

正如上面的评论所指出的,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'))

Here是绘制横截面的一个很好的示例

相关问题 更多 >