12、【opencv入门】线性滤波:方框滤波、均值滤波与高斯滤波

一、理论与概念讲解

1、关于平滑处理

  “平滑处理“(smoothing)也称“模糊处理”(bluring),是一项简单且使用频率很高的图像处理方法。平滑处理的用途有很多,最常见的是用来减少图像上的噪点或者失真。在涉及到降低图像分辨率时,平滑处理是非常好用的方法。

2、图像滤波与滤波器

  图像滤波,即在尽量保留图像细节特征的条件下对目标图像的噪声进行抑制,是图像预处理中不可缺少的操作,其处理效果的好坏将直接影响到后续图像处理和分析的有效性和可靠性。
  消除图像中的噪声成分叫作图像的平滑化或滤波操作。信号或图像的能量大部分集中在幅度谱的低频和中频段是很常见的,而在较高频段,感兴趣的信息经常被噪声淹没。因此一个能降低高频成分幅度的滤波器就能够减弱噪声的影响。

图像滤波的目的有两个:

  一是抽出对象的特征作为图像识别的特征模式;

  一个是为适应图像处理的要求,消除图像数字化时所混入的噪声。
而对滤波处理的要求也有两条:

  一是不能损坏图像的轮廓及边缘等重要信息;

  二是使图像清晰视觉效果好。
平滑滤波是低频增强的空间域滤波技术。它的目的有两类:

  一类是模糊;

  一类是消除噪音。
  空间域的平滑滤波一般采用简单平均法进行,就是求邻近像元点的平均亮度值。邻域的大小与平滑的效果直接相关,邻域越大平滑的效果越好,但邻域过大,平滑会使边缘信息损失的越大,从而使输出的图像变得模糊,因此需合理选择邻域的大小。
  关于滤波器,一种形象的比喻法是:我们可以把滤波器想象成一个包含加权系数的窗口,当使用这个滤波器平滑处理图像时,就把这个窗口放到图像之上,透过这个窗口来看我们得到的图像。
  滤波器的种类有很多, 在新版本的OpenCV中,提供了如下五种常用的图像平滑处理操作方法,且他们分别被封装在单独的函数中,使用起来非常方便:

1     方框滤波——boxblur函数
2     均值滤波(邻域平均滤波)——blur函数
3     高斯滤波——GaussianBlur函数
4     中值滤波——medianBlur函数
5     双边滤波——bilateralFilter函数

3、线性滤波器的简介

线性滤波器:线性滤波器经常用于剔除输入信号中不想要的频率或者从许多频率中选择一个想要的频率。

几种常见的线性滤波器:

1     允许低频率通过的低通滤波器。
2     允许高频率通过的高通滤波器。
3     允许一定范围频率通过的带通滤波器。
4     阻止一定范围频率通过并且允许其它频率通过的带阻滤波器。
5     允许所有频率通过、仅仅改变相位关系的全通滤波器。
6     阻止一个狭窄频率范围通过的特殊带阻滤波器,陷波滤波器(Band-stop filter)。

4、关于滤波和模糊

关于滤波和模糊,

  滤波是将信号中特定波段频率滤除的操作,是抑制和防止干扰的一项重要措施。

  滤波可分低通滤波和高通滤波两种。而高斯滤波是指用高斯函数作为滤波函数的滤波操作,至于是不是模糊,要看是高斯低通还是高斯高通,低通就是模糊,高通就是锐化。

5、邻域算子与线性邻域滤波

  邻域算子(局部算子)是利用给定像素周围的像素值的决定此像素的最终输出值的一种算子。而线性邻域滤波是一种常用的邻域算子,像素的输出值取决于输入像素的加权和,具体过程如下图。

