OpenCV学习(五)平滑滤波

原创 2013年12月03日 14:22:36

在图像处理中,会用到对图像的平滑处理。平滑常常被用来消除图像的噪声,下面介绍几个常用于图像平滑处理的函数。

(1)blur

平滑处理常常用到的是线性滤波器。线性滤波器的数学基础建立在卷积的概念之上。对于一个线性滤波器,变换之后图像的输出值其实就是该点的某个邻域内各点输入像素值的加权求和。假如我们现在有一幅M*N的图像要进行线性滤波,我们采用的卷积核(kernel)的大小为m*n,那么处理之后的图像每一点的输出像素值由下面这个式子确定:

上面这个式子的分子其实就是源图像f(x,y)与我们所用的卷积核w做卷积运算,而分母是w中各个元素的求和结果,这样就为邻域内不同位置分配了不同的权重。需要注意的是,一旦w被确定,那么分母其实是一个常数,所以我们应该提前算出分母的值,避免直接代入上式造成的重复计算和不必要的时间浪费。

而OpenCV中的blur函数就是上面所提到的线性滤波中的一种。它和其他线性滤波算法的共性点自然在于它也是一个线性滤波器,满足叠加原理。

个性的地方在于这种线性滤波所用的卷积核的各个元素都是相等的,都是1,也就是它对该点邻域内的各个其他点“一视同仁”,“不分彼此”。所以这时候,分母的大小是m*n,分子呢?自然就是邻域内各点像素值求和。而这种“大锅饭”的坏处,便是消除噪声的同时也使得源图像的某些细节变得模糊不清。具体模糊的程度要看我们所选用的卷积核的大小(也就是m,n的大小)和图像细节处的像素范围。请看下面的例子(取自冈萨雷斯数字图像处理一书),这个例子分别展示了原图和经过blur函数处理后的图像,所用的卷积核的大小依次为3*3,5*5,9*9,15*15,35*35,从字母a和噪声区域以及几个圆形点的变化可以比较容易地看出图像模糊是如何受到卷积核大小和细节本身所影响的:


	string name="image.tif";
	Mat src=imread (name,CV_LOAD_IMAGE_GRAYSCALE);
	imshow ("Original",src);
	imwrite ("Original Image.tif",src);
	Mat dst;
	blur (src,dst,Size(3,3));
	imshow ("size=3",dst);
	imwrite ("size=3.tif",dst);

(2)高斯滤波器

由于blur的“大锅饭”政策使得图像的模糊比较严重,所以人们又提出了其他的分配权重的方法,其中之一便是高斯滤波。高斯滤波的思想便是根据邻域内的点到中心点的距离来按照高斯函数分配权重。还记得那个钟形函数的曲线吗??


“假设图像是1维的,那么观察上图,不难发现中间像素的加权系数是最大的, 周边像素的加权系数随着它们远离中间像素的距离增大而逐渐减小。”

二维的Guass函数可以写成如下的形式:


值得一提的是,高斯函数作为一种很典型的函数,被广泛地应用于权重分配等一些场合,比如局部加权回归(locally wighted regression)。

同时要注意到,虽然高斯函数并不是一个线性函数,但是由于它的思想仍然是给不同区域的像素点分配不同的权重。本质上和blur并没有什么区别,它仍然是一种线性滤波。它的非线性体现在了权重的分配过程中,而这个分配过程是距离的函数,并不是像素值的函数!

同时标准差的不同,也会对滤波结果产生比较大的影响。实际过程中如何选择标准差,也是应该考虑的问题之一。可以知道,标准差越大,高斯函数的图像会越扁平,也就是说权重的下降速度会变慢,高斯滤波也就越接近于blur;标准差越小,高斯函数的图像会越高瘦,也就是认为数据集中程度较好,权重会随着距离的增大迅速衰减,这样模糊化程度会减弱,但是消除噪声的能力会下降。所以说,标准差的大小实际是对权重分配的再调整。

	string name="image.tif";
	Mat src=imread (name,CV_LOAD_IMAGE_GRAYSCALE);
	imshow ("Original",src);
	Mat dst;
	GaussianBlur (src,dst,Size(3,3),1);   //x方向的标准差选为1
	imshow ("高斯滤波",dst);
	imwrite ("高斯滤波结果.tif",dst);


(3)中值滤波

中值滤波是一种非线性滤波方法。顾名思义,这种滤波方法会将邻域内的像素值排序,选出中值赋给输出像素点。这种滤波方法模糊化较轻,尤其适用于被椒盐噪声污染的图像。

