Ітерація через масивний масив


135

Чи є менш дослідна альтернатива цьому:

for x in xrange(array.shape[0]):
    for y in xrange(array.shape[1]):
        do_stuff(x, y)

Я придумав це:

for x, y in itertools.product(map(xrange, array.shape)):
    do_stuff(x, y)

Що зберігає один відступ, але все ще досить некрасиво.

Я сподіваюся на те, що буде схоже на цей псевдокод:

for x, y in array.indices:
    do_stuff(x, y)

Чи існує щось подібне?


Я в python 2.7 і використовую ваше рішення з itertools; я читав у коментарях, що використання itertools буде швидше. однак, (можливо, тому, що я в 2.7), мені також довелося розпакувати карту в циклі for. for x, y in itertools.product(*map(xrange, array.shape)):
ALM

У довідці NumPy є сторінка під назвою "Ітерація над масивами": docs.scipy.org/doc/numpy/reference/arrays.nditer.html
Кейсі

Відповіді:


187

Я думаю, ти шукаєш нечисленний номер .

>>> a =numpy.array([[1,2],[3,4],[5,6]])
>>> for (x,y), value in numpy.ndenumerate(a):
...  print x,y
... 
0 0
0 1
1 0
1 1
2 0
2 1

Щодо виступу. Це трохи повільніше, ніж розуміння списку.

X = np.zeros((100, 100, 100))

%timeit list([((i,j,k), X[i,j,k]) for i in range(X.shape[0]) for j in range(X.shape[1]) for k in range(X.shape[2])])
1 loop, best of 3: 376 ms per loop

%timeit list(np.ndenumerate(X))
1 loop, best of 3: 570 ms per loop

Якщо ви турбуєтеся про ефективність, ви можете трохи оптимізувати, переглянувши реалізацію ndenumerate, яка робить 2 речі, перетворившись на масив та циклічно. Якщо ви знаєте, що у вас є масив, ви можете викликати .coordsатрибут плоского ітератора.

a = X.flat
%timeit list([(a.coords, x) for x in a.flat])
1 loop, best of 3: 305 ms per loop

1
Зауважте, це працює, але це надзвичайно повільно. Вам краще повторити вручну.
Марті

43

Якщо вам потрібні лише індекси, ви можете спробувати numpy.ndindex:

>>> a = numpy.arange(9).reshape(3, 3)
>>> [(x, y) for x, y in numpy.ndindex(a.shape)]
[(0, 0), (0, 1), (0, 2), (1, 0), (1, 1), (1, 2), (2, 0), (2, 1), (2, 2)]

15

див. nditer

import numpy as np
Y = np.array([3,4,5,6])
for y in np.nditer(Y, op_flags=['readwrite']):
    y += 3

Y == np.array([6, 7, 8, 9])

y = 3не буде працювати, використовувати y *= 0і y += 3замість цього.


2
або скористайтеся y [...] = 3
Дональд Хобсон
Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.