>>> 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)
我需要在e
中找到值为大于15的真行或行。我可以使用for循环迭代,但是,我想知道numpy是否有办法可以做这更有效率?
答案 0 :(得分:50)
获取至少有一个项目大于15的行号:
>>> np.where(np.any(e>15, axis=1))
(array([1, 2], dtype=int64),)
答案 1 :(得分:23)
您可以使用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]
...
答案 2 :(得分:1)
一种简单明了的方法::使用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)]
或a[(a>4) != 0]
而不是a[np.nonzero(a>4)]
,因为它们可以正确处理0维数组。有关更多详细信息,请参见documentation。如以下示例所示,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]])
>>> [print(i,j) for i,j in a]
0 0
0 4
3 4
4 0
4 3
4 4