将不规则数据经纬度添加到 Xarray 数据集进行插值

Question

所以我有一个 Xarray 数据集ds ，其中包含 3 个变量：（示例数据）

Data variables:
    CO2      (time, data) float32 400.0 400.0 ... 400.0
    ...

其中data维度是一维数组，我有匹配长度的一维数组，其中包含data每个条目的纬度和经度。 纬度和经度不定期采样。 理想情况下，我会添加纬度和经度，因此它显示如下：

Data variables:
    CO2      (time, data, lat, lon) float32 400.0 400.0 ... 400.0
    ...

如何将纬度/经度添加为维度，以便我可以使用 Xarray interp 函数对其进行插值？ （或者还有其他更好的方法来插入这些数据吗？）

例如： ds_interp = ds.interp(time=dt, lat=dlat, lon=dlon, assume_sorted=False)

Answer 1

在本节底部的 xarray 文档中，您可以看到一条注释，说明插值仅适用于常规网格。 我将不得不对数据进行网格化并对其进行插值，我希望避免这种情况，或者使用支持不规则网格的 scipy.interpolate 手动将其插入为 numpy 数组。