CN115715852A - 一种基于蒙特卡洛模拟的囊结构优化方法 - Google Patents

一种基于蒙特卡洛模拟的囊结构优化方法 Download PDF

Info

Publication number
CN115715852A
CN115715852A CN202211444518.5A CN202211444518A CN115715852A CN 115715852 A CN115715852 A CN 115715852A CN 202211444518 A CN202211444518 A CN 202211444518A CN 115715852 A CN115715852 A CN 115715852A
Authority
CN
China
Prior art keywords
capsule
capsule structure
optimization
target parameters
diameter
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
Application number
CN202211444518.5A
Other languages
English (en)
Inventor
李东洁
梁雨
王雪莹
姚钢
徐东昊
杨柳
张宇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin University of Science and Technology
Original Assignee
Harbin University of Science and Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Harbin University of Science and Technology filed Critical Harbin University of Science and Technology
Priority to CN202211444518.5A priority Critical patent/CN115715852A/zh
Publication of CN115715852A publication Critical patent/CN115715852A/zh
Pending legal-status Critical Current

Links

Images

Landscapes

  • Radiation-Therapy Devices (AREA)

Abstract

一种基于蒙特卡洛模拟的囊结构优化方法,涉及肿瘤放化疗植入囊结构优化领域。目的是解决传统放化疗计划中囊结构优化需要依靠大量迭代试错的问题。具体包括以下步骤:首先,选择一组目标参数作为囊结构初始值,并针对囊结构辐照特点确定一个治疗评判标准;然后,结合评判标准,不断改变目标参数,循环将囊结构进行蒙特卡洛模拟,直至得到一组优化的目标参数;最后,对优化目标参数的囊结构执行剂量场计算,并检查计算结果是否符合优化要求,若符合,则结束流程,若不符合,对该组优化目标参数进行微调,返回蒙特卡洛模拟。本发明实现了囊结构最优,可以最大程度放化疗肿瘤,为放射治疗前手术方案的制定和调整提供依据。

Description

