三维医学图像配准在图像引导放疗中应用.docxVIP

三维医学图像配准在图像引导放疗中应用 摘要：为了确定病人的摆位误差，实现精确放疗, 提出一种改进的Demons弹性配准算法。采用FDK算法对锥 形束CT(CBCT)图像进行三维重建，利用可视化工具包(VTK) 体绘制法可视化重建结果；在分割与配准工具包(ITK)基 础上实现Demons算法，并基于对称梯度的思想，将参考图 像和浮动图像的梯度场信息加入到Demons算法中，给出新 的Demons形变力公式。分别使用单模态和多模态医学图像 进行配准实验，结果显示改进的Demons算法与原始Demons 算法相比，配准速度更快、精度更髙。基于对称梯度的Demons 算法更适用于图像引导放射治疗中CBCT重建图像与CT计划 图像间的配准。 关键词：图像引导放射治疗；图像配准；Demons算法; 医学影像；三维重建；对称梯度 中图分类号：TP391 文献标志码：A 0引言 图像引导放射治疗(ImageGuided Radiotherapy, IGRT) 将放射治疗机与成像设备结合在一起，在患者治疗前、治疗 中利用各种先进的影像设备对肿瘤及正常器官进行实时的 监控，釆集图像信息，确定治疗靶区和重要组织在治疗过程 中的运动和分次治疗间的位移误差，并在必要时进行位置和 剂量分布的校正，减少放疗间靶区位移误差和摆位误差，监 测和校正放疗时肿瘤和正常组织运动引起的误差，已成为实 现精确放射治疗的有效方法。 医学图像配准可以匹配治疗前和计划时的单模或多模 图像数据，根据病人体位位移调整病人摆位，是实现IGRT 平台的关键技术。3D/3D配准保证图像在三个维度位置的一 致性，是图像配准临床应用的关键。按照空间变换，3D/3D 图像配准分为基于刚性变换(Rigid Transformation)的 3D/3D配准和基于弹性变换(Elastic Transformation)的 3D/3D配准。并且随着计算机技术的发展，国内外学者对 3D/3D图像配准的研究也越来越多。例如Jung等［1］提出基 于互信息(Mutual Information)的三维配准并采用GPU加 速配准速度，Pennec等［2］基于光流场模型实现了三维图像 配准，Nam等［3］基于B样条的弹性变化实现MR和CT 3D多 模图像间的配准等。刚性配准方法不适合处理大形变配准 ［7］。但是在放疗过程中，由于病人的呼吸、脏器的蠕动和 体位的变动，软组织器官难免产生较大形变。因此，本文重 点研究IGRT系统中的3D/3D图像弹性配准。 Demons算法［8］是一种基于图像灰度的弹性配准方 法。该算法使用参考图像的灰度梯度信息来决定浮动图像每 个像素的移动，具有鲁棒性强、误差小的特点；但仅靠梯度 信息往往是不充足的，特别是当梯度非常小的时候，图像变 形的方向不能确定，容易导致错误的配准变换，因此一些学 者针对此问题对Demons算法进行了改进［9-12］ o例如Wang 等［9］提出了 一种改进算法Active Demons算法，将浮动图 像的梯度信息也加入到光流场方程中，但作者并没有对双向 力的分配问题进行更深入研究；Vercauteren等［10］提出了 一种基于灰度的配准模型，将图像灰度的均方差作为相似性 测度，没有考虑图像灰度的梯度场对配准过程的影响；周付 根等［11］将图像间的对称梯度和待配准图像的梯度引入到 Demons算法中，生成一种合力并计算图像形变力，但两个单 向力仅是简单的叠加，没有分析两个作用力相对重要性；目 前弹性配准理论还有待进一步发展，真正意义上的应用仍不 多，尤其针对图像引导放疗领域的医学图像3D/3D弹性配准 研究较少，其计算效率和稳定性也有待提高。 FDS团队自主研发了精确放射治疗系统 (Advanced/Accurate Radiation Therapy System, ARTS) ［15-22］,并针对图像引导放射治疗，在ARTS基础上进一步 开发了 IGRT原型系统(ARTSIGRT) ［23-25］ o本文的研究工 作是ARTSIGRT 1.0的核心部分。首先在ITK基础上实现了 基于Demons的弹性配准算法，结果显示该配准算法鲁棒性 强、误差小；其次，分析了 Demons配准算法的不足，在此 基础上引入对称梯度概念，提出一种带权重的对称Demons (Weighted Symmetric Demons, WSD)配准算法，以希望 应用于IGRT中的多模态图像配准。最后通过单模态CT/CT 图像配准和多模态CT/CBCT图像配准实验，分析WSD配准算 法的精度，结果表明该算法可以提高配准精度，加快收敛速 度。 2实验与结果 本文利用WSD算法分别配准不同模态的三维临床医学图 像，并与原始Demons算法的配准结果进行比较，从而验证 改进算法的有效