example7_10.m

来自「This complete matlab for neural network」· M 代码 · 共 26 行

%装载并显示原始图像
load flujet;
subplot(1,2,1);
image(X);colormap(map);
title('原始图像');
axis square;
%==============================
%首先利用小波db3对图像X进行2层分解
[c,l]=wavedec2(X,2,'db3');
%==============================
%全局阈值
[thr,sorh,keepapp]=ddencmp('cmp','wv',X);
%==============================
%进行压缩处理：对所有高频系数进行同样的阈值量化处理
[Xcmp,cxc,lxc,perf0,perfl2]=wdencmp('gbl',c,l,'db3',2,thr,sorh,keepapp);
%==============================
%图示压缩结果
subplot(1,2,2);
image(Xcmp);colormap(map);
title('压缩后的图像');
axis square;
disp('小波分解系数中置0的系数个数百分比：');
perf0
disp('压缩后保留能量百分比：');
perfl2