一种基于动态半影模型的动态调强装置及优化方法制造方法及图纸

本发明专利技术公开了一种基于动态半影模型的动态调强装置及优化方法。涉及医学物理领域，具体优化方法的步骤如下：首先采用叶片速度时序函数，在叶片允许的最大速度和最大加速度的前提下，建立不同时刻下的叶片速度的时序模型；其次，在动态调强优化过程中，基于已建立的不同时刻下的叶片速度的时序模型以及动态半影模型，采用直接机器参数优化方法快速优化生成动态调强计划。本发明专利技术能够快速优化得到满足临床要求的动态调强计划。床要求的动态调强计划。床要求的动态调强计划。

【技术实现步骤摘要】
一种基于动态半影模型的动态调强装置及优化方法


[0001]本专利技术属于医学物理领域，涉及了一种用于动态调强优化问题模型，具体的，是涉及一种基于动态半影模型的动态调强装置及优化方法。

技术介绍

[0002]随着放射治疗技术的蓬勃发展，其放射治疗计划的装置和方法也在相继提出并付诸使用，尤其近些年快速发展的动态调强技术。
[0003]动态调强即在射线照射时多叶准直器叶片处于运动状态，利用多叶准直器的相对应的每对叶片在射线照射时在计算机控制下对靶区进行扫描以实现对射野强度的调节。动态模式依靠对每个叶片运动的速度和照射的剂量率的调整来实现强度的调试。对于动态调强，由于在治疗时，叶片在运动过程中治疗机一直在出束，所以要想缩短总的治疗实施时间就应该尽量缩短治疗出束的时间；对于目前的动态调强，存在以下缺陷：第一，照射的时间较长，效率较低；第二，没有考虑叶片半影，造成实际的精度有误差。

技术实现思路

