一种放疗计划强度分布优化方法、装置及放射治疗设备与流程

1.本发明属于放射治疗技术领域，具体涉及一种放疗计划强度分布优化方法、装置及放射治疗设备。


背景技术：

2.逆向计划调强放射治疗技术是现在医疗机构放疗科采用的主流技术，它的优化方法有：两步法、一步法(直接子野优化)。两步法指的是：强度分布优化、光栅子野分割。强度分布优化是根据患者影像结合临床要求，求解出能满足临床要求的最优射线强度分布的过程；光栅子野分割是将最优强度分布转化为物理上多叶光栅(multi
‑
layerscollimator，mlc)可执行的一系列子野的过程。在强度分布优化 (fluence map optimization，fmo)中，每个连续的射束被分成很多细小的束(beamlets)，分别考虑每个小束在介质中的剂量沉积，使用梯度下降法等最优化方法求解每个小束的强度值。最后对于特定的某个点，将每个射束对它的剂量贡献叠加，得到该点的总剂量。
3.fmo的一个重要步骤是计算射线强度在人体组织内的剂量沉积，一般软件均采用笔形束法(pencil beam，pb)。pb算法的剂量值是笔形束核与强度分布的二维卷积结果，笔形束核描述的是一束入射在无限大均匀介质中的能量扩散模式。求解笔形束核的方法有：mc模拟、直接从测量数据中得到、反卷积从测量的宽束射线剂量分布中得到、有限笔形束法。针对非均匀介质的情况，考虑了介质的3d密度校正、射线在到达模体前的散射和电子污染。尽管pb算法可以较好地处理不规则射野的情况，但在人体非均匀组织和体表不规则轮廓上计算精度仍然不够，尤其对于气管、喉这样的空腔组织，剂量差可达 100％，如图2所示。为了更好地模拟射线在非均匀组织中的剂量沉积，许多研究提出了原射线注量校正、核函数校正等等方法，但这些本质都是数据拟合的方法，精度远不及mc方法。
4.mc算法又称计算机随机模拟方法，其原理基于随机数的生成，通过模拟相当大数目的粒子在介质中发生相互作用的过程，得到通量、能谱、剂量在整个介质中的分布。由于mc是逼真地模拟大量粒子(多达上亿个)能量沉积的过程，所以模拟的结果是公认最精确的。
5.商用软件的fmo大多考虑到运算速度而使用pb算法，由于pb 算法自身的精度不高，优化得到的最优强度分布是伪最优的，由此强度分布经最终mc算法计算得到的剂量分布不一定符合预期的临床要求，其结果与临床要求设置的初衷是相悖的。因此，在计划优化中采用mc算法是一个重要的研究方向。
6.多数研究者使用egs4等工具，预先生成所有beamlets的mc沉积矩阵，在强度分布优化中直接调用数值即可。直接在cpu上执行 mc剂量计算的效率很低，为了加快运算速度，mc剂量计算在gpu 上的应用逐渐增多，发展出了基于gpu的mc剂量计算软件库，如gdpm，gpumcd。基于cpu的mc优化提供了解决方法，在gpu 上的应用可以有效地提高计算速度，但对于强度分布优化而言，速度依旧太慢，不能满足商用软件的要求。此外，这两种方法需要将沉积矩阵存储在计算上，占用大量硬盘和内存空间，且每个病人都需要计算一次，时间开
销很大，不能满足每日有大量病人的医疗机构。


技术实现要素：

