资源描述:
《小波去噪主程序资料.doc》由会员上传分享,免费在线阅读,更多相关内容在教育资源-天天文库。
1、小波去噪主程序%xiaoboquzao.m%对彩色图像进行去噪I=imread('菊花.png','png');%读入图像X=im2double(I);%转换成双精度类型x_noise=imnoise(X,'gaussian',0.01);%加入高斯噪声%提取三个通道信息xr=x_noise(:,:,1);%R通道xg=x_noise(:,:,2);%G通道xb=x_noise(:,:,3);%B通道%估计三个通道的阈值[Cr,Sr]=wavedec2(xr,2,'sym4');[Cg,Sg]=wavedec2(xg,2,'sym
2、4');[Cb,Sb]=wavedec2(xb,2,'sym4');thr_r=Donoho(xr);%R通道全局阈值thr_g=Donoho(xg);%G通道全局阈值thr_b=Donoho(xb);%B通道全局阈值%对三个通道分别进行去噪%Donoho全局阈值软阈值公式x_soft_r=wdenoise(xr,'gbl','s',thr_r,'sym4',2);x_soft_g=wdenoise(xg,'gbl','s',thr_g,'sym4',2);x_soft_b=wdenoise(xb,'gbl','s',thr_b,
3、'sym4',2);%Donoho全局阈值硬阈值公式----------------------------------------------x_hard_r=wdenoise(xr,'gbl','h',thr_r,'sym4',2);x_hard_g=wdenoise(xg,'gbl','h',thr_g,'sym4',2);x_hard_b=wdenoise(xb,'gbl','h',thr_b,'sym4',2);%Birge-Massart策略软阈值公式----------------------------------
4、------------thr_lvd_r=thr_rx_soft_lvd_r=wdenoise(xr,'lvd','s',thr_lvd_r,'sym4',1);thr_lvd_g=thr_gx_soft_lvd_g=wdenoise(xg,'lvd','s',thr_lvd_g,'sym4',1);thr_lvd_b=thr_bx_soft_lvd_b=wdenoise(xb,'lvd','s',thr_lvd_b,'sym4',1);%半软阈值-----------------------------------------
5、----------------------x1_r=den1(xr,'sym4',2,thr_r);x1_g=den1(xg,'sym4',2,thr_g);x1_b=den1(xb,'sym4',2,thr_b);%半软阈值+均值滤波----------------------------------------------------x1_5_r=den1_5_1(xr,'sym4',2,thr_r,0.5*thr_r);x1_5_g=den1_5_1(xg,'sym4',2,thr_g,0.5*thr_g);x1_5_b=
6、den1_5_1(xb,'sym4',2,thr_b,0.5*thr_b);%自适应阈值-------------------------------------------------------------x4_r=den4(xr,'sym4',2);x4_g=den4(xg,'sym4',2);x4_b=den4(xb,'sym4',2);%恢复去噪后的图像x_soft=cat(3,x_soft_r,x_soft_g,x_soft_b);%Donoho软阈值x_hard=cat(3,x_hard_r,x_hard_g,x_h
7、ard_b);%Donoho硬阈值x1=cat(3,x1_r,x1_g,x1_b);%半软阈值x4=cat(3,x4_r,x4_g,x4_b);%自适应阈值x1_5=cat(3,x1_5_r,x1_5_g,x1_5_b);%半软阈值+均值滤波%计算去噪图像与原图像峰值信噪比psnr_soft=PSNR_color(x_soft,X)psnr_hard=PSNR_color(x_hard,X)psnr1=PSNR_color(x1,X)psnr1_5=PSNR_color(x1_5,X)psnr4=PSNR_color(x4,X)%
8、显示去噪后的图像figure;imshow(X);title('原图像');figure;imshow(x_noise);title('带噪声图像');figure;imshow(x_soft);title('Donoho软阈值')figure;i