def uwt_inv(a, d, wavelet, boundary, level, initial_signal=None, axis=None):
hkr = nputils.atrou(get_wavelet_obj(wavelet).get_rec_hk(), pow(2, level))
gkr = nputils.atrou(get_wavelet_obj(wavelet).get_rec_gk(), pow(2, level))
c1 = nputils.convolve(a, hkr, boundary, axis=axis, mode="valid")
c2 = nputils.convolve(d, gkr, boundary, axis=axis, mode="valid")
return 1 / 2. * (c1 + c2)
def uwt(signal, wavelet, boundary, level, initial_signal=None, axis=None):
hkd = nputils.atrou(get_wavelet_obj(wavelet).get_dec_hk(), pow(2, level))
gkd = nputils.atrou(get_wavelet_obj(wavelet).get_dec_gk(), pow(2, level))
a = nputils.convolve(signal, hkd, boundary, axis=axis)
d = nputils.convolve(signal, gkd, boundary, axis=axis)
return (a, d)
def uwt_inv(a, d, wavelet, boundary, level, initial_signal=None, axis=None):
hkr = nputils.atrou(get_wavelet_obj(wavelet).get_rec_hk(), pow(2, level))
gkr = nputils.atrou(get_wavelet_obj(wavelet).get_rec_gk(), pow(2, level))
c1 = nputils.convolve(a, hkr, boundary, axis=axis, mode="valid")
c2 = nputils.convolve(d, gkr, boundary, axis=axis, mode="valid")
return 1 / 2. * (c1 + c2)
def uwt(signal, wavelet, boundary, level, initial_signal=None, axis=None):
hkd = nputils.atrou(get_wavelet_obj(wavelet).get_dec_hk(), pow(2, level))
gkd = nputils.atrou(get_wavelet_obj(wavelet).get_dec_gk(), pow(2, level))
a = nputils.convolve(signal, hkd, boundary, axis=axis)
d = nputils.convolve(signal, gkd, boundary, axis=axis)
return (a, d)
def uiwt(signal, wavelet, boundary, level, initial_signal=None, axis=None):
hkd = nputils.atrou(get_wavelet_obj(wavelet).get_dec_hk(), pow(2, level))
a = nputils.convolve(signal, hkd, boundary, axis=axis, mode='same')
d = signal - a
return (a, d)
def uiwt(signal, wavelet, boundary, level, initial_signal=None, axis=None):
hkd = nputils.atrou(get_wavelet_obj(wavelet).get_dec_hk(), pow(2, level))
a = nputils.convolve(signal, hkd, boundary, axis=axis, mode='same')
d = signal - a
return (a, d)