图注:邻域滤波(卷积):左边图像与中间图像的卷积产生右边图像。目标图像中蓝色标记的像素是利用原图像中红色标记的像素计算得到的。

  邻域算子除了用于局部色调调整以外,还可以用于图像滤波,实现图像的平滑和锐化,图像边缘增强或者图像噪声的去除。本篇文章,我们介绍的主角是线性邻域滤波算子,即用不同的权重去结合一个小邻域内的像素,来得到应有的处理效果。

线性滤波处理的输出像素值是输入像素值的加权和 :

其中的加权和我们称其为“核”,滤波器的加权系数,即滤波器的“滤波系数”。

上面的式子可以简单写作:

 

其中f表示输入像素值,h表示加权系数“核“,g表示输出像素值

在新版本的OpenCV中,提供了如下三种常用的线性滤波操作,他们分别被封装在单独的函数中,使用起来非常方便: 

1     方框滤波——boxblur函数
2     均值滤波——blur函数
3     高斯滤波——GaussianBlur函数

6、方框滤波(box Filter)

  方框滤波(box Filter)被封装在一个名为boxblur的函数中,即boxblur函数的作用是使用方框滤波器(box filter)来模糊一张图片,从src输入,从dst输出。函数原型如下:

1 C++: void boxFilter(InputArray src,OutputArray dst, int ddepth, Size ksize, Point anchor=Point(-1,-1), boolnormalize=true, int borderType=BORDER_DEFAULT )

参数详解:

1     第一个参数,InputArray类型的src,输入图像,即源图像,填Mat类的对象即可。该函数对通道是独立处理的,且可以处理任意通道数的图片,但需要注意,待处理的图片深度应该为CV_8U, CV_16U, CV_16S, CV_32F 以及 CV_64F之一。
2     第二个参数,OutputArray类型的dst,即目标图像,需要和源图片有一样的尺寸和类型。
3     第三个参数,int类型的ddepth,输出图像的深度,-1代表使用原图深度,即src.depth()。
4     第四个参数,Size类型(对Size类型稍后有讲解)的ksize,内核的大小。一般这样写Size( w,h )来表示内核的大小( 其中,w 为像素宽度, h为像素高度)。Size(3,3)就表示3x3的核大小,Size(5,5)就表示5x5的核大小
5     第五个参数,Point类型的anchor,表示锚点(即被平滑的那个点),注意他有默认值Point(-1,-1)。如果这个点坐标是负值的话,就表示取核的中心为锚点,所以默认值Point(-1,-1)表示这个锚点在核的中心。
6     第六个参数,bool类型的normalize,默认值为true,一个标识符,表示内核是否被其区域归一化(normalized)了。
7     第七个参数,int类型的borderType,用于推断图像外部像素的某种边界模式。有默认值BORDER_DEFAULT,我们一般不去管它。

boxFilter()函数方框滤波所用的核为:

  其中f表示原图,h表示核,g表示目标图,当normalize=true的时候,方框滤波就变成了我们熟悉的均值滤波。也就是说,均值滤波是方框滤波归一化(normalized)后的特殊情况。其中,归一化就是把要处理的量都缩放到一个范围内,比如(0,1),以便统一处理和直观量化。

  而非归一化(Unnormalized)的方框滤波用于计算每个像素邻域内的积分特性,比如密集光流算法(dense optical flow algorithms)中用到的图像倒数的协方差矩阵(covariance matrices of image derivatives)

  如果我们要在可变的窗口中计算像素总和,可以使用integral()函数。

7、均值滤波

  均值滤波,是最简单的一种滤波操作,输出图像的每一个像素是核窗口内输入图像对应像素的像素的平均值( 所有像素加权系数相等),其实说白了它就是归一化后的方框滤波。通过源码剖析会发现,blur函数内部中其实就是调用了一下boxFilter。

