CN103083821B - 一种直接优化静态调强子野形状及其机器跳数的***和方法 - Google Patents
一种直接优化静态调强子野形状及其机器跳数的***和方法 Download PDFInfo
- Publication number
- CN103083821B CN103083821B CN201310029736.7A CN201310029736A CN103083821B CN 103083821 B CN103083821 B CN 103083821B CN 201310029736 A CN201310029736 A CN 201310029736A CN 103083821 B CN103083821 B CN 103083821B
- Authority
- CN
- China
- Prior art keywords
- ziye
- dose
- optimization
- launched field
- shape
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Withdrawn - After Issue
Links
Landscapes
- Radiation-Therapy Devices (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种直接优化静态调强子野形状及其机器跳数的***和方法。该***包括输入模块、数据处理模块、优化模块及优化结果输出模块,本该发明主要是根据放疗计划制定者要求,自动建立相应的目标函数,然后采用遗传算法进行子野形状的调整,采用共轭梯度算法进行子野权重的优化,两种算法交替优化,最后得到满足要求子野形状及其对应机器跳数权重。该方法可以快速的得到满足要求的每个射野方向的各个子野形状及其机器跳数,优化结果就是直接可以实施的子野序列。本发明不仅可以有效的解决传统方法存在的精度的问题,而且操作简单,可以提高制定计划和实施计划效率。
Description
技术领域
本发明涉及到精确放疗中,采用多叶光栅进行调强时,最小化多叶光栅对应的子野个数及总的机器跳数的***和方法。
背景技术
静态调强放疗中得到可以实施的子野形状及其机器跳数的一般步骤是首先借助优化算法优化得到每个方向上满足期望剂量分布的最优强度,然后对优化后的强度进行子野分割,转换为MLC可以实施的子野,这种方式简称为“二步法”。由于在强度分布优化时,并没有考虑多叶光栅本身所包含的一些约束限制,这些限制只有在子野分割时考虑,因此导致子野分割会产生很多个子野形状和及其对应的权重。长期的临床应用表明,“二步法”存在的缺陷,可以总结为:(1)每个照射野的子野数目多,总的照射跳数(MU)往往是常规/适形照射野的数倍,这种情况不仅会增加照射的时间,还将增加产生漏射、机头散射和tongue-groove效应的概率;(2)可能会生成一些小的、形状很不规则的子野,这些子野在实施过程中产生的误差较大,射线利用率低,因此可能影响患者受照剂量的准确度;(3)在将理想的强度分布转换为可以实施的强度分布时,两个强度会有差别,这种差别往往使计划质量变差,本已满足的临床要求的计划可能不再满足。
发明内容
本发明技术解决问题:为了克服“二步法”的缺陷,提供一种根据期望的剂量分布,直接优化每个射野方向上各个子野的形状及其机器跳数的***和方法,可以快速得到满足要求的每个射野方向的各个子野形状及其机器跳数,不仅可以有效的解决传统“二步法”存在的精度问题,而且操作简单,还提高了计划制定和计划实施的效率。
本发明的技术方案是:一种直接优化静态调强子野形状及其机器跳数的***,如图1所示,包括:
(1)数据导入模块,用于导入从放疗计划***中输出的人体影像CT数据、靶区及感兴趣区域定义信息、射野方向信息、多叶光栅信息及计划设计者对计划的要求信息;
(2)数据处理模块,根据数据导入模块导入的靶区信息,射野方向信息,计算靶区在射野方向上的投影,向外扩一定的大小得到对应射野方向的射野形状,然后根据多叶光栅类型,对每个射野形状进行网格化,在多叶光栅运动方向上网格大小为0.5cm,在多叶光栅运动垂直方向,网格大小为多叶光栅叶片在等中心平面厚度,最后采用基于蒙卡的快速有限笔形束剂量计算方法计算射野内每个单元野对各个计算采样点的剂量影响;
(3)优化模块,主要是进行优化模型的建立及采用优化算法进行优化;根据导入模块导入的计划设计者对计划的要求,通过二次目标函数形式转化为优化的目标函数,然后供步骤优化过程使用;采用基于物理剂量或剂量-体积约束的目标函数及约束,表达为:
此时,xm k为待优化的各子野的权重,aim为第m个子野对第i个采样点的剂量影响,由剂量计算方法计算得到;第一个式子为优化的目标函数,优化的目标是使该目标函数值最小,为靶区对目标函数的贡献,为第j个危及器官对目标函数的贡献;wPTV和分别是靶区和第j个危及器官的权重归一化因子;式(2)中di是通过(3)式计算出来的第i个计算采样点的剂量,DPTV是靶区的处方剂量,是第j个危及器官的剂量限值,nPTV和分别为靶区和第j个危及器官的计算采样点数,NOAR是优化中考虑的危及器官的个数;式(3)中,Nray为子野的总个数,aim为单位强度的第m个子野对第i个采样点的剂量贡献,xm k为向量的第m个分量即第m个子野的权重,因此第i个采样点的剂量是所有子野对该点剂量贡献的总和;所有的aim构成的矩阵事先由剂量计算模块计算好供优化计算时直接调用;xm k即为优化时待调整的参数,优化的目标即为找到最优的xm k使得各点的剂量满足医生给定的剂量要求,其中K为当前迭代的次数;
根据数据处理模块中得到的射野形状对应的多叶光栅位置来初始化遗传算法群体,遗传算法通过交叉变异得到新的个体,即每个个体对应一组子野形状,采用共轭梯度算法进行每个子野的权重优化,共轭梯度算法优化的目标函数为根据计划设计者的要求建立的优化模型,如公式(1)~(3)表达式,然后评估遗传算法群体中每个个体,计算目标函数的值,根据目标函数值采用联赛选择方式选择好的个体进入下一代进化过程中;重复上述步骤直到遗传算法群体最优个体的目标函数值满足要求,终止优化计算;
(4)数据输出模块,将遗传算法优化得到的最优个体通过解码得到对应的最优的子野形状,最后输出子野形状、机器跳数及根据子野形状及机器跳数计算出的剂量分布结果。
一种直接优化静态调强子野形状及其机器跳数的方法,实现步骤如下:
(1)按照Dicom-RT格式导入人体CT影像数据、靶区及感兴趣区域定义信息、射野方向信息,根据靶区在射野方向上的投影并向外扩一定的大小作为对应射野方向的射野形状,然后供步骤(4)使用;根据CT影像数据、靶区及感兴趣区域定义信息对各区域进行采样得到各区域的计算点信息供步骤(4)和步骤(5)使用;
(2)导入每个方向的子野个数、采用的多叶光栅(Multi-LeafCollimator,MLC)类型供步骤(5)中遗传算法初始化个体时使用,多叶光栅类型信息还供步骤(4)中射野网格化使用;
(3)导入计划设计者对计划的要求,通过二次目标函数形式转化为优化的目标函数,然后供步骤(5)使用;采用基于物理剂量或剂量-体积约束的目标函数及约束,表达为:
此时,xm k为待优化的各子野的权重,aim为第m个子野对第i个采样点的剂量影响,由剂量计算方法计算得到;第一个式子为优化的目标函数,优化的目标是使该目标函数值最小,为靶区对目标函数的贡献,为第j个危及器官对目标函数的贡献;wPTV和分别是靶区和第j个危及器官的权重归一化因子;式(5)中di是通过(6)式计算出来的第i个计算采样点的剂量,DPTV是靶区的处方剂量,是第j个危及器官的剂量限值,nPTV和分别为靶区和第j个危及器官的计算采样点数,NOAR是优化中考虑的危及器官的个数;式(6)中,Nray为子野的总个数,aim为单位强度的第m个子野对第i个采样点的剂量贡献,xm k为向量的第m个分量即第m个子野的权重,因此第i个采样点的剂量是所有子野对该点剂量贡献的总和;所有的aim构成的矩阵事先由剂量计算模块计算好供优化计算时直接调用;xm k即为优化时待调整的参数,优化的目标即为找到最优的xm k使得各点的剂量满足医生给定的剂量要求,其中K为当前迭代的次数;
(4)根据步骤(1)中的射野形状及步骤(2)中的多叶光栅类型,对每个射野进行网格化,在多叶光栅运动方向上网格大小为0.5cm,在多叶光栅运动垂直方向,网格大小为多叶光栅叶片在等中心平面厚度,然后采用基于蒙卡的快速有限笔形束剂量计算方法计算射野内每个单元野,即笔形束对各个计算采用点的剂量影响;此步骤是根据步骤(1)中射野方向信息将射野划分成网格,然后每个网格即为本步骤中的每个单元野,即笔形束;
(5)根据步骤(1)中得到的射野形状对应的多叶光栅位置来初始化遗传算法群体,遗传算法通过交叉变异得到新的个体,即每个个体对应一组子野形状,采用共轭梯度算法进行每个子野的权重优化,共轭梯度算法优化步骤(3)中的目标函数,评估遗传算法群体中每个个体,计算目标函数的值,根据目标函数值采用联赛选择方式选择好的个体进入下一代进化过程中;重复上述步骤直到遗传算法群体最优个体的目标函数值满足要求,终止优化计算;
(6)最后输出最优的子野形状、机器跳数及剂量分布结果。
所述步骤(5)中遗传算法群体初始化,每个个体对应的子野序列中每个子野初始形状为根据靶区当前射野方向上投影外扩0.1cm~1cm对应的MLC形状。
所述遗传算法的编码方式采用二维二进制编码,且每行中的编码为1的位置连续;交叉和变异采用常规的二进制编码方式下的交叉变异方式,不满足每行中1的位置连续的交叉变异被放弃。
所述共轭梯度算法在迭代过程中,迭代次数小于20,当迭代出现负值时,直接将其赋值为0,并继续优化。
本发明与现有技术相比的优点在于:
(1)本发明通过优化可以直接得到满足计划设计者要求的子野形状及其机器跳数,从而提高了优化后计划的可实施性;
(2)本发明采用遗传算法和共轭梯度算法相结合的方法进行优化,利用两个算法的优点,优化速度快,比原来的两步法得到最终结果的速度及效果好。
附图说明
图1为本发明***的组成框图;
图2为本发明方法实现流程图。
具体实施方式
如图2所示,本发明一种直接优化静态调强子野形状及其机器跳数的方法,具体实现如下:
(1)按照Dicom-RT格式导入CT影像数据、靶区及感兴趣区域定义信息、射野方向信息,根据靶区在射野方向上的投影并向外扩一定的大小作为对应射野方向的射野形状,根据CT影像数据靶区及感兴趣区域定义信息对计算点进行采样得到各区域的计算点信息;
比如,针对鼻炎测试模型,导入57张CT影像数据及对应各张勾画好的感兴趣区域(左右眼睛、脊髓、脑干和左右腮腺)及靶区信息(GTV、PTV)、射野方向信息采用常规的360度范围内均匀分布9个野的方式。
(2)计划设计者设定每个方向上的子野个数为7、采用的多叶光栅(Multi-LeafCollimatorMLC)类型为瓦里安40对多叶光栅,本***将多叶光栅的叶片对数、每个叶片的厚度及约束限制输入给优化模块使用;
(3)导入计划设计者对计划的要求,针对鼻炎测试模型,靶区PTV的处方剂量为6000cGy,GTV的处方剂量为7000cGy,脊髓最大剂量为4500cGy,脑干最大剂量不大于5400cGy,左右腮腺超过35cGy的体积不大于50%,左右眼睛的最大剂量为1000cGy。并将其通过二次目标函数形式转化为优化的目标函数,详见权利要求1中的步骤3,然后供步骤5使用;
(4)根据步骤1中的射野形状及步骤2中的多叶光栅类型,对每个射野进行网格化,在多叶光栅运动方向上网格大小为0.5cm,在多叶光栅运动垂直方向,网格大小为多叶光栅叶片在等中心平面厚度,然后采用基于有限笔形束剂量计算方法计算射野内每个单元野对各个计算采用点的剂量影响;根据步骤1中得到的射野形状对应的多叶光栅位置来初始化遗传算法群体,遗传算法通过交叉变异得到新的个体(每个个体对应一组子野形状),采用共轭梯度算法进行每个子野的权重优化,共轭梯度算法优化步骤3中的目标函数,评估遗传算法群体中每个个体,计算目标函数的值,根据目标函数值采用联赛选择方式选择好的个体进入下一代进化过程中。重复上述步骤直到遗传算法群体最优个体的目标函数值满足要求,终止优化计算;最后输出最优的子野形状、机器跳数及剂量分布结果。
表1为针对同一个测试模型,设置相同的射野大小及方向、靶区和各器官的约束剂量或剂量体积约束,用本发明确定出的子野个数及其子野权重,与传统“二步法”确定出的子野个数及权重进行对比。对本例检验表明,得到几乎同样的剂量分布情况下,本发明计算得到的子野个数及总的MU分别为63和2833,而传统“两步法”得到的子野个数和总的MU分别为327和4369。可以看出本发明比传统方法子野个数减少了81%,MU减少了35%。如果两个计划均采用剂量率为400MU/min的加速器照射,那么“二步法”与本发明对应的机器出束时间分别为10.92min和7.08min。如果每个子野间的切换时间为7s,那么“二步法”切换子野形状所需要的时间总共为38.15min,而本发明切换子野只需要7.35min,因此本发明将大大缩短实施时间。
表1常规两步法与直接子野所用子野个数和MU对比表
总之,本发明可以快速的得到满足要求的每个射野方向的各个子野形状及其机器跳数,优化结果就是直接可以实施的子野序列。因此,不仅可以有效的解决传统方法存在的精度的问题,而且操作简单,可以提高制定计划和实施计划效率。
本发明未详细阐述部分属于本领域公知技术。
以上所述,仅为本发明部分具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本领域的人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。
Claims (5)
1.一种直接优化静态调强子野形状及其机器跳数的***,其特征在于包括:
(1)数据导入模块,用于导入从放疗计划***中输出的人体影像CT数据、靶区及感兴趣区域定义信息、射野方向信息、多叶光栅信息及计划设计者对计划的要求信息;
(2)数据处理模块,根据数据导入模块导入的靶区信息,射野方向信息,计算靶区在射野方向上的投影,向外扩一定的大小得到对应射野方向的射野形状,然后根据多叶光栅类型,对每个射野形状进行网格化,在多叶光栅运动方向上网格大小为0.5cm,在多叶光栅运动垂直方向,网格大小为多叶光栅叶片在等中心平面厚度,最后采用基于蒙卡的快速有限笔形束剂量计算方法计算射野内每个单元野对各个计算采样点的剂量影响;
(3)优化模块,根据计划设计者对计划的要求进行优化模型的建立及采用优化算法进行优化;根据导入模块导入的计划设计者对计划的要求,通过二次目标函数形式转化为优化的目标函数,然后供步骤优化过程使用;采用基于物理剂量或剂量-体积约束的目标函数及约束,表达为:
此时,xm k为待优化的各子野的权重,aim为第m个子野对第i个采样点的剂量影响,由剂量计算方法计算得到;式(1)为优化的目标函数,优化的目标是使该目标函数的值最小fPTV(xk)为靶区内各个采样点剂量与期望剂量差别的平均值,为第j个危及器官各个采样点剂量与期望剂量差别的平均值;wPTV和分别是靶区和第j个危及器官的权重归一化因子;式(2)中di是通过(3)式计算出来的第i个计算采样点的剂量,DPTV是靶区的处方剂量,是第j个危及器官的剂量限值,nPTV和分别为靶区和第j个危及器官的计算采样点数,NOAR是优化中考虑的危及器官的个数;式(3)中,Nray为子野的总个数,aim为单位强度的第m个子野对第i个采样点的剂量贡献,xm k为xk向量的第m个分量即第m个子野的权重,因此第i个采样点的剂量是所有子野对该点剂量贡献的总和;所有的aim构成的矩阵事先由剂量计算模块计算好供优化计算时直接调用;xm k即为优化时待调整的参数,优化的目标即为找到最优的xm k使得各点的剂量满足医生给定的剂量要求,其中K为当前迭代的次数;
根据数据处理模块中得到的射野形状对应的多叶光栅位置来初始化遗传算法群体,遗传算法通过交叉变异得到新的个体,即每个个体对应一组子野形状,采用共轭梯度算法进行每个子野的权重优化,共轭梯度算法优化的目标函数为根据计划设计者的要求建立的优化模型,如公式(1)~(3)表达式,然后评估遗传算法群体中每个个体,计算目标函数的值,根据目标函数值采用联赛选择方式选择好的个体进入下一代进化过程中;重复上述步骤直到遗传算法群体最优个体的目标函数值满足要求,终止优化计算;
(4)数据输出模块,将遗传算法优化得到的最优个体通过解码得到对应的最优的子野形状,最后输出子野形状、机器跳数及根据子野形状及机器跳数计算出的剂量分布结果。
2.一种直接优化静态调强子野形状及其机器跳数的方法,其特征在于实现步骤如下:
(1)按照Dicom-RT格式导入人体CT影像数据、靶区及感兴趣区域定义信息、射野方向信息,根据靶区在射野方向上的投影并向外扩一定的大小作为对应射野方向的射野形状,然后供步骤(4)使用;根据CT影像数据、靶区及感兴趣区域定义信息对各区域进行采样得到各区域的计算点信息供步骤(4)和步骤(5)使用;
(2)导入每个方向的子野个数、采用的多叶光栅(Multi-LeafCollimator,MLC)类型供步骤(5)中遗传算法初始化个体时使用,多叶光栅类型信息还供步骤(4)中射野网格化使用;
(3)导入计划设计者对计划的要求,通过二次目标函数形式转化为优化的目标函数,然后供步骤(5)使用;采用基于物理剂量或剂量-体积约束的目标函数及约束,表达为:
此时,xm k为待优化的各子野的权重,aim为第m个子野对第i个采样点的剂量影响,由剂量计算方法计算得到;公式(4)为优化的目标函数,优化的目标是使该目标函数的值最小,fPTV(xk)为靶区对目标函数的贡献,为第j个危及器官对目标函数的贡献;wPTV和分别是靶区和第j个危及器官的权重归一化因子;式(4)中di是通过(5)式计算出来的第i个计算采样点的剂量,DPTV是靶区的处方剂量,是第j个危及器官的剂量限值,nPTV和分别为靶区和第j个危及器官的计算采样点数,NOAR是优化中考虑的危及器官的个数;式(6)中,Nray为子野的总个数,aim为单位强度的第m个子野对第i个采样点的剂量贡献,xm k为xk向量的第m个分量即第m个子野的权重,因此第i个采样点的剂量是所有子野对该点剂量贡献的总和;所有的aim构成的矩阵事先由剂量计算模块计算好供优化计算时直接调用;xm k即为优化时待调整的参数,优化的目标即为找到最优的xm k使得各点的剂量满足医生给定的剂量要求,其中K为当前迭代的次数;
(4)根据步骤(1)中的射野形状及步骤(2)中的多叶光栅类型,对每个射野进行网格化,在多叶光栅运动方向上网格大小为0.5cm,在多叶光栅运动垂直方向,网格大小为多叶光栅叶片在等中心平面厚度,然后采用基于蒙卡的快速有限笔形束剂量计算方法计算射野内每个单元野,即笔形束对各个计算采用点的剂量影响;此步骤是根据步骤(1)中射野方向信息将射野划分成网格,然后每个网格即为本步骤中的每个单元野,即笔形束;
(5)根据步骤(1)中得到的射野形状对应的多叶光栅位置来初始化遗传算法群体,遗传算法通过交叉变异得到新的个体,即每个个体对应一组子野形状,采用共轭梯度算法进行每个子野的权重优化,共轭梯度算法优化步骤(3)中的目标函数,评估遗传算法群体中每个个体,计算目标函数的值,根据目标函数值采用联赛选择方式选择好的个体进入下一代进化过程中;重复上述步骤直到遗传算法群体最优个体的目标函数值满足要求,终止优化计算;
(6)最后输出最优的子野形状、机器跳数及剂量分布结果。
3.根据权利要求2所述的直接优化静态调强子野形状及其机器跳数的方法,其特征在于:所述步骤(5)中遗传算法群体初始化,每个个体对应的子野序列中每个子野初始形状为根据靶区当前射野方向上投影外扩0.1cm~1cm对应的MLC形状。
4.根据权利要求2所述的直接优化静态调强子野形状及其机器跳数的方法,其特征在于:所述遗传算法的编码方式采用二维二进制编码,且每行中的编码为1的位置连续;交叉和变异采用常规的二进制编码方式下的交叉变异方式,不满足每行中1的位置连续的交叉变异被放弃。
5.根据权利要求2所述的直接优化静态调强子野形状及其机器跳数的方法,其特征在于:所述共轭梯度算法在迭代过程中,迭代次数小于20,当迭代出现负值时,直接将其赋值为0,并继续优化。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310029736.7A CN103083821B (zh) | 2013-01-27 | 2013-01-27 | 一种直接优化静态调强子野形状及其机器跳数的***和方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310029736.7A CN103083821B (zh) | 2013-01-27 | 2013-01-27 | 一种直接优化静态调强子野形状及其机器跳数的***和方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103083821A CN103083821A (zh) | 2013-05-08 |
CN103083821B true CN103083821B (zh) | 2015-12-23 |
Family
ID=48197274
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310029736.7A Withdrawn - After Issue CN103083821B (zh) | 2013-01-27 | 2013-01-27 | 一种直接优化静态调强子野形状及其机器跳数的***和方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103083821B (zh) |
Families Citing this family (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104636832B (zh) * | 2015-02-27 | 2018-09-21 | 上海联影医疗科技有限公司 | 剂量分布变化估算方法及装置、直接子野优化方法及*** |
CN105031820B (zh) * | 2015-09-01 | 2018-04-03 | 上海联影医疗科技有限公司 | 调强放疗逆向优化方法和装置 |
CN105617536A (zh) * | 2015-12-24 | 2016-06-01 | 上海联影医疗科技有限公司 | 旋转逆向调强放疗优化方法及装置 |
CN105709341B (zh) * | 2016-01-15 | 2018-03-27 | 中国科学院合肥物质科学研究院 | 一种基于梯度法和漫水填充法的调强子野优化方法 |
US10143859B2 (en) * | 2016-09-22 | 2018-12-04 | Varian Medical Systems International Ag | Controlling and shaping the dose distribution outside treatment targets in external-beam radiation treatments |
CN108066902B (zh) * | 2016-11-14 | 2019-12-20 | 上海东软医疗科技有限公司 | 一种子野信息确定方法和装置 |
CN106730411B (zh) * | 2016-12-21 | 2019-11-19 | 上海联影医疗科技有限公司 | 放射治疗计划优化方法及*** |
CN107297033B (zh) * | 2017-07-13 | 2019-02-26 | 中科超精(安徽)科技有限公司 | 冻结笔形束的放射治疗笔形束权重优化加速方法和*** |
CN107998519A (zh) * | 2017-12-26 | 2018-05-08 | 荣科科技股份有限公司 | 一种用于imrt的分子动力学再优化算法 |
CN108211136B (zh) * | 2017-12-29 | 2021-01-15 | 苏州雷泰医疗科技有限公司 | 基于动态光栅的旋转调强优化方法及放射治疗设备 |
CN108175953B (zh) * | 2017-12-29 | 2020-07-21 | 苏州雷泰医疗科技有限公司 | 一种基于动态光栅的旋转调强优化方法及放射治疗设备 |
CN110246562B (zh) * | 2018-03-07 | 2021-08-06 | 有方(合肥)医疗科技有限公司 | 确定放疗***中的子射束强度的方法、装置及计算机*** |
CN110237439B (zh) * | 2018-12-21 | 2020-12-18 | 苏州雷泰医疗科技有限公司 | 用于正交双层光栅装置的静态调强子野优化方法 |
CN110020710B (zh) * | 2019-03-08 | 2023-10-27 | 华南理工大学 | 一种基于人工蜂群算法的射束方向及权重多目标优化方法 |
CN110211664B (zh) * | 2019-04-25 | 2022-11-04 | 安徽大学 | 一种基于机器学习自动设计放射治疗方案的*** |
CN110215623B (zh) * | 2019-06-11 | 2020-11-10 | 苏州雷泰智能科技有限公司 | 基于子野优化的正交双层光栅旋转调强的实现方法及装置 |
CN111986778B (zh) * | 2020-07-31 | 2024-02-20 | 上海联影医疗科技股份有限公司 | 一种调强计划优化***、装置及存储介质 |
CN112999531B (zh) * | 2021-02-25 | 2022-06-24 | 中科超精(南京)科技有限公司 | 一种基于呼吸运动的动态优化*** |
CN112926156B (zh) * | 2021-02-25 | 2024-04-26 | 中科超精(南京)科技有限公司 | 一种基于位置可变半影模型的直接子野优化装置及其优化方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101422640A (zh) * | 2008-11-25 | 2009-05-06 | 中国科学院等离子体物理研究所 | 优化放射治疗射束强度分布的多目标优化方法和*** |
CN102136041A (zh) * | 2011-04-18 | 2011-07-27 | 深圳市海博科技有限公司 | 治疗计划逆向规划方法和治疗计划*** |
CN102247660A (zh) * | 2011-04-18 | 2011-11-23 | 深圳市海博科技有限公司 | 治疗计划逆向规划方法和治疗计划*** |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7529339B2 (en) * | 2003-12-12 | 2009-05-05 | University Of Western Ontario | Method and system for optimizing dose delivery of radiation |
-
2013
- 2013-01-27 CN CN201310029736.7A patent/CN103083821B/zh not_active Withdrawn - After Issue
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101422640A (zh) * | 2008-11-25 | 2009-05-06 | 中国科学院等离子体物理研究所 | 优化放射治疗射束强度分布的多目标优化方法和*** |
CN102136041A (zh) * | 2011-04-18 | 2011-07-27 | 深圳市海博科技有限公司 | 治疗计划逆向规划方法和治疗计划*** |
CN102247660A (zh) * | 2011-04-18 | 2011-11-23 | 深圳市海博科技有限公司 | 治疗计划逆向规划方法和治疗计划*** |
Also Published As
Publication number | Publication date |
---|---|
CN103083821A (zh) | 2013-05-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103083821B (zh) | 一种直接优化静态调强子野形状及其机器跳数的***和方法 | |
US9067064B2 (en) | Optimization process for volumetric modulated arc therapy | |
CN102921115B (zh) | 一种建立基于测量数据的医用直线加速器简便照射源模型的方法 | |
CN110237439B (zh) | 用于正交双层光栅装置的静态调强子野优化方法 | |
Tian et al. | A large-scale combinatorial many-objective evolutionary algorithm for intensity-modulated radiotherapy planning | |
CN107551411B (zh) | 针对射程不确定性的质子重离子调强放疗鲁棒优化方法 | |
CN115829972B (zh) | 一种放疗计划三维剂量分布与通量的同步预测方法及装置 | |
CN112086172A (zh) | 一种三维剂量计算方法、计算机设备以及可读介质 | |
CN105561485A (zh) | 放射治疗计划的优化方法和装置 | |
CN106204749A (zh) | 稀疏低秩特征表达的三维模型共分割方法 | |
Yousefi et al. | Novel Method for More Efficient Optimizing the Knowledge-Based Planning: Specific Voxels of each Structure Influenced by Dominant Beamlets (SVSIDB) | |
CN111494815B (zh) | 基于混合变尺度模型的三维剂量计算方法、装置及介质 | |
CN112086173B (zh) | 三维剂量计算方法、装置、计算机设备及可读介质 | |
CN105477789A (zh) | 基于二次规划模型抑制总出束时间的动态调强放疗方法 | |
CN102136041A (zh) | 治疗计划逆向规划方法和治疗计划*** | |
CN116341182A (zh) | 基于多目标进化和传统优化方法协同的质子调强多目标优化方法及其*** | |
CN116870377A (zh) | 一种基于深度学习的肝脏肿瘤放疗剂量预测方法 | |
CN113178242B (zh) | 一种基于耦合式生成对抗网络的自动计划优化*** | |
CN113516233B (zh) | 一种面向vmat放疗计划的神经网络预测装置 | |
CN115310373A (zh) | 一种制氢电解槽剩余寿命预测方法 | |
US20230264048A1 (en) | Method and Apparatus for Controlling Multileaf Collimator in Radiotherapy | |
CN112419439B (zh) | 一种治疗计划中质子重离子能量快速反演的方法 | |
Liu et al. | NeuralDAO: Incorporating neural network generated dose into direct aperture optimization for end‐to‐end IMRT planning | |
CN110289075B (zh) | 一种基于模糊熵的直接子野优化方法及*** | |
CN110876839B (zh) | 不均匀网格分布模拟直线加速器治疗计划的剂量计算方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C41 | Transfer of patent application or patent right or utility model | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20150924 Address after: Hefei City, Anhui province Jinzhai road 230061 No. 310 tiger building 6 floor Applicant after: Hefei Chaoan Medical Technology Co., Ltd. Address before: Daoxiang road in Shushan District of Hefei City, Anhui province 230031 No. 9 Building Room 618 Applicant before: Hefei Chaoan Medical Technology Co., Ltd. |
|
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
AV01 | Patent right actively abandoned | ||
AV01 | Patent right actively abandoned |
Granted publication date: 20151223 Effective date of abandoning: 20190130 |