CN108415058A - 放射线的剂量计算方法及系统 - Google Patents

Abstract

本发明适用于医学肿瘤放射治疗技术领域，提供了一种放射线的剂量计算方法及系统。所述方法包括：A、将笔束核分为主核和散射核两部份，分别构建主剂量模型和散射剂量模型；B、利用预先建立的接受放射线照射治疗的模体的模型及预设的放射治疗深度算法，计算所述放射线照射的深度；C、根据所述放射线照射的深度对所述主剂量模型和散射剂量模型分别进行蒙特卡罗仿真模拟，并根据所述蒙特卡罗仿真模拟的结果计算所述笔束核的笔束核值；D、根据所述放射线照射的深度和所述笔束核值计算放射线照射所述放射线照射治疗的模体的剂量。借此，本发明的放射线的剂量计算方法可大大降低计算复杂度，提高剂量计算速度，在保证速度的同时也提高了准确度。

Description

放射线的剂量计算方法及系统

技术领域

本发明涉及医学肿瘤放射治疗技术领域，尤其涉及一种放射线的剂量计算方法及系统。

背景技术

癌症(即恶性肿瘤)是一种严重威胁人类健康的常见多发的重大疾病。70年代以来，癌症发病数以年均3％到5％的速度递增，已成为人类第二大死因。

癌症治疗的目的是治愈癌症患者，延长生命和提高生活质量。放射治疗作为世界卫生组织提出当前癌症治疗的三大手段之一，其使用光子、电子等射线照射肿瘤区域，利用放射生物学效应最大程度地杀伤癌细胞、破坏癌组织并使癌肿瘤缩小，同时尽可能使肿瘤周围的正常组织受到的损伤最小化，在治疗肿瘤的同时以最小的代价换取患者长期高质量的存活。目前70％左右的肿瘤病人在病程的不同阶段需要接受放射治疗。

肿瘤放射治疗的剂量模拟的实现主要有以下三部分：图像信息处理、人体建模和剂量计算。

(1)图像信息处理

首先要对来自医学影像设备的数据，例如CT(Computed Tomography，计算机断层成像)、MRI(Magnetic Resonance Imaging，磁共振成像)、PET(positron emissiontomography PET,正电子发射计算机断层显像)等进行图像的预处理，包括配准、融合、校正、滤波和增强等，然后对皮肤边界、病灶靶区、危及器官等进行分割。

(2)人体建模

将医学影像数据(如CT数据)转化为不同剂量计算方法所需的模型，包含三维几何(解剖结构、位置属性)和物理信息(材料组成、密度信息)，可真实再现人体复杂解剖结构及元素组成，为高精度剂量计算提供精确的模型基础。人体建模模块可以根据剂量计算方法的不同需求建立相应的计算模型。对于解析剂量计算方法，提供人体几何材料分布信息数据。

(3)剂量计算

传统的笔形束算法具有快速和均匀区域精度较高的优势，但在组织非均匀区域存在较大误差；蒙特卡罗方法在所有区域均可保持较高精度，但是极为耗时。传统的笔形束算法使用的剂量沉积核来自于经验的百分深度剂量和百分离轴剂量分布，不能准确描述各个点的沉积核。

剂量计算是强调放射治疗计划系统的核心内容，其计算精度直接影响放射物理师制定治疗计划的质量，以致最终影响放疗的精度和放疗疗效。考虑临床实际需要，计划系统对剂量计算速度及时间有着严格的要求，特别是在逆向计划中，优化功能需要进行大量反复的剂量计算试算，因此对计算速度要求更加苛刻。快速高精度剂量计算方法的研究始终是强调放射治疗计划系统研究的主要内容之一。

解析算法中的笔形束算法因其形象的描述了宽广照射野离散而成的“基本束”特征，而广泛用于临床治疗计划，它的基本思路是将照射束约束设备，如多叶光栅、挡板等形成的不规则野离散成一个个截面无限小的基本束，也就是零野条件下从源点射出的一束无穷小尺寸的射线，每个基本束都会在模体产生剂量分布，这个分布被称为笔束核，人体中任意一点的剂量可以看作为各基本束在此点形成剂量的叠加。

笔形束算法计算简单、所需的计算硬件要求也较低，但对于人体组织的不均匀性和体表不规则轮廓的处理精度不高，如果加上对不同介质的不均匀修正，则可以提高不同介质剂量的准确性。

