小波變換-軟硬閾值半軟閾值圖像去噪matlab程序復(fù)習(xí)進(jìn)程

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;收集于網(wǎng)絡(luò),如有侵權(quán)請(qǐng)聯(lián)系治理員刪除精品文檔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