OpenCV之cvSmooth函数平滑滤波

1、cvSmooth函数用法

定义原型

<span style="font-size:12px;">  void cvSmooth( const CvArr* src, CvArr* dst,int smoothtype=CV_GAUSSIAN,int param1, int param2, double param3, double param4 );</span>

    src:输入图像.
    dst:输出图像.
    smoothtype平滑方法:

CV_BLUR_NO_SCALE(简单不带尺度变换的模糊)-对每个象素的 param1×param2 领域求和。如果邻域大小是变化的,可以事先利用函数 cvIntegral 计算积分图像. 但是该选项没有计算平均值的操作，所以比 CV_BLUR 的运算快，但是容易溢出，因此输入输出图像必须有不同的数值精度，才能保证不会溢出。不缩放比例的simple blur 支持8位的图像输入，但是结果必须是 IPL_DEPTH_16S(CV_16S) 或者 IPL_DEPTH_32S(CV_32S) 简单无缩放的滤波不适合 IN_PLACE 方式：因为输入与输出类型必须不同。

CV_BLUR (simple blur)- -对每个象素param1×param2邻域求和并做尺度变换 1/(param1×param2)。
CV_GAUSSIAN(gaussian blur) - -对图像进行核大小为 param1×param2 的高斯卷积。
CV_MEDIAN(median blur) - -对图像进行核大小为param1×param1 的中值滤波 (邻域是方的)。
CV_BILATERAL(双向滤波) - -应用双向 3x3 滤波，彩色 sigma=param1，空间 sigma=param2.。

    关于双向滤波，可参考http://www.dai.ed.ac.uk/CVonline/LOCAL_COPIES/MANDUCHI1/Bilateral_Filtering.html。
      param1 平滑操作的第一个参数.
      param2 平滑操作的第二个参数. 对于简单/非尺度变换的高斯模糊的情况，如果param2的值为零，则表示其被设定为param1。
      param3 对应高斯参数的 Gaussian sigma (标准差). 如果为零，则标准差由下面的核尺寸计算：
      sigma = (n/2 - 1)*0.3 + 0.8, 其中 n=param1 对应水平核,n=param2 对应垂直核 . 对小的卷积核 (3×3 to 7×7) 使用如上公式所示的标准 sigma 速度会快。如果 param3 不为零，而 param1 和 param2 为零，则核大小由sigma 计算 (以保证足够精确的操作).

         函数 cvSmooth 可使用上面任何一种方法平滑图像。每一种方法都有自己的特点以及局限。
         没有缩放的图像平滑仅支持单通道图像，并且支持8位到16位的转换(与cvSobel和cvaplace相似)和32位浮点数到32位浮点数的变换格式。
   简单模糊和高斯模糊支持 1- 或 3-通道, 8-比特和 32-比特浮点图像。这两种方法可以（in-place）方式处理图像。
   中值和双向滤波工作于 1- 或 3-通道， 8-位图像，但是不能以 in-place 方式处理图像.

中值滤波中值滤波法是一种非线性平滑技术，它将每一象素点的灰度值设置为该点某邻域窗口内的所有象素点灰度值的中值。实现方法：

通过从图像中的某个采样窗口取出奇数个数据进行排序
用排序后的中值取代要处理的数据即可
中值滤波法对消除椒盐噪音非常有效，在光学测量条纹图象的相位分析处理方法中有特殊作用，但在条纹中心分析方法中作用不大。中值滤波在图像处理中,常用于用来保护边缘信息,是经典的平滑噪声的方法

中值滤波原理

中值滤波是基于排序统计理论的一种能有效抑制噪声的非线性信号处理技术，中值滤波的基本原理是把数字图像或数字序列中一点的值用该点的一个拎域中各点值的中值代替，让周围的像素值接近的值，从而消除孤立的噪声点。方法是去某种结构的二维滑动模板，将板内像素按照像素值的大小进行排序，生成单调上升（或下降）的为二维数据序列。二维中值滤波输出为g（x,y）=med{f(x-k,y-l),(k,l∈W)} ，其中,f(x,y)，g(x,y)分别为原始图像和处理后图像。W为二维模板，通常为2*2，3*3区域，也可以是不同的的形状，如线状，圆形，十字形，圆环形等。

高斯滤波

