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

Question

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

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

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

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 文档at the bottom of this section 中，您可以看到一条注释说插值仅适用于常规网格。我将不得不对数据进行网格化并对其进行插值，这是我希望避免的，或者使用支持不规则网格的 scipy.interpolate 手动将其插入为 numpy 数组。