CN101477205B - 基于多算法的放射源反演方法 - Google Patents
基于多算法的放射源反演方法 Download PDFInfo
- Publication number
- CN101477205B CN101477205B CN2009101161160A CN200910116116A CN101477205B CN 101477205 B CN101477205 B CN 101477205B CN 2009101161160 A CN2009101161160 A CN 2009101161160A CN 200910116116 A CN200910116116 A CN 200910116116A CN 101477205 B CN101477205 B CN 101477205B
- Authority
- CN
- China
- Prior art keywords
- radioactive source
- calculation
- photon
- information
- power spectrum
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Measurement Of Radiation (AREA)
Abstract
本发明公开了一种放射源反演方法。即利用水箱或者其他辐射剂量学测量设备测得PDD曲线,通过本方法反演得到高精度的放射源信息——光子能谱,或者质子能谱,或者电子能谱与光子污染等放射源信息。使用本发明具有如下优点:可同时实现重建放射源的光子能谱、电子能谱、质子能谱与光子污染等多种放射源信息;克服传统能谱计算方法不能获得精确多种放射源信息的缺点;获得精确的放射源能谱的同时,还可以获得光子污染等信息;为剂量计算提供正确的放射源信息;可以获得离散的或者连续的能谱与光子污染能谱信息;采用多种算法实现精确的放射源,克服传统单一算法不能保证能谱反演有解的问题。
Description
技术领域
本发明属于放射源反演与重建方法领域。
背景技术
利用体外测量的剂量信息可以反推源的信息就是源反演。源反演方法保证得到正确的放射源信息——包括电子能谱、质子能谱、光子能谱、光子污染等信息,这为了精确了解放射源的射线成分构成提供了重要保证。
精确的放射源信息是提高辐射剂量计算精度的重要参数,不管是蒙特卡罗模拟或者解析剂量计算,放射源信息都起到重要作用。目前,获取放射源信息的途径主要有三种,其一是蒙特卡罗模拟源,但要根据具体源建立准确几何和物理模型,缺乏灵活性;其二采用康普顿散射仪(Compton Scatter Spectrometer)等实验测量,但耗费大不便于测量,并且所得结果因不能得到光子污染情况而不能应用于剂量计算;最后是利用体外测量剂量数据——如测量水箱中的百分深度剂量(Percentage DepthDose,PDD)曲线反演能谱,这种方法简单,可操作性强;只要建立适当的数学模型就可以实现,而且PDD曲线的测量,普通的辐射测量设备都能实现。在利用PDD反演能谱中,本专利的方法测试表明了本专利的正确性与可靠性,并同时表明忽略了光子污染的能谱反演所得到的能谱是错误的。传统方法忽略了光子污染,并且方法单一,不能保证所得放射源信息的正确性。本专利将通过发展新方法实现精确多算法的多种放射源反演,主要实现重建放射源的光子能谱、质子能谱、电子能谱与光子污染等多种放射源信息。
发明内容
本发明的目的在于提供一种基于多算法的放射源反演方法。
本发明的技术方案如下:
基于多算法的放射源反演方法,其特征在于通过如下步骤实现:
第一步,获得单能深度剂量曲线PPD,单能深度剂量曲线是通过国际公开使用的蒙特卡罗(例如EGSnrc)等程序模拟单向垂直入射水箱或者其他测量介质的放射源——放射源包括:放射性同位素产生的α、β、γ射线,以及各类X射线机或由加速器产生的X射线、电子线、质子束等,在10cm×10cm射野或者其他射野大小能量沉积获得——由于电子平衡的影响,射野较大会更好,而小野会带来一定的计算误差;能量间隔从0到60MeV或者更大能量;
第二步,利用三维水箱或者其他测量介质测量获取百分深度剂量曲线:把放射源机头和机架角度调整到0度,水箱或者其他测量介质放置在水平地面,中心点与机头中心点连线垂直于水平地面,机头开野大小设置为10cm×10cm射野或者其他射野大小——由于电子平衡的影响,射野较大会更好,而小野会带来一定的计算误差;利用放射源的电离室探头或者其他测量仪器按照一定的深度间隔扫描;
第三步,放射源反演;
(3)采用如下的数学模型实现:
其中,
D′(zj)为根据反演出来的能谱重建出来的PDD曲线;
σ为均方根误差(Root Mean Squared Error);
m为用于拟合时的测量PDD数据个数;
N+4为需要拟合的系数个数;N为能群的个数;
D(Ei,z)为能量为Ei的深度剂量曲线;
能谱为Ф(E)=|ai|;其中ai可以是常系数也可以是带有待定系数的经验公式;
Dc为归一化系数;
μp为伴随射线入射的高能散射光平均衰减系数;
μe为在轫致辐射光子和次级光子衰减系数;
v为与入射表面剂量有关的系数。
(4)放射源信息反演:通过利用传统成熟的非线性反演算法包括:C1:Levenberg-Marquardt;C2:Quasi-Newton;C3:Gradient;C4:Conjugate-Gradient;C5:Newton;C6:Principal-Axis;C7:NMinimize等算法,根据已知的测量PDD曲线D(z)和单能PDD曲线D(E,z),进行数据拟合求得能谱以及光子污染分布函数
第四步,多算法结果比较;根据以上提到的反演算法的计算结果:能谱与光子污染,以及对应的均方根误差,综合评价最优的计算结果;综合评价方法如下:
测量PDD与根据反演所得能谱反演计算D(z)对比;选取计算D(z)与测量PDD在“建成”区附近的均方根误差比较小,采用公式(1)计算误差,并且总的均方根误差比较小的计算结果;
第五步,最优结果输出与显示;将以上比较的最优结果同时以图像以及数据的形式输出。
所述的计算结果,其特征在于具有同时实现重建多种放射源能谱与光子污染——放射源包括:放射性同位素产生的α、β、γ射线,以及各类X射线机或由加速器产生的X射线、电子线、质子束等。
所述的计算结果,其特征在于采用多种算法实现,并能够选取最优的放射源。
所述的计算结果,其特征在于能重建多种放射源能谱的同时,还可以获得光子污染等信息。
所述的计算结果,其特征在于可以获得离散的或者连续的能谱与光子污染信息。
经过实验表明,本专利可以得到高精度的放射源信息——光子能谱、质子能谱、电子能谱与光子污染等多种放射源信息。
使用本发明具有如下优点:
(1)可同时实现重建放射源的光子能谱、电子能谱、质子能谱与光子污染等多种放射源信息;
(2)克服传统能谱计算方法不能获得精确多种放射源信息的缺点;
(3)获得精确的放射源能谱的同时,还可以获得光子污染等信息。为剂量计算提供正确的放射源信息;
(4)可以获得离散的或者连续的能谱与光子污染能谱信息;
(5)采用多种算法实现精确的放射源,克服传统单一算法不能保证能谱反演有解的问题。
附图说明
图1是本发明结构示意图;
图2是测量PDD的测量几何模型;
图3是电子能谱反演输入:测量电子PDD曲线;
图4是电子能谱反演结果与原能谱(点线)的对比;
图5是光子能谱反演输入:测量光子PDD曲线
图6是光子能谱反演结果与原能谱(点线)的对比。
具体实施方式
基于多算法的放射源反演方法
实例1:以Varian2100放射源(电子加速器)输出的20MeV电子源为例。
第一步,获得单能深度剂量曲线PPD,作为单能PPD数据库;单能深度剂量曲线是通过国际公开使用的蒙特卡罗程序EGSnrc模拟单向垂直入射的射线在10cm×10cm射野能量沉积获得,能量间隔从0到60MeV,能量间隔为0.25MeV;
利用三维水箱测量获取测量PDD,如图2为测量模型。利用自动水箱测量Varian加速器PDD曲线:把放射源的机头和机架角度调整到0度,水箱放置在水平地面,水箱中心点与机头中心点连线垂直于水平地面。机头射野大小设置为10cm×10cm。利用电离室探头按照取0.125cm深度间隔扫描,得到如图3的测量PDD曲线。
第二步,放射源反演。运行基于公式(1)的数学模型的程序(计算步骤如图1),反演得到能谱、光子污染函数以及对应均方根误差。并通过选择不同的算法,得到不同算法对应的能谱、光子污染函数以及对应的均方根误差。
(1)放射源的能谱与光子污染信息,采用如下的数学模型实现:
其中,
D′(zj)为根据反演出来的能谱重建出来的PDD曲线;
σ为均方根误差(Root Mean Squared Error);
m为用于拟合时的测量PDD数据个数;
N+4为需要拟合的系数个数;N为能群的个数;
D(Ei,z)为能量为Ei的深度剂量曲线;
能谱为Ф(E)=|ai|;其中ai可以是常系数也可以是带有待定系数的经验公式,本例子采用常系数;
Dc为归一化系数;
μp为伴随射线入射的高能散射光平均衰减系数;
μe为轫致辐射光子和次级光子衰减系数;
v为与入射表面剂量有关的系数。
(2)放射源信息反演:本例通过利用传统成熟的非线性反演算法包括:C1:Levenberg-Marquardt;C2:Quasi-Newton;C3:Gradient;C4:Conjugate-Gradient;C5:Newton;C6:Principal-Axis;C7:NMinimize算法,根据已知的测量PDD曲线D(z)和单能PDD曲线D(E,z),进行数据拟合求得能谱以及光子污染分布函数
第三步,根据均方根误差,选择最优的计算结果。
第四步,最优结果输出与显示。如图表1与图4所示,其中C1算法的结果为最优结果。
表1是反演参数。
表1
其中,第一栏的Dc、μp、μe和v值是根据原光子污染函数的拟合结果
实例2:
以Varian2100放射源(电子加速器打靶产生的光子束)输出的6MV光子源为例。
第一步,在实例1的基础上,把图3的平行源改为光子点源入射,用户输入的测量PDD曲线见图5,其他模型参数不变。第二步到第四步与实例1的步骤一致。输出的光子能谱结果见图6,其中C1/C2/C3/C6的算法结果较好,而C1算法的结果为最优结果。
制作原理:用户利用三维水箱测量获取PDD。利用自动水箱或者其他介质测量放射源PDD曲线:把放射源机头和机架角度调整到0度,水箱或者其他介质放置在水平地面,测量装置中心点与放射源机头中心点连线垂直于水平地面。放射源机头开野大小设置为10cm×10cm或者其他射野大小——由于电子平衡的影响,射野较大会更好,而小野会带来一定的计算误差。利用电离室探头或者其他探测器按照一定的深度间隔扫描——深度间隔可以根据用户具备的电离室条件调整。用户将所得的PDD数据按照“第一列是坐标值,第二列是对应的PDD值”格式作为输入文件;然后,运行本方法分析程序,即可输出最优的放射源信息——光子能谱、或者质子能谱、或者电子能谱与光子污染等还有其他放射源信息。
Claims (5)
1.基于多算法的放射源反演方法,其特征在于通过如下步骤实现:
第一步,获得单能深度剂量曲线,单能深度剂量曲线是通过国际公开使用的蒙特卡罗程序模拟单向垂直入射水箱或者其他测量介质的放射源--放射源包括:放射性同位素产生的α、β、γ射线,以及各类X射线机或由加速器产生的X射线、电子线、质子束,在10cm×10cm射野或者其他射野大小能量沉积获得--由于电子平衡的影响,射野较大会更好,而小野会带来一定的计算误差;能量间隔从0到60MeV或者更大能量;
第二步,利用三维水箱或者其他测量介质测量获取百分深度剂量曲线:把放射源机头和机架角度调整到0度,水箱或者其他测量介质放置在水平地面,中心点与机头中心点连线垂直于水平地面,机头开野大小设置为10cm×10cm射野或者其他射野大小--由于电子平衡的影响,射野较大会更好,而小野会带来一定的计算误差;利用放射源的电离室探头或者其他测量仪器按照一定的深度间隔扫描;
第三步,放射源反演;
(1)采用如下的数学模型实现:
其中,
D′(zj)为根据反演出来的能谱重建出来的百分深度剂量曲线曲线;
σ为均方根误差(Root Mean Squared Error);
m为用于拟合时的测量百分深度剂量曲线数据个数;
N+4为需要拟合的系数个数;N为能群的个数;
D(Ei,z)为能量为Ei的深度剂量曲线;
能谱为Φ(E)=|ai|;其中ai是常系数或者是带有待定系数的经验公式;
Dc为归一化系数;
μp为伴随射线入射的高能散射光平均衰减系数;
μe为在轫致辐射光子和次级光子衰减系数;
v为与入射表面剂量有关的系数。
(2)放射源信息反演:通过利用传统成熟的非线性反演算法包括:C1:Levenberg-Marquardt;C2:Quasi-Newton;C3:Gradient;C4:Conjugate-Gradient;C5:Newton;C6:Principal-Axis;C7:NMinimize,根据已知的测量百分深度剂量曲线D(z)和单能百分深度剂量曲线D(E,z),进行数据拟合求得能谱以及光子污染分布函数
第四步,多算法结果比较;根据以上提到的反演算法的计算结果:能谱与光子污染,以及对应的均方根误差,综合评价最优的计算结果;综合评价方法如下:
测量百分深度剂量曲线与根据反演所得能谱反演计算D(z)对比;选取计算D(z)与测量百分深度剂量曲线在“建成”区附近的均方根误差比较小,采用公式(1)计算误差,并且总的均方根误差比较小的计算结果;
第五步,最优结果输出与显示;将以上比较的最优结果同时以图像以及数据的形式输出。
2.根据权利要求1所述的基于多算法的放射源反演方法,其特征在于所述的计算结果,具有同时实现重建多种放射源能谱与光子污染--放射源包括:放射性同位素产生的α、β、γ射线,以及各类X射线机或由加速器产生的X射线、电子线、质子束。
3.根据权利要求1所述的基于多算法的放射源反演方法,其特征在于所述的计算结果,采用多种算法实现,并能够选取最优的放射源。
4.根据权利要求1所述的基于多算法的放射源反演方法,其特征在于所述的计算结果,能重建多种放射源能谱的同时,还获得光子污染信息。
5.根据权利要求1所述的基于多算法的放射源反演方法,其特征在于所述的计算结果,能够获得离散的或者连续的能谱与光子污染信息。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009101161160A CN101477205B (zh) | 2009-01-22 | 2009-01-22 | 基于多算法的放射源反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009101161160A CN101477205B (zh) | 2009-01-22 | 2009-01-22 | 基于多算法的放射源反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101477205A CN101477205A (zh) | 2009-07-08 |
CN101477205B true CN101477205B (zh) | 2012-07-18 |
Family
ID=40837950
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009101161160A Expired - Fee Related CN101477205B (zh) | 2009-01-22 | 2009-01-22 | 基于多算法的放射源反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101477205B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102033239B (zh) * | 2009-09-28 | 2013-07-10 | 同方威视技术股份有限公司 | 加速器x射线能量测量系统 |
CN102353975B (zh) * | 2011-08-31 | 2013-06-12 | 珠海和佳医疗设备股份有限公司 | 碘[125i]密封源水中吸收剂量的测量及修正方法 |
CN102426377A (zh) * | 2011-09-06 | 2012-04-25 | 中国科学院合肥物质科学研究院 | 一种三维剂量反演方法 |
JP2016061584A (ja) * | 2014-09-16 | 2016-04-25 | キヤノン株式会社 | 放射線計測装置 |
CN105079984B (zh) * | 2015-07-06 | 2017-10-24 | 上海联影医疗科技有限公司 | 一种源压缩和解压缩的方法及装置 |
CN106199672B (zh) * | 2016-06-30 | 2019-01-08 | 中国科学院合肥物质科学研究院 | 一种基于蒙特卡罗光子模拟的卷积叠加剂量计算方法 |
CN106932810B (zh) * | 2017-04-01 | 2018-02-23 | 西安一体医疗科技有限公司 | 一种伽玛射线剂量的卷积计算方法 |
CN108333621A (zh) * | 2018-02-01 | 2018-07-27 | 中国计量科学研究院 | 加速器辐射质测量方法 |
CN116134434A (zh) * | 2020-08-13 | 2023-05-16 | 上海联影医疗科技股份有限公司 | 辐射源建模系统和方法 |
CN112231899B (zh) * | 2020-09-30 | 2023-02-03 | 天津医科大学 | X射线连续能谱的确定方法、装置、电子设备及存储介质 |
CN112419439B (zh) * | 2020-12-04 | 2023-11-21 | 兰州大学 | 一种治疗计划中质子重离子能量快速反演的方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6301329B1 (en) * | 1998-02-09 | 2001-10-09 | The University Of Southampton | Treatment planning method and apparatus for radiation therapy |
CN201051145Y (zh) * | 2007-05-11 | 2008-04-23 | 北京大北华光科技发展有限责任公司 | 三维剂量测量装置 |
CN100432700C (zh) * | 2006-12-29 | 2008-11-12 | 成都川大奇林科技有限责任公司 | 医用电子加速器能谱测量方法 |
-
2009
- 2009-01-22 CN CN2009101161160A patent/CN101477205B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6301329B1 (en) * | 1998-02-09 | 2001-10-09 | The University Of Southampton | Treatment planning method and apparatus for radiation therapy |
CN100432700C (zh) * | 2006-12-29 | 2008-11-12 | 成都川大奇林科技有限责任公司 | 医用电子加速器能谱测量方法 |
CN201051145Y (zh) * | 2007-05-11 | 2008-04-23 | 北京大北华光科技发展有限责任公司 | 三维剂量测量装置 |
Non-Patent Citations (2)
Title |
---|
Ernesto Mainegra-Hing et al.Calculation of photon energy deposition kernels and electron dose point kernels in water.《Medical Physics》.2005,第32卷(第3期),685-699. * |
江海燕 等.BNCT中载能粒子对肿瘤细胞损伤效应的MonteCarlo模拟及分析.《核技术》.2004,第27卷(第9期),687-692. * |
Also Published As
Publication number | Publication date |
---|---|
CN101477205A (zh) | 2009-07-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101477205B (zh) | 基于多算法的放射源反演方法 | |
Dudouet et al. | Double-differential fragmentation cross-section measurements of 95 MeV/nucleon 12 C beams on thin targets for hadron therapy | |
Sterpin et al. | Analytical computation of prompt gamma ray emission and detection for proton range verification | |
Yamaguchi et al. | Imaging of monochromatic beams by measuring secondary electron bremsstrahlung for carbon-ion therapy using a pinhole x-ray camera | |
Taddeucci et al. | Multiple-scattering corrections to measurements of the prompt fission neutron spectrum | |
Zarifi et al. | Characterization of prompt gamma-ray emission with respect to the Bragg peak for proton beam range verification: A Monte Carlo study | |
CN101477202B (zh) | 一种能量沉积核的获取方法 | |
Haj-Heidari et al. | Method for developing HPGe detector model in Monte Carlo simulation codes | |
CN103091699A (zh) | 散射法测量强γ射线能谱的装置及方法 | |
CN101477203A (zh) | 一种解析蒙特卡罗剂量计算方法 | |
Thomas et al. | What is neutron metrology and why is it needed? | |
Wrońska et al. | Prompt-gamma emission in GEANT4 revisited and confronted with experiment | |
Birattari et al. | A Bonner Sphere Spectrometer with extended response matrix | |
Bedogni et al. | Testing a newly developed single-sphere neutron spectrometer in reference monochromatic fields from 147 keV to 14.8 MeV | |
CN112906193A (zh) | 一种探测器γ响应函数的高斯展宽及解谱方法 | |
Wu et al. | Experimental validation of a BNCT epithermal neutron flux detector using 55Mn (n, γ) 56Mn reaction at an accelerator-based BNCT facility | |
Xie et al. | Calibration and unfolding of the pulse height spectra of liquid scintillator-based neutron detectors using photon sources | |
Kulisek et al. | Analysis of an indirect neutron signature for enhanced UF6 cylinder verification | |
CN113466915A (zh) | 一种基于高纯锗探测器的辐射剂量率测量系统及方法 | |
MacMullin et al. | Measurement of the elastic scattering cross section of neutrons from argon and neon | |
Shu et al. | Novel method exploration of monitoring neutron beam using Cherenkov photons in BNCT | |
Liu et al. | Performance evaluation of the source description of the THOR BNCT epithermal neutron beam | |
Wernli et al. | Neutron dosimetry around accelerators in Switzerland | |
Wang et al. | Experimental and theoretical study of long counters on the departure of “point” assumption and scattering background influence | |
Röttger et al. | The PTB neutron reference fields (PIAF)-quasi-monoenergetic neutron reference fields in the energy range from thermal to 200 MeV |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120718 Termination date: 20150122 |
|
EXPY | Termination of patent right or utility model |