%裝載并圖示原始圖像
load gatlin2;
subplot(2,2,1);
image(X);
colormap(map);
title('原始圖像');
axis square;
%==============================
%生成含噪圖像并圖示
init=2055615866;
randn('seed',init);
XX=X+8*randn(size(X));
subplot(2,2,2);
image(XX);
colormap(map);
title('含噪圖像');
axis square;
%==============================
%對圖像進行消噪處理
%用小波函數coif2對圖像XX進行2層分解
[c,l]=wavedec2(XX,2,'coif2');
n=[1,2];%設置尺度向量
p=[10.28,24.08];%設置閾值向量
%==============================
%對三個高頻系數進行閾值處理
nc=wthcoef2('h',c,l,n,p,'s');
nc=wthcoef2('v',c,l,n,p,'s');
nc=wthcoef2('d',c,l,n,p,'s');
X1=waverec2(nc,l,'coif2');
subplot(2,2,3);
image(X1);
colormap(map);
title('第一次消噪后的圖像');
axis square;
%==============================
%再次對三個高頻系數進行閾值處理
mc=wthcoef2('h',nc,l,n,p,'s');
mc=wthcoef2('v',nc,l,n,p,'s');
mc=wthcoef2('d',nc,l,n,p,'s');
%==========================================
%對更新后的小波分解結構進行重構并圖示結果
X2=waverec2(mc,l,'coif2');
subplot(2,2,4);
image(X2);
colormap(map);
title('第二次消噪后的圖像');
axis square;