忽略 nan 值并执行 numpy.polyval 的函数

Question

我正在尝试编写一个函数，在其中提供一个 mxn 数组，它的目标是遍历数组并为每个值执行 numpy.polyval。数组将包含一些 nan 值，因此我希望函数以与原始数组相同的顺序返回 nan 值。 到目前为止，这是我的功能：

def polyval_NaN(p,data):
    m,n = data.shape
    data2 = data
    for i in np.arange(1,m):
        for j in np.arange(1,n):
            if np.isnan(data2[i,j]) == False:
                data2[i,j] = np.polyval(p,data[i,j])
    return data2

感谢您的帮助

Answer 1

像这样使用布尔索引：

>>> import numpy as np
>>> 
>>> p = np.array((4,5,3))
>>> a = np.arange(12.).reshape(3,4)
>>> a.ravel()[[2,5,7]] = np.nan
>>> a
array([[ 0.,  1., nan,  3.],
       [ 4., nan,  6., nan],
       [ 8.,  9., 10., 11.]])
>>> 
>>> result = np.full_like(a, np.nan)
>>> good = ~np.isnan(a)
>>> result[good] = np.polyval(p, a[good])
>>> 
>>> result
array([[  3.,  12.,  nan,  54.],
       [ 87.,  nan, 177.,  nan],
       [299., 372., 453., 542.]])