%?Nlayer_denoise.m
clear
clc
%?load?sample?image?and?its?noise?plus?edge?image;
%load?sample_circle_30.mat;
%load?noise_edge_circle.mat;
T=128;
%x=noise_edge;
X=imread（‘source.bmp‘）;
x=imnoise（X‘gaussian‘00.05）;
%?分層閾值處理
for?i=1:6?
????[cs]=wavedec2（xi‘sym6‘）;?%?i層小波分解
????[thr3nkeep3]=wdcbm2（cs2）;%?提取i層閾值
????Duo（::i）=wdencmp（‘lvd‘cs‘sym6‘ithr3‘h‘）;%?i層小波去噪
????sny（i）=10*log10（sum（sum（double（X）.^2））/sum（sum（（double（X）-double（Duo（::i）））.^2）））;
????%?i層小波去噪后信噪比
end
%?顯示各層去噪后圖像
for?i=1:6
????subplot（42i+2）
????imshow（Duo（::i）[]）
end
%?顯示原圖像及加噪后圖像
subplot（421）
imshow（X）
title（‘原圖像‘）
subplot（422）
imshow（x[]）
title（‘添加噪聲后的圖像‘）
%?顯示不同層次去噪后圖像信噪比走勢圖
figure（2）
plot（sny）
axis（[0?6?5?13]）
title（‘不同層次去噪后圖像信噪比走勢圖‘）
????

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件??????66616??2000-07-25?11:21??YuZhiQuZao\Miss.bmp

?????文件????????873??2008-10-06?14:49??YuZhiQuZao\Nlayer_denoise.m

?????文件???????2052??2008-10-06?17:26??YuZhiQuZao\selfAdaptable.m

?????文件????????146??2008-10-06?14:02??YuZhiQuZao\yuzhicl.m

?????文件???????2709??2008-10-06?16:56??YuZhiQuZao\zheZhongYuZhi.m

?????文件????????742??2008-10-06?17:24??YuZhiQuZao\zishy.m

?????文件?????144718??2008-10-06?16:58??YuZhiQuZao\圖像\noise.bmp

?????文件?????144718??2008-10-06?16:58??YuZhiQuZao\圖像\ruanYuZhi.bmp

?????文件?????144718??2008-10-06?16:58??YuZhiQuZao\圖像\yingYuZhi.bmp

?????文件?????144718??2008-10-06?16:58??YuZhiQuZao\圖像\zheZhong.bmp

?????文件?????144718??2008-10-06?17:27??YuZhiQuZao\圖像\ziShiYing.bmp

?????目錄??????????0??2008-10-06?17:27??YuZhiQuZao\圖像

?????目錄??????????0??2008-10-06?17:27??YuZhiQuZao

-----------?---------??----------?-----??----

???????????????796728????????????????????13