python - 不挂在对齐点上的 scipy.interpolate.griddata 的替代品

Question

我有一个点数据集，我试图在网格上插值。这些点以网格方式对齐，缺少一些点，如下所示：

更复杂的是，其他输入点集可能无法在网格上对齐，因此我尝试使用scipy.interpolate.griddata将这些值插入到常规网格中。但是，有时我的底层网格与输入点数据集的采样率完全对齐并griddata挂起。根据这个问题 scipy.interpolate.griddata，如果样本出现在完全对齐的 3 个点上，这对我来说有时会发生。

是否有将点数据插值到 Python 中的常规网格的高性能替代方案？

Answer 1

我的解决方案是随机扰动每个输入点的 x 和 y 坐标，随机因子介于1e-6 * grid_cell_size. 1e-6根据经验，我选择比 32 位浮点数中的最小 delta 差大几个数量级，但又足够小，以至于它引入的误差仍然被插值方案引入的误差所掩盖。

我仍然对其他想法持开放态度，但这样我可以scipy通过引入一个非常小的数值误差来使用股票插值函数。