numpy获取索引值为true

Question

>>> ex=np.arange(30)
>>> e=np.reshape(ex,[3,10])
>>> e
array([[ 0,  1,  2,  3,  4,  5,  6,  7,  8,  9],
       [10, 11, 12, 13, 14, 15, 16, 17, 18, 19],
       [20, 21, 22, 23, 24, 25, 26, 27, 28, 29]])
>>> e>15
array([[False, False, False, False, False, False, False, False, False,
        False],
       [False, False, False, False, False, False,  True,  True,  True,
         True],
       [ True,  True,  True,  True,  True,  True,  True,  True,  True,
         True]], dtype=bool)

我需要找到具有true或行的行,e其值大于15.我可以使用for循环进行迭代,但是,我想知道numpy是否可以更有效地执行此操作？

Answer 1

要获取至少一个项目大于15的行号:

>>> np.where(np.any(e>15, axis=1))
(array([1, 2], dtype=int64),)

Answer 2

一种简单而干净的方法：使用np.argwhere按元素而不是维度对索引进行分组np.nonzero(a)（即，np.argwhere为每个非零元素返回一行）。

>>> a = np.arange(10)
>>> a
array([0, 1, 2, 3, 4, 5, 6, 7, 8, 9])
>>> np.argwhere(a>4)
array([[5],
       [6],
       [7],
       [8],
       [9]])

np.argwhere(a)与 相同np.transpose(np.nonzero(a))。

注意：您不能使用a(np.argwhere(a>4))来获取a. 推荐的方法是使用a[(a>4).astype(bool)]ora[(a>4) != 0]而不是a[np.nonzero(a>4)]因为它们正确处理 0-d 数组。有关更多详细信息，请参阅文档。如在下面的例子中可以看出，a[(a>4).astype(bool)]并且a[(a>4) != 0]可以被简化为a[a>4]。

另一个例子：

>>> a = np.array([5,-15,-8,-5,10])
>>> a
array([  5, -15,  -8,  -5,  10])
>>> a > 4
array([ True, False, False, False,  True])
>>> a[a > 4]
array([ 5, 10])
>>> a = np.add.outer(a,a)
>>> a
array([[ 10, -10,  -3,   0,  15],
       [-10, -30, -23, -20,  -5],
       [ -3, -23, -16, -13,   2],
       [  0, -20, -13, -10,   5],
       [ 15,  -5,   2,   5,  20]])
>>> a = np.argwhere(a>4)
>>> a
array([[0, 0],
       [0, 4],
       [3, 4],
       [4, 0],
       [4, 3],
       [4, 4]])
>>> for i,j in a: print(i,j)
... 
0 0
0 4
3 4
4 0
4 3
4 4
>>>

Answer 3

您可以使用非零功能.它返回给定输入的非零索引.

简单的方法

>>> (e > 15).nonzero()

(array([1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2]), array([6, 7, 8, 9, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9]))

看到指数更清洁,使用transpose方法:

>>> numpy.transpose((e>15).nonzero())

[[1 6]
 [1 7]
 [1 8]
 [1 9]
 [2 0]
 ...

不坏的方式

>>> numpy.nonzero(e > 15)

(array([1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2]), array([6, 7, 8, 9, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9]))

或干净的方式:

>>> numpy.transpose(numpy.nonzero(e > 15))

[[1 6]
 [1 7]
 [1 8]
 [1 9]
 [2 0]
 ...

Answer 4

当您只需要行 idx 时，我更喜欢np.flatnonzero(arr)该选项。有效，但它返回一个元组而不是数组。相当于.nonzero()arr.nonzero()flatnonzero()np.nonzero(np.ravel(arr))[0]

正如评论中提到的，np.where()NumPy 文档不鼓励这样做。