自主学习周 MATLAB图像添加噪声和降噪学习报告Word文档格式.docx

自主学习周 MATLAB图像添加噪声和降噪学习报告Word文档格式.docx

《自主学习周 MATLAB图像添加噪声和降噪学习报告Word文档格式.docx》由会员分享，可在线阅读，更多相关《自主学习周 MATLAB图像添加噪声和降噪学习报告Word文档格式.docx（8页珍藏版）》请在冰豆网上搜索。

它在数学类科技应用软件中在数值计算方面首屈一指。

MATLAB可以进行矩阵运算、绘制函数和数据、实现算法、创建用户界面、连matlab开发工作界面接其他编程语言的程序等，主要应用于工程计算、控制设计、信号处理与通讯、图像处理、信号检测、金融建模设计与分析等领域。

图像噪声

目前大多数数字图像系统中，输入图像都是采用先冻结再扫描方式将多维图像变成一维电信号，再对其进行处理、存储、传输等加工变换。

最后往往还要在组成多维图像信号，而图像噪声也将同样受到这样的分解和合成。

在这些过程中电气系统和外界影响将使得图像噪声的精确分析变得十分复杂。

另一方面图像只是传输视觉信息的媒介，对图像信息的认识理解是由人的视觉系统所决定的。

不同的图像噪声，人的感觉程度是不同的，这就是所谓人的噪声视觉特性课题。

图像噪声在数字图像处理技术中的重要性越来越明显，如高放大倍数航片的判读，Ｘ射线图像系统中的噪声去除等已经成为不可缺少的技术步骤。

3.高斯噪声

　　所谓高斯噪声是指它的概率密度函数服从高斯分布（即正态分布）的一类噪声。

如果一个噪声，它的幅度分布服从高斯分布，而它的功率谱密度又是均匀分布的，则称它为高斯白噪声。

高斯白噪声的二阶矩不相关，一阶矩为常数，是指先后信号在时间上的相关性。

高斯白噪声包括热噪声和散粒噪声。

编辑本段高斯分布

高斯分布的概念

高斯分布，也称正态分布，又称常态分布，记为N（μ，σ^2），其中μ，σ^2为分布的参数，分别为高斯分布的期望和方差。

当有确定值时，p（x）也就确定了，特别当μ=0，σ^2=1时，X的分布为标准正态分布。

高斯噪声的抑制

在电子技术测控系统中，对干扰与噪声抑制方法主要有屏蔽、合理接地、隔离、合理布线、净化电源、滤波、采用专用器件等等措施。

　　除了采用通用的噪声抑制方法外，对高斯噪声的抑制方法常常采用数理统计方法。

4.椒盐噪声

椒盐噪声是由图像传感器，传输信道，解码处理等产生的黑白相间的亮暗点噪声。

椒盐噪声往往由图像切割引起。

去除脉冲干扰级椒盐噪声最常用的算法是中值滤波。

路面图像属于结构光图像，使用区域分割技术中的阀值分割法消除白噪声及部分椒盐噪声，而不能使用中值滤波对白噪声及椒盐噪声进行滤波，因为滤波模板在图像中漫游时会改变光条中像素的真实灰度分布，给随后的重心法细化过程带来负面影响。

大量的实验研究发现，由摄像机拍摄得到的图像受离散的脉冲、椒盐噪声和零均值的高斯噪声的影响较严重。

噪声给图像处理带来很多困难，对图像分割、特征提取、图像识别等具有直接影响。

因此，实时采集的图像需进行滤波处理。

消除图像中的噪声成份叫做图像的平滑化或滤波操作。

滤波的目的有两个:

一是抽出对象的特征作为图像识别的特征模式;

二是为适应计算机处理的要求，消除图像数字化时所混入的噪声。

对滤波处理的要求有两条:

一是不能损坏图像轮廓及边缘等重要信息;

二是使图像清晰，视觉效果好。

