%裝載并顯示原始圖像
load belmont2;
subplot(1,2,1);image(X);colormap(map);title('原始圖像');
%首先利用db3小波對(duì)圖像X進(jìn)行2層分解
[c,l]=wavedec2(X,2,'db3');
%全局閾值
[thr,sorh,keepapp]=ddencmp('cmp','wv',X);
%壓縮處理：對(duì)所有高頻系數(shù)進(jìn)行同樣的閾值量化處理
[Xcmp,cxc,lxc,perf0,perfl2]=wdencmp('gbl',c,l,'db3',2,thr,sorh,keepapp);
%將壓縮后的圖像與原始圖像相比較
subplot(1,2,2);
image(Xcmp);colormap(map);
title('壓縮后的圖像');
%顯示相關(guān)參數(shù)
disp('小波分解系數(shù)中為0的系數(shù)個(gè)數(shù)百分比：');
perf0
disp('壓縮后保留能量百分比：');
perfl2