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

三维医学图像配准在图像引导放疗中的应用 　　摘要：为了确定病人的摆位误差，实现精确放疗，提出一种改进的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]。例如Wang等[9]提出了一种改进算法Active Demons 算法，将浮动图像的梯度信息也加入到光流场方程中，但作者并没有对双向力的分配问题进行更深入研究；Vercauteren等[10]提出了一种基于灰度的配准模型，将图像灰度的均方差作为相似性测度，没有考虑图像灰度的梯度场对配准过程的影响；周付根等[11]将图像间的对称梯度和待配准图像的梯度引入到Demons算法中，生成一种合力并计算图像形变力，但两个单向力仅是简单的叠加，没有分析两个作用力相对重要性；目前弹性配准理论还有待进一步发展，真正意义上的应用仍不多，尤其针对图像引导放疗领域的医学图像3D/3D弹性配准研究较少，其计算效率和稳定性也有待提高。 　　FDS团队自主研发了精确放射治疗系统（Advanced/Accurate Radiation Therapy System， ARTS）[15-22]，并针对图像引导放射治疗，在ARTS基础上进一步开发了IGRT原型系统（ARTSIGRT）[23-25]。本文的研究工作是ARTSIGRT 1.0的核心部分。首先在ITK基础上实现了基于Demons的弹性配准算法，结果显示该配准算法鲁棒性强、误差小；其次，分析了Demons配准算法的不足，在此基础上引入对称梯度概念，提出一种带权重的对称Demons（Weighted Symmetric Demons， WSD）配准算法，以希望应用于IGRT中的多模态图像配准。最后通过单模态CT/CT图像配准和多模态CT/CBCT图像配准实验，分析WSD配准算法的精度，结果表明该算法可以提高配准精度，加快收敛速度。 　　2实验与结果 　　本文利用WSD算法分别配准不同模态的三维临床医学图像，并与原始Demons算法的配准结果进行比较，从而验证改进算法的有效性。下面通过两组实验验证WSD算法的精度和效率。 　　第