谐波均值滤波器

Halcom

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

复制代码

im1 = harmonymean_filter(im,3,3);  % 应用谐波均值滤波