六抽头滤波

Halcom

六抽头滤波:
filter_coef = [1 -5 20 20 -5 1]/32;  %(6-tap filter)系数
im_inter = six_tap_filter(im, filter_coef); % 六抽头滤波

1. function A_interp_pic = six_tap_filter(A, filter_coef)
   
2. % 六抽头插值滤波
   
3. % 函数输入:    A: 输入图像
   
4. %             filter_coef: 滤波器系数
   
5. % 函数输出:   
   
6. %             A_interp_pic:  六抽头插值滤波图像
   
7. 
   
8. if (length(size(A)) == 3)
   
9. % 如果输入图像为3D数组，则重复插值滤波3次
   
10.     for i=1:3
    
11.         A_interp_pic(:,:,i) = Bilinear_Filter_interp(A(:,:,i),filter_coef);
    
12.     end
    
13.    
    
14. else
    
15.     [m,n] = size(A);       % 求行列
    
16.     % 初始化
    
17.     A1 = [];           % 初始化
    
18.     A_interp_pic = []; % 初始化
    
19.     % 列插值
    
20.     A_col= filter2(filter_coef,A);       % 滤波
    
21.     for i=1:n
    
22.         A1 = [A1,A(:,i),A_col(:,i)];
    
23.     end
    
24.     A1(:,end) = [];    % 去边缘
    
25.     % 行插值
    
26.     A_rows = filter2(filter_coef,A1')';  % 滤波
    
27.     for i=1:m
    
28.         A_interp_pic = [A_interp_pic; A1(i,:); A_rows(i,:)];
    
29.     end
    
30.     A_interp_pic(end,:) = []; % 去边缘
    
31.    
    
32. end

复制代码