小波对图像的二次压缩

suffer

1. load flujet;    %装载并显示原始图像
   
2. subplot(1,2,1);
   
3. image(X);
   
4. colormap(map);
   
5. title('原始图像');
   
6. axis square;
   
7. %首先利用小波db3对图像X进行2层分解
   
8. [c,l]=wavedec2(X,2,'db3');
   
9. %全局阈值
   
10. [thr,sorh,keepapp]=ddencmp('cmp','wv',X);
    
11. %进行压缩处理：对所有高频系数进行同样的阈值量化处理
    
12. [Xcmp,cxc,lxc,perf0,perfl2]=wdencmp('gbl',c,l,'db3',2,thr,sorh,keepapp);
    
13. %图示压缩结果
    
14. subplot(1,2,2);
    
15. image(Xcmp);colormap(map);
    
16. title('压缩后的图像');
    
17. axis square;
    
18. disp('小波分解系数中置0的系数个数百分比：');
    
19. perf0
    
20. disp('压缩后保留能量百分比：');
    
21. perfl2

复制代码