对比一下就可以发现中值滤波很好地去除了噪声。

	string name="CircuitBoard.tif";
	Mat src=imread (name,CV_LOAD_IMAGE_GRAYSCALE);
	//线性滤波
	Mat blurFilter;
	blur(src,blurFilter,Size(3,3));
	//中值滤波
	Mat medianFilter;
	medianBlur (src,medianFilter,3);

	imshow("椒盐噪声污染的图像",src);
	imshow("线性滤波结果",blurFilter);
	imshow("中值滤波结果",medianFilter);
	imwrite ("线性滤波结果.tif",blurFilter);
	imwrite ("中值滤波结果.tif",medianFilter);

(4)双边滤波

双边滤波也是一种非线性滤波方法。出于减轻图象去噪后模糊的问题,双边滤波不仅考虑权重随着距离衰减,同样也认为像素值相差越多,所占的权重应该越小。“双边滤波器也给每一个邻域像素分配一个加权系数。 这些加权系数包含两个部分, 第一部分加权方式与高斯滤波一样,第二部分的权重则取决于该邻域像素与当前像素的灰度差值。”

由于权重的分配过程中不但涉及到距离信息,更涉及到了像素值,所以容易想到这是一种非线性滤波方法。

对于这种滤波方法,这个网址给出了一份比较详细的数学论证:双边滤波的数学原理。要注意一点,双边滤波不一定要靠高斯函数来实现,但是高斯函数是一种很常用的分配权重的方法。而且我个人认为,由于“中心极限定理”,高斯函数应该是一种比较好的分配方法。

这种方法的好处便是能够在去噪和抗模糊化中取得一种比较好的折中效果。

这个函数的调用格式为(这里参考了邻域滤波:方框、高斯、中值、双边滤波这篇博客和图像平滑处理这个教程的内容):

    bilateralFilter(InputArray src, OutputArray dst, int d, double sigmaColor, double sigmaSpace,  
                          int borderType=BORDER_DEFAULT );  

d 表示滤波时像素邻域直径,d为负时由 sigaColor计算得到;d>5时不能实时处理。

sigmaColor、sigmaSpace非别表示颜色空间和坐标空间的滤波系数sigma。可以简单的赋值为相同的值。<10时几乎没有效果;>150时为油画的效果。
borderType可以不指定。

版权声明:本文为博主原创文章,未经博主允许不得转载。

相关文章推荐

OpenCV五种平滑处理(模糊处理)

cvSmooth()函数原型: CVAPI(void) cvSmooth( const CvArr* src, CvArr* dst,                       int smoot...

OpenCV学习(五)平滑滤波

在图像处理中,会用到对图像的平滑处理。平滑常常被用来消除图像的噪声,下面介绍几个常用于图像平滑处理的函数。 (1)blur 平滑处理常常用到的是线性滤波器。线性滤波器的数学基础建立在卷积的...

【OpenCV入门教程之九】 非线性滤波专场:中值滤波、双边滤波

正如我们上一篇文章中讲到的,线性滤波可以实现很多种不同的图像变换。然而非线性滤波,如中值滤波器和双边滤波器,有时可以达到更好的实现效果。邻域算子的其他一些例子还有对二值图像进行操作的形态学算子,用于计...
  • zhmxy555
  • zhmxy555
  • 2014年04月08日 15:37
  • 118372

【OpenCV学习笔记 009】图像滤波

  • 2016年10月01日 22:56
  • 2.36MB
  • 下载

openCV学习笔记(二)------图像平滑处理和边缘检测

平滑处理和边缘检测

opencv学习(7)——图像平滑

#include #include #include "opencv2/imgproc/imgproc.hpp" #include "opencv2/imgcodecs.hpp" #include ...

Opencv学习笔记(六):图像的平滑

图像的平滑 Last Edit 2013/12/26 图像的平滑操作就是对图像进行滤波。这部分主要学习如何使用4种线性滤波器。下面是Opencv2提供的4个函数。 函数一:blur:对应的掩模是一...

OpenCv学习笔记之图像平滑处理

在介绍平滑处理之前先来说说与平滑有关的滤波基础知识。   滤波是将信号中特定波段频率滤除的操作,是抑制和防止干扰的一项重要措施。滤波器是建立的一个数学模型,通过这个模型来将图像数据进行能量转化,能量低...
  • wzy8862
  • wzy8862
  • 2015年12月17日 20:15
  • 491

【学习opencv】IplImage平滑滤波 cvSmooth

显示时用到的函数及头文件: #include "stdafx.h" #include "opencv2/opencv.hpp" #include #include #include #inclu...

Opencv学习笔记(四)--图像处理平滑,锐化操作

图像平滑算法图像平滑与图像模糊是同一概念,主要用于图像的去噪。平滑要使用滤波器,为不改变图像的相位信息,一般使用线性滤波器。几种不同的平滑方法:1. 归一化滤波器 Blurs an image us...
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:OpenCV学习(五)平滑滤波
举报原因:
原因补充:

(最多只允许输入30个字)