分段迭代B-snake模型.docVIP

分段迭代B-snake模型 第22卷第5期 2009年l0月 四川理工学院(自然科学版) JournalofSichuanUniversityofScienceamp;Engineering(NaturalScienceEdition) Vo1.22No.5 0ct.2009 文章编号:1673-1549(2009)05-0096-05 分段迭代B—snake模型 王文哲一,唐克伦,牟宗魁一,郑彬 (1.~1)ll理工学院机械工程学院,四川自贡643000; 2.四川省高校过程装备与控制工程重点实验室,四川自贡643000) 摘要:传统B—snake模型以整体迭代为基础,已经逼近到目标的轮廓曲线仍要参加迭代计算,影 响轮廓曲线的收敛速度.提出了一种改进的B—snake模型——分段迭代B—snake模型,它以分段迭代 为基础,在迭代过程中只计算尚未逼近的轮廓曲线,直至所有的轮廓曲线均逼近目标轮廓.实验证明, 与传统B—snake模型相比,分段迭代B—snake模型可以快速收敛到目标轮廓. 关键词:分段迭代;B—snake;轮廓提取;整体迭代 中图分类号:TP391文献标识码:A 引言 主动轮廓模型,又称为snake模型,是Kass等人… 在1987年提出来的,snake模型对噪声和对比度不敏 感,能将目标从复杂背景中分割出来,并能有效地跟踪 目标的形变和非刚体的复杂运动而被广泛用于如边缘 提取,图像分割和分类,物体跟踪,3D重建,立体视觉匹 配等图像处理领域. 在当前的CAD/CAM系统中,B样条(B—spline)曲 线曲面拟合技术,已经成为几何造型的核心部分. Menet将B样条曲线与主动轮廓模型结合,提出了B 样条主动轮廓模型,也称为B—snake模型;该模型以B 样条曲线的控制点为未知量,将主动轮廓表示为B样条 曲线的形式,然后采用有限差分法对欧拉方程进行求 解,得到使模型能量极小化的控制点坐标.鲁爱东等 提出了一种用户交互操作与B—snake相结合的半自动 分割方法,在超声图像中提取肝瘤轮廓,并成功地应用 于肝肿瘤手术仿真系统.Cheng等人成功采用自适 应有限元法求解snake模型,并将其用于MRI图像的三 维重构. 在传统的B—snake模型中,主动轮廓模型迭代是以 整体迭代的方式进行,即所有控制点均参与每一次迭代 计算过程,已经逼近到目标轮廓的曲线依然要参加进一 步的迭代计算,这无形之中增加了迭代计算的工作量, 影响了曲线收敛速度.在传统B—snake模型的基础上, 本文提出了一种改进的B—snake模型——分段迭代B — snake模型.分段迭代B—snake模型的基本原理是一 次迭代之后通过检测找出已经逼近到目标的曲线单元, 并将其分离出来,在进一步的迭代中只计算尚未逼近目 标的曲线单元;当所有曲线单元均逼近目标之后,对各 个部分曲线进行组合生成整体轮廓曲线.与传统B— snake模型相比,分段迭代B—snake模型能够快速,稳定 收敛到目标轮廓,节约了迭代时间,提高了snake模型的 迭代收敛速度. l传统B—snake模型 Kass等提出的经典snake模型可以表示为定义在 s[0,1]上的参数曲线,即:(s)=((s),Y(s)),其 中s为弧长;以(s)为变量的snake模型能量函数E 可以表示为: 11 Eh:【÷(l()l+I(s)I)+ s((s))(1) 式中参数和操纵着模型的物理行为和局部连续性. snake模型寻找此能量函数局部最小值的过程,是轮廓 曲线由初始位置向真实轮廓逐渐逼近的过程,即通过能 量函数的动态优化来逼近目标的真实轮廓. B样条主动轮廓模型结合了B样条曲线的特性,在 收稿日期:2009-04-17 作者简介:王文哲(1983.),男,湖北天门人,硕士生,主要从事计算机辅助设计方面的研究. 第22卷第5期王文哲等:分段迭代B—snake模型97 实际使用过程中,常采用三次均匀B样条曲线来表示主 动轮廓模型的轮廓曲线一. 由于四个连续的控制点可以确定一段三次均匀B 样条曲线,即有限元计算单元,而参与迭代计算的轮廓 曲线就是由这样的曲线单元顺序连接而成,所以可以将 有限元单元e(s)在几何上表示为形函数Ⅳ与单元节点 向量:的乘积: e(s)=N? = [日..3(S),日(),B:3(),B3.3()] ? [,+,+,+,(2) 其中B(s),12=0,l,2,3为三次B样条的基函数, 为B样条曲线的控制点.结合式(1),进一步从8E=0 可以得到单元e(s)的有限元方程 K?P=F(3) 其中单元刚度矩阵为 ONrON +J8(4) 单元载荷向量 F:fNq(e(s))dsJ0 (5) 按照各单元在整体轮廓中所占位置对单元刚度矩 阵,载荷向量及节点向量进行组集,即可得到模型的