噪声滤波方法均值滤波

一、理论与概念讲解

1、非线性滤波概述

之前我们说的线性滤波，即两个信号之和的响应和他们各自响应之和相等。换句话说，每个像素的输出值是一些输入像素的加权和，线性滤波器易于构造，并且易于从频率响应角度来进行分析。

　　其实在很多情况下，使用邻域像素的非线性滤波也许会得到更好的效果。比如在噪声是散粒噪声而不是高斯噪声，即图像偶尔会出现很大的值的时候。在这种情况下，用高斯滤波器对图像进行模糊的话，噪声像素是不会被去除的，它们只是转换为更为柔和但仍然可见的散粒。这就到了中值滤波登场的时候了。

2、中值滤波

　　中值滤波（Median filter）是一种典型的非线性滤波技术，基本思想是用像素点邻域灰度值的中值来代替该像素点的灰度值，该方法在去除脉冲噪声、椒盐噪声的同时又能保留图像边缘细节。

　　中值滤波是基于排序统计理论的一种能有效抑制噪声的非线性信号处理技术，其基本原理是把数字图像或数字序列中一点的值用该点的一个邻域中各点值的中值代替，让周围的像素值接近的真实值，从而消除孤立的噪声点，对于斑点噪声（speckle noise）和椒盐噪声（salt-and-pepper noise）来说尤其有用，因为它不依赖于邻域内那些与典型值差别很大的值。中值滤波器在处理连续图像窗函数时与线性滤波器的工作方式类似，但滤波过程却不再是加权运算。

中值滤波在一定的条件下可以克服常见线性滤波器如最小均方滤波、方框滤波器、均值滤波等带来的图像细节模糊，而且对滤除脉冲干扰及图像扫描噪声非常有效，也常用于保护边缘信息, 保存边缘的特性使它在不希望出现边缘模糊的场合也很有用，是非常经典的平滑噪声处理方法。

中值滤波与均值滤波器比较

中值滤波器与均值滤波器比较的优势：在均值滤波器中，由于噪声成分被放入平均计算中，所以输出受到了噪声的影响，但是在中值滤波器中，由于噪声成分很难选上，所以几乎不会影响到输出。因此同样用3x3区域进行处理，中值滤波消除的噪声能力更胜一筹。中值滤波无论是在消除噪声还是保存边缘方面都是一个不错的方法。 

中值滤波器与均值滤波器比较的劣势：中值滤波花费的时间是均值滤波的5倍以上。

　　顾名思义，中值滤波选择每个像素的邻域像素中的中值作为输出，或者说中值滤波将每一像素点的灰度值设置为该点某邻域窗口内的所有像素点灰度值的中值。

例如，取3 x 3的函数窗，计算以点[i,j]为中心的函数窗像素中值步骤如下：

1      按强度值大小排列像素点．
2      选择排序像素集的中间值作为点[i,j]的新值．

这一过程如图下图所示．

　　一般采用奇数点的邻域来计算中值，但如果像素点数为偶数时，中值就取排序像素中间两点的平均值。

但是对一些细节多，特别是线、尖顶等细节多的图像不宜采用中值滤波。

3、双边滤波

　　双边滤波（Bilateral filter）是一种非线性的滤波方法，是结合图像的空间邻近度和像素值相似度的一种折衷处理，同时考虑空域信息和灰度相似性，达到保边去噪的目的。具有简单、非迭代、局部的特点。

双边滤波器的好处是可以做边缘保存（edge preserving），一般过去用的维纳滤波或者高斯滤波去降噪，都会较明显地模糊边缘，对于高频细节的保护效果并不明显。双边滤波器顾名思义比高斯滤波多了一个高斯方差sigma－d，它是基于空间分布的高斯滤波函数，所以在边缘附近，离的较远的像素不会太多影响到边缘上的像素值，这样就保证了边缘附近像素值的保存。但是由于保存了过多的高频信息，对于彩色图像里的高频噪声，双边滤波器不能够干净的滤掉，只能够对于低频信息进行较好的滤波。

在双边滤波器中，输出像素的值依赖于邻域像素值的加权值组合：

 

而加权系数w(i,j,k,l)取决于定义域核和值域核的乘积。其中定义域核表示如下（如图）：

 

定义域滤波对应图示：

值域核表示为：

 

值域滤波：

两者相乘后，就会产生依赖于数据的双边滤波权重函数：

 

二、示例

1、中值滤波——medianBlur函数

　medianBlur函数使用中值滤波器来平滑（模糊）处理一张图片，从src输入，而结果从dst输出。且对于多通道图片，每一个通道都单独进行处理，并且支持就地操作（In-placeoperation）。函数原型如下：

1 C++: void medianBlur(InputArray src,OutputArray dst, int ksize)

参数详解：

1     第一个参数，InputArray类型的src，函数的输入参数，填1、3或者4通道的Mat类型的图像；当ksize为3或者5的时候，图像深度需为CV_8U，CV_16U，或CV_32F其中之一，而对于较大孔径尺寸的图片，它只能是CV_8U。
2     第二个参数，OutputArray类型的dst，即目标图像，函数的输出参数，需要和源图片有一样的尺寸和类型。我们可以用Mat::Clone，以源图片为模板，来初始化得到如假包换的目标图。
3     第三个参数，int类型的ksize，孔径的线性尺寸（aperture linear size），注意这个参数必须是大于1的奇数，比如：3，5，7，9 ...