蒙特卡罗方法对放射线与模体每一个粒子的随机相互作用进行模拟，获得粒子在人体组织中的剂量的分布。当前蒙特卡罗方法及其仿真程序在放射治疗剂量计算及验证方面得到了广泛的应用，同时蒙特卡罗剂量计算方法的致命弱点--收敛速度缓慢、计算时间漫长，使得长期以来一直无法在临床放射治疗上应用。目前它在以下方面起到重要的作用：①用来仿真各种复杂条件下的剂量分布；②作为其他剂量算法精度的验证工具；③为其他解析模型提供剂量沉积核；④分离原射线和散射线。

综上可知，现有技术在实际使用上，显然存在不便与缺陷，所以有必要加以改进。

发明内容

针对上述的缺陷，本发明的目的在于提供一种放射线的剂量计算方法及系统，将蒙特卡罗仿真模拟和笔形束模型相结合的进行剂量计算，可大大降低计算复杂度，提高剂量计算速度。

为了实现上述目的，本发明提供一种放射线的剂量计算方法，包括：

A、将笔束核分为主核和散射核两部份，分别构建主剂量模型和散射剂量模型；

B、利用预先建立的接受放射线照射治疗的模体的模型及预设的放射治疗深度算法，计算所述放射线照射的深度；

C、根据所述放射线照射的深度对所述主剂量模型和散射剂量模型分别进行蒙特卡罗仿真模拟，并根据所述蒙特卡罗仿真模拟的结果计算所述笔束核的笔束核值；

D、根据所述放射线照射的深度和所述笔束核值计算放射线照射所述放射线照射治疗的模体的剂量。

根据本发明所述的放射线的剂量计算方法，在所述步骤A之前或之后包括步骤：

E、采集所述需要进行放射治疗的人体部位的图像数据，并根据所述图像数据建立所述模型；或者

F、根据仿真照射预定物质的设置，建立所述模体的模型。

根据本发明所述的放射线的剂量计算方法，在所述步骤A中，根据光子沉积理论，将笔束核分为主核和散射核两部份；其中，所述主核为单位质量内，一次带电粒子和第一次跟介质交互的光子所释放的二次带电粒子所沉积的能量；所述散射核为单位质量内，多次跟介质交互的、来自韧致辐射和湮灭产生的光子所释放的带电粒子所沉积的能量。

根据本发明所述的放射线的剂量计算方法，所述步骤B包括：

B1、通过等效深度法将所述模体的物理深度转换为等效水体深度，所述等效水体深度为所述放射线照射的深度；

B2、将笔形束按照相同尺寸划分为多个单元野，并根据加速器的附件信息获得每个所述单元野的强度信息，产生并获得每个单元野的位置信息。

根据本发明所述的放射线的剂量计算方法，所述步骤C包括：

C1、通过公式(1)对所述主剂量模型和散射剂量模型分别进行蒙特卡罗仿真模拟，拟合出所述放射线照射的深度z处A_z、a_z、B_z、b_z四个拟合参数；

其中，所述R为圆形的放射野的半径；

C2、通过公式(2)以及所述A_z、a_z、B_z、b_z四个参数计算所述笔束核的笔束核值；

其中，为主核，为散射核，ρ表示介质密度，P表示单位质量的能量分布，P/ρ表示单位体积的能量分布，r是所述笔形束剂量沉积核中任意点到所述笔形束的垂直距离，A_z、a_z、B_z、b_z是深度z上的拟合参数。

根据本发明所述的放射线的剂量计算方法，所述步骤D包括：

D1、根据所述模型内任意一点与所述单元野的位置关系调用对应的所述笔束核的笔束核值，并乘以所述单元野的强度，获得所述单元野对该点的剂量贡献值，再卷积叠加所有所述单元野的剂量贡献值，获得该点的放射线照射的剂量；

D2、每个所述单元野根据所述放射线照射的跳数信息确定相应权重，将同一方向的所有所述单元野的权重叠加，形成一个所述放射线照射的方向的光通分布，所述光通分布和所述模型内所述的点的剂量沉积核卷积叠加，获得不规则野的剂量分布；所述不规则野的剂量分布为所述放射线照射治疗的模体的剂量。

根据本发明所述的放射线的剂量计算方法，所述步骤D1中，所述模体内任意一点的放射线照射的剂量的计算公式为：

其中，r为所述模型内的任意一点，s表示射野，Φ_E表示能量为E的粒子的能量通量，Π是笔核值。

根据本发明所述的放射线的剂量计算方法，在所述步骤E中，所述需要进行放射治疗的人体部位仅包括放疗靶区和危及器官；在根据所述图像数据建立所述模型后，将所述模型置入放疗系统等中心上；

在所述步骤F中，所述模型为水箱模型，水骨模型，水肺模型；并且在建立所述模体的模型后，将所述模体放在等中心上。

