Wavelet packet decomposition

From Infogalactic: the planetary knowledge core
Jump to: navigation, search

Lua error in package.lua at line 80: module 'strict' not found.

Originally known as Optimal Subband Tree Structuring (SB-TS) also called Wavelet Packet Decomposition (WPD) (sometimes known as just Wavelet Packets or Subband Tree) is a wavelet transform where the discrete-time (sampled) signal is passed through more filters than the discrete wavelet transform (DWT).

Introduction

In the DWT, each level is calculated by passing only the previous wavelet approximation coefficients (cAj) through discrete-time low and high pass quadrature mirror filters.[1] However, in the WPD, both the detail (cDj (in the 1-D case), cHj, cVj, cDj (in the 2-D case)) and approximation coefficients are decomposed to create the full binary tree.[2] [3][4][5][6][7]

File:Wavelets - WPD.png
Wavelet Packet decomposition over 3 levels. g[n] is the low-pass approximation coefficients, h[n] is the high-pass detail coefficients

For n levels of decomposition the WPD produces 2n different sets of coefficients (or nodes) as opposed to (3n + 1) sets for the DWT. However, due to the downsampling process the overall number of coefficients is still the same and there is no redundancy.

From the point of view of compression, the standard wavelet transform may not produce the best result, since it is limited to wavelet bases that increase by a power of two towards the low frequencies. It could be that another combination of bases produce a more desirable representation for a particular signal. The best basis algorithm by Coifman and Wickerhauser[1] finds a set of bases that provide the most desirable representation of the data relative to a particular cost function (e.g. entropy).

There were relevant studies in signal processing and communications fields to address the selection of subband trees (orthogonal basis) of various kinds, e.g. regular, dyadic, irregular, with respect to performance metrics of interest including energy compaction (entropy), subband correlations and others. [3] [4] [5] [6] [7]

Discrete wavelet transform theory (continuous in the variable(s)) offers an approximation to transform discrete (sampled) signals. In contrast, the discrete subband transform theory provides a perfect representation of discrete signals.[5]

Gallery

Applications

Wavelet packet was successfully applied in preclinical diagnosis.[8]

References

  1. 1.0 1.1 Coifman RR & Wickerhauser MV, 1992. Entropy-Based Algorithms for Best Basis Selection, IEEE Transactions on Information Theory, 38(2).
  2. Daubechies, I. (1992), Ten lectures on wavelets, SIAM
  3. 3.0 3.1 A.N. Akansu and Y. Liu, On Signal Decomposition Techniques, (Invited Paper), Optical Engineering Journal, special issue Visual Communications and Image Processing, vol.30, pp. 912-920, July 1991.
  4. 4.0 4.1 H. Caglar, Y. Liu and A.N. Akansu, Statistically Optimized PR-QMF Design, Proc. SPIE Visual Communications and Image Processing, vol. 1605, pp. 86-94, 1991.
  5. 5.0 5.1 5.2 A.N. Akansu and R.A. Haddad, Multiresolution Signal Decomposition: Transforms, Subbands, and Wavelets. Boston, MA: Academic Press, ISBN 978-0-12-047141-6, 1992.
  6. 6.0 6.1 A. Benyassine and A.N. Akansu, Performance Analysis and Optimal Structuring of Subchannels for Discrete Multitone Transceivers , Proc. IEEE Proc. IEEE International Symposium on Circuits and Systems (ISCAS), pp. 1456-1459, April 1995.
  7. 7.0 7.1 M.V. Tazebay and A.N. Akansu, Adaptive Subband Transforms in Time-frequency Excisers for DSSS Communications Systems, IEEE Trans. Signal Processing, vol. 43, pp. 2776-2782, Nov. 1995.
  8. Lua error in package.lua at line 80: module 'strict' not found.

External links

  • An implementation of wavelet packet decomposition can be found in MATLAB wavelet toolbox: [1].
  • An implementation for R can be found in the wavethresh package: [2].
  • An illustration and implementation of wavelet packets along with its code in C++ can be found at [3].
  • JWave: An implementation in Java for 1-D and 2-D wavelet packets using Haar, Daubechies, Coiflet, and Legendre wavelets.