是否有一个numpy函数有效地允许nwise迭代?

# http://seriously.dontusethiscode.com/2013/04/28/nwise.html
from itertools import tee, islice
nwise = lambda xs, n=2: zip(*(islice(xs, idx, None) for idx, xs in enumerate(tee(xs, n))))


例如。将均值应用于nwise元素?要获得移动平均线?

最佳答案

一般用途:

import numpy as np
from numpy.lib.stride_tricks import as_strided

def moving_slice(a, k):
    a = a.ravel()
    return as_strided(a, (a.size - k + 1, k), 2 * a.strides)

均线更好:
def moving_avg(a, k):
    ps = np.cumsum(a)
    return (ps[k-1:] - np.r_[0, ps[:-k]]) / k

例:
a = np.arange(10)

moving_avg(a, 4)
# array([ 1.5,  2.5,  3.5,  4.5,  5.5,  6.5,  7.5])

ms = moving_slice(a, 4)
ms
# array([[0, 1, 2, 3],
#        [1, 2, 3, 4],
#        [2, 3, 4, 5],
#        [3, 4, 5, 6],
#        [4, 5, 6, 7],
#        [5, 6, 7, 8],
#        [6, 7, 8, 9]])

# no data are copied:
a[4] = 0
ms
# array([[0, 1, 2, 3],
#        [1, 2, 3, 0],
#        [2, 3, 0, 5],
#        [3, 0, 5, 6],
#        [0, 5, 6, 7],
#        [5, 6, 7, 8],
#        [6, 7, 8, 9]])

关于python - python 。 nwise numpy数组迭代,我们在Stack Overflow上找到一个类似的问题:https://stackoverflow.com/questions/47506663/

10-12 13:05