自适应中值滤波算法滤除医学图像脉冲噪声要点

课程设计说明书学生姓名：学号：学院：专业：题目：自适应中值滤波算法滤除医学图像脉冲噪声指导教师：职称:2013年月日课程设计任务书12/13学年第二学期学院：专业：学生姓名：学号：课程设计题目：自适应中值滤波算法滤除医学图像脉冲噪声起迄日期：2013年6月7日～2013年6月28日课程设计地点：指导教师：系主任：下达任务书日期:2013年6月7日课程设计任务书1．设计目的：综合运用医学图像处理的相关知识，基于理论推导，设计一种自适应中值滤波算法，利用MATLAB作为编程工具进行计算机实现，从而加深对所学知识的理解。2．设计内容和要求（包括原始数据、技术参数、条件、设计要求等）：设计一种自适应的中值滤波算法，并将其应用在医学图像的脉冲噪声抑制中，达到抑制噪声与保持边缘的良好平衡。3．设计工作任务及工作量的要求〔包括课程设计计算说明书(论文)、图纸、实物样品等〕：（1）查阅文献资料，了解相关知识；（2）确定设计方法；（3）编写MATLAB程序，并上机调试；（4）撰写课程设计说明书。课程设计任务书4．主要参考文献：1韩利竹，王华，MATALB电子仿真与应用(第2版[M])，北京：国防工业版社，2003：98～238.2RafaclC.Gonzalez,RichardE.Woods著，阮秋琦，阮宇智等译。数字图版[M]北京：电子工业出版社，2007：55～193.3赵广元著，matlab与控制系统仿真实践[M]，北京：北京航空航天大学出版社，2009：56～201.4SorinZoican,Improvedmedianfilterforimpulsenoiseremoval[J].TELSIKSSerbiaandMotenegra,Nis,2003,10(1-3):681～684.5RvanNee.OFDMWirelessMultimediaCommunications[J].RrasadR.ArtechHouse，1999，33：29861S.Haykin著.郑宝玉等译.自适应滤波器原理社[M]，北京：第四版.电子工业出版社.2003：35～356.7何振亚，自适应信号处理[M]，北京：科学出版社，2002：23～99.8邹国良，自适应滤波理论及应用[M]，石家庄：河北大学出版社，1997：55～345.5．设计成果形式及要求：提供课程设计说明书一份；MATLAB程序；6．工作计划及进度：6月7日~6月12日：查资料；6月13日~6月19日：确定设计方案；6月20日~6月25日：编写程序，上机调试；6月26日~6月27日：完成课程设计说明书；6月28日：答辩。系主任审查意见：签字：年月日Abstract：Anadaptivemedianfilterthatcanremoveimpulsenoiseinmedicalimagesisimplemented，anditiscorn—paredwithtraditionalmedianfilterandaveragingfilterbytwoobjectiveevaluationcriteria，root-mean-squareerroranddestroynessdegreeofnoisetooriginalimage．Thepixelatthecenteroffilteringwindowisdetectedtodeterminewhetheritisanimpulsenoisebasedongivencondition．Medianfilterwithadaptivewindowisappliedtoelim—inatethenoise．Thealgorithmovercomestheshortcomingofmedianfilterthatreplacesthevalueofeachpixelbythemedianofthegraylevelsintheneighborhoodofthatpixel，andreducesunnecessarylossofimagedetail．Simu—lationresultsbasedonMATLABshowthat，forlargerimpulsenoisedensity,itnotonlyeliminatesnoiseeffectively,butalsopreservestheedgesanddetailinformationoftheobjects．Thealgorithmhasbeenusedinthevirtualendo—scopesystem，andthefilteringperformanceisverysatisfactory．Keywords：virtualendoscopesystem；medicalimagefiltering；adaptivealgorithm；medianfilter；impulsenoise摘要：实现了一种滤除医学图像脉冲噪声的自适应中值滤波算法，用均方根误差和噪声对原图像的毁坏程度两个客观评价指标对该算法及传统均值、中值滤波方法进行了比较与评价。根据设定条件检测滤波窗口中心像素是否为脉冲噪声，采取滤波窗口自适应的算法来滤除脉冲噪声，去除了传统中值滤波对所有像素均用中值代替的弊端，减少了不必要的图像细节损失。基于MATLAB的仿真实验表明，对于较大密度的脉冲噪声，该算法在有效抑制噪声的同时，能较好地保护边缘和细节信息。该算法已应用于虚拟内窥镜系统中，取得了令人满意的效果。关键词：虚拟内窥镜系统；医学图像去噪；自适应算法；中值滤波；脉冲噪声目录1.引言··························································12.设计方案简介··················································23.设计条件及主要参数表··········································33.1中值滤波的去噪过程·········································33.2自适应中值滤波算法·········································33.2.1脉冲噪声特点···········································33.2.2一种自适应中值滤波·····································33.3主要参数····················································44.设计结果·······················································45.设计评述及体会·················································8参考文献·························································9附录·····························································10中北大学课程设计说明书11.引言计算机断层扫描(CT)、核磁共振(MRI)和正电子放射断层成像(PET)等高精度影像设备是虚拟内窥镜技术(VirtualEndoscope)所需的医学图像数据的来源。由于受到强电磁场的干扰或是图像传输过程中比特位发生了错误，这些影像设备所获取的医学图像数据常含有脉冲噪声，从而导致图像丢失重要的细节信息，影响了组织的分割与提取，严重地影响三维重建结果的质量和人们对图像的正确分析。因此，虚拟内窥镜系统中医学图像数据预处理模块的第一步就是图像滤波。通过滤波处理在实现噪声抑制的同时保持边缘和细节信息，增强图像的特征。中值滤波器作为非线性滤波的代表在抑制脉冲噪声方面表现出的简单有效性，使得其在实际处理过程中得到了广泛的应用。但是，传统的中值滤波对图像中的所有点均进行滤波，改变了图像中未被脉冲噪声污染的像素点的灰度值，从而影响了滤波图像的质量。因此出现了很多基于标准中值滤波器的改进算法，例如加权中值滤波（WMF）、自适应的中值滤波器（AMF）、方向中值滤波器（DWM）、开关中值滤波器（SMF）等。这些改进算法在去噪和保持图像细节方面较标准中值滤波器均有了很大程度的提高，但在实际应用中都有各自的局限性。WMF通过对滤波窗口内的像素赋不同的权值，降低了细节的损失程度，但同时其去噪性能也有所下降。DWM在滤波时，引入了图像的方向信息，采用多个不同的滤波窗口检测脉冲噪声的走向，从而进一步提高了算法的细节信息，特别是边缘信息保护能力。但其对每一个像素点均采用多个滤波窗口进行滤波，增加了算法的计算量，延长了计算时间。SMF增加了噪声检测机制以区分噪声点与信号点，在噪声密度低时能更好地保留图像细节，但其性能随着输入图像信噪比的降低接近于标准中值滤波且其计算时间随噪声密度的增大而增加。与这些方法相比较，AMF具有较为优秀的滤波性能，但随着噪声密度的增加，其保护细节的能力下降较陕。通过分析脉冲噪声特点，结合虚拟内窥镜系统实时图像处理的要求，实现了一种自适应中值滤波算法。该算法根据设定条件检测滤波窗口中心像素是否为脉冲噪声，然后自适应地采用不同大小的窗口对脉冲噪声点进行中值滤波处理，具有思想简单、速度快等特点。对较大密度的脉冲噪声，该算法在有效抑制噪声的同时，能较好地保护边缘和细节信息。中北大学课程设计说明书22.设计方案简介滤波器被广泛地用于图象的预处理,抑制图象噪声,增强对比度,以及强化图象的边沿特征.运用较为广泛的线性滤波器如平均值滤波器,能较好地抑制图象中的加性噪声.但是,线性滤波器会引起图象的钝化或模糊,使得图象中物体边界产生位移.特别是,在图象受到乘性噪声或脉冲噪声的干扰,如超声波及雷达成像中普遍存在的斑点噪声,线性滤波器就不能取得预期的效果.中值滤波器，就像其名字一样，是用该像素的相邻像素的灰度中值来代替该像素的值，是一种非线性滤波器.例如滤波窗口由3×3个象素组成,则其中5个象素的灰度值会小于等于该滤波器的输出灰度值,同时5个象素的灰度值会大于等于滤波器的输出.由此可见,对于离散的脉冲噪声,当其出现的次数小于窗口尺寸的一半时,将被抑制掉同时也能较好地保证图象的边沿特征,而且易于实现.因此它被广泛地应用于图象处理,尤其是医学图象处理,如超声波图象.但由于其使用的滤窗大小是固定不变的,当窗中噪声像素数超过有用像素之半时（噪声密度较大时）,中值滤波滤波作用大大降低。在脉冲噪声强度大于0.2时,中值滤波效果就显得不是令人满意。而自适应中值滤波器会根据一定的设定条件改变滤窗的大小,即当噪声面积较大时,通过增加滤窗的大小将噪声予以去除,同时当判断滤窗中心的像素不是噪声时,不改变其当前像素值,即不用中值代替。这样,自适应中值滤波器可以处理噪声概率更大的脉冲噪声,同时在平滑非脉冲噪声图像时能够更好地保持图像细节,这是传统中值滤波器做不到的。实验框图如下：图2.1自适应中值滤波器的流程图开始载入有噪图像检测像素点是否在设置的范围内中值滤波输出滤波后的图像结束是中北大学课程设计说明书33.设计条件及主要参数表3.1中值滤波的去噪过程给定的图像为二维信号，在信号中加入指定的椒盐噪声,然后利用中值滤波进行去噪。虽然有关中值滤波的函数是在matlab函数库中已经提供，但在图像处理中利用中值滤波去除图像中的椒盐噪声噪声却是一种有效的方法，利用中值滤波函数去除图像中的噪声过程如下：（1）使用imread()读入原始的彩色图像。（2）因为使用中值滤波器只能对灰度图像进行处理，所以利用rgb2gray()将彩色图像转化为灰度图像。（3）用imnoise()在灰度图像中加入椒盐噪声。（4）利用medfilt2()函数进行中值滤波,并在matlab环境下运行。3.2自适应中值滤波算法3.2.