一种基于蒙特卡洛模拟的囊结构优化方法
技术领域
本发明涉及放化疗植入囊结构优化领域。具体涉及一种基于蒙特卡洛模拟的囊结构优化设计方法。
背景技术
随着放射性治疗在医学领域的进一步发展,大约有70%的癌症病人在治疗过程中会使用到放疗。在放射治疗前,医生往往会根据CT诊断的结果来判定肿瘤组织的大小进而选择放射性治疗中使用的植入囊的大小,以此达到最佳的放射治疗效果。当根据肿瘤组织大小确定了植入囊大小时,植入囊结构设计决定了辐照剂量场范围大小和肿瘤治疗的效果好坏。如何优化囊结构,是放化疗囊植入的重中之重。在放化疗计划中使用传统技术对植入囊结构优化不仅需要依靠人工大量迭代试错,还需要诊断过程主治医师有较高的治疗水平,同时也增加了诊断时间和工作量。为了解决这一问题,本发明提出了一种基于蒙特卡洛模拟的囊结构优化方法,实现上述过程自动优化,摆脱主治医师人工深度参与,并可得到最优化的植入囊结构,可以最大程度放化疗肿瘤,促进治疗向个体化、精准化治疗发展。本发明适用于放射性治疗术前为放射治疗前手术方案的制定和调整提供依据。
发明内容
1.本发明为解决植入囊结构最优化的问题,提出一种基于蒙特卡洛模拟的囊结构优化方法,其特征在于,该基于蒙特卡洛模拟的囊结构优化方法包括如下步骤:
步骤一:确定囊整体尺寸结构,根据优化要求,选择一组目标参数作为囊结构初始值,并针对囊结构辐照特点确定一个治疗评判标准;
步骤二:结合评判标准,按照预先设定的步长,不断改变已选择的目标参数的取值范围,循环将囊结构进行蒙特卡洛模拟,直至得到一组最优化的目标参数;
步骤三:对最优化目标参数的囊结构执行剂量场计算,并检查计算结果是否符合优化要求,若符合,则结束流程,得到最优化囊结构,若不符合,则进入步骤四;
步骤四:对该组优化目标参数进行微调,返回步骤二。
2.作为进一步优化说明,步骤一中,常根据CT诊断的结果来判定肿瘤组织的大小进而选择放射性治疗中使用的植入囊的大小,选择现有尺寸的囊结构或者适合对应肿瘤尺寸的囊结构作为初始值。
3.作为进一步优化说明,步骤一中,目标参数包括螺距、囊直径、微腔直径、薄膜宽度、囊壁最小厚度、囊材料;
囊螺距:囊螺距与囊螺旋圈数相关,在囊整体尺寸确定的情况下,囊螺距越大,囊螺旋圈数越少;囊螺距影响囊直径,当螺距较小时,囊直径不能过大;通常,囊螺距取值范围在0-5mm之间;
囊直径:囊直径为微腔所在圆心面的直径,囊直径决定微腔直径,囊直径要始终大于微腔直径;通常,囊直径取值范围在0.1-1mm之间;
微腔直径:微腔分为化疗微腔与放疗微腔,微腔直径要始终小于囊直径;通常,微腔直径取值范围为囊直径尺寸的30%-100%之间;
薄膜宽度:两个微腔之间由薄膜相隔,薄膜宽度影响微腔之间的距离;通常,薄膜宽度取值范围在0-0.1mm之间;
囊壁最小厚度:腔壁与囊壁之间的最小距离为囊壁最小厚度;通常,囊壁最小厚度取值范围在0-0.7mm之间;
囊材料:囊材料选用可以植入人体的常见医用材料,常用于作为植入囊的材料有聚四氟乙烯、聚醚醚酮、氮化硅、硅橡胶和钛合金等。
4.作为进一步优化说明,步骤一中治疗评判标准由主治医师根据病种特点设定,由辐照病灶处肿瘤和组织器官的辐射敏感性来定义约束条件;最简单评判标准为肿瘤区域有效剂量和非肿瘤区域有效剂量比值,值越大表示囊结构参数越合理。
5.作为进一步优化说明,步骤二中按照预先设定的步长,改变目标参数为螺距、囊直径、微腔直径、薄膜宽度、囊壁最小厚度、囊材料中各参数预先设定的取值范围。
6.作为进一步优化说明,步骤二中预先设定的步长、取值范围考虑制造工艺、病人对不同材料过敏程度等因素影响。
7.作为进一步优化说明,步骤二中将所有组进行蒙特卡洛模拟得到的结果结合评价标准得到对应评价标准值,对所有组的评价标准值排序,得到一组最优化的目标参数。
8.作为进一步优化说明,步骤三中检查计算结果为各目标参数之间是否存在冲突,若发生参数冲突,进入步骤四。
9.作为进一步优化说明,步骤三中计算结果不存在冲突,对最优化目标参数的囊结构执行剂量场计算,与初始囊结构剂量场对比,是否符合优化要求,若符合,则结束流程,得到最优化囊结构,若不符合,则进入步骤四。
10.作为进一步优化说明,步骤四中对冲突的优化目标参数进行微调,直至各目标参数取值合理。
11.作为进一步优化说明,步骤四中微调方法为针对冲突的优化目标参数,取值范围和步长进行小范围的调整,返回步骤二再对冲突的优化目标参数进行少量几次优化,直至得到优化目标参数合理的囊结构。
本发明与现有技术相比,具有如下的优点和有益效果:
本文提出来一种基于蒙特卡洛模拟的囊结构优化方法。目的是解决传统技术中放化疗计划植入囊结构优化需要依靠人工大量迭代试错的问题。本发明能够获得更加客观准确的不同囊结构的放射治疗计划优化结果,解决了人工优化过程中选择和修改目标参数非常依赖主治医师的经验,实现了囊结构的自动优化,找到植入囊最优结构,可以最大程度放化疗肿瘤,促进治疗向个体化、精准化治疗发展。
附图说明
图1为本发明中基于蒙特卡洛模拟的囊结构优化方法的流程图;
图2为本发明中肿瘤放化治疗植入囊结构图。
具体实施方式:
为了使本发明的目的、技术方案及优点更加的清楚明白,以下结合附图及实施例,对本发明进一步详细说明。
囊结构参数影响空间剂量场和肿瘤辐照效果。一般情况下,一组优化目标参数有很多,每个优化目标参数取值范围也较大,使得使用传统技术对囊结构优化需要依靠人工大量迭代试错,依赖主治医师的经验。本发明使用蒙特卡洛模拟方法对上述过程自动优化,摆脱主治医师人工深度参与,并可得到最优化的植入囊结构,可以最大程度放化疗肿瘤。
本实例中,结合图2说明,双腔囊整体直径为20mm,在保持双腔囊整体直径不变的情况下,对双腔囊结构优化。初始双腔囊结构参数为:螺距为2mm,囊直径为1mm,微腔直径为0.4mm,薄膜宽度为0.1mm,囊壁最小厚度为0.05mm,囊材料为聚四氟乙烯。
结合图1说明,本实例基于蒙特卡洛模拟的囊结构优化方法按照以下步骤进行:
步骤一:确定囊整体尺寸结构,根据优化要求,选择一组目标参数作为囊结构初始值,并针对囊结构辐照特点确定一个治疗评判标准;
本步骤中,选择现有或者医学上针对肿瘤大小常用的囊结构作为初始值,一般脑胶质瘤大小最大至5cm左右,选取整体直径为20mm作为双腔囊结构;所选择的目标参数有螺距、囊直径、微腔直径、薄膜宽度、囊壁最小厚度、囊材料;选取螺距区间为1-3mm之间;选取囊直径区间为0.8-1mm之间;选取微腔直径为囊直径的30%-40%之间;选取囊壁最小厚度为0-0.2mm;选取薄膜宽度为0-0.4mm之间;选取囊材料为聚四氟乙烯或聚醚醚酮;治疗评判标准由主治医师根据病种特点设定,由辐照病灶处肿瘤和组织器官的辐射敏感性来定义约束条件,以最简单的肿瘤区域有效剂量和非肿瘤区域有效剂量比值作为评判标准,值越大表示囊结构参数越合理。
步骤二:结合评判标准,按照预先设定的步长,不断改变已选择的目标参数的取值范围,循环将囊结构进行蒙特卡洛模拟,直至得到一组最优化的目标参数;
本步骤中,选取螺距预先设定的步长为1mm,其余预先设定的步长为0.1mm;根据对应目标参数预先设定的取值范围,改变目标参数,取值范围考虑制造工艺、病人对材料过敏程度等因素影响;将改变参数后的囊结构进行蒙特卡洛模拟,将结果结合评判标准,得出具体评价标准值;不断改变不同参数范围,循环将不同囊结构参数进行蒙特卡洛模拟,得到评价标准值;根据评价标准值大小,得到一组最优的目标参数。
步骤三:对最优化目标参数的囊结构执行剂量场计算,并检查计算结果是否符合优化要求,若符合,则结束流程,得到最优化囊结构,若不符合,则进入步骤四;
本步骤中,检查最优化囊结构各项参数之间是否存在冲突,并执行剂量场计算,与初始囊结构对比优化效果是否符合优化要求,若符合要求,则结束流程,若不符合,则进入步骤四。
步骤四:对该组优化目标参数进行微调,返回步骤二。
本步骤中,对冲突的优化目标参数进行微调,针对冲突的优化目标参数,取值范围和步长进行小范围的调整,直至各目标参数取值合理;返回步骤二再对冲突的优化目标参数进行少量几次优化,直至得到优化目标参数合理的囊结构。最终得到的最优化囊结构参数为螺距为2mm,囊直径为0.8mm,囊微腔直径为0.32mm,囊壁最小厚度为0.04mm,薄膜宽度为0.08mm,囊材料为聚醚醚酮。

