kuwahara_filter滤波器

Halcom

kuwahara_filter滤波器
im1 = kuwahara(im,5);

1. function [Y,Xpad] = kuwahara(X,WINSZ)
   
2. % kuwahara_filter滤波器
   
3. % 对图像边界轮廓有较强的保护作用
   
4. % kuwahara nonlinear edge-preserving filtering
   
5. % 函数输入：
   
6. %          X：二维图像矩阵
   
7. %          WINSZ： window size
   
8. % 函数输出：
   
9. %          Y：滤波图像
   
10. %          Xpad：点扩展矩阵
    
11. % The Kuwahara filter 4块，(最直观的如下所示，5x5 pixels).
    
12. %
    
13. %    ( a  a  ab   b  b)
    
14. %    ( a  a  ab   b  b)
    
15. %    (ac ac abcd bd bd)
    
16. %    ( c  c  cd   d  d)
    
17. %    ( c  c  cd   d  d)
    
18. 
    
19. if nargin < 2
    
20.     WINSZ = 5;  % 默认值
    
21. end
    
22. 
    
23. if ~isa(X,'double')
    
24.     X = im2double(X);  % double类型
    
25. end
    
26. 
    
27. Pad = floor(WINSZ/2);  % 向下取整
    
28. 
    
29. Xpad = padarray(X,[Pad, Pad],'replicate');  % 点扩展
    
30. % A = [ 1 2;
    
31. %       3 4 ];
    
32. % B = padarray(A,[3 2],'replicate','post')
    
33. % B =
    
34. %      1     2     2     2
    
35. %      3     4     4     4
    
36. %      3     4     4     4
    
37. %      3     4     4     4
    
38. %      3     4     4     4
    
39. [padRows,padCols] = size(Xpad);                    % 求维数
    
40. Y = zeros(size(X));
    
41. 
    
42. nRowIters = length((Pad+1):(padRows-Pad));
    
43. count = 1;
    
44. for i = (Pad+1):(padRows-Pad)
    
45.     for j = (Pad+1):(padCols-Pad)
    
46.         % 分成每一个小窗，Q1~Q4
    
47.         w = Xpad((i-Pad):(i+Pad),(j-Pad):(j+Pad));
    
48.         wnw = w(1:(Pad+1),1:(Pad+1));
    
49.         wne = w(1:(Pad+1),(Pad+1):WINSZ);
    
50.         wsw = w((Pad+1):WINSZ,1:(Pad+1));
    
51.         wse = w((Pad+1):WINSZ,(Pad+1):WINSZ);
    
52.         % 计算方差
    
53.         s = var([wnw(:),wne(:),wsw(:),wse(:)]);
    
54.         m = mean([wnw(:),wne(:),wsw(:),wse(:)]);
    
55.         [y,k] = min(s);  
    
56.         % 计算均值
    
57.         Y(i,j) = m(k);
    
58.     end
    
59. 
    
60. end

复制代码