%讀原圖?
X=imread（‘cameraman.tif‘）;?

subplot（221）;?
imshow（X）;?
title（‘原圖‘）;?
axis?square;?
sigma=0.01;
M=9;

Y=imnoise（X‘gaussian‘0sigma）;?
subplot（222）;?
imshow（Y）;?
title（‘加噪圖像‘）;?
axis?square;?

%%小波分解
Y=double（Y）;
[CACHCVCD]=?dwt2（Y‘db1?‘）;
%?[cs]=wavedec2（j1‘haar‘）;
y?=[CACH;CVCD];
y=uint8（y）;
subplot（223）;?
title（‘小波分解‘）;
imshow（y）;?

%%對不同的高頻進行維納濾波，可選
y=CV;???????????????
sz=size（y）;
%cacualte?q（ij）
temp=ones（sqrt（M））;
y=double（y）;
temp=double（temp）;
Q=conv2（ytemp‘same‘）/M;
k=1+sqrt（2/M）;
Q（find（Q>k*sigma））=y（find（Q>k*sigma））;
Q（find（Q<=k*sigma?））=0;

Q=double（Q）;
Q=conv2（Qtemp‘same‘）/M;
Q1=Q;
Q（find（Q<0））=0;
a=Q./（Q1+eps）;
xx=a.*y;

%%%重構
XX?=?IDWT2（CACHCVxx‘db1‘）;
XX=uint8（XX）;
subplot（224）;?
imshow（XX）;?
title（‘重構以后的圖像‘）;

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

?????文件????????886??2009-01-04?19:44??小波去噪\tian.m

?????文件?????322095??2008-12-24?20:47??小波去噪\圖像小波域維納濾波去噪算法的改進（1）.pdf

?????目錄??????????0??2009-01-12?19:39??小波去噪

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

???????????????322981????????????????????3