2

我正在尝试按照屏蔽广播的方式做一些事情，其中​​只广播某些值。

假设我有一个更大的数组更大的数组和一个更小的数组smaller_array：

import numpy as np
import numpy.ma as ma

bigger_array = np.zeros((4,4), dtype=np.int32)
smaller_array = np.ones((2,2), dtype=np.int32)

现在，我只想要较小数组的前三个值来替换较大数组的某个部分的值，但是掩码并没有达到我希望的效果：

masked_smaller_array = ma.masked_array(smaller_array, mask=[(0, 0), (0, 1)])
bigger_array[2:4, 2:4] = masked_smaller_array

这只是返回与常规广播相同的内容，即：

[[0 0 0 0]
 [0 0 0 0]
 [0 0 1 1]
 [0 0 1 1]]

而不是我所希望的

[[0 0 0 0]
 [0 0 0 0]
 [0 0 1 1]
 [0 0 1 0]]

在通过覆盖之前剥离屏蔽值

bigger_array[2:4, 2:4] = masked_smaller_array[~masked_smaller_array.mask]

也没有用，因为这会使阵列变平，使广播不兼容。

有没有其他方法可以达到同样的效果？

4

2 回答 2

1

当你有——

bigger_array[2:4, 2:4] = masked_smaller_array[~masked_smaller_array.mask]

您只需要索引等式的左侧即可使用相同的掩码。因此，解决它的一种方法是修复它，就像这样 -

# Mask corresponding to smaller array from where elements are to be taken
select_mask = ~masked_smaller_array.mask

# Use the mask on source (smaller array) to select specific elements from it 
# and update sliced and masked (with same mask) places in bigger array
bigger_array[2:4, 2:4][select_mask] = smaller_array[select_mask]

样品运行 -

In [59]: bigger_array = np.zeros((4,4), dtype=np.int32)
    ...: smaller_array = np.ones((2,2), dtype=np.int32)
    ...: masked_smaller_array =ma.masked_array(smaller_array,mask=[(0, 0),(0, 1)])
    ...: 

In [60]: select_mask = ~masked_smaller_array.mask

In [61]: select_mask
Out[61]: 
array([[ True,  True],
       [ True, False]], dtype=bool)

In [62]: bigger_array[2:4, 2:4][select_mask] = smaller_array[select_mask]

In [63]: bigger_array
Out[63]: 
array([[0, 0, 0, 0],
       [0, 0, 0, 0],
       [0, 0, 1, 1],
       [0, 0, 1, 0]], dtype=int32)
于 2015-12-19T21:35:24.447 回答
0

实际上可以将您的掩码用作条件，并在条件为真时numpy.where填充它。smaller_array顺便说一句，您可以使用 IIUCmake_mask代替masked_array：

your_mask = [(0, 0), (0, 1)]
mask = ma.make_mask(your_mask)
np.where(~mask, smaller_array, bigger_array[2:4, 2:4])


In [106]: mask
Out[106]: 
array([[False, False],
       [False,  True]], dtype=bool)

In [108]: np.where(~mask, smaller_array, bigger_array[2:4, 2:4])
Out[108]: 
array([[1, 1],
       [1, 0]], dtype=int32)

然后你可以为你的变量分配它

于 2015-12-19T21:55:27.873 回答