一种基于图的彩色图像分割算法-沃焱.pdf
《一种基于图的彩色图像分割算法-沃焱.pdf》由会员分享,可在线阅读,更多相关《一种基于图的彩色图像分割算法-沃焱.pdf(8页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、第44卷第9期2016年9月华南理工大学学报(自然科学版)Joumal of South ChiIla UIliVersi哆of Technology(Namral Science Edition)V0144 No9S印tember 2016文章编号:1000565x(2016)09-0001-08一种基于图的彩色图像分割算法水沃焱 金璇(华南理工大学计算机科学与工程学院,广东广州510006)摘要:为了解决图像分割中容易出现的过分割问题,提出了一种基于图的彩色图像分割算法该算法在区域合并的基础上,首先用Mean shift方法进行预处理,得到初始过分割区域后对其构造邻接图,然后计算邻接区域的
2、颜色、纹理及边缘特征相似性以判断区域是否需要合并直到所有满足条件的区域都被合并为了保持图像的全局属性,文中通过查找最优合并成本的方式进行区域合并实验结果表明:即使在图像目标和背景区域颜色比较相似时,文中算法也能较好地实现对目标区域的完整分割;与其他4种算法相比,文中算法具有更好的分割性能关键词:图像分割;Gabor滤波器;纹理;边缘;区域相似性中图分类号:TP39141 doi:103969jissn1000565x201609001图像分割的目的是将整个图像分成若干个互不相交的非空子区域,每个子区域的内部是连通的,同一区域在亮度、颜色等特征上具有相似性,而相邻区域之间不具有相似性分割问题是图
3、像处理、计算机视觉中一个基础性和挑战性的问题,在图像识别、目标检测等技术中有重要的作用该问题在20世纪70年代就受到了广泛的关注,人们提出了多种图像分割方法,如阈值方法J、聚类方法拉。J、图论方法Hj J、基于区域的方法帕。1等而基于图的分割方法旧。纠因能考虑到所有像素点携带的特征和空间信息而备受关注在构造的图中,图顶点表示图像的像素点或区域,顶点之间用边连接,边的权值可以表示像素点或区域之间的不相似性最早提出基于图的分割方法是Zahn旧。的基于最小生成树方法该方法的分割标准是消除最小生成树中权值较大的边,权值为像素点之间的亮度差异,由于突变区域亮度大于稳定区域亮度,采用固定阈值和局部测量方式
4、进行分割可能使突变区域被分成多个区域,或亮度较低区域被合并为解决这个问题,Urcluhanl91提出了用最小权值对边进行归一化的分割方法,此方法容易导致错误分割,并且难以保持图像的全局属性Felzensw舶等叫提出了通过计算区域内部和区域之间的颜色差异来预测边界是否存在,从而进行区域合并在这些以像素点为顶点的的图方法中,图像越大,顶点越多,且每个顶点需与周围的多个顶点建立邻接关系,使得图像处理复杂度较高为解决这个问题,Rezvani陆等1先用Mean shift方法进行预处理,然后结合Ncut和AverageCut进行图像分割B0等副先用分水岭分割,再结合连续概率比测试和最大可能性标准进行合并
5、预测,并动态合并满足条件收稿日期:20151120。基金项目:国家自然科学基金资助项目(61472145);广东省自然科学基金资助项目(2016A030313472);华南理工大学中央高校基本科研业务费专项资金资助项目(2015zZ031);广东省一教育部产学研合作专项(20138090500015)Fo哪da6伽items:Supported by the National Natural Science Foundation of China(61472145),the Natural Science Foundation ofGuangdong Pmvince(2016。A0303 l
6、3472)龃d the Special fbnds for UniversityIndustry Co叩eration of Guangdong Province and theMinistry of Education of China(20138090500015)作者简介:沃焱(1975一),女,博士,教授,主要从事多媒体应用技术研究Email:woy龃scuteducn万方数据2 华南理工大学学报(自然科学版) 第44卷的区域这些方法中,图的顶点表示一个小的区域,每个顶点只需与周围少部分顶点建立连接,提高了图像处理速度但这些方法仅根据颜色进行分割,当图像目标和背景区域颜色相似时,难以得
7、到好的分割结果针对这些问题,文中结合图像的颜色、纹理和边缘特征,提出了一种基于图的彩色图像分割算法为提高图像的分割速度,文中先使用Mean sh“”1进行预处理得到初始过分割区域,并利用这些区域构造相应的邻接图为保持图像的分割完整性,文中提取了图像的边缘轮廓,并计算邻接区域的边缘相似性为减少图像目标和背景颜色相似时的过分割现象,文中利用Gabor滤波提取图像的纹理特征,并通过统计邻接区域的纹理分布来计算其纹理相似性为了保持图像的全局属性,文中根据查找最优合并成本的方式进行区域合并1 区域相似性度量文中利用区域合并来分割图像由于区域比像素点携带更多的信息,并且基于区域的处理方法能提高图像分割效率
8、,因此文中先用Mean shift预处理得到初始过分割区域集合JR=R。,尺:,氐,然后通过度量区域相似性将相似的区域进行合并文中从颜色、纹理和边缘3个特征来度量区域相似性11 颜色相似性为度量区域颜色相似性,先将图像分成R、G、B三个通道,邻接区域R。与Ri(R。,尺iR)的颜色差异D珊(Ri,Ri)可表示为Diff(R,弓)=(1)式中,p。,Ri,p,Rf,cf(m)、cf(p,)分别为像素点n、p,在第f通道的颜色值,、i分别为区域Ri、Ri的像素点个数Diff(尺。,尺f)越大,R、尺f颜色相似的概率越小;反之,其颜色相似的概率越大则邻接区域R。、Ri颜色相似性条件可表示为cl(屯驴
9、告警如鹰鼠 (2)占,为阈值通过比较占,与区域颜色差异大小来判断颜色相似性,因此s。的取值应与图像颜色差异及分布相关当区域颜色差异较大时,图像颜色分布比较分散,为避免过分割现象,s,需相应地增大;当区域差异较小时,为避免欠分割现象,s,需相应地减小由于颜色的标准差能反映颜色的分布情况,因此文中利用图像颜色的标准差来计算占,:81=2(stdR+stdG+stdB)3max(Dif)(3)式中,maDc(Diff)为邻接区域对的最大颜色差异,stdR、stdG、stdB分别为图像在R、G、B通道归一化后的标准差stdR、stdG、stdB越小,说明图像颜色分布越集中,颜色差异越小,对应的占,值也
10、越小;反之,stdR、stdG、stdB越大时,说明颜色分布越分散,对应的颜色差异越大,s,值也越大12边缘相似性为了较好地保证分割图像的完整性,合并前需考虑区域间是否有较强的边缘轮廓Arbelaez等4 o利用局部多尺度信息和全局边缘信息提出了一种有效的提取边缘轮廓的方法局部多尺度边缘信息和全局边缘信息的结合,既能提取出图像的边缘轮廓,又能凸显出显著的边缘,有效地对图像进行描述基于这些优点,文中使用该方法来提取图像的边缘轮廓121边缘轮廓提取对于局部边缘提取,该方法通过计算在LAB颜色空间中的L、a、b通道以及纹理通道的方向梯度信号G(并,y,p)来预测像素点在p方向存在边缘的概率对于L、a
11、、b通道的方向梯度,该方法先在像素点p 7(z,y)处放置一圆盘,并在p方向将圆盘划分成两半,然后求出两个半圆盘的亮度或者颜色直方图分布两个半圆盘直方图、的z2距离就是像素点p 7的方向梯度: 疋2(m,)-专耕(4)式中,。、:分别为两个半圆盘直方图和中第n个桶的频数每个桶表示量化后的亮度或者颜色强度对纹理通道进行相应处理后可用以上方式求出像素点的纹理方向梯度为了检测到更多轮廓,该方法从盯2,盯,2盯3个尺度分别计算每个通道的方向梯度信号多尺度多通道的方向梯度信号的结合,形成局部多尺度边缘信息mPb:mPb(石,y,口)=a ch,。G ch出)(x,y,p)(5)式中:mPb(戈,y,日)
12、为像素点在(戈,y)处日方向的局部多尺度边缘概率;s为尺度;G chm)(x,y,p)是颜色通道ch上以(戈、),)为圆心、盯(ch,s)为圆盘半径,从角度p上进行切分的两个半圆盘直方图的差万方数据第9期 沃焱等:一种基于图的彩色图像分割算法异;am;为权值,用来衡量不同通道尺度下的概率分布像素点的局部多尺度边缘能量可表示为Empb(x,),)=ma)【mPb(z,y,p) (6)对于全局边缘检测,可用谱聚类将局部多尺度信息结合到全局框架首先根据局部多尺度边缘概率定义相似矩阵彬,l,=exp(一呻蜂Empb(p)p) (7)式中,驴为连接两个像素点i、,的线段,Empb(p 7)为线段经过的像
13、素点p 7的Empb,常数p取值为01对角矩阵D。=乞Wi,根据(Dw),=加v获得特征向量v。,y1,一,l,把每个特征向量yt看成一个图像,在p方向上与高斯方向导数滤波器作卷积,可获得方向梯度信号V。v。(戈,y)不同特征向量的方向梯度的结合,构成谱聚类下的边缘梯度概率:1sPb(戈,y,咿)=_lIV。V。(戈,y) (8)=1、女式中,l历为权值,A。为第后个特征向量的特征值结合多尺度局部梯度信息和谱聚类全局梯度信息可获得图像的边缘轮廓gPb:gPb(戈,y,p)=卢“;G。h,(。,)(戈,),臼)+s ch“sPb(z,臼) (9)式中,卢。和p为权值对于像素点的全局边缘能量,可以
14、用像素点的最大gPb边缘概率来表示:Egpb(戈,y)=m9xgPb(戈,y,p) (10)122边缘相似性判断由于gPb用边缘梯度信号描述,梯度越大,存在边缘的可能性越大而邻接区域的边缘相似性主要在于检测邻接区域边界之间是否存在边缘,因此文中根据邻接区域间相邻边界点的梯度概率来度量边缘相似性在图1所示的示例中,尺。、Ri为邻接区域Ri、R,的相邻边界点集合P=p。l 1mP,P为Ri、足之间的相邻边界点个数根据边界点的梯度幅值可将咒、尺i间的边缘梯度概率Grad(月i,月i)表示为4吣墨图1 邻接K域分布同Fig1 Dist曲ution of adjacent regions边界点集合P的边
15、缘概率平均值:Grad(尺。,弓)=Egpb(p。(z),Jp。(y)P(11)式中,p。(戈)、p。(y)表示第m个边界点的戈、),坐标Gmd(Ri,R,)越大,尺i、尺,之间存在边缘的可能性越大,两区域边缘特征越不相似;反之,Grad(R,Rj)越小,Ri、R,之间存在边缘的可能性越小,边缘相似的概率越大因此,Ri、尺,的边缘相似性条件可表示为c2(屯耻矗罴汨n剐妇2 (12)c2一B卜ii其他一一 02式中,占:为阈值当区域边界点的梯度概率小于s:时,邻接区域具有边缘相似性;反之,则不具有边缘相似性13纹理相似性纹理是图像处理中的常用特征在文献1012方法中,仅用颜色特征来进行图像分割,
16、当背景和目标区域颜色相似时,难以得到较好的分割结果,因此在合并时需要考虑区域间的纹理相似性Gabor滤波15 3是一种常用的提取纹理的方法,能有效提取不同方向和尺度的纹理信息因此,文中用Gabor滤波器来提取纹理特征二维Gabor滤波器表示为 g(戈,),;痧,8,妒,6,7)=exp(一等)exp(i(2丌芳叫) (13)式中,zcosp+),sinp,y=一戈sinp+),cosp,毋为扬乏长,9 7为Gabor核函数方向,沙为相位偏移,6为高斯标准差,y为空间纵横比由于纹理有较强的方向性,因此文中从8个方向(o,詈,詈,荨,詈,警,莩,)提取纹理先利用式(13)得到8个Gabor模板,然
17、后用Gabor模板对图像进行滤波,可获得各方向的的纹理特征:=J,0踟j (14)式中,L,为p方向的纹理特征,O为卷积运算符,为原图,骱为p方向的Gabor模板由于像素点在各方向上的纹理梯度不同,为度量区域纹理相似性,需要统计区域中各像素点的纹理分布像素点在p方向上的纹理梯度越大,该方向的纹理越强因此,可以把像素点的最大纹理梯度所在方向作为该点的主方向,记为p:(戈,y),纹理梯度作为该点的能量值E(戈,_y):万方数据4 华南理工大学学报(自然科学版) 第44卷E(z,y)2m戕m(x,y) (15)式中,L,(戈,y)为像素点(z,),)在p 7方向上的纹理梯度选出所有点的主方向和能量值
18、后,再用方向直方图进行纹理统计初始时,直方图在各方向上的纹理能量为0对区域足的所有像素点,为直方图的各个方向投票,每票都带有权值,当像素点的主方向酢(戈,),)对应于直方图的p方向时,权值为该像素点的能量值E(戈,y),否则为0则像素点在p 7方向上的权值E(戈,y,臼7)可表示为Ecz,y,p,=:jxy :三:!:; c-6,当Ri中所有像素点投完票,可获得图2(a)所示的纹理分布直方图在p 7方向上的能量分布可表示为。t尺i(目)=E 7(舅(露),y(七),p) (17)式中,E 7(戈(后),),(后),p)为R。中第1|个像素点在p方向上的权值同理,对区域Ri纹理统计后,可得到图2
19、(b)所示的纹理分布讪蛰-_2 3 4 5 6 7 8方向索引l 2 3 4 5 6 7 8方向索引(a)片纹唯分布 (h)t?Z理分仿图2图l中R。、置的区域纹理分布Fig2 Texture distribution of re舀ons R。and足in Fig1统计邻接区域足、兄i的纹理分布后,其纹理差异可以表示为1厂r一tDist(尺i,R)2寺荟tRi(9,)-t弓(口,)2(18)tDist(尺。,尺,)越大,R。与尺,纹理相似的概率越小,反之纹理相似的概率越大因此,区域Ri、尺i的纹理相似条件可表示为c3(驴盅巍一即妇3 (19)c舭i,弓)2 ii其他一1。 (19)占,为参数当
20、两区域之问的纹理差异小于占,时,两个区域在纹理上有相似性,反之不具有相似性14 区域合并条件合并条件是区域合并的关键文中度量了图像颜色、纹理以及边缘特征的相似性,如果仅把颜色相似性作为合并条件,当目标和背景区域颜色相似时,容易将目标和背景区域融合同理,如果仅把纹理或者边缘相似性作为最终合并条件,也难以获得较好的分割结果若利用颜色和纹理或者颜色和边缘相似性作为合并条件,与之前单一特征相比,合并效果应该要好,但如果将颜色、纹理和边缘特征相结合,使得合并条件更加严格,既能使相似的区域被合并,又能够有效地减少过分割现象因此,文中结合颜色、纹理以及边缘特征相似性,将构造的区域合并条件表示为C 7(尺i,
21、R,)=C。(尺。,Ri)C2(尺。,尺,)C3(R。,尺j)(20)仅当C(R。,R,)=1时才进行合并,此时c。(Ri,尺i)、C:(Ri,尺,)、C,(R。,尺,)均为1,表明区域Ri、Rj在颜色、纹理和边缘特征上都具有相似性2 基于图的区域合并21 邻接图的构造为了更好地描述图像,文中对Mean sh甜131产生的过分割区域集合R构建邻接图令G 7=(y,e)表示一个无向图,y是图的顶点集合,对图进行初始化时y=R(秽i,a,)e,e是边集合,(秽。,)表示连接顶点秽i和秽i的边,该边的上的权值伽(口i,口,)用于衡量秽i、移i的不相似性结合区域在颜色、纹理和边缘特征上的差异,可将连接
22、口;、Fi的边的权值甜(耽,秽表示为埘(。,秽,)=Diff(秽。,秽,)+Grad(秽。,可i)+tDist(创i,q) (21)22 区域合并过程在区域合并过程中,合并顺序的优劣能反映算法能否保持图像全局属性由于区域之间的差异性,合并时会产生合并成本,记为costF,差异越大,costF越大为使合并顺序最优,每次合并时costF应最小由于边的权值能反映区域间的差异,因此区域Ri和Ri之间的costF可用权值表示,即costF(Rj,尺,)=埘(尺。,R,) (22)在合并过程中,选择costF最小的边(Ri,尺i),若R。、尺,满足合并条件就将Ri、尺i合并,并更新与R。、Ri相关联的边若
23、尺i、Ri不满足合并条件,因其他边的costF比costF(R。,R,)大,表明其他区域间的颜色、纹理以及边缘上的特征差异更大,因而不会满足合并条件,此时区域合并可以终止万方数据第9期 沃焱等:一种基于图的彩色图像分割算法 5S翡菁霎溉邻接图中所有边的合并3实验分析1)利用式(22)计算出邻接图中所有边的合并 。 死但卅1,I成本costF2)在边集合e中选择合并成本costF最小的边(尺i,Ri),利用式(20)计算c 7(尺i,R小3)若C(Ri,尺i)=1,则将Ri合并人Ri,并对邻接图作以下更新(a)删除Ri与Ri间的边(Ri,R山(b)对于所有从Ri出发的边(尺i,R触),删除边(R
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 一种 基于 彩色 图像 分割 算法 沃焱
限制150内