Demons算法在四维CT图像配准中的应用研究.pdf

中国体视学学会cT理论与应用分会、中国工程物理研究院 2008年7月 应用电子学研究所、(CT理论与应用研究》编辑部论文集 Demons算法在四维CT图像配准中的应用1 段卓镭，周付根 (北京航空航天大学图像中心，北京100083) 摘要：由呼吸运动引起的器官的形变是4D精确放射治疗中最主要的误差来源。利用4D图像配准 技术来分析不同呼吸时相中器官的形变是目前研究的热点，而研究的方法是基于形变模型的3D配 准算法。基于光流场的Demons算法具有速度快和精度高的特点．本文将其应用到4DcT图像配准 中，并在原始图像形变力的计算中引入了对称梯度。实验结果表明．结合对称梯度的Demons方法 具有更快的收敛速度，而且匹配误差小． 关键词：Demons算法；对称梯度；4DCT图像配准 1 引言 在放射治疗过程中，由于呼吸运动的存在，肺部、胰腺、肝脏和其他胸腹部的肿瘤， 均可随着呼吸运动产生位移和形变，从而造成了治疗计划中的剂量分布和目标实际接受的 剂量分布之间存在差异。 传统的3DCT影像只采集了呼吸周期中某一个时刻的信息，不能准确代表患者接受治 疗时器官在呼吸周期中的状态。随着4DCT在临床上的应用，放射治疗中补偿目标位置随 时间的变化成为可能。4DcT提供了一个呼吸周期内不同时相的数据，目前已经有机构开展 了4DcT放射治疗计划的研究，例如u-21。制定4D放射治疗计划最关键的技术是4DcT图 像配准。把不同呼吸时相的CT图像都相对于一个参考时相图像进行3D形变图像配准。通 过配准后得到不同时相之间的空间变换参数，可以准确有效地计算和分析器官在呼吸周期 内的位移和形变情况。 图像配准实质上是在3D空间寻找一个最佳变换．可以视为一个优化问题嘲，通过对空 间变换参数的调整达到三维体素之间的一一对应。三维体素形变配准算法可以归纳为基于 空间变换的方法和基于物理模型的方法。基于空间变换的形变配准主要有基于多项式、基 函数和样条函数的方法。而基于物理模型的形变方法认为图像之间的差异是由某种物理形 变引起的，其思想是构造能够拟合这种形变的物理模型，目前主要的物理模型有弹性模型、 粘性流体模型和光流场模型H。。 CT图像具有灰度上的连续性，故采用基于灰度的配准方法比较合适。而基于光流模型 的Demons配准算法正是一种基于体素灰度的全自动、稳定性强的匹配方法。文献n们中通过 对比三种典型的形变配准算法显示Demons算法更适合于临床应用。因此，本文采用基于光 流场的Demons配准算法，并将对称梯度引入到原始算法中图像形变力的计算。实验结果表 明，结合对称梯度的Demons方法具有更快的收敛速度，而且匹配误差小。 基金项目：广东省科技计划项目“四维精确放射治疗计划系统”(200723一E0431)． 段卓镭等：Demons算法在四维cT图像配准中的应用 33 2原始Demons算法及其改进算法 2．1 Thirion’s Demons配准算法 的方法。他在图像中引入了“Demons”力，即图像的形变力。每一个体素上的“Demons” 力驱动待配准图像向参考图像方向变形，以达到两幅图像间的匹配。参考图像的梯度Vf是 驱动“Demons”的内力(internal force)‘引。 力(external 参考图像F上的灰度，m是待配准图像M上的灰度。P点在两幅图像上的偏移可以由下式得 出。 ；：一 !Z二12 (1) ：×可f M+与譬。 吖是参考图像的梯度，它代表了参考图像上相邻点之间的联系，所以把它称为来自参 考图像的内力(internal force)。U一所)是两幅图像对应像素点的灰度差，被称为形变 的外力(external force)。为了在全局范围内使该变换连续，通常的做法是在每一个迭代 后，使用高斯滤波来平滑所得到的偏移。此算法是一个循环迭代过程，