-
-
Notifications
You must be signed in to change notification settings - Fork 524
Open
Labels
Description
The output of waverec2 is not always consistent with the input of wavedec2. For example, the following
import numpy as np
import pywt
nx = 129
ny = 255
data = np.random.randn(nx, ny)
alpha = pywt.wavedec2(data, 'db1', mode='zero', level=1)
xrec = pywt.waverec2(alpha, 'db1', mode='zero')
print(data.shape, xrec.shape)
print(np.abs(data - xrec[0:data.shape[0],0:data.shape[1]]).max())produces
(129, 255) (130, 256)
1.7763568394002505e-15They do at least seem to be consistent within the overlap region. Is this the expected behaviour?