% 轉載原始圖像
load belmont2;
% X包含原始圖像
x=X(100:200,100:200);
nbc=size(map,1);
subplot(221),image(x),colormap(map);
axis square;
title('原始圖像');
% 使用wdencmp進行圖像壓縮
% 得到X的小波分解結構
n=5;w='sym2';
[c,l]=wavedec2(x,n,w);
% 小波系數的閾值處理
thr=20;
[xd,cxd,lxd,perf0,perfl2]=wdencmp('gbl',c,l,w,n,thr,'h',l);
% 畫出壓縮后的圖像
subplot(222),image(xd),colormap(map);
title('全局閾值壓縮圖像');
axis square;
subplot(223),image(x),colormap(map);
axis square;
title('原始圖像');
% 定義三個方向的各層不同的閾值
thr_h=[17 18];
thr_d=[19 20];
thr_v=[21 22];
thr=[thr_h;thr_d;thr_v];
% 進行圖像壓縮
[xd,cxd,lxd,perf0,perfl2]=wdencmp('lvd',x,'sym8',2,thr,'h');
% 畫出壓縮后的圖像
subplot(224),image(xd),colormap(map);
title('獨立閾值壓縮圖像');
axis square;