program_13_10.m

《小波分析理论与Matlab R2007实现》源代码

% 装载原始图像
load belmont2; 
nbc = size(map,1);
% 使用sym4小波对图像进行3层分解
wname = 'sym4'; lev = 3;
[c,s] = wavedec2(X,lev,wname);
% 图像压缩，使用推荐参数选择各层不同的阈值
alpha = 1.5; m = 2.7*prod(s(1,:));
[thr,nkeep] = wdcbm2(c,s,alpha,m)
% 使用上面的阈值和硬阈值处理方式进行图像压缩
[xd,cxd,sxd,perf0,perfl2] =wdencmp('lvd',c,s,wname,lev,thr,'h');
% 画出原始图像和压缩后的图像
colormap(pink(nbc));
subplot(121), image(wcodemat(X,nbc)),
title('原始图像')
subplot(122), image(wcodemat(xd,nbc)),
title('压缩后的图像')
disp('恢复百分比: ');
num2str(perfl2)
disp('压缩百分比: ');
num2str(perf0)