[发明专利]一种伽玛射线剂量的卷积计算方法

说明书

技术领域

本发明属于放射治疗计划剂量并行计算方法技术领域，涉及一种伽玛射线剂量的卷积计算方法。

背景技术

伽玛射线立体定向回转聚焦放疗机是一种融立体定向技术和放射外科技术于一体的立体定向放射外科治疗设备，它采用伽玛射线几何聚焦方式，通过精确的立体定向将经过规划的大剂量伽玛射线集中照射于体内的预选靶区，一次性、致死性地摧毁靶点内的组织，以达到外科手术切除或损毁的治疗效果，立体定向放射外科与普通外科有本质的不同，它避免了传统外科开放式手术所带来的术后出血、感染以及可能损伤重要功能器官的危险，创建了一种无创伤的手术方法。

目前，伽玛射线立体定向回转聚焦放疗机治疗计划的剂量计算，采用机器出厂时，测量各个尺寸准直器在单个伽玛放射源下在水中沉积剂量的OAR(离轴比)和TMR(深度剂量比)，作为剂量计算的基础数据；治疗计划的剂量分布是由受照区空间内所有单个体数据离散点的剂量组成，计算单个体数据点的剂量时，首先根据单个伽玛放射源的焦点剂量和放射源的空间分布以及OAR和TMR数据计算单个伽玛放射源在此点的剂量贡献，然后对所有放射源的剂量贡献进行数值积分，再结合该点的电子密度计算沉积剂量，结果即为该点的综合沉积剂量；最后遍历体数据的所有离散点，这样就完成了单个靶点的剂量计算，多个靶点的剂量由单个靶点的剂量叠加而成。

然而，每台机器出厂时，准直器尺寸因为工艺原因都有差异，而这个差异导致了OAR和TMR的差异，OAR和TMR的不准确会直接影响剂量计算的精度。所以，就每台机器出厂时都要进行单源测试。

发明内容

本发明的目的是提供一种伽玛射线剂量的卷积计算方法，解决了现有技术中每台机器出厂时都要进行单源测试的问题。

本发明所采用的技术方案是，一种伽玛射线剂量的卷积计算方法，具体按照以下步骤实施：

步骤1，采用蒙特卡罗方法模拟生成伽玛放射源点源在水中的剂量沉积分布K(x,y,z)；

步骤2，依据准直器编号对应的机械尺寸，采用光线追踪法计算由准直器形成的单个伽玛放射源在水中的原射线沉积分布T_n(x,y,z)，n为准直器的编号，共有N个准直器，n＝1,2,3,...,N，(x,y,z)为坐标点；

步骤3，将步骤1生成的放射源点源在水中的剂量沉积分布K(x,y,z)和步骤2中形成的单个伽玛放射源在水中的原射线沉积分布T_n(x,y,z)进行卷积计算，生成单个伽玛放射源在水中的剂量沉积分布：

G_n(x,y,z)＝T_n(x,y,z)*K(x,y,z)＝∫∫∫T_n(x,y,z)·K(X-x,Y-y,Z-z)dxdydz，

其中，n为准直器的编号；

步骤4，给伽玛射线立体定向回转聚焦放疗机的放疗计划系统中导入病人信息和影像体数据，导入或勾画皮肤轮廓、危及器官轮廓和靶区；

步骤5，根据步骤4导入的信息和数据，设计和调整靶点位置、权重、准直器编号、拉弧角度，制定放疗布靶方案；

步骤6，对步骤4导入体数据中的每个离散点采用GPU并行计算，计算单个靶点在体数据的吸收剂量分布：D_T＝(x,y,z；m)，m为靶点的编号，m＝1,2,3,...,M；

步骤7，重复步骤6，对所有靶点在体数据的吸收剂量分布进行叠加，得到靶区的吸收剂量分布：式中，m为靶点的编号，M为靶点的个数；

步骤8，对经步骤7得到靶区的吸收剂量分布和治疗计划进行评价，如果不满足临床要求，则转入步骤5；如果满足临床要求，则输出计划，计算结束。

本发明的特点还在于，

步骤6中计算单个靶点在体数据的吸收剂量分布，具体按照以下步骤实施：

步骤6.1，初始化GPU，进入并行计算；

步骤6.2，把治疗头的旋转角度离散化，每隔0.5～5度采样一次，计算各个采样角度处每个放射源的源皮距；