在多维numpy数组[重复项]的每个索引处切片不同的范围

问题描述 投票:0回答:1

我有一个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,但是性能对于这段代码很重要,因此我希望将其向量化。

我已经编辑了问题,以便更清楚地看到所需的输出和错误。

python arrays numpy slice
1个回答
1
投票
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]])
© www.soinside.com 2019 - 2024. All rights reserved.