7.为了解决上述技术问题，本发明提出了一种放疗计划强度分布优化方法、装置及放射治疗设备。
8.为了达到上述目的，本发明的技术方案如下：
9.一方面，本发明公开一种放疗计划强度分布优化方法，包括以下步骤：
10.s1：通过神经网络，将pb剂量映射到mc剂量上；
11.s2：求解mc剂量相对于pb剂量的梯度，得到基于mc算法的沉积矩阵；
12.s3：将基于mc算法的沉积矩阵带入fmo中，实现放疗计划强度分布优化。
13.在上述技术方案的基础上，还可做如下改进：
14.作为优选的方案，s1具体包括以下步骤：
15.s1.1：数据准备，选取n1个病人ct，每个病人添加n2个射野，分别计算每个射野的pb剂量和mc剂量，形成n1*n2组数据，每组数据均包括对应的ct、pb剂量和mc剂量；
16.s1.2：数据划分，将n1*n2组数据按比例划分为训练集测试集；
17.s1.3：创建卷积神经网络；
18.s1.4：采用训练集数据训练创建好的卷积神经网络；
19.s1.5：在训练好的卷积神经网络上测试测试集数据，得到预测剂量，计算预测剂量和mc剂量的gamma通过率。
20.作为优选的方案，卷积神经网络包括：输入层、融合层、卷积层、激活层、池化层、上采样层、输出层。
21.作为优选的方案，
[0022][0023][0024]
其中，a为基于mc算法的沉积矩阵；
[0025]
d
pb
与d
mc
分别为基于pb算法的剂量与基于mc算法的剂量；
[0026]
a
jm
为基于mc算法的沉积矩阵的每一个值；
[0027]
代表该体素点(第j个)的mc剂量相对于第k个pb 剂量的导数；
[0028]
dd
pb,k
/dx
m
代表第k个pb剂量相对于第m个强度值的导数；
[0029]
n
dose
代表剂量个数；
[0030]
由上式可得，通过对pb剂量到mc剂量映射的卷积神经网络的梯度求解，能够得到基于mc算法的沉积矩阵。
[0031]
作为优选的方案，fmo具体包括以下步骤：
[0032]
s3.1：将射野分割成笔形束；
[0033]
s3.2：引入基于mc算法的沉积矩阵；
[0034]
s3.3：设置初始强度值，将其作为优化的初始值；
[0035]
s3.4：利用基于mc算法的沉积矩阵求解目标函数的导数；
[0036]
s3.5：求解当前梯度下的最优步长；
[0037]
s3.6：更新强度值、剂量和目标函数；
[0038]
s3.7：重复步骤s3.4
‑
s3.6，直至满足停止条件后，迭代停止。
[0039]
作为优选的方案，还包括以下步骤：
[0040]
s4：根据放疗计划强度分布优化方法计算人体内的剂量和各个器官的dvh，评估计划效果。
[0041]
另一方面，本发明还公开一种放疗计划强度分布优化装置，包括：
[0042]
神经网络映射模块，用于通过神经网络，将pb剂量映射到mc 剂量上；
[0043]
沉积矩阵获取模块，用于求解mc剂量相对于pb剂量的梯度，得到基于mc算法的沉积矩阵；
[0044]
强度分布优化模块，用于将基于mc算法的沉积矩阵带入fmo 中，实现放疗计划强度分布优化。
[0045]
另一方面，本发明还公开放射治疗设备，用于执行上述任一种放疗计划强度分布优化方法，或，包括：上述放疗计划强度分布优化装置。
[0046]
本发明一种放疗计划强度分布优化方法、装置及放射治疗设备具有以下有益效果：
[0047]
本发明将mc算法应用在fmo中可以确保优化过程得到的强度分布是真实可靠的，使用这样的强度分布进行多叶光栅子野分割，得到的子野走位是更加合理的，保证了最终得到的计划剂量分布更接近于临床要求，计划效果有显著提升。
附图说明
[0048]
为了更清楚地说明本发明实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本发明的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。
[0049]
图1为本发明实施例提供的放疗计划强度分布优化方法的流程图。
[0050]
图2为mc剂量与pb剂量的对比图；
[0051]
图2(a)为mc剂量图；
[0052]
图2(b)为pb剂量图。
[0053]
图3为本发明实施例提供的训练集中的一层数据图；
[0054]
图3(a)为ct图；
[0055]
图3(b)为pb剂量数据图；
[0056]
图3(c)为mc剂量数据图。
[0057]
图4为本发明实施例提供的强度分布图优化流程图。
[0058]
图5为本发明实施例提供的pb剂量到mc剂量的映射方法流程图。
[0059]
图6为本发明实施例提供的卷积神经网络的结构示意图。
[0060]
图7为本发明实施例提供的卷积神经网络预测的neural剂量 gamma通过率与pb剂量和gamma通过率对比图；
[0061]
图7(a)为mc剂量图；
[0062]
图7(b)为mc剂量分布图；
[0063]
图7(c)为groundtruth图；
[0064]
图7(d)为pb剂量图；
[0065]
图7(e)为pb剂量分布图；
[0066]
图7(f)为pb剂量gamma通过率图；
[0067]
图7(g)为卷积神经网络预测的neural剂量图；
[0068]
图7(h)为卷积神经网络预测的neural剂量分布图；
[0069]
图7(i)为卷积神经网络预测的neural剂量gamma通过率图。
[0070]
图8为本发明实施例提供的3*3卷积核对4*4矩阵(外侧补0) 卷积过程图。
[0071]
图9为本发明实施例提供的反池化上采样过程图。
[0072]
图10为本发明实施例提供的插值法上采样过程图。
[0073]
图11为本发明实施例提供的反卷积上采样过程图。
[0074]
图12为本发明实施例提供的基于卷积神经网络的优化dvh与基于mc的优化dvh对比图。
具体实施方式
[0075]
下面结合附图详细说明本发明的优选实施方式。
[0076]
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0077]
为了达到本发明的目的，一种放疗计划强度分布优化方法、装置及放射治疗设备的其中一些实施例中，如图1所示，该放疗计划强度分布优化方法，包括以下步骤：
[0078]
s1：通过神经网络，将pb剂量映射到mc剂量上；
[0079]
s2：求解mc剂量相对于pb剂量的梯度，得到基于mc算法的沉积矩阵；
[0080]
s3：将基于mc算法的沉积矩阵带入fmo中，实现放疗计划强度分布优化。
[0081]
值得注意的是，s1中，pb剂量和mc剂量可以但不限于通过计划系统tps计算得到。
[0082]
其中pb剂量的计算方法可根据公式或试验结果计算主射线和散射线在组织内的剂量沉积，mc剂量的计算方法是模拟数十亿个粒子在组织内的剂量沉积。
[0083]
为了进一步地优化本发明的实施效果，在另外一些实施方式中，其余特征技术相同，不同之处在于，s1具体包括以下步骤：
[0084]
s1.1：数据准备，选取n1个病人ct，每个病人添加n2个射野，分别计算每个射野的pb剂量和mc剂量，形成n1*n2组数据，每组数据均包括对应的ct、pb剂量和mc剂量；
[0085]
s1.2：数据划分，将n1*n2组数据按比例划分为训练集测试集；
[0086]
s1.3：创建卷积神经网络；
[0087]
s1.4：采用训练集数据训练创建好的卷积神经网络；
[0088]
s1.5：在训练好的卷积神经网络上测试测试集数据，得到预测剂量，计算预测剂量和mc剂量的gamma通过率。
[0089]
进一步，卷积神经网络包括：输入层、融合层、卷积层、激活层、池化层、上采样层、
输出层，如图6所示。
[0090]
在具体地实施例中，pb剂量到mc剂量的映射方法，具体可以分为如下步骤，如图5所示；
[0091]
s1.1：数据准备：
[0092]
选取50例病人，病人ct涵盖头颈、胸腹、盆腔等部位；每个病人添加72个射野(间隔0.5
°
)，射野的等中心位置尽可能分布在不同 ct层，使得射野可以覆盖更多的ct范围；
[0093]
分别计算每个射野的pb剂量和mc剂量，值得注意的是，这里剂量网络保持统一，可以采用3mm网格；
[0094]
整体来看，共有3600组ct、pb剂量、mc剂量的数据，如图3 所示；
[0095]
由于ct数据的采样间距、大小与剂量网格不一致，因此需要将 ct数据插值到与剂量网格一致的大小和间隔。
[0096]
s1.2：数据划分：将3600组数据按7：2：1的比例划分为训练集、验证集和测试集。
[0097]
s1.3：创建卷积神经网络：卷积神经网络由输入层、融合层、卷积层、激活层、池化层、上采样层、输出层组成。
[0098]
对ct与pb剂量图像执行随机裁剪(crop)，从而提高神经网络的泛化能力，处理后的图像大小为128
×
128，两者经融合层形成一组输入数据(128
×
128
×
2)，而后经64个卷积核处理，形成128
×
128
×
64 的张量，继续卷积处理，形成128
×
128
×
64的张量，经最大池化处理后形成64
×
64
×
64的张量，以此类推，直至上采样、卷积处理后得到 128
×
128
×
64的张量，再次卷积后得到128
×
128的输出图像。
[0099]
卷积神经网络的输出真值为mc剂量，损失函数定义为mc剂量和预测剂量之间的均方误差(mean square error,mse)，是目标变量和预测值之间的平方距离之和的均值：
[0100][0101]
其中：n为128
×
128。
[0102]
s1.4：训练卷积神经网络；
[0103]
网络训练的过程就是在不断迭代求解卷积核，是特征提取的过程；
[0104]
训练的具体步骤为：
[0105]
i.随机选取数据(batch size＝512)，对输入图像进行在线数据增广，方法包含平移、2d旋转，目的是增大数据量，提高卷积神经网络的泛化能力；
[0106]
ii.使用随机梯度下降法反向传播得到网络所有的权重和偏移的梯度，使用学习率10
‑5更新权重和偏移，在此基础上，得到新的预测值，与mc剂量对比，得到均方误差，重复此过程200次，完成一个epoch；
[0107]
iii.重复步骤i),ii)，直至均方误差小于阈值或达到最大epoch(如： 5000)，迭代停止，卷积神经网络训练完毕；
[0108]
上述计算过程在nvidiartx 2080ti上执行。
[0109]
在该实施例中，将数据划分为训练集、验证集和测试集。在训练的过程中，一般在训练时，几个epoch结束后跑一次验证集看下效果。从而可以及时发现模型或者参数的问题，及时终止训练，重新调参或者调整模型。此外，可以验证模型的泛化能力，如果在验证集
上的效果比训练集上差很多，就应考虑模型是否过拟合了。
[0110]
s1.5：测试效果；
[0111]
在训练好的卷积神经网络上测试，输入测试集的ct和pb剂量数据，得到预测剂量，将所有的预测剂量数据保存在本地，使用剂量验证软件(如：dvs,linatech)计算预测剂量和mc剂量的gamma 通过率(3mm，3％，10％)，通过率在95％以上即为合格。
[0112]
gamma通过率是衡量“测量”和“计算”间差异的定量方法，是吻合距离(distance to agreement,dta)和剂量差异(dose difference,dd) 的组合：
[0113][0114][0115]
r(r
m
,r
c
)＝|r
c
‑
r
m
|；
[0116]
δ(r
m
,r
c
)＝d
c
(r
c
)
‑
d
m
(r
m
)；
[0117]
其中：
[0118]
r(r
m
,r
c
)代表对于测量点r
m
与计算点r
c
间的距离；
[0119]
δ(r
m
,r
c
)代表两者之间的剂量差；
[0120]
δd
m
代表距离阈值(可以但不限于取3mm)；
[0121]
δd
m
代表剂量的阈值(可以但不限于取最大剂量的3％)；
[0122]
当γ(r
m
)≤1时，该点通过，否则不过。
[0123]
若只考虑距离因素(δ(r
m
,r
c
)＝0)，那么测量点通过的依据是：在距离阈值(如：3mm)内，存在一个计算点的剂量值等于测量点；
[0124]
若只考虑剂量因素(r(r
m
,r
c
)＝0)，那么测量点通过的依据是：在同一位置，计算剂量和测量剂量之间的差不大于剂量阈值(如：最大值的3％)。
[0125]
如图7所示，其显示某一层pb剂量和neural剂量(基于神经网络的剂量，即为卷积神经网络的预测剂量)的gamma通过率，mc剂量为真值，可见neural剂量更接近于mc剂量，统计结果显示，neural 剂量的gamma通过率达到99％，而pb剂量的gamma通过率仅有62％。
[0126]
为了进一步地优化本发明的实施效果，在另外一些实施方式中，其余特征技术相同，不同之处在于，获取基于mc算法的沉积矩阵，其原理为：
[0127][0128][0129]
其中，a为基于mc算法的沉积矩阵；
[0130]
d
pb
与d
mc
分别为基于pb算法的剂量与基于mc算法的剂量；
[0131]
a
jm
为基于mc算法的沉积矩阵的每一个值；
[0132]
代表该体素点(第j个)的mc剂量相对于第k个pb 剂量的导数；
[0133]
dd
pb,k
/dx
m
代表第k个pb剂量相对于第m个强度值的导数；
[0134]
n
dose
代表剂量个数；
[0135]
由上式可得，通过对pb剂量到mc剂量映射的卷积神经网络的梯度求解，能够得到
基于mc算法的沉积矩阵。
[0136]
值得注意的是，上述公式将问题“使用mc算法计算沉积矩阵”转化为“求解mc剂量对pb剂量的导数”，其中：pb剂量对强度值x的导数 dd
pb,k
/dx
m
是原算法得到的沉积矩阵，即tps软件里使用主射线和散射线叠加的方法计算出的矩阵。
[0137]
而“求解mc剂量对pb剂量的导数”的本质就是“pb剂量到mc 剂量映射的卷积神经网络”的梯度求解，求解方法可参考反向传播算法。
[0138]
1)求解各个类型网络的梯度，下面分别描述卷积层(带激活层)、池化层、上采样层的梯度求解。
[0139]
i.带激活层的卷积层梯度：如图8所示，记输入为x，卷积层的权重和偏置分别为w和b，输出为y，激活函数为那么有：
[0140][0141]
conv2(w,x)代表w与x的离散二维卷积，也可以写作：
[0142][0143]
因此有
[0144][0145][0146]
其中：
[0147]
k反映了卷积核大小(2*k 1＝kernel size)，当采用3
×
3卷积核时，k为1；
[0148]
激活函数若为relu函数，那么它的导数为：
[0149][0150]
激活函数若为sigmoid函数，那么它的导数为：
[0151][0152]
依照上述公式可以得到卷积层(带激活函数)的导数。
[0153]
ii.池化层的梯度：以最大池化(max pooling)，2
×
2的卷积核为例，该层的输出等于对应输入的最大值，即 y(i,j)＝max{x(p,q),p＝i
‑
1,i,q＝j
‑
1,j}，那么梯度可表示为：
[0154][0155]
iii.上采样层的梯度：上采样(up sampling)有三种方法：插值法、反卷积和反池化；
[0156]
a)如图10所示，插值法：将输入图像放大指某一尺寸，使用插值法计算每一点的像素值，以双线性插值为例，上采样后的图像某一点像素值表示为下式，其中上采样图像i,j处对应了原图的i，j:，输出对输入的梯度记为插值权重w；
[0157][0158]
b)如图11所示，反卷积：卷积可以表示为两个矩阵的乘积，反卷积就是卷积核矩阵的转置与原矩阵相乘，是小尺寸到大尺寸的一种转换方法，本质也是一种卷积，它的梯度计算方法可以参考卷积的方法。
[0159]
c)如图9所示，反池化：以最大池化为例，记录下max
‑
pooling 在对应kernel中的坐标，在反池化过程中，将一个元素根据kernel进行放大，根据之前的坐标将元素填写进去，其它位置补0。
[0160]
它的梯度计算类似于池化层的方法：最大值处的梯度直接传递，其他位置梯度为0。
[0161]
2)依次反向计算梯度：考虑包含一个输入层a、隐藏层b、输出层c的网络，信息传递方向为a
→
b
→
c，c中每一个点c
ij
相对于 a中每一个点a
st
的梯度可通过b传递，写作:
[0162][0163]
考虑整个卷积神经网络，从最末尾的输出层出发，依次计算出相对于前一层、前两层、前三层的梯度，不断向前递归，直至计算得到相对于输入层的梯度。
[0164]
为了进一步地优化本发明的实施效果，在另外一些实施方式中，其余特征技术相同，不同之处在于，fmo具体包括以下步骤：
[0165]
s3.1：将射野分割成笔形束；
[0166]
s3.2：引入基于mc算法的沉积矩阵；
[0167]
s3.3：设置初始强度值，将其作为优化的初始值；
[0168]
s3.4：利用基于mc算法的沉积矩阵求解目标函数的导数；
[0169]
s3.5：求解当前梯度下的最优步长；
[0170]
s3.6：更新强度值、剂量和目标函数；
[0171]
s3.7：重复步骤s3.4
‑
s3.6，直至满足停止条件后，迭代停止。
[0172]
在具体的实施例中，强度分布优化是计划优化两步法的第一步，使用梯度下降法执行此优化的步骤见图4，包括：
[0173]
s3.1：射野分割成笔形束(beamlets)：将设置的射野钨门范围划分成n
×
m个矩形网络，n代表光栅的宽度方向，m代表光栅运动方向，通常情况下n方向的网格按照光栅宽度划分，便于执行下一步的光栅子野分割，m方向按照等间距分割；
[0174]
s3.2：引入基于mc算法的沉积矩阵a；
[0175]
s3.3：设置初始强度值：作为优化的初始值；
[0176]
s3.4：求解目标函数的导数：这里定义的优化目标函数为：
[0177][0178]
[0179]
其中：
[0180]
是目标函数，表示为靶区(ptv)与危及器官(oar)的剂量二次型误差的加权；
[0181]
是危机器官的剂量误差值，α是它的加权；
[0182]
是危机器官的剂量误差值，β是它的加权；
[0183]
是待优化变量,即为所有射野离散成的笔形束的强度通量图；为oar的二次型误差，是所有器官中所有体素点计算剂量和理想剂量之差的二范数；
[0184]
n
oar
为危及器官的个数；
[0185]
n
i
为第i个危机器官的体素点个数；
[0186]
w
j
为第j个体素点的权重；
[0187]
δ
j
为惩罚因子，当计算剂量满足预设的约束时，即为0，否则记为1；
[0188]
为根据射野强度通量图计算的j点剂量值，是计算剂量；
[0189]
p
j
为理想剂量值；
[0190]
同理，为靶区的二次型误差，其中n
ptv
为靶区的体素点个数；
[0191]
a
jm
为第m个笔形束对第j个体素的剂量贡献，组成的矩阵就是基于mc算法的沉积矩阵a。
[0192]
因此，目标函数的导数为：
[0193][0194]
s3.5：求解最优步长：求解当前梯度下的最优步长；
[0195]
s3.6：更新强度值、剂量和目标函数；
[0196]
s3.7：判断是否满足停止条件：当自变量(强度值x)或目标函数j
obj
(x)的绝对变化值或相对变化值小于阈值ε时，迭代停止，否则重复步骤s3.4
‑
s3.4。
[0197]
为了进一步地优化本发明的实施效果，在另外一些实施方式中，其余特征技术相同，不同之处在于，该放疗计划强度分布优化方法还包括以下步骤：
[0198]
s4：根据放疗计划强度分布优化方法计算人体内的剂量和各个器官的dvh，评估计划效果。
[0199]
在具体的实施例中，直接使用gdpm计算所有笔形束(beamlets) 的沉积矩阵(总粒子数约为2
×
109)，将其代入优化流程，得到最优强度分布，进一步地获得剂量分布，作为强度分布优化的基准(baseline)。
[0200]
使用本发明公开的方法得到基于mc算法的沉积矩阵，加入优化流程中去，得到优化强度分布，使用mc算法计算人体内的剂量，并统计出器官的dvh。
[0201]
对比两者的dvh，验证本发明的计划效果。
[0202]
两者的结果图可见图12，图12中：neural dose代表基于神经网络的结果(即本发明的结果)，mc剂量代表基于随机计算的结果，可见两者结果几乎完全一致，验证了本发明的可行性。
[0203]
另一方面，本发明实施例还公开一种放疗计划强度分布优化装置，包括：
[0204]
神经网络映射模块，用于通过神经网络，将pb剂量映射到mc 剂量上；
[0205]
沉积矩阵获取模块，用于求解mc剂量相对于pb剂量的梯度，得到基于mc算法的沉积矩阵；
[0206]
强度分布优化模块，用于将基于mc算法的沉积矩阵带入fmo 中，实现放疗计划强度分布优化。
[0207]
另一方面，本发明实施例还公开放射治疗设备，用于执行上述任一实施例公开的放疗计划强度分布优化方法，或，包括：上述实施例公开的放疗计划强度分布优化装置。
[0208]
本发明一种放疗计划强度分布优化方法、装置及放射治疗设备具有以下有益效果：
[0209]
本发明将mc算法应用在fmo中可以确保优化过程得到的强度分布是真实可靠的，使用这样的强度分布进行多叶光栅子野分割，得到的子野走位是更加合理的，保证了最终得到的计划剂量分布更接近于临床要求，计划效果有显著提升。
[0210]
上述实施例只为说明本发明的技术构思及特点，其目的在于让本领域普通技术人员能够了解本发明的内容并加以实施，并不能以此限制本发明的保护范围，凡根据本发明精神实质所作的等效变化或修饰，都应涵盖在本发明的保护范围内。