根据本发明所述的放射线的剂量计算方法，所述笔形束是从放射源射出的无穷小尺寸射线，所述笔束核是单位时间照射的所述笔形束在无限大均匀介质中的单位能量分布沉积。

本发明还提供一种实现所述放射线的剂量计算方法的系统，所述系统包括：

模型构建模块，用于将笔束核分为主核和散射核两部份，分别构建主剂量模型和散射剂量模型；

放射线照射的深度计算模块，用于利用预先建立的接受放射线照射治疗的模体的模型及预设的放射治疗深度算法，计算所述放射线照射的深度；

仿真计算模块，用于根据所述放射线照射的深度对所述主剂量模型和散射剂量模型分别进行蒙特卡罗仿真模拟，并根据所述蒙特卡罗仿真模拟的结果计算所述笔束核的笔束核值；

剂量计算模块，根据所述放射线照射的深度和所述笔束核值计算放射线照射所述放射线照射治疗的模体的剂量。

本发明通过使用高精度的将蒙特卡罗仿真模拟和笔形束模型相结合的剂量计算算法，基于一种笔束核的解析描述，阐述了一种半解析光子笔形束积分方法，首先，将笔束核分为主核和散射核两部份，分别构建主剂量模型和散射剂量模型；利用预先建立的接受放射线照射治疗的模体的模型及预设的放射治疗深度算法，计算所述放射线照射的深度；然后，根据所述放射线照射的深度对所述主剂量模型和散射剂量模型分别进行蒙特卡罗仿真模拟，并根据所述蒙特卡罗仿真模拟的结果计算所述笔束核的笔束核值；最后，根据所述放射线照射的深度和所述笔束核值计算放射线照射所述放射线照射治疗的模体的剂量。蒙特卡罗仿真模拟和笔形束模型相结合的剂量计算算法克服了传统的直接测量的方法获取单元小射束的能量沉积核。它也没有像传统的蒙特卡罗模拟一样保留模体各个点的蒙卡模拟结果，形成庞大的沉积核数据文件，而是只需模拟射线在单元野垂直入射到模体的情况，并且最后只需保留模体的中心轴和离轴的剂量信息为笔束核的数据。借此，本发明的放射线的剂量计算方法可大大降低计算复杂度，提高剂量计算速度，在保证速度的同时也提高了准确度。

附图说明

图1是本发明放射线的剂量计算系统的结构示意图；

图2是本发明优选放射线的剂量计算系统的结构示意图；

图3是本发明放射线的剂量计算方法的流程图；

图4是本发明放射线的剂量计算方法优选实施例的流程图之一；

图5是本发明放射线的剂量计算方法优选实施例的流程图之二。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。

参见图1和图2，所述放射线的剂量计算系统100主要包括模型构建模块10、放射线照射的深度计算模块20、仿真计算模块30、剂量计算模块40，其中：

模型构建模块10，根据光子沉积理论，用于将笔束核分为主核和散射核两部份，分别构建主剂量模型和散射剂量模型；主要依据是带电粒子与物质相互作用的主要过程包括与核外电子发生非弹性碰撞产生的电离、激发以及与原子核发生非弹性碰撞产生的韧致辐射以及正负电子对湮没辐射；带电粒子的能量损耗一部分发射光学光谱和X线与光能的形式释放，即产生光子。光子本身是电磁辐射过程的表现，与带电粒子相比，光子与物质相互作用表现出不同的特点：光子不直接沉积能量，而是先把能量传递给碰撞生成的带电粒子，带电粒子与模体多次作用后被模体吸收，将能量沉积在模体中。所述主核为单位质量内，一次带电粒子和第一次跟介质交互的光子所释放的二次带电粒子所沉积的能量；所述散射核为单位质量内，多次跟介质交互的、来自韧致辐射和湮灭产生的光子所释放的带电粒子所沉积的能量。

放射线照射的深度计算模块20，用于利用预先建立的接受放射线照射治疗的模体的模型及预设的放射治疗深度算法，计算所述放射线照射的深度。

仿真计算模块30，用于根据所述放射线照射的深度对所述主剂量模型和散射剂量模型分别进行蒙特卡罗仿真模拟，并根据所述蒙特卡罗仿真模拟的结果计算所述笔束核的笔束核值；

剂量计算模块40，根据所述放射线照射的深度和所述笔束核值计算放射线照射所述放射线照射治疗的模体的剂量。

图2是本发明优选放射线的剂量计算系统的结构示意图，所述系统100包括：模型构建模块10、放射线照射的深度计算模块20、仿真计算模块30、剂量计算模块40、模体模型建立模块50、图片预处理模块60、人体建模模块70。

