从numpy数组中取出具有给定步幅/步长的子arrays

可以说我有一个Python Numpy数组数组a。

numpy.array([1,2,3,4,5,6,7,8,9,10,11].) 

我想从这个长度为5的数组中创build一个子序列matrix,结果matrix因此将如下所示:

 numpy.array([[1,2,3,4,5],[4,5,6,7,8],[7,8,9,10,11]]) 

一种可能的实现方法是使用for循环。

 result_matrix = np.zeros((3, 5)) for i in range(0, len(a), 3): result_matrix[i] = a[i:i+5] 

有没有一个更清洁的方式来实现这个是Numpy?

方法1:使用broadcasting

 def broadcasting_app(a, L, S ): # Window len = L, Stride len/stepsize = S nrows = ((a.size-L)//S)+1 return a[S*np.arange(nrows)[:,None] + np.arange(L)] 

方法2:使用更高效的NumPy strides

 def strided_app(a, L, S ): # Window len = L, Stride len/stepsize = S nrows = ((a.size-L)//S)+1 n = a.strides[0] return np.lib.stride_tricks.as_strided(a, shape=(nrows,L), strides=(S*n,n)) 

样品运行 –

 In [143]: a Out[143]: array([ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11]) In [144]: broadcasting_app(a, L = 5, S = 3) Out[144]: array([[ 1, 2, 3, 4, 5], [ 4, 5, 6, 7, 8], [ 7, 8, 9, 10, 11]]) In [145]: strided_app(a, L = 5, S = 3) Out[145]: array([[ 1, 2, 3, 4, 5], [ 4, 5, 6, 7, 8], [ 7, 8, 9, 10, 11]])