Opencv 图片的初步处理——平滑处理

明天上午还有HP的二面  本来要睡了  想了下  还是又写了一篇

关键函数：

Smooth

各种方法的图像平滑

void cvSmooth( const CvArr* src, CvArr* dst,
               int smoothtype=CV_GAUSSIAN,
               int param1=3, int param2=0, double param3=0, double param4=0 );

src

输入图像.

dst

输出图像.

smoothtype

平滑方法:

• CV_BLUR_NO_SCALE (简单不带尺度变换的模糊) - 对每个象素的 param1×param2 领域求和。如果邻域大小是变化的，可以事先利用函数 cvIntegral 计算积分图像。
• CV_BLUR (simple blur) - 对每个象素param1×param2邻域 求和并做尺度变换 1/(param1•param2).
• CV_GAUSSIAN (gaussian blur) - 对图像进行核大小为 param1×param2 的高斯卷积
• CV_MEDIAN (median blur) - 对图像进行核大小为param1×param1 的中值滤波 (i.e. 邻域是方的).
• 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 方式处理图像.

中值滤波

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

1. 通过从图像中的某个采样窗口取出奇数个数据进行排序
2. 用排序后的中值取代要处理的数据即可

中值滤波法对消除椒盐噪音非常有效，在光学测量条纹图象的相位分析处理方法中有特殊作用，但在条纹中心分析方法中作用不大。中值滤波在图像处理中,常用于用来保护边缘信息,是经典的平滑噪声的方法

中值滤波原理

中值滤波是基于排序统计理论的一种能有效抑制噪声的非线性信号处理技术，中值滤波的基本原理是把数字图像或数字序列中一点的值用该点的一个拎域中各 点值的中值代替，让周围的像素值接近的值，从而消除孤立的噪声点。方法是去某种结构的二维滑动模板，将板内像素按照像素值的大小进行排序，生成单调上升 （或下降）的为二维数据序列。二维中值滤波输出为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平滑滤波，剔除噪声，然后求二阶导矢，用二阶导的过零点确定边缘，在计算时也是频域乘积=>空 域卷积。

滤波器就是建立的一个数学模型，通过这个模型来将图像数据进行能量转化，能量低的就排除掉，噪声就是属于低能量部分

其实编程运算的话就是一个模板运算，拿图像的八连通区域来说，中间点的像素值就等于八连通区的像素值的均值，这样达到平滑的效果

若使用理想滤波器，会在图像中产生振铃现象。采用高斯滤波器的话，系统函数是平滑的，避免了振铃现象。

+++++++++++++++++++++++++++++++++++++++++++++++++++割一割+++++++++++++++++++++++++++++++++++++++++++++++++++++++

cvCreateImage

cvCreateImage： 　　

创建头并分配数据 　　

IplImage* cvCreateImage( CvSize size, int depth, int channels ); 　　

参数说明： 　　

size 图像宽、高. 　　

depth 图像元素的位深度，可以是下面的其中之一： 　　

IPL_DEPTH_8U - 无符号8位整型 　　

IPL_DEPTH_8S - 有符号8位整型 　　

IPL_DEPTH_16U - 无符号16位整型 　　

IPL_DEPTH_16S - 有符号16位整型 　　

IPL_DEPTH_32S - 有符号32位整型 　　

IPL_DEPTH_32F - 单精度浮点数　　

IPL_DEPTH_64F - 双精度浮点数 　　

channels： 　　每个元素（像素）通道号.可以是 1, 2, 3 或 4.通道是交叉存取的，

例如通常的彩色图像数据排列是：b0 g0 r0 b1 g1 r1 ...

虽然通常 IPL 图象格式可以存贮非交叉存取的图像，并且一些OpenCV 也能处理他, 但是这个函数只能创建交叉存取图像. 　　

函数 cvCreateImage 创建头并分配数据，这个函数是下列的缩写型式： 　　

header = cvCreateImageHeader(size,depth,channels); 　　cvCreateData(header);

+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++割一割++++++++++++++++++++++++++++++++++++++++++++++++++++

代码：

#include "stdafx.h"

#include <cv.h>
#include <cxcore.h>
#include <highgui.h>


int _tmain(int argc, _TCHAR* argv[])
{
    //载入图片
    IplImage *in = cvLoadImage("funny-pictures.jpg");

    //创建显示窗口
    cvNamedWindow("In",CV_WINDOW_AUTOSIZE);
    cvNamedWindow("Out",CV_WINDOW_AUTOSIZE);

    //生成输出图片
    IplImage *out=cvCreateImage(cvGetSize(in), IPL_DEPTH_8U, 3);

    //对输出图像平滑处理
    cvSmooth( in, out, CV_GAUSSIAN, 5,5 );

    //显示图像
    cvShowImage("In",in);
    cvShowImage("Out",out);

    cvWaitKey();
    cvDestroyWindow("In");
    cvDestroyWindow("Out");
    cvReleaseImage(&in);
    cvReleaseImage(&out);

    return 0;
}

完成

结果：

注意   其中

复制图像:

IplImage* img1=cvCreateImage(cvSize(640,480),IPL_DEPTH_8U,1); 
IplImage* img2;
img2=cvCloneImage(img1);  // 注意通过cvCloneImage得到的图像
                      // 也要用 cvReleaseImage 释放，否则容易产生内存泄漏