图片预处理模块60，用于根据来自医学图像的数据进行图像预处理，包括配准、融合、校正、滤波和增强等，然后对皮肤边界、病灶靶区、危及器官进行精确的分割。

人体建模模块70，用于将医学图像的数据转化为不同剂量计算方法所需的模型，包括三维几何(解剖结构、位置属性)和物理信息(材料组成、密度信息)，真实再现人体复杂解剖结构及元素组成。

模体模型建立模块50，用于采集所述需要进行放射治疗的人体部位的图像数据，并根据所述图像数据建立所述模型，具体的是，所述需要进行放射治疗的人体部位仅包括放疗靶区和危及器官；在根据所述图像数据建立所述模型后，将所述模型置入放疗系统等中心上；或者

模体模型建立模块50，用于根据仿真照射预定物质的设置，建立所述模体的模型，所述模型优选为水箱模型，水骨模型，水肺模型；并且在建立所述模体的模型后，将所述模体放在等中心上。

模型构建模块10，根据光子沉积理论，用于将笔束核分为主核和散射核两部份，分别构建主剂量模型和散射剂量模型。

放射线照射的深度计算模块20还包括：

深度等效子模块21，用于通过等效深度法将所述模体的物理深度转换为等效水体深度，所述等效水体深度为所述放射线照射的深度；

笔形束划分子模块22，用于将笔形束按照相同尺寸划分为多个单元野，并根据加速器的附件信息获得每个所述单元野的强度信息，产生并获得每个单元野的位置信息。所述笔形束是从放射源射出的无穷小尺寸射线。

仿真计算模块30，用于根据所述放射线照射的深度对所述主剂量模型和散射剂量模型分别进行蒙特卡罗仿真模拟，并根据所述蒙特卡罗仿真模拟的结果计算所述笔束核的笔束核值；仿真计算模块30通过公式(1)对所述主剂量模型和散射剂量模型分别进行蒙特卡罗仿真模拟，拟合出所述放射线照射的深度z处A_z、a_z、B_z、b_z四个拟合参数；公式(1)如下：

其中，所述R为圆形的放射野的半径，在设定照射参数时设定；并根据所述蒙特卡罗仿真模拟的结果计算所述笔束核的笔束核值；公式(1)物理意义为，把半径为R的圆形放射野上所有无穷小笔形束在中心轴处的剂量积分，最终积分结果表示笔形束在中心轴所沉积的剂量。只需要给出多个R值，记录每次仿真得到的中心轴剂量，就可以拟合出深度z处A_z、a_z、B_z、b_z四个参数。由于公式(1)把中心轴剂量表达为R的函数，使得当笔束核被用于剂量计算时增加了拟合准确度。同时也避免了当r趋向于0时所带来的奇异性问题。仿真计算模块30通过公式(2)以及所述A_z、a_z、B_z、b_z四个参数计算所述笔束核的笔束核值；兆伏级光子束产生的笔束核可以用公式(2)所示的两个沿半径方向的指数分布高精度解析，式中有四个参数，不同能量光子束产生的笔束核参数不同，通过一些仿真平台得到具体加速器的具体能量的笔束核参数。进行蒙特卡罗仿真模拟的时候，主核和散射核被分开模拟，公式(2)如下：其中，为主核，为散射核，ρ表示介质密度，P表示单位质量的能量分布，P/ρ表示单位体积的能量分布，r是笔形束剂量沉积核(笔束核)中任意点到所述笔形束的垂直距离，A_z、a_z、B_z、b_z是深度z上的拟合参数。由于通过蒙特卡罗仿真模拟的结果计算所述笔束核的笔束核值包含了加速器能谱的信息，从而保证剂量计算的准确性。

剂量计算模块40还包括：

第一剂量计算子模块41，用于根据所述模型内任意一点与所述单元野的位置关系调用对应的所述笔束核的笔束核值，并乘以所述单元野的强度，获得所述单元野对该点的剂量贡献值，再卷积叠加所有所述单元野的剂量贡献值，获得该点的放射线照射的剂量；

第二剂量计算子模块42，用于将每个所述单元野根据所述放射线照射的跳数信息确定相应权重，将同一方向的所有所述单元野的权重叠加，形成一个所述放射线照射的方向的光通分布，所述光通分布和所述模型内所述的点的剂量沉积核卷积叠加，获得不规则野的剂量分布；所述不规则野的剂量分布为所述放射线照射治疗的模体的剂量。剂量计算模块40具体的剂量计算公式如下所示：所述模体内任意一点的放射线照射的剂量的计算公式为：

