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.]])

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

1 个答案: