CN104361253A - 一种利用规则化方法确定后装源驻留时间的方法 - Google Patents
一种利用规则化方法确定后装源驻留时间的方法 Download PDFInfo
- Publication number
- CN104361253A CN104361253A CN201410713531.5A CN201410713531A CN104361253A CN 104361253 A CN104361253 A CN 104361253A CN 201410713531 A CN201410713531 A CN 201410713531A CN 104361253 A CN104361253 A CN 104361253A
- Authority
- CN
- China
- Prior art keywords
- residence time
- formula
- reference point
- dosage
- vector
- 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.)
- Granted
Links
Landscapes
- Radiation-Therapy Devices (AREA)
- Nuclear Medicine (AREA)
- Measuring And Recording Apparatus For Diagnosis (AREA)
Abstract
本发明公开了一种利用规则化方法确定后装源驻留时间的方法,包括以下步骤;(I)根据病人的CT图像,确定靶区及紧要器官的轮廓,设置驻留点位置和参考点位置及处方剂量;(II)根据式(1)计算参考点的计量剂量:(III)根据计算剂量与处方剂量的偏差、相邻驻留位的时间差构建评价函数式(2):O(T)=(D0-AT)2+αTT·T;(IV)评价函数对驻留时间向量T进行微分,得到式(3);(V)简化式(3)得到一个正规方程组,并用向量表示如下:(B+αI)T=F;(VI)通过改变不同的α值,得到不同的驻留时间解集。本发明利用规则化方法获得后装驻留时间,同时,引入时间平滑因子解决了相邻驻留位中驻留时间相差大的问题。
Description
技术领域
本发明涉及一种确定后装驻留时间的方法,具体的说,是涉及一种利用规则化方法确定后装源驻留时间的方法。
背景技术
后装放射治疗是指把不带放射源的施源器置于治疗部位,由电脑遥控步进电机将放射源送入容器进行放射治疗,如此可避免放置治疗容器过程中医务人员因放射受伤。可以将放射源准确安全地输送到患者需要治疗的部位进行放射治疗。由于放置位置准确、距病体组织近等优点,在治疗妇科、鼻咽、食道、支气管、直肠、膀胱、乳腺及胰腺等肿瘤中取得了明显的临床治疗效果。
三维近距放射治疗计划系统的应用也越来越广泛,基本思路是根据临床对放射源周围剂量分布的特殊要求,控制步进源在不同驻留位置停留不同的时间,提高治疗质量。根据临床所规定的靶区处方剂量(即参考剂量)值,对步进源的驻留时间进行优化处理。然而,在优化求解的过程中由于临床情况的复杂性,如驻留位较多、剂量参考点不能与驻留位一一对应等,尤其是相邻驻留位中驻留时间相差大,使得计算结果中出现负值驻留时间,这在临床和物理学上都是没有意义的;同时,解的结果中,相邻驻留位的驻留时间可能相差很大,临床上也很难接受,因为它会造成某些临近驻留位的剂量过低或者过高。
发明内容
本发明的目的在于克服上述缺陷,提供一种利用规则化方法确定后装源驻留时间的方法。
为了实现上述目的,本发明采用技术方案如下:
一种利用规则化方法确定后装驻留时间的方法,包括以下步骤:
(I)根据病人的CT图像,确定靶区及紧要器官的轮廓,设置驻留点位置和参考点位置及处方剂量;
(I I)根据式(1)计算参考点的计算剂量:
(III)根据计算剂量与处方剂量的偏差、驻留位的时间构建评价函数式(2):
O(T)=(D0-AT)2+αTT·T (2)
(IV)要使得评价函数值最小,评价函数对驻留时间向量T进行微分,得到式(3):
(V)简化式(3)得到一个正规方程组,并用向量表示如下:
(B+αI)T=F (4)
(VI)将设置的α值代入式(4),即可求得驻留时间解集;
式中,dj表示第j个参考点的计算剂量,aij表示i驻留点在j参考点处的剂量率,驻留点数目为N,参考点数目为M。ti表示第i个驻留点的驻留时间,D表示参考点的计算剂量列向量,表示M个参考点的参考剂量列向量,α表示驻留时间平滑因子,B=ATA,F=ATD0,A为剂量率的M×N矩阵,AT是A的转置矩阵,I为单位矩阵,T是驻留时间的N维列向量。其中,D是通过式(1)计算得到的,D0是在计算之前设置的一个参考剂量值。
将式(1)用向量表示如下:
D=AT
其中,D=[d1,d2,…,dm]T
T=[t1,t2,…,tn]T。
进一步的,若计算得出的驻留时间的解为负数,则取值为0。
与现有技术相比,本发明具有以下有益效果:
本发明利用规则化方法获得后装驻留时间,同时,引入时间平滑因子解决了相邻驻留位中驻留时间相差大的问题,另一方面,有效地解决了因临床情况的复杂性,而导致驻留位较多、剂量参考点不能与驻留位一一对应等问题。
附图说明
图1是本发明中α=0时的驻留时间分布图。
图2是本发明中α=5时的驻留时间分布图。
图3是本发明中α=10时的驻留时间分布图。
具体实施方式
下面结合附图和实施例对本发明作进一步说明,本发明的实施方式包括但不限于下列实施例。
实施例
本实施例提供了一种利用规则化方法确定后装驻留时间的方法,包括以下步骤;
(一)根据病人的CT图像,确定靶区及紧要器官的轮廓,设置驻留点位置和参考点位置及处方剂量;(二)根据式(1)计算参考点的计算剂量:将上式用向量表示如下:D=AT,其中
D=[d1,d2,…,dm]T,T=[t1,t2,…,tn]T。
(三)根据计算剂量与处方剂量的偏差、驻留位的时间构建评价函数式(2):O(T)=(D0-AT)2+αTT·T;(四)评价函数对驻留时间向量T进行微分,得到式(3):(五)简化式(3)得到一个正规方程组,并用向量表示如下:(B+αI)T=F;(六)将设置的α值代入式(4),即可求得驻留时间解集。代入不同的α值,即可得到不同的驻留时间解集,如图1至3所示。
需要说明的是,因驻留时间不能为负,故若计算得出的驻留时间的解为负数,则取值为0。
测试例
采用本技术方案,一个单管施源器,设置16个驻留点,参考点位于靶区边界距离施源器1cm处,参考剂量设置为500cGy,分别设置α=0,α=5,α=10得到3组不同的驻留时间,如图1至3所示,由图可知,α=0时相邻位置的时间偏差最大为28.8s,α=5时相邻位置的时间偏差最大为12.5s,α=10时相邻位置的时间偏差最大为9.5s。
按照上述实施例,便可很好地实现本发明。值得说明的是,基于上述设计原理的前提下,为解决同样的技术问题,即使在本发明所公开的结构基础上做出的一些无实质性的改动或润色,所采用的技术方案的实质仍然与本发明一样,故其也应当在本发明的保护范围内。
Claims (3)
1.一种利用规则化方法确定后装源驻留时间的方法,其特征在于,包括以下步骤;
(I)根据病人的CT图像,确定靶区及紧要器官的轮廓,设置驻留点位置和参考点位置及处方剂量;
(II)根据式(1)计算参考点的计算剂量:
(III)根据计算剂量与处方剂量的偏差、驻留时间构建评价函数式(2):
O(T)=(D0-AT)2+αTT·T (2)
(IV)使得评价函数值最小,则评价函数对驻留时间向量T进行微分,得到式(3):
(V)简化式(3)得到一个正规方程组,并用向量表示如下:
(B+αI)T=F (4)
(VI)将设置的α值代入式(4),即可求得驻留时间解集;
式中,dj表示第j个参考点的计算剂量,aij表示i驻留点在j参考点处的剂量率,驻留点数目为N,参考点数目为M,ti表示第i个驻留点的驻留时间,D表示参考点的计算剂量列向量,表示M个参考点的参考剂量列向量,α表示驻留时间平滑因子,B=ATA,F=ATD0,A为剂量率的M×N矩阵,AT是A的转置矩阵,I为单位矩阵,T是驻留时间的N维列向量。
2.根据权利要求1所述的一种利用规则化方法确定后装源驻留时间的方法,其特征在于,将式(1)用向量表示如下:
D=AT
其中,D=[d1,d2,…,dm]T
T=[t1,t2,…,tn]T。
3.根据权利要求1所述的一种利用规则化方法确定后装源驻留时间的方法,其特征在于,若计算得出的驻留时间的解为负数,则取值为0。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410713531.5A CN104361253B (zh) | 2014-11-28 | 2014-11-28 | 一种利用规则化方法确定后装源驻留时间的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410713531.5A CN104361253B (zh) | 2014-11-28 | 2014-11-28 | 一种利用规则化方法确定后装源驻留时间的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104361253A true CN104361253A (zh) | 2015-02-18 |
CN104361253B CN104361253B (zh) | 2017-08-25 |
Family
ID=52528512
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410713531.5A Active CN104361253B (zh) | 2014-11-28 | 2014-11-28 | 一种利用规则化方法确定后装源驻留时间的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104361253B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106215334A (zh) * | 2016-08-23 | 2016-12-14 | 广州医科大学附属肿瘤医院 | 一种后装放疗剂量的逆向调强优化计算方法及系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100998497A (zh) * | 2006-12-29 | 2007-07-18 | 四川大学 | 确定电子束剂量分布的方法 |
WO2013112309A1 (en) * | 2012-01-24 | 2013-08-01 | Spectro, Inc. | Integrated, portable sample analysis system and method |
CN203436708U (zh) * | 2013-09-12 | 2014-02-19 | 中国医学科学院北京协和医院 | 一种柱状施源器治疗装置 |
-
2014
- 2014-11-28 CN CN201410713531.5A patent/CN104361253B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100998497A (zh) * | 2006-12-29 | 2007-07-18 | 四川大学 | 确定电子束剂量分布的方法 |
WO2013112309A1 (en) * | 2012-01-24 | 2013-08-01 | Spectro, Inc. | Integrated, portable sample analysis system and method |
CN203436708U (zh) * | 2013-09-12 | 2014-02-19 | 中国医学科学院北京协和医院 | 一种柱状施源器治疗装置 |
Non-Patent Citations (4)
Title |
---|
CHENG KONG ET AL.: "A new index comparable to BED for evaluating the biological efficacy of hypofractionated radiotherapy schemes in early stage non-small cell lung cancer:Analysis of data from the literature", 《LUNG CANCER》 * |
余辉等: "后装腔内治疗剂量计算的方法及临床应用", 《中国医疗设备》 * |
勾成俊等: "模拟动力学方法应用于调强放射治疗的研究", 《四川大学学报(自然科学版)》 * |
缪洪波: "中子刀治疗系统的软件研制", 《中国优秀硕士学位论文全文数据库 医药卫生科技辑》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106215334A (zh) * | 2016-08-23 | 2016-12-14 | 广州医科大学附属肿瘤医院 | 一种后装放疗剂量的逆向调强优化计算方法及系统 |
CN106215334B (zh) * | 2016-08-23 | 2019-04-16 | 广州医科大学附属肿瘤医院 | 一种后装放疗剂量的逆向调强优化计算方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN104361253B (zh) | 2017-08-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111068187B (zh) | 用于执行适应性放射疗法阶段的自动化工作流的系统 | |
CN104117151B (zh) | 一种在线自适应放疗计划优化方法 | |
EP3157627B1 (en) | System and method for automatic treatment planning | |
Pawlowski et al. | Reduction of dose delivered to organs at risk in prostate cancer patients via image-guided radiation therapy | |
CN104107062B (zh) | 一种评估放射治疗效果方法及系统 | |
CN114206438B (zh) | 使用投影图像预测放射治疗控制点 | |
Lee et al. | American Brachytherapy Society consensus guidelines for locally advanced carcinoma of the cervix. Part III: low-dose-rate and pulsed-dose-rate brachytherapy | |
CN104225809B (zh) | 一种带有呼吸补偿的4d放疗计划的实施方法及其设备 | |
EP3632508B1 (en) | Systems, methods and devices for automated target volume generation | |
CN109499014A (zh) | 一种妇科肿瘤后装手术辅助装置的制作方法 | |
Hyer et al. | The dosimetric impact of heterogeneity corrections in high-dose-rate 192Ir brachytherapy for cervical cancer: Investigation of both conventional Point-A and volume-optimized plans | |
EP3750595B1 (en) | Method and system for robust radiotherapy treatment planning for biological uncertainties | |
WO2015194554A1 (ja) | 放射線治療シミュレーション装置 | |
Petit et al. | Dose recalculation in megavoltage cone-beam CT for treatment evaluation: removal of cupping and truncation artefacts in scans of the thorax and abdomen | |
CN104361253A (zh) | 一种利用规则化方法确定后装源驻留时间的方法 | |
WO2021011499A1 (en) | Independent stereotactic radiotherapy dose calculation and treatment plan verification | |
Zhang et al. | Three dimensional expansion of margins for single-fraction treatments: Stereotactic radiosurgery brain cases | |
CN118574659A (zh) | 基本剂量计算 | |
US20230144962A1 (en) | Computer-implemented method for radiotherapy treatment planning, computer program product and computer system for performing the method | |
US20230141197A1 (en) | Computer-implemented method for radiotherapy treatment planning, computer program product and computer system for performing the method | |
Sonier et al. | Implementation of a volumetric modulated arc therapy treatment planning solution for kidney and adrenal stereotactic body radiation therapy | |
CN109999365B (zh) | 一种用于实现放疗逆向优化方法的计算设备和存储介质 | |
JP2020178935A (ja) | 線量分布判定システム、深層学習装置、線量分布判定方法及びコンピュータプログラム | |
Marqa et al. | A dose verification tool for high-dose-rate interstitial brachytherapy treatment planning in accelerated partial breast irradiation | |
Khamfongkhruea et al. | Dosimetric Evaluation of U-Net Deep Learning-Based Synthetic Computed Tomography for Adaptive Radiation Treatment Planning |
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 |