小波变换-软硬阈值半软阈值图像去噪matlab程序复习进程

1、小 波 变 换 - 软 硬 阈 值 半 软 阈 值 图 像 去 噪 mat l ab 程 序精品文档%软阈值硬阈值半软阈值巴特沃斯滤波%均值为 0 方差 502clcclose allclear allmap=gray256; x=imreadhzw.jpg; x=rgb2grayx; subplot2,3,1; imagex; colormapmap; title 原始图片 ;axis square; init=2055615866; randnseed,init;x1=50.*randnsizex; x=doublex nx=x+x1; subplot2,3,2; imagenx; col

2、ormapmap;title 加噪后的图片 ;axis square;c=num2strc;text100,100,PSNR:text300,100,c;%硬阈值thr,sorh,keepapp=ddencmpnx1=wdencmp gbl ,nx,sym5den ,2,thr, wv h,nx; ,keepapp;subplot2,3,3;imagenx1;title 2.D. ￥ .o .axis squarea1=psnrnx1,x;a1=num2stra1;text100,100,PSNR:;text300,100,a1;%软阈值nx2=wdencmpgbl,nx,sym5,2,thr

3、,s,keepapp;subplot2,3,4;imagenx2;.D. ￥ .o .title axis square;c=psnrnx2,x; PSNR:;c=num2strc;text100,100,text300,100,c;%半软阈值nx3=hsoftnx, subplot2,3,5;sym5,2,0.5,thr;收集于网络,如有侵权请联系治理员删除精品文档imagenx3;.D. ￥ .o .title . axis squarea4=psnrnx3,x; ;PSNR:;a4=num2stra4; text100,100,text300,100,a4;%巴特沃斯g=fft2nx;

4、g=fftshiftg; m,n=sizeg; N=3;d0=60; n1=floorm/2; n2=floorn/2; for i=1:m for j=1:n d=sqrti-n12+j-n22; h=1/1+d/d02*N; gi,j=h*gi,j; end end g=ifftshiftg; g=realifft2g; subplot2,3,6; imageg; title .1 .2 axis square ; a5=psnrg,x; a5=num2stra5; text100,100, PSNR: ; text300,100,a5;function X=hsoftx,wname,n,thr,thrl C,S=wavedec2x,n,wname; dcoef=CprodS1,:+1:end; ind=findabsdcoef=thrl&absdcoefthr+prodS1,:; Cind=signCind.*thr/thr-thrl*absCind-thrl; A1=wrcoef2 a ,C,S,wname,n; H1=wrcoef2 h ,C,S,wname,n; V1=wrcoef2 v