Python Numpy 2-dimensional array iteration

Question

Would like to build a list of indices into a 2 dimensional bool_ array, where True.

import numpy
arr = numpy.zeros((6,6), numpy.bool_)
arr[2,3] = True
arr[5,1] = True
results1 = [[(x,y) for (y,cell) in enumerate(arr[x].flat) if cell] for x in xrange(6)]
results2 = [(x,y) for (y,cell) in enumerate(arr[x].flat) if cell for x in xrange(6)]

results 1:

[[], [], [(2, 3)], [], [], [(5, 1)]]

results 2 is completely wrong

Goal:

[(2, 3),(5, 1)]

Any way to do this without flattening the list afterwards, or any better way to do this in general?

Answer 1

I think the function you're looking for is numpy.where. Here's an example:

>>> import numpy
>>> arr = numpy.zeros((6,6), numpy.bool_)
>>> arr[2,3] = True
>>> arr[5,1] = True
>>> numpy.where(arr)
(array([2, 5]), array([3, 1]))

You can turn this back into an index like this:

>>> numpy.array(numpy.where(arr)).T
array([[2, 3],
       [5, 1]])

Answer 2

>>> import numpy as np
>>> arr = np.zeros((6,6), np.bool_)
>>> arr[2,3] = True
>>> arr[5,1] = True
>>> np.argwhere(arr)
array([[2, 3],
       [5, 1]])