同态滤波

Halcom

同态滤波
rL = 0.3999;  % 低频增益
rH = 0.71;   % 高频增益
D0 = 1;   % 截止频率

im_e = Homom_filter(im,D0,rL,rH);   % 同态滤波

1. function im_e = Homom_filter(im,d,rL,rH)
   
2. % 同态滤波
   
3. %函数输入:
   
4. %         im:  输入的二维图像矩阵
   
5. %         d:   截止频率
   
6. %         rL   低频增益
   
7. %         rH   高频增益
   
8. %函数输出:
   
9. %         im_e: 重构滤波图像
   
10. 
    
11. if ~isa(im,'double')
    
12.     im = double(im);
    
13. end
    
14. 
    
15. [r c]=size(im); % 输入图像维数
    
16. % 高斯高通滤波
    
17. A=zeros(r,c);
    
18. for i=1:r
    
19.     for j=1:c
    
20.         A(i,j)=(((i-r/2).^2+(j-c/2).^2)).^(.5);
    
21.         B(i,j)=A(i,j)*A(i,j);
    
22.         H(i,j)=(1-exp(-((B(i,j)).^2/d.^2)));  % 传递函数
    
23.     end
    
24. end
    
25. 
    
26. % 同态滤波器传递函数：
    
27. H=((rH-rL).*H)+rL;
    
28. % 取对数
    
29. im_l=log2(im + 1e-5);
    
30. % 离散傅里叶变换
    
31. im_f=fft2(im_l);
    
32. % 滤波
    
33. im_nf=H.*im_f;
    
34. % DFT反变换
    
35. im_n=abs(ifft2(im_nf));
    
36. % 指数变换，消除取对数
    
37. im_e = exp(im_n);  % 滤波矩阵
    
38. im_e = uint8(im_e);

复制代码