recursive_cc.py 文件源码

def Gaussian1D(signal, sigma, padding=0):
    n = signal.shape[0]
    tmp = np.zeros(n + padding)
    if tmp.shape[0] < 4:
        raise ValueError('Signal and padding too short')
    tmp[:n] = signal
    B, A = __gausscoeff(sigma)
    tmp = lfilter(B, A, tmp)
    tmp = tmp[::-1]
    tmp = lfilter(B, A, tmp)
    tmp = tmp[::-1]
    return tmp[:n]