matlab小波分析在图像消噪中的应用实例

ttksd

　　例一：利用小波分析给定一个二维含噪图像进行消噪处理
　　源程序如下：

1. 　　clear all; %清楚所有变量
   
2. 
   
3. 　　RGB=imread('5.bmp'); %读取MATLAB目录下的测试图片 名为：5.bmp
   
4. 
   
5. 　　[x,map]=rgb2ind(RGB,128); %把真彩色图像RGB转换为索引图像
   
6. 
   
7. 　　subplot(2,2,1);
   
8. 
   
9. 　　image(x);
   
10. 
    
11. 　　colormap(map);
    
12. 
    
13. 　　title('原始的图像');
    
14. 
    
15. 　　init=2055615866; %加入噪声
    
16. 
    
17. 　　rand('seed',init);
    
18. 
    
19. 　　X=double(x);
    
20. 
    
21. 　　XX=X+randn(size(x))/10;
    
22. 
    
23. 　　subplot(2,2,2);
    
24. 
    
25. 　　image(XX);
    
26. 
    
27. 　　colormap(map);
    
28. 
    
29. 　　title('含噪声的图像');
    
30. 
    
31. 　　[c,l]=wavedec2(XX,2,'sym5'); %进行小波分析
    
32. 
    
33. 　　a1=wrcoef2('a',c,l,'sym5',1);
    
34. 
    
35. 　　a2=wrcoef2('a',c,l,'sym5',2);
    
36. 
    
37. 　　subplot(2,2,3);
    
38. 
    
39. 　　image(a1);
    
40. 
    
41. 　　colormap(map);
    
42. 
    
43. 　　title('第一层的重构图像');
    
44. 
    
45. 　　subplot(2,2,4);
    
46. 
    
47. 　　image(a2);
    
48. 
    
49. 　　colormap(map);
    
50. 
    
51. 　　title('第二层的重构图像');
    

复制代码

　　试验结果的图片：


　　例一：利用小波分析给定一个二维含噪图像进行消噪处理例二：利用二维小波变换给定图像进行小波消噪处理

　　例二：利用二维小波变换给定图像进行小波消噪处理

1. 　　clear all
   
2. 
   
3. 　　RGB=imread('5.bmp');
   
4. 
   
5. 　　[x,map]=rgb2ind(RGB,128);
   
6. 
   
7. 　　subplot(2,2,1);
   
8. 
   
9. 　　image(x);
   
10. 
    
11. 　　colormap(map);
    
12. 
    
13. 　　title('原始图像');
    
14. 
    
15. 　　init=2055615866;
    
16. 
    
17. 　　rand('seed',init);
    
18. 
    
19. 　　X=double(x);
    
20. 
    
21. 　　XX=X+randn(size(x))/10;
    
22. 
    
23. 　　subplot(2,2,2);
    
24. 
    
25. 　　image(XX);
    
26. 
    
27. 　　colormap(map);
    
28. 
    
29. 　　title('含噪图像');
    
30. 
    
31. 　　[c,l]=wavedec2(XX,3,'coif2');
    
32. 
    
33. 　　n=[1,2];
    
34. 
    
35. 　　p=[10.28,24.08];
    
36. 
    
37. 　　%nc=wthcoef2('h',c,l,n,p,'s');
    
38. 
    
39. 　　%nc=wthcoef2('v',c,l,n,p,'s');
    
40. 
    
41. 　　nc=wthcoef2('d',c,l,n,p,'s');
    
42. 
    
43. 　　X1=waverec2(nc,l,'coif2');
    
44. 
    
45. 　　subplot(2,2,3);
    
46. 
    
47. 　　image(X1);
    
48. 
    
49. 　　colormap(map);
    
50. 
    
51. 　　%mc=wthcoef2('h',nc,l,n,p,'s');
    
52. 
    
53. 　　mc=wthcoef2('v',nc,l,n,p,'s');
    
54. 
    
55. 　　%mc=wthcoef2('d',nc,l,n,p,'s');
    
56. 
    
57. 　　X2=waverec2(mc,l,'coif2');
    
58. 
    
59. 　　subplot(2,2,4);
    
60. 
    
61. 　　image(X2);
    
62. 
    
63. 　　colormap(map);
    
64. 
    
65. 　　title('第二次消噪后的图像');
    

复制代码

转自：http://www.zhixing123.cn/matlab/23150.html