一种改进的活动区域轮廓模型——无需水平集重新初始化
- 格式:pdf
- 大小:306.27 KB
- 文档页数:5
8 2008.44(2) Computer Engineering and Applications计算机工程与应用
一种改进的活动区域轮廓模型
——无需水平集重新初始化
孙文杰 ,陈允杰 ,汤 杨 ,韦志辉 ,王平安 ,夏德深
SUN Wen-jie .CHEN Yun-jie1,TANG YangI,WEI Zhi-huP,WANG Ping—an2,XIA De—shen
1.南京理工大学计算机科学与技术学院,南京2 10094
2.香港中文大学计算机科学与工程学系,香港
1.School of Computer Science&Technology,Nanjing University of Science and Technology,Nanjing 210094,China
2.Department of Computer Science&Engineering,Chinese University of Hong Kong,Hong Kong,China
E—mail:wenjie—s@163.coln
SUN Wen—jie,CHEN Yun-jie,TANG Yang,et a1.New active region contour model without re—initialization.Computer
Engineering and Applications,2008。44(2): 1 l_
Abstract ̄Active region contour model has been successfully applied in image segmentation,object tracking,etc,which has many pros compared with active contour model based on edge.But,during the evolution of level set,it’s numerically necessary to keep
the evolving level set function close to a signed distance function.Re—initialization,a technique for periodically re—initializing the
level set function to a signed distance function during the evolution,results in slowing speed of evolution,increasing complex of implement.To overcome the re—initialization,the author proposes to add a panelized energy into active region contour model to
keep level set close to a signed distance function;it thus speeds up the CHIVe evolution and the segmentation.The auhtor applies
the model to segment the texture image,brain MRI and track objects;the experiments prove the model is effective. Key words:level set;re-initialization;CV model;active region contour model;image segmentation
摘要:基于区域的活动区域模型已经成功应用在图像分割、目标跟踪等领域,较之基于梯度的活动轮廓模型具有很多优点。但
是,这些水平集模型在演化过程中,为了保持为符号距离函数,必须对其重新初始化,降低了曲线演化速度,增加了实现复杂度。为
了解决重新初始化问题,在测地活动区域模型的能量函数中,加入惩罚项来约束水平集保持为符号距离函数,无需再重新初始化,
极大地提高了演化速度。将其运用在纹理图像、脑MR图像分割以及视频跟踪中,实验证明该模型是有效的。
关键词:水平集;重新初始化;CV模型;活动区域轮廓模型;图像分割
文章编号:1002—8331(2008)02—0008—04 文献标识码:A 中图分类号:TP391.4
l引言
水平集方法(Level Set Methods)是s.Osher和J.A.Sethian 于1988年在研究曲线以曲率相关的速度演化时提出来的1]1,用
于描述曲线的演化过程。其基本思想是将当前正在演化的闭合 曲线(曲面)的问题转化为更高维空间中的水平集函数曲线(曲
面)的隐式解。与其他曲线演化方法相比,它的最大优势在于它
的稳定性和拓扑无关性。
活动轮廓模型与水平集理论相结合的曲线演化方法,是目
前广为关注的一种分割方法。该方法利用闭合曲线或曲面的几
何特性,建立曲线变形的能量函数,将这个能量函数最小化,使
闭合曲线逐渐逼近图像中目标边界埘。曲线演化方法又可分为
边界活动轮廓模型口1、区域活动轮廓模型f 1和混合活动轮廓 模型 3种。基于边界的活动轮廓图像分割技术优于传统的基
于边缘检测的分割,但由于它也是建立在图像边缘检测基础上
的,因而与传统方法具有相同的局限性,即对噪声较敏感。而基
于区域的主动轮廓图像分割技术则由于直接使用主动轮廓内
部和外部的全部灰度信息,没有使用基于梯度的边缘检测函
数,故对噪声较大的图像和弱边界图像仍能工作得很好。
无论是基于梯度的还是区域的模型,为了保持水平集接近
符号距离函数(Signed Distance Function,SDF),都要进行水平
集重新初始化运算,大大增大了计算量和实现的复杂度。
Gomes和Faugems 曾指出,重新初始化使得水平集理论和应 用存在不一致现象,而且当前许多初始化的方法在初始化的时
候会产生偏离真实水平集现象。水平集重新初始化还仍是一个
基金项日:香港特区政府研究资助局研究项目(No.CUHK/4185/00E);香港中文大学研究基金(No.2050345 o 作者简介:孙文杰(1977一),女,博士研究生,主要研究方向为图像处理、模式识别、数值分析与算法;陈允杰(1980一),男,博士研究生,主要研究方 向为图像处理、模式识别、数值分析与算法;汤杨(1980一),男,博士,研究方向为图像处理、模式识别、数值分析与算法;韦志辉(1963一), 男,博士,教授,研究方向为变分偏微分方程在图像处理中的应用、小波分析、模式识别;王平安(1961一),男,博士,教授,主要研究方向
为交互式科学计算可视化、三维医学图像、虚拟现实等;夏德深(1941一),男,博士,教授,博士生导师,主要研究方向为图像处理、卫星遥 感、模式识别。 一
维普资讯 http://www.cqvip.com 孙文杰,陈允杰,汤杨,等:一种改进的活动区域轮廓模型——无需水平集重新初始化2008,44(2) 9
难题。 ¨在文献[7]中因此提出了新的无需重新初始化的水平集
模型,实验证明¨方法实现简单,效率较高,具有很大的优越
性。但¨的方法是基于梯度的边界活动轮廓模型,对于弱边界 和内部有空洞时将无能为力,另外对初始曲线的要求也较高。
为了克服以上缺点,将n的无需初始化方法进行推广,本文将
其应用到活动区域轮廓模型中。 本文第2章介绍区域活动轮廓模型框架,第3章介绍b
的无需初始化思想,第4章将无需初始化和区域的活动轮摩 模
型结合得到改进的分割模型,第5章给出实验结果和分析,最
后进行总结。
2区域活动轮廓模型框架
假定图像n中有目标 和背景 ,两部分构成, = u
,1-1 n 1-12=0。根据最大后验概率准则(Maximum A—posterior
Probability,MAP),图像分割的任务就是根据数据D寻找模式
分类 ,最大化后验概率P(MID)
P(MID)=—P(D IM )P(M) (1)
由于P(D)是一定值,因此最大化P(MID)即最大化: 、
E(M)=P(DIM)P(M) (2) 也即最小化:
E( )一log(P(DIM))一log(P(M)) (3)
其中,log是以e为底的自然对数。应用到图像中,就是模式
就是当前点属于n 还是 ,数据D就是图像,,最大化后验概
率P(MID)可写成
11 ̄II㈩ = ㈩
记P,(s)=P(1(x)=sIx∈n.),根据轮廓曲线的长度估计当前点的
先验概率的p(x n,),假设最后分割后的轮廓曲线长度最短
并假设图像中每个点是相互独立的,式(2)可写成
np ( )np:( )e (5) ∈{ 】 ∈{ 式(3)可写成
E(,)=一J logp ( ) 一J 1ogpz( )dx+v J,ds (6)
,为分割曲线,应用到水平集框架下得到
E((b)=一J H((b)l。gP ( ) 一J (1一日((b))logpz( )(k+
I IVH((b)ldx (7)
可以看到,若内外区域的分布满足高斯分布,采用高斯函数估计P
. P(s)=— 一e (8) 、/2订 采用梯度下降法求解式(7),得到
,( (1。gP1--10gPz T )) (9)
若采用高斯模型来估计概率,每次更新式(9)后,和CV模型一
样需要更新内外区域的均值 , :,标准差 , :
1= 坍(西)dx H(西)dx , 2= I ,(1-H(+))dx J n
l(1-H(西))dx J n (10) 当 = 时,带入式(7)得到的就是CV模型的能量函数式:
E((b)=J ( )2H((b)dx+J ( z) (1一日((b)) +
Il VH((b)l J n (12)
因此,CV模型可看作是式(7)的特例。很明显CV模型是假设 图像是分段常量的函数,当内外区域的均值一样时,将得不到
正确的分割结果。如图1所示:图1(a)是一幅内外区域具有相 同均值,不同方差的实验图,图1(b)是CV模型得到的错误分
割结果,图1(c)是采用式(7),采用高斯概率估计得到的正确
分割结果,较之CV模型更具实用性。
■一一
(a)原图 (b)CV分割结果(c
图l合成图像分割
但由于高斯分布只能处理均值不同或方差不同的情形,因
此当两类目标具有相同均值和相同方差时,将导致分割失败。
并且由于图像的分布不一定满足高斯分布如图7(d)的斑马图 像,采用高斯分布来进行概率估计,和CV模型一样也有局限
性。图2(a)是一幅内外区域具有相同均值,相同方差的实验
图。图2(b)是真实的分割结果,图2(c)是采用高斯分布得到的
错误分割结果,图2(d)是本文模型得到的分割结果。
■■
(a)原图 (b)真实分割结果
●■
c)采用高斯概率结果 (d)本文模型结果
图2合成图像分割
文献[8]中采用调整窗口尺寸参数的Parzen窗来估计概率