中值滤波

Halcom

1. function b = med_filt2(a,mn)
   
2. % 中值滤波
   
3. % 输入：
   
4. %      a：输入二维图像局长呢
   
5. %      mn：为[m,n]滤波模板
   
6. % 输出：
   
7. %      b:中值滤波图像
   
8. 
   
9. domain = ones(mn);          % 模板
   
10. if (rem(prod(mn), 2) == 1)
    
11.     order = (prod(mn)+1)/2; % 中值
    
12.     b = ordfilt2(a, order, domain, 'zeros'); % 排序滤波
    
13. else
    
14.     order1 = prod(mn)/2;
    
15.     order2 = order1+1;
    
16.     b = ordfilt2(a, order1, domain, 'zeros');  % 排序滤波
    
17.     b2 = ordfilt2(a, order2, domain, 'zeros'); % 排序滤波
    
18.         if islogical(b)  % 逻辑运算
    
19.                 b = b | b2;  % 或运算
    
20.         else
    
21.                 b =        imlincomb(0.5, b, 0.5, b2); % 0.5*b + 0.5*b2
    
22.         end
    
23. end
    

复制代码