我有一个像这样的numpy数组:
>>> I
array([[ 1., 0., 2., 1., 0.],
[ 0., 2., 1., 0., 2.]])
这样的阵列A:
>>> A = np.ones((2,5,3))
我想获得以下矩阵:
>>> result
array([[[ False, False, True],
[ False, True, True],
[ False, False, False],
[ False, False, True],
[ False, True, True]],
[[ False, True, True],
[ False, False, False],
[ False, False, True],
[ False, True, True],
[ False, False, False]]], dtype=bool)
最好通过示例进行说明:
I[0,0] = 1 -> result[0,0,:2] = False
和result[1,1,2:] = True
I[1,0] = 0 -> result[1,1,0] = False
和result[1,1,1:] = True
这是我目前的实施(正确):
result = np.empty((A.shape[0], A.shape[1], A.shape[2]))
r = np.arange(A.shape[2])
for i in xrange(A.shape[0]):
result[i] = r > np.vstack(I[i])
print result.astype(np.bool)
有没有办法以更快的方式实现(避免for循环)?
谢谢!
答案 0 :(得分:1)
您只需要在I
上添加其他维度,以便正确地广播r
:
result = r > I.reshape(I.shape[0],I.shape[1],1)
e.g。
In [41]: r>I.reshape(2,5,1)
Out[41]:
array([[[False, False, True],
[False, True, True],
[False, False, False],
[False, False, True],
[False, True, True]],
[[False, True, True],
[False, False, False],
[False, False, True],
[False, True, True],
[False, False, False]]], dtype=bool)