5.中值滤波　　

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

实现方法：

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区域，也可以是不同的的形状，如线状，圆形，十字形，圆环形等。

参考资料：

二、程序源代码

M=imread（'

'

）%打开图像“”

subplot（2,3,1）

imshow（M）%在一行一列显示原始图像

title（'

original'

）

P1=imnoise（M,'

gaussian'

%加入高斯躁声

subplot（2,3,2）

imshow（P1）%在一行二列显示加入高斯躁声后的图像

gaussiannoise'

）;

P2=imnoise（M,'

salt&

pepper'

%加入椒盐躁声

subplot（2,3,3）

imshow（P2）%在一行三列显示加入椒盐躁声后的图像

peppernoise'

P12=im2bw（P1,

P22=im2bw（P2,%将添加噪音后的图像转为二进制

P01=im2bw（M,

subplot（2,3,4）

imshow（P01）

original-binary'

）%将原图像转为二进制并于二行一列显示

P3=medfilt2（P12）%对高斯躁声中值滤波

subplot（2,3,5）

imshow（P3）%在二行二列显示高斯躁声中值滤波后的图像

medfiltergaussian'

P4=medfilt2（P22）%对椒盐躁声中值滤波

subplot（2,3,6）

imshow（P4）

medfiltersalt&

）%在二行三列显示对椒盐躁声中值滤波

三、运行结果

运行环境：

MATLABR2010b

四、相关函数

IMNOISE

Addnoisetoimage.

J=IMNOISE（I,TYPE,...）AddnoiseofagivenTYPEtotheintensityimage

I.TYPEisastringthatcanhaveoneofthesevalues:

'

Gaussianwhitenoisewithconstant

meanandvariance

localvar'

Zero-meanGaussianwhitenoise

withanintensity-dependentvariance

poisson'

Poissonnoise

"

OnandOff"

pixels

speckle'

Multiplicativenoise

DependingonTYPE,youcanspecifyadditionalparameterstoIMNOISE.All

numericalparametersarenormalized;

theycorrespondtooperationswith

imageswithintensitiesrangingfrom0to1.

IM2BW

Convertimagetobinaryimagebythresholding.

IM2BWproducesbinaryimagesfromindexed,intensity,orRGBimages.Todo

this,itconvertstheinputimagetograyscaleformat（ifitisnotalready

anintensityimage）,andthenconvertsthisgrayscaleimagetobinaryby

thresholding.TheoutputbinaryimageBWhasvaluesof1（white）forall

pixelsintheinputimagewithluminancegreaterthanLEVELand0（black）

forallotherpixels.（NotethatyouspecifyLEVELintherange[0,1],

regardlessoftheclassoftheinputimage.）

BW=IM2BW（I,LEVEL）convertstheintensityimageItoblackandwhite.

BW=IM2BW（X,MAP,LEVEL）convertstheindexedimageXwithcolormapMAPto

blackandwhite.

BW=IM2BW（RGB,LEVEL）convertstheRGBimageRGBtoblackandwhite.

NotethatthefunctionGRAYTHRESHcanbeusedtocomputeLEVEL

automatically.

MEDFILT2

2-Dmedianfiltering.

B=MEDFILT2（A,[MN]）performsmedianfilteringofthematrix

Aintwodimensions.Eachoutputpixelcontainsthemedian

valueintheM-by-Nneighborhoodaroundthecorresponding

pixelintheinputimage.MEDFILT2padstheimagewithzeros

ontheedges,sothemedianvaluesforthepointswithin

[MN]/2oftheedgesmayappeardistorted.

B=MEDFILT2（A）performsmedianfilteringofthematrixA

usingthedefault3-by-3neighborhood.

B=MEDFILT2（...,PADOPT）controlshowthematrixboundaries

arepadded.PADOPTmaybe'

zeros'

（thedefault）,

symmetric'

or'

indexed'

.IfPADOP