Claims (1)

1.一种基于蒙特卡洛模拟的囊结构优化方法,其特征在于,该方法包括如下步骤:
步骤一:确定囊整体尺寸结构,根据优化要求,选择一组目标参数作为囊结构初始值,并针对囊结构辐照特点确定一个治疗评判标准;
所述步骤一中:根据CT诊断的结果来判定肿瘤组织的大小进而选择放射性治疗中使用的植入囊的大小,选择现有尺寸的囊结构或者适合对应肿瘤尺寸的囊结构作为初始值;目标参数包括:螺距、囊直径、微腔直径、薄膜宽度、囊壁最小厚度、囊材料;囊螺距取值范围在0-5mm之间;囊直径取值范围在0.1-1mm之间;微腔直径取值范围为囊直径尺寸的30%-100%之间;薄膜宽度取值范围在0-0.1mm之间;囊壁最小厚度取值范围在0-0.7mm之间;常用于作为植入囊的材料有聚四氟乙烯、聚醚醚酮、氮化硅、硅橡胶和钛合金等;治疗评判标准由主治医师根据病种特点设定,由辐照病灶处肿瘤和组织器官的辐射敏感性来定义约束条件;
步骤二:结合评判标准,按照预先设定的步长,不断改变已选择的目标参数的取值范围,循环将囊结构进行蒙特卡洛模拟,直至得到一组最优化的目标参数;
所述步骤二中:按照预先设定的步长,改变目标参数为螺距、囊直径、微腔直径、薄膜宽度、囊壁最小厚度、囊材料中各参数预先设定的取值范围;预先设定的步长、取值范围需要考虑制造工艺、病人对不同材料过敏程度等因素影响;
步骤三:对最优化目标参数的囊结构执行剂量场计算,并检查计算结果是否符合优化要求,若符合,则结束流程,得到最优化囊结构,若不符合,则进入步骤四;
所述步骤三中:检查计算结果为各目标参数之间是否存在冲突,若发生参数冲突,进入步骤四;若计算结果不存在冲突,对最优化目标参数的囊结构执行剂量场计算,与初始囊结构剂量场对比,是否符合优化要求,若符合,则结束流程,得到最优化囊结构,若不符合,则进入步骤四;
步骤四:对该组优化目标参数进行微调,返回步骤二;
所述步骤四中:微调方法为针对冲突的优化目标参数,取值范围和步长进行小范围的调整,返回步骤二再对冲突的优化目标参数进行少量几次优化,直至得到优化目标参数合理的囊结构。
CN202211444518.5A 2022-11-18 2022-11-18 一种基于蒙特卡洛模拟的囊结构优化方法 Pending CN115715852A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211444518.5A CN115715852A (zh) 2022-11-18 2022-11-18 一种基于蒙特卡洛模拟的囊结构优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211444518.5A CN115715852A (zh) 2022-11-18 2022-11-18 一种基于蒙特卡洛模拟的囊结构优化方法

