毕业设计论文matlab图像分割.docx
- 文档编号:8101786
- 上传时间:2023-01-28
- 格式:DOCX
- 页数:18
- 大小:28.77KB
毕业设计论文matlab图像分割.docx
《毕业设计论文matlab图像分割.docx》由会员分享,可在线阅读,更多相关《毕业设计论文matlab图像分割.docx(18页珍藏版)》请在冰豆网上搜索。
毕业设计论文matlab图像分割
毕业设计(论文)matlab图像分割
毕业设计(论文)matlab图像分割
数字图像的多分辨率分析处理方法研究
―基于小波变换的医学图像分割的研究
电信学院电子信息工程专业
摘要
图像分割是一种重要的图像分析技术。
对图像分割的研究一直是图像技术研究中的热点和焦点。
医学图像分割是图像分割的一个重要应用领域,也是一个经典难题,至今已有上千种分割方法,既有经典的方法也有结合新兴理论的方法。
本论文首先介绍了双峰法以及最大类方差自动阈值法,然后重点介绍一种基于小波变换的图像分割方法,该方法先对图像的灰度直方图进行小波多尺度变换,然后从较大的尺度系数到较小的尺度系数逐步定位出灰度阈值。
最后,对这几种算法的分割效果进行了比较。
实验结果表明,本设计能够实时稳定的对目标分割提取,分割效果良好。
医学图像分割是医学图像处理中的一个经典难题。
图像分割能够自动或半自动描绘出医学图像中的解剖结构和其它感兴趣的区域,从而有助于医学诊断。
关键词:
小波变换;图像分割;阈值
Abstract
Theimagesegmentationisanimportanttechnologyofimageprocessing.Itisstillahotpointandfocusofimageprocessing.Medicalimagesegmentationisanimportantapplicationinthefieldofimagesegmentation,anditisalsoaclassicaldifficultproblemforresearchers.Thousandsofmethodshavebeenputforwardtomedicalimagesegmentation.Someuseclassicalmethodsandothersusenewmethods.
Inthispaper,firstintroducedthepetronasmethodandimumbetweenclassvariance.Thenfocusintroducedamethodofimagesegmentationbasedonwavelettransformisdiscussed.Inthismethod,thewaveletmultiscaletransformofimagegrayhistogramisdonefirst.Moreover,thegraythresholdisgraduallyfoundoutfromlargescalecoefficientstosmallscalecoefficients.Finally,theeffectsofthemethodsinsegmentationarecompared.Theexperimentalresultsindicatethatthesystemcanobtainagoodperformanceofimagesegmentation.
Medicalimagesegmentationisaclassicalpuzzleforresearchers.Imagesegmentationisthemethodtodelineateanatomicstructuresorotherinterestedregionsautomaticallyorsemi-automatically,whichishelpfultodiagnosisandplaysacrucialroleinmanymedicalimagingapplications.
Keywords:
WaveletTransform;ImageSegmentation;threshold
第一章绪论1
1.1图像分割技术的现状和发展情况1
1.2图像分割主要研究方法1
边缘检测法2
区域提取法2
阈值分割法3
结合特定理论工具的分割法31.3论文的内容与结构安排...4第二章图像分割预处理52.1图像平滑5
中值滤波原理5
平滑效果分析6
2.2灰度调整7
灰度调整原理7
灰度调整效果分析7
2.3本章小结8
第三章基于阈值的图像分割技术93.1阈值分割原理9
3.2图像分割方法10
图像二值化10
双峰法10
最大方差自动取阈值法12
3.3本章小结13
第四章基于小波图像阈值分割技术144.1基于小波阈值分割技术简述144.2小波分析14
小波变换14
小波分割算法及步骤15
4.3阈值选取以及实验分析16
直方图分辨率的小波表示16
多分辨率阈值选取17
实验分析18
4.4本章小结20
第五章总结与展望21
5.1工作总结................................................21
5.2工作展望............................................21
致谢22
参考文献23
附录24
第一章绪论
本章对论文涉及的研究领域进行了较为详细的综述。
简要介绍了医学图像分割的研究目的和意义,给出了医学图像分割的基本方法及步骤。
在对医学图像分割问题的起源、发展和研究现状进行简要综述的基础上,介绍了该领域当前的研究热点及论文的主要研究内容。
1.1图像分割技术的现状和发展情况
图像分割算法的研究已有几十年的历史,一直以来都受到人们的高度重视。
关于图像分割的原理和方法国内外已有不少的论文发表,但一直以来没有一种分割方法适用于所有图像分割处理。
传统的图像分割方法存在着不足,不能满足人们的要求,为进一步的图像分析和理解带来了困难。
随着计算机技术的迅猛发展,及其相关技术的发展和成熟,结合图像增强等技术,能够在计算机上实现图像分割处理。
其中最主要的技术是图像分割技术,从图像中将某个特定区域与其它部分进行分离并提取出来的处理。
图像分割的方法有许多种,有阈值分割方法,边界分割方法,区域提取方法,结合特定理论工具的分割方法等。
早在1965年就有人提出检测边缘算子,边缘检测已产生不少经典算法。
越来越多的学者开始将数学形态学、模糊理论、遗传算法理论、分形理论和小波变换理论等研究成果运用到图像分割中,产生了结合特定数学方法和针对特殊图像分割的先进图像分割技术。
尤其是近年来迅速发展起来的小波理论为图像处理带来了新的理论和方法。
小波变换具有良好局部特性,当小波函数尺度较大时,抗噪声的能力强,当小波函数尺度较小时,提取图像细节的能力强,这样就可以很好地解决抑制噪声和提取图像边缘细节之间的矛盾。
1.2图像分割主要研究方法
图像分割是图像处理中的一项关键技术,自20世纪70年代起一直受到人们的高度重视,至今已提出了上千种各种类型的分割算法,现提出的分割算法大都是针对具体问题的,并没有一种适合于所有图像的通用分割算法,而且近年来每年都有上百篇相关研究报道发表。
然而,还没有制定出选择合适分割算法的标准,这给图像分割技术的应用带来许多实际问题。
因此,对图像分割的研究还在不断深入之中,是目前图像处理中研究的热点之一。
图像分割在图像工程中的位置它起着承上启下的作用,可以认为是介于低层次处理和高层次处理的中间层间。
最近几年又出现了许多新思路、新方法、或改进算法。
下面对一些经典传统方法作简要的概述。
多年来人们对图像分割提出了不同的解释和表述,借助集合概念对图像分割可给出如下定义:
令集合R代表整个图像区域,对R的图像分割可以看做是将R分成N个满足以下条件的非空子集R1,R2,R3,„,RN;1)在分割结果中,每个区域的像素有着相同的特性在分割结果中,不同子区域具有不同的特性,并且它们没有公共性分割的所有子区域的并集就是原来的图像各个子集是连通的区域边缘检测的实质是采用某种算法来提取出图像中对象与背景问的交界线。
我们将边缘定义为图像中灰度发生急剧变化的区域边界。
图像灰度的变化情况可以用图像灰度分布的梯度来反映,因此我们可以用局部图像微分技术来获得边缘检测算子。
经典的边缘检测方法,是通过对原始图像中像素的某小邻域构造边缘检测算子来达到检测边缘这一目的。
第二章图像分割预处理
由于受多种因素(光学系统失真、系统噪声、曝光不足或过量等)条件的影响,得到的图片往往信息微弱,无法辨识,需要进行增强处理。
增强的目的,在于提供一个满足一定要求的图像,或对图像进行变换,以进行分析。
2.1图像平滑
图像平滑的目的是为了减少图像噪声。
图像的噪声来自于多方面,有来自于系统外部干扰,如电磁波或经电源串进系统内部而引起的外部噪声,也有来自于系统内部的干扰,如摄像机的热噪声,电器机械运动而产生的抖动噪声等内部噪声。
因此,去除噪声,恢复原始图像是图像处理中的一个重要内容。
噪声主要来自下面三个方面:
(1)光电子噪声:
主要由光的统计本质和图像传感器的光电转换过程引起的(如光电管的光量子噪声和电子起伏噪声);
(2)电子噪声:
主要来自电子元器件(如电阻引起的热噪声);
(3)光学噪声:
主要由光学现象产生的(如胶片的粒状结构产生的颗粒噪声);
图像在生成和传输过程中受到这些噪声的干扰和影响,使图像处理结果变差。
因此,抑制或消除这些噪声从而改善图像质量,在图像处理过程中是一个重要的预处理,也称为对图像的平滑滤波过程。
中值滤波原理
中值滤波是一种非线性滤,由于它在实际运算过程中并不需要图像的统计特性,所以比较方便。
中值滤波首先是被应用在一维信号处理技术中,后来被二维图像信号处理技术所应用。
在一定的条件下,可以克服线性滤波器所带来的图像
细节模糊,而且对滤除脉冲干扰及图像扫描噪声最为有效。
(2-1)
对二维序列Xi,j的中值滤波,滤波窗口也是二维的,但这种二维窗口可以有各种不同的形状,如线状、方形、圆形、十字形、圆环形等。
二维中值滤波可表示为:
(2-2)
在实际使用窗口时,窗口的尺寸一般先用再取逐渐增大,直到其滤波效果满意为止G为输入信号频谱,F为输出信号频谱,定义为中值滤波器的频率响应特性,实现表明H是与G有关,呈不规则波动不大的曲线,其均值比较平坦,可以认为信号经中值滤波后,传输函数近似为1,即中值滤波对信号的频域影响不大,频谱基本不变。
平滑效果分析
图2-3所示的是图像中值滤波前后的效果比较,其中图2-3(a)是含有噪声的原图,图2-3(b)是用中值滤波处理后的图像,滤波窗口为3×3,可见,中值滤波后的图像不仅滤去了椒盐类噪声,而且边缘得到了较好的保护。
a带噪声图像b消噪后图像
图2-3带噪声图像与中值滤波后图像比较
2.2灰度调整
在成像过程中,扫描系统、光电转换系统中的很多的因素,如光照强弱、感光部件灵敏度、光学系统不均匀性、元器件特性不稳定等均可造成图像亮度分布的不均匀,导致某些部分亮,某些部分暗。
灰度调整就是在图像采集系统中对图像像素进行修正,使整幅图像成像均匀。
灰度调整原理
灰度调整可使图像动态范围增大,图像对比度扩展,图像变清晰,特征明显,是图像增强的重要手段之一。
在曝光不足或过度的情况下,图像灰度可能会局限在一个很小的范围内。
这时在显示器上看到的将是一个模糊不清、似乎没有灰度层次的图像。
采用线性灰度调整对图像每一个像素灰度作线性拉伸,将有效地改善图像视觉效果。
灰度调整效果分析
(a)灰度调整前b灰度调整后
c原始图像直方图(d)调整后直方图
图2-4灰度调整前后直方图比较
由图2-4可以看出b视觉效果较a明显,灰度调整前后直方图的比较可以看出,调整后直方图d去除了原始直方图c的噪声直方图,灰度调整后图像明显清晰了。
2.3本章小结
为了操作直观,本论文直接对灰度图像进行处理,要是彩色图像,须在操作前将其转换为灰度图像,本章通过对图像的平滑以及灰度调整,达到了比较明显的去噪效果。
在图像的直方图上也表现得十分明显,这样为后面进行基于直方图的操作提供了较好的条件。
第三章基于阈值的图像分割技术
当非灰度图像转换为灰度图像后,图像中各目标区域的灰度值会不一样,如果图像的灰度直方图具有明显的双峰值或多峰值特征,可以利用阈值化方法求取最佳阈值,然后对图像进行合理分割。
3.1阈值分割原理
阈值化图像分割是一种最基本的图像分割方法,经过半个多世纪的研究,现已提取了大量的算法。
其基本原理就是选取一个或多个处于灰度图像范围之中的灰度阈值,然后将图像中各个像素的灰度值与阈值比较,并根据比较的结果将图像中的对应像素分成两类或多类,从而把图像划分成互不重叠的区域集合,达成图像分割的目的。
采用阈值化图像分割时通常需要对图像作一定的模型假设。
利用图像模型尽可能了解图像有几个不同的区域组成。
基于图像分割模型经常采用这样一种假设:
目标或背景内相邻像素间的灰度值是相似的,但不同目标或背景的像素在灰度上存有差异。
设原始图像为fx,y,按照一定准则在fx,y中找到某种特征值,该特征值便是进行分割时的阈值T,或者找到某个合适的区域空间Ω,将图像分割成两个部分,分割后的图像为
(3-1)
对于有多种阈值情况,分割后的图像可以表示为:
其中是一组分割阈值,是经分割后对应不同区域的图像灰度值,K为分割后的区域或目
标数,。
无论是单阈值分割还是多阈值分割,都是选取一个比较合理的阈值,以确定图像中每个像素点应该属于目标区域还是背景区域,从而产生相应的二值图像。
))1,2,3,„,i,„L)i的像素数设为ni,则图像的全部像素数为
3-3
将其标准化后,像素数为,其中,i?
S,pi?
0,(3-4)
设有某一图像灰度直方图,t为分离两区域的阈值。
由直方图统计可被t分
离后的区域1、区域2占整图像的面积比以及整幅图像、区域1、区域2的平均灰度为:
区域1的面积比:
;区域2的面积比(3-5)
或者
整幅图像平均灰度;区域1的平均灰度;
区域2的平均灰度3-6
式中,G为图像的灰度级数。
整图像平均灰度与区域1、区域2平均灰度值之间的关系为
(3-7)
同一区域常常具有灰度相似特性,而不同区域之间则表现为明显的灰度差异,当被阈值t分离的两个区域之间灰度差较大时,两个区域的平均灰度u1,u2与整图像平均灰度u之差也较大,区域间的方差就是描述这种差异的有效参数,其表达式为:
(3-8)
式中,表示了图像被阈值t分割后的两个阈值之间的方差。
显然不同的t值,就会得到不同的区域方差,也就是说,区域方差、区域1均值、区域2均值、区域面积比、区域面积比都是阈值t的函数,因此式(3-8)可写为:
(3-9)
经数学推导,区域间的方差可表示为:
(3-10)
被分割的两区域间的方差达最大时,被认为是两区域的最佳分离状态,由此确定阈定值T:
,以最大方差决定阈值不需要认为设定其他参数,是一种自动选
择阈值的方法,它不仅适用于两区域的单阈值选择,也可以扩展到多区域的多阈值选择中去。
a原始图像b最大方差法分割后图像
图3-11最大方差自动取阈值法
该方法将图像分成两个类,当类间方差与类内方差的分离度最大时即为最佳阈值(由图3-11表明,该方法能够准确而快速地对图像进行二值化,特别是当对象物和背景的灰度值的差具有一定大小的时候,效果更明显。
3.3本章小结
双峰法和最大类方差法区域分割技术,是图像分割中最重要而且有效的技术之一,在实际的图像处理系统中得到了广泛应用。
特别是在需要实时性较强的图像处理系统中,快速而准确的图像阈值化方法就成为非常重要的研究目标。
对一幅具体的图像,选用何种算法,要进行对比实验,不存在一种通用的图像分割算法。
而这些都只是传统的分割方法,只是全局阈值分割方法中较好的方法,但是对图像进行局部阈值分割上面方法显得尤为困难,下面将介绍小波变换的分割技术,也是本论文的重点。
第四章基于小波图像阈值分割技术
小波变换是近年来得到广泛应用的数学工具,与傅里叶变换、窗口傅里叶变换相比,小波变换是空间时间和频率的局域变换,能有效地从信号中提取信息。
4.1基于小波阈值分割技术简述
本论文利用小波变换对含噪图像的直方图进行多尺度分解,先在较大的尺度下找出图像分割阈值的粗略值,然后逐渐减小尺度,精确定位分割阈值,算法采用MATLAB编程仿真。
基于小波变换的阈值法图像分割技术则能够有效地避免噪声的影响。
该方法的基本思想是首先由二进制小波变换将图像的直方图分解为不同层次的小波系数,然后依据给定的分割准则和小波系数选择阈值门限,最后利用阈值标出图像分割的区域。
整个分割过程是从粗到细,由尺度变化来控制,即起始分割由粗略的L2(R)子空间上投影的直方图来实现,如果分割不理想,则利用直方图在精细的子空间上的小波系数逐步细化图像分割。
4.2小波分析
基于小波变换的阈值法图像分割技术能有效地弥补传统的图像阈值法分割技术的不足,具有较强的抗噪声性能,同时,对于直方图为多峰值的情况,可以利用小波的多分辨率分解,对灰度阈值进行合理地选择,实现对图像的分割处理。
换
由于图像的直方图可以看作是一维信号,而直方图上的突变点(波峰点和波谷点),往往可以代表图像灰度变化的特征。
因此Jean-ChristopheOlivo提出了用小波变换对直方图进行处理的方法实现自动阈值提取。
Olivo通过检测直方图小波变换的奇异点和区域极值点给出直方图峰值点的特性。
而小波变换的波峰
和波谷点可以代表图像中灰度代表值和阈值点。
利用小波变换多尺度特性实现对图像的阈值分割。
又由于小波变换具有多分辨率的特性,因此可以通过对医学图像直方图的小波变换,实现由粗到细的多层次结构的阈值分割。
首先在最低分辨率一层进行,然后逐渐向高层推进。
小波变换的零交叉点表示了在分辨率2j时低通信号的局部跳变点。
当尺度2j减小时,信号的局部微小细节逐渐增多,因此,能够检测出各微小细节的灰度突变点;当尺度2j增大时,信号的局部细节逐渐消失,而结构较大的轮廓却能清晰地反映出来,因而能检测出该结构较大的灰度突变点。
因此,可以选择小波为光滑函数的二阶导数,对图像的一维直方图信号进行小波变换,检测出直方图信号的突变点,由此搜索出两峰之问的谷点作为分割阈值点。
这就是小波变换用于图像分割的基本原理。
对图像的直方图来说,它的各层的小波分解系数表示不同分辨率下的细节信号,它与小波近似信号联合构成直方图的多分辨率小波分解表示。
给定直方图,考虑其多分辨率小波分解表示的零交叉点和极值点来确定直方图的峰值点和谷点。
小波分割算法及步骤
分割算法的计算量与图像尺寸大小呈线性变化,本论文介绍直方图的多分辨率分析。
对于每个整数j?
Z(Z整数集合),表示在j分辨率下的二进制有理数。
因此,对于任何j?
Z,是一组在实数轴上的等间隔采样点集合,如果ij,则表示低分辨率(较粗)的采样点;反之,ij,则表示高分辨率(较细)的采样点。
假定f表示为一幅图像,g是图像f中最大灰度,则直方图表示为
(4-1)
式中表示计数操作,是离散函数。
令,离散函数表示成连续函数,看作是由几个分段常数函数组成。
对于j?
Z,按采样点采样,则表示在j分辨率下的直
方图。
进一步可以用Haar尺度函数的平移与伸缩表示,即
(4-2)
(4-3)
多层表达曲线可以表示如下
(4-4)
对直方图进行小波分解,利用小波系数,按式(4-4)重建直方图,从近似直方图中选择阈值,完成分割阈值。
步骤如下:
步骤1:
预设分割区域为M,分解级数,L为图像最大灰度值;
步骤2:
小波分解曲线,得到,令j0,;
步骤3:
,将大于j分解层次的系数设置为0,用式(4-4)重建,在重建直方图中,找出满足和条件的标号l(灰度),并且统计标号l的个数n;
步骤4:
如果nM,则jj+1,当jJ时,转向步骤:
3;
步骤5:
从重建直方图中,找到阈值;
步骤6:
像素值与阈值比较,标出所在区域。
4.3阈值选取以及实验分析
本论文所采用波谷点确定为图像分割的阈值点,两阈值平均点作为后一阈值和前一阈值之间区间灰度的代表值。
直方图分辨率的小波表示
设图像的灰度范围为0,1,2,„,N-1,灰度值x0xN-1对应的像素为n0,则一幅图像总像素为M:
(4-5)
灰度值x出现的概率为:
(4-6)
由上式可以建立该图像的直方图,它反映了该图HxPx,xO,1,„,N-1上灰度分布的统计特性,是基于像素灰度的图像分割方法的基础。
为了建立小波变换的多分辨率分解表示,引入尺度函数x,其傅立叶变换
满足条件:
(4-7)
可见,x相当于低通滤波器,这样图像直方图Hx的低通分量为:
(4-8)
设原始图像直方图信号各尺度之间的各阶小波变换。
可以证明:
信号在在尺度为时被平滑掉的高频成分,可以用尺度的小波变换来恢复,我们称集合为图像直方图信号的多分辨率小波分解表示。
直方图信号多分辨小波分解由一个最低分辨率下的近似信号和一组分辨率的细节信号所组成。
这是一种介于频域和时域的表示。
为图像分析提供了一个由粗到细的分层框架。
多分辨率阈值选取
基于直方图和小波变换的图像分割技术由以下几个步骤组成:
首先由粗分辨率下的图像直方图细节信息确定分割区域类数;其次,在相邻峰之间自动确定最优阈值;最后用求出的最优阈值分割原图像。
由于图像的原始直方图一般不够平滑或含有一定的噪声,因此,有必要对原始直方图进行平滑处理,以利于分割目标。
其方法为:
在空间域中采用保护边缘平滑方法平滑直方图,它既能保留原直方图基本变化特性,又能消除小峰的跳动。
或者选取大尺度下的小波变换系数对直方图进行处理,也可以减小噪声的影响。
在分辨率为2j时,由小波分解后的直方图近似信号的极大值确定初始区域
类数,即确定峰的数目。
对于灰度级数不多的原始影像,一个区域类通常对应直方图中的一个峰,然而,对于一幅复杂图像,经小波分解后平滑直方图中的每个峰则不一定都对应一个区域类,它也可能从属于邻近的一个峰,因而有必要通过检查认定哪些峰对应于分割区域类。
峰的独立性判断是为了消除不能成为一类的那些峰,独立峰应满足三个条件:
1应具有一定的灰度范围;
2应具有一定峰下面积;
3应具有一定的峰谷差。
峰与峰之间的谷点选取首先在最低分辨率层进行,然后逐层推进,直到信号的最高分辨率层,并在每一层的阈值选取中,采用前一层的选取结构为引导。
这种由粗到精的控制策略,能有效地选取最优分割阈值,同时较好地克服了噪声干扰和搜索空间大的问题。
设最优阈值分别为为T1,T2,T3,„,Tkk为正整数,即可用这些最优阈值分割原始图像fx,y,得到分割结果gx,y,公式如下:
(4-9)
其中CkkO,1,(((,K表示分割后的类别代码。
实验分析
(a)))
一、全局二值化方法由于采用的是用一个固定门限值来分割,因此门限值的选取十分重要。
虽然众多学者提出了许多种选取“最优门限”的方法,但这种分割方法在光照不
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 毕业设计 论文 matlab 图像 分割