[0004]专利技术目的：本专利技术的目的是提供了一种基于动态半影模型的动态调强装置及优化方法，建立基于不同时刻下的叶片速度的时序模型，并基于此实施直接机器参数优化模型，实现动态调强计划的快速、精确优化。
[0005]技术方案：本专利技术所述的一种基于动态半影模型的动态调强装置，包括用户输入装置，叶片关于速度的时序模型装置，动态调强优化装置及用户输出装置；
[0006]进一步的，所述用户输入装置包括用户输入器官组织、叶片最大运动速度、叶片最大加速度、剂量变化率、机架角及计划优化参数等基本参数信息；
[0007]进一步的，所述叶片关于速度的时序模型装置包括使用叶片速度
‑
时间函数在不同时刻上建立叶片速度的时序模型；具体的，在射线源照射过程中，基于光栅叶片的相对运动实现剂量调节，在叶片允许的最大速度和最大加速度的前提下，考虑叶片的非匀速运动，建立不同时刻下的叶片速度的时序模型；
[0008]进一步的，所述动态调强优化装置包括通过叶片速度的时序模型得到目标函数对叶片运动时间的偏导数，基于建立的叶片速度的时序模型以及位置的可变半影模型，采用直接机器参数优化生成动态调强计划；
[0009]进一步的，所述用户输出装置，包括叶片通过优化装置得到的叶片运动时间以及所对应的运动位置、叶片的运动速度及叶片的加速度等信息；
[0010]本专利技术还提出一种基于动态半影模型的动态调强优化方法，由叶片速度的时序模型建立模块，动态调强优化模块共两个模块组成；
[0011]其中，所述叶片速度
‑
时间建立模块由机器参数作为输入参数，采用叶片可变速度时间函数建立叶片速度的时序模型；
[0012]所述动态调强优化模块，基于建立的叶片速度的时序模型，采用基于半影模型的
直接机器参数优化生成动态调强计划。
[0013]进一步的，所述叶片速度
‑
时间建立模块包括：
[0014]使用叶片速度
‑
时间函数在不同时刻上建立叶片速度的时序模型：在射线源照射过程中，基于光栅叶片的相对运动实现剂量调节，在叶片允许的最大速度和最大加速度的前提下，考虑叶片的非匀速运动，建立不同时刻下的叶片速度的时序模型，计算公式如下所示：
[0015]式(1)表示叶片位置x与叶片加速度a(t)以及速度v(t)之间的关系：
[0016]x＝γ(a(t),v(t))
ꢀꢀꢀ
(1)
[0017]其中，a(t)表示叶片t时刻点的加速度，v(t)表示叶片t时刻点的速度，γ表示x与叶片加速度a(t)以及速度v(t)之间的映射关系。
[0018]式(2)表示叶片t时刻点的加速度约束：
[0019]|a(t)|<a
max
ꢀꢀꢀ
(2)
[0020]其中，a
max
是叶片运动的最大加速度。
[0021]式(3)表示的是叶片t时刻点的速度约束：
[0022]|v(t)|<v
max
ꢀꢀꢀ
(3)
[0023]其中，v
max
是叶片运动的最大速度。
[0024]进一步的，所述动态调强优化模块的执行步骤具体如下：
[0025]步骤1：用户输入叶片参数，构建叶片模型，使用蒙特卡罗方法，模拟不同位置下的半影参数：从叶片运动的初始化的位置到最终位置之间进行离散化处理，对于其中的每一个位置进行蒙特卡罗方法模拟，用此方式，可以得到不同位置下的半影参数；
[0026]步骤2：根据不同位置的半影参数，建立基于位置的可变半影模型，半影模型即考虑到不同位置的射束的散射，修正位置的偏差；
[0027]步骤3：使用可变半影模型的直接机器参数优化方法，优化子野形状；当满足结束条件时，则结束，否则继续使用可变半影模型的直接机器参数装置和方法，优化子野形状；
[0028]进一步的，在步骤2中，所述建立基于位置的可变半影模型是：针对不同位置的半影参数，在叶片和放射源的不同相对位置分别进行计算，当放射源被叶片遮挡的程度不同的时候，形成一组不同的强度曲线，通过数据拟合，可以得到可变半影模型；
[0029]所述步骤3中，直接机器参数优化方法，使用直接机器参数优化调整叶片位置，利用目标函数对时间的偏导数，采用共轭梯度法计算目标函数极小值，目标函数计算公式如下：
[0030]F＝F(x)＝min∑
i
[∑
j
M
xij
*M
j
(x)
‑
D
i
]2ꢀꢀꢀ
(4)
[0031][0032]其中，v(t),a(t)分别来自于式(2)与式(3)，F为目标函数值，x为叶片位置，1≤i≤S，1≤j≤N，其中S与N均为正整数，S表示器官点总数，N表示粒子束总数，M
xij
是当叶片位置为x时第j个粒子束对第i个器官点的剂量贡献，D
i
为第i个器官点的处方剂量，M
j
(x)是第j个粒子束穿过叶片x后的剂量强度分布，t为叶片运动时间，是关于叶片在位置x处速度、加速度与时间的函数。
[0033]对上述公式目标函数F求偏导数，通过叶片速度的时序模型得到目标函数对光栅时间偏导数，偏导数计算公式如下：
[0034][0035]其中，来自于式(1)，F为目标函数值，t为叶片运动时间，x为叶片位置，D为器官点的剂量贡献矩阵，M是剂量强度分布。
[0036]使用直接机器参数优化调整叶片位置，利用目标函数对时间的偏导数，采用共轭梯度法计算目标函数极小值，快速、精确优化动态调强计划。
[0037]有益效果：本专利技术与现有技术相比，本专利技术通过叶片速度
‑
时间参数，建立叶片速度的时序模型，初创性的考虑了基于位置的可变半影模型，采用直接机器参数约束优化，能够快速交付精确的动态调强计划。
附图说明<本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于动态半影模型的动态调强装置，其特征在于，包括用户输入装置，叶片关于速度的时序模型装置，动态调强优化装置及用户输出装置；所述用户输入装置包括用户输入器官组织、叶片最大运动速度、叶片最大加速度、剂量变化率、机架角及计划优化参数；所述叶片关于速度的时序模型装置包括使用叶片速度
‑
时间函数建立叶片速度的时序模型；所述动态调强优化装置包括通过叶片速度的时序模型得到目标函数对叶片运动时间的偏导数；所述用户输出装置包括叶片通过优化装置得到的叶片运动时间以及所对应的运动位置、叶片的运动速度及叶片的加速度信息。2.根据权利要求1所述的一种基于动态半影模型的动态调强装置的优化方法，其特征在于，包括模型建立模块及动态调强优化模块。3.根据权利要求2所述的一种基于动态半影模型的动态调强装置的优化方法，其特征在于，所述模型建立模块，是根据优化的参数信息，使用叶片速度
‑
时间函数在不同时刻上建立叶片速度的时序模型：具体的，在射线源照射过程中，基于光栅叶片的相对运动实现剂量调节，在叶片允许的最大速度和最大加速度的前提下，考虑叶片的非匀速运动，建立不同时刻下的叶片速度的时序模型，计算公式如下所示：式(1)表示叶片位置x与叶片加速度a(t)以及速度v(t)之间的关系：x＝γ(a(t),v(t))
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)其中，a(t)表示叶片t时刻点的加速度，v(t)表示叶片t时刻点的速度，γ表示叶片位置x与叶片加速度a(t)、速度v(t)之间的映射关系；式(2)表示叶片t时刻点的加速度约束：|a(t)|&lt;a
max
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(2)其中，a
max
表示叶片运动的最大加速度；式(3)表示叶片t时刻点的速度约束：|v(t)|&lt;v
max
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(3)其中，v
max
表示叶片运动的最大速度。4.根据权利要求2所述的一种基于动态半影模型的动态调强装置的优化方法，其特征在于，所述动...

【专利技术属性】
技术研发人员：ꢀ七四专利代理机构，
申请(专利权)人：中科超精南京科技有限公司，
类型：发明
国别省市：