(1)均值滤波的理论简析

  均值滤波是典型的线性滤波算法,主要方法为邻域平均法,即用一片图像区域的各个像素的均值来代替原图像中的各个像素值。一般需要在图像上对目标像素给出一个模板(内核),该模板包括了其周围的临近像素(比如以目标像素为中心的周围8(3x3-1)个像素,构成一个滤波模板,即去掉目标像素本身)。再用模板中的全体像素的平均值来代替原来像素值。即对待处理的当前像素点(x,y),选择一个模板,该模板由其近邻的若干像素组成,求模板中所有像素的均值,再把该均值赋予当前像素点(x,y),作为处理后图像在该点上的灰度个g(x,y),即个g(x,y)=1/m ∑f(x,y) ,其中m为该模板中包含当前像素在内的像素总个数。

 (2)均值滤波的缺陷

  均值滤波本身存在着固有的缺陷,即它不能很好地保护图像细节,在图像去噪的同时也破坏了图像的细节部分,从而使图像变得模糊,不能很好地去除噪声点。

 (3)在OpenCV中使用均值滤波——blur函数

blur函数的原型:

1 C++: void blur(InputArray src, OutputArraydst, Size ksize, Point anchor=Point(-1,-1), int borderType=BORDER_DEFAULT )
1     第一个参数,InputArray类型的src,输入图像,即源图像,填Mat类的对象即可。该函数对通道是独立处理的,且可以处理任意通道数的图片,但需要注意,待处理的图片深度应该为CV_8U, CV_16U, CV_16S, CV_32F 以及 CV_64F之一。
2     第二个参数,OutputArray类型的dst,即目标图像,需要和源图片有一样的尺寸和类型。比如可以用Mat::Clone,以源图片为模板,来初始化得到如假包换的目标图。
3     第三个参数,Size类型(对Size类型稍后有讲解)的ksize,内核的大小。一般这样写Size( w,h )来表示内核的大小( 其中,w 为像素宽度, h为像素高度)。Size(3,3)就表示3x3的核大小,Size(5,5)就表示5x5的核大小
4     第四个参数,Point类型的anchor,表示锚点(即被平滑的那个点),注意他有默认值Point(-1,-1)。如果这个点坐标是负值的话,就表示取核的中心为锚点,所以默认值Point(-1,-1)表示这个锚点在核的中心。
5     第五个参数,int类型的borderType,用于推断图像外部像素的某种边界模式。有默认值BORDER_DEFAULT,我们一般不去管它。

8、高斯滤波

 (1)高斯滤波的理论简析

  高斯滤波是一种线性平滑滤波,适用于消除高斯噪声,广泛应用于图像处理的减噪过程。通俗的讲,高斯滤波就是对整幅图像进行加权平均的过程,每一个像素点的值,都由其本身和邻域内的其他像素值经过加权平均后得到。高斯滤波的具体操作是:用一个模板(或称卷积、掩模)扫描图像中的每一个像素,用模板确定的邻域内像素的加权平均灰度值去替代模板中心像素点的值。

  大家常常说高斯滤波最有用的滤波操作,虽然它用起来,效率往往不是最高的。

  高斯模糊技术生成的图像,其视觉效果就像是经过一个半透明屏幕在观察图像,这与镜头焦外成像效果散景以及普通照明阴影中的效果都明显不同。高斯平滑也用于计算机视觉算法中的预先处理阶段,以增强图像在不同比例大小下的图像效果(参见尺度空间表示以及尺度空间实现)。从数学的角度来看,图像的高斯模糊过程就是图像与正态分布做卷积。由于正态分布又叫作高斯分布,所以这项技术就叫作高斯模糊。

  图像与圆形方框模糊做卷积将会生成更加精确的焦外成像效果。由于高斯函数的傅立叶变换是另外一个高斯函数,所以高斯模糊对于图像来说就是一个低通滤波操作。

   高斯滤波器是一类根据高斯函数的形状来选择权值的线性平滑滤波器。高斯平滑滤波器对于抑制服从正态分布的噪声非常有效。一维零均值高斯函数为:

  其中,高斯分布参数Sigma决定了高斯函数的宽度。对于图像处理来说,常用二维零均值离散高斯函数作平滑滤波器。二维高斯函数为:

 

