基于ITK的数字重建放射影像重建算法与应用.pdfVIP

基于ITK的数字重建放射影像重建算法与应用 闫锋1，2 罗月童1，2龙鹏程2李佳2黄善清2孔维华1，2 曹瑞芬2 李贵2胡丽琴2吴宜灿2，3FDS团队 1合肥工业大学计算机与信息学院VCC研究室，合肥230009 2中国科学院等离子体物理研究所，安徽合肥230031 3中国科学技术大学核科学技术学院，安徽合肥230027 Reconstlllcted 摘要：数字重建放射影像(Digitally 中有广泛应用，有多种实现方法。ITK是开源的图像分析和处理的开发包，在医学图像处理等领域有 着广泛的应用。本文在讨论了DRR图像的基本原理的基础上，介绍了一种基于ITK的数字重建影像成 像算法，并将该算法应用于自主开发的精确放疗系统，效果良好。 关键词：数字重建放射影像ITK放疗计划系统 1引 言 Recons仃ucted 数字重建放射影像(Digital 来越多的取代传统CT模拟定位机所用的胶片图像。DI汛广泛应用于CT模拟定位、图像引导放射治疗 的叠加情况，还能用于肿瘤靶区的定义、验证和确认，以此确保病人在实际治疗时的摆位与在治疗计划中 的摆位相同。 目前，DRR重建算法主要采用光线投射法实现。刘鹏，高军等人通过对光场重现理论进行延伸，用光 踪算法14J，无需解方程求交点，只包含加减运算，计算速度快，缩减了成像时间。汪俊，吴章文等采用基 于射线追踪的方法，对密度进行指数变换，通过调整衰减系数和密度指数来实现高密度组织与低密度组织 间对比度的提高pJ，节省了内存空间。陈为采用基于纹理映射的光线投射算法来加速生成DRR图像161。 ITK是一个功能非常强大、免费开源的开发包17J，实现了大量的图像处理算法，在医学图像处理领域 有着广泛应用。目前国内关于ITK在数字重建影像方面的应用很少。因为ITK中的算法都经过了大量测试， 鲁棒性强，所以基于ITK开发不仅能节约大量时间，而且能保证算法的鲁棒性，所以本文讨论了基于ITK 工具包实现了数字重建影像的成像的方法。 2开发．环境 and ITK(Insi曲tSegmentationregistration 开发包系统。IK以免费开源的形式开发，形成了一个不断增长的、可以自己维护与开发的用户社区及开发 人员社区，使ITK能够不断的进步。ITK主要用于图像分析和处理，实现了当前大部分图像分割和配准算 法，且对图像的预处理有很好的效果，但是缺乏对可视化方面的支持，需要与其他可视化工具配合使用。 Filter中相关算法处理读入的图像数据，最后产生新的lIIlageFile。 3Dl汛算法实现 DI汛图像的生成过程是x光学衰减的过程，这个过程符合光学吸收模型规律，因此DRR中x光线的 心 衰减过程可用表达式，(s)=厶eXp(一【f(f)出)描述。其中： W ◆ s为光学投射方向的长度参数； ◆J(s)为距离s处的光学强度； ◆ f(f)是光线强度的衰减系数，即沿光线投射方向f处的光线吸收率； ◆ 厶是光线进入三维数据场时的光线强度。 DRR生成是一个模拟X光线穿透CT体元，经过衰减和吸收后投影到成像平面后累加的过程。DRR 生成的过程如图1所示，基本可以分为五个步骤： (1)建立一个CT图像组的三维体素矩阵。 (2)沿虚源向CT图像组发出若干条投影线，射线数量与DRR平面的像素个数相同。 (3)获得每条投影线通过此体素矩阵的交点，将这些点的电子密度值累加。 (4)求投影线通过体素矩阵的有效射线长度。 (5)将电子密度累加值与射线长度相乘，求出的值按灰度显示，即为DRR图像。 图lDRR图像生成过程示意图 在ITK中，可以利用其中的一些函数来实现Dl讯算法设计的五个步骤，最终实现DRR算法。实现过 程如下： (1)通过ITK中的Im