2010年第1期中图分类号:TI田11.73文献标识码:A文章编号:1009—2552(2010)01—0065—03图像平滑处理方法的比较研究平丽(河南理工大学电气学院,焦作454000)摘要:图像平滑处理在数字图像处理中有着重要的作用。系统分析目前具有代表性的图像平滑处理方法,并用MATLAB7.0软件实现了各种算法。实验结果表明,各种方法均有各自的优缺点和适用条件,在做图像平滑处理之前,应对图像进行分析,针对图像的特点和应用需求选用合适的方法。关键词:图像增强;图像平滑J~Comparison…Og一●imagesmoottnngprocessingmethodsPINGU(ElectricalSchod,HmP‘灯懈:tI疵Urdver蛳,J'momo4544300,China)Abstract:ImagesmoothingprocessingpkIysallimportantroleindigitalimageprocessing.11lispapercomprehensivelyanalyzestherepresentativemethodsofimagesmoothingprocessingatpresent.and砌嘲eachalgorithmwiththebe删selectedMATLAB7.0software.ResultsindicatethateachmethodhasSOllleadvantagesandlimitations.Itshouldacc砌illgtothecharacteristicsoftheimageaswellasapplicationm州sbeforeconductingimage咖tlliJlgprocessing.Keywords:imageenhancement;imagesmoothiIlgO引言的含义均相同。]为含噪图像信号,则平滑处理后的图像增强的方法比较多,有加噪音,平滑,锐化,图像信号g(m,n)为:伪彩色增强等多种方法。本文重点对其中的平滑处射M理方法进行分析。图像的简单平滑是图像增强处理g(m,n)=∑默HI,一蠡,,I—f)r(一后.一z)=中最基本的方法之一,它利用卷积运算对图像邻域“m,n)*丁(一m,一厅)(1)的像素灰度进行平均化,从而达到减少图像中杂点式(1)可视为模板运算【副,其中(2M+1)×(2肘+1)影响、降低图像对比度的目的。阵列丁称为乘积模板。其算法可分解为:本文分别选取了空域平滑处理中的均值法,领①将模板r在待处理图像中平移,并将模板中域平均法,中值滤波法,巴特沃失低通滤波这几种平心与图像中某个像素位置重合。滑处理方法进行详细分析,并对不同的方法进行比②将模板r上的系数(矩阵元素)与模板下对较总结出各种方法的优缺点和适用范围,从而为人应的像素灰度值相乘。们在实际应用中选择合适的图像平滑方法提供实验③将所有乘积相加。依据。文中所有方法的实践均在MATLAB'/.0环境④将和(模板的输出响应)赋给处理后图像中中实现。对应模板中心位置的像素。1平滑处理方法⑤重复上述过程,直至模板丁的中心遍历待处1.1平滑处理的模板运算理图像中所有像素。设厂(//'t,凡)[符号.厂(m,n)不仅可以表示二元离散函数厂(m,n),而且还可以表示矩阵厂中的元素收稿日期:2009—07—07作者简介:平丽(1986一),女,硕士研究生,主要研究方向为图像处“m,n),以下g(m,n)、T(m,n)和h(m,rg)等符号理,故障诊断。一65—万 方数据要达到平滑的目的,模板r的系数首先应具有非负性,即对于所有m,n,有T(m,n)E0。因此该处理算法实际上就是进行(2M+1)×(2M+1)邻域的加权平均,各权值(模板系数)不同,平滑效果就不一样。此外,还应该满足归一性,即三(m,rt)T(m,11,)=l,这是为了使图像的平均灰度值在处理前后不发生变化,以及输出图像仍在原来的灰度范围内。尽管在空域算法中使用乘积模板的概念简洁清晰,但要进行频域分析,还须利用卷积模板的概念,为此首先引入二维离散卷积的定义。二维离散信号^(m,rt)和以(m,n)的卷积八rrt,17,)定义为八m,凡)=^(玑,n)*五(m,n)=∑∑^(k,z)五(m—k,n—z)=七=一∞l=一∞∑∑^(m—k,11,一z)A(k,z)(2)这样,式(1)可以写为Mg(m,,t)=∑∑以m—k,n—z)丁(一k,一z)=.厂(111,,,1)*r(一m,一n)(3)于是上述平滑处理系统可看作一个离散滤波器,且是线性移不变的,其冲激响应阵列h(m,n)=T(一m,一n),如图l所示。冲激响应阵列h又称为卷积模板,当乘积模板丁满足对称条件T(rrt,n)=丁(一m,一n)时,有h=T,这时h和丁可以不加区别地统称为模板。1.2邻域平均算法及其改进在式(1)中,若模板各系数均取常值:T(m,凡)=1/N,其中Ⅳ=(2M+1)2,这时输出图像信号与原图像信号之间的关系为:MMg(m,凡)=面1∑2f(rrt+k,凡+z)(4)或表示为:g(m,n)=丙1∑八k,f)(5)式中,S是以点(m,n)为中心的(2肘+1)×(2M+1)邻域。在这种方法中,对邻域中的各点平等对待(权值相同),算法简单,运算速度快,并且仍然满足三(rrt,凡)丁(m,n)=1。可以证明,平滑后的噪声标准差降低为原来的1/佃,因此邻域尺寸越大,抑制噪声的效果越显著。若要突出中心点的核心地位,以及强调临近点的作用,可将上述算法加以改进,一种改进算法的模板与邻域平均算法的模板对比如下(以3×3尺寸为.『Il11例):邻域平均h,=百1ll1了Lll和一种改进算法1l1j一66一万 方数据1211h2=素l242…LI。l21J1.3中值滤波在中值滤波算法中,对于孤立像素的属性并不非常关注,而是认为图像中的每个像素都跟邻域内其他像素有着密切的关系,对于每一个邻域,算法都会在采样得到的若干像素中,选择一个最有可能代表当前邻域特征的像素的灰度作为中心像素灰度,这样就有效避免了离散型杂点对图像的影响。明白了中值滤波算法的思路后,如何有效地选择一个最有可能代表邻域特征的像素灰度就成了算法的核心问题。直观地,邻域像素中灰度大小居中的点通常能很好地描述邻域的属性,因此在中值滤波算法中,中心像素的灰度是通过对邻域内像素灰度排序然后取中值来确定的。1.4巴特沃失低通滤波1.4.1分频对一幅图像来说,能量主要集中在低频分量中,而噪声和图像细节信息主要集中在高频分量中。在图像处理过程中,将两种分量分离开,进行不同的操作与处理,可以取得更佳的处理效果。分频的目的就是将输人图像的高频分量和低频分量分开。进行分频时,既可以选择高通滤波器又可以选择低通滤波器,都能达到分频的目的。本文实验的过程中使用的是低通滤波器。在诸多低通滤波器中,Buter-worth低通滤波器没有“振铃”现象,而且能够提高图像的细节清晰度。Buterworth低通滤波器(BLPF)传递函数:1爿【it,口J2订T丙再瓣式中,Do为截止频率,D(Ⅱ.t,)= ̄/“2+口2,在D(Ⅱ.t,)=仇处,日=告日。;为阶数,取正整数,用来控制衰减速度。对~幅输入图像“菇,Y),可将其分解为高频分量五(髫,Y)和低频分量^(菇,Y),如下式所示:以菇,Y)=厶(z,Y)+^(菇,Y)将图像.厂(戈,Y)通过2阶Buterworth低通滤波器g(戈,Y),即八髫,Y)与g(戈,Y)进行卷积,可得到图像的低频分量“菇,Y),如下:工(髫,),):/.(x,Y)*g(戈,,,)其中,g(茗,),)2r_n灭;瓣(D。为截止频率)在图像∥戈,Y)中减去低频分量^(茗,Y)就得到高频分量:五(石,Y)=八戈,Y)一^(茗,Y)1.4.2低频分量均衡化图像的能量主要集中在低频分量.^(石,Y)中,对图像的整体视觉效果起着决定性的作用,必然要2仿真结果与分析利用MATLAB对图像平滑处理算法进行仿真,结果如图1所示,其中图l(a)是原始图像,与受到噪声污染的图像,图1(b)是原图像和均值法后的图像,图l(c)是原图像和领域平均法后的图像,图l(d)是原图像和中值滤波法后的图像,图l(e)是原图像和巴特沃失低通滤波法后的图像。由仿真结果可知,简单平滑的结果中不难发现,平滑的邻域边长越大,可以消除的杂点影响就越多,同时图像的对比度下降也越多。使用简单平滑的方法进行杂点消除往往以大幅降低图像清晰度为代价。另外,使用简单平滑处理实质上只能将杂点无限放大而使人眼不易察觉,但不能真正消除杂点,这是图像简单平滑对该分量进行直方图均衡。对低频分量八菇,Y)利用传统算法进行直方图均衡,得到t,(石,Y):t,(石,Y)=艘(^(菇,Y))噪声存在于高频分量中,如果直接进行直方图均衡,必将使噪声得到增强,影响图像质量。因此,将以(石,Y)乘上一个加权系数口得到t。(菇,Y):th(戈,Y)=口*五(菇,Y)口的取值因具体情况而定:若图像受噪声影响比较严重,D的取值应小于l,以免放大噪声;若图像受噪声影响比较小,口的取值应大于Z。输出图像八茹,Y)为高频分量五(菇,Y)和低频分量^(戈,Y)之和:处理的一大缺陷。图像的中值滤波很好地解决了图像杂点的问题,然而对于细节较多的图像,中值滤波处理常常导致图像信息的丢失。巴特沃失低通滤波算法的处理结果较传统算法与原图更为接近,通过以戈,Y)=^(菇,Y)+^(石,Y)所以,输出图像以菇,Y)为五(菇,Y)和HE(f,(菇,Y))的线性组合,即:增强,更好地将退化图像恢复出原始图像的信息,并且还有了更为丰富的信息。采用Buterworth低通滤t(戈,Y)=a*五(菇,Y)+嬲(^(茗,,,))t(x,Y)就是希望得到的图像。通过分频,输人图像被分解为高频分量和低频分量,对两分量进行不同的操作,再经过线性组合得到最终的增强图像,图像清晰度高,细节信息丰富,增强效果明显。波器将图像分为高频分量和低频分量,低频分量进行HE算法直方图均衡,然后再将高频分量和低频分量进行融合得到最终处理结果,这样可以突出图像的细节信息,使图像整体视觉效果得到改善,有利于后续图像处理工作的进行。图1仿真后的结果(下转第70页)万方数据 一67—用共享MFCDLL),设工程名为Serialcomm。retum1;建立Serialcomm.DLL的工程,在定义串口GPS类CGPSSefialConrn的class后面加关键字AFX—EXT—CLASS,在CGPSSerialCormn类文件中添加:撑include”}将JSerialComm.h复制到VC工程目录中,可以将jni.h和jni—md.h复制到工程目录中,并在JSeri—JSeriMcomm.h,,,实现冠以关键字JNIEXPORT的导出函数如下:JNIEXPORTjintJNICALLJava—JSerialcomm—alcomm.h头文件中的加入#include”jni.h『,,然后编译生成SerialComm.dll,即完成了既可被Java访问也可被VC++访问的VC+DLL设计。2.3WatchSerialcomm(JN[Env*,jclass){Java调用vC++DLL将2.2节生成的Serialconma.dll文件复制到CGPSSerialcommComm;//声明串口GPS类JSerialComm.java所在的目录中,在JSerialCoram.java去掉最后两行代码,然后再编译生成JSerialComm.bStartWatch;class文件,这样在java中运行.class,可以看调用对象CommBOOLbStartWatch=Comm.StartWatch();,/if(bStartWatch==TRUE)retum1;elserelllnlSetcomm(),WatchSerialComm()开始检测串口。当然我们可以去掉JSerialcomm.java中的main方法再编译生成JSerialcomm.class,这个Java类结合动态连接库Serialconma.dll可供任何Java访问串口GPS。O:3结束语虽然Java也提供串口编程的JavaCormnunieation}JNIEXPORTvoidJNICALLJava—JSerialComm—API开发包,但利用现成的VC++串口GPS类,设计供Java访问的VC++DLL总比重新学习使用JavaCommunicationSetComm(JNIEnv*,jclass){API开发包再设计Java串口GPs类要CGPSSefialcommComm;简单的多,而且也做到了不同程序设计语言代码的可重用性。参考文献:[1]苏德富,陈志琦,霍林.基于JNl的Java与其他高级语言之间通信的研究与应用[J].广西大学学报:自然科学版,2005,30(1):6-9.[2]任俊伟,林东岱.JN!技术实现跨平台开发的研究[J].计算机应用研究,2005,22(7):180—184.if(Conch.Init(”COMl”,CBR一1200,NOPAR一1TY,7,ONESTOPBIT)!=TRuE)//初始化和设置AfxMessageBox(”串El初始化失败!”);elseAfxMessageBox(”串El初始化成功!”);return;[3]原佳丽,杨仁广.盂祥增.Java调用vc++的动态链接库[J].电脑知识与技术,2007(7):169—170.}JNIEXPORTjintJNICALLJava—JSerialcomm—JS—[4]胡乃乎,赵振.基于VC,6.0串口通信动态链接库(DLL)的设计与应用[J].微计算机信息,2007。23(3—3):214—216.[5]张永,李其申,李镍岚.基于Vc++的动态链接库的创建与调用方法[J].南昌航空工业学院学报:自然科学版。2005.19(3):81—84.topWatch(JNIEnv*,jclass){CGPSSerialCommComm;inti;[6]王敏。邵定宏,陆金桂.动态链接库技术及其应用实例[J].微计算机信息,2006,22(9~3):272—274.’t,㈣●Ⅷ●…(上接第67页)i=Comm.StopWatch();/,停止检测串口,…,…,㈣,…,…,-o●,●o●苏:常州工学院学报,2004.[8]137.责任编辑:肖滨[73李俊生.一种图像平滑处理方法的空域算法与频域分析[J].江冈萨霄斯.数字图像处理[M].北京:电子工业出版社,2004:132—参考文献:[1]朱秀昌,刘峰,胡栋.数字图像处理与图像通信[M].北京:北京邮电大学出版杜。2002:67—84.[2]陈廷标,夏良正.数字图像处理[M].北京:人民邮电出版社,1990:255—271.[9]【10jKennethRCattlemen.I)i西talimageprocessing[M].E叫酬Cliffs,NewJersey,USA:PrenticeHallinternational,Ine.,1998.TsumotoS[3]张晓黎,陆新华,吴冲龙.模板运算的传递甬数分析[J].计算机与现代化.2001(6):20—23.[4]阮秋琦.数字图像处理[M].北京:电子工业出版杜,2004.f5]高浩军,杜字人.中值滤波在圈像处理中的应用[J].电子J=程师,2004(8):35一=16.[6]夏良正.数字图像处理[M].南京:东南大学出版社,1999:43—83.[C]/l㈣ngsKotheirH.Extraction0fknowledgefromdatabasesbasedInternationalRSofIEEEC0耐mmeFuzzySys—terns,NewOrleans,Imuisiana.USA,1996:748—754.[1】]SungoJ髓,L∞YongHtM.Centerwe磷tedapplications∞imagemedla/ifiltersandenhancement[C].IEEETransactionsCircuitsandSystema,1991,38(9):984—993.责任编辑:肖滨一70一万方数据 图像平滑处理方法的比较研究
作者:作者单位:刊名:英文刊名:年,卷(期):被引用次数:
平丽, PING Li
河南理工大学电气学院,焦作,454000信息技术
INFORMATION TECHNOLOGY2010(1)1次
1.Ko Sung-Jea;Lee Yong Hoon Center weighted median filters and their applications on imageenhancement 1991(09)
2.Tsumeto S H Extraction of knowledge from databases based on RS 19963.Kenneth R Cattlemen Digital image processing 19984.冈萨雷斯 数字图像处理 2004
5.李俊生 一种图像平滑处理方法的空域算法与频域分析[期刊论文]-常州工学院学报 2004(2)6.夏良正 数字图像处理 1999
7.高浩军;杜宇人 中值滤波在图像处理中的应用[期刊论文]-电子工程师 2004(08)8.阮秋琦 数字图像处理 2004
9.张晓黎;陆新华;吴冲龙 模板运算的传递函数分析[期刊论文]-计算机与现代化 2001(06)10.陈廷标;夏良正 数字图像处理 1990
11.朱秀昌;刘峰;胡栋 数字图像处理与图像通信 2002
1.黄身桂.段念.黄辉.徐西鹏 基于图像拼接的砂轮表面三维形貌图重构[期刊论文]-金刚石与磨料磨具工程2010(6)
2.黄身桂.段念.黄辉.徐西鹏 基于图像拼接的砂轮表面三维形貌图重构[期刊论文]-金刚石与磨料磨具工程2010(6)
本文链接:http://d.g.wanfangdata.com.cn/Periodical_xxjs201001019.aspx