%裝載并圖示原始圖像
load flujet;
subplot(2,2,1);
image(X);
colormap(map);
title('原始圖像');
axis square;
%==============================
%生成含噪圖像并圖示
init=2055615866;
randn('seed',init);
X1=X+20*randn(size(X));
subplot(2,2,2);
image(X1);
colormap(map);
title('含噪圖像');
axis square;
%==============================
%用小波sym2對圖像X1進行一層小波包分解
T=wpdec2(X1,1,'sym2');
%==============================
%設置全局閾值
thr=8.342;
%==============================
%對圖像的小波包分解系數進行軟閾值量化以實現圖像消噪
NT=wpthcoef(T,0,'s',thr);
%==============================
%由于所含的是高頻白噪聲，因此僅對低頻系數進行重構
X2=wprcoef(NT,1);
%==============================
%畫出結果圖像
subplot(2,2,3);
image(X2);
colormap(map);
title('消噪后的圖像');
axis square;