基于Snake模型的图像目标轮廓自动跟踪方法_精品文档.pdf

1、第 33 卷第 2 期2003 年 3 月 东 南 大学 学 报(自 然 科 学 版)JOURNAL OF SOUTHEAST UNIVERSITY(Natural Science Edition)Vol133 No12Mar.2003 基于 Snake 模型的图像目标轮廓自动跟踪方法王立功1 于甬华2 姜晓彤1 罗立民1(1东南大学生物科学与医学工程系,南京 210096)(2山东省肿瘤医院放疗科,济南 250117)摘要:针对实验数据像素灰度的分布特点,提出了一种对目标轮廓线进行有效和可靠的搜索和跟踪策略.由于数据中病变组织与其邻近组织像素灰度差别相对明显,首先通过采用一种改进的轮廓自动跟

2、踪方法对目标轮廓进行跟踪,将得到的轮廓线经采样得到其离散控制点作为 Snake轮廓搜索和跟踪算法的输入,既克服了 Snake 方法对初始轮廓线控制点分布的局限性,又避免了采用单一轮廓跟踪方法跟踪目标轮廓线的不确定性,提高了分割病变组织的速度和准确性,此方法具有较高实用意义.关键词:Snake 模型;轮廓;目标跟踪;方向编码;交互分割中图分类号:TP391 文献标识码:A 文章编号:1001-0505(2003)02-0215-04Snake model based automatic tracking methodfor image object contourWang Ligong1 Yu

3、Yonghua2 Jiang Xiaotong1 Luo Limin1(1Department of Biological Science and Medical Engineering,Southeast University,Nanjing 210096,China)(2Department of Radiotherapy,Shandong Province TumorHospital,Jinan 250117,China)Abstract:Aiming at the distribution characteristic of pixel gray level in tested dat

4、a,this paper proposes avalid and reliable searching and tracking policy for the object contour.Because of the relatively distinct con-trast of gray level between unhealthy tissues and neighboring ones in data,an improved method is used totrack the object contour firstly,then the contour is sampled t

5、o get the discrete controlling points of it and in-put them into the Snake contour tracking algorithm subroutine,thus avoiding the limitation of the initial con-tour controlling points of the Snake method and increasing the speed and accuracy of the segmentation for theunhealthy tissues.The proposed

6、 method is significant for practical application.Key words:Snake model;contour;object tracking;orientation coding;interactive segmentation 收稿日期:2002-08-27.作者简介:王立功(1970),男,博士生,bastion ;罗立民(联系人),男,博士,教授,博士生导师,luo.list .医学图像方面的许多应用都依赖于对图像中目标轮廓的准确跟踪,而目标轮廓可以通过很多途径得到.在对图像的研究和应用中,人们往往仅对图像中的某些部分感兴趣.这些部分常称为

7、目标或前景,其余部分则称为背景,它们一般对应于图像中特定的、具有独特性质的区域.为了辨识和分析目标,需要将它们分离出来,提取出感兴趣区域,在此基础上才有可能对研究对象做进一步的分析和处理.图像分割是由图像处理到图像分析的关键步骤,也是一种基本的计算机视觉技术,这是因为图像的分割、目标的分离、特征的提取和参数的测量是将原始图像转化为更抽象、更紧凑的形式,使得更高层次的分析和理解成为可能的关键步骤.图像分割的主要任务就是要把图像分割成各具特性的区域并提取出感兴趣目标,而区域的特性可以是像素的灰度、颜色和纹理等特征,预先定义的目标可以对应单个区域,也可以对应多个区域.由于图像分割几乎是图像处理和分析

8、以及理解应用领域的第一步,因而多年来对该领域的研究成果也颇为丰富,许多学者针对各自特定研究对象提出了不少有价值的新方法和新观点1 4.然而就大多数方法而言,真正能投入到临床应用的并不多,很多方法都只是停留在实验室研究阶段.1 问题的提出实现图像分割的完全自动化几乎可以被认为是该研究领域的最高目标,然而在大多数场合下全自动分割需要人工干预才能得到较为理想的分割效果.目前临床上绝大多数图像分割工作仍然是由人工来完成,由经验丰富的临床医生在原始胶片或计算机屏幕显示的图像上直接勾画出有关组织的边界或感兴趣区域,尽管到目前为止,人工分割的精度相对较高,但此法费时、费力,需要专业知识,且对操作者的勾画技能

9、有较高要求,如欲分割一个完整的图像数据集需要花费临床医生很长时间才能完成,同时,分割结果也因人而异,可再现性较差,在已实现商品化的大多数应用中均采用此法.本文试图对人工分割方法做一些改进,以提高分割速度和准确性.在放疗计划系统应用中,精确提取出肿瘤组织的边界是整个系统得以顺利进行的最根本保证,若勾画出的病变组织边界多于实际边界则会在随后的放疗中伤及周围的健康组织,对病人术后的生命质量产生很大影响.反之,若少于实际边界则会造成照射不彻底,同样会造成病人术后可能的疾病复发.基于此,采用M.Kass 提出的主动轮廓线模型/蛇0方法5对原始病变组织轮廓进行跟踪,以提高放疗计划系统的自动化程度.Snak

10、e方法及其改进方法在静态图像处理方面主要用于检测图像边缘、曲率、边界以及区域分割和骨架提取等,该法比较适合于复杂的生物医学组织结构的分割.基于 Snake 的算法可分为两类:一类是变分方法;另一类是离散网格法.变分法寻找全局最优解,要求对目标的先验知识知道得足够多,计算量很大,时间开销也很大.离散网格法存在如下一些问题,限制了其实际应用:为使算法收敛,需酌情调整收敛值,否则可能导致不收敛;代表轮廓线的点数固定不变,不能随目标大小作相应改变,而这一点在实际应用中几乎不可行.文献 6 提出了一种变型的 Snake 方法,但其实际上同样也存在着上述的一些问题.M.Kass的 Snake 模型用参数来

