自主学习周 MATLAB图像添加噪声和降噪学习报告.docx

1、自主学习周 MATLAB图像添加噪声和降噪学习报告自主学习周 MATLAB图像添加噪声和降噪学习与实践报告学生姓名： 班 级： 学 号： 导 师： 时 间： 2011年9月一、背景资料MATLABMATLAB是由美国mathworks公司发布的主要面对科学计算、可视化以及交互式程序设计的高科技计算环境。它将数值分析、矩阵计算、科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的视窗环境中，为科学研究、工程设计以及必须进行有效数值计算的众多科学领域提供了一种全面的解决方案，并在很大程度上摆脱了传统非交互式程序设计语言（如C、Fortran）的编辑模式，代表了当今国际科学

2、计算软件的先进水平。 MATLAB和Mathematica、Maple并称为三大数学软件。它在数学类科技应用软件中在数值计算方面首屈一指。MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连 matlab开发工作界面接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。图像噪声目前大多数数字图像系统中，输入图像都是采用先冻结再扫描方式将多维图像变成一维电信号，再对其进行处理、存储、传输等加工变换。最后往往还要在组成多维图像信号，而图像噪声也将同样受到这样的分解和合成。在这些过程中电气系统和外界影响将使得图像噪声的

3、精确分析变得十分复杂。另一方面图像只是传输视觉信息的媒介，对图像信息的认识理解是由人的视觉系统所决定的。不同的图像噪声，人的感觉程度是不同的，这就是所谓人的噪声视觉特性课题。 图像噪声在数字图像处理技术中的重要性越来越明显，如高放大倍数航片的判读，射线图像系统中的噪声去除等已经成为不可缺少的技术步骤。3.高斯噪声所谓高斯噪声是指它的概率密度函数服从高斯分布（即正态分布）的一类噪声。如果一个噪声，它的幅度分布服从高斯分布，而它的功率谱密度又是均匀分布的，则称它为高斯白噪声。高斯白噪声的二阶矩不相关，一阶矩为常数，是指先后信号在时间上的相关性。 高斯白噪声包括热噪声和散粒噪声。编辑本段高斯分布高斯

4、分布的概念高斯分布，也称正态分布，又称常态分布，记为N（，2），其中，2为分布的参数，分别为高斯分布的期望和方差。当有确定值时，p(x)也就确定了，特别当=0，2=1时，X的分布为标准正态分布。高斯噪声的抑制在电子技术测控系统中，对干扰与噪声抑制方法主要有屏蔽、合理接地、隔离、合理布线、净化电源、滤波、采用专用器件等等措施。 除了采用通用的噪声抑制方法外，对高斯噪声的抑制方法常常采用数理统计方法。4.椒盐噪声 椒盐噪声是由图像传感器，传输信道，解码处理等产生的黑白相间的亮暗点噪声。椒盐噪声往往由图像切割引起。去除脉冲干扰级椒盐噪声最常用的算法是中值滤波。路面图像属于结构光图像，使用区域分割技术

5、中的阀值分割法消除白噪声及部分椒盐噪声，而不能使用中值滤波对白噪声及椒盐噪声进行滤波，因为滤波模板在图像中漫游时会改变光条中像素的真实灰度分布，给随后的重心法细化过程带来负面影响。 大量的实验研究发现，由摄像机拍摄得到的图像受离散的脉冲、椒盐噪声和零均值的高斯噪声的影响较严重。噪声给图像处理带来很多困难，对图像分割、特征提取、图像识别等具有直接影响。因此，实时采集的图像需进行滤波处理。消除图像中的噪声成份叫做图像的平滑化或滤波操作。滤波的目的有两个:一是抽出对象的特征作为图像识别的特征模式;二是为适应计算机处理的要求，消除图像数字化时所混入的噪声。对滤波处理的要求有两条:一是不能损坏图像轮廓及

