xarray多维坐标插值提取时间序列

Question

我有一个多维坐标数据集，其中lon和lat取决于(x,y)维度：

<xarray.Dataset>
Dimensions:  (y: 310, x: 225, time: 52560, level: 10)
Coordinates:
    lat      (y, x) float32 ...
    lon      (y, x) float32 ...
  * time     (time) datetime64[ns] 2017-01-01 ... 2017-12-31T23:50:00
  * x        (x) int16 0 1 2 3 4 5 6 7 8 ... 216 217 218 219 220 221 222 223 224
  * y        (y) int16 0 1 2 3 4 5 6 7 8 ... 301 302 303 304 305 306 307 308 309
  * level    (level) float32 40.0 60.0 80.0 100.0 ... 170.0 200.0 250.0 300.0
Data variables:
    wspd     (time, level, y, x) float32 dask.array<chunksize=(100, 10, 310, 225), meta=np.ndarray>
Attributes: (12/130)
    TITLE:                           D3E5
    START_DATE:                      2011-07-01_00:00:00
    SIMULATION_START_DATE:           2010-01-01_00:00:00
    WEST-EAST_GRID_DIMENSION:        [256]
    SOUTH-NORTH_GRID_DIMENSION:      [352]
    BOTTOM-TOP_GRID_DIMENSION:       [50]
    ...                              ...
    ISOILWATER:                      [14]
    history:                         Thu Aug  8 14:07:48 2019: ncrcat -7 -L 3...
    NCO:                             4.4.2
    title:                           D3E5
    nco_openmp_thread_number:        [1]
    actual_range:                    [100. 100.]

在局部坐标系中， x的值范围为 0-309， y的范围为 0-224。

为了将我的wspd变量的时间序列的数据过滤到给定的lon, lat坐标，我按照这篇文章找到了最近的点。 这对我很有用。

现在感兴趣的问题：

最近的邻居对我来说是不够的，我想在高级插值下的xarray中进行插值，方法是使用xr.DataArray.interp(x, y, method="cubic")为我的观点获得更具代表性的时间序列兴趣。 只是我不能这样做，因为.interp()不适用于多维坐标，例如lon(x,y)和lat(x,y) 。 一个想法是找到与我给定的(lon,lat)坐标完全相同的(x_actual,y_actual)坐标，而不仅仅是最近的坐标，并使用.interp(x_actual,y_actual, method="cudic") 。 但我不知道这是否可行以及如何去做。

有人知道如何为多维坐标中给定的点获取插值数据变量吗？

感谢您的帮助，如果我能提供更多信息，请告诉我。