import numpy as np
A = np.arange(27).reshape([3, 3, 3])
B = np.zeros(A.size)
i, j, k = np.ogrid[0:3, 0:3, 0:3]
B[i + j + k] = A
print(B)
Эффективный способ получения подмассива массива, в котором отображаются индексы