6、边缘等重要信息;二是使图像清晰，视觉效果好。5.中值滤波中值滤波法是一种非线性平滑技术，它将每一像素点的灰度值设置为该点某邻域窗口内的所有像素点灰度值的中值. 实现方法： 1：通过从图像中的某个采样窗口取出奇数个数据进行排序 2：用排序后的中值取代要处理的数据即可 中值滤波法对消除椒盐噪音非常有效，在光学测量条纹图象的相位分析处理方法中有特殊作用，但在条纹中心分析方法中作用不大. 中值滤波在图像处理中,常用于用来保护边缘信息,是经典的平滑噪声的方法 中值滤波原理。中值滤波是基于排序统计理论的一种能有效抑制噪声的非线性信号处理技术，中值滤波的基本原理是把数字图像或数字序列中一点的值用该点的一个邻

7、域中各点值的中值代替，让周围的像素值接近的真实值，从而消除孤立的噪声点。方法是去某种结构的二维滑动模板，将板内像素按照像素值的大小进行排序，生成单调上升（或下降）的为二维数据序列。二维中值滤波输出为g（x,y）=medf(x-k,y-l),(k,lW) ，其中,f(x,y)，g(x,y)分别为原始图像和处理后图像。W为二维模板，通常为2*2，3*3区域，也可以是不同的的形状，如线状，圆形，十字形，圆环形等。 参考资料：二、程序源代码M=imread() % 打开图像“”subplot(2,3,1) imshow(M) %在一行一列显示原始图像title(original)P1=imnoise(

8、M,gaussian, %加入高斯躁声 subplot(2,3,2) imshow(P1) %在一行二列显示加入高斯躁声后的图像 title(gaussian noise);P2=imnoise(M,salt & pepper, %加入椒盐躁声subplot(2,3,3)imshow(P2) %在一行三列显示加入椒盐躁声后的图像 title(salt & pepper noise);P12=im2bw(P1, P22=im2bw(P2, %将添加噪音后的图像转为二进制 P01=im2bw(M,subplot(2,3,4)imshow(P01)title(original-binary) %将原

9、图像转为二进制并于二行一列显示P3=medfilt2(P12) %对高斯躁声中值滤波subplot(2,3,5)imshow(P3) %在二行二列显示高斯躁声中值滤波后的图像title(medfilter gaussian)P4=medfilt2(P22) %对椒盐躁声中值滤波subplot(2,3,6)imshow(P4)title(medfilter salt & pepper noise) %在二行三列显示对椒盐躁声中值滤波三、运行结果 运行环境：MATLAB R2010b四、相关函数IMNOISE Add noise to image. J = IMNOISE(I,TYPE,.) Ad

10、d noise of a given TYPE to the intensity image I. TYPE is a string that can have one of these values: gaussian Gaussian white noise with constant mean and variance localvar Zero-mean Gaussian white noise with an intensity-dependent variance poisson Poisson noise salt & pepper On and Off pixels speck

11、le Multiplicative noise Depending on TYPE, you can specify additional parameters to IMNOISE. All numerical parameters are normalized; they correspond to operations withimages with intensities ranging from 0 to 1.IM2BW Convert image to binary image by thresholding. IM2BW produces binary images from i

12、ndexed, intensity, or RGB images. To do this, it converts the input image to grayscale format (if it is not already an intensity image), and then converts this grayscale image to binary by thresholding. The output binary image BW has values of 1 (white) for all pixels in the input image with luminan

13、ce greater than LEVEL and 0 (black) for all other pixels. (Note that you specify LEVEL in the range 0,1, regardless of the class of the input image.) BW = IM2BW(I,LEVEL) converts the intensity image I to black and white. BW = IM2BW(X,MAP,LEVEL) converts the indexed image X with colormap MAP to black

14、 and white. BW = IM2BW(RGB,LEVEL) converts the RGB image RGB to black and white. Note that the function GRAYTHRESH can be used to compute LEVELautomatically.MEDFILT22-D median filtering. B = MEDFILT2(A,M N) performs median filtering of the matrix A in two dimensions. Each output pixel contains the m

15、edian value in the M-by-N neighborhood around the corresponding pixel in the input image. MEDFILT2 pads the image with zeros on the edges, so the median values for the points within M N/2 of the edges may appear distorted. B = MEDFILT2(A) performs median filtering of the matrix A using the default 3-by-3 neighborhood. B = MEDFILT2(.,PADOPT) controls how the matrix boundaries are padded. PADOPT may be zeros (the default), symmetric, or indexed. If PADOP