伽玛刀TPS中放射线剂量分布的三维可视化方法.docxVIP

医学 TPS 中放射线剂量分布的三维可视化方法 摘 要: 针对放射治疗计划系统 GyroRTPS 需要 ,给出了 一个放射线剂量分布三维可视化的实现方法 , 以帮助医师多 角度查看病灶组织任意剖面上的放射线剂量分布情况。首先 论述系统中对体数据和剂量数据的坐标变换、插值计算等预 处理 ;然后根据 VTK 工具包中可处理的数据类型对数据进行 转换 ,给出了基于 VTK 实现剂量分布三维显示及其与体数据 融合显示的方法。 实例测试结果和临床使用表明 ,基于该方法 开发的软件系统运行可靠 , 可更好地辅助医师进行放射治疗 决策。 关键词 : 医学治疗计划系统 ;放射线 ;剂量分布 ; 可视化方 法 0 引言 治疗计划系统 (Treatment Planning System, TPS) 是 伽玛刀治疗的配套软件。在临床放射治疗前 , 通过对满足 DICOM(Digital Imaging and Communications inMedicine) 的医学序列 CT ( Computed Tomography) 或 MRI (Magnetic ResonanceImaging) 等图像的分析 , 实现 对病灶区域的精确定位 , 同时由治疗计划系统模拟射线照射 时的剂量分布情况 [1-2], 辅助医生制订准确的放疗方案 ,使放 射治疗最大限度杀死癌变细胞的同时 , 避免或减少对正常组 织和重要器官的伤害。 简单的二维 TPS 依据某一方向的 CT 切片数据进行模拟 放射 ,只能显示单张切片上的剂量分布情况 ,不能全面地反映 病灶区的受剂量情况 , 无法模拟真实的立体放射治疗。 一个完善的放射治疗 TPS[3], 需要精确描述病灶区域和 正常组织的空间关系 , 建立准确的剂量分布数据场 [4], 最终能 够对人体患病组织的剂量空间分布进行多角度多层面的可 视化显示处理 ,帮助医生制订准确详尽的放射治疗方案 , 以期 达到提高放射治疗精度、减少放射损伤的目的。对剂量分布 采用三维立体方式的可视化 ,将成为伽玛刀放射治疗中一个 不可缺少的功能。 1 二维剂量分布显示 剂量通常是指一次给药后产生药物治疗效果的数量 , 这 里则用于表征射线强度。 在伽玛刀放射治疗的配套软件 TPS 中,剂量分布是指伽玛射线的粒子照射到人体中能量吸收的 空间分布。剂量计算一般由 TPS 中的单独模块完成 ,剂量计 算完成后 , 空间某点的剂量值存放在一个三维数组中供显示 模块调用。传统的剂量显示按照等值线绘制方式得到平面等 剂量线图 ,如图 1 所示;或者显示空间各层的等剂量云图 ,如图 2 所示。 图1等剂倉找屋示效采 图2分层等剂趾祸显彌果 虽然平面显示可方便地查看某一切片上的剂量分布情 况,但是无法提供一个三维 （立体）的剂量分布显示[5],医生也 无法从整体上查看受照射病灶组织中剂量分布情况 ，所以仍 然需要医生依靠经验，由多幅二维图像估计病灶组织的大小 及形状，“构思”剂量分布与病灶组织的三维几何关系。 2体数据和剂量数据的预处理 随着计算机图形学和科学计算可视化技术的发展 ，可以 很方便地利用一些可视化工具包 （例如VTK等）对体数据进行 三维重建，并对绘制结果提供任意面剖切显示 [6]。本文依据 作者所在项目组从事研发的“陀螺刀”放射治疗计划系统 GyroRTPS（Gyro Radiotherapy Treatment Planning System）来探索利用VTK进行剂量分布的三维可视化研究。 本文中涉及的体数据为符合 DICOM协议的CT或MRI图像 数据,剂量则依据“陀螺刀”设备的制造原理由 GyroRTPS 计算得到。 2. 1 坐标转换 在GyroRTPS中，由剂量计算模块计算得到的剂量数据 存放在一个float型的三维数组中，其计算是在如图3所示的 坐标系中进行。为了使处理的数据量尽可能小，并不是对每张 图片的所有像素都做处理，只计算兴趣图像区内的剂量分布 ’ 兴趣图像区的大小和位置由医师勾勒得出。 所以，剂量数据的 分布原点也不是坐标系的原点。 利就分布原点 圏3刑呈坐标椒休数据坐标木盘:阁 如图3所示,剂量分布原点所在的坐标系为计算剂量数 据时所用的坐标系，其外层的坐标系则为原始图像所在坐标 系,两者所用并不是一个坐标系。 而在对DICOM 图像序列进 行三维重建时的坐标系（在此只考虑XY轴方向的二维坐标系） 则是以原始图像坐标系的点（0,N）为原点（N为原始DICOM 图像的宽度 ), 以向右为 X 轴正方向 ,以向上为 Y 轴正方向 ,X 轴与原始 DICOM 图像下边重合的坐标系。 所以在对原始剂 量数据进行三维重建之前 ,必须进行坐标转换 , 将剂量坐标系 下的剂量数据转换为三维重建坐标系下的数据。 设原始图像所在坐标系为