% 當前延拓模式是補零
% 裝載原始圖像
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',1);
% 畫出壓縮后圖像
subplot(222),image(xd),colormap(map);
title('全局閾值壓縮圖像');
axis square;
% 定義三個方向各層不同的閾值
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(223),image(xd),colormap(map);
title('獨立閾值壓縮后圖像');
axis square;