从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]])