其中，r为所述模型内的任意一点，s表示射野，Φ_E表示能量为E的粒子的能量通量，Π是上述公式(2)预先计算的笔核值。

参见图3，所述方法可由如图1或图2所示的放射线的剂量计算系统100实现，所述放射线的剂量计算方法包括：

步骤S301，将笔束核分为主核和散射核两部份，分别构建主剂量模型和散射剂量模型，优选的是，根据光子沉积理论，将笔束核分为主核和散射核两部份；其中，所述主核为单位质量内，一次带电粒子和第一次跟介质交互的光子所释放的二次带电粒子所沉积的能量；所述散射核为单位质量内，多次跟介质交互的、来自韧致辐射和湮灭产生的光子所释放的带电粒子所沉积的能量；所述笔束核是单位时间照射的所述笔形束在无限大均匀介质中的单位能量分布沉积。本步骤可由模型构建模块10实现。

步骤S302，利用预先建立的接受放射线照射治疗的模体的模型及预设的放射治疗深度算法，计算所述放射线照射的深度，模体是指置入接受治疗的模体，本步骤可由放射线照射的深度计算模块20实现。

步骤S303，根据所述放射线照射的深度对所述主剂量模型和散射剂量模型分别进行蒙特卡罗仿真模拟，并根据所述蒙特卡罗仿真模拟的结果计算所述笔束核的笔束核值，本步骤可由仿真计算模块30实现.

步骤S304，根据所述放射线照射的深度和所述笔束核值计算放射线照射所述放射线照射治疗的模体的剂量，本步骤可由剂量计算模块40实现。

参见图4，所述方法可由如图1或图2所示的放射线的剂量计算系统100实现，所述方法包括步骤如下：

步骤S401，根据来自医学图像的数据进行图像预处理，包括配准、融合、校正、滤波和增强等，然后对皮肤边界、病灶靶区、危及器官进行精确的分割。本步骤可由图片预处理模块60实现。

步骤S402，将医学图像的数据转化为不同剂量计算方法所需的模型，包括三维几何(解剖结构、位置属性)和物理信息(材料组成、密度信息)，真实再现人体复杂解剖结构及元素组成。本步骤可由人体建模模块70实现。

步骤S403，如果是仿真照射预定物质供研究使用，则根据仿真照射预定物质的设置，可以通过一些平台直接制作模体的模型，比如水箱模型，水骨模型，水肺模型等，并且在建立所述模体的模型后，将所述模体放在等中心上。或者对于人体仿真来看说，则采集所述需要进行放射治疗的人体部位的图像数据，并根据所述图像数据建立所述模型，所述模型优选为三维模型，所述需要进行放射治疗的人体部位仅包括放疗靶区和危及器官；在根据所述图像数据建立所述模型后，将所述模型置入放疗系统等中心上。本步骤可由模体模型建立模块50在步骤S404之前或之后实现。

步骤S404，将笔束核分为主核和散射核两部份，分别构建主剂量模型和散射剂量模型。

步骤S405，通过等效深度法将所述模体的物理深度转换为等效水体深度，所述等效水体深度为所述放射线照射的深度；所谓等效放射深度法，就是通过把介质中的物理深度按照相对于水的衰减系数比做缩放，得到等效放射深度。本步骤可由深度等效子模块21实现。

步骤S406，将笔形束按照相同尺寸划分为多个单元野，并根据加速器的附件信息获得每个所述单元野的强度信息，产生并获得每个单元野的位置信息，供后续仿真计算模块30和剂量计算模块40调用。所述笔形束是从放射源射出的无穷小尺寸射线。本步骤可由笔形束划分子模块22实现。

步骤S407，根据所述放射线照射的深度对所述主剂量模型和散射剂量模型分别进行蒙特卡罗仿真模拟，优选的是，通过公式(1)对所述主剂量模型和散射剂量模型分别进行蒙特卡罗仿真模拟，拟合出所述放射线照射的深度z处A_z、a_z、B_z、b_z四个拟合参数；公式(1)如下：

