mythre.m
上传用户:laurel
上传日期:2013-06-05
资源大小:2k
文件大小:1k
源码类别:

波变换

开发平台:

Matlab

  1. hzsx1=imread('f:/pictures/hzsx1.bmp');
  2. sx=rgb2gray(hzsx1);
  3. subplot(2,2,1),imshow(sx);
  4. subplot(2,2,1),image(sx);
  5. title('原始图像');
  6. axis square;
  7. [c1,l1]=wavedec2(sx,2,'sym4');
  8. a1=wrcoef2('a',c1,l1,'sym4',1);
  9. a2=wrcoef2('a',c1,l1,'sym4',2);
  10. subplot(2,2,2),image(a2);
  11. title('强制消噪处理');
  12. axis square;
  13. [c,l]=wavedec2(sx,5,'sym4');
  14. thr=50;
  15. [sxd,cxd,lxd,perf0,perf12]=wdencmp('gbl',sx,'sym2',2,thr,'h',1);
  16. subplot(2,2,3),image(sxd);
  17. title('全局阈值去噪');
  18. axis square;
  19. thr_h=[75 70];
  20. thr_d=[74 72];
  21. thr_v=[78 68];
  22. thr2=[thr_h;thr_d;thr_v]
  23. [yxd,cxd,lxd,perf0,perf12]=wdencmp('lvd',sx,'sym2',2,thr2,'h');
  24. subplot(2,2,4),image(yxd);
  25. title('给定阈值去噪');
  26. axis square;