开通博客赚积分
发布资源赚积分
分类

源码开发语言/平台
当前位置：首页 > 源码/资料 > 其他 > 查看源码
denoise2.m
资源名称：小波去噪.rar　[点击查看]
上传用户：sla11nk8
上传日期：2013-03-09
资源大小：21k
文件大小：4k
源码类别：

其他

开发平台：

Matlab

denoise2.m：源码内容
  1. % denoise2.m
  2. %---本程序是基于LeventCode平台作的。
  3. %---包括以下方法：
  4. % BivaShrink方法、模型1、模型2、模型3（TrivaShrink方法）、BayesShrink方法、LAWMLShrink方法的DWT实现和DT_CWT实现。
  5. %---运行方法：
  6. % 将某方法对应行的代码前的“%”去掉，而其他方法前一律都加“%”屏蔽掉。
  7. % 噪声方差大小可以自行设定
  9. I=imread('lena.bmp');%读取图像数据
  10. figure;
  11. imshow(I,256);%显示原图像
  13. %产生噪声图像
  14. theta_noise=20;%噪声方差(可设为其他值）------------------------------------------
  15. noise=GWN2(length(I),theta_noise);
  16. Inoise=double(I)+noise;
  17. figure;
  18. imshow(Inoise,[]);%显示带噪图像
  20. %%%%%%%%%%%%%%%%%%%%%%--------BivaShrink方法---------%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  21. %----------Levent原创Bivariate Shrinkage方法：denoising_dwt-----------------------
  22. Idenoise=denoising_dwt(Inoise);
  23. %---------------------------------------------------------------------------------
  25. %----------复数小波改进-----------------------------------------------------------
  26. %Idenoise=denoising_dtdwt(Inoise);
  27. %---------------------------------------------------------------------------------
  28. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  31. %%%%%%%%%%%%%%%%%%%%%%%--------Cai:模型1方法---------%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  32. %----------Cai:这里用当前系数y1和领域系数y3:BivaShrink13(模型1)-------------------
  33. %Idenoise=BivaShrink13(Inoise);
  34. %---------------------------------------------------------------------------------
  36. %----------复数小波改进-----------------------------------------------------------
  37. %Idenoise=dt_BivaShrink13(Inoise);
  38. %---------------------------------------------------------------------------------
  39. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  42. %%%%%%%%%%%%%%%%%%%%%%%--------Cai:模型2方法---------%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  43. %----------Cai:这里用父层系数y2和邻域系数y3:BivaShrink23--------------------------
  44. %Idenoise=BivaShrink23(Inoise);
  45. %---------------------------------------------------------------------------------
  47. %----------复数小波改进-----------------------------------------------------------
  48. %Idenoise=dt_BivaShrink23(Inoise);
  49. %---------------------------------------------------------------------------------
  50. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  53. %%%%%%%%%%%%%%%%--------Cai:模型3(TrivaShrink方法)---------%%%%%%%%%%%%%%%%%%%%%%%
  54. %----------Cai:这里用当前系数y1、父层系数y2和邻域系数y3---------------------------
  55. %----------是一种三维的萎缩函数:BivaShrink123-------------------------------------
  56. %Idenoise=BivaShrink123(Inoise);
  57. %---------------------------------------------------------------------------------
  59. %----------复数小波改进-----------------------------------------------------------
  60. %Idenoise=dt_BivaShrink123(Inoise);
  61. %---------------------------------------------------------------------------------
  62. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  65. %%%%%%%%%%%%%%%%%%%%%%%%--------BayesShrink---------%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  66. %-------------BayesShrink的方差局部自适应版本－-----------------------------------
  67. %Idenoise=denoising_BayesShrink(Inoise);
  68. %---------------------------------------------------------------------------------
  70. %----------复数小波改进-----------------------------------------------------------
  71. %Idenoise=dt_BayesShrink(Inoise);
  72. %---------------------------------------------------------------------------------
  73. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  76. %%%%%%%%%%%%%%%%%%%%%%%%--------LAWMLShrink---------%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  77. %--------------LAWMLShrink--------------------------------------------------------
  78. %Idenoise=denoising_LAWMLShrink(Inoise);
  79. %---------------------------------------------------------------------------------
  81. %----------复数小波改进-----------------------------------------------------------
  82. %Idenoise=dt_LAWMLShrink(Inoise);
  83. %---------------------------------------------------------------------------------
  84. %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
  86. %输出去噪后的MSE和PSNR值
  87. MSE2_=MSE2(double(I),Idenoise)
  88. PSNR2=PSNR(MSE2_)
  89. figure;
  90. imshow(Idenoise,[]);%显示恢复图像