其中，所述R为圆形的放射野的半径，在设定照射参数时设定；并根据所述蒙特卡罗仿真模拟的结果计算所述笔束核的笔束核值；公式(1)物理意义为，把半径为R的圆形放射野上所有无穷小笔形束在中心轴处的剂量积分，最终积分结果表示笔形束在中心轴所沉积的剂量。只需要给出多个R值，记录每次仿真得到的中心轴剂量，就可以拟合出深度z处A_z、a_z、B_z、b_z四个参数。由于公式(1)把中心轴剂量表达为R的函数，使得当笔束核被用于剂量计算时增加了拟合准确度。同时也避免了当r趋向于0时所带来的奇异性问题。并根据所述蒙特卡罗仿真模拟的结果计算所述笔束核的笔束核值；具体的是，通过公式(2)以及所述A_z、a_z、B_z、b_z四个参数计算所述笔束核的笔束核值；兆伏级光子束产生的笔束核可以用公式(2)所示的两个沿半径方向的指数分布高精度解析，式中有四个参数，不同能量光子束产生的笔束核参数不同，通过一些仿真平台得到具体加速器的具体能量的笔束核参数。进行蒙特卡罗仿真模拟的时候，主核和散射核被分开模拟，公式(2)如下：其中，为主核，为散射核，ρ表示介质密度，P表示单位质量的能量分布，P/ρ表示单位体积的能量分布，r是笔形束剂量沉积核(笔束核)中任意点到所述笔形束的垂直距离，A_z、a_z、B_z、b_z是深度z上的拟合参数。由于通过蒙特卡罗仿真模拟的结果计算所述笔束核的笔束核值包含了加速器能谱的信息，从而保证剂量计算的准确性。本步骤可由仿真计算模块30实现。

步骤S408，根据所述模型内任意一点与所述单元野的位置关系调用对应的所述笔束核的笔束核值，并乘以所述单元野的强度，获得所述单元野对该点的剂量贡献值，再卷积叠加所有所述单元野的剂量贡献值，获得该点的放射线照射的剂量；剂量计算模型主要由两部分组成：笔束核入射粒子能量通量。笔束核描述了截面无限小的射束入射到半无限大均匀介质中的能量分布，可用蒙特卡罗仿真模拟得到；入射粒子能量通量是指模体或病人体内相应空间点所接受的粒子能量通量，两者进行卷积得到任意空间点的剂量。

步骤S409，将每个所述单元野根据所述放射线照射的跳数信息确定相应权重，将同一方向的所有所述单元野的权重叠加，形成一个所述放射线照射的方向的光通分布，所述光通分布和所述模型内所述的点的剂量沉积核卷积叠加，获得不规则野的剂量分布；所述不规则野的剂量分布为所述放射线照射治疗的模体的剂量。则剂量计算公式如下所示：所述模体内任意一点的放射线照射的剂量的计算公式为：其中，r为所述模型内的任意一点，s表示射野，Φ_E表示能量为E的粒子的能量通量，Π是上述公式(2)预先计算的笔核值。

解析剂量算法认为能量局部沉积(次级电子不输运)，或者采用某种简单的方法近似处理次级电子的输运过程，存在明显局限性。此类解析剂量算法基本上都是假设介质为半无穷大的前提下，得到均匀模体中的剂量分布特点表达式或数据库，同时采用半经验方法处理不均匀介质的存在，计算光子或电子剂量分布会带来显著误差。而蒙特卡罗仿真模拟通过对粒子与物质相互作用进行随机模拟，来获得粒子在人体组织中沉积能量的分布。虽然目前蒙特卡罗仿真模拟已在医学剂量计算及验证方面得到广泛的应用，并被公认为当前所有剂量算法中最精确的一种，被称为医学领域剂量计算的“金标准”。但是该方法存在收敛速度慢、计算时间长的致命弱点，长期以来一直无法为临床放射治疗所接受，只能作为其他算法的验证基准，或者是用来为卷积/迭加模型计算能量沉积核——只需要一次计算获得数据库，在实际计算中不需要蒙卡计算等。考虑到蒙特卡罗计算时间过长，临床使用的放疗计划系统100几乎没有采用蒙特卡罗。而卷积/叠加算法在带电粒子失衡条件下精度不高，常规的卷积/叠加算法也存在计算速度慢的弱点。规则束模型及其修正模型速度上能够满足，但精度是一个不可忽视的问题。传统笔形束算法使用的笔束核则来自于经验的百分深度剂量和百分离轴剂量分布。本发明基于将蒙特卡罗仿真模拟和笔形束模型相结合的剂量计算算法研究，通过Gate仿真平台高精度拟合笔形束剂量沉积核，即笔束核。通过将不规则射野划分为若干个有限截面大小的笔形束，将笔束核有限截面的笔形束卷积计算患者人体三维剂量分布。为进一步提高剂量计算速度，仅计算放疗靶区和危及器官剂量分布，避免全身三维卷积计算。为了计算非均匀人体中的剂量分布，采样水等效放射深度沿射线路径单方向进行卷积。

参见图5，所述方法可由如图1或图2所示的放射线的剂量计算系统100实现，具体实施例实施过程包括：

