小波去噪

Halcom

小波去噪:
[thr,sorh,keepapp] = ddencmp_thr('den','wv',im);  
% ‘gb1’表示使用全局门限进行去噪
% 'sym4'小波变换函数
N = 4; % 小波变换的尺度
im1 = wdencmp('gbl',im,'sym4',N,thr,sorh,keepapp);  % 小波滤波

1. function [thr,sorh,keepapp] = ddencmp_thr(dorc,worwp,x)
   
2. % 函数输入：
   
3. %         dorc: 'den'还是‘cmp’
   
4. %         worwp： 'wv'小波还是‘wp’小波包
   
5. %         x：输入图像二维矩阵
   
6. % 函数输出：
   
7. %         thr：阈值
   
8. %         sorh：软阈值‘s’还是硬阈值'h'
   
9. %         keepapp：常数1；
   
10. % 默认值
    
11. keepapp = 1; % keepapp=1时，表示保持低频图像系数不变，keepapp=0时，表示可以改变。
    
12. if isequal(dorc,'den') && isequal(worwp,'wv')
    
13.     sorh = 's';   % 软阈值门限
    
14. else
    
15.     sorh = 'h';   % 硬阈值门限
    
16. end
    
17. 
    
18. % 最著名的阈值形式
    
19. n = numel(x);  % 一般为65536
    
20. % nominal threshold.
    
21. switch dorc
    
22.   case 'den'
    
23.     switch worwp
    
24.       case 'wv' , thr = sqrt(2*log(n));               % wavelets.
    
25.       case 'wp' , thr = sqrt(2*log(n*log(n)/log(2))); % wavelet packets.
    
26.     end
    
27.   case 'cmp' ,  thr = 1;
    
28. end

复制代码

小波包去噪:
[thr,sorh,keepapp] = ddencmp_thr('den','wp',im);  
im1 = wpdencmp(im,sorh,4,'sym4','threshold',thr,keepapp);