如何追加掩码数组

Question

我有两个数组，我想把它们合并成一个新的数组，但我需要保留被遮罩的信息。我试过使用numpy.append()，但是它会丢失这些被遮罩的信息。

>>> maska
masked_array(data = [-- 1 3 2 1 -- -- 3 6],
     mask = [ True False False False False  True  True False False],
     fill_value = 0)

>>> b
masked_array(data = [-- 1 3 2],
     mask = [ True False False False], fill_value = 0)

>>> np.append(maska,b)
masked_array(data = [0 1 3 2 1 0 0 3 6 0 1 3 2],
     mask = False, fill_value = 999999)

Answer 1

试试使用带掩码的类追加功能：np.ma.append

# make first Array
In[1]: maska = np.ma.array([999, 1,3,2,1, 999, 999, 3, 6],
                    mask = [True, False, False, False, False, True, True, False, False])
In[2]: maska
Out[2]: masked_array(data=[--, 1, 3, 2, 1, --, --, 3, 6],
                       mask=[ True, False, False, False, False,  True,  True, False, False], fill_value=999999)

# make 2nd array
In[3]: b = np.ma.array([999, 1, 3, 2], mask=[True, False, False, False])
In[4]: b
Out[4]: masked_array(data=[--, 1, 3, 2],
                     mask=[ True, False, False, False], fill_value=999999)

# Now append the two arrays 
In[5]: np.ma.append(maska, b)
Out[5]: masked_array(data=[--, 1, 3, 2, 1, --, --, 3, 6, --, 1, 3, 2],
                     mask=[ True, False, False, False, False,  True,  True, False,
               False,  True, False, False, False], fill_value=999999)

Answer 2

这确实很奇怪，连 np.ma.hstack 都不能用。不过，你可以通过手动组合这些遮罩来实现你想要的效果：

In [1]: import numpy as np

In [2]: def masked_hstack(tup):
   ...:     return np.ma.masked_array(np.hstack(tup),
   ...:            mask=np.hstack([arr.mask for arr in tup]))
   ...:

In [3]: a, b = [0, 1, 3, 2, 1, 0, 0, 3, 6], [0, 1, 3, 2]

In [4]: maska, maskb = [np.ma.masked_equal(arr, 0) for arr in a, b]

In [5]: masked_hstack((maska, maskb))
Out[5]:
masked_array(data = [-- 1 3 2 1 -- -- 3 6 -- 1 3 2],
             mask = [ True False False False False  True  True False False  True False False
 False],
       fill_value = 999999)

如果需要的话，你还可以把 fill_value 设置为 0。