巴特沃斯低通滤波器

Halcom

1. function im5 = freqfilter_btw_lp(im,D0,n)
   
2.     if ~isa(im,'double')
   
3.         im1 = double(im)/255;
   
4.     end
   
5. im2 = fft2(im1);      % 傅里叶变换
   
6. im3 = fftshift(im2);  % 中心化
   
7. 
   
8. [N1, N2] = size(im3);
   
9. n1 = fix(N1 / 2);
   
10. n2 = fix(N2 / 2);
    
11. for i = 1:N1
    
12.     for j = 2:N2
    
13.         d = sqrt((i-n1)^2+(j-n2)^2);
    
14.         h = 1/(1 + 0.414 * (d / D0)^(2*n));  % 巴特沃斯低通滤波器
    
15.         result(i,j) = h * im3(i,j);
    
16.     end
    
17. end
    
18. result = ifftshift(result);    % 反中心化
    
19. im4 = ifft2(result);           % 反变换
    
20. im5 = im2uint8(real(im4));     % 滤波图像
    
21. 
    
22. end

复制代码