Filip Wasilewski wrote: > Actually you have not. The algorithm you presented gives completely > wrong results. Have a look at quick&dirty(TM) implementation bellow.
God grief. I followed the implementation in Ingrid Daubechies' and Wim Sweldens' original wavelet lifting paper (J. Fourier Anal. Appl., 4: 247-269, 1998). If you look at the factorized polyphase matrix for D4 (which gives the inverse transform), their implementation of the lifting steps the original paper is incorrect. That is scary. A free PDF or PS of the paper is available here: http://cm.bell-labs.com/who/wim/papers/factor/index.html Thanks a lot! Sturla Molden -- http://mail.python.org/mailman/listinfo/python-list