我有一个m x n numpy数组arr
,对于arr
的每一列,我都有要访问的给定行范围。我有一个n x 1数组vec
,它描述此范围何时开始,并且该范围具有一定的持续时间d
。如何有效地提取此d x n感兴趣的数组?可以通过巧妙的切片来完成吗?我最初的想法是尝试类似的方法:
arr = np.tile(np.arange(10),(4,1)).T
vec = np.array([3,4,5,4])
d = 3
vec_2 = vec+d
out = arr[vec:vec2,np.arange(n)]
但是这会产生以下错误:TypeError: only integer scalar arrays can be converted to a scalar index
。所需的输出将是以下数组:
array([[3, 4, 5, 4],
[4, 5, 6, 5],
[5, 6, 7, 6],
[6, 7, 8, 7])
我可以遍历d
,但是性能对于这段代码很重要,因此我希望将其向量化。
我已经编辑了问题,以便更清楚地看到所需的输出和错误。
In [489]: arr=np.arange(24).reshape(6,4)
In [490]: vec=np.array([0,2,1,3])
利用linspace
的最新扩展来生成多个数组:
In [493]: x = np.linspace(vec,vec+2,3).astype(int)
In [494]: x
Out[494]:
array([[0, 2, 1, 3],
[1, 3, 2, 4],
[2, 4, 3, 5]])
In [495]: arr[x, np.arange(4)]
Out[495]:
array([[ 0, 9, 6, 15],
[ 4, 13, 10, 19],
[ 8, 17, 14, 23]])
列迭代方法:
In [498]: np.stack([arr[i:j,k] for k,(i,j) in enumerate(zip(vec,vec+3))],1)
Out[498]:
array([[ 0, 9, 6, 15],
[ 4, 13, 10, 19],
[ 8, 17, 14, 23]])