Documentation

This is machine translation

Translated by Microsoft
Mouseover text to see original. Click the button below to return to the English verison of the page.

Note: This page has been translated by MathWorks. Please click here
To view all translated materals including this page, select Japan from the country navigator on the bottom of this page.

lwt2

2-D lifting wavelet transform

Syntax

[CA,CH,CV,CD] = lwt2(X,W)
X_InPlace = lwt2(X,LS)
lwt2(X,W,LEVEL)
X_InPlace = lwt2(X,W,LEVEL,'typeDEC',typeDEC)
[CA,CH,CV,CD] = LWT2(X,W,LEVEL,'typeDEC',typeDEC)

Description

lwt2 performs a 2-D lifting wavelet decomposition with respect to a particular lifted wavelet that you specify.

[CA,CH,CV,CD] = lwt2(X,W) computes the approximation coefficients matrix CA and detail coefficients matrices CH, CV, and CD, obtained by a lifting wavelet decomposition, of the matrix X. W is a lifted wavelet name (see liftwave).

X_InPlace = lwt2(X,LS) computes the approximation and detail coefficients. These coefficients are stored in place:

  • CA = X_InPlace(1:2:end,1:2:end)

  • CH = X_InPlace(2:2:end,1:2:end)

  • CV = X_InPlace(1:2:end,2:2:end)

  • CD = X_InPlace(2:2:end,2:2:end)

lwt2(X,W,LEVEL) computes the lifting wavelet decomposition at level LEVEL.

X_InPlace = lwt2(X,W,LEVEL,'typeDEC',typeDEC) or [CA,CH,CV,CD] = LWT2(X,W,LEVEL,'typeDEC',typeDEC) with typeDEC = 'w' or 'wp' computes the wavelet or the wavelet packet decomposition using lifting, at level LEVEL.

Instead of a lifted wavelet name, you may use the associated lifting scheme LS: lwt2(X,LS,...) instead of LWT2(X,W,...).

For more information about lifting schemes, see lsinfo.

Examples

% Start from the Haar wavelet and get the
% corresponding lifting scheme.
lshaar = liftwave('haar');

% Add a primal ELS to the lifting scheme.
els = {'p',[-0.125 0.125],0};
lsnew = addlift(lshaar,els);

% Perform LWT at level 1 of a simple image.
x = reshape(1:16,4,4);
[cA,cH,cV,cD] = lwt2(x,lsnew)

cA =

    5.7500   22.7500
   10.0000   27.0000


cH =

    1.0000    1.0000
    1.0000    1.0000


cV =

    4.0000    4.0000
    4.0000    4.0000


cD =

     0     0
     0     0

% Perform integer LWT of the same image.
lshaarInt = liftwave('haar','int2int');
lsnewInt = addlift(lshaarInt,els);
[cAint,cHint,cVint,cDint] = lwt2(x,lsnewInt)

cAint =

     3    11
     5    13


cHint =

     1     1
     1     1


cVint =

     4     4
     4     4


cDint =

     0     0
     0     0

More About

collapse all

Tips

When X represents an indexed image, X, as well as the output arrays cA,cH,cV,cD, or X_InPlace are m-by-n matrices. When X represents a truecolor image, it is an m-by-n-by-3 array, where each m-by-n matrix represents a red, green, or blue color plane concatenated along the third dimension.

For more information on image formats, see the image and imfinfo reference pages .

Algorithms

This function implements the polyphase algorithm.

lwt reduces to dwt with zero-padding extension mode and without extra-coefficients.

References

Strang, G.; T. Nguyen (1996), Wavelets and filter banks, Wellesley-Cambridge Press.

Sweldens, W. (1998), "The Lifting Scheme: a Construction of Second Generation of Wavelets," SIAM J. Math. Anal., 29 (2), pp. 511–546.

See Also

Introduced before R2006a

Was this topic helpful?