高斯高通滤波器

Halcom

1. function im5 = freqfilter_gaussian_Hp(im,D0)
   
2. % 高斯高通滤波器gaussian
   
3. % input:
   
4. %     M,N:频域滤波器的尺寸
   
5. %     D0：带阻滤波器的截止频率
   
6. %     n ：阶次
   
7. % output:
   
8. %       H：M x N的矩阵，表示频域滤波器矩阵，数据类型为double，
   
9.     if ~isa(im,'double')
   
10.         im1 = double(im)/255;
    
11.     end
    
12. im2 = fft2(im1);      % 傅里叶变换
    
13. im3 = fftshift(im2);  % 中心化
    
14. 
    
15. [N1, N2] = size(im3);
    
16. n1 = fix(N1 / 2);
    
17. n2 = fix(N2 / 2);
    
18. for i = 1:N1
    
19.     for j = 2:N2
    
20.         D = sqrt((i-n1)^2+(j-n2)^2);
    
21.         h = 1- exp(-D.^2./2./D0./D0);  % 高斯高通滤波器
    
22.         result(i,j) = h * im3(i,j);
    
23.     end
    
24. end
    
25. result = ifftshift(result);    % 反中心化
    
26. im4 = ifft2(result);           % 反变换
    
27. im5 = im2uint8(real(im4));     % 滤波图像
    
28. end

复制代码