(2)在OpenCV中使用高斯滤波——GaussianBlur函数

   GaussianBlur函数的作用是用高斯滤波器来模糊一张图片,对输入的图像src进行高斯滤波后用dst输出。它将源图像和指定的高斯核函数做卷积运算,并且支持就地过滤(In-placefiltering)。

1  C++: void GaussianBlur(InputArray src,OutputArray dst, Size ksize, double sigmaX, double sigmaY=0, intborderType=BORDER_DEFAULT )
1     第一个参数,InputArray类型的src,输入图像,即源图像,填Mat类的对象即可。它可以是单独的任意通道数的图片,但需要注意,图片深度应该为CV_8U,CV_16U, CV_16S, CV_32F 以及 CV_64F之一。
2     第二个参数,OutputArray类型的dst,即目标图像,需要和源图片有一样的尺寸和类型。比如可以用Mat::Clone,以源图片为模板,来初始化得到如假包换的目标图。
3     第三个参数,Size类型的ksize高斯内核的大小。其中ksize.width和ksize.height可以不同,但他们都必须为正数和奇数。或者,它们可以是零的,它们都是由sigma计算而来。
4     第四个参数,double类型的sigmaX,表示高斯核函数在X方向的的标准偏差。
5     第五个参数,double类型的sigmaY,表示高斯核函数在Y方向的的标准偏差。若sigmaY为零,就将它设为sigmaX,如果sigmaX和sigmaY都是0,那么就由ksize.width和ksize.height计算出来。
6     为了结果的正确性着想,最好是把第三个参数Size,第四个参数sigmaX和第五个参数sigmaY全部指定到。
7     第六个参数,int类型的borderType,用于推断图像外部像素的某种边界模式。有默认值BORDER_DEFAULT,我们一般不去管它。

二、示例

 1、方框滤波

 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 
10 int main( )
11 {
12        //载入原图
13        Mat image=imread("1.jpg");
14 
15        //创建窗口
16        namedWindow("方框滤波【原图】" );
17        namedWindow("方框滤波【效果图】");
18        namedWindow("方框滤波【效果图1】");
19 
20        //显示原图
21        imshow("方框滤波【原图】", image );
22 
23        //进行滤波操作
24        Mat out;
25        boxFilter(image, out, -1,Size(3, 3));
26 
27        Mat out1;
28        boxFilter(image, out1, -1, Size(5, 5));
29 
30        //显示效果图
31        imshow("方框滤波【效果图】" ,out );
32 
33        imshow("方框滤波【效果图1】", out1);
34 
35        waitKey(0 );
36 }

 效果显示:

                              原图                                                    效果图                                                效果图1

       

 (2)均值滤波

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

效果展示:

                              原图                                                 效果图                                                   效果图1

       

(3)高斯滤波

 1 #include "opencv2/core/core.hpp"
 2 #include"opencv2/highgui/highgui.hpp"
 3 #include"opencv2/imgproc/imgproc.hpp"
 4 
 5 
 6 using namespace cv;
 7 
 8 //控制台应用程序的入口函数,我们的程序从这里开始
 9 int main( )
10 {
11     //载入原图
12     Mat image=imread("1.jpg");
13 
14     //创建窗口
15     namedWindow("高斯滤波【原图】" );
16     namedWindow("高斯滤波【效果图】");
17     namedWindow("高斯滤波【效果图1】");
18     //显示原图
19     imshow("高斯滤波【原图】", image );
20 
21     //进行均值滤波操作
22     Mat out;
23     GaussianBlur(image, out, Size( 3, 3 ), 0, 0 );
24 
25     Mat out1;
26     GaussianBlur(image, out1, Size( 5, 5 ), 0, 0 );
27 
28     //显示效果图
29     imshow("高斯滤波【效果图】" ,out );
30     imwrite("效果图.jpg", out);
31 
32     imshow("高斯滤波【效果图1】", out1);
33     imwrite("效果图1.jpg", out1);
34 
35     waitKey(0 );
36 }