示例：

1        //载入原图
2        Mat image=imread("1.jpg");
3        //进行中值滤波操作
4        Mat out;
5        medianBlur( image, out, 7);

完整示例：

1 #include "opencv2/core/core.hpp"
 2 #include"opencv2/highgui/highgui.hpp"
 3 #include"opencv2/imgproc/imgproc.hpp"
 4 #include <iostream>
 5 
 6 using namespace std;
 7 using namespace cv;
 8 
 9 int main( )
10 {
11        //载入原图
12        Mat image=imread("1.jpg");
13        if(!image.data)
14        {
15            cout << "载入错误！" << endl;
16            return false;
17        }
18 
19        //创建窗口
20        namedWindow("中值滤波【原图】" );
21        namedWindow("中值滤波【效果图】");
22 
23        //显示原图
24        imshow("中值滤波【原图】", image );
25 
26        //进行中值滤波操作
27        Mat out;
28        medianBlur( image, out, 7);
29 
30        //显示效果图
31        imshow("中值滤波【效果图】" ,out );
32        imwrite("效果图.jpg", out);
33 
34        waitKey(0 );
35 }

效果图：

2、双边滤波——bilateralFilter函数

用双边滤波器来处理一张图片，由src输入图片，结果于dst输出。函数原型如下

1 C++: void bilateralFilter(InputArray src, OutputArraydst, int d, double sigmaColor, double sigmaSpace, int borderType=BORDER_DEFAULT)

1     第一个参数，InputArray类型的src，输入图像，即源图像，需要为8位或者浮点型单通道、三通道的图像。
2     第二个参数，OutputArray类型的dst，即目标图像，需要和源图片有一样的尺寸和类型。
3     第三个参数，int类型的d，表示在过滤过程中每个像素邻域的直径。如果这个值我们设其为非正数，那么OpenCV会从第五个参数sigmaSpace来计算出它来。
4     第四个参数，double类型的sigmaColor，颜色空间滤波器的sigma值。这个参数的值越大，就表明该像素邻域内有更宽广的颜色会被混合到一起，产生较大的半相等颜色区域。
5     第五个参数，double类型的sigmaSpace坐标空间中滤波器的sigma值，坐标空间的标注方差。他的数值越大，意味着越远的像素会相互影响，从而使更大的区域足够相似的颜色获取相同的颜色。当d>0，d指定了邻域大小且与sigmaSpace无关。否则，d正比于sigmaSpace。
6     第六个参数，int类型的borderType，用于推断图像外部像素的某种边界模式。注意它有默认值BORDER_DEFAULT。

示例：

1 //载入原图
2 Mat image=imread("1.jpg");
3 //进行双边滤波操作
4 Mat out;
5 bilateralFilter( image, out, 25, 25*2, 25/2 );

完整示例：

1 #include "opencv2/core/core.hpp"
 2 #include "opencv2/highgui/highgui.hpp"
 3 #include "opencv2/imgproc/imgproc.hpp"
 4 #include <iostream>
 5 
 6 using namespace std;
 7 using namespace cv;
 8 
 9 int main( )
10 {
11     //载入原图
12     Mat image=imread("1.jpg");
13     if(!image.data)
14     {
15        cout << "载入出错！" << endl;
16        return false;
17     }
18 
19     //创建窗口
20     namedWindow("双边滤波【原图】" );
21     namedWindow("双边滤波【效果图】");
22 
23     //显示原图
24     imshow("双边滤波【原图】", image );
25 
26     //进行双边滤波操作
27     Mat out;
28     bilateralFilter( image, out, 25, 25*2, 25/2 );
29 
30     //显示效果图
31     imshow("双边滤波【效果图】" ,out );
32     imwrite("效果图.jpg", out);
33 
34     waitKey(0 );
35 }

效果图：

综合示例（涵盖放宽滤波、均值滤波、高斯滤波、中值滤波和双边滤波）：

