CN112231899B - X射线连续能谱的确定方法、装置、电子设备及存储介质 - Google Patents
X射线连续能谱的确定方法、装置、电子设备及存储介质 Download PDFInfo
- Publication number
- CN112231899B CN112231899B CN202011069868.9A CN202011069868A CN112231899B CN 112231899 B CN112231899 B CN 112231899B CN 202011069868 A CN202011069868 A CN 202011069868A CN 112231899 B CN112231899 B CN 112231899B
- Authority
- CN
- China
- Prior art keywords
- energy spectrum
- pdd
- continuous
- continuous energy
- matrix data
- 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.)
- Active
Links
- 238000001228 spectrum Methods 0.000 title claims abstract description 276
- 238000000034 method Methods 0.000 title claims abstract description 40
- 239000011159 matrix material Substances 0.000 claims abstract description 83
- 230000006870 function Effects 0.000 claims abstract description 71
- 238000004088 simulation Methods 0.000 claims abstract description 15
- 238000005259 measurement Methods 0.000 claims abstract description 14
- 238000012545 processing Methods 0.000 claims description 9
- 238000004422 calculation algorithm Methods 0.000 claims description 8
- 231100000987 absorbed dose Toxicity 0.000 claims description 6
- 238000004590 computer program Methods 0.000 claims description 5
- 239000012528 membrane Substances 0.000 claims description 3
- 230000003595 spectral effect Effects 0.000 claims 2
- 230000003287 optical effect Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 5
- 238000000342 Monte Carlo simulation Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 150000002500 ions Chemical class 0.000 description 2
- 239000013307 optical fiber Substances 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- WFKWXMTUELFFGS-UHFFFAOYSA-N tungsten Chemical compound [W] WFKWXMTUELFFGS-UHFFFAOYSA-N 0.000 description 2
- 229910052721 tungsten Inorganic materials 0.000 description 2
- 239000010937 tungsten Substances 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 230000005281 excited state Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000008707 rearrangement Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- High Energy & Nuclear Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Measurement Of Radiation (AREA)
Abstract
本发明公开了一种X射线连续能谱的确定方法、装置、电子设备及存储介质,包括:获取X射线能谱数据与能量的对应关系,基于所述对应关系建立初始连续能谱模型;对所述初始连续能谱模型进行离散,获取预设数量的离散能谱函数;基于水下不同深度序列,获取所述预设数量单能的模拟PDD矩阵数据和预设范围连续能谱的测量PDD矩阵数据;基于所述离散能谱函数、所述模拟PDD矩阵数据以及测量PDD矩阵数据确定所述初始连续能谱模型中各能谱系数,得到目标连续能谱模型;获取待估算的连续能量范围,将所述能量范围输入至所述目标连续能谱函数中,得到所述连续能量范围对应的连续能谱。实现了简便快速的获取连续能谱,提高了能谱获取的速率。
Description
技术领域
本发明实施例涉及计算机成像技术领域,尤其涉及一种X射线连续能谱的确定方法、装置、电子设备及存储介质。
背景技术
X射线光子束的能谱信息在CT图像定量分析、能谱CT图像重建及CT辐射剂量计算等领域都发挥着重要的作用。
目前X射线能谱的获取常用间接估计方法,可分为经验型、半经验型和蒙特卡罗模拟型。其中,经验型能谱估计是基于对实际能谱信息进行多项式拟合方法,运算速度较快但缺乏对物理机制的分析;半经验型能谱估计是利用体模装置测量X射线透射数据间接估算,需要对物理过程进行推导,工作量大,时间成本较高;基于蒙特卡罗模拟的能谱模型需要用到在临床实际中并不容易获得的X射线产生装置等内部结构信息,工作难度大。
发明内容
本发明提供X射线连续能谱的确定方法、装置、电子设备及存储介质,实现了简便快速的获取连续能谱,提高了能谱获取的速率。
第一方面,本发明实施例提供了一种X射线连续能谱的确定方法,该方法包括:
获取X射线能谱数据与能量的对应关系,基于所述对应关系建立初始连续能谱模型;
对所述初始连续能谱模型进行离散,获取预设数量的离散能谱函数;
基于水下不同深度序列,获取所述预设数量单能的模拟PDD矩阵数据和预设范围连续能谱的测量PDD矩阵数据;
基于所述离散能谱函数、所述模拟PDD矩阵数据以及测量PDD矩阵数据确定所述初始连续能谱模型中各能谱系数,得到目标连续能谱模型;
获取待估算的连续能量范围,将所述能量范围输入至所述目标连续能谱函数中,得到所述连续能量范围对应的连续能谱。
第二方面,本发明实施例还提供了一种X射线连续能谱的确定装置,该装置包括:
初始连续能谱模型确定模块,用于获取X射线能谱数据与能量的对应关系,基于所述对应关系建立初始连续能谱模型;
离散能谱函数确定模块,用于对所述初始连续能谱模型进行离散,获取预设数量的离散能谱函数;
PDD矩阵数据获取模块,用于基于水下不同深度序列,获取所述预设数量单能的模拟PDD矩阵数据和预设范围连续能谱的测量PDD矩阵数据;
目标连续能谱模型获取模块,用于基于所述离散能谱函数、所述模拟PDD矩阵数据以及测量PDD矩阵数据确定所述初始连续能谱模型中各能谱系数,得到目标连续能谱模型;
连续能谱获取模块,用于获取待估算的连续能量范围,将所述能量范围输入至所述目标连续能谱函数中,得到所述连续能量范围对应的连续能谱。
第三方面,本发明实施例还提供了一种电子设备,所述电子设备包括:
一个或多个处理器;
存储装置,用于存储一个或多个程序,
当所述一个或多个程序被所述一个或多个处理器执行,使得所述一个或多个处理器实现如本发明任意实施例提供的X射线连续能谱的确定方法。
第四方面,本发明实施例还提供了一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现本发明任意实施例提供的X射线连续能谱的确定方法。
本发明通过获取X射线能谱数据与能量的对应关系,基于所述对应关系建立初始连续能谱模型;对所述初始连续能谱模型进行离散,获取预设数量的离散能谱函数;基于水下不同深度序列,获取所述预设数量单能的模拟PDD矩阵数据和预设范围连续能谱的测量PDD矩阵数据;基于所述离散能谱函数、所述模拟PDD矩阵数据以及测量PDD矩阵数据确定所述初始连续能谱模型中各能谱系数,得到目标连续能谱模型;获取待估算的连续能量范围,将所述能量范围输入至所述目标连续能谱函数中,得到所述连续能量范围对应的连续能谱,解决了现有技术中获取连续能谱需要测量大量数据、工作量大的导致效率低下的问题,实现了简便快速的获取连续能谱,提高了能谱获取的速率。
附图说明
图1为本发明实施例一提供的X射线连续能谱的确定方法的流程图;
图2为本发明实施例二提供的X射线连续能谱的确定装置的结构示意图;
图3为本发明实施例三提供的电子设备的结构示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步的详细说明。可以理解的是,此处所描述的具体实施例仅仅用于解释本发明,而非对本发明的限定。另外还需要说明的是,为了便于描述,附图中仅示出了与本发明相关的部分而非全部结构。
实施例一
图1为本发明实施例一中的X射线连续能谱的确定方法的流程图,本实施例可以适用于对连续能谱进行获取的情况,该方法可以由连续能谱获取的装置执行,该连续能谱获取装置可以采用软件和/或硬件的方式实现,本实施例对能谱获取装置的型号和设备不加以限制。如图1所示,该方法具体包括:
S110、获取X射线能谱数据与能量的对应关系,基于所述对应关系建立初始连续能谱模型。
X射线的产生是入射电子使原子内壳层电子激发而发生的现象。即内壳层电子被激发后跳到比费米能高的能级上,电子轨道内出现的空位被外壳层轨道的电子填入时,作为多余能量放出的就是X射线。以高能X射线为试样,将试样原子内层轨道的电子激发出来,形成缺少内层电子的激发态离子。该离子极不稳定,其外层电子迅速向内层空穴跃迁,同时发射另一能量的X射线。此激发过程持续进行可得到不同能量线系的X射线,通过能量色散方式可获得X射线强度随能量变化的能量谱。为了获取X射线连续能谱,可选的,可以利用X光机以及钨能量或者铅能量过滤器产生X射线。具体的,在本发明实施例中,可以对X光机设置不同的电压以便采用不同速度的电子能量轰击钨靶滤波器,即获取产生的不同管电压下的X射线能谱数据,并形成各管电压以及当前管电压下对应的X射线能谱数据的映射关系,并利用非线性拟合软件对各管电压以及当前管电压下对应的X射线能谱数据进行拟合,例如Curve Expert、ORIGIN、MATLAB等非线性曲线拟合软件。
具体的,基于对X射线能谱数据与对应能量进行拟合建立X射线的初始连续能谱模型,得到初始的连续能谱函数表达式为:
其中,E为待估算的连续能谱的电子能量范围,a,b,c和d为所述初始连续能谱模型中待确定的能谱系数,m是X射线特征峰个数,δ()是狄克拉函数,Ei是X射线特征能量。
可选的,该初始连续能谱函数的表达式还可以为:
其中,E为待估算的连续能谱的电子能量范围,a,b,c和d为所述初始连续能谱模型中待确定的能谱系数,m是X射线特征峰个数,δ()是狄克拉函数,Ei是X射线特征能量。
S120、对所述初始连续能谱模型进行离散,获取预设数量的离散能谱函数。
在本发明实施例中,为了求解初始连续能谱模型中的连续能谱函数,可以根据能谱的电子能量范围将连续的能谱函数进行离散,获得以1keV为能量间隔的离散能谱函数。其中,离散能谱的个数与连续能谱的电子能量范围相等。具体的,离散能谱函数的表达式为:
s=[S(E1),S(E2),...,S(EN),]N
其中,s()为离散能谱函数,S()为连续能谱函数,E为待估算的连续能谱的电子能量范围,N为离散能谱的个数。
可选的,将连续能谱函数离散为离散能谱函数,离散能谱函数可以表达为s=F(a,b,c,d),即是函数为s,参数为(a,b,c,d)的函数。通过求解参数(a,b,c,d)进而确定连续能谱函数。
S130、基于水下不同深度序列,获取所述预设数量单能的模拟PDD矩阵数据和预设范围连续能谱的测量PDD矩阵数据。
在本发明实施例中,PDD(percentage depth dose,百分深度剂量)为体膜内射线中心轴上某一深度d处的吸收剂量Dd与参考深度d0处吸收剂量D0之比的百分数,是描述射线中心轴不同深度处相对剂量分布的物理量。具体的,PDD的表达式为:
其中,Dd为射线中心轴上深度为d处的点的吸收剂量,D0为射线中心轴上某一深度d0处的吸收剂量。
可选的,可以采用预设模拟算法获取单能X射线扇形束在特定深度的模拟单能PDD曲线。具体的,预设算法可以是蒙特卡罗模拟算法(EGSnrc/DOSXYZnrc程序包),也可以采用其他模拟算法,本实施例对此不加以限制。通过设置预设模拟算法下模拟射束特征参数,基于能量间隔对预设数量的单能进行模拟,得到预设能量范围内在水下特定深度下的模拟单能谱PDD曲线。可选的,模拟射束参数包括能量范围,能量间隔,射束特征扇形角和能量过滤器厚度。可选的,预设数量为能量范围与能量间隔的比值数据。为了获得较好的模拟结果,可选的,可设置各模拟射束特征参数为:能量范围为1keV~150keV,能量间隔为1keV,射束特征扇形角为55度,能量过滤器厚度为2cm,上述各模拟射束特征参数只是作为可选实施例,也可以根据实际情况进行具体设置,本实施例对此不加以限制。
可选的,基于预设间隔设置水下不同深度序列D(D1,D2,...,DN,),以及预设数量的各个单能X射线水介质中的模拟PDD曲线,计算在水下不同深度序列的各单能模拟PDD曲线,所形成的模拟PDD矩阵数据记为PDDS。具体的,PDDS的表达式为:
其中,n是每个单能PDD曲线的数据个数,N是预设数量单能(E1,E2,...,EN,)的个数。
可选的,可以控制CT设备在水下特定深度序列下进行预设能量范围内连续能谱的PDD数据的测量,得到水下特定深度序列下预设范围内连续能谱PDD曲线。具体的,根据CT设备的型号不同,预设能量可以是80kV、90kV、100kV、120kV和140kV。
可选的,基于预设间隔设置水下不同深度序列获取在水下不同深度序列的各连续能谱的测量PDD曲线,所形成的测量PDD矩阵数据记为PDDm。具体的,PDDm的表达式为:
PDDm(D)=[PDDm(D1),PDDm(D2),...,PDDm(Dn)]n
其中,n是水下深度序列个数。
S140、基于所述离散能谱函数、所述模拟PDD矩阵数据以及测量PDD矩阵数据确定所述初始连续能谱模型中各能谱系数,得到目标连续能谱模型。
在本发明实施例中,为了得到目标连续能谱模型,需要求解初始连续能谱模型中的连续能谱函数中的各能谱系数。可选的,基于连续能谱函数离散得到的离散能谱函数求解各能谱系数。具体的,基于离散能谱函数、模拟PDD矩阵数据以及测量PDD矩阵数据与模型误差之间的预设能谱规则确定模型误差。可选的,预设能谱规则可以表达为:
PDDm=s*PDDs+ε
其中,s为离散能谱函数,PDDs为模拟PDD矩阵数据,PDDm为测量PDD矩阵数据,ε为模型误差。可选的,可以采用跌倒算法,求解模型误差的最小值。可选的,基于所述模型误差的最小值原则,对所述离散能谱函数、所述模拟PDD矩阵数据以及测量PDD矩阵数据进行解析,确定所述连续能谱模型中各能谱系数。具体的,最小原则的表达式为:
(a,b,c,d)=argminF(a,b,c,d)=argmin((PDDm-s*PDDs)2)
可选的,基于上述表达式求解各模型系数,并将各模型系数带入至初始连续能谱模型中确定目标连续能谱模型,得到模型中的目标连续能谱函数。
S150、获取待估算的连续能量范围,将所述能量范围输入至所述目标连续能谱函数中,得到所述连续能量范围对应的连续能谱。
在本发明实施例中,通过求解初始连续能谱模型中的各能谱系数确定目标连续能谱模型。获取待估算的连续能谱的能量范围。具体的,能量范围包括但不限于1keV~80keV、1keV~90keV、1keV~100keV、1keV~120keV、1keV~140keV和1keV~150keV。将待估算的连续能谱的能量范围输入至目标连续能谱模型中,通过计算目标连续能谱函数,得到能量范围内的连续能谱,替代了复杂的连续能谱计算过程,提高了连续能谱的确定效率。
本发明实施例提供了一种X射线连续能谱的确定方法,通过获取X射线能谱数据与能量的对应关系,基于所述对应关系建立初始连续能谱模型;对所述初始连续能谱模型进行离散,获取预设数量的离散能谱函数;基于水下不同深度序列,获取所述预设数量单能的模拟PDD矩阵数据和预设范围连续能谱的测量PDD矩阵数据;基于所述离散能谱函数、所述模拟PDD矩阵数据以及测量PDD矩阵数据确定所述初始连续能谱模型中各能谱系数,得到目标连续能谱模型;获取待估算的连续能量范围,将所述能量范围输入至所述目标连续能谱函数中,得到所述连续能量范围对应的连续能谱,解决现有技术中获取连续能谱工作量大且工作难度大的问题,实现了简便快速的获取连续能谱,提高了能谱获取的速率。
实施例二
图2是本发明实施例二中的X射线连续能谱的确定装置的结构示意图。如图2所示,所述装置包括:
初始连续能谱模型确定模块210,用于获取X射线能谱数据与能量的对应关系,基于所述对应关系建立初始连续能谱模型;
离散能谱函数确定模块220,用于对所述初始连续能谱模型进行离散,获取预设数量的离散能谱函数;
PDD矩阵数据获取模块230,用于基于水下不同深度序列,获取所述预设数量单能的模拟PDD矩阵数据和预设范围连续能谱的测量PDD矩阵数据;
目标连续能谱模型获取模块240,用于基于所述离散能谱函数、所述模拟PDD矩阵数据以及测量PDD矩阵数据确定所述初始连续能谱模型中各能谱系数,得到目标连续能谱模型;
连续能谱获取模块250,用于获取待估算的连续能量范围,将所述能量范围输入至所述目标连续能谱函数中,得到所述连续能量范围对应的连续能谱。
可选的,所述初始连续能谱模型为:
其中,所述E为待估算的连续能量范围,所述a,b,c和d为所述初始连续能谱模型中待确定的能谱系数,m是X射线特征峰个数,δ()是狄克拉函数,Ei是X射线特征能量。
可选的,所述初始连续能谱模型基于对X射线能谱数据与对应能量进行拟合得到。
可选的,所述PDD矩阵数据获取模块230,包括第一PDD矩阵数据获取单元,
所述第一PDD矩阵数据获取单元包括:
第一PDD曲线获取子单元,用于设置预设模拟算法下模拟射束特征参数,基于能量间隔对预设数量的单能进行模拟,得到预设能量范围内在水下不同深度序列下的模拟单能谱PDD曲线。
第一PDD矩阵数据获取单元,用于基于得到的所述模拟单能谱PDD曲线形成模拟PDD矩阵数据。
可选的,所述PDD矩阵数据获取模块230,还包括第二PDD矩阵数据获取单元,
所述第二PDD矩阵数据获取单元包括:
第二PDD曲线获取子单元,用于控制CT设备在水下不同深度序列下,进行预设能量范围内连续能谱的PDD数据的测量,得到水下不同深度序列下,预设范围内连续能谱PDD曲线;
第二PDD矩阵数据获取子单元,用于接收所述CT设备发送的所述水下不同深度序列下预设范围内连续能量的PDD曲线,形成所述测量PDD矩阵数据。
可选的,预设数量的单能是预设能量范围内的离散单能。
可选的,目标连续能谱模型获取模块240,包括:
模型误差确定单元,用于基于所述离散能谱函数、所述模拟PDD矩阵数据以及测量PDD矩阵数据与模型误差之间的预设能谱规则确定模型误差;
能谱系数确定单元,用于基于所述模型误差的最小值原则,对所述离散能谱函数、所述模拟PDD矩阵数据以及测量PDD矩阵数据进行解析,确定所述连续能谱模型中各能谱系数。
本发明实施例提供了一种X射线连续能谱的确定装置,通过获取X射线能谱数据与能量的对应关系,基于所述对应关系建立初始连续能谱模型;对所述初始连续能谱模型进行离散,获取预设数量的离散能谱函数;基于水下不同深度序列,获取所述预设数量单能的模拟PDD矩阵数据和预设范围连续能谱的测量PDD矩阵数据;基于所述离散能谱函数、所述模拟PDD矩阵数据以及测量PDD矩阵数据确定所述初始连续能谱模型中各能谱系数,得到目标连续能谱模型;获取待估算的连续能量范围,将所述能量范围输入至所述目标连续能谱函数中,得到所述连续能量范围对应的连续能谱,解决现有技术中获取连续能谱工作量大且工作难度大的问题,实现了简便快速的获取连续能谱,提高了能谱获取的速率。
实施例三
图3是本发明实施例三中的电子设备的结构示意图。图3示出了适于用来实现本发明实施方式的示例性电子设备312的框图。图3显示的电子设备312仅仅是一个示例,不应对本发明实施例的功能和使用范围带来任何限制。
如图3所示,电子设备312以通用计算设备的形式表现。电子设备312的组件可以包括但不限于:一个或者多个处理器或者处理单元316,系统存储器328,连接不同系统组件(包括系统存储器328和处理单元316)的总线318。
总线318表示几类总线结构中的一种或多种,包括存储器总线或者存储器控制器,外围总线,图形加速端口,处理器或者使用多种总线结构中的任意总线结构的局域总线。举例来说,这些体系结构包括但不限于工业标准体系结构(ISA)总线,微通道体系结构(MAC)总线,增强型ISA总线、视频电子标准协会(VESA)局域总线以及外围组件互连(PCI)总线。
电子设备312典型地包括多种计算机系统可读介质。这些介质可以是任何能够被电子设备312访问的可用介质,包括易失性和非易失性介质,可移动的和不可移动的介质。
系统存储器328可以包括易失性存储器形式的计算机系统可读介质,例如随机存取存储器(RAM)330和/或高速缓存存储器332。电子设备312可以进一步包括其它可移动/不可移动的、易失性/非易失性计算机系统存储介质。仅作为举例,存储系统334可以用于读写不可移动的、非易失性磁介质(图3未显示,通常称为“硬盘驱动器”)。尽管图3中未示出,可以提供用于对可移动非易失性磁盘(例如“软盘”)读写的磁盘驱动器,以及对可移动非易失性光盘(例如CD-ROM,DVD-ROM或者其它光介质)读写的光盘驱动器。在这些情况下,每个驱动器可以通过一个或者多个数据介质接口与总线318相连。存储器328可以包括至少一个程序产品,该程序产品具有一组(例如至少一个)程序模块,这些程序模块被配置以执行本发明各实施例的功能。
具有一组(至少一个)程序模块326的程序/实用工具336,可以存储在例如存储器328中,这样的程序模块326包括但不限于操作系统、一个或者多个应用程序、其它程序模块以及程序数据,这些示例中的每一个或某种组合中可能包括网络环境的实现。程序模块326通常执行本发明所描述的实施例中的功能和/或方法。
电子设备312也可以与一个或多个外部设备314(例如键盘、指向设备、显示器324等)通信,还可与一个或者多个使得用户能与该电子设备312交互的设备通信,和/或与使得该电子设备312能与一个或多个其它计算设备进行通信的任何设备(例如网卡,调制解调器等等)通信。这种通信可以通过输入/输出(I/O)接口322进行。并且,电子设备312还可以通过网络适配器320与一个或者多个网络(例如局域网(LAN),广域网(WAN)和/或公共网络,例如因特网)通信。如图所示,网络适配器320通过总线318与电子设备312的其它模块通信。应当明白,尽管图3中未示出,可以结合电子设备312使用其它硬件和/或软件模块,包括但不限于:微代码、设备驱动器、冗余处理单元、外部磁盘驱动阵列、RAID系统、磁带驱动器以及数据备份存储系统等。
处理单元316通过运行存储在系统存储器328中的程序,从而执行各种功能应用以及数据处理,例如实现本发明实施例所提供的X射线连续能谱的确定方法,该方法包括:
获取X射线能谱数据与能量的对应关系,基于所述对应关系建立初始连续能谱模型;
对所述初始连续能谱模型进行离散,获取预设数量的离散能谱函数;
基于水下不同深度序列,获取所述预设数量单能的模拟PDD矩阵数据和预设范围连续能谱的测量PDD矩阵数据;
基于所述离散能谱函数、所述模拟PDD矩阵数据以及测量PDD矩阵数据确定所述初始连续能谱模型中各能谱系数,得到目标连续能谱模型;
获取待估算的连续能量范围,将所述能量范围输入至所述目标连续能谱函数中,得到所述连续能量范围对应的连续能谱。
实施例四
本发明实施例四还提供了一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现如本发明实施例所提供的X射线连续能谱的确定方法,该方法包括:
获取X射线能谱数据与能量的对应关系,基于所述对应关系建立初始连续能谱模型;
对所述初始连续能谱模型进行离散,获取预设数量的离散能谱函数;
基于水下不同深度序列,获取所述预设数量单能的模拟PDD矩阵数据和预设范围连续能谱的测量PDD矩阵数据;
基于所述离散能谱函数、所述模拟PDD矩阵数据以及测量PDD矩阵数据确定所述初始连续能谱模型中各能谱系数,得到目标连续能谱模型;
获取待估算的连续能量范围,将所述能量范围输入至所述目标连续能谱函数中,得到所述连续能量范围对应的连续能谱。
本发明实施例的计算机存储介质,可以采用一个或多个计算机可读的介质的任意组合。计算机可读介质可以是计算机可读信号介质或者计算机可读存储介质。计算机可读存储介质例如可以是——但不限于——电、磁、光、电磁、红外线、或半导体的系统、装置或器件,或者任意以上的组合。计算机可读存储介质的更具体的例子(非穷举的列表)包括:具有一个或多个导线的电连接、便携式计算机磁盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、光纤、便携式紧凑磁盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。在本文件中,计算机可读存储介质可以是任何包含或存储程序的有形介质,该程序可以被指令执行系统、装置或者器件使用或者与其结合使用。
计算机可读的信号介质可以包括在基带中或者作为载波一部分传播的数据信号,其中承载了计算机可读的程序代码。这种传播的数据信号可以采用多种形式,包括但不限于电磁信号、光信号或上述的任意合适的组合。计算机可读的信号介质还可以是计算机可读存储介质以外的任何计算机可读介质,该计算机可读介质可以发送、传播或者传输用于由指令执行系统、装置或者器件使用或者与其结合使用的程序。
计算机可读介质上包含的程序代码可以用任何适当的介质传输,包括——但不限于无线、电线、光缆、RF等等,或者上述的任意合适的组合。
可以以一种或多种程序设计语言或其组合来编写用于执行本发明操作的计算机程序代码,所述程序设计语言包括面向对象的程序设计语言—诸如Java、Smalltalk、C++,还包括常规的过程式程序设计语言—诸如”C”语言或类似的程序设计语言。程序代码可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或电子设备上执行。在涉及远程计算机的情形中,远程计算机可以通过任意种类的网络——包括局域网(LAN)或广域网(WAN)—连接到用户计算机,或者,可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。
注意,上述仅为本发明的较佳实施例及所运用技术原理。本领域技术人员会理解,本发明不限于这里所述的特定实施例,对本领域技术人员来说能够进行各种明显的变化、重新调整和替代而不会脱离本发明的保护范围。因此,虽然通过以上实施例对本发明进行了较为详细的说明,但是本发明不仅仅限于以上实施例,在不脱离本发明构思的情况下,还可以包括更多其他等效实施例,而本发明的范围由所附的权利要求范围决定。
Claims (10)
1.一种X射线连续能谱的确定方法,其特征在于,包括:
获取X射线能谱数据与能量的对应关系,基于所述对应关系建立初始连续能谱模型;
对所述初始连续能谱模型进行离散,获取预设数量的离散能谱函数;
基于水下不同深度序列,获取所述预设数量单能的模拟PDD矩阵数据和预设范围连续能谱的测量PDD矩阵数据;
基于所述离散能谱函数、所述模拟PDD矩阵数据以及测量PDD矩阵数据确定所述初始连续能谱模型中各能谱系数,得到目标连续能谱模型;
获取待估算的连续能量范围,将所述能量范围输入至目标连续能谱函数中,得到所述连续能量范围对应的连续能谱;
其中,所述PDD为体膜内射线中心轴上深度d处的吸收剂量Dd与参考深度d0处吸收剂量D0之比的百分数,是描述射线中心轴不同深度处相对剂量分布的物理量。
3.根据权利要求1所述的方法,其特征在于,所述初始连续能谱模型基于对X射线能谱数据与对应能量进行拟合得到。
4.根据权利要求1所述的方法,其特征在于,所述基于水下不同深度序列,获取所述预设数量单能的模拟PDD矩阵数据,包括:
设置预设模拟算法下模拟射束特征参数,基于能量间隔对预设数量的单能进行模拟,得到预设能量范围内在水下不同深度序列下的模拟单能谱PDD曲线;
基于得到的所述模拟单能谱PDD曲线形成模拟PDD矩阵数据。
5.根据权利要求4所述的方法,其特征在于,所述基于水下不同深度序列,获取预设范围连续能谱的测量PDD矩阵数据,包括:
控制CT设备在水下不同深度序列下,进行预设能量范围内连续能谱的PDD数据的测量,得到水下不同深度序列下,预设范围内连续能谱PDD曲线;
接收所述CT设备发送的所述水下不同深度序列下预设范围内连续能谱的PDD曲线,形成所述测量PDD矩阵数据。
6.根据权利要求5所述的方法,其特征在于,预设数量的单能是预设能量范围内的离散单能。
7.根据权利要求1所述的方法,所述基于所述离散能谱函数、所述模拟PDD矩阵数据以及测量PDD矩阵数据确定所述初始连续能谱模型中各能谱参数,得到目标连续能谱模型,包括:
基于所述离散能谱函数、所述模拟PDD矩阵数据以及测量PDD矩阵数据与模型误差之间的预设能谱规则确定模型误差;
基于所述模型误差的最小值原则,对所述离散能谱函数、所述模拟PDD矩阵数据以及测量PDD矩阵数据进行解析,确定所述连续能谱模型中各能谱系数。
8.一种X射线连续能谱的确定装置,其特征在于,包括:
初始连续能谱模型确定模块,用于获取X射线能谱数据与能量的对应关系,基于所述对应关系建立初始连续能谱模型;
离散能谱函数确定模块,用于对所述初始连续能谱模型进行离散,获取预设数量的离散能谱函数;
PDD矩阵数据获取模块,用于基于水下不同深度序列,获取所述预设数量单能的模拟PDD矩阵数据和预设范围连续能谱的测量PDD矩阵数据;
目标连续能谱模型获取模块,用于基于所述离散能谱函数、所述模拟PDD矩阵数据以及测量PDD矩阵数据确定所述初始连续能谱模型中各能谱系数,得到目标连续能谱模型;
连续能谱获取模块,用于获取待估算的连续能量范围,将所述能量范围输入至目标连续能谱函数中,得到所述连续能量范围对应的连续能谱;
其中,所述PDD为体膜内射线中心轴上深度d处的吸收剂量Dd与参考深度d0处吸收剂量D0之比的百分数,是描述射线中心轴不同深度处相对剂量分布的物理量。
9.一种电子设备,其特征在于,所述电子设备包括:
一个或多个处理器;
存储装置,用于存储一个或多个程序,
当所述一个或多个程序被所述一个或多个处理器执行,使得所述一个或多个处理器实现如权利要求1-7中任一所述的X射线连续能谱的确定方法。
10.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,该程序被处理器执行时实现如权利要求1-7中任一所述的X射线连续能谱的确定方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011069868.9A CN112231899B (zh) | 2020-09-30 | 2020-09-30 | X射线连续能谱的确定方法、装置、电子设备及存储介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011069868.9A CN112231899B (zh) | 2020-09-30 | 2020-09-30 | X射线连续能谱的确定方法、装置、电子设备及存储介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112231899A CN112231899A (zh) | 2021-01-15 |
CN112231899B true CN112231899B (zh) | 2023-02-03 |
Family
ID=74119511
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011069868.9A Active CN112231899B (zh) | 2020-09-30 | 2020-09-30 | X射线连续能谱的确定方法、装置、电子设备及存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112231899B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116304519B (zh) * | 2023-05-12 | 2023-07-28 | 苏州益腾电子科技有限公司 | X射线管的实时热量计算方法、装置、系统和存储介质 |
CN118131302B (zh) * | 2024-05-07 | 2024-08-30 | 迈胜医疗设备有限公司 | 医用粒子加速器初级能谱的计算方法及相关设备 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101477205A (zh) * | 2009-01-22 | 2009-07-08 | 中国科学技术大学 | 基于多算法的放射源反演方法 |
CN103135125A (zh) * | 2011-12-01 | 2013-06-05 | 中国辐射防护研究院 | 一种过滤x射线参考辐射的能谱分析方法 |
CN103336295A (zh) * | 2013-06-19 | 2013-10-02 | 南京航空航天大学 | 一种医用电子直线加速器光子束能谱获取方法 |
CN105866821A (zh) * | 2016-03-25 | 2016-08-17 | 上海联影医疗科技有限公司 | 一种获取射束能谱的方法、装置及剂量分布计算方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102015023B (zh) * | 2008-05-13 | 2014-01-22 | 三菱电机株式会社 | 粒子射线治疗装置及粒子射线治疗方法 |
-
2020
- 2020-09-30 CN CN202011069868.9A patent/CN112231899B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101477205A (zh) * | 2009-01-22 | 2009-07-08 | 中国科学技术大学 | 基于多算法的放射源反演方法 |
CN103135125A (zh) * | 2011-12-01 | 2013-06-05 | 中国辐射防护研究院 | 一种过滤x射线参考辐射的能谱分析方法 |
CN103336295A (zh) * | 2013-06-19 | 2013-10-02 | 南京航空航天大学 | 一种医用电子直线加速器光子束能谱获取方法 |
CN105866821A (zh) * | 2016-03-25 | 2016-08-17 | 上海联影医疗科技有限公司 | 一种获取射束能谱的方法、装置及剂量分布计算方法 |
Non-Patent Citations (2)
Title |
---|
Photon Energy Spectrum Reconstruction Based on Monte Carlo and Measured Percentage Depth Dose in Accurate Radiotherapy;Gui LI et al.;《Progress in NUCLEAR SCIENCE and TECHNOLOGY》;20111231;第2卷;第160-164页 * |
基于模拟退火法的医用电子加速器6MV X射线能谱重建;刘娟等;《中国生物医学工程学报》;20130831;第32卷(第04期);第385-389页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112231899A (zh) | 2021-01-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112231899B (zh) | X射线连续能谱的确定方法、装置、电子设备及存储介质 | |
Bourque et al. | A stoichiometric calibration method for dual energy computed tomography | |
Bert et al. | Geant4-based Monte Carlo simulations on GPU for medical applications | |
US20160095561A1 (en) | Photon counting apparatus | |
Collins-Fekete et al. | A maximum likelihood method for high resolution proton radiography/proton CT | |
Sterpin et al. | Analytical computation of prompt gamma ray emission and detection for proton range verification | |
Marcotulli et al. | NuSTAR perspective on high-redshift MeV blazars | |
Omar et al. | A model for the energy and angular distribution of x rays emitted from an x‐ray tube. Part I. Bremsstrahlung production | |
López-Corredoira et al. | Cosmological test with the QSO Hubble diagram | |
O'Connell et al. | fastCAT: Fast cone beam CT (CBCT) simulation | |
Traoré Dubuis et al. | Total cross section of furfural by electron impact: Experiment and theory | |
CN111001097B (zh) | 一种放疗剂量评估系统、装置及存储介质 | |
Shunhavanich et al. | Fluid‐filled dynamic bowtie filter: Description and comparison with other modulators | |
Krah et al. | Scattering proton CT | |
Bontempi et al. | Monte Carlo simulator of realistic x‐ray beam for diagnostic applications | |
Cajgfinger et al. | Fixed forced detection for fast SPECT Monte-Carlo simulation | |
Timoshenko et al. | Simulation of radiation field inside interplanetary spacecraft | |
Somasundaram et al. | Development and validation of an open source Monte Carlo dosimetry model for wide‐beam CT scanners using Fluka | |
Lakshmanan et al. | Simulations of nuclear resonance fluorescence in Geant4 | |
Badali et al. | Characterization of an x-ray source with a partitioned diamond-tungsten target for electronic brachytherapy with 3D beam directionality | |
JP7092921B2 (ja) | 制御装置 | |
Scot et al. | The Monte Carlo code MCSHAPE: Main features and recent developments | |
RU2527137C2 (ru) | Способ определения изотопного отношения делящегося вещества, содержащегося в камере деления | |
Kirkpatrick et al. | Comparing the output of measured and GEANT4 simulated X-ray tubes | |
US11872071B2 (en) | Method for correcting a spectral image |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
OL01 | Intention to license declared | ||
OL01 | Intention to license declared |