Opencv 图像添加噪声

出处:http://blog.sina.com.cn/s/blog_7071fc6f0100zu0b.html

大牛帖:http://www.opencv.org.cn/forum/viewtopic.php?t=5299

[cpp]

// Xoo.cpp : Defines the entry point for the console application.   

1.   
2. #include "cv.h"   
3. #include "highgui.h"   
4. int main( int argc, char** argv )  
5. {  
6.     // 结构中载入图像：图像也是BMP图像（cvLoadImage）或者其它格式   
7.     // XML/YAML (cvLoad)   
8.     CvImage img(argc > 1 ? argv[1] : "lena.jpg", 0, CV_LOAD_IMAGE_COLOR),  
9.         img_yuv, y, noise;  
10.     //  RNG   
11.     //初始化随机数生成器状态   
12.     //CvRNG cvRNG( int64 seed=-1 );   
13.     //seed   
14.     //64-bit 的值用来初始化一个随机序列   
15.     //函数 cvRNG 初始化随机数生成器并返回其状态。   
16.     //指向这个状态的指针可以传递给函数 cvRandInt,    
17.     //cvRandReal 和 cvRandArr .    
18.     //在通常的实现中使用一个 multiply-with-carry generator 。   
19.     CvRNG rng = cvRNG(-1);  
20.     if( !img.data() ) // 检查图像是否被载入   
21.         return -1;  
22.     img_yuv = img.clone(); // 克隆图像   
23.     // 色彩空间转换   
24.     //void cvCvtColor( const CvArr* src, CvArr* dst, int code );   
25.     //src==source   
26.     //dst=destination   
27.     //int code==转换的具体格式   
28.     cvCvtColor( img, img_yuv, CV_BGR2YCrCb );   
29.     // 创建图像   
30.     //CvImage::create   
31.     //void CvImage::create(CvSize size, int depth, int channels);   
32.     //创建一个影象。   
33.     //size   
34.     //图像大小   
35.     //depth   
36.     //像素深度   
37.     //channels   
38.     //通道数   
39.     y.create( img.size(), IPL_DEPTH_8U, 1 );   
40.     noise.create( img.size(), IPL_DEPTH_32F, 1 );  
41.     //cvSplit 函数，分解图像到单个色彩通道上，然后单独处理。   
42. //void cvSplit( const CvArr* src, CvArr* dst0, CvArr* dst1,   
43. //               CvArr* dst2, CvArr* dst3 );   
44. //#define cvCvtPixToPlane cvSplit   
45. //src   
46. //原数组.   
47. //dst0...dst3   
48. //目标通道   
49. //函数 cvSplit 分割多通道数组成分离的单通道数组d。   
50. //可获得两种操作模式 . 如果原数组有N通道且前N输出数组非NULL,    
51. //所有的通道都会被从原数组中提取，如果前N个通道只有一个通道非NULL函数只提取该指定通道，   
52. //否则会产生一个错误，余下的通道（超过前N个通道的以上的）必须被设置成NULL，   
53. //对于设置了COI的IplImage 结使用cvCopy 也可以从图像中提取单通道。   
54.     cvSplit( img_yuv, y, 0, 0, 0 );   
55.     // 正态分布的随机数组   
56.     //  用随机数填充数组并更新 RNG 状态   
57.     //void cvRandArr( CvRNG* rng, CvArr* arr, int dist_type, CvScalar param1, CvScalar param2 );   
58.     //rng   
59.     //被 cvRNG 初始化的 RNG 状态.   
60.     //arr   
61.     //输出数组   
62.     //dist_type   
63.     //分布类型:   
64.     //CV_RAND_UNI - 均匀分布   
65.     //CV_RAND_NORMAL - 正态分布 或者 高斯分布   
66.     //param1   
67.     //分布的第一个参数。如果是均匀分布它是随机数范围的闭下边界。如果是正态分布它是随机数的平均值。   
68.     //param2   
69.     //分布的第二个参数。如果是均匀分布它是随机数范围的开上边界。如果是正态分布它是随机数的标准差。   
70.     //函数 cvRandArr 用均匀分布的或者正态分布的随机数填充输出数组。   
71.     //在下面的例子中该函数被用来添加一些正态分布的浮点数到二维数组的随机位置。   
72.     cvRandArr( &rng, noise, CV_RAND_NORMAL, cvScalarAll(0), cvScalarAll(20) );   
73.     //各种方法的图像平滑   
74.     //void cvSmooth( const CvArr* src, CvArr* dst,   
75.     //               int smoothtype=CV_GAUSSIAN,   
76.     //               int param1=3, int param2=0, double param3=0, double param4=0 );   
77.     //src   
78.     //输入图像.   
79.     //dst   
80.     //输出图像.   
81.     //smoothtype   
82.     //平滑方法:   
83.     //CV_BLUR_NO_SCALE (简单不带尺度变换的模糊) - 对每个象素的 param1×param2 领域求和。   
84.     //如果邻域大小是变化的，可以事先利用函数 cvIntegral 计算积分图像。   
85.     //CV_BLUR (simple blur) - 对每个象素param1×param2邻域 求和并做尺度变换 1/(param1•param2).  
86.     //CV_GAUSSIAN (gaussian blur) - 对图像进行核大小为 param1×param2 的高斯卷积   
87.     //CV_MEDIAN (median blur) - 对图像进行核大小为param1×param1 的中值滤波 (i.e. 邻域是方的).  
88.     //CV_BILATERAL (双向滤波) - 应用双向 3x3 滤波，彩色 sigma=param1，空间 sigma=param2. 关于双向滤波，   
89.     //可参考 http://www.dai.ed.ac.uk/CVonline/LOCAL_COPIES/MANDUCHI1/Bilateral_Filtering.html  
90.     //param1   
91.     //平滑操作的第一个参数.   
92.     //param2   
93.     //平滑操作的第二个参数. 对于简单/非尺度变换的高斯模糊的情况，如果param2的值 为零，则表示其被设定为param1。   
94.     //param3   
95.     //对应高斯参数的 Gaussian sigma (标准差). 如果为零，则标准差由下面的核尺寸计算：   
96.     //sigma = (n/2 - 1)*0.3 + 0.8, 其中 n=param1 对应水平核,   
97.     //                                 n=param2 对应垂直核.   
98.     //对小的卷积核 (3×3 to 7×7) 使用如上公式所示的标准 sigma 速度会快。如果 param3 不为零，  
99.     //而 param1 和 param2 为零，则核大小有 sigma 计算 (以保证足够精确的操作).   
100.     //函数 cvSmooth 可使用上面任何一种方法平滑图像。每一种方法都有自己的特点以及局限。   
101.     //没有缩放的图像平滑仅支持单通道图像，并且支持8位到16位的转换(与cvSobel和cvaplace相似)和32位浮点数到32位浮点数的变换格式。   
102.     //简单模糊和高斯模糊支持 1- 或 3-通道, 8-比特 和 32-比特 浮点图像。这两种方法可以（in-place）方式处理图像。   
103.     //中值和双向滤波工作于 1- 或 3-通道， 8-位图像，但是不能以 in-place 方式处理图像.   
104.     //   
105.     //中值滤波   
106.     //中值滤波法是一种非线性平滑技术，它将每一象素点的灰度值设置为该点某邻域窗口内的所有象素点灰度值的中值。实现方法：   
107.     //通过从图像中的某个采样窗口取出奇数个数据进行排序   
108.     //用排序后的中值取代要处理的数据即可   
109.     //中值滤波法对消除椒盐噪音非常有效，在光学测量条纹图象的相位分析处理方法中有特殊作用，但在条纹中心分析方法中作用不大。   
110.     //中值滤波在图像处理中,常用于用来保护边缘信息,是经典的平滑噪声的方法   
111.     //中值滤波原理   
112.     //中值滤波是基于排序统计理论的一种能有效抑制噪声的非线性信号处理技术，中值滤波的基本原理是把数字图像   
113.     //或数字序列中一点的值用该点的一个拎域中各点值的中值代替，让周围的像素值接近的值，从而消除孤立的噪声点。   
114.     //方法是去某种结构的二维滑动模板，将板内像素按照像素值的大小进行排序，生成单调上升（或下降）的为二维数据序列。   
115.     //二维中值滤波输出为g（x,y）=med{f(x-k,y-l),(k,l∈W)} ，其中,f(x,y)，g(x,y)分别为原始图像和处理后图像。   
116.     //W为二维模板，通常为2*2，3*3区域，也可以是不同的的形状，如线状，圆形，十字形，圆环形等。  
117.     //高斯滤波   
118.     //高斯滤波实质上是一种信号的滤波器，其用途是信号的平滑处理，我们知道数字图像用于后期应用，其噪声是最大的问题，   
119.     //由于误差会累计传递等原因，很多图像处理教材会在很早的时候介绍Gauss滤波器，用于得到信噪比SNR较高的图像（反应真实信号）。   
120.     //于此相关的有Gauss-Lapplace变换，其实就是为了得到较好的图像边缘，先对图像做Gauss平滑滤波，   
121.     //剔除噪声，然后求二阶导矢，用二阶导的过零点确定边缘，在计算时也是频域乘积=>空域卷积。   
122.     //滤波器就是建立的一个数学模型，通过这个模型来将图像数据进行能量转化，能量低的就排除掉，噪声就是属于低能量部分   
123.     //其实编程运算的话就是一个模板运算，拿图像的八连通区域来说，中间点的像素值就等于八连通区的像素值的均值，这样达到平滑的效果   
124.     //若使用理想滤波器，会在图像中产生振铃现象。采用高斯滤波器的话，系统函数是平滑的，避免了振铃现象。   
125.     cvSmooth( noise, noise, CV_GAUSSIAN, 5, 5, 1, 1 ); // GAUSSIAN滤波做平衡   
126.     //将帧叠加到累积器（accumulator）中   
127.     //void cvAcc( const CvArr* image, CvArr* sum, const CvArr* mask=NULL );   
128.     //image   
129.     //输入图像, 1- 或 3-通道, 8-比特或32-比特浮点数. (多通道的每一个通道都单独处理).   
130.     //sum   
131.     //同一个输入图像通道的累积，32-比特或64-比特浮点数数组.   
132.     //mask   
133.     //可选的运算 mask.   
134.     // noise = noise + y   
135.     cvAcc( y, noise );    
136.     //ConvertScale   
137.     //使用线性变换转换数组   
138.     //void cvConvertScale( const CvArr* src, CvArr* dst, double scale=1, double shift=0 );   
139.     //#define cvCvtScale cvConvertScale   
140.     //#define cvScale  cvConvertScale   
141.     //#define cvConvert( src, dst )  cvConvertScale( (src), (dst), 1, 0 )   
142.     //src   
143.     //输入数组.   
144.     //dst   
145.     //输出数组   
146.     //scale   
147.     //比例因子.   
148.     //shift   
149.     //该加数被加到输入数组元素按比例缩放后得到的元素上   
150.     //函数 cvConvertScale 有多个不同的目的因此就有多个同义函数（如上面的#define所示）。    
151.     //该函数首先对输入数组的元素进行比例缩放，然后将shift加到比例缩放后得到的各元素上，   
152.     //即： dst(I)=src(I)*scale + (shift,shift,...)，最后可选的类型转换将结果拷贝到输出数组。   
153.     //多通道的数组对各个通道是独立处理的。   
154.     //类型转换主要用舍入和溢出截断来完成。   
155.     //也就是如果缩放+转换后的结果值不能用输出数组元素类型值精确表达，   
156.     //就设置成在输出数组数据轴上最接近该数的值。   
157.     //如果 scale=1, shift=0 就不会进行比例缩放. 这是一个特殊的优化，   
158.     //相当于该函数的同义函数名：cvConvert 。如果原来数组和输出数组的类型相同，   
159.     //这是另一种特殊情形，可以被用于比例缩放和平移矩阵或图像，   
160.     //此时相当于该函数的同义函数名：cvScale。   
161.     cvConvert( noise, y ); // y = noise * 1 + 0   
162.     //Merge   
163.     //从几个单通道数组组合成多通道数组或插入一个单通道数组   
164.     //void cvMerge( const CvArr* src0, const CvArr* src1,   
165.     //              const CvArr* src2, const CvArr* src3, CvArr* dst );   
166.     //#define cvCvtPlaneToPix cvMerge   
167.     //src0... src3   
168.     //输入的通道.   
169.     //dst   
170.     //输出数组.   
171.     //函数cvMerge 是前一个函数的反向操作。   
172.     //如果输出数组有N个通道并且前N个输入通道非NULL，   
173.     //就拷贝所有通道到输出数组，如果在前N个通道中只有一个单通道非NULL ，   
174.     //只拷贝这个通道到输出数组，否则 就会产生错误。   
175.     //除前N通道以外的余下的通道必须置NULL。   
176.     //对于设置了COI的 IplImage结构使用 cvCopy也可以实现向图像中插入一个通道 。   
177.     cvMerge( y, 0, 0, 0, img_yuv );  // 图层合并   
178.     cvCvtColor( img_yuv, img, CV_YCrCb2BGR ); // 图像色彩空间转换   
179.     //cvNamedWindow   
180.     //   
181.     //创建窗口   
182.     //int cvNamedWindow( const char* name, int flags=CV_WINDOW_AUTOSIZE );   
183.     //name    
184.     //窗口的名字，它被用来区分不同的窗口，并被显示为窗口标题。   
185.     //flags    
186.     //窗口属性标志，为1时表示会根据图像自动调整窗口大小。   
187.     //目前唯一支持的标志是CV_WINDOW_AUTOSIZE。   
188.     //当这个标志被设置后，用户不能手动改变窗口大小，   
189.     //窗口大小会自动调整以适合被显示图像（参考cvShowImage）。   
190.     //函数cvNamedWindow创建一个可以放置图像和trackbar的窗口。   
191.     //被创建的窗口可以通过它们的名字被引用。   
192.     //如果已经存在这个名字的窗口，这个函数将不做任何事情。   
193.     cvNamedWindow( "image with grain", CV_WINDOW_AUTOSIZE );  
194.     //cvShowImage   
195.     //   
196.     //在指定窗口中显示图像   
197.     //void cvShowImage( const char* name, const CvArr* image );   
198.     //name   
199.     //窗口的名字。   
200.     //image   
201.     //被显示的图像。   
202.     //函数cvShowImage 在指定窗口中显示图像。   
203.     //如果窗口创建的时候被设定标志CV_WINDOW_AUTOSIZE，   
204.     //那么图像将以原始尺寸显示；否则，图像将被伸缩以适合窗口大小。   
205.     img.show( "image with grain" ); // cvShowImage的另外一种形式   
206.     //cvWaitKey   
207.     //   
208.     //等待按键事件   
209.     //int cvWaitKey( int delay=0 );   
210.     //delay    
211.     //延迟的毫秒数。   
212.     //函数cvWaitKey无限制的等待按键事件（delay<=0时）；或者延迟"delay"毫秒。   
213.     //返回值为被按键的值，如果超过指定时间则返回-1。   
214.     //注释：这个函数是HighGUI中唯一能够获取和操作事件的函数，   
215.     //所以在一般的事件处理中，它需要周期地被调用，   
216.     //除非HighGUI被用在某些能够处理事件的环境中。   
217.       
218.     cvWaitKey();  
219.     return 0;  
220.     // 所有图像自动释放，这是使用C++类比较方便的地方   
221. }