python小波分解与重构_python - 使用pyWavelets进行多级局部小波重构 - 堆栈内存溢出...

我设法编写了自己的wrcoef函数版本，该版本似乎可以正常工作：

import pywt

import numpy as np

def wrcoef(X, coef_type, coeffs, wavename, level):

N = np.array(X).size

a, ds = coeffs[0], list(reversed(coeffs[1:]))

if coef_type =='a':

return pywt.upcoef('a', a, wavename, level=level)[:N]

elif coef_type == 'd':

return pywt.upcoef('d', ds[level-1], wavename, level=level)[:N]

else:

raise ValueError("Invalid coefficient type: {}".format(coef_type))

level = 4

X = [0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17]

coeffs = pywt.wavedec(X, 'db1', level=level)

A4 = wrcoef(X, 'a', coeffs, 'db1', level)

D4 = wrcoef(X, 'd', coeffs, 'db1', level)

D3 = wrcoef(X, 'd', coeffs, 'db1', 3)