效果展示:

                              原图                                                  效果图                                                       效果图1

       

【综合示例】

 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;//存储图片的Mat类型
10 int g_nBoxFilterValue=3;  //方框滤波参数值
11 int g_nMeanBlurValue=3;  //均值滤波参数值
12 int g_nGaussianBlurValue=3;  //高斯滤波参数值
13 
14 
15 //轨迹条的回调函数
16 static void on_BoxFilter(int, void *);            //方框滤波
17 static void on_MeanBlur(int, void *);           //均值滤波
18 static void on_GaussianBlur(int, void *);                    //高斯滤波
19 
20 int main(  )
21 {
22     //载入原图
23     g_srcImage= imread( "1.jpg", 1 );
24     if(!g_srcImage.data )
25     {
26         cout << "载入错误!" << endl;
27         return false;
28     }
29 
30     //克隆原图到三个Mat类型中
31     g_dstImage1= g_srcImage.clone( );
32     g_dstImage2= g_srcImage.clone( );
33     g_dstImage3= g_srcImage.clone( );
34 
35     //显示原图
36     namedWindow("【<0>原图窗口】", 1);
37     imshow("【<0>原图窗口】",g_srcImage);
38 
39 
40     //方框滤波】
41     //创建窗口
42     namedWindow("【<1>方框滤波】", 1);
43     //创建轨迹条
44     createTrackbar("内核值:", "【<1>方框滤波】",&g_nBoxFilterValue, 40,on_BoxFilter );
45     on_MeanBlur(g_nBoxFilterValue,0);
46     //显示窗口
47     imshow("【<1>方框滤波】", g_dstImage1);
48     imwrite("方框滤波.jpg", g_dstImage1);
49 
50 
51     //均值滤波
52     //创建窗口
53     namedWindow("【<2>均值滤波】", 1);
54     //创建轨迹条
55     createTrackbar("内核值:", "【<2>均值滤波】",&g_nMeanBlurValue, 40,on_MeanBlur );
56     on_MeanBlur(g_nMeanBlurValue,0);
57     //显示窗口
58     imshow("【<2>均值滤波】", g_dstImage2);
59     imwrite("均值滤波.jpg", g_dstImage2);
60 
61     //高斯滤波
62     //创建窗口
63     namedWindow("【<3>高斯滤波】", 1);
64     //创建轨迹条
65     createTrackbar("内核值:", "【<3>高斯滤波】",&g_nGaussianBlurValue, 40,on_GaussianBlur );
66     on_GaussianBlur(g_nGaussianBlurValue,0);
67     //显示窗口
68     imshow("【<3>高斯滤波】", g_dstImage3);
69     imwrite("高斯滤波.jpg", g_dstImage3);
70     //等待按键输入
71     waitKey(0);
72 
73     return 0;
74 }
75 
76 static void on_BoxFilter(int, void *)
77 {
78     //方框滤波操作
79     boxFilter(g_srcImage, g_dstImage1, -1,Size( g_nBoxFilterValue+1, g_nBoxFilterValue+1));
80 }
81 
82 static void on_MeanBlur(int, void *)
83 {
84     //均值滤波操作
85     blur(g_srcImage, g_dstImage2, Size( g_nMeanBlurValue+1, g_nMeanBlurValue+1),Point(-1,-1));
86 }
87 
88 static void on_GaussianBlur(int, void *)
89 {
90     //高斯滤波操作
91     GaussianBlur(g_srcImage, g_dstImage3, Size( g_nGaussianBlurValue*2+1,g_nGaussianBlurValue*2+1 ), 0, 0);
92 }

效果展示:

                       方框滤波                                                 均值滤波                                                  高斯滤波

       

原文地址:https://www.cnblogs.com/Long-w/p/9663190.html