%裝載并圖示原始圖像
load gatlin2;
subplot(2,2,1);
image(X);
colormap(map);
title('原始圖像');

%生成含噪圖像并圖示
init=2055615866;
randn('seed',init);
XX=X+8*randn(size(X));
subplot(2,2,2);
image(XX);
colormap(map);
title('含噪圖像');

%對圖像進行消噪處理
%用小波函數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('第一次消噪后的圖像');

%再次對高頻小波系數進行閾值處理
%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('第二次消噪后的圖像');