在 Python 中插值网格时忽略 NaN

Question

我有一个网格速度场，我想在 Python 中进行插值。目前我正在使用 scipy.interpolate 的RectBivariateSpline 来执行此操作，但我希望能够通过将网格中的某些值设置为 NaN 来定义我的字段的边缘。但是，当我这样做时，它会打乱整个网格的插值，从而有效地将其变为 NaN。显然这是 scipy fitpack 中的一个错误，那么解决这个问题的最佳方法是什么？我希望能够将 NaN 保留在网格中，以便稍后处理边缘和越界，但我不希望它影响网格其余部分的插值。

Answer 1

所有正确实现浮点的语言（包括 python）都允许您通过将数字与其自身进行比较来测试 NaN。

x 不等于 x 当且仅当 x 为 NaN。

您将能够使用它来相应地过滤您的数据集。

Answer 2

样条拟合/插值是全局的，因此即使是单个 nan 也可能会弄乱整个网格。