Publications (1)

Publication Number Publication Date
CN115715852A true CN115715852A (zh) 2023-02-28

Family

ID=85255466

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211444518.5A Pending CN115715852A (zh) 2022-11-18 2022-11-18 一种基于蒙特卡洛模拟的囊结构优化方法

Country Status (1)

Country Link
CN (1) CN115715852A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116312955A (zh) * 2023-03-29 2023-06-23 中国科学院苏州生物医学工程技术研究所 用于控制和调节膝骨关节炎光治疗剂量的系统和方法

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1537648A (zh) * 2003-10-22 2004-10-20 高春平 球囊导管型肿瘤近距离内放射治疗装置
CN109248385A (zh) * 2018-08-20 2019-01-22 北京东方瑞云科技有限公司 基于蒙特卡洛树搜索的放射治疗计划优化系统及方法
CN208823772U (zh) * 2017-09-18 2019-05-07 王娟 一种可短期插值治疗肿瘤的球囊粒子链
CN110141770A (zh) * 2019-04-02 2019-08-20 成都真实维度科技有限公司 一种瘤体内放射性粒子源植入针道路径规划及针道生成方法
US20200126683A1 (en) * 2017-04-24 2020-04-23 Advanced Accelerator Applications Accelerator-driven neutron activator for brachytherapy
CN211536247U (zh) * 2019-12-12 2020-09-22 戴建荣 用于放疗的囊状施照器
CN112587808A (zh) * 2020-12-08 2021-04-02 北京航空航天大学 粒子支架剂量分布计算方法
CN113181563A (zh) * 2021-04-30 2021-07-30 珠海横乐医学科技有限公司 粒子植入肿瘤内放疗剂量规划方法、系统及介质
CN115317809A (zh) * 2022-08-12 2022-11-11 华中科技大学同济医学院附属协和医院 一种采用结果引导实现放射性粒子植入的方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1537648A (zh) * 2003-10-22 2004-10-20 高春平 球囊导管型肿瘤近距离内放射治疗装置
US20200126683A1 (en) * 2017-04-24 2020-04-23 Advanced Accelerator Applications Accelerator-driven neutron activator for brachytherapy
CN208823772U (zh) * 2017-09-18 2019-05-07 王娟 一种可短期插值治疗肿瘤的球囊粒子链
CN109248385A (zh) * 2018-08-20 2019-01-22 北京东方瑞云科技有限公司 基于蒙特卡洛树搜索的放射治疗计划优化系统及方法
CN110141770A (zh) * 2019-04-02 2019-08-20 成都真实维度科技有限公司 一种瘤体内放射性粒子源植入针道路径规划及针道生成方法
CN211536247U (zh) * 2019-12-12 2020-09-22 戴建荣 用于放疗的囊状施照器
CN112587808A (zh) * 2020-12-08 2021-04-02 北京航空航天大学 粒子支架剂量分布计算方法
CN113181563A (zh) * 2021-04-30 2021-07-30 珠海横乐医学科技有限公司 粒子植入肿瘤内放疗剂量规划方法、系统及介质
CN115317809A (zh) * 2022-08-12 2022-11-11 华中科技大学同济医学院附属协和医院 一种采用结果引导实现放射性粒子植入的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
万斌、姚进: "近距离放射治疗系统的研究进展", 《中国医疗设备》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116312955A (zh) * 2023-03-29 2023-06-23 中国科学院苏州生物医学工程技术研究所 用于控制和调节膝骨关节炎光治疗剂量的系统和方法
CN116312955B (zh) * 2023-03-29 2023-10-20 中国科学院苏州生物医学工程技术研究所 用于控制和调节膝骨关节炎光治疗剂量的系统和方法

