A faster way than "for loop" - Indexing 3d numpy array by 2d array

Question

I have a 3d (3, 2, 3) array and first dimension(3) is flexible it can be any size.

arr = np.array(
[[[56, 24, 32],
  [56, 24, 32]],
 [[51, 27, 72],
  [51, 27, 72]],
 [[36, 14, 49],
  [36, 14, 49]]])

Indexing array is (2,3):

idxs = np.array(
[[1, 0, 2],
 [2, 1, 0]])

I wanna index arr by idxs. Expected result is :

[[[24 56 32]
  [32 24 56]]
 [[27 51 72]
  [72 27 51]]
 [[14 36 49]
  [49 14 36]]])

If i use for loop like below it takes a lot of time :

for i, arr2d in enumerate(arr):
    for j, (arr1d, idx) in enumerate(zip(arr2d, idxs)):
        arr[i, j] = arr1d[idx]

So my question is : How can i speed up this process ?

arr has a redundant axis: np.squeeze(arr[:,::2])[:,idxs]. — Michael Szczesny
– Michael Szczesny, Commented Apr 22, 2022 at 10:59

Daniel F · Accepted Answer · 2022-04-22 14:46:36Z

1

Use np.take_along_axis

np.take_along_axis(arr, idxs[None, ...], 2)
Out[]: 
array([[[24, 56, 32],
        [32, 24, 56]],

       [[27, 51, 72],
        [72, 27, 51]],

       [[14, 36, 49],
        [49, 14, 36]]])

answered Apr 22, 2022 at 14:46

Daniel F

14.5k2 gold badges34 silver badges59 bronze badges

Sign up to request clarification or add additional context in comments.

Collectives™ on Stack Overflow

A faster way than "for loop" - Indexing 3d numpy array by 2d array

1 Answer 1

Comments

Your Answer

Hot Network Questions

Collectives™ on Stack Overflow

1 Answer 1

Comments

Your Answer

Sign up or log in

Post as a guest

Related