高斯滤波实质上是一种信号的滤波器，其用途是信号的平滑处理，我们知道数字图像用于后期应用，其噪声是最大的问题，由于误差会累计传递等原因，很多图像处理教材会在很早的时候介绍Gauss滤波器，用于得到信噪比SNR较高的图像（反应真实信号）。于此相关的有Gauss-Lapplace变换，其实就是为了得到较好的图像边缘，先对图像做Gauss平滑滤波，剔除噪声，然后求二阶导矢，用二阶导的过零点确定边缘，在计算时也是频域乘积=>空域卷积。

滤波器就是建立的一个数学模型，通过这个模型来将图像数据进行能量转化，能量低的就排除掉，噪声就是属于低能量部分。其实编程运算的话就是一个模板运算，拿图像的八连通区域来说，中间点的像素值就等于八连通区的像素值的均值，这样达到平滑的效果。若使用理想滤波器，会在图像中产生振铃现象。采用高斯滤波器的话，系统函数是平滑的，避免了振铃现象。

2、实例(OpenCV 2.48+VS2012)

 
#include<cv.h>  
#include<highgui.h>  
#include<stdio.h>  
void main()
{IplImage* img_src = 0;IplImage* img_dst = 0;char* pictureImage = (char*)("11.bmp");img_src = cvLoadImage(pictureImage, 1);if (!img_src){printf("Could not load image file:%s", pictureImage);exit(0);}img_dst = cvCreateImage(cvGetSize(img_src), IPL_DEPTH_8U, 3);cvNamedWindow("原始图像:", CV_WINDOW_AUTOSIZE);cvMoveWindow("原始图像:", 100, 100);cvShowImage("原始图像:", img_src);cvNamedWindow("简单模糊:", CV_WINDOW_AUTOSIZE);cvMoveWindow("简单模糊:", 200, 100);//(simple blur)简单模糊- -对每个象素param1×param2邻域求和并做尺度变换 1/(param1×param2)。cvSmooth(img_src, img_dst, CV_BLUR, 3, 3, 0, 0);cvShowImage("简单模糊:", img_dst);//注意使用CV_BLUR_NO_SCALE时，原始图像和目标图像,他们的类型不能是一致的。IplImage* img_dst2 = 0;img_dst2 = cvCreateImage(cvGetSize(img_src), IPL_DEPTH_16S, 3);cvNamedWindow("简单不带尺度变换的模糊:", CV_WINDOW_AUTOSIZE);cvMoveWindow("简单不带尺度变换的模糊:", 300, 100);//(简单不带尺度变换的模糊)-对每个象素的 param1×param2 领域求和。如果邻域大小是变化的,可以事先利用函数 cvIntegral 计算积分图像cvSmooth(img_src, img_dst2, CV_BLUR_NO_SCALE, 3, 3, 0, 0);cvShowImage("简单不带尺度变换的模糊:", img_dst2);cvNamedWindow("高斯卷积:", CV_WINDOW_AUTOSIZE);cvMoveWindow("高斯卷积:", 400, 100);//对图像进行核大小为 param1×param2 的高斯卷积。cvSmooth(img_src, img_dst, CV_GAUSSIAN, 3, 3, 0, 0);cvShowImage("高斯卷积:", img_dst);cvNamedWindow("中值滤波:", CV_WINDOW_AUTOSIZE);cvMoveWindow("中值滤波:", 500, 100);//对图像进行核大小为param1×param1 的中值滤波 (邻域是方的)。cvSmooth(img_src, img_dst, CV_MEDIAN, 3, 3, 0, 0);cvShowImage("中值滤波:", img_dst);//  cvSmooth(pFrameMat,pFrameMat,CV_BILATERAL,3,0,0);//双向滤波  //参数不对，改成  //cvSmooth(pFrameMat,pFrameMat,CV_BILATERAL,3,3,11,11);  //试试，后面两个参数分别代表了位置上和颜色上的sigma值，不能设为0的  cvNamedWindow("双向滤波:", CV_WINDOW_AUTOSIZE);cvMoveWindow("双向滤波:", 600, 100);//(双向滤波) - -应用双向 3x3 滤波，彩色 sigma=param1，空间 sigma=param2.。cvSmooth(img_src, img_dst, CV_BILATERAL, 3, 3, 11, 11);cvShowImage("双向滤波:", img_dst);cvWaitKey(0);cvReleaseImage(&img_src);cvReleaseImage(&img_dst);cvDestroyWindow("原始图像:");cvDestroyWindow("简单不带尺度变换的模糊:");cvDestroyWindow("高斯卷积:");cvDestroyWindow("中值滤波:");cvDestroyWindow("双向滤波:");
}