巴特沃斯高通滤波器

Halcom

1. function im5 = freqfilter_btw_Hp(im,D0,n)
   
2. % 巴特沃斯高通滤波器
   
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.         if d==0
    
22.             h=0;
    
23.         else
    
24.             h = 1/(1 + 0.414 * (D0 / d)^(2*n));  % 巴特沃斯高通滤波器
    
25.         end
    
26.         result(i,j) = h * im3(i,j);
    
27.     end
    
28. end
    
29. result = ifftshift(result);    % 反中心化
    
30. im4 = ifft2(result);           % 反变换
    
31. im5 = im2uint8(real(im4));     % 滤波图像
    
32. end

复制代码