基于小波变换的分解重构（图像压缩）

Halcom

    采用Bior算子，对图像进行小波分解、重构，对图像的第一层和第二层信息进行压缩，MATLAB程序如下：

1. clc,clear,close all                    % 清屏、清工作区、关闭窗口
   
2. warning off                         % 消除警告
   
3. feature jit off                       % 加速代码执行
   
4. load wbarb;
   
5. subplot(2,2,1);image(X);
   
6. colormap(map);
   
7. xlabel('(a) 原始图像');
   
8. axis square   
   
9. disp('压缩前图像X的大小');       % 显示文字
   
10. whos('X')                        % 显示图像属性
    
11. 
    
12. %对图像用小波进行层小波分解
    
13. [c,s]=wavedec2(X,2,'bior3.7');    % 提取小波分解结构中的一层的低频系数和高频系数
    
14. cal=appcoef2(c,s,'bior3.7',1);    % 水平方向
    
15. ch1=detcoef2('h',c,s,1);         % 垂直方向
    
16. cv1=detcoef2('v',c,s,1);         % 斜线方向
    
17. cd1=detcoef2('d',c,s,1);
    
18. %各频率成份重构
    
19. a1=wrcoef2('a',c,s,'bior3.7',1);
    
20. h1=wrcoef2('h',c,s,'bior3.7',1);
    
21. v1=wrcoef2('v',c,s,'bior3.7',1);
    
22. d1=wrcoef2('d',c,s,'bior3.7',1);
    
23. c1=[a1,h1;v1,d1];
    
24. 
    
25. %显示分频信息
    
26. subplot(2,2,2);image(c1);
    
27. colormap(jet)                     % 设置色彩索引图
    
28. axis square;                      % 设置显示比例
    
29. xlabel ('(b) 分解后低频和高频信息');
    
30. 
    
31. ca1=appcoef2(c,s,'bior3.7',1);
    
32. ca1=wcodemat(ca1,440,'mat',0);
    
33. %改变图像高度并显示
    
34. ca1=0.5*ca1;
    
35. 
    
36. subplot(2,2,3);image(ca1);
    
37. colormap(map);                    % 设置色彩索引图
    
38. axis square;                       % 设置显示比例
    
39. xlabel('(c) 第一次压缩图像');
    
40. disp('第一次压缩图像的大小为：');   % 显示文字
    
41. whos('ca1')                       % 显示图像属性
    
42. 
    
43. %保留小波分解第二层低频信息进行压缩
    
44. ca2=appcoef2(c,s,'bior3.7',2);
    
45. %首先对第二层信息进行量化编码
    
46. ca2=wcodemat(ca2,440,'mat',0);
    
47. %改变图像高度并显示
    
48. ca2=0.25*ca2;
    
49. 
    
50. subplot(2,2,4);image(ca2);  
    
51. colormap(map);                    % 设置色彩索引图
    
52. axis square;                       % 设置显示比例
    
53. xlabel('(d) 第二次压缩图像');  
    
54. disp('第二次压缩图像的大小为：');   % 显示文字
    
55. whos('ca2')                       % 显示图像属性

复制代码