[发明专利]不均匀网格分布模拟直线加速器治疗计划的剂量计算方法

权利要求书

1.一种不均匀网格分布模拟直线加速器治疗计划剂量计算方法，适于在计算设备中执行，其特征在于：包括如下步骤：

(1)设置不均匀网格参数：

(i)设置剂量计算的二维通量网格的边界；

(ii)根据多叶准直器叶片的厚度初步分割二维通量网格y方向的初始网格间距；

(iii)根据多叶准直器片间耦合位置二次分割二维通量网格y方向的网格间距，得到y方向的二次分割网格间距；

(iv)设置二维通量网格x方向的网格间距；

(2)根据片间漏射透射测量结果设定对应网格的初始通量值；

(3)根据不均匀网格参数和各网格对应的初始通量值调用剂量计算方法计算得到每个网格对总剂量分布的贡献，最后将各网格的剂量计算结果叠加，得到整个射野在选定模体中的剂量分布；

(4)测量选定模体在该射野下的真实剂量；

(5)对比真实剂量与剂量计算的整体差异，判断真实剂量与剂量计算的差异是否在预设的计算精度范围内；

若差异在预设的计算精度范围内，接受当前网格分割参数和通量值，用于病人治疗计划的剂量计算；

若差异超过预设的计算精度范围时，重复步骤(6)～(7)；直到剂量计算结果与真实剂量的差异满足预设精度；

(6)进一步调整二维通量网格y方向的二次分割网格间距，得到y方向的优化网格间距；以及进一步调整各网格对应的初始通量值，得到优化通量值；

(7)根据优化后的网格参数和各网格对应的优化通量值，调用剂量计算方法计算得到每个网格对总剂量分布的贡献，最后将各网格的剂量计算结果叠加，得到整个射野在选定模体中的剂量。

2.根据权利要求1所述的不均匀网格分布模拟直线加速器治疗计划剂量计算方法，其特征在于：步骤(i)中，剂量计算的二维通量网格的边界由每个治疗计划中钨门运动形成的最大外轮廓确定。

3.根据权利要求1所述的不均匀网格分布模拟直线加速器治疗计划剂量计算方法，其特征在于：步骤(ii)中，初步分割的二维通量网格y方向的网格间距与和其对应的多叶准直器叶片的厚度一致。

4.根据权利要求1所述的不均匀网格分布模拟直线加速器治疗计划剂量计算方法，其特征在于：步骤(iii)中，根据多叶准直器中突出和刻槽的位置对二维通量网格y方向的网格间距进行二次分割。

5.根据权利要求1所述的不均匀网格分布模拟直线加速器治疗计划剂量计算方法，其特征在于：步骤(iv)中，二维通量网格x方向的网格分割为等间距分割。

6.根据权利要求1所述的不均匀网格分布模拟直线加速器治疗计划剂量计算方法，其特征在于：步骤(iv)中，二维通量网格x方向相邻网格的间距小于一对叶片闭合时的间距值。

7.根据权利要求1所述的不均匀网格分布模拟直线加速器治疗计划剂量计算方法，其特征在于：步骤(3)中，所述的剂量计算方法为蒙特卡罗剂量计算方法、笔形束计算模型、NeuralRad剂量计算方法或RBM计算方法。

8.一种计算设备，包括：

一个或多个处理器；

存储器；以及

一个或多个程序，其中所述一个或多个程序存储在所述存储器中并被配置为由一个或多个处理器执行，所述的一个或多个程序包括用于实现上述如权利要求1-7中任一项所述的不均匀网格分布模拟直线加速器治疗计划剂量计算方法的指令。

9.一种存储一个或多个程序的计算机可读存储介质，所述的一个或多个程序包括指令，所述指令适于由存储器加载并执行上述权利要求1-7中任一所述的不均匀网格分布模拟直线加速器治疗计划剂量计算方法。