141
>>> 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 是否可以更有效地做到这一点?

4

4 回答 4

128

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

>>> np.where(np.any(e>15, axis=1))
(array([1, 2], dtype=int64),)
于 2013-04-18T23:38:20.983 回答
78

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

简单的方法

>>> (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]
 ...
于 2014-02-16T18:34:49.347 回答
50

一种简单而干净的方法:用于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)),但它为 0 维数组生成正确形状的结果。

注意:您不能使用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
>>>
于 2019-09-01T04:23:21.150 回答
0

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

如评论中所述,np.where()NumPy 文档不鼓励这样做。

于 2021-12-07T20:15:00.840 回答