www.gusucode.com > 《matlab图像处理与界面编程宝典》秦襄培 编著，每章的MATLAB源代码程序 > 第20章/代码20-2.txt

    
clear;                                   % 清理工作空间                 
load wbarb;                              % 装入图形信号
figure;                                  % 新建窗口
image(X);                                % 显示图像
colormap(map);                           % 设置色彩索引图
axis square;                             % 设置显示比例
title('原始图像');                       % 设置图像标题
disp('压缩前图像的大小');                % 显示文字
whos('X');                               % 显示图像属性
axis square;                             % 设置显示比例
%对图像用db3小波进行二层小波分解
[c,s]=wavedec2(X,5,'db3');
[thr,sorh,keepapp]=ddencmp('cmp','wv',X);
[Xcomp,cxc,lxc,perf0,perfl2]=wdencmp('gbl',c,s,'db3',5,thr,sorh,keepapp);
figure;                                  % 新建窗口
image(Xcomp);                            % 显示图像
colormap(map);                           % 设置色彩索引图
axis square;                             % 设置显示比例
title(' 压缩后的图像');                  % 设置图像标题
disp('压缩后图像的大小');                % 显示文字
whos('Xcomp')                            % 显示图像属性
disp('小波分解系数中值为0的系数个数百分比'); % 显示文字
disp(perf0);                             % 显示数值
disp('压缩后剩余能量百分比');            % 显示文字
disp(perfl2);                            % 显示数值