在 numpy masked_array 中用 nan 替换掩码

Question

>> masks = [[1,1],[0,0]]    
>> [np.ma.masked_array(data=np.array([1.0,2.0]), mask=m, fill_value=np.nan).mean() for m in masks]
   [masked, 1.5]

我想用 nan 替换 masked 结果。有没有办法直接用 numpy 的 masked_array?

Answer 1

我觉得你可以np.ones

masks=np.array([[1,1],[0,0]])

np.ma.masked_array(data=np.array([1.0,2.0])*np.ones(masks.shape), mask=masks, fill_value=np.nan).mean(axis=1)
Out[145]: 
masked_array(data=[--, 1.5],
             mask=[ True, False],
       fill_value=1e+20)

Answer 2

In [232]: M = np.ma.masked_array(data=np.array([1.0,2.0]),mask=[True, False])

filled 方法用填充值替换屏蔽值：

In [233]: M.filled()                                                         
Out[233]: array([1.e+20, 2.e+00])
In [234]: M.filled(np.nan)         # or with a value of your choice.                                                   
Out[234]: array([nan,  2.])

或者像你一样，在定义数组的时候指定填充值：

In [235]: M = np.ma.masked_array(data=np.array([1.0,2.0]),mask=[True, False],
     ...:  fill_value=np.nan)                                                
In [236]: M                                                                  
Out[236]: 
masked_array(data=[--, 2.0],
             mask=[ True, False],
       fill_value=nan)
In [237]: M.filled()                                                         
Out[237]: array([nan,  2.])

屏蔽均值方法跳过填充值：

In [238]: M.mean()                                                           
Out[238]: 2.0
In [239]: M.filled().mean()                                                  
Out[239]: nan
In [241]: np.nanmean(M.filled())    # so does the `nanmean` function
In [242]: M.data.mean()             # mean of the underlying data                                                      
Out[242]: 1.5

在 numpy masked_array 中用 nan 替换掩码

replace masked with nan in numpy masked_array

python

numpy

masked-array