CN116500665A - 傅里叶尺度变换特性的微剂量探测器组织等效换算方法 - Google Patents

傅里叶尺度变换特性的微剂量探测器组织等效换算方法 Download PDF

Info

Publication number
CN116500665A
CN116500665A CN202310561795.2A CN202310561795A CN116500665A CN 116500665 A CN116500665 A CN 116500665A CN 202310561795 A CN202310561795 A CN 202310561795A CN 116500665 A CN116500665 A CN 116500665A
Authority
CN
China
Prior art keywords
deposition energy
micro
energy
dose detector
equivalent conversion
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
CN202310561795.2A
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.)
Chengdu Technological University CDTU
Original Assignee
Chengdu Technological University CDTU
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 Chengdu Technological University CDTU filed Critical Chengdu Technological University CDTU
Priority to CN202310561795.2A priority Critical patent/CN116500665A/zh
Publication of CN116500665A publication Critical patent/CN116500665A/zh
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/29—Measurement performed on radiation beams, e.g. position or section of the beam; Measurement of spatial distribution of radiation
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/16—Measuring radiation intensity
    • G01T1/161—Applications in the field of nuclear medicine, e.g. in vivo counting
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00—Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/36—Measuring spectral distribution of X-rays or of nuclear radiation spectrometry
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01T—MEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T7/00—Details of radiation-measuring instruments
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E30/00—Energy generation of nuclear origin
    • Y02E30/30—Nuclear fission reactors

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • General Physics & Mathematics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Molecular Biology (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Biomedical Technology (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • Measurement Of Radiation (AREA)

Abstract

本发明公开了基于傅里叶尺度变换特性的微剂量探测器组织等效换算方法，属于辐射剂量测量技术领域，包括：获取给定条件下入射辐射在微剂量探测器中的沉积能量分布D(x)，并计算入射辐射在微剂量探测器和生物组织中的平均沉积能量比；将D(x)进行尺度变换，并将D(x)、进行尺度变换后得到的H(x)进行傅里叶变换，得到D(x)、H(x)分别在傅里叶空间中的模MD、MH；根据平均沉积能量比对模MD、MH差值进行修正，然后将修正后的差值进行傅里叶逆变换、归一化处理；最后将归一化处理的结果与H(x)进行卷积，得到生物组织中的沉积能谱T(x)。本发明能够准确可靠的将微剂量探测器测得的沉积能量谱转换为生物组织中沉积能量谱，有效改善当前组织等效换算方法在低能区换算时存在偏差的问题，实现测量目标。

Description

傅里叶尺度变换特性的微剂量探测器组织等效换算方法
技术领域
本发明涉及辐射剂量测量与辐射防护技术领域，尤其涉及一种基于傅里叶尺度变换特性的微剂量探测器组织等效换算方法。
背景技术
固体微剂量探测器由于能直接将探测器灵敏体积制成生物组织大小，可以准确表征入射辐射在微体积内的相互作用，在近年来引发了广泛的关注。不过，由于探测器材料(主要为硅和金刚石)与生物组织材料之间存在差值，如物理特性、密度、电离能等，致使放射性粒子或射线在微剂量探测器中的沉积能量分布与生物组织中的沉积能量分布并不一致。因此，为获得生物组织微观体积内的沉积能量分布，需要对微剂量探测器的测量结果进行组织等效换算。
当前，关于微剂量探测器测得的沉积能谱向同等大小生物组织的换算，有学者提出了一种基于谱图分布变换的换算方法。由于在换算过程中未对低能沉积事件的特点加以考虑，致使低能区的换算结果并不理想，即无法有效将微剂量探测器测得的沉积能量谱转换为生物组织中沉积能量谱，特别在低能区换算时存在较大误差。
发明内容
本发明的目的在于克服现有技术的问题，提供了一种基于傅里叶尺度变换特性的微剂量探测器组织等效换算方法。
本发明的目的是通过以下技术方案来实现的：傅里叶尺度变换特性的微剂量探测器组织等效换算方法，该方法包括以下步骤：
获取给定条件(一定能量射线粒子辐照场)下入射辐射在微剂量探测器中的沉积能量分布D(x)；
计算当前条件下入射辐射在微剂量探测器和生物组织中的平均沉积能量比λ；
根据平均沉积能量比λ对测得的沉积能量分布D(x)进行尺度变换，得到尺度变换后的沉积能量分布H(x)；
将D(x)、H(x)进行傅里叶变换，得到D(x)、H(x)分别在傅里叶空间中的模MD、MH；
计算模MD和MH之间的差值，并根据平均沉积能量比对差值进行修正得到ME；
将修正后的差值ME进行傅里叶逆变换得到沉积能量分布差E(x)，再对E(x)进行归一化处理得到C(x)；
将归一化处理的结果C(x)与尺度变换后的沉积能量分布H(x)进行卷积，得出同等条件下与微剂量探测器大小相同的生物组织中的沉积能谱T(x)。
在一示例中，所述平均沉积能量比λ的计算式为：
式中，ε表示入射辐射在介质中的沉积能量；ε(Ein,微剂量探测器)表示入射能量为Ein时，微剂量探测器中的沉积能量；ε(Ein,组织)表示入射能量为Ein时，生物组织中的沉积能量。
在一示例中，粒子平均沉积能量ε的计算式为：
ΔE＝Ei-Ei+1
ε＝E0-En
式中，L表示入射辐射在介质(微剂量探测器或生物组织)中穿过的距离；ΔE表示能量间隔；E0表示入射辐射进入介质之前的能量；En表示入射辐射离开介质时的能量；i表示第i个能量间隔点；S(Ei)表示能量大小为Ei时，介质对入射粒子的线性阻止本领。
在一示例中，对测得的沉积能量分布D(x)进行尺度变换的计算式为：
H(x)＝λD(λx)
式中，x表征粒子沉积能量。
在一示例中，计算模MD和MH之间的差值，并根据平均沉积能量比对差值进行修正的计算式为：
在一示例中，沉积能谱T(x)的计算式为：
T(x)＝C(x)*H(x)
式中，x表征粒子沉积能量。
在一示例中，傅里叶变换、傅里叶空间中模的计算公式为：
式中，1≤k≤N，表示粒子沉积能量数据点x在傅里叶空间对应的数据点，N表示数据长度，x表征粒子沉积能量。
在一示例中，傅里叶逆变换计算公式为：
式中，1≤k≤N，表示粒子沉积能量数据点x在傅里叶空间对应的数据点，N表示数据维度，x表征粒子沉积能量。
在一示例中，剔除沉积能量分布差E(x)中的负值。
在一示例中，剔除E(x)负值后进行归一化处理的计算式为：
式中，E1(x)表示剔除E(x)分布中负值后得到的分布；C(x)表示归一化后的分布，x表征粒子沉积能量。
需要进一步说明的是，上述方法各示例对应的技术特征可以相互组合或替换构成新的技术方案。
与现有技术相比，本发明有益效果是：
本发明方法能够广泛应用于固体微剂量探测器的组织等效修正中，使探测器适用于粒子全射程内测量；结合仿真结果可知，本发明方法能够准确可靠的将微剂量探测器测得的沉积能量谱转换为生物组织中沉积能量谱，有效改善当前组织等效换算方法在低能区换算时存在偏差的问题，实现测量目标；整个过程操作简便，在一定程度上降低测量成本。
附图说明
下面结合附图对本发明的具体实施方式作进一步详细的说明，此处所说明的附图用来提供对本申请的进一步理解，构成本申请的一部分，在这些附图中使用相同的参考标号来表示相同或相似的部分，本申请的示意性实施例及其说明用于解释本申请，并不构成对本申请的不当限定。
图1为50MeV质子在骨骼体膜中的百分深度剂量曲线图；
图2为本发明公开的一种基于傅里叶尺度变换特性的微剂量探测器组织等效换算方法的流程图；
图3为百分深度剂量曲线坪区金刚石微剂量探测器的等效换算结果与骨骼中沉积能谱的对比图；
图4为百分深度剂量曲线布拉格峰前端处金刚石微剂量探测器的等效换算结果与骨骼中沉积能谱的对比图；
图5为百分深度剂量曲线布拉格峰后端处金刚石微剂量探测器的等效换算结果与骨骼中沉积能谱的对比图。
具体实施方式
下面结合附图对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
在本发明的描述中，需要说明的是，属于“中心”、“上”、“下”、“左”、“右”、“竖直”、“水平”、“内”、“外”等指示的方向或位置关系为基于附图所述的方向或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。此外，使用序数词(例如，“第一和第二”、“第一至第四”等)是为了对物体进行区分，并不限于该顺序，而不能理解为指示或暗示相对重要性。
在本发明的描述中，需要说明的是，除非另有明确的规定和限定，属于“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本发明中的具体含义。
此外，下面所描述的本发明不同实施方式中所涉及的技术特征只要彼此之间未构成冲突就可以相互结合。
本实施例以10×10×10μm3大小金刚石微剂量探测器在50MeV质子辐照场中测得的能谱向同等大小生物组织中沉积能量等效换算为例，采用本发明方法对金刚石微剂量探测器在骨骼百分深度剂量曲线坪区、布拉格峰前端、布拉格峰后端测得的能量沉积分布进行组织等效换算。其中，50MeV质子照射下骨骼的百分深度剂量曲线如图1所示，它表示射线粒子(50MeV质子)在骨骼中某一深度处的吸收剂量与参考深度处吸收剂量之比的百分数，是描述射线中心轴不同深度处相对剂量分布的物理量。
此时，如图2所示，本发明一种基于傅里叶尺度变换特性的微剂量探测器组织等效换算方法包括以下步骤：
S1：获取金刚石微剂量探测器在百分深度剂量曲线坪区、布拉格峰前端、布拉格峰后端处的沉积能量分布D1(x)、D2(x)、D3(x)；
本实施例中，选取的布拉格曲线坪区、布拉格峰前端、布拉格峰后端测量点的三个测量点在体膜中的具体深度参数如表1所示：
表1选取的测量点在骨骼体膜中的深度位置
S2：计算在测量点1、测量点2、测量点3处入射能量下金刚石微剂量探测器与骨骼之间的平均沉积能量比，λ1、λ2、λ3；
S3：根据计算得到的平均沉积能量比将测量点1、测量点2、测量点3处测得的沉积能量分布D1(x)、D2(x)、D3(x)进行尺度变换，得到H1(x)、H2(x)、H3(x)；
S4：将分布D1(x)、D2(x)、D3(x)和H1(x)、H2(x)、H3(x)进行傅里叶变换并求得变换后结果在傅里叶空间里的模MD1、MD2、MD3和MH1、MH2、MH3；
S5：计算模MD1和MH1、MD2和MH2、MD3和MH3之间的差异并根据平均沉积能量比对差异进行修正得到ME1、ME2、ME3；
S6：将ME1、ME2、ME3进行傅里叶逆变换得到E1(x)、E2(x)、E3(x)，并对E1(x)、E2(x)、E3(x)进行归一化处理得到C1(x)、C2(x)、C3(x)；
S7：将归一化处理的结果C1(x)、C2(x)、C3(x)分别与尺度变换后的沉积能量分布H1(x)、H2(x)、H3(x)进行卷积计算，得出测量点1、测量点2、测量点3处骨骼中的沉积能量分布T1(x)、T2(x)、T3(x)。
为验证本发明方法的组织等效换算效果，通过软件模拟获取各测量点处质子在相同大小骨骼中的沉积能量分布B1(x)、B2(x)、B3(x)，将运用本发明提出的方法得到的沉积能量分布T1(x)、T2(x)、T3(x)与B1(x)、B2(x)、B3(x)进行对比，结果如图3-5所示，图3-5中横坐标ε表示骨骼中沉积的能量，纵坐标F(ε)表示沉积能量的概率分布。其中图3为测量点1处金刚石微剂量探测器的等效换算结果与骨骼中沉积能量分布的对比；图4为测量点2处金刚石微剂量探测器的等效换算与骨骼中沉积能量分布的对比；图5为测量点3处金刚石微剂量探测器的等效换算结果与骨骼中沉积能量分布的对比。综合上述图3-5可以看出，通过本发明方法对微剂量探测器进行组织等效换算后得到的沉积能量分布与相同条件下骨骼中的沉积能量分布吻合度较好，尤其是在测量点1、2处低能区，换算结果与实际结果之间的重合度较高，即：本发明方法能将探测器测得的沉积能量分布准确转换为相同条件下生物组织体中的沉积能量分布，且与现有换算方法相比，该方法可以有效消除低能区的换算偏差，提高换算结果与实际结果之间的吻合度。
以上具体实施方式是对本发明的详细说明，不能认定本发明的具体实施方式只局限于这些说明，对于本发明所属技术领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干简单推演和替代，都应当视为属于本发明的保护范围。

Claims (10)

1.一种傅里叶尺度变换特性的微剂量探测器组织等效换算方法，其特征在于：其包括以下步骤：
获取给定条件下入射辐射在微剂量探测器中的沉积能量分布D(x)；
计算当前条件下入射辐射在微剂量探测器和生物组织中的平均沉积能量比λ；
根据平均沉积能量比λ对测得的沉积能量分布D(x)进行尺度变换，得到尺度变换后的沉积能量分布H(x)；
将D(x)、H(x)进行傅里叶变换，得到D(x)、H(x)分别在傅里叶空间中的模MD、MH；
计算模MD和MH之间的差值，并根据平均沉积能量比对差值进行修正得到ME；
将修正后的差值ME进行傅里叶逆变换得到沉积能量分布差E(x)，再对E(x)进行归一化处理得到C(x)；
将归一化处理的结果C(x)与尺度变换后的沉积能量分布H(x)进行卷积，得出同等条件下与微剂量探测器大小相同的生物组织中的沉积能谱T(x)。
2.根据权利要求1所述的傅里叶尺度变换特性的微剂量探测器组织等效换算方法，其特征在于：所述平均沉积能量比λ的计算式为：
式中，ε表示入射辐射在介质中的沉积能量；ε(Ein,微剂量探测器)表示入射能量为Ein时，微剂量探测器中的沉积能量；ε(Ein,组织)表示入射能量为Ein时，生物组织中的沉积能量。
3.根据权利要求1所述的傅里叶尺度变换特性的微剂量探测器组织等效换算方法，其特征在于：粒子平均沉积能量ε的计算式为：
ΔE＝Ei-Ei+1
ε＝E0-En
式中，L表示入射辐射在介质中穿过的距离；ΔE表示能量间隔；E0表示入射辐射进入介质之前的能量；En表示入射辐射离开介质时的能量；i表示第i个能量间隔点；S(Ei)表示能量大小为Ei时，介质对入射粒子的线性阻止本领。
4.根据权利要求1所述的傅里叶尺度变换特性的微剂量探测器组织等效换算方法，其特征在于：对测得的沉积能量分布D(x)进行尺度变换的计算式为：
H(x)＝λD(λx)
式中，x表征粒子沉积能量。
5.根据权利要求1所述的傅里叶尺度变换特性的微剂量探测器组织等效换算方法，其特征在于：计算模MD和MH之间的差值，并根据平均沉积能量比对差值进行修正的计算式为：
6.根据权利要求1所述的傅里叶尺度变换特性的微剂量探测器组织等效换算方法，其特征在于：沉积能谱T(x)的计算式为：
T(x)＝C(x)*H(x)
式中，x表征粒子沉积能量。
7.根据权利要求1所述的傅里叶尺度变换特性的微剂量探测器组织等效换算方法，其特征在于：傅里叶变换、傅里叶空间中模的计算公式为：
式中，1≤k≤N，表示粒子沉积能量数据点x在傅里叶空间对应的数据点，N表示数据长度，x表征粒子沉积能量。
8.根据权利要求1所述的傅里叶尺度变换特性的微剂量探测器组织等效换算方法，其特征在于：傅里叶逆变换计算公式为：
式中，1≤k≤N，表示粒子沉积能量数据点x在傅里叶空间对应的数据点，N表示数据维度，x表征粒子沉积能量。
9.根据权利要求1所述的傅里叶尺度变换特性的微剂量探测器组织等效换算方法，其特征在于：剔除沉积能量分布差E(x)中的负值。
10.根据权利要求9所述的傅里叶尺度变换特性的微剂量探测器组织等效换算方法，其特征在于：剔除E(x)负值后进行归一化处理的计算式为：
式中，E1(x)表示剔除E(x)分布中负值后得到的分布；C(x)表示归一化后的分布，x表征粒子沉积能量。
CN202310561795.2A 2023-05-17 2023-05-17 傅里叶尺度变换特性的微剂量探测器组织等效换算方法 Pending CN116500665A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310561795.2A CN116500665A (zh) 2023-05-17 2023-05-17 傅里叶尺度变换特性的微剂量探测器组织等效换算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310561795.2A CN116500665A (zh) 2023-05-17 2023-05-17 傅里叶尺度变换特性的微剂量探测器组织等效换算方法

Publications (1)

Publication Number Publication Date
CN116500665A true CN116500665A (zh) 2023-07-28

Family

ID=87318230

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310561795.2A Pending CN116500665A (zh) 2023-05-17 2023-05-17 傅里叶尺度变换特性的微剂量探测器组织等效换算方法

Country Status (1)

Country Link
CN (1) CN116500665A (zh)

Similar Documents

Publication Publication Date Title
Fuchey et al. Exclusive neutral pion electroproduction in the deeply virtual regime
Muir et al. Monte Carlo calculations of electron beam quality conversion factors for several ion chamber types
Motz et al. K-ionization cross sections for relativistic electrons
Colautti et al. Miniaturized microdosimeters as LET monitors: First comparison of calculated and experimental data performed at the 62 MeV/u 12C beam of INFN-LNS with four different detectors
Suchowerska et al. Perturbation of radiotherapy beams by radiographic film: measurements and Monte Carlo simulations
Zecca et al. An experimental and theoretical investigation into positron and electron scattering from formaldehyde
Casnati et al. Measurements of the total atomic differential cross section of elastic scattering of 59.54-keV photons
CN116500665A (zh) 傅里叶尺度变换特性的微剂量探测器组织等效换算方法
Freier et al. Interaction between D 2 and He 3 in the Neighborhood of the 18.6-Mev Level of Li 5
Abudra’a et al. Dosimetry formalism and calibration procedure for electronic brachytherapy sources in terms of absorbed dose to water
Roberts et al. Kilovoltage energy imaging with a radiotherapy linac with a continuously variable energy range
CN114034721B (zh) 一种微剂量探测器的组织等效修正方法
Brown et al. K-shell ionization cross sections for Ti, Cr, Co, Ni, Cu, and Zn targets induced by protons in the energy range 300—2400 keV
Straticiuc et al. New setup for basic radiobiology studies using a 3 MV TandetronTM: Design and developments
Miklukho et al. Study of nucleon correlations in nuclei by (p, p’) inelastic reaction at 1 GeV
Yu et al. Experimental measurement and Monte Carlo simulation the correction factor for the medium-energy X-ray free-air ionization chamber
CN118447968A (zh) 非组织等效材料的组织等效修正方法
Furre et al. Measurement of dose rate at the interface of cell culture medium and glass dishes by means of ESR dosimetry using thin films of alanine
King et al. Development of Standard X-Ray Beams for Calibration of Radiobiology Cabinet and Conformal Irradiators
Shaw et al. The Clatterbridge high-energy neutron therapy facility: measurements of beam parameters for clinical use
Plagnard et al. X-ray spectrometry for calculating conversion coefficients from air kerma to operational quantities in radiation protection
CN116956600A (zh) 一种高能电子束束流密度校准方法
Rindhatayathon et al. Development of the Primary Standard for IAEA TRS-457 at OAP, Thailand
Liang et al. Establishment of 60-200kV X-ray reference radiation of high air-kerma rate
CN118571384A (zh) 一种基于蒙特卡洛模拟建立高能x射线参考辐射质的方法

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