逆谐波均值滤波

Halcom

1. function im = conharmmean_filter(x,m,n,q)
   
2. % 逆谐波均值滤波
   
3. % 函数输入:
   
4. %         x：输入二维图像矩阵
   
5. %         m,n：滤波掩膜尺寸
   
6. % 函数输出:
   
7. %         im：输出图像矩阵，数据类型与输人相同
   
8. if ~isa(x,'double')
   
9.     x = double(x)/255;
   
10. end
    
11. 
    
12. im = imfilter(x.^(q+1),ones(m,n),'replicate');        % 滤波算子
    
13. im = im./(imfilter(x.^q,ones(m,n),'replicate')+eps);   % 逆谐波均值滤波
    
14. im = im2uint8(im);                              % 类型转换
    
15. end
    

复制代码

im1 = conharmmean_filter(im,3,3,3/2);   % 应用逆谐波均值滤波