main.asv

用提升方法实现小波任意级数分解

clear all
close all

[X]=imread('lena.bmp');
map=gray(256);
m=length(X);
figure(1)
subplot(1,2,1);
imshow(X);title('原始图像');axis square;

tic
%整数小波变换
de_x=liftwavedec2(X,m,8);
de_x=uint8(de_x);
data =de_x;
% compress data
fprintf('Compresing data ... ')
[zipped,info] = norm2huff(data);
fprintf('Done!\n')

% decompress data
fprintf('Decompressing data ... ')
unzipped = huff2norm(zipped,info);
fprintf('Done!\n')

whos data zipped unzipped

[a,b]=size(data);
unz=reshape(unzipped,a,b);

% test it
isOK = isequal(de_x(:),unz(:))

%整数小波反变换
re_x=liftwaverec2(,m,8);
re_x=uint8(re_x);
% test it
isOK = isequal(X(:),re_x(:))

cmp=gray(256);
subplot(1,2,2);colormap(cmp);imshow(uint8(re_x));title('重建图');axis square;