小波變換-軟硬閾值半軟閾值圖像去噪matlab程序

小波變換-軟硬閾值半軟閾值圖像去噪matlab程序小波變換-軟硬閾值半軟閾值圖像去噪matlab程序小波變換-軟硬閾值半軟閾值圖像去噪matlab程序xxx公司小波變換-軟硬閾值半軟閾值圖像去噪matlab程序文件編號：文件日期：修訂次數(shù)：第1.0次更改批準審核制定方案設(shè)計，管理制度%軟閾值硬閾值半軟閾值巴特沃斯濾波clccloseallclearallmap=gray(256);x=imread('');x=rgb2gray(x);subplot(2,3,1);image(x);colormap(map);title('原始圖片');axissquare;init=66;randn('seed',init);x1=50.*randn(size(x));%均值為0方差50^2x=double(x)nx=x+x1;subplot(2,3,2);image(nx);colormap(map);title('加噪后的圖片');axissquare;c=num2str(c);text(100,100,'PSNR:');text(300,100,c);%硬閾值[thr,sorh,keepapp]=ddencmp('den','wv',nx);nx1=wdencmp('gbl',nx,'sym5',2,thr,'h',keepapp);subplot(2,3,3);image(nx1);title('ó2?D?μè￥??oóí???');axissquare;a1=psnr(nx1,x);a1=num2str(a1);text(100,100,'PSNR:');text(300,100,a1);%軟閾值nx2=wdencmp('gbl',nx,'sym5',2,thr,'s',keepapp);subplot(2,3,4);image(nx2);title('èí?D?μè￥??oóí???');axissquare;c=psnr(nx2,x);c=num2str(c);text(100,100,'PSNR:');text(300,100,c);%半軟閾值nx3=hsoft(nx,'sym5',2,,thr);subplot(2,3,5);image(nx3);title('°?èí?D?μè￥??oóí???');axissquare;a4=psnr(nx3,x);a4=num2str(a4);text(100,100,'PSNR:');text(300,100,a4);%巴特沃斯g=fft2(nx);g=fftshift(g);[m,n]=size(g);N=3;d0=60;n1=floor(m/2);n2=floor(n/2);fori=1:mforj=1:nd=sqrt((i-n1)^2+(j-n2)^2);h=1/(1+(d/d0)^(2*N));g(i,j)=h*g(i,j);endendg=ifftshift(g);g=(real(ifft2(g)));subplot(2,3,6);image(g);title('°íì????1μíí¨??2¨');axissquare;a5=psnr(g,x);a5=num2str(a5);text(100,100,'PSNR:');text(300,100,a5);functionX=hsoft(x,wname,n,thr,thrl)[C,S]=wavedec2(x,n,wname);dcoef=C(prod(S(1,:))+1:end);ind=find(abs(dcoef)<thrl)+prod(S(1,:));C(ind)=0;ind=find(abs(dcoef)>=thrl&abs(dcoef)<thr)+prod(S(1,:));C(ind)=sign(C(ind)).*((thr/(thr-thrl)*(abs(C(ind))-thrl)));A1=wrcoef2('a',C,S,wname,n);H1