步骤S501，根据来自医学图像的数据进行图像预处理，包括配准、融合、校正、滤波和增强等，然后对皮肤边界、病灶靶区、危及器官进行精确的分割。

步骤S502，将医学图像的数据转化为不同剂量计算方法所需的模型，包括三维几何(解剖结构、位置属性)和物理信息(材料组成、密度信息)，真实再现人体复杂解剖结构及元素组成。

步骤S503，设置笔形束照射参数，包括圆形的放射野的半径R。

步骤S504，将笔束核分为主核和散射核两部份，分别构建主剂量模型和散射剂量模型。

步骤S505，通过等效深度法将所述模体的物理深度转换为等效水体深度，所述等效水体深度为所述放射线照射的深度。

步骤S506，根据步骤S503设置的笔形束照射参数，获得笔形束信息；将笔Ψ_i形束按照相同尺寸划分为多个单元野，调用对应的蒙特卡罗笔束核库进行蒙特卡罗仿真模拟获得笔束核值，并根据加速器的附件信息获得每个所述单元野的强度信息，产生并获得每个单元野的位置信息。

步骤S507，根据步骤S506的每个所述单元野的强度信息以及每个单元野的位置信息调用对应的笔束核值，并乘以所述单元野的强度，获得所述单元野对该点的剂量贡献值，再卷积叠加所有所述单元野的剂量贡献值，获得该点的放射线照射的剂量。

该有限笔束剂量计算公式如公式(3)所示，第一项表示光通分布，第二项表示沉积核在有限截面s(i，z)上积分所形成的有限笔束，两者相乘叠加得到计算点的剂量分布。i表示不同的有限笔束，Ψ_i表示从光通分布中得到的单元野权重，s(i，z)表示单元野确定的有限笔束在深度z上的截面，为公式(2)所示的主剂量核和散射剂量核。

步骤S508，输出剂量分布计算结果。

在本发明的一个实施例中，通过计算一例鼻咽癌病人强调放疗计划的剂量分布，对比分析该病人强调放疗结构集(靶区和危及器官)中定义的感兴趣区域剂量。

区域	实际剂量Gy	计算剂量Gy	误差
				PTV	62.663	62.795	0.040
脑干	14.988	14.107	0.029
				左腮	29.389	28.618	0.027
视交叉	58.999	60.030	0.022

表1

表1为计划靶区和危及器官误差统计表，其中的第一列为剂量统计区域，分别统计了医生勾画的PTV和危及器官(脑干、左腮、视交叉)；第2、3列为对应区域剂量的平均值；最后一列是剂量计算综合误差。

从表中可以看出，无论是靶区还是危及器官区域，综合误差都控制在5％之内，符合临床放疗对于剂量计算的精度要求。

综上所述，本发明通过使用高精度的将蒙特卡罗仿真模拟和笔形束模型相结合的剂量计算算法，基于一种笔束核的解析描述，阐述了一种半解析光子笔形束积分方法，首先，将笔束核分为主核和散射核两部份，分别构建主剂量模型和散射剂量模型；利用预先建立的接受放射线照射治疗的模体的模型及预设的放射治疗深度算法，计算所述放射线照射的深度；然后，根据所述放射线照射的深度对所述主剂量模型和散射剂量模型分别进行蒙特卡罗仿真模拟，并根据所述蒙特卡罗仿真模拟的结果计算所述笔束核的笔束核值；最后，根据所述放射线照射的深度和所述笔束核值计算放射线照射所述放射线照射治疗的模体的剂量。蒙特卡罗仿真模拟和笔形束模型相结合的剂量计算算法克服了传统的直接测量的方法获取单元小射束的能量沉积核。它也没有像传统的蒙特卡罗模拟一样保留模体各个点的蒙卡模拟结果，形成庞大的沉积核数据文件，而是只需模拟射线在单元野垂直入射到模体的情况，并且最后只需保留模体的中心轴和离轴的剂量信息为笔束核的数据。借此，本发明的放射线的剂量计算方法可大大降低计算复杂度，提高剂量计算速度，在保证速度的同时也提高了准确度。

当然，本发明还可有其它多种实施例，在不背离本发明精神及其实质的情况下，熟悉本领域的技术人员当可根据本发明作出各种相应的改变和变形，但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。

Claims (10)

1.一种放射线的剂量计算方法，其特征在于，包括：

A、将笔束核分为主核和散射核两部份，分别构建主剂量模型和散射剂量模型；

B、利用预先建立的接受放射线照射治疗的模体的模型及预设的放射治疗深度算法，计算所述放射线照射的深度；

