% 裝載原始圖像
load belmont2 
nbc = size(map,1);
% 使用coif2執行圖像的3層小波分解
wname = 'coif2'; lev = 3;
[c,s] = wavedec2(X,lev,wname);
% 由第1層的高頻系數估計噪聲標準差
det1 = detcoef2('compact',c,s,1);
sigma = median(abs(det1))/0.6745;
% 圖像降噪時，使用wbmpen函數選擇全局閾值
alpha = 1.2;
thr = wbmpen(c,l,sigma,alpha)
% 使用軟閾值和保存的低頻信號，進行圖像降噪
keepapp = 1;
xd = wdencmp('gbl',c,s,wname,lev,thr,'s',keepapp);
% 畫出原始圖像和降噪圖像
colormap(pink(nbc));
subplot(121), image(wcodemat(X,nbc))
title('原始圖像')
subplot(122), image(wcodemat(xd,nbc))
title('降噪圖像')