1 #include <opencv2/core/core.hpp>
  2 #include<opencv2/highgui/highgui.hpp>
  3 #include<opencv2/imgproc/imgproc.hpp>
  4 #include <iostream>
  5 
  6 using namespace std;
  7 using namespace cv;
  8 
  9 Mat g_srcImage,g_dstImage1,g_dstImage2,g_dstImage3,g_dstImage4,g_dstImage5;
 10 int g_nBoxFilterValue=6;  //方框滤波内核值
 11 int g_nMeanBlurValue=10;  //均值滤波内核值
 12 int g_nGaussianBlurValue=6;  //高斯滤波内核值
 13 int g_nMedianBlurValue=10;  //中值滤波参数值
 14 int g_nBilateralFilterValue=10;  //双边滤波参数值
 15 
 16 //轨迹条回调函数
 17 static void on_BoxFilter(int, void *);            //方框滤波
 18 static void on_MeanBlur(int, void *);           //均值块滤波器
 19 static void on_GaussianBlur(int, void *);                    //高斯滤波器
 20 static void on_MedianBlur(int, void *);               //中值滤波器
 21 static void on_BilateralFilter(int, void*);                    //双边滤波器
 22 
 23 int main(  )
 24 {
 25     //载入原图
 26     g_srcImage= imread( "1.jpg", 1 );
 27     if(!g_srcImage.data )
 28     {
 29        cout << "载入出错！" << endl;
 30        return false;
 31     }
 32 
 33     //克隆原图到四个Mat类型中
 34     g_dstImage1= g_srcImage.clone( );
 35     g_dstImage2= g_srcImage.clone( );
 36     g_dstImage3= g_srcImage.clone( );
 37     g_dstImage4= g_srcImage.clone( );
 38     g_dstImage5= g_srcImage.clone( );
 39 
 40     //显示原图
 41     namedWindow("【<0>原图窗口】", 1);
 42     imshow("【<0>原图窗口】",g_srcImage);
 43 
 44 
 45     //=================【<1>方框滤波】=========================
 46     //创建窗口
 47     namedWindow("【<1>方框滤波】", 1);
 48     //创建轨迹条
 49     createTrackbar("内核值：", "【<1>方框滤波】",&g_nBoxFilterValue, 50,on_BoxFilter );
 50     on_MeanBlur(g_nBoxFilterValue,0);
 51     imshow("【<1>方框滤波】", g_dstImage1);
 52     imwrite("方框.jpg", g_dstImage1);
 53 
 54 
 55     //=================【<2>均值滤波】==========================
 56     //创建窗口
 57     namedWindow("【<2>均值滤波】", 1);
 58     //创建轨迹条
 59     createTrackbar("内核值：", "【<2>均值滤波】",&g_nMeanBlurValue, 50,on_MeanBlur );
 60     on_MeanBlur(g_nMeanBlurValue,0);
 61     imshow("【<2>均值滤波】",g_dstImage2);
 62     imwrite("均值.jpg", g_dstImage2);
 63 
 64     //=================【<3>高斯滤波】===========================
 65     //创建窗口
 66     namedWindow("【<3>高斯滤波】", 1);
 67     //创建轨迹条
 68     createTrackbar("内核值：", "【<3>高斯滤波】",&g_nGaussianBlurValue, 50,on_GaussianBlur );
 69     on_GaussianBlur(g_nGaussianBlurValue,0);
 70     imshow("<3>高斯滤波", g_dstImage3);
 71     imwrite("高斯.jpg", g_dstImage3);
 72 
 73     //=================【<4>中值滤波】===========================
 74     //创建窗口
 75     namedWindow("【<4>中值滤波】", 1);
 76     //创建轨迹条
 77     createTrackbar("参数值：", "【<4>中值滤波】",&g_nMedianBlurValue, 50,on_MedianBlur );
 78     on_MedianBlur(g_nMedianBlurValue,0);
 79     imshow("<4>中值滤波", g_dstImage4);
 80     imwrite("中值.jpg", g_dstImage4);
 81 
 82     //=================【<5>双边滤波】===========================
 83     //创建窗口
 84     namedWindow("【<5>双边滤波】", 1);
 85     //创建轨迹条
 86     createTrackbar("参数值：", "【<5>双边滤波】",&g_nBilateralFilterValue, 50,on_BilateralFilter);
 87     on_BilateralFilter(g_nBilateralFilterValue,0);
 88     imshow("<5>双边滤波", g_dstImage5);
 89     imwrite("双边.jpg", g_dstImage5);
 90 
 91     waitKey(0);
 92 
 93     return 0;
 94 }
 95 //            描述：方框滤波操作的回调函数
 96 
 97 static void on_BoxFilter(int, void *)
 98 {
 99     //方框滤波操作
100     boxFilter(g_srcImage, g_dstImage1, -1,Size( g_nBoxFilterValue+1, g_nBoxFilterValue+1));
101 }
102 
103 //            描述：均值滤波操作的回调函数
104 
105 static void on_MeanBlur(int, void *)
106 {
107     blur(g_srcImage, g_dstImage2, Size( g_nMeanBlurValue+1, g_nMeanBlurValue+1),Point(-1,-1));
108 }
109 
110 //            描述：高斯滤波操作的回调函数
111 
112 static void on_GaussianBlur(int, void *)
113 {
114     GaussianBlur(g_srcImage, g_dstImage3, Size( g_nGaussianBlurValue*2+1,g_nGaussianBlurValue*2+1 ), 0, 0);
115 }
116 
117 //            描述：中值滤波操作的回调函数
118 
119 static void on_MedianBlur(int, void *)
120 {
121     medianBlur( g_srcImage, g_dstImage4, g_nMedianBlurValue*2+1 );
122 }
123 
124 //            描述：双边滤波操作的回调函数
125 
126 static void on_BilateralFilter(int, void *)
127 {
128     bilateralFilter( g_srcImage, g_dstImage5, g_nBilateralFilterValue, g_nBilateralFilterValue*2,g_nBilateralFilterValue/2 );
129 }