C、根据所述放射线照射的深度对所述主剂量模型和散射剂量模型分别进行蒙特卡罗仿真模拟，并根据所述蒙特卡罗仿真模拟的结果计算所述笔束核的笔束核值；

D、根据所述放射线照射的深度和所述笔束核值计算放射线照射所述放射线照射治疗的模体的剂量。

2.根据权利要求1所述的放射线的剂量计算方法，其特征在于，在所述步骤A之前或之后包括步骤：

E、采集所述需要进行放射治疗的人体部位的图像数据，并根据所述图像数据建立所述模型；或者

F、根据仿真照射预定物质的设置，建立所述模体的模型。

3.根据权利要求1所述的放射线的剂量计算方法，其特征在于，在所述步骤A中，根据光子沉积理论，将笔束核分为主核和散射核两部份；其中，所述主核为单位质量内，一次带电粒子和第一次跟介质交互的光子所释放的二次带电粒子所沉积的能量；所述散射核为单位质量内，多次跟介质交互的、来自韧致辐射和湮灭产生的光子所释放的带电粒子所沉积的能量。

4.根据权利要求3所述的放射线的剂量计算方法，其特征在于，所述步骤B包括：

B1、通过等效深度法将所述模体的物理深度转换为等效水体深度，所述等效水体深度为所述放射线照射的深度；

B2、将笔形束按照相同尺寸划分为多个单元野，并根据加速器的附件信息获得每个所述单元野的强度信息，产生并获得每个单元野的位置信息。

5.根据权利要求3所述的放射线的剂量计算方法，其特征在于，所述步骤C包括：

C1、通过公式(1)对所述主剂量模型和散射剂量模型分别进行蒙特卡罗仿真模拟，拟合出所述放射线照射的深度z处A_z、a_z、B_z、b_z四个拟合参数；

其中，所述R为圆形的放射野的半径；

C2、通过公式(2)以及所述A_z、a_z、B_z、b_z四个参数计算所述笔束核的笔束核值；

其中，为主核，为散射核，ρ表示介质密度，P表示单位质量的能量分布，P/ρ表示单位体积的能量分布，r是所述笔形束剂量沉积核中任意点到所述笔形束的垂直距离，A_z、a_z、B_z、b_z是深度z上的拟合参数。

6.根据权利要求5所述的放射线的剂量计算方法，其特征在于，所述步骤D包括：

D1、根据所述模型内任意一点与所述单元野的位置关系调用对应的所述笔束核的笔束核值，并乘以所述单元野的强度，获得所述单元野对该点的剂量贡献值，再卷积叠加所有所述单元野的剂量贡献值，获得该点的放射线照射的剂量；

D2、每个所述单元野根据所述放射线照射的跳数信息确定相应权重，将同一方向的所有所述单元野的权重叠加，形成一个所述放射线照射的方向的光通分布，所述光通分布和所述模型内所述的点的剂量沉积核卷积叠加，获得不规则野的剂量分布；所述不规则野的剂量分布为所述放射线照射治疗的模体的剂量。

7.根据权利要求6所述的方法，其特征在于，所述步骤D1中，所述模体内任意一点的放射线照射的剂量的计算公式为：

其中，r为所述模型内的任意一点，s表示射野，Φ_E表示能量为E的粒子的能量通量，Π是笔核值。

8.根据权利要求2所述的方法，其特征在于，在所述步骤E中，所述需要进行放射治疗的人体部位仅包括放疗靶区和危及器官；在根据所述图像数据建立所述模型后，将所述模型置入放疗系统等中心上；

在所述步骤F中，所述模型为水箱模型，水骨模型，水肺模型；并且在建立所述模体的模型后，将所述模体放在等中心上。

9.根据权利要求4所述的方法，其特征在于，所述笔形束是从放射源射出的无穷小尺寸射线，所述笔束核是单位时间照射的所述笔形束在无限大均匀介质中的单位能量分布沉积。

10.一种包括如权利要求1～9任一项所述的放射线的剂量计算方法的系统，其特征在于，所述系统包括：

模型构建模块，用于将笔束核分为主核和散射核两部份，分别构建主剂量模型和散射剂量模型；

放射线照射的深度计算模块，用于利用预先建立的接受放射线照射治疗的模体的模型及预设的放射治疗深度算法，计算所述放射线照射的深度；

仿真计算模块，用于根据所述放射线照射的深度对所述主剂量模型和散射剂量模型分别进行蒙特卡罗仿真模拟，并根据所述蒙特卡罗仿真模拟的结果计算所述笔束核的笔束核值；

剂量计算模块，根据所述放射线照射的深度和所述笔束核值计算放射线照射所述放射线照射治疗的模体的剂量。