Similar Documents

Publication Publication Date Title
CN101184527B (zh) 并入功能性成像信息的放射治疗方案设备
Kalkman et al. Current focussing in cochlear implants: an analysis of neural recruitment in a computational model
Solin et al. Three-dimensional photon treatment planning of the intact breast
CN115715852A (zh) 一种基于蒙特卡洛模拟的囊结构优化方法
Guidolin et al. Stereotactic image-guided neoadjuvant ablative single-dose radiation, then lumpectomy, for early breast cancer: the SIGNAL prospective single-arm trial of single-dose radiation therapy
CA2515244A1 (en) Arrangement of seeds and spacers for brachytherapy
US9486646B2 (en) System and method for control of external beam radiation
CN106169028A (zh) 调强放疗中危及器官吸收剂量预测数学模型的构建方法
Morrison et al. Advanced collapsed cone engine dose calculations in tissue media for COMS eye plaques loaded with I‐125 seeds
US20030095923A1 (en) System for imaging malignant tumors using carbon 13 with MRI
CN117545529A (zh) 促进向异质身体施用治疗辐射的方法和装置
Morrison et al. Delivered dose uncertainty analysis at the tumor apex for ocular brachytherapy
CN113262390B (zh) 一种基于放射性粒子植入的肿瘤放疗效果评价方法
Paschal et al. Monte Carlo modeling of the Elekta Versa HD and patient dose calculation with EGSnrc/BEAMnrc
Chand et al. Effect of CT number to relative electron density curves acquired at different tube voltage and current on radiotherapy dose calculation
Ur Rab et al. Giant phyllodes tumor of breast post fibroadenoma excision: a rare case report and review of literature
Emberru Assessment of a Treatment Planning Protocol for the Reduction of Dosimetry Calculation Errors in Radiotherapy for Head and Neck Patients With Dental Implants
CN217612537U (zh) 一种用于宫颈癌的后装精准近距离治疗装置
Fuss et al. Intensity-modulated radiosurgery: improving dose gradients and maximum dose using post inverse-optimization interactive dose shaping
Tien et al. Dosimetric study of Varian Universal Multi-Channel Cylinder System for High-Dose-Rate 192Ir Brachytherapy
Craig et al. Comparison of correction protocols for image-guided radiation therapy
Barten et al. Artificial intelligence based planning of HDR prostate brachytherapy: first clinical experience
Rak et al. Einfluss cochleärer Parameter auf die aktuelle Cochleaimplantatversorgung: Entwicklung eines Konzepts für die personalisierte Medizin
Fijałkowski et al. Three-dimensional (3D) real-time conformal brachytherapy-a novel solution for prostate cancer treatment. Part I. Rationale and method
Yadav et al. Development of Anthropomorphic Heterogeneous Pelvic Phantom and Its Comparison with Homogeneous Phantom in Advance Radiation Therapy: Dosimetry Analysis

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination