CN103977499A - 粒子植入治疗计划设计方法 - Google Patents
粒子植入治疗计划设计方法 Download PDFInfo
- Publication number
- CN103977499A CN103977499A CN201410216372.8A CN201410216372A CN103977499A CN 103977499 A CN103977499 A CN 103977499A CN 201410216372 A CN201410216372 A CN 201410216372A CN 103977499 A CN103977499 A CN 103977499A
- Authority
- CN
- China
- Prior art keywords
- treatment plan
- particles
- particle
- msub
- target area
- 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.)
- Pending
Links
- 239000002245 particle Substances 0.000 title claims abstract description 109
- 238000011282 treatment Methods 0.000 title claims abstract description 70
- 238000000034 method Methods 0.000 title claims abstract description 45
- 238000002513 implantation Methods 0.000 title claims abstract description 23
- 238000011156 evaluation Methods 0.000 claims abstract description 17
- 210000000920 organ at risk Anatomy 0.000 claims description 14
- 239000000203 mixture Substances 0.000 claims description 2
- 238000005457 optimization Methods 0.000 abstract description 27
- 238000011221 initial treatment Methods 0.000 abstract description 2
- 238000013461 design Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 206010028980 Neoplasm Diseases 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 230000005855 radiation Effects 0.000 description 3
- 230000002285 radioactive effect Effects 0.000 description 3
- 230000002776 aggregation Effects 0.000 description 2
- 238000004220 aggregation Methods 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 210000002307 prostate Anatomy 0.000 description 2
- 238000002922 simulated annealing Methods 0.000 description 2
- 238000002560 therapeutic procedure Methods 0.000 description 2
- 210000001519 tissue Anatomy 0.000 description 2
- 230000004931 aggregating effect Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000002725 brachytherapy Methods 0.000 description 1
- 230000007012 clinical effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 210000003128 head Anatomy 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 239000007943 implant Substances 0.000 description 1
- 210000004185 liver Anatomy 0.000 description 1
- 210000004072 lung Anatomy 0.000 description 1
- 238000002595 magnetic resonance imaging Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 210000003739 neck Anatomy 0.000 description 1
- 210000000496 pancreas Anatomy 0.000 description 1
- 230000008092 positive effect Effects 0.000 description 1
- 230000000693 radiobiological effect Effects 0.000 description 1
- 238000001959 radiotherapy Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 208000024891 symptom Diseases 0.000 description 1
- 230000001225 therapeutic effect Effects 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
Landscapes
- Radiation-Therapy Devices (AREA)
Abstract
本发明是一种粒子植入治疗计划设计方法,针对治疗计划对剂量分布的要求,提出了一种有效的粒子判优准则,合理地利用了粒子和靶区以及OARs的相对位置关系的先验知识,能够快速地得出治疗计划。本方法首先获得所有可能位置的粒子,根据判优准则确定各粒子的评价值,评价值最小的粒子为当前最优粒子;不断将当前最优粒子植入靶区,直到靶区处方剂量覆盖度达到预定值,得到初步的治疗计划;当能减少穿刺针数目时,移除对应粒子,重新选取当前最优粒子植入,得到新的治疗计划;直到穿刺针数目不能减少时得到最终的治疗计划。本发明结合启发性优化方法的特点,能够在保证剂量分布的前提下,有效地减少穿刺针数,提高了粒子植入治疗计划的可行性。
Description
技术领域
本发明涉及计算机技术,具体涉及一种粒子植入治疗计划设计方法。
背景技术
放射性粒子植入治疗也被称为近距离放射治疗,是一种用于治疗肿瘤的放射治疗手段。粒子植入通过模板引导,经皮穿刺等介入手段将放射性粒子(125I,103Pd)植入癌变肿瘤的区域,在对肿瘤进行持续性低剂量照射的同时,能有效地减小危及器官(OARs)辐射剂量,具有较好的放射生物效应,取得了较好的临床效果。此外,粒子植入治疗还具有副作用小,使用安全、费用低、微创等特点,被广泛地应用于头颈部、肺部、肝部、胰腺、前列腺等部位肿瘤的治疗。
放射性粒子植入治疗成功的关键是合理规划粒子在靶区的排列分布,确保靶区的辐射剂量达到处方剂量且尽可能的均匀分布,同时减小周围危及器官承受的剂量。为了达到治疗效果,临床治疗中需要通过超声、计算机断层扫描(CT)或者核磁共振成像(MRI)等成像设备确定靶区的形状和位置,以及靶区周围的临近组织以及相对位置关系,然后制定治疗计划,确定粒子的数目以及分布位置。
传统的治疗计划设计是一个不断试错修改(trial and error)的过程:物理师根据临床经验设计出初步的治疗计划,然后计算剂量,根据靶区和危及器官内的剂量分布调整粒子数目和位置,重复这个步骤直到得出一个可行的治疗计划。然而,治疗计划设计通常需要综合权衡考虑多个目标,设计的过程繁重且耗时较长,而且在很大程度上依赖于物理师的临床经验。
现有针对治疗计划设计的优化方法针对前列腺粒子植入提出,根据优化方法求解的手段可以分为:随机性(stochastic)优化方法,确定性(deterministic)优化方法和启发性(heuristic)优化方法。随机性优化方法包括:模拟退火方法(simulated annealing,SA)和遗传方法(geneticalgorithm,GA)。随机性方法的优化模型中将靶区和OAR等效为椭圆或者圆柱,通过随机的确定植入粒子的位置得到治疗计划,然后计算剂量并评估模型中的剂量分布,得出当前计划的评价度,并应用评价度指导优化过程。基于简化模型的剂量分布评价规则在一定程度上并不能体现真实的剂量分布情况,实际应用中很可能得出不可行的治疗计划。同时随机性优化方法针对穿刺针道已经确定的情况提出,优化问题相对简单。对于针道尚未确定、靶区体积较大的病例,优化问题较为复杂,应用随机性优化方法有可能陷入局部最优解的死循环。确定性优化方法将所有可能位置是否植入粒子编码为1和0。根据对靶区最大最小剂量值、剂量分布均整度,OAR最大剂量以及穿刺针数量的要求确定多目标函数和约束条件,通过分支定界(branch and bound)方法求解混合整数规划(mixed-integer programming)。由于优化模型比较复杂,计算量较大,计算时间仍然长达数小时。此外,目标函数中各项的权重也直接影响着优化结果,在实际应用中需要根据不同的病症进行调整。现有的启发性优化方法以粒子对靶区和OAR的剂量贡献度为评价准则,选择当前最优粒子直到靶区处方剂量覆盖度达到预定值。由于评价准则没有考虑靶区剂量的均整度,启发性优化方法采用一种基于等剂量线的约束规则防止粒子聚集。当针数达到预定值之后,将可选粒子的位置限制在已有穿刺针上,如果当前没有可用粒子,同时靶区剂量未达到预定要求,增加针数阈值,移除所有粒子重新优化。相对于随机性优化方法和确定性优化方法,启发性优化方法能够较快地得出治疗计划,但是,该方法针对避免粒子聚集和减少穿刺针提出的机制还存在一定的问题,实际应用中很可能得到穿刺针较少但剂量分布较差的治疗计划,同时迭代次数较多,计算时间较长,不能满足术中计划设计的临床需要。
发明内容
本发明针对治疗计划对剂量分布的要求,提供了一种粒子植入治疗计划设计方法,本发明方法提出了一种有效的粒子判优准则,合理地利用了粒子和靶区以及OARs的相对位置关系的先验知识,能够快速地得出治疗计划。本发明提供的粒子植入治疗计划设计方法结合启发性优化方法的特点,提出用于减少穿刺针数目的双重迭代策略,能够在保证剂量分布的前提下,有效地减少了穿刺针数,提高了粒子植入治疗计划计划的可行性。
本发明提供的粒子植入治疗计划设计方法,包括如下步骤:
步骤1:设所有可能位置的粒子组成集合S,计算各粒子对靶区和OARs的平均剂量;
步骤2:按照判优准则确定集合S中各粒子的评价值,根据评价值对粒子进行排序,评价值最小的粒子为当前最优粒子;
判优准则Cs定义为:
其中,为粒子对靶区的平均剂量;为粒子对不同危及器官的平均剂量,w1,w2为对应的权重;D为粒子植入之后靶区内剂量分布的标准差;
步骤3:设治疗计划中包含的粒子组成集合Z,集合Z初始为空;处方剂量覆盖靶区的体积比定义为靶区处方剂量覆盖度。如果当前靶区处方剂量覆盖度小于预定值Ct,将当前最优粒子从集合S中删除,加入治疗计划集合Z中,重新计算靶区内剂量分布,然后执行步骤2;否则,当前治疗计划集合Z就是一个初步的治疗计划,转步骤4执行。
步骤4:若第一次执行到本步骤,直接执行步骤5;否则,判断当前的治疗计划与上一次治疗计划相比,穿刺针的针数是否减少,若是,执行步骤5,否则执行步骤6;
步骤5:找到当前的治疗计划中只包含一个粒子的穿刺针,从集合Z中移除对应粒子,并将移除的粒子加入粒子集合S中,然后重新计算靶区内剂量分布,继续转步骤2执行;
步骤6:将当前的治疗计划作为最终治疗计划。
本发明的优点与积极效果在于:
1、本发明根据粒子植入治疗的临床目标提出了一种有效的判优准则,同时合理地利用了粒子和靶区以及OARs的相对位置关系的先验知识;能够在保证靶区处方剂量覆盖度的前提下,使靶区内的剂量尽可能均匀地分布,同时减小对OARs的辐射剂量;
2、利用启发性优化方法的特点,本发明方法是一种双重迭代策略,利用启发性优化方法的特点,在保证剂量分布满足要求的前提下,减少了穿刺针的数目,提高了粒子植入治疗计划的可行性。
附图说明
图1是本发明提供的粒子植入治疗计划设计方法的流程示意图。
具体实施方式
下面将结合附图和实施例对本发明作进一步的详细说明。
如图1所示,本发明提供的粒子植入治疗计划设计方法包括如下步骤:
步骤1:确定所有可能位置的粒子,并计算各粒子对靶区和OARs的平均剂量,记为集合S。
步骤2:按照判优准则对集合S中的粒子进行排序。判优准则Cs定义为:
其中,为粒子对靶区的平均剂量,为粒子对不同OAR组织的平均剂量,w1,w2为对应的权重,实际应用中可以根据OAR的敏感程度进行调整。D为将粒子植入靶区之后靶区内剂量分布的标准差。当前靶区内剂量分布根据不断植入最优粒子或者移除粒子后有更新,所以对应D值也有更新,需要重新计算。
根据式(1)获得各粒子的评价值,然后再根据评价值从小达到的顺序将粒子进行排序,评价值最小的粒子为当前最优粒子。
步骤3:治疗计划中包含的粒子组成集合Z。处方剂量覆盖靶区的体积比定义为靶区处方剂量覆盖度。如果当前靶区处方剂量覆盖度小于预定值Ct,将当前最优粒子从集合S中删除,加入治疗计划集合Z中,重新计算靶区内剂量分布,然后执行步骤2;否则,当前治疗计划集合Z就是一个初步的治疗计划,转步骤4执行。
预定值Ct设定为大于等于95%的值。治疗计划集合Z初始为空。
重复步骤2到步骤3,直到靶区处方覆盖度达到预定值,从而得出一个初步的治疗计划。
步骤4:若第一次执行到此步骤,直接执行步骤5;否则,判断当前的治疗计划与上一次治疗计划相比,穿刺针的针数是否减少,若是,执行步骤5,否则执行步骤6。
步骤5:找到当前的治疗计划中只包含一个粒子的穿刺针,从集合Z中移除对应粒子,并将移除的粒子加入粒子集合S中,然后重新计算靶区内剂量分布,转步骤2继续执行。
步骤6:得到最终的治疗计划,当前治疗计划就是最终治疗计划。
在选取粒子的过程中没有考虑穿刺针数的限制,因此初步的治疗计划很可能使用了较多的穿刺针,降低了治疗计划的可行性。本发明方法采用双重迭代策略,利用启发性优化方法的特点,在保证剂量分布满足要求的前提下,减少穿刺针的数目。
首先移除初步治疗计划中只包含一个粒子的穿刺针(单个粒子),然后更新剂量场,重新计算靶区内剂量分布,根据粒子判优准则选取当前最优粒子,重复该过程直到靶区覆盖度达到预定的要求,得到一个更新的治疗计划。继续移除更新之后的治疗计划中的单个粒子,重复选取粒子的过程,直到治疗计划使用的穿刺针数不再减少。
按照评价准则搜索得到的是当前粒子分布情况下的最优粒子,并不一定是全局最优粒子。被移除的单个粒子很可能不会被再次选中。对于其他粒子,一旦被选中就一直被保留下来。双重迭代过程利用启发性搜索的特点,不断的移除治疗计划中的单个粒子,没有对穿刺针数做强制的限制,能够在保证剂量分布不退化的前提下减少穿刺针的数目。
Cs中的平均剂量比值应用了粒子与靶区和OARs的相对位置关系的先验知识,D考虑了粒子植入之后对靶区剂量均整度的影响,在避免粒子聚集的前提下,优先选择靠近靶区中心且远离OARs的粒子。
Claims (2)
1.一种粒子植入治疗计划设计方法,其特征在于,包括如下步骤:
步骤1:设所有可能位置的粒子组成集合S,确定各粒子对靶区和危及器官OARs的平均剂量;
步骤2:按照判优准则确定集合S中各粒子的评价值Cs,根据评价值对粒子进行排序,评价值最小的粒子为当前最优粒子;
定义为:
其中,为粒子对靶区的平均剂量;为粒子对不同危及器官的平均剂量,w1,w2为对应的权重;D为粒子植入当前靶区之后靶区内剂量分布的标准差;
步骤3:设治疗计划中包含的粒子组成集合Z,集合Z初始为空;处方剂量覆盖靶区的体积比定义为靶区处方剂量覆盖度;如果当前靶区处方剂量覆盖度小于预定值Ct,将当前最优粒子从集合S中删除,加入治疗计划集合Z中,重新计算靶区内剂量分布,然后执行步骤2;否则,当前治疗计划集合Z就是一个初步的治疗计划,转步骤4执行;
步骤4:若第一次执行到本步骤,直接执行步骤5;否则,判断当前的治疗计划与上一次治疗计划相比,穿刺针的针数是否减少,若是,执行步骤5,否则执行步骤6;
步骤5:找到当前的治疗计划中只包含一个粒子的穿刺针,从集合Z中移除对应粒子,并将移除的粒子加入粒子集合S中,然后重新计算靶区内剂量分布,继续转步骤2执行;
步骤6:将当前的治疗计划作为最终治疗计划。
2.根据权利要求1所述的一种粒子植入治疗计划设计方法,其特征在于,步骤3中所述的预定值Ct设定为大于等于95%的值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410216372.8A CN103977499A (zh) | 2014-05-21 | 2014-05-21 | 粒子植入治疗计划设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410216372.8A CN103977499A (zh) | 2014-05-21 | 2014-05-21 | 粒子植入治疗计划设计方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN103977499A true CN103977499A (zh) | 2014-08-13 |
Family
ID=51269815
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410216372.8A Pending CN103977499A (zh) | 2014-05-21 | 2014-05-21 | 粒子植入治疗计划设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103977499A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105797284A (zh) * | 2016-04-25 | 2016-07-27 | 张宏涛 | 一种放射性粒子放射剂量验证方法及系统 |
CN111862022A (zh) * | 2020-07-13 | 2020-10-30 | 中山大学 | 全身多部位放射治疗危及器官自动勾画方法 |
CN112933428A (zh) * | 2021-01-28 | 2021-06-11 | 苏州普能医疗科技有限公司 | 插植导板设计方法和插植导板 |
CN112957620A (zh) * | 2021-01-28 | 2021-06-15 | 苏州普能医疗科技有限公司 | 近距离放射治疗计划设计系统和方法 |
CN113181563A (zh) * | 2021-04-30 | 2021-07-30 | 珠海横乐医学科技有限公司 | 粒子植入肿瘤内放疗剂量规划方法、系统及介质 |
CN113877072A (zh) * | 2020-07-02 | 2022-01-04 | 南京大学 | 一种近距离放射治疗中对粒子源分布进行优化的动态规划算法 |
-
2014
- 2014-05-21 CN CN201410216372.8A patent/CN103977499A/zh active Pending
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105797284A (zh) * | 2016-04-25 | 2016-07-27 | 张宏涛 | 一种放射性粒子放射剂量验证方法及系统 |
CN113877072A (zh) * | 2020-07-02 | 2022-01-04 | 南京大学 | 一种近距离放射治疗中对粒子源分布进行优化的动态规划算法 |
CN111862022A (zh) * | 2020-07-13 | 2020-10-30 | 中山大学 | 全身多部位放射治疗危及器官自动勾画方法 |
CN111862022B (zh) * | 2020-07-13 | 2022-08-16 | 广州柏视医疗科技有限公司 | 全身多部位放射治疗危及器官自动勾画方法 |
CN112933428A (zh) * | 2021-01-28 | 2021-06-11 | 苏州普能医疗科技有限公司 | 插植导板设计方法和插植导板 |
CN112957620A (zh) * | 2021-01-28 | 2021-06-15 | 苏州普能医疗科技有限公司 | 近距离放射治疗计划设计系统和方法 |
CN112933428B (zh) * | 2021-01-28 | 2022-09-27 | 苏州普能医疗科技有限公司 | 插植导板设计方法和插植导板 |
CN113181563A (zh) * | 2021-04-30 | 2021-07-30 | 珠海横乐医学科技有限公司 | 粒子植入肿瘤内放疗剂量规划方法、系统及介质 |
CN113181563B (zh) * | 2021-04-30 | 2023-01-03 | 珠海横乐医学科技有限公司 | 粒子植入肿瘤内放疗剂量规划方法、系统及介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Byrne et al. | Varian ethos online adaptive radiotherapy for prostate cancer: Early results of contouring accuracy, treatment plan quality, and treatment time | |
US10328282B2 (en) | System and method for novel chance-constrained optimization in intensity-modulated proton therapy planning to account for range and patient setup uncertainties | |
US10850122B2 (en) | Optimization methods for radiation therapy planning | |
CN103977499A (zh) | 粒子植入治疗计划设计方法 | |
CN103282967B (zh) | 用于放射疗法的自动化治疗计划 | |
CN110327554B (zh) | 基于预测剂量分布引导的调强放疗计划优化方法及应用 | |
US11738211B2 (en) | Dose-distribution estimation in proton therapy | |
JP2013529481A (ja) | 放射線量を推定し、推定された放射線量を操作するためのシステムおよび方法 | |
JP2016511117A (ja) | 同位体源外部ビーム放射線治療のためのシステム及び方法 | |
Pokharel et al. | Evaluation of hybrid inverse planning and optimization (HIPO) algorithm for optimization in real‐time, high‐dose‐rate (HDR) brachytherapy for prostate | |
US20180111005A1 (en) | Method of selecting beam geometries | |
Anbumani et al. | Treatment planning methods in high dose rate interstitial brachytherapy of carcinoma cervix: a dosimetric and radiobiological analysis | |
Betts et al. | Optimised robust treatment plans for prostate cancer focal brachytherapy | |
Yang et al. | Stereotactic magnetic resonance-guided adaptive radiation therapy (SMART) for abdominopelvic oligometastases | |
Luong et al. | Efficient, effective, and insightful tackling of the high-dose-rate brachytherapy treatment planning problem for prostate cancer using evolutionary multi-objective optimization algorithms | |
CN109248385B (zh) | 基于蒙特卡洛树搜索的放射治疗计划优化系统 | |
Niedzielski et al. | Analysis of esophageal‐sparing treatment plans for patients with high‐grade esophagitis | |
Sachpazidis et al. | Investigating the role of constrained CVT and CVT in HIPO inverse planning for HDR brachytherapy of prostate cancer | |
Shaw et al. | Equivalence of Gyn GEC-ESTRO guidelines for image guided cervical brachytherapy with EUD-based dose prescription | |
CN114588561A (zh) | 射野强度分布优化方法、装置及放射治疗设备 | |
Sarkar et al. | Prostate implant dosimetric outcomes and migration patterns between bio-absorbable coated and uncoated brachytherapy seeds | |
CN103977498B (zh) | 一种放射性粒子锥形布源的三维构建方法 | |
Shen et al. | Automatic dose prediction using deep learning and plan optimization with finite‐element control for intensity modulated radiation therapy | |
Galuppi et al. | EP-1254: Proctoscopic evaluation of rectal toxicity after radiohormone therapy in prostate cancer | |
Lausch et al. | iCONE‐SRS: Development of inverse treatment planning for cone‐based stereotactic radiosurgery |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20140813 |