我有一个网格化的速度场,我想在Python中插值。目前我使用scipy.interpolate的RectBivariateSpline来实现这一点,但我希望能够通过将网格中的某些值设置为NaN来定义字段的边缘。然而,当我这样做时,它会搞乱整个网格的插值,有效地使其到处都是NaN。显然这是scipy fitpack中的一个错误,那么解决这个问题的最佳方法是什么呢?我希望能够保持网格中的NaN,以便稍后处理边缘和边界外,但我不希望它影响网格其余部分的插值。
RectBivariateSpline
lsmd5eda1#
所有正确实现浮点的语言(包括python)都允许您通过将数字与其本身进行比较来测试NaN。x不等于x,当且仅当x是NaN。您可以使用它来相应地过滤数据集。
x
ippsafx72#
样条线拟合/插值是全局的,因此即使是一个nan也可能会扰乱整个网格。
nan
2条答案
按热度按时间lsmd5eda1#
所有正确实现浮点的语言(包括python)都允许您通过将数字与其本身进行比较来测试NaN。
x
不等于x
,当且仅当x
是NaN。您可以使用它来相应地过滤数据集。
ippsafx72#
样条线拟合/插值是全局的,因此即使是一个
nan
也可能会扰乱整个网格。