CN105930636B - 一种自动确定目标函数权重的放射治疗方案优化系统 - Google Patents
一种自动确定目标函数权重的放射治疗方案优化系统 Download PDFInfo
- Publication number
- CN105930636B CN105930636B CN201610217049.1A CN201610217049A CN105930636B CN 105930636 B CN105930636 B CN 105930636B CN 201610217049 A CN201610217049 A CN 201610217049A CN 105930636 B CN105930636 B CN 105930636B
- Authority
- CN
- China
- Prior art keywords
- information
- optimization
- weight
- module
- scheme optimization
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/20—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Medical Informatics (AREA)
- Public Health (AREA)
- Data Mining & Analysis (AREA)
- Databases & Information Systems (AREA)
- Pathology (AREA)
- Epidemiology (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Radiation-Therapy Devices (AREA)
Abstract
本发明公开一种自动确定目标函数权重的放射治疗方案优化系统,包括信息输入模块、预处理模块、方案优化模块、权重调整模块、方案优化结果输出模块;信息输入模块获取病人三维信息、器官勾画信息、治疗头信息,目标函数信息、以及目标函数中所使用器官的DV限制参数;预处理模块通过剂量计算引擎,计算不同照射方向的剂量沉积矩阵;方案优化模块通过方案优化引擎,进行方案优化;权重调整模块根据方案优化模块的输入信息,计算各器官的剂量分布和各子目标函数值,并判断各子目标函数是否都满足优化的终止条件,方案优化结果输出模块将方案优化信息输出。本发明自动方案优化是一个迭代的过程,不需要输入理想DVH曲线。
Description
技术领域
本发明属于放射治疗方案优化技术领域,涉及自动确定目标函数权重的放射治疗方案优化系统。
背景技术
在放射治疗的方案优化过程中,方案优化的总目标函数一般表示为各子目标函数与对应权重乘积的总和。子目标函数的权重表示各子目标函数在优化过程中的重要性。在方案优化之前,各子目标函数所使用的权重值是预先由物理师根据经验预先设定的,且该权重值在优化过程中是固定不变的。在方案优化结束后,物理师对优化结果进行方案优化质量评价。如果优化结果不满足临床放射治疗要求,物理师根据优化结果,适当地调整各子目标函数对应的权重值,并重新进行方案优化。该操作循环进行,直到方案优化结果满足临床放射治疗的要求为止。在这种手动调整权重的方法中,优化过程需要循环执行多次,增加了制作放射治疗方案所需的时间;同时还需要物理师根据经验对各子目标函数的权重值进行修改,增加了物理师的工作量。针对各子目标函数的权重取值问题,已有学者展开了多目标优化问题的研究。对于传统多目标优化,首先求出目标函数的帕累托最优边界,然后在帕累托最优边界上由物理师或医生选取最优的优化结果。当总目标函数中包含较多的子目标函数,求解帕累托最优边界是非常耗时的。而且,需要物理师或医生在帕累托边界上选取最优的优化结果,同样增加了物理师或医生的工作量。因此,需要提出一种快速的自动放射治疗方案优化技术,在自动优化的过程中,根据优化结果自动调整子目标函数的权重值。
在现有自动放射治疗方案优化技术中,兰(Nan)提出了一中基于已有剂量-体积直方图(dose-volume histogram,DVH)的自适应放射治疗重新优化方法(具体参见Nan L,Masoud Z,Andres U-S,et al.Automatic treatment plan re-optimization foradaptive radiotherapy guided with the initial plan DVHs.Phys.Med.Biol.,2013,58:8725-8738)。该方法的优点是可以在优化过程中自动的调整子目标函数的权重值;缺点是需要预先知道理想DVH分布。在自适应放射治疗过程中,可以用同一个病人前一次分次照射方案优化的DVH分布近似地代替该病人当前分次照射的理想DVH分布。然而在一般的方案优化中,由于病人解剖结构之间的差异非常大,用这种近似代替的方法是不可行的。因此,需要发明一种不需要理想DVH的自动方案优化方法和系统。
一般从三维剂量分布、DVH曲线以及生物准则等方面对放射治疗方案优化结果的质量评价。查看三维剂量分布,主要为了确认靶区高剂量分布的覆盖情况,以及周围正常组织接受到的高剂量照射。传统的优化质量评价主要是根据美国肿瘤放射治疗协作组织(theRadiation Therapy Oncology Group,RTOG)或临床定量分析正常组织效应准则(theQuantitative Analysis of Normal Tissue Effects in the Clinic,QUANTEC)给出的评价标准对DVH进行定量的评价。通常方案优化中物理子目标函数的设置就是参考这些评价标准进行设定的,并且满足这些评价标准的优化结果被认为是满足临床要求的。因此,在优化过程中,可以根据实际剂量分布与这些评价标准自动的差值自动地修改对应的子目标函数权重,使得优化结果满足所有的评价标准。在临床放射治疗计划系统中,生物准则已经被作为子目标函数参与放射治疗方案优化,以及被用来定量地评价方案优化结果。使用生物准则对方案优化结果进行方案优化和质量评价是放射治疗未来发展的一个趋势。同样,相对于物理子目标函数,生物子目标函数也可以使用类似的方法进行自动权重调整。
发明内容
本发明主要是针对传统放射治疗方案制作过程中物理师循环调整子目标函数并进行重新优化的技术问题,提出了一种自动确定目标函数权重的放射治疗方案优化系统。
本发明通过以下技术方案实现:
一种自动确定目标函数权重的放射治疗方案优化系统,包括信息输入模块、预处理模块、方案优化模块、权重调整模块、方案优化结果输出模块;
所述的信息输入模块包括3个功能单元:病人三维信息获取功能单元、器官勾画功能单元和其它信息输入功能单元,在病人三维信息获取功能单元中,病人通过硬件设备进行成像,得到病人的三维密度信息,在器官勾画功能单元中进行自动勾画而获取器官勾画信息,从而获得病人的各器官勾画信息,在其它信息输入功能单元中输入治疗头信息,目标函数信息,以及目标函数中所使用器官的DV限制参数;
所述的预处理模块包含一个剂量计算功能单元,在使用当前治疗头的情况下,通过剂量计算引擎,计算不同照射方向的剂量沉积矩阵,首先,根据获取的病人的各器官勾画信息、等中心点和照射方向确定每个方向需要计算的射束;然后,根据病人的三维密度信息和治疗头信息,通过剂量计算引擎计算每个射束的剂量分布,从而得到剂量沉积矩阵;
所述的方案优化模块包括利用已知信息输入模块的其它信息输入功能单元的目标函数和预处理模块的剂量沉积矩阵,通过方案优化引擎,进行方案优化,首先,在所有的照射方向中选择一个最佳的子野形状;其次,通过子野权重优化算法重新计算已有子野的权重值;最后,删除权重值为零的子野,并对剩余子野形状进行微调,以达到子野形状最优,如果总的子野个数超过预设的上限值,则在迭代过程中不再增加新的子野,只做子野权重优化,删除子野和子野形状微调,方案优化进行一段时间或迭代一定次数之后停止优化,并输出优化结果到权重调整模块;
所述的权重调整模块根据方案优化模块的输入信息,计算各器官的剂量分布和各子目标函数值,并判断各子目标函数是否都满足优化的终止条件,如果方案优化结果满足临床要求,方案优化结束并将方案优化结果输出到方案优化结果输出模块;否则修正各子目标函数的权重,继续方案优化,权重修改后,将按照新的子目标函数权重继续之前停止的方案优化,方案优化和权重调整循环进行,直到方案优化的终止条件都满足而停止。
所述的方案优化结果输出模块包括两个功能单元:放射治疗功能单元和优化结果显示功能单元,放射治疗功能单元将方案优化所得子野形状、子野权重、射线源等信息输入到放射治疗设备,以控制放射治疗设备对病人进行放射治疗,其中子野形状和子野权重分别用来控制多叶准直器的开口形状和在对应开口形状下的射线照射时间长短,优化结果显示功能单元将优化结果中的三维剂量分布和DVH曲线分别输入到计算机和打印机,以显示方案优化的结果和生成报表。
进一步,治疗头信息包括完整的治疗头形状结构和参数信息,以及照射方向,等中心的位置等。目标函数信息为针对当前具体的病例设置的目标函数,包括各子目标函数的类型、所使用的参数值。
进一步,在器官勾画功能单元中进行自动勾画而获取器官勾画信息,可以通过物理师在所获取的三维密度信息上进行勾画。
进一步,在器官勾画功能单元中进行自动勾画而获取器官勾画信息,可以通过自动勾画软件进行自动勾画而获取病人的各器官勾画信息。
进一步,剂量计算引擎所采用的算法可以是笔形束算法、点核算法、蒙特卡罗算法和确定性方法等。
进一步,所述的子目标函数可以分为7类,分别为最大剂量子目标函数、最小剂量子目标函数、均匀剂量子目标函数、DVH子目标函数、等效均一剂量子目标函数、肿瘤控制率子目标函数和正常组织并发症概率子目标函数。
本发明与现有技术相比,优点如下:
(1)本发明可以进行自动方案优化,自动方案优化是一个迭代的过程,在迭代过程中,根据当前方案优化的结果与给定的处方剂量的差值自动修正各子目标函数的权重,对不满足条件的目标函数的增加惩罚,从而保证迭代结束后所有目标函数都满足条件,实现自动方案优化,相比Nan提出的方法,本发明不需要输入理想DVH曲线。
(2)本发明可以应用于传统的方案优化,也可以应用于自适应放射治疗的方案优化中;
(3)本发明可以利用生物准则进行自动方案优化。
附图说明:
下面结合附图对本发明的具体实施方式作进一步详细说明。
图1是本发明的自动确定目标函数权重的放射治疗方案优化系统的具体实施流程图;
图2是本发明的最大DVH限制的权重修正示意图;
图3是本发明的最小DVH限制的权重修正示意图;
图4是本发明的信息输入模块的硬件构成示意图;
图5是本发明的方案优化结果输出模块构成示意图;
图6是本发明的自动确定目标函数权重的放射治疗方案优化系统组成框图。
具体实施方式
为使本发明实施例的目的、特点和优点能够更为明显易懂,下面结合附图对本发明的具体实施作详细的说明。
如图6所示,本发明提供了一种自动确定目标函数权重的放射治疗方案优化系统1,包括信息输入模块11、预处理模块12、方案优化模块13、权重调整模块14、方案优化结果输出模块15;所述的信息输入模块11的硬件构成如图4所示,该模块包括3个功能单元:病人三维信息获取功能单元111、器官勾画功能单元112和其它信息输入功能单元113。在病人三维信息获取功能单元111中,病人通过硬件设备进行成像,得到病人的三维密度信息,硬件常用的设备是螺旋CT。为了增加螺旋CT成像的质量,采用功能成像设备(如核磁共振成像设备)对病人进行成像,以增强螺旋CT成像的质量或直接用来表示病人的三维密度信息。在器官勾画功能单元112中进行自动勾画而获取器官勾画信息,可以通过物理师在所获取的三维密度信息上进行勾画,从而获得病人的各器官勾画信息,也可以通过自动勾画软件进行自动勾画而获取病人的各器官勾画信息。在其它信息输入功能单元113中输入治疗头信息,目标函数信息,以及目标函数中所使用器官的DV限制参数等。治疗头信息包括完整的治疗头形状结构和参数信息,以及照射方向,等中心的位置等。目标函数信息为针对当前具体的病例设置的目标函数,包括各子目标函数的类型、所使用的参数值等。
所述的器官勾画功能单元112和其它信息输入功能单元113都在计算机上实现。以上所述硬件设备以及器官勾画的软件可以在大型医疗设备公司获取得到。
所述的预处理模块12包含一个剂量计算功能单元(图未视出),在使用当前治疗头的情况下,通过剂量计算引擎,计算不同照射方向的剂量沉积矩阵。实施的技术方案是,首先,根据获取的病人的各器官勾画信息、等中心点和照射方向确定每个方向需要计算的射束;然后,根据病人的三维密度信息和治疗头信息,通过剂量计算引擎计算每个射束的剂量分布,从而得到剂量沉积矩阵。剂量计算引擎所采用的算法可以是笔形束算法、点核算法、蒙特卡罗算法和确定性方法等。在这个模块中,我们使用经典笔形束算法计算剂量分布。为了提高笔形束算法的计算速度,该算法被嵌入到硬件设备中,形成单独的剂量计算硬件设备。
通常,该剂量计算功能单元在计算机上实现,同样其软件可以在大型医疗设备公司都获取得到。由于该功能单元的执行无需用户的参与,因此该功能单元可以由专门的计算剂量分布的硬件设备实现。
所述的方案优化模块13包括利用已知信息输入模块11的其它信息输入功能单元113的目标函数和预处理模块12的剂量沉积矩阵,通过方案优化引擎,进行方案优化。实施的技术方案是,在已知目标函数和剂量沉积矩阵的情况下,就可以进行方案优化来确定最终各子野的形状和权重。方案优化算法可以采用两步法(注量图优化和射野分割)或直接子野优化方法。在本发明具体实施方式中,以直接子野优化方法为例进行说明。子野优化是一个迭代的过程。在每次迭代中,首先,在所有的照射方向中选择一个最佳的子野形状;其次,通过子野权重优化算法重新计算已有子野的权重值;最后,删除权重值为零的子野,并对剩余子野形状进行微调,以达到子野形状最优。如果总的子野个数超过预设的上限值,则在迭代过程中不再增加新的子野,只做子野权重优化,删除子野和子野形状微调。方案优化进行一段时间或迭代一定次数之后停止优化,并输出优化结果到权重调整模块14,如每迭代十次,停止优化并将结果输出到权重调整模块14。
所述的权重调整模块14根据方案优化模块13的输入信息,计算各器官的剂量分布和各子目标函数值(限制参数),并判断各子目标函数是否都满足优化的终止条件,如果方案优化结果满足临床要求,方案优化结束并将方案优化结果输出到方案优化结果输出模块15;否则修正各子目标函数的权重,继续方案优化,权重修改后,将按照新的子目标函数权重继续之前停止的方案优化,方案优化和权重调整循环进行,直到方案优化的终止条件都满足而停止。
所述的子目标函数可以分为7类,分别为最大剂量子目标函数、最小剂量子目标函数、均匀剂量子目标函数、DVH子目标函数、等效均一剂量(equivalent uniform dose,EUD)子目标函数、肿瘤控制率(tumor control probability,TCP)子目标函数和正常组织并发症概率(normal tissue complication probability,NTCP)子目标函数。所有子目标函数权重的采用相似的修正方法,当前子目标函数的权重值为wold,修改后的子目标函数权重值为wnew,二者的关系为
wnew=wold×factor,
或者
wnew=wold+factor,
其中factor为修正因子。在优化过程中,需要根据方案优化模块13中的优化时间或迭代次数来决定选择使用乘法或加法对权重进行修正。如果系统设置的优化时间较短或迭代次数比较少,这时需要使用加法进行校正,以避免因权重变化太大而引起的优化结果反复振荡现象。如果系统设置的优化时间较长或迭代次数比较多,这时需要使用乘法进行校正,以缩短优化所需的时间。
以下是7种不同的子目标函数权重修正因子的具体计算方法:
1、最大剂量子目标函数:
factor=D2/D1,
其中D1为处方剂量,D2为该子目标函数对应器官的最大剂量值;
2、最小剂量子目标函数:
factor=D1/D2,
其中D1为处方剂量,D2为该子目标函数对应器官的最小剂量值;
3、均匀剂量子目标函数:
factor=max(D1,D2)/min(D1,D2),
其中D1为处方剂量,D2为该子目标函数对应器官的均匀剂量值;
4、DVH子目标函数:对于最大DVH限制,如果对应DVH曲线从该DVH限制点的上方通过,则认为该子目标函数所对应器官的剂量分布不满足临床需求,如图2所示。在图2中,DVH限制为点S201,其对应的剂量-体积值分别为D1和V1。而在当前的DVH曲线上,体积为V1的点在S202处,对应剂量值为D2;剂量为D1的点在S203处,对应体积为V2。在点S201、S202和S203所围成的多边形范围内的剂量分布是不满足当前DVH限制的。因此,需要增加子目标函数的权重值,以期望在下次方案优化后该子目标函数满足要求。此时,修正因子可以根据剂量条件获得:
factor=D2/D1;
也可以根据体积条件计算得到:
factor=V2/V1。
对于最小DVH限制,如果DVH曲线从该限制点的下方通过,则认为该DVH所对应器官的剂量分布不满足临床需求,如图3所示,DVH限制为点S301,其对应的剂量-体积值分别为D1和V1。而在当前的DVH曲线上,体积为V1的点在S302处,对应剂量值为D2;剂量为D1的点在S303处,对应体积为V2。在点S301、S302和S303所围成的多边形范围内的剂量分布是不满足当前DVH限制的,需要修正对应子目标函数的权重。修正方法与上述最大DVH限制的修正方法相似,其中根据剂量条件和体积条件计算所得修正因子分别为
factor=D1/D2
和
factor=V1/V2;
5、EUD子目标函数:当EUD的参数a值大于等于1时,该子目标函数被用来控制正常组织或靶区上的高剂量分布,此时修正因子的计算方法与最大子目标函数修正因子的计算方法相似。当参数a值小于1时,该子目标函数被用来控制靶区的低剂量分布,此时修正因子的计算方法与最小子目标函数修正因子的计算方法相似;
6、TCP子目标函数:计算方法与最小剂量子目标函数修正因子的计算方法相似;
7、NTCP子目标函数:计算方法与最大剂量子目标函数修正因子的计算方法相似。
所述的方案优化模块13和权重调整模块14是本发明的主要部分。这两个模块分别包含了方案优化功能单元(图未视出)和权重调整功能单元(图未视出),由于这两个功能单的运行无需用户的参与,这两个功能单元可以在计算机上执行。在实现过程中,由于方案优化模块13和权重调整模块14被循环调用,因此,为了提高运算速度,这两个模块的实现也上可以被嵌入到一个硬件设备上,以快速地实现方案优化功能和权重调整功能,从而形成独立的方案优化和权重调整硬件设备。
所述的方案优化结果输出模块15包括两个功能单元:放射治疗功能单元151和优化结果显示功能单元152(如图5所示)。放射治疗功能单元151将方案优化所得子野形状、子野权重、射线源等信息输入到放射治疗设备,以控制放射治疗设备对病人进行放射治疗,其中子野形状和子野权重分别用来控制多叶准直器的开口形状和在对应开口形状下的射线照射时间长短。常用的放射治疗设备可以在大型医疗设备公司获取得到。优化结果显示功能单元152将优化结果中的三维剂量分布和DVH曲线分别输入到计算机和打印机,以显示方案优化的结果和生成报表。
如图1、图6所示,本发明的自动确定目标函数权重的放射治疗方案优化方法:
步骤S101为信息输入,在信息输入模块11中,输入病人的三维密度信息、器官勾画信息、治疗头信息,目标函数信息,以及目标函数中所使用器官的DV限制参数等。其中病人的三维密度信息可以是CT图像、MR图像或其它方式获取的病人密度信息。器官勾画信息可以通过物理师在三维密度信息上进行勾画而获得的,也可以通过自动勾画软件进行自动勾画而获取。治疗头信息包括完整的治疗头形状结构和参数信息,以及照射方向,等中心的位置等。目标函数信息为针对当前具体的病例设置的目标函数,包括各子目标函数的类型、所使用的参数等。
步骤S102为预处理,在预处理模块12中,首先,根据获取的病人的各器官勾画信息、等中心点和照射方向确定每个方向需要计算的射束;然后,根据病人的三维密度信息和治疗头信息,通过剂量计算引擎计算每个射束的剂量分布,从而得到剂量沉积矩阵。剂量计算引擎所采用的算法可以是笔形束算法、点核算法、蒙特卡罗算法和确定性方法等。在这个模块中,我们使用经典笔形束算法计算剂量分布。为了提高笔形束算法的计算速度,该算法被嵌入到硬件设备中,形成单独的剂量计算硬件设备。
步骤S103为方案优化,在方案优化模块13中,在已知信息输入模块11的其它信息输入功能单元113的目标函数和预处理模块12的剂量沉积矩阵的情况下,就可以进行方案优化来确定最终各子野的形状和权重。方案优化算法可以采用两步法(注量图优化和射野分割)或直接子野优化方法。在本发明具体实施方式中,以直接子野优化方法为例进行说明。子野优化是一个迭代的过程。在每次迭代中,首先,在所有的照射方向中选择一个最佳的子野形状;其次,通过子野权重优化算法重新计算已有子野的权重值;最后,删除权重值为零的子野,并对剩余子野形状进行微调,以达到子野形状最优。每迭代十次,停止优化并将结果输出到权重调整模块。如果总的子野个数超过预设的上限值,则在迭代过程中不再增加新的子野,只做子野权重优化,删除子野和子野形状微调。
步骤S104为权重调整,在权重调整模块14中,根据方案优化模块13的输入信息,计算各器官的剂量分布和各子目标函数的值,并判断各子目标函数是否都满足优化的终止条件。如果所有子目标函数的条件都满足,方案优化结束,输出优化结果,否则,按照发明内容中的权重修正方式对各子目标函数的权重就行修正。权重修改后,将按照新的子目标函数权重继续之前停止的方案优化。方案优化和权重调整循环进行,直到方案优化的终止条件都满足而停止。
在实现过程中,方案优化模块13和权重调整模块14被循环调用。因此,为了提高运算速度,这两个模块的实现被嵌入到一个硬件设备上,以快速地实现方案优化功能和权重调整功能。
步骤S105为输出方案优化结果,在方案优化结果输出模块15中,方案优化结果输出模块输出三维剂量分布、DVH曲线、子野形状和子野权重等。在这个模块中子野形状、子野权重以及其它关于射线源的信息(参数)被输入到放射治疗设备,由这些输入参数来控制放射治疗设备,对病人进行放射治疗,其中子野形状和子野权重分别用来控制多叶准直器的开口形状和在对应开口形状下的射线照射时间长。三维剂量分布和DVH曲线被输入到计算机和打印机,分别用来显示方案优化结果和打印报表。
本发明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干的改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (6)
1.一种自动确定目标函数权重的放射治疗方案优化系统,其特征在于:包括信息输入模块、预处理模块、方案优化模块、权重调整模块、方案优化结果输出模块;
所述的信息输入模块包括3个功能单元:病人三维信息获取功能单元、器官勾画功能单元和其它信息输入功能单元,在病人三维信息获取功能单元中,病人通过硬件设备进行成像,得到病人的三维密度信息,在器官勾画功能单元中进行自动勾画而获取器官勾画信息,从而获得病人的各器官勾画信息,在其它信息输入功能单元中输入治疗头信息,目标函数信息,以及目标函数中所使用器官的DV限制参数;
所述的预处理模块包含一个剂量计算功能单元,在使用当前治疗头的情况下,通过剂量计算引擎,计算不同照射方向的剂量沉积矩阵,首先,根据获取的病人的各器官勾画信息、等中心点和照射方向确定每个方向需要计算的射束;然后,根据病人的三维密度信息和治疗头信息,通过剂量计算引擎计算每个射束的剂量分布,从而得到剂量沉积矩阵;
所述的方案优化模块包括利用已知信息输入模块的其它信息输入功能单元的目标函数和预处理模块的剂量沉积矩阵,通过方案优化引擎,进行方案优化,首先,在所有的照射方向中选择一个最佳的子野形状;其次,通过子野权重优化算法重新计算已有子野的权重值;最后,删除权重值为零的子野,并对剩余子野形状进行微调,以达到子野形状最优,如果总的子野个数超过预设的上限值,则在迭代过程中不再增加新的子野,只做子野权重优化,删除子野和子野形状微调,方案优化进行一段时间或迭代一定次数之后停止优化,并输出优化结果到权重调整模块;
所述的权重调整模块根据方案优化模块的输入信息,计算各器官的剂量分布和各子目标函数值,并判断各子目标函数是否都满足优化的终止条件,如果方案优化结果满足临床要求,方案优化结束并将方案优化结果输出到方案优化结果输出模块;否则修正各子目标函数的权重,继续方案优化,权重修改后,将按照新的子目标函数权重继续之前停止的方案优化,方案优化和权重调整循环进行,直到方案优化的终止条件都满足而停止;
所述的方案优化结果输出模块包括两个功能单元:放射治疗功能单元和优化结果显示功能单元,放射治疗功能单元将方案优化所得子野形状、子野权重、射线源信息输入到放射治疗设备,以控制放射治疗设备对病人进行放射治疗,其中子野形状和子野权重分别用来控制多叶准直器的开口形状和在对应开口形状下的射线照射时间长短,优化结果显示功能单元将优化结果中的三维剂量分布和DVH曲线分别输入到计算机和打印机,以显示方案优化的结果和生成报表。
2.如权利要求1所述的一种自动确定目标函数权重的放射治疗方案优化系统,其特征在于:治疗头信息包括完整的治疗头形状结构和参数信息,以及照射方向,等中心的位置;目标函数信息为针对当前具体的病例设置的目标函数,包括各子目标函数的类型、所使用的参数值。
3.如权利要求1所述的一种自动确定目标函数权重的放射治疗方案优化系统,其特征在于:在器官勾画功能单元中进行自动勾画而获取器官勾画信息,通过物理师在所获取的三维密度信息上进行勾画。
4.如权利要求1所述的一种自动确定目标函数权重的放射治疗方案优化系统,其特征在于:在器官勾画功能单元中进行自动勾画而获取器官勾画信息,通过自动勾画软件进行自动勾画而获取病人的各器官勾画信息。
5.如权利要求1所述的一种自动确定目标函数权重的放射治疗方案优化系统,其特征在于:剂量计算引擎所采用的算法可以是笔形束算法、点核算法、蒙特卡罗算法和确定性方法。
6.如权利要求1所述的一种自动确定目标函数权重的放射治疗方案优化系统,其特征在于:所述的子目标函数分为7类,分别为最大剂量子目标函数、最小剂量子目标函数、均匀剂量子目标函数、DVH子目标函数、等效均一剂量子目标函数、肿瘤控制率子目标函数和正常组织并发症概率子目标函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610217049.1A CN105930636B (zh) | 2016-03-29 | 2016-03-29 | 一种自动确定目标函数权重的放射治疗方案优化系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610217049.1A CN105930636B (zh) | 2016-03-29 | 2016-03-29 | 一种自动确定目标函数权重的放射治疗方案优化系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105930636A CN105930636A (zh) | 2016-09-07 |
CN105930636B true CN105930636B (zh) | 2019-07-09 |
Family
ID=56840553
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610217049.1A Active CN105930636B (zh) | 2016-03-29 | 2016-03-29 | 一种自动确定目标函数权重的放射治疗方案优化系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105930636B (zh) |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
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 |
CN108057177A (zh) * | 2017-12-11 | 2018-05-22 | 太原工业学院 | 一种基于模糊推理的混合准则放疗参数自动优化方法 |
CN108364678B (zh) * | 2018-02-10 | 2022-02-01 | 福建省肿瘤医院(福建省肿瘤研究所、福建省癌症防治中心) | 一种基于机器学习的鼻咽癌放疗计划自动化设计方法 |
CN109125952B (zh) * | 2018-07-18 | 2020-07-31 | 中北大学 | 基于核模型的卷积叠加能量沉积计算方法 |
CN109248385B (zh) * | 2018-08-20 | 2022-05-13 | 北京东方瑞云科技有限公司 | 基于蒙特卡洛树搜索的放射治疗计划优化系统 |
CN109192273B (zh) * | 2018-09-10 | 2021-09-28 | 东莞东阳光高能医疗设备有限公司 | 硼中子俘获治疗的束流评估方法、装置、设备及存储介质 |
CN110404184A (zh) * | 2019-06-13 | 2019-11-05 | 苏州同调医学科技有限公司 | 一种测算放疗射线剂量分布和剂量目标函数的方法和系统 |
WO2021036366A1 (zh) * | 2019-08-29 | 2021-03-04 | 北京连心医疗科技有限公司 | 标准化的人工智能自动放疗计划方法和系统 |
CN110415785A (zh) * | 2019-08-29 | 2019-11-05 | 北京连心医疗科技有限公司 | 人工智能引导放疗计划的方法及系统 |
CN113744385A (zh) * | 2020-05-29 | 2021-12-03 | 美敦力(上海)管理有限公司 | 器官三维数字模型的显示方法和装置 |
WO2022104497A1 (zh) * | 2020-11-17 | 2022-05-27 | 西安大医集团股份有限公司 | 靶点的权重确定方法、装置及放射治疗系统 |
CN113117253B (zh) * | 2021-04-20 | 2022-09-16 | 中北大学 | 一种基于各向异性核的剂量计算系统 |
CN112972912B (zh) * | 2021-04-22 | 2021-08-17 | 北京大学第三医院(北京大学第三临床医学院) | 非共面放疗射束入射路径优化方法 |
CN114712733A (zh) * | 2022-04-29 | 2022-07-08 | 上海联影医疗科技股份有限公司 | 一种自动确定射野信息的方法、系统和装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104866928A (zh) * | 2015-06-04 | 2015-08-26 | 上海联影医疗科技有限公司 | 用于自适应放射治疗的优化方法及系统 |
CN105031819A (zh) * | 2015-08-25 | 2015-11-11 | 上海联影医疗科技有限公司 | 一种剂量优化方法与系统 |
-
2016
- 2016-03-29 CN CN201610217049.1A patent/CN105930636B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104866928A (zh) * | 2015-06-04 | 2015-08-26 | 上海联影医疗科技有限公司 | 用于自适应放射治疗的优化方法及系统 |
CN105031819A (zh) * | 2015-08-25 | 2015-11-11 | 上海联影医疗科技有限公司 | 一种剂量优化方法与系统 |
Non-Patent Citations (2)
Title |
---|
放射治疗计划的自动优化及再优化关键技术研究;李楠;《万方数据》;20141209;第1-104页 |
精确放射治疗剂量计算及方案优化方法研究;张鹏程;《万方数据》;20150730;第1-120页 |
Also Published As
Publication number | Publication date |
---|---|
CN105930636A (zh) | 2016-09-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105930636B (zh) | 一种自动确定目标函数权重的放射治疗方案优化系统 | |
Sibolt et al. | Clinical implementation of artificial intelligence-driven cone-beam computed tomography-guided online adaptive radiotherapy in the pelvic region | |
US20210085997A1 (en) | Three-dimensional radiotherapy dose distribution prediction | |
Hellebust et al. | Dosimetric impact of interobserver variability in MRI-based delineation for cervical cancer brachytherapy | |
CN105930637A (zh) | 一种自动确定目标函数权重的放射治疗方案优化方法 | |
Voet et al. | Fully automated volumetric modulated arc therapy plan generation for prostate cancer patients | |
US7817778B2 (en) | Interactive treatment plan optimization for radiation therapy | |
EP2686068B1 (en) | Studying dosimetric impact of motion to generate adaptive patient-specific margins in ebrt planning | |
US7809107B2 (en) | Method for controlling modulation strength in radiation therapy | |
US8952346B2 (en) | Systems and methods for isotopic source external beam radiotherapy | |
JP2016525427A (ja) | 治療計画の自動化 | |
CN108711447A (zh) | 基于体素权重因子的自动调强多目标优化方法及其应用 | |
EP3573715B1 (en) | Estimation of the achievability of treatment goals in radiation therapy | |
Van Haveren et al. | Lexicographic extension of the reference point method applied in radiation therapy treatment planning | |
EP3634578B1 (en) | Treatment plan generation for radiation therapy treatment | |
Nguyen et al. | Computerized triplet beam orientation optimization for MRI‐guided Co‐60 radiotherapy | |
Scaggion et al. | Updating a clinical Knowledge-Based Planning prediction model for prostate radiotherapy | |
US7206377B2 (en) | Predictive organ dynamics database and code | |
CN107998519A (zh) | 一种用于imrt的分子动力学再优化算法 | |
Hong et al. | Clinical experience of automated SBRT paraspinal and other metastatic tumor planning with constrained hierarchical optimization | |
CN113870976A (zh) | 一种基于剂量预评估的自适应放疗剂量调强优化计算方法 | |
Baker et al. | Real-time interactive planning for radiotherapy of head and neck cancer with volumetric modulated arc therapy | |
Torelli | Spatiotemporal Fractionation in Radiotherapy | |
Cilla et al. | Automated VMAT treatment planning for complex cancer cases: A Feasibility Study | |
Cavus et al. | Optimizing volumetric modulated arc therapy prostate planning using an automated Fine-Tuning process through dynamic adjustment of optimization parameters |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |