获取索引映射的数组的子数组的有效方法

Question

我有一个矩阵说a。我需要得到它的一个子矩阵，它的索引基本上来自主矩阵索引的映射（这个映射不一定是 1-1）。我有以下代码来生成子矩阵，这里的映射被认为是sum。

import numpy as np
def transform(A):
    B=np.zeros(A.flatten().shape[0])
    for i in range(A.flatten().shape[0]):
        multi_idx=np.unravel_index(i,A.shape)
        B[np.sum(multi_idx)]=A[multi_idx] #the mapping applied on the indices: B[np.sum(multi_idx)]
    return B       
A=np.arange(27).reshape([3,3,3])
print A
print transform(A)        

有输出：

[[[ 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]]]
[  0.   9.  18.  21.  24.  25.  26.   0.   0.   0.   0.   0.   0.   0.   0.
   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.]

Answer 1

np.ogrid 可以方便地根据数组中的索引进行表达式。例如，

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)

产量

[  0.   9.  18.  21.  24.  25.  26.   0.   0.   0.   0.   0.   0.   0.   0.
   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.   0.]

注意赋值

B[X] = A

等价于

B[X.ravel()] = A.ravel()

作业是按从左到右的顺序完成的。因此，如果X 有许多重复值，则只有最后一个值最终会影响B。这具有以您希望的方式处理地图的非一对一性的效果。