基于小波变换的CT图像去噪增强研究.docx
《基于小波变换的CT图像去噪增强研究.docx》由会员分享,可在线阅读,更多相关《基于小波变换的CT图像去噪增强研究.docx(46页珍藏版)》请在冰豆网上搜索。
基于小波变换的CT图像去噪增强研究
本
科
毕
业
论
文
作者:
杨鹏专业:
电气工程指导教师:
完成日期:
2013.5.25
原创性声明
本人声明:
所呈交的论文是本人在导师指导下进行的研究成果。
除了文中特别加以标注和致谢的地方外,论文中不包含其他人已发表或撰写过的研究成果。
参与同一工作的其他同志对本研究所做的任何贡献均已在论文中作了明确的说明并表示了谢意。
签名:
日期:
本论文使用授权说明
本人完全了解南通大学有关保留、使用学位论文的规定,即:
学校有权保留论文及送交论文复印件,允许论文被查阅和借阅;学校可以公布论文的全部或部分内容。
(保密的论文在解密后应遵守此规定)
学生签名:
指导教师签名:
日期:
南通大学毕业设计(论文)立题卡
注:
1、此表一式三份,学院、教研室、学生档案各一份。
南通大学
毕业设计(论文)任务书
题目基于小波变换的CT图像增强去噪研究
学生姓名杨鹏学院电气工程学院专业电气工程及其自动化班级电094学号0912002084起讫日期2013.01.06-2013.06.06指导教师王建平职称讲师
发任务书日期2013年01月06日
注:
此表为参考表格,学院可根据专业特点,对该表格进行适当的修改。
南通大学本科生毕业设计(论文)开题报告
注:
1、学院可根据专业特点,可对该表格进行适当的修改。
南通大学电气工程学院
2013年06月
摘要
随着科技的迅速发展,医学CT图像已成为医生诊断和治疗病人的重要工具,而CT技术由于其图像分辨率较高,对人体危害小,已被作为主要手段来研究病理以及解剖等方面。
但在实际应用中,图像的传输、接收、显示等都会造成图像质量的降低,从而导致图像的可读性差,甚至于医生无法对病变部位进行有效观察并做出诊断。
因此,图像的处理就显的尤为重要。
基于小波变换的图像去噪与增强方法有很多,本文主要讨论了各种传统的小波去噪与增强方法,并对它们的图像去噪与增强效果进行了分析。
通过对传统方法的分析研究,本文提出了小波局部阈值软硬函数折中去噪和改进过的小波变换图像子带增强方法对图像进行去噪和增强,并与传统方法进行比较。
实验结果表明,本文提出的方法不仅增强了图像的细节特征以及提高图像清晰度,而且视觉效果更佳,具有一定的价值。
关键词:
医学CT图像,小波变换,图像去噪,图像增强
ABSTRACT
Withtherapiddevelopmentofscienceandtechnology,medicalimaginghasbecomeanimportanttoolfordoctorstodiagnoseandtreatpatients.Becauseofitshighimageresolution,smallharmtohumanbody,CTtechnologyhasbeenastheprimarymeanstostudythepathologyandanatomy.However,inpracticalapplications,imagetransmission,reception,display,etc.willcauseadecreaseinimagequality,resultinginpoorreadabilityoftheimage,evendoctorscannotbeeffectiveobservationofthelesiontomakeadiagnosis.Accordingly,theimageprocessingisparticularlyimportant.
Therearemanymeansbasedonwavelettransformimagedenoisingandenhancement,thispapermainlydiscussesthevarioustraditionalwaveletdenoisingandenhancementmethods,andtheirimagedenoisingandenhancementeffectswereanalyzed.Throughtheanalysisoftraditionalmethods,thispaperpresentsthelocalthresholdwaveletdenoisingcompromisehardnessandimprovedfunctionoversub-bandwavelettransformimageenhancementmethodforimagedenoisingandenhancement,andcomparedwiththetraditionalmethod.Experimentalresultsshowthattheproposedmethodnotonlyenhancesimagedetailfeaturesandimproveimageclarityandbettervisualeffects,hasacertainvalue.
Keywords:
medicalCTimage,wavelettransform,imagedenoising,imageenhancement
目录
摘要....................................................IABSTRACT...................................................II
第一章绪论.................................................1
1.1课题研究的背景及意义...............................................................................1
1.2本课题的研究内容.......................................................................................2
1.3基于小波变换的图像去噪发展与现状.......................................................2
1.4基于小波变换的图像增强发展与现状.......................................................3
1.5本文的主要工作...........................................................................................3
第2章基于小波变换的图像增强和去噪的理论基础...............4
2.1小波变换简介..............................................4
2.2经典傅里叶变换............................................4
2.3短时傅里叶变换............................................5
2.4连续小波变换..............................................5
2.5离散小波变换..............................................6
2.6二进制小波变换............................................7
2.7多分辨率分析..............................................8
第三章基于小波变换的CT图像去噪研究.........................11
3.1小波去噪原理.............................................11
3.2小波去噪方法.............................................13
3.2.1小波萎缩法..........................................14
3.2.2投影方法............................................16
3.2.3相关法..............................................16
3.3去噪图像的评价指标.......................................16
3.4基于小波变换的图像去噪算法...............................17
3.6实验结果分析.............................................20
第四章基于小波变换的医学图像增强研究........................22
4.1基于小波变换的图像增强的原理.............................23
4.2传统的直方图均衡图像增强方法.............................23
4.3基于小波变换的CT图像子图像带增强........................24
4.4本文图像增强方法的代码实现...............................25
4.5实验结果与分析...........................................25
第五章总结与展望.............................................27
5.1总结.....................................................27
5.2展望.....................................................27
附录.......................................................28
参考文献...................................................31
致谢...................................................33
第一章绪论
1.1课题研究的背景及意义
70年代初期,X-CT的发明在医学影像领域掀起了一场巨大的变革,自此之后生物医学影像得到了飞速发展,各种不同的数字成像技术被发现。
它们包括数字X线摄影(DR)、计算机X线摄影(CR)、超声成像、直接数字X线摄影(DDR0、磁共振成像技术(MR)、计算机断层成像(CT)、磁共振血管造影术(MRA)、正电子发射断层成像(PET)等。
CT图像处理技术是一种非常有用的技术,由于它可以让人类的视觉透过事物表面往内延伸的特性,该技术被广泛应用于生物医学。
但是,CT图像在采集过程中常常会受到噪声的干扰,使得实际得到的图像中含有噪声成分。
噪声的存在使得图像在结构上的相关性被破坏,从而影响了特征的提取。
例如,对人体病灶部位拍摄的CT图像会出现用肉眼无法分辨的情况,因为病灶部位在灰度以及形状上和其他部分相似。
因此,用计算机对CT图像做预处理就非常必要,否则就很难对图像做出评估。
为了改善CT图像的视觉效果,通常需要使用一些图像处理技术,其中图像去噪以及图像增强技术是最常用的两种图像处理技术。
CT图像在处理过程中主要有两大问题:
(1)去除噪声的同时要尽量保存图像的边缘和细节信息。
由于医生需要通过观察CT图像的细节特征来确定某部位是否发生病变,因此,在对CT图像的去噪的同时要尽可能的保留图像的细节特征。
常用的去噪方法很多,如领域平均法、维纳滤波法、低通滤波法等。
但这些方法都很容易破坏图像的细节特征,并且只要是对图像的锐化处理都会使得细节特征变模糊。
(2)如何增强图像以及提高图像的清晰度。
图像经过去噪处理后,不可避免的会使一些边缘和细节变模糊,这就需要对图像进行增强处理。
对于以上两个问题,本文在研究了小波变换的基础上提出了自己的图像去噪和增强方法。
小波变换最近几十年来的新兴理论,被广泛应用于各个领域,例如,信号处理、语音识别与合成、自动控制、图像处理与分析、天体物理学。
传统的图像处理方法在时域和频域上无法做到同时具有高分辨率,但小波变换却在时域和频域同时具有很好的局部化特征。
而且由于它可以聚焦到分析对象的所有细节上,它也被称为“数字显微镜”。
它的这
个特性,使得它非常实用于图像的去噪。
本课题的研究意义在于通过对CT图像的去噪和增强以提高图像对比度和清晰度,从而改善图像的质量,为医疗诊断提供帮助。
1.2本课题的研究内容
几十年来,小波变换技术一直在迅速发展,然而在图像处理方面,兼顾去噪和增强的算法研究一直是该领域的一个难点。
本课题的工作是利用小波变换的方法来实现对医学CT图像的去噪和增强。
1.3基于小波变换的图像去噪发展与现状
通常图像去噪方法主要分为两类:
一类基于频域,如低通滤波、Wiener滤波等,它们是对整幅图像的全局处理。
另一类基于空间域,如中值滤波、统计滤波等,他们对图像中某一像素中邻域的局部处理。
这两类去噪方法只能在频率域或只能在空间域展开,因此,这两类图像去噪方法都会造成图像边缘以及细节信息的损失,造成图像模糊[1]。
为此,人们就需要一种新的技术来处理图像,小波变换应运而生。
小波变换技术自问世以来被广泛应用于信号和图像处理领域,作为一种有效的时域和频域分析方法,使得它在信号和图像处理领域有很好的前景。
凭借其良好的时频局部化能力和多分辨率分析能力,小波变换可以有效地把信号和噪声区分开以及有效处理短时瞬态信号和含宽带噪声信号等。
因此,小波变换非常实用于图像去噪的同时兼顾图像细节特征。
自70年代初小波变换技术问世到今天,该技术得到了飞速发展,不同研究者基于小波变换消除噪声提出了各种不同的方法。
早期的小波去噪与有损压缩类似,首先对含噪信号或图像进行正交小波变换,然后选定一个固定的阈值与小波系数进行比较,低于该值的小波系数置零,高于该值的小波系数保留。
后来,小波收缩新号去噪方法被发现,它的去噪效果相对以前的方法更好。
后来,该方法的提出者Standford大学的Donohol和johnstone又在此基础上提出了硬阈值和软阈值准则,并对该理论做了不断完善。
但阈值法消噪会对图像产生一些不好的改变,针对这些问题,平移不变小波去噪法被发现了。
它既能有效抑制阈值去噪法产生的伪Gibbs现象,又能减小原始信号和估计信号之间的MSE和提高了原始信号和估计信号之间的SNR,从而避免了特征模糊化现象的产生。
1992年,Mallat&zhong提出了小波模极大值去噪法[2],该方法根据有用信号和噪声的小波变换在奇异点的模极大值的不同传播特性来去噪。
HuaXie等人通过马尔科夫模型来实现对图像的去噪[3]。
一直到现在,有关小波去噪的新方法不断出现在各种文献上,然而在实际应用中并不能得到很好的应用,一些问题长期得不到解决,还有待研究。
因此,该领域还有很大的探索价值,期
待越来越简单,去噪效果更好的去噪方法被发现。
1.4基于小波变换的图像增强发展与现状
图像增强技术发展的初期,经典的图像增强技术可分为三类:
第一类是空域变换增强技术,如灰度变换等;第二类是频域变换增强技术,如同态滤波等;第三类是其他增强技术。
随着时代的进步,经典的图像增强技术已不能满足现代人的要求,越来越多的新的图像增强技术被研究者们提出,以达到更好的增强效果。
这些新提出的方法大体可分为以下几类:
对经典图像增强方法的改进,如对反锐化掩膜法[4]的改进等;基于遗传算法、神经网络、模拟退火法等;基于数字形态学的图像增强技术:
还有就是基于小波变换的图像增强技术。
很多传统的图像增强技术在增强图像的同时使得噪声也被增强,使得图像的一些细节特征被破坏。
而新提出的一些图像增强技术虽然达到了图像增强的效果,但实用性差,还需要做进一步改进。
鉴于这种现状,如何增强图像的同时,不会造成噪声增强以及图像细节损失已经成为了该技术的研究方向。
目前来说,医学CT图像主要采用基于小波变换的多尺度分析法来增强图像,但这种方法容易造成边缘特征的丢失,因此基于小波变换的图像增强技术还需要探索发展。
国内外有不少的学者致力于对小波变换的图像增强研究,例如A.Laine和S.Mallat,。
A.Laine就曾提出过乳房X光片的多尺度分析法。
基于小波变换的图像增强技术已被广泛应用于很多领域,但是由于其算法的复杂性,还有很多领域并没有使用该方法,随着算法的不断改进,基于小波变换的图像增强技术必定会在未来越来越受到人们的青睐。
1.5本文的主要工作
本文在小波变换的理论基础上对CT图像的去噪和增强算法进行了研究,提出了软硬阈值折中去噪法实现对CT图像的去噪处理,并改进了子带图像增强法用作对图像的增强处理。
并通过MATLAB软件用本文的方法对人脑CT图像进行去噪与增强实验,而且与传统的方法做了对比和分析。
本文分为5章:
第一章阐述了本课题的研究背景和意义,并且分析了国内外的图像去噪和增强现状。
第二章介绍了基于小波变换的图像与增强的理论知识。
第三章和第四章分别了图像的去噪原理和增强原理以及对各种方法做了介绍,顺势提出了自己的图像去噪与增强方法,并通过实验将本文方法与传统方法进行了对比。
第五章是本人对全文的总结以及对未来的展望。
第2章基于小波变换的图像增强和去噪的理论基础
小波变换属于时频分析的一种。
传统的信号分析建立在傅里叶变换的基础上,但由于傅里叶变换在时频分析上的局限性,致使其无法处理非平稳信号。
为此,人们提出了短时傅里叶变换来克服该问题。
但短时傅里叶变换本质是对单一分辨率信号的分析方法,针对该方法的缺陷,小波变换应运而生。
小波变换又可分为连续小波变换、离散小波变换、二进制小波变换等,它们之间也各有各的优缺点。
连续小波变换在所有可能的缩放和平移上操作,而离散小波变换采用所有缩放和平移值的特定子集。
可以说,小波变换技术的出现是历史必然的选择。
传统的变换技术在处理信号方面已经无法跟上时代的脚步,因此,在国内外各位学者的探索下,小波变换诞生了。
小波变换技术不仅继承了传统变换方法的优点,而且克服了传统方法的弊端,具有非常广泛的应用前景。
2.1小波变换简介
小波变换是一种信号的时间—尺度(时间—频率)分析方法,它具有多分辨率分析的特点,而且在时频两域都具有表征信号局部特征的能力,是一种窗口大小固定不变但形状可改变,时间窗和频率窗都可以改变的时频局部化分析方法。
即在低频部分具有较高的频率分辨率和较低的时间分辨率,在高频部分具有较高的时间分辨率和较低的频率分辨率,很适合于探测正常信号中夹带的瞬态反常现象并展示其成分,所以被誉为分析信号的数字显微镜,利用连续小波变换进行动态系统检测与诊断具有良好的效果。
2.2经典傅里叶变换
傅里叶变换是众多学科领域里重要的应用工具之一。
从实用的观点看,当人们考虑傅里叶分析的时候通常是指傅里叶变换和傅里叶级数。
函数f(t∈1(LR的连续傅里叶变换定义为
((iwtFeftdt+∞
--∞=
⎰(2-1)
反变换定义为:
1((2iwtFteFwdwπ+∞-∞
=⎰(2-2)
根据公式可以发现,f(t在空间域(,-∞+∞的大小决定了(Fω在各个频点的值,而(Fω在整个频率域上的值也决定了f(t在各个时间点的大小。
由于傅里叶分析只能反应信号时域上的频谱,而不能体现信号的局部特征[5],因此傅里叶变换不具有时域局部化分析
的能力。
函数f(t的波形经过傅里叶变换后得到很多个正弦波,分析这些正弦波的标准基可得出傅里叶变换在频域具有出色的局部化的分析功能。
2.3短时傅里叶变换
为了解决标准傅里叶变换的局部分析能力只限于频域的问题。
1946年,短时傅里叶变换被提出。
短时傅里叶变换的原理是:
对信号进行划分,将其分成许多小的时间段,用傅里叶变换对每个时间段进行分析,以得到每个时间段各自的频率。
表达式表示为
*=((itR
Sftgedtωωττ--⎰(,)(2-3)
其中“*”是复共轭,g(t表示有紧支集的函数,f(t表示开始分析的信号。
在该变换中,iteω起着限频的作用,f(t的作用是作为时限。
如果时间τ变换,g(t的“时间窗”在t轴上
移动,使f(t缓缓的进行分析。
所以g(t又被称作窗口函数,而Sωτ(,)大体显示出了f(t
在时刻为τ时、频率为ω的“信号”的含有度。
而信号在[],τδτδ-+、[],ωεωε-+区域内的状态可以表示其在窗函数上的展开,该区域被称作窗口,δ和ε分别称为窗口的时宽和频宽,表示了视频分析中的分辨率,窗宽越小则分辨率就越高。
很显然,希望δ和ε都非常小,从而时频分析效果更佳,但海森堡测不准原则[6]指出δ和ε是相互制约的,不可能同时都非常小。
由此可见,虽然短时傅里叶变换具有了一些时域上的局部分析能力,但它本身却有个无法解决的问题:
矩形窗口的形状根据窗函数g(t的改变而变化。
因此,短时傅里叶变换过于单一,分辨率的改变会使得原窗函数g(t也必须做出相应改变。
所以短时傅里叶变换适用于稳定波形的分析研究,而不适合对急剧变化的波形分析研究,对急剧变化的波形的高频部分的分析需要具有良好的时间分辨率,而对急剧变化的波形的低频部分的分析则需要好的频率分辨率,而短时傅里叶变换不能满足这样的要求。
2.4连续小波变换
设2(tLRψ∈(),其傅里叶变换为(ωψ,若(ωψ满足以下的允许条件(完全重构条件或恒等分辨条件)[7]
2(RCωωωψψ=<∞⎰(2-4)
,(tψ被称为基本小波。
对母函数(tψ作伸缩变换和平移变换可得
(abtbta-⎛⎫ψ=
⎪⎝
⎭,;0abR
a∈≠(2-5)称之为小波序列。
a被称为伸缩变量,b被称为平移变量。
函数组2((ftLR∈的连续小波变换可写为
(,,((
fabR
tb
Wabfaftdta
--=ψ=ψ⎰
(2-6)
逆变换公式为
211((,(ftb
ftabdadbCaa+∞+∞
ψ
-∞-∞
-=
ψ⎰⎰(2-7)在小波变换中,由于(tψ产生的小波,(abtψ对被分析的信号起着观测窗的作用,所以
(tψ还应该遵守条件:
(tdt+∞
-∞
ψ<∞⎰
(2-8)
故函数(ωψ是连续变换的。
换句话