小波变换软硬阈值半软阈值图像去噪matlab程序.doc

%软阈值　硬阈值　半软阈值 巴特沃斯滤波
ｃlc
close ａll
cｌｅar　aｌl

maｐ=grａy(256)；　
x＝imｒead(’hzw。jｐg＇); 　
x＝rｇｂ2gray（x)；
subploｔ（２，3,１)；
iｍage(x); 　 　
coloｒｍａｐ（maｐ)；　
title（’原始图片'）；
ａxｉs sｑuａｒｅ; 　
inｉｔ=２05５615866; 　 　
rａnｄn(’sｅｅd',ｉnｉt)；
x１=50.*randn（ｓｉｚe(ｘ))；　　　　％均值为0　方差50^2

x=double（x)
nx＝x+x1；
ｓubplot(２,3,2)；
imａge(nx);
ｃoloｒmap(mａp);
tiｔlｅ（’加噪后的图片');
axｉs ｓqｕare；
c=nｕｍ2str(c);
tｅxt(1００，100,＇PSＮＲ：’)；
tｅxt(３0０，１0０，c）；

％硬阈值
［thｒ,sｏrh，ｋｅｅpapp］=ddencmｐ（'den','wv'，nx)；
ｎx1＝wdenｃmｐ(’gｂl’,ｎx，'ｓym５＇，２，tｈr,'ｈ’，keeｐapp);
suｂｐlot(2,3,３)；
imａgｅ（ｎx1)；
title(’Ó²ãÐֵȥÔëºóͼÏñ');
axiｓ sqｕare；
a１=ｐsｎr(nx１,x);
ａ１＝ｎum２str（a1）;
text(１00,１00，’ＰSNR:');
text(300，10０，a1)；
%软阈值
ｎx2=wｄｅncmp（’gｂl’,ｎｘ,’ｓym5’,2,thr，'s’，ｋeｅpapｐ)；
ｓubpｌｏt（2,3,4）；
imaｇe(ｎx2）；
titlｅ('ÈíãÐֵȥÔëºóͼÏñ’）;
ａxis squaｒｅ;
c=pｓnr（nx2,x); 　　
c=nｕm2ｓｔr(ｃ）;
ｔeｘt(100,１00，＇PSＮR:＇）;
text（300,100,c)；

%半软阈值
nx３=ｈｓoｆt（nx,’ｓｙm５’,２,0．5，ｔhr)；
subpｌｏt(２,3，５);
iｍage(nx3）；
ｔｉtle（’°ëÈíãÐֵȥÔëºóͼƬ＇)；
axｉs sｑｕare;
ａ4＝ｐsｎr（ｎx3，x); 　　
a4=nuｍ２ｓｔｒ(ａ4);
tｅxt(１00,100,'PSNR：'）;
ｔext(３0０，100，a４)；
　
%巴特沃斯
ｇ＝ffｔ２(nx)；
g=fftsｈｉｆt（g）；　
[m,n］=sｉze(g)；
N=３;d0=６0;
ｎ1=flｏｏｒ（ｍ/2)；
ｎ2＝fｌoｏr(n/２）；
for i=1:m
　　 fｏｒ j=１：ｎ
　 　 d=sｑrt((i—n1）^2+（j—n2)^2)；
　 　 h=１/（１+（d/d0)＾（2＊N)）;
　 　　 g(ｉ，ｊ)=h＊g(ｉ，j);
　 end
end
g=ifftｓｈｉft(g