CN111714147A - 能量刻度曲线获取方法、装置、计算机设备和存储介质 - Google Patents
能量刻度曲线获取方法、装置、计算机设备和存储介质 Download PDFInfo
- Publication number
- CN111714147A CN111714147A CN202010604518.1A CN202010604518A CN111714147A CN 111714147 A CN111714147 A CN 111714147A CN 202010604518 A CN202010604518 A CN 202010604518A CN 111714147 A CN111714147 A CN 111714147A
- Authority
- CN
- China
- Prior art keywords
- energy
- decay characteristic
- gamma decay
- adc
- values
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 56
- 230000005253 gamme decay Effects 0.000 claims abstract description 75
- 239000013078 crystal Substances 0.000 claims abstract description 46
- 238000011088 calibration curve Methods 0.000 claims abstract description 45
- 230000005855 radiation Effects 0.000 claims abstract description 20
- 239000002245 particle Substances 0.000 claims abstract description 19
- 238000012216 screening Methods 0.000 claims abstract description 16
- 238000004590 computer program Methods 0.000 claims description 10
- 238000006243 chemical reaction Methods 0.000 claims description 8
- 238000012937 correction Methods 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 4
- 230000002285 radioactive effect Effects 0.000 abstract description 5
- 238000001514 detection method Methods 0.000 abstract description 4
- 238000002600 positron emission tomography Methods 0.000 description 30
- 230000008569 process Effects 0.000 description 10
- 238000004891 communication Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 8
- BPQQTUXANYXVAA-UHFFFAOYSA-N Orthosilicate Chemical compound [O-][Si]([O-])([O-])[O-] BPQQTUXANYXVAA-UHFFFAOYSA-N 0.000 description 4
- ANDNPYOOQLLLIU-UHFFFAOYSA-N [Y].[Lu] Chemical compound [Y].[Lu] ANDNPYOOQLLLIU-UHFFFAOYSA-N 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000001228 spectrum Methods 0.000 description 4
- 239000000126 substance Substances 0.000 description 4
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 3
- 230000004060 metabolic process Effects 0.000 description 3
- 229910052710 silicon Inorganic materials 0.000 description 3
- 239000010703 silicon Substances 0.000 description 3
- WQZGKKKJIJFFOK-GASJEMHNSA-N Glucose Natural products OC[C@H]1OC(O)[C@H](O)[C@@H](O)[C@@H]1O WQZGKKKJIJFFOK-GASJEMHNSA-N 0.000 description 2
- 230000002776 aggregation Effects 0.000 description 2
- 238000004220 aggregation Methods 0.000 description 2
- 238000002591 computed tomography Methods 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 235000014113 dietary fatty acids Nutrition 0.000 description 2
- 229930195729 fatty acid Natural products 0.000 description 2
- 239000000194 fatty acid Substances 0.000 description 2
- 150000004665 fatty acids Chemical class 0.000 description 2
- 239000008103 glucose Substances 0.000 description 2
- 238000009206 nuclear medicine Methods 0.000 description 2
- 102000039446 nucleic acids Human genes 0.000 description 2
- 108020004707 nucleic acids Proteins 0.000 description 2
- 150000007523 nucleic acids Chemical class 0.000 description 2
- 102000004169 proteins and genes Human genes 0.000 description 2
- 108090000623 proteins and genes Proteins 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 229910052797 bismuth Inorganic materials 0.000 description 1
- JCXGWMGPZLAOME-UHFFFAOYSA-N bismuth atom Chemical compound [Bi] JCXGWMGPZLAOME-UHFFFAOYSA-N 0.000 description 1
- OKTJSMMVPCPJKN-BJUDXGSMSA-N carbon-11 Chemical compound [11C] OKTJSMMVPCPJKN-BJUDXGSMSA-N 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 229940079593 drug Drugs 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000005251 gamma ray Effects 0.000 description 1
- 238000007429 general method Methods 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000004020 luminiscence type Methods 0.000 description 1
- BOYZAERJCXIRAX-UHFFFAOYSA-N lutetium(3+);trisilicate Chemical compound [Lu+3].[Lu+3].[Lu+3].[Lu+3].[O-][Si]([O-])([O-])[O-].[O-][Si]([O-])([O-])[O-].[O-][Si]([O-])([O-])[O-] BOYZAERJCXIRAX-UHFFFAOYSA-N 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000002503 metabolic effect Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000007493 shaping process Methods 0.000 description 1
- FVAUCKIRQBBSSJ-UHFFFAOYSA-M sodium iodide Chemical compound [Na+].[I-] FVAUCKIRQBBSSJ-UHFFFAOYSA-M 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/58—Testing, adjusting or calibrating thereof
- A61B6/582—Calibration
- A61B6/585—Calibration of detector units
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/037—Emission tomography
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/42—Arrangements for detecting radiation specially adapted for radiation diagnosis
- A61B6/4208—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector
- A61B6/4241—Arrangements for detecting radiation specially adapted for radiation diagnosis characterised by using a particular type of detector using energy resolving detectors, e.g. photon counting
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Medical Informatics (AREA)
- Engineering & Computer Science (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Physics & Mathematics (AREA)
- High Energy & Nuclear Physics (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Measurement Of Radiation (AREA)
Abstract
本申请涉及一种能量刻度曲线获取方法、装置、计算机设备和存储介质,其中,该能量刻度曲线获取方法包括:获取探测器接收到的来自于探测器晶体自身辐射粒子的本底事件;对所述本底事件进行甄别,得到多个ADC值以及对应的核素的伽玛衰变特征能峰值;根据多个所述ADC值以及对应的核素的伽玛衰变特征能峰值,得到能量刻度曲线。通过探测器晶体的本底固有放射现象,避免了使用放射源,减少了医师受到辐射的剂量,减少了医院的使用成本,并且相对于有源检测,其操作简单,可以在任何时期获取能量刻度曲线。
Description
技术领域
本申请涉及医疗技术领域,特别是涉及一种能量刻度曲线获取方法、装置、计算机设备和存储介质。
背景技术
正电子发射计算机断层扫描(PET,Positron Emission Computed Tomography),是核医学领域比较先进的临床检查影像技术。其大致方法是,将某种物质,一般是生物生命代谢中必须的物质,如:葡萄糖、蛋白质、核酸、脂肪酸,标记上短寿命的放射性核素(如F18,碳11等),注入人体后,通过对于该物质在代谢中的聚集,来反映生命代谢活动的情况,从而达到诊断的目的。由于正电子发射计算机断层扫描具有灵敏度高、特异性高以及安全性好等特点,因此,正电子发射计算机断层扫描被更为广泛的应用。本申请中的ADC是指模数转换器,是指将连续变化的模拟信号转换成离散的数字信号的器件;ADC值是指模数转换器的采样值。PET探测器在接收到光子之后,通过康普顿散射或光电作用沉积能量到相应的探测器晶体中,并激发出与所沉积的能量成比例的可见光光子,这些光子通过光电倍增管(PMT)或硅光电倍增器(SiPM)转换为相应的电信号,再通过后续电路进行增益、成形、数字化后形成ADC值,通过已知的能量刻度曲线可算出已知ADC值所表达的光子能量。
目前的传统技术中,使用Ge-68以及Na-22等正电子放射桶源或棒源进行数据采集,并得到相应的能量刻度曲线。但由于存在放射源,因此会对医师造成一定程度的放射伤害。增加了医师受到辐射的剂量,并且操作流程复杂。
发明内容
本申请实施例提供了一种能量刻度曲线获取方法、装置、计算机设备和存储介质,以至少解决相关技术中医师受辐射剂量高以及操作流程复杂的问题。
第一方面,本申请实施例提供了一种能量刻度曲线获取方法,包括:获取探测器接收到的来自于探测器晶体自身辐射粒子的本底事件;对所述本底事件进行甄别,得到多个ADC值以及对应的核素的伽玛衰变特征能峰值;根据多个所述ADC值以及对应的核素的伽玛衰变特征能峰值,得到能量刻度曲线。
在其中一个实施例中,所述对所述本底事件进行甄别,得到多个ADC值以及对应的核素的伽玛衰变特征能峰值包括:根据所述本底事件,得到多种核素的伽玛衰变特征能峰;对多种所述核素的伽玛衰变特征能峰进行模数转换,得到相应特征能峰的ADC值;选取至少三个伽玛衰变特征能峰以及相应伽玛衰变特征能峰对应的ADC值。
在其中一个实施例中,所述选取至少三个伽玛衰变特征能峰以及相应伽玛衰变特征能峰对应的ADC值包括:选取0keV、307keV和596keV三个伽玛衰变特征能峰以及相应伽玛衰变特征能峰对应的ADC值。
在其中一个实施例中,所述根据多个所述ADC值以及对应的核素的伽玛衰变特征能峰值,得到能量刻度曲线包括:对多个所述ADC值以及对应的核素的伽玛衰变特征能峰值进行插值拟合,得到能量刻度曲线。
在其中一个实施例中,所述根据多个所述ADC值以及对应的核素的伽玛衰变特征能峰值,得到能量刻度曲线之后包括:根据所述能量刻度曲线,确定能量刻度曲线上任意能峰值对应的ADC值。
在其中一个实施例中,所述根据多个所述ADC值以及对应的核素的伽玛衰变特征能峰值,得到能量刻度曲线之后包括:根据所述能量刻度曲线,确定511keV能峰值对应的ADC值。
在其中一个实施例中,所述探测器接收来自于探测器晶体自身辐射粒子的本底事件的时长不小于30分钟。
第二方面,本申请实施例提供了一种能量刻度曲线获取装置,包括:获取模块,用于获取探测器接收到的来自于探测器晶体自身辐射粒子的本底事件;甄别模块,用于对所述本底事件进行甄别,得到多个ADC值以及对应的核素的伽玛衰变特征能峰值;曲线计算模块,用于根据多个所述ADC值以及对应的核素的伽玛衰变特征能峰值,得到能量刻度曲线。
第三方面,本申请实施例提供了一种PET探测器晶体能量校正方法,利用上述任一种能量刻度曲线获取方法获取能量刻度曲线,确定511keV能峰值对应的ADC值,根据所述确定后的511keV能峰值对应的ADC值对PET探测器晶体能量进行校正。
第四方面,本申请实施例提供了一种计算机设备,包括存储器、处理器以及存储在所述存储器上并可在所述处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现如上述第一方面所述的能量刻度曲线获取方法或第三方面所述的PET探测器晶体能量校正方法。
第五方面,本申请实施例提供了一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时实现如上述第一方面所述的能量刻度曲线获取方法或第三方面所述的PET探测器晶体能量校正方法。
相比于相关技术,本申请实施例提供的能量刻度曲线获取方法,通过获取探测器接收到的来自于探测器自身辐射粒子的本底事件;再对接收到的本底事件进行甄别,得到多个ADC值以及对应的核素的伽玛衰变特征能峰值;最后根据多个ADC值以及与对应核素的伽玛衰变特征能峰值,得到能量刻度曲线。通过探测器晶体的本底固有放射现象,避免了使用放射源,减少了医师受到辐射的剂量,减少了医院的使用成本,并且相对于有源检测,其操作简单,可以在任何时期获取能量刻度曲线。
本申请的一个或多个实施例的细节在以下附图和描述中提出,以使本申请的其他特征、目的和优点更加简明易懂。
附图说明
此处所说明的附图用来提供对本申请的进一步理解,构成本申请的一部分,本申请的示意性实施例及其说明用于解释本申请,并不构成对本申请的不当限定。在附图中:
图1为一个实施例中能量刻度曲线获取方法的流程图;
图2为一个实施例中能量刻度曲线获取装置的结构框图;
图3为一个实施例中计算机设备的硬件结构示意图。
具体实施方式
为了使本申请的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本申请进行描述和说明。应当理解,此处所描述的具体实施例仅仅用以解释本申请,并不用于限定本申请。基于本申请提供的实施例,本领域普通技术人员在没有作出创造性劳动的前提下所获得的所有其他实施例,都属于本申请保护的范围。
显而易见地,下面描述中的附图仅仅是本申请的一些示例或实施例,对于本领域的普通技术人员而言,在不付出创造性劳动的前提下,还可以根据这些附图将本申请应用于其他类似情景。此外,还可以理解的是,虽然这种开发过程中所作出的努力可能是复杂并且冗长的,然而对于与本申请公开的内容相关的本领域的普通技术人员而言,在本申请揭露的技术内容的基础上进行的一些设计,制造或者生产等变更只是常规的技术手段,不应当理解为本申请公开的内容不充分。
在本申请中提及“实施例”意味着,结合实施例描述的特定特征、结构或特性可以包含在本申请的至少一个实施例中。在说明书中的各个位置出现该短语并不一定均是指相同的实施例,也不是与其它实施例互斥的独立的或备选的实施例。本领域普通技术人员显式地和隐式地理解的是,本申请所描述的实施例在不冲突的情况下,可以与其它实施例相结合。
除非另作定义,本申请所涉及的技术术语或者科学术语应当为本申请所属技术领域内具有一般技能的人士所理解的通常意义。本申请所涉及的“一”、“一个”、“一种”、“该”等类似词语并不表示数量限制,可表示单数或复数。本申请所涉及的术语“包括”、“包含”、“具有”以及它们任何变形,意图在于覆盖不排他的包含;例如包含了一系列步骤或模块(单元)的过程、方法、***、产品或设备没有限定于已列出的步骤或单元,而是可以还包括没有列出的步骤或单元,或可以还包括对于这些过程、方法、产品或设备固有的其它步骤或单元。本申请所涉及的“连接”、“相连”、“耦接”等类似的词语并非限定于物理的或者机械的连接,而是可以包括电气的连接,不管是直接的还是间接的。本申请所涉及的“多个”是指两个或两个以上。“和/或”描述关联对象的关联关系,表示可以存在三种关系,例如,“A和/或B”可以表示:单独存在A,同时存在A和B,单独存在B这三种情况。字符“/”一般表示前后关联对象是一种“或”的关系。本申请所涉及的术语“第一”、“第二”、“第三”等仅仅是区别类似的对象,不代表针对对象的特定排序。
正电子发射型计算机断层显像(Positron Emission Computed Tomography,PET),是核医学领域比较先进的临床检查影像技术。是将某种物质,一般是生物生命代谢中必须的物质,如:葡萄糖、蛋白质、核酸、脂肪酸,标记上短寿命的放射性核素(如18F,11C等),注入人体后,放射性核素在衰变过程中释放出正电子,一个正电子在行进十分之几毫米到几毫米后遇到一个电子后发生湮灭,从而产生方向相反的一对能量为511KeV的光子。这对光子,通过高度灵敏的照相机捕捉,并经计算机进行散射和随机信息的校正。经过对不同的正电子进行相同的分析处理,我们可以得到在生物体内聚集情况的三维图像,从而达到诊断的目的。
目前常用于PET设备的探测器主要有碘化钠(NaI)、锗酸铋(BGO)、硅酸镥(LSO)、硅酸钇镥(LYSO)等闪烁晶体探测器。本实施例中以硅酸钇镥(LYSO)闪烁晶体探测器为例进行说明。硅酸钇镥闪烁晶体(LYSO)以其高光输出、快发光衰减、有效原子序数多、密度大、物化性质稳定以及对γ射线探测效率高等特点被用作PET设备的闪烁晶体。由于LYSO中存在的Lu176,因此具有固有放射现象,Lu176在衰变过程中释放出能量分别为88kev、202kev和307kev的三种能量的γ光子以及能量随机范围为0kev~593kev的β粒子。
能量刻度,是将PET***采集到的指代粒子能量的ADC值映射为实际能量的转换操作。PET探测器在接收到光子之后,激发出可见光光子,这些光子通过光电倍增管(PMT)或硅光电倍增器(SiPM)转换为相应的电信号,该电信号经过模数转换生成ADC值,通过对已知能量与相应ADC值进行映射,完成能量刻度,得到能量刻度曲线。其中,ADC值即为采集到的电信号进行模数转换后得到的数值。
本实施例提供的能量刻度曲线获取方法,通过探测器晶体的固有放射现象,采集伽玛特征能峰的本底事件,可在晶体所得到的本底事件能谱中得到307kev特征能峰和596kev特征能峰,由已知的三个特征能峰对实际应用中的511kev特征能峰进行插值拟合,得到ADC值与实际光子能量的对应关系。
本实施例提供了一种能量刻度曲线获取方法。图1为一个实施例中能量刻度曲线获取方法的流程图,如图1所示,该流程包括如下步骤:
步骤S102,获取探测器接收到的来自于探测器晶体自身辐射粒子的本底事件。
具体地,探测器晶体为硅酸钇镥闪烁晶体(LYSO),由于LYSO中存在的Lu176,在Lu176衰变时,产生β粒子以及γ光子。其中β粒子在发生衰变的探测器晶体被立即吸收,而γ光子需要通过整个视野到达对面的晶体才能被相应的探测器晶体吸收,接收β粒子和γ光子的探测器晶体的连线称为响应线,探测器接收到β粒子和γ光子的事件称为本底事件。
在其中一个实施例中,所述探测器接收来自于探测器晶体自身辐射粒子的本底事件的时长不小于30分钟。具体地,PET设备在空闲时间,采集无源状态下的本底事件,由于本底中各类散射事件的比例较高,很容易降低本底能谱伽玛特征能峰的信噪比,因此,对于本底事件理论上采集的时间越长越好。由于,实际需要,在无源情况下至少采集30分钟的本底事件,才能保证每个探测器晶体显示出更好的能谱结构。
在其中一个实施例中,可以预先设置PET***的时间窗和能量窗;基于所述时间窗和能量窗获取本底事件。其中,所述时间窗的值不小于所述PET***临床阈值;所述能量窗的值不小于所述PET***临床阈值。时间窗和能量窗用于体现PET设备的灵敏度。其中,时间窗的选定主要依据时间分辨率进行设定,时间窗的选择应比时间分辨率稍大;时间分辨率的定义为对已知符合事件中相对的两个探测器相应的时间差分布的半高宽。能量窗的选定主要依据能量分辨率进行设定,好的能量分辨率可以选择较小的能量窗。在本实施例中,为了更好的适应LYSO本底固有的放射性现象,需要设置合适的时间窗以及能量窗。优选的,时间窗的选定应该不小于PET***临床时间窗阈值,能量窗的选定应该不小于PET***临床能量窗阈值。其中时间窗阈值以及能量窗阈值,均为在PET***实际应用过程中,通过大量数据统计得到的时间窗阈值以及能量窗阈值。
步骤S104,对所述本底事件进行甄别,得到多个ADC值以及对应的核素的伽玛衰变特征能峰值。
具体地,根据所述本底事件,得到多种核素的伽玛衰变特征能峰;对多种所述核素的伽玛衰变特征能峰进行模数转换,得到相应特征能峰的ADC值;选取至少三个伽玛衰变特征能峰以及相应伽玛衰变特征能峰对应的ADC值。更具体的,探测器晶体在接收到β粒子或γ光子之后,沉积能量到相应的探测器晶体中,并激发出与所沉积的能量成比例的可见光光子,这些光子通过光电倍增管(PMT)或硅光电倍增器(SiPM)转换为相应的电信号,根据电信号生成本底事件的能谱图。对该能谱图进行能量甄别,得到多种核素的伽玛衰变特征能峰。其中,该核素包括Lu176等探测器晶体中存在的具有放射性的核素。以核素为Lu176为例进行举例说明,Lu176在衰变过程中释放出能量分别为88kev、202kev和307kev的三种能量的γ光子以及能量随机范围为0kev~593kev的β粒子,因此,Lu176对应的伽玛衰变特征能峰包括:88kev光子的伽玛衰变特征能峰、202kev光子的伽玛衰变特征能峰、307kev光子的伽玛衰变特征能峰、596kev光子的伽玛衰变特征能峰以及各种能量值对应的伽玛衰变特征能峰。对获取到的伽玛衰变特征能峰进行模数转换,得到相应伽玛衰变特征能峰的ADC值。在得到所有伽玛衰变特征能峰以及相应伽玛衰变特征能峰的ADC值之后,从中选取部分伽玛衰变特征能峰以及与其对应的ADC值。由于在能量刻度的过程中,越多的伽玛衰变特征能峰和与其对应ADC值被标定出来,能量刻度就越准确。在本实施例中,优选的,选取0keV、307keV和596keV三个伽玛衰变特征能峰以及相应伽玛衰变特征能峰对应的ADC值。
步骤S106,根据多个所述ADC值以及对应的核素的伽玛衰变特征能峰值,得到能量刻度曲线。
具体地,对多个所述ADC值以及对应的核素的伽玛衰变特征能峰值进行插值拟合,得到能量刻度曲线。更具体的,根据选取的0keV、307keV和596keV三个伽玛衰变特征能峰以及相应伽玛衰变特征能峰对应的ADC值,进行插值拟合,得到能量刻度曲线。该能量刻度曲线能够覆盖所有能量范围,其能量范围为0到2000keV。
上述能量刻度曲线获取方法,通过获取探测器接收到的来自于探测器自身辐射粒子的本底事件;再对接收到的本底事件进行甄别,得到多个ADC值以及对应的核素的伽玛衰变特征能峰值;最后根据多个ADC值以及与对应核素的伽玛衰变特征能峰值,得到能量刻度曲线。通过探测器晶体的本底固有放射现象,避免了使用放射源,减少了医师受到辐射的剂量,减少了医院的使用成本,并且相对于有源检测,其操作简单,可以在任何时期获取能量刻度曲线。
在其中一个实施例中,根据多个所述ADC值以及对应的核素的伽玛衰变特征能峰值,得到能量刻度曲线之后,根据所述能量刻度曲线,确定能量刻度曲线上任意能峰值对应的ADC值。具体的,根据能量刻度曲线,得到所有能峰值以及对应的ADC值。在实际的使用过程中,每一种特征峰核素的药物在PET实际使用中都有应用存在,因此在实际使用过程中可能会碰到任一种光子能量,需要提前获取所有能峰值以及对应的ADC值,以便于在后期实际PET扫描过程中,找到实际的光子能量。
在其中一个实施例中,根据多个所述ADC值以及对应的核素的伽玛衰变特征能峰值,得到能量刻度曲线之后,根据所述能量刻度曲线,确定511keV能峰值对应的ADC值。511keV能峰是目前比较常用的PET放射性核素注射剂所产生的特征峰,因此,需要提前根据能量刻度曲线,得到511keV能峰值对应的ADC值。
本实施例提供的能量刻度曲线获取方法,利用探测器晶体中存在的Lu-176核素自然衰变所产生的307keV和596keV,通过长期累积探测器晶体自然衰变产生的本地数据,对两个特征能峰进行拟合甄别,最终依据0keV、307keV和596keV三个特征峰以及对应的ADC值,得到能量刻度曲线。再根据能量刻度曲线得到实际使用中关注的511keV对应的ADC值。上述能量刻度曲线获取方法避免了对放射源的使用,提高了操作的安全性和可行性,有效的利用无患者扫描情况下PET设备的本底数据,更有效及时的对探测器的性能刻度进行更新。
在其中一个实施例中,还提供了一种PET探测器晶体能量校正方法,利用上述任一个实施例提供的能量刻度曲线获取方法获取能量刻度曲线,确定511keV能峰值对应的ADC值,根据所述确定后的511keV能峰值对应的ADC值对PET探测器晶体能量进行校正。具体的,在实际的使用过程中,首先根据获取到的能量刻度曲线,找到能量刻度曲线中511keV能峰值对应的ADC值。根据511keV能峰值对应的ADC值,对探测器晶体接收到的能量进行校正,使探测器对接收到的能量判别更加的准确,进一步的提高PET设备的成像精准度。
需要说明的是,在上述流程中或者附图的流程图中示出的步骤可以在诸如一组计算机可执行指令的计算机***中执行,并且,虽然在流程图中示出了逻辑顺序,但是在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤。
本实施例还提供了一种能量刻度曲线获取装置,该装置用于实现上述实施例及优选实施方式,已经进行过说明的不再赘述。如以下所使用的,术语“模块”、“单元”、“子单元”等可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。
图2为一个实施例中能量刻度曲线获取装置的结构框图,如图2所示,该装置包括:获取模块100、甄别模块200以及曲线计算模块300。
获取模块100,用于获取探测器接收到的来自于探测器晶体自身辐射粒子的本底事件;
甄别模块200,用于对所述本底事件进行甄别,得到多个ADC值以及对应的核素的伽玛衰变特征能峰值;
曲线计算模块300,用于根据多个所述ADC值以及对应的核素的伽玛衰变特征能峰值,得到能量刻度曲线。
甄别模块200,还用于根据所述本底事件,得到多种核素的伽玛衰变特征能峰;对多种所述核素的伽玛衰变特征能峰进行模数转换,得到相应特征能峰的ADC值;选取至少三个伽玛衰变特征能峰以及相应伽玛衰变特征能峰对应的ADC值。
甄别模块200,还用于选取0keV、307keV和596keV三个伽玛衰变特征能峰以及相应伽玛衰变特征能峰对应的ADC值。
曲线计算模块300,还用于对多个所述ADC值以及对应的核素的伽玛衰变特征能峰值进行插值拟合,得到能量刻度曲线。
能量刻度曲线获取装置还包括:提取模块。
提取模块,用于根据所述能量刻度曲线,确定能量刻度曲线上任意能峰值对应的ADC值。
提取模块,还用于根据所述能量刻度曲线,确定511keV能峰值对应的ADC值。
需要说明的是,上述各个模块可以是功能模块也可以是程序模块,既可以通过软件来实现,也可以通过硬件来实现。对于通过硬件来实现的模块而言,上述各个模块可以位于同一处理器中;或者上述各个模块还可以按照任意组合的形式分别位于不同的处理器中。
另外,结合图1描述的本申请实施例能量刻度曲线获取方法以及PET探测器晶体能量校正方法可以由计算机设备来实现。图3为根据本申请实施例的计算机设备的硬件结构示意图。
计算机设备可以包括处理器31以及存储有计算机程序指令的存储器32。
具体地,上述处理器31可以包括中央处理器(CPU),或者特定集成电路(Application Specific Integrated Circuit,简称为ASIC),或者可以被配置成实施本申请实施例的一个或多个集成电路。
其中,存储器32可以包括用于数据或指令的大容量存储器。举例来说而非限制,存储器32可包括硬盘驱动器(Hard Disk Drive,简称为HDD)、软盘驱动器、固态驱动器(SolidState Drive,简称为SSD)、闪存、光盘、磁光盘、磁带或通用串行总线(Universal SerialBus,简称为USB)驱动器或者两个或更多个以上这些的组合。在合适的情况下,存储器32可包括可移除或不可移除(或固定)的介质。在合适的情况下,存储器32可在数据处理装置的内部或外部。在特定实施例中,存储器32是非易失性(Non-Volatile)存储器。在特定实施例中,存储器32包括只读存储器(Read-Only Memory,简称为ROM)和随机存取存储器(RandomAccess Memory,简称为RAM)。在合适的情况下,该ROM可以是掩模编程的ROM、可编程ROM(Programmable Read-Only Memory,简称为PROM)、可擦除PROM(Erasable ProgrammableRead-Only Memory,简称为EPROM)、电可擦除PROM(Electrically Erasable ProgrammableRead-Only Memory,简称为EEPROM)、电可改写ROM(Electrically Alterable Read-OnlyMemory,简称为EAROM)或闪存(FLASH)或者两个或更多个以上这些的组合。在合适的情况下,该RAM可以是静态随机存取存储器(Static Random-Access Memory,简称为SRAM)或动态随机存取存储器(Dynamic Random Access Memory,简称为DRAM),其中,DRAM可以是快速页模式动态随机存取存储器(Fast Page Mode Dynamic Random Access Memory,简称为FPMDRAM)、扩展数据输出动态随机存取存储器(Extended Date Out Dynamic RandomAccess Memory,简称为EDODRAM)、同步动态随机存取内存(Synchronous Dynamic Random-Access Memory,简称SDRAM)等。
存储器32可以用来存储或者缓存需要处理和/或通信使用的各种数据文件,以及处理器31所执行的可能的计算机程序指令。
处理器31通过读取并执行存储器32中存储的计算机程序指令,以实现上述实施例中的任意一种能量刻度曲线获取方法。
在其中一些实施例中,计算机设备还可包括通信接口33和总线30。其中,如图3所示,处理器31、存储器32、通信接口33通过总线30连接并完成相互间的通信。
通信接口33用于实现本申请实施例中各模块、装置、单元和/或设备之间的通信。通信端口33还可以实现与其他部件例如:外接设备、图像/数据采集设备、数据库、外部存储以及图像/数据处理工作站等之间进行数据通信。
总线30包括硬件、软件或两者,将计算机设备的部件彼此耦接在一起。总线30包括但不限于以下至少之一:数据总线(Data Bus)、地址总线(Address Bus)、控制总线(Control Bus)、扩展总线(Expansion Bus)、局部总线(Local Bus)。举例来说而非限制,总线30可包括图形加速接口(Accelerated Graphics Port,简称为AGP)或其他图形总线、增强工业标准架构(Extended Industry Standard Architecture,简称为EISA)总线、前端总线(Front Side Bus,简称为FSB)、超传输(Hyper Transport,简称为HT)互连、工业标准架构(Industry Standard Architecture,简称为ISA)总线、无线带宽(InfiniBand)互连、低引脚数(Low Pin Count,简称为LPC)总线、存储器总线、微信道架构(Micro ChannelArchitecture,简称为MCA)总线、***组件互连(Peripheral Component Interconnect,简称为PCI)总线、PCI-Express(PCI-X)总线、串行高级技术附件(Serial AdvancedTechnology Attachment,简称为SATA)总线、视频电子标准协会局部(Video ElectronicsStandards Association Local Bus,简称为VLB)总线或其他合适的总线或者两个或更多个以上这些的组合。在合适的情况下,总线80可包括一个或多个总线。尽管本申请实施例描述和示出了特定的总线,但本申请考虑任何合适的总线或互连。
该计算机设备可以基于获取到的计算机指令,执行本申请实施例中的能量刻度曲线获取方法,从而实现结合图1描述的能量刻度曲线获取方法以及PET探测器晶体能量校正方法。
另外,结合上述实施例中的能量刻度曲线获取方法以及PET探测器晶体能量校正方法,本申请实施例可提供一种计算机可读存储介质来实现。该计算机可读存储介质上存储有计算机程序指令;该计算机程序指令被处理器执行时实现上述实施例中的任意一种能量刻度曲线获取方法以及PET探测器晶体能量校正方法。
以上所述实施例的各技术特征可以进行任意的组合,为使描述简洁,未对上述实施例中的各个技术特征所有可能的组合都进行描述,然而,只要这些技术特征的组合不存在矛盾,都应当认为是本说明书记载的范围。
以上所述实施例仅表达了本申请的几种实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本申请构思的前提下,还可以做出若干变形和改进,这些都属于本申请的保护范围。因此,本申请专利的保护范围应以所附权利要求为准。
Claims (10)
1.一种能量刻度曲线获取方法,其特征在于,包括:
获取探测器接收到的来自于探测器晶体自身辐射粒子的本底事件;
对所述本底事件进行甄别,得到多个ADC值以及对应的核素的伽玛衰变特征能峰值;
根据多个所述ADC值以及对应的核素的伽玛衰变特征能峰值,得到能量刻度曲线。
2.根据权利要求1所述的能量刻度曲线获取方法,其特征在于,所述对所述本底事件进行甄别,得到多个ADC值以及对应的核素的伽玛衰变特征能峰值包括:
根据所述本底事件,得到多种核素的伽玛衰变特征能峰;
对多种所述核素的伽玛衰变特征能峰进行模数转换,得到相应特征能峰的ADC值;
选取至少三个伽玛衰变特征能峰以及相应伽玛衰变特征能峰对应的ADC值。
3.根据权利要求2所述的能量刻度曲线获取方法,其特征在于,所述选取至少三个伽玛衰变特征能峰以及相应伽玛衰变特征能峰对应的ADC值包括:
选取0keV、307keV和596keV三个伽玛衰变特征能峰以及相应伽玛衰变特征能峰对应的ADC值。
4.根据权利要求1所述的能量刻度曲线获取方法,其特征在于,所述根据多个所述ADC值以及对应的核素的伽玛衰变特征能峰值,得到能量刻度曲线包括:
对多个所述ADC值以及对应的核素的伽玛衰变特征能峰值进行插值拟合,得到能量刻度曲线。
5.根据权利要求1所述的能量刻度曲线获取方法,其特征在于,所述根据多个所述ADC值以及对应的核素的伽玛衰变特征能峰值,得到能量刻度曲线之后包括:
根据所述能量刻度曲线,确定能量刻度曲线上任意能峰值对应的ADC值。
6.根据权利要求1所述的能量刻度曲线获取方法,其特征在于,所述根据多个所述ADC值以及对应的核素的伽玛衰变特征能峰值,得到能量刻度曲线之后包括:
根据所述能量刻度曲线,确定511keV能峰值对应的ADC值。
7.根据权利要求1所述的能量刻度曲线获取方法,其特征在于,
所述探测器接收来自于探测器晶体自身辐射粒子的本底事件的时长不小于30分钟。
8.一种能量刻度曲线获取装置,其特征在于,包括:
获取模块,用于获取探测器接收到的来自于探测器晶体自身辐射粒子的本底事件;
甄别模块,用于对所述本底事件进行甄别,得到多个ADC值以及对应的核素的伽玛衰变特征能峰值;
曲线计算模块,用于根据多个所述ADC值以及对应的核素的伽玛衰变特征能峰值,得到能量刻度曲线。
9.一种PET探测器晶体能量校正方法,其特征在于,利用权利要求1至7中任一种能量刻度曲线获取方法获取能量刻度曲线,确定511keV能峰值对应的ADC值,根据所述确定后的511keV能峰值对应的ADC值对PET探测器晶体能量进行校正。
10.一种计算机设备,包括存储器、处理器以及存储在所述存储器上并可在所述处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现如权利要求1至7中任一项所述的能量刻度曲线获取方法或如权利要求9所述的PET探测器晶体能量校正方法。
Priority Applications (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010604518.1A CN111714147A (zh) | 2020-06-29 | 2020-06-29 | 能量刻度曲线获取方法、装置、计算机设备和存储介质 |
PCT/CN2021/095840 WO2021238929A1 (zh) | 2020-05-25 | 2021-05-25 | 一种成像设备的状态检测方法和*** |
EP21814520.9A EP4144296A4 (en) | 2020-05-25 | 2021-05-25 | CONDITION DETECTION METHOD AND SYSTEM FOR AN IMAGING DEVICE |
JP2022572797A JP7512440B2 (ja) | 2020-05-25 | 2021-05-25 | 画像形成装置の状態の検出方法及びシステム |
US18/058,785 US20230102139A1 (en) | 2020-05-25 | 2022-11-25 | Systems and methods for state detection of imaging devices |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010604518.1A CN111714147A (zh) | 2020-06-29 | 2020-06-29 | 能量刻度曲线获取方法、装置、计算机设备和存储介质 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111714147A true CN111714147A (zh) | 2020-09-29 |
Family
ID=72569733
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010604518.1A Pending CN111714147A (zh) | 2020-05-25 | 2020-06-29 | 能量刻度曲线获取方法、装置、计算机设备和存储介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111714147A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2021238929A1 (zh) * | 2020-05-25 | 2021-12-02 | 上海联影医疗科技股份有限公司 | 一种成像设备的状态检测方法和*** |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109259786A (zh) * | 2018-09-19 | 2019-01-25 | 明峰医疗***股份有限公司 | 基于lyso闪烁体pet***的能量自刻度方法 |
CN109709596A (zh) * | 2018-11-23 | 2019-05-03 | 南京航空航天大学 | 一种放射性能谱无源自动能量刻度方法 |
-
2020
- 2020-06-29 CN CN202010604518.1A patent/CN111714147A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109259786A (zh) * | 2018-09-19 | 2019-01-25 | 明峰医疗***股份有限公司 | 基于lyso闪烁体pet***的能量自刻度方法 |
CN109709596A (zh) * | 2018-11-23 | 2019-05-03 | 南京航空航天大学 | 一种放射性能谱无源自动能量刻度方法 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2021238929A1 (zh) * | 2020-05-25 | 2021-12-02 | 上海联影医疗科技股份有限公司 | 一种成像设备的状态检测方法和*** |
JP7512440B2 (ja) | 2020-05-25 | 2024-07-08 | 上海聯影医療科技股▲ふん▼有限公司 | 画像形成装置の状態の検出方法及びシステム |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Surti et al. | Advances in time-of-flight PET | |
Ter-Pogossian et al. | Super PETT I: a positron emission tomograph utilizing photon time-of-flight information | |
Gu et al. | NEMA NU-4 performance evaluation of PETbox4, a high sensitivity dedicated PET preclinical tomograph | |
Gu et al. | Performance evaluation of HiPET, a high sensitivity and high resolution preclinical PET tomograph | |
Del Guerra et al. | Performance evaluation of the fully engineered YAP-(S) PET scanner for small animal imaging | |
EP1688092A2 (en) | Radiological inspection apparatus and radiological inspection method | |
US9606199B2 (en) | MR-compatible blood sampling system for PET imaging applications in combined PET/MR imaging system | |
CN110772274A (zh) | 一种多核素的伽马成像***与方法 | |
CN108474862B (zh) | 具有lu谱减除的能量校准 | |
CN111568452A (zh) | Pet***状态检测方法、装置、计算机设备和存储介质 | |
CN111685785B (zh) | Pet晶***置查找表的校正方法、装置以及计算机设备 | |
CN111568453A (zh) | 能量校正状态检测方法、装置、计算机设备和存储介质 | |
King et al. | Introduction to the physics of molecular imaging with radioactive tracers in small animals | |
Livieratos | Basic principles of SPECT and PET imaging | |
CN113057667B (zh) | Pet探测器信号采样方法、装置、电子装置和存储介质 | |
KR101174485B1 (ko) | Pet 장치 및 pet 장치의 신호 처리 방법 | |
KR101595929B1 (ko) | 양전자방출단층촬영장치용 반응 위치 판별 방법 및 시스템 | |
Dahlbom | PET Imaging: Basic and New Trends | |
CN111714147A (zh) | 能量刻度曲线获取方法、装置、计算机设备和存储介质 | |
CN112932517B (zh) | 本底放射性实现时间自校正方法及装置 | |
Surti et al. | PET instrumentation | |
CN111728625B (zh) | Pet***时间同步方法、装置、计算机设备和存储介质 | |
Kijewski | Positron emission tomography (PET) and single-photon emission computed tomography (SPECT) physics | |
CN113712584A (zh) | Pet扫描方法、装置、电子装置和存储介质 | |
Raylman et al. | Corrections for the effects of accidental coincidences, Compton scatter, and object size in positron emission mammography (PEM) imaging |
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 | ||
CB02 | Change of applicant information | ||
CB02 | Change of applicant information |
Address after: 201807 Shanghai City, north of the city of Jiading District Road No. 2258 Applicant after: Shanghai Lianying Medical Technology Co., Ltd Address before: 201807 Shanghai City, north of the city of Jiading District Road No. 2258 Applicant before: SHANGHAI UNITED IMAGING HEALTHCARE Co.,Ltd. |