11、表示轮廓线:M(s)=(x(s),y(s),定义轮廓线的能量函数 E*Snake=Q10ESnake(M(s)ds=Q10 Eint(M(s)+Eimg(M(s)+Eext(M(s)ds(1)式中,Eint为主动轮廓线的内部能量;Eimg为图像自身作用力产生的能量;Eext为外部限制力产生的能量,由此可见 Snake 的行为是由外部力和内部力来控制的,内部力起平滑性限制的作用,外部力引导初始轮廓线向图像特征逼近,初始轮廓线在内、外力的作用下寻找其能量函数的极小值.式(1)中的Eimg可进一步描述为Eimg=ElEl+EedgEedg+EterEter(2)式(2)是直线、边缘以及界线能量项的加

12、权线性组合,可分别从图像 I(x,y)计算出.El=I(x,y);Eedg=-I(x,y)2,Eter为用高斯函数平滑过的图像中各级轮廓线的曲率;Eint由一阶项和二阶项组成Eint(s)=(A(s)M c(s)2+B(s)M d(s)2)2(3)式中,一阶项当轮廓线不封闭时有较大值,二阶项当轮廓线的曲率变大时有较大值.某个点的 A值和B值决定了轮廓线在该点的延伸和弯曲程度,如A=0,则该点不连续;B=0,则该点为角点.从求取能量函数极小值的推导过程可以看出,M.Kass及其改进的算法(变分法)均要求外力必须可微,这一点使算法的应用受到限制.Snake 方法针对不同作用对象及相应系数取不同的值

13、,作者在实验中取 A=110,B为 0 或 110(取决于该点是否为角点),采用直接来自于临床的 T1 加权 MRI 图像数据,数据层内空间分辨率 256 256,层厚 5 mm,层间距 6 mm,挑选了含有病变组织的 5 幅空间序列MRI 图像,数据灰度对比相对明显,病变组织清晰可见.采用 HPP450 微机,平均跟踪时间约0135 s.图 1 中第 1 行数据为采集自临床实际数据,图 1 采用 Snake 方法跟踪肿瘤组织的运行结果第 2行为通过人机交互产生的初始轮廓线,为对比实验,轮廓线基本上没有封闭,第 3 行为采用 Snake方法跟踪出的病变组织的轮廓线.由跟踪过程可知,/蛇0能很快

14、滑向病变组织的边缘,然而,从轮廓216东南大学学报(自然科学版)第33 卷线包络结果来看并不令人满意,跟踪出的轮廓线能够完全贴住病变组织边界的出现概率并不高,有时需要反复设定轮廓线的初始控制点并改变预先设定的 A和B值才能取得比较理想的结果,而大部分结果则是出现轮廓线的局部漂移,即跟踪出的轮廓线往往会将病变组织周围的正常健康组织也纳入其中,这种情形对于临床应用来说是不可行的.2 改进的方法作者认为,上述问题的出现,根本原因是 Snake方法对初始轮廓线控制点的设定条件及 A和 B值的取值不具普适性.为此,提出如下改进方法:首先对图像进行滤波处理以提高信噪比,尽可能地去除噪声干扰,然后用阈值法对

15、图像进行二值化处理,阈值法对图像进行分割的关键是选取合适的阈值,采用熵值法7求取阈值,其选取准则基于使背景和目标的灰度分布的有效信息量最大,一幅图像的有效信息量用图像的熵来衡量,调节阈值 S,使得 S 处的熵函数 U(S)取得最大值,此时就取得了最佳图像分割阈值.其算法为,计算灰度直方图,并定义熵H(A)=-ESi=1PiPSlnPiPSH(B)=-ESi=1Pi1-PSlnPi1-PS式中,Pi=fiN,PS=ESi=1Pi,fi为图像中灰度为i 的点数;N 为图像的总点数.定义图像的熵函数为:U(S)=H(A)+H(B),求得熵函数 U(S)最大时对应的值 S,此值即为最佳分割阈值T.图2

16、为任选含病变组织的一幅图像所对应的直方图,以及采用熵函数法进行二值化后的结果.由图示结果可见,图 2 采用熵阈值选取法进行二值化的结果采用熵阈值法分割图像能够比较理想地保留病变组织区域,并且该组织边界与周围组织间的界限对比明显,组织区域大小与其实际大小从分割结果来看大致相当,尽管如此,仍然无法肯定分割出的组织边界恰处于实际病变组织的边界上.基于此,首先采用一种适应性目标轮廓跟踪算法,该算法源于边界方向链码轮廓跟踪算法的思想.链码是对边界点的一种编码表示方法,其特点是利用一系列具有特定长度和方向的相连的直线段来表示目标的边界.因为每个线段的长度固定而方向数目取为有限,所以只有边界的起始点需用边界点的绝对坐标表示,其余点都可只用接续方向来代表偏移量.常用的有 4和 8 方向链码,分别表示为:上、下、左、右和上、下、左、右、左上、左下、右上、右下,如图 3 所示.这样,对于图像中任一目标轮廓只要用处于轮廓线上任一点的坐标和相应的方向编码即可表示,很明显,此轮廓线表示方法具有目标平移和尺度变换不变性,但不具目标旋转的不变性.基于类似思想,采用一种可称之为/摸索法0 的轮廓跟踪方法来跟踪目标轮廓