CN111736207A - 一种pet时间标定方法 - Google Patents

一种pet时间标定方法 Download PDF

Info

Publication number
CN111736207A
CN111736207A CN202010469508.1A CN202010469508A CN111736207A CN 111736207 A CN111736207 A CN 111736207A CN 202010469508 A CN202010469508 A CN 202010469508A CN 111736207 A CN111736207 A CN 111736207A
Authority
CN
China
Prior art keywords
time
channel
threshold
value
pet
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
Application number
CN202010469508.1A
Other languages
English (en)
Other versions
CN111736207B (zh
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.)
Guangdong Mingfeng Medical Technology Co ltd
Original Assignee
Guangdong Mingfeng Medical Technology Co ltd
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 Guangdong Mingfeng Medical Technology Co ltd filed Critical Guangdong Mingfeng Medical Technology Co ltd
Priority to CN202010469508.1A priority Critical patent/CN111736207B/zh
Publication of CN111736207A publication Critical patent/CN111736207A/zh
Application granted granted Critical
Publication of CN111736207B publication Critical patent/CN111736207B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01TMEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T7/00Details of radiation-measuring instruments
    • G01T7/005Details of radiation-measuring instruments calibration techniques
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01TMEASUREMENT OF NUCLEAR OR X-RADIATION
    • G01T1/00Measuring X-radiation, gamma radiation, corpuscular radiation, or cosmic radiation
    • G01T1/29Measurement performed on radiation beams, e.g. position or section of the beam; Measurement of spatial distribution of radiation
    • G01T1/2914Measurement of spatial distribution of radiation
    • G01T1/2985In depth localisation, e.g. using positron emitters; Tomographic imaging (longitudinal and transverse section imaging; apparatus for radiation diagnosis sequentially in different planes, steroscopic radiation diagnosis)

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Molecular Biology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Nuclear Medicine (AREA)
  • Measurement Of Radiation (AREA)

Abstract

本发明公开了一种PET时间标定方法,涉及医学影像设备技术领域,基于动态阈值的双阈值比较器实现时间信号的甄别,包括一个高阈值比较器和一个动态低阈值比较器,高阈值比较器对每个探测通道标定能量参考电压,低阈值比较器对每个探测通道标定时间参考电压。采用高阈值标定各探测通道的能量参考电压,并获取各通道增益比例关系,基于外部源统计各通道的时间柱状图,采用步进迭代的方式获取各通道的时间分辨率,并将各通道时间分辨率最优的阈值作为该通道的最佳低阈值,最终实现对各探测通道时间的精准标定,从而达到了提升TOF‑PET系统时间分辨率及图像质量的目的。

Description

一种PET时间标定方法
技术领域
本发明涉及医学影像设备技术领域,尤其是一种PET时间标定方法。
背景技术
正电子发射计算机断层显像(Positron Emission Computed Tomography,简称PET)是目前最先进的医疗诊断设备之一。其工作原理是将同位素标记药物注入体内,该放射性核素会释放出正电子,并与体内的负电子快速发生湮灭辐射,并产生两个能量相同,方向相反的γ光子。
PET系统通过采用环绕人体的核探测器装置,可以获取光子被探测器探测到的时间、位置和能量信息,并根据该信息重建图像。常用的核探测器包括由多个闪烁晶体组成的晶体阵列(Crystal array)和光电探测器。其中,晶体阵列用于接收病患体内释放出的光子(例如γ光子)并将其转换成可见光,光电探测器用于将可见光转换成脉冲信号,脉冲信号可以为电压信号或电流信号,其对应的脉冲波形在时间轴上呈现线性上升和指数衰减的特性,因此,通过设置阈值比较器,将脉冲波形上超过阈值的时刻确定为脉冲信号的时间信息,即:脉冲信号的过阈时刻。
PET成像的关键在于γ射线的脉冲信号采集,如何精确获取光子的能量和时间信息,是决定PET成像质量好坏的关键。尤其对于飞行时间(TOF,TimeofFlight)PET的信号采集,精确的时间信息识别是至关重要,可以大幅度降低噪声信号的干扰,为图像重建提供更加精确的信息。
但是现有技术中,通过单一固定阈值鉴别出脉冲信号的时间信息(时刻)与γ光子击中闪烁体的实际时刻有较大差异,这就会导致最终图像信噪比降低。因此,需要一种更为精准的时间标定方法,实现对TOF-PET系统各通道的时间标定,降低原始信号的噪声干扰,最终实现TOF-PET系统时间分辨率的提升。
例如,公布号为CN107843914A的发明专利申请,公布了一种PET时间标定方法和PET系统,通过根据各个探测通道的通道基线值,确定时间标定参考电压和能量标定参考电压,进而确定单位电压延时,可以根据不同通道之间的基线偏差,结合单位电压延时来校正单事件时间,从而使得单事件时间的测定更加准确。
又如,公布号为CN109893154A的发明专利申请,公布了一种基于低秩约束的PET时间校正方法,通过在PET时间校正线性模型的基础上加入低秩约束,提高了对大噪声的鲁棒性,同时再解决问题时采用稀疏最小二乘的方式,可以解决因为大量数据所带来的计算时间问题。同时,本发明使PET系统能够获得更好的时间分辨率,从而利用TOF信息获取更好的空间分辨率,使PET系统能够在医学成像中提供更为准确的信息,为临床诊断提供更好的帮助。
本发明旨在提供一种新的PET时间标定方法,基于动态阈值的双阈值比较器实现时间信号的甄别,采用高阈值标定各探测通道的能量参考电压,并获取各通道增益比例关系,基于外部源统计各通道的时间柱状图,采用步进迭代的方式获取各通道的时间分辨率,并将各通道时间分辨率最优的阈值作为该通道的最佳低阈值,最终实现对各探测通道时间的精准标定,从而达到了提升TOF-PET系统时间分辨率及图像质量的目的。
发明内容
本发明要解决的技术问题是:解决通过单一固定阈值鉴别出脉冲信号的时间信息与γ光子击中闪烁体的实际时刻有较大差异的问题。
本发明为解决上述技术问题所采用的技术方案是:提供一种PET时间标定方法,采用动态阈值的双阈值比较器实现时间信息的甄别,包括以下步骤:
S1.获取TOF-PET系统各探测通道增益值;
S2.选取block内任意一个通道的增益值作为参考基准,获取block内各通道的增益值比例关系;
S3.增益值比例关系作为对应block内各通道的动态阈值比例关系,更快的找到合适的阈值;
S4.根据动态阈值比例关系,获取各通道的动态阈值的初始阈值;
S5.根据初始阈值,获取各通道的时间分布柱状图;
S6.根据时间分布柱状图,对各通道进行时间标定;
S7.根据时间标定值获取对应通道时间标定后的时间分布柱状图,并获取相应的时间分辨率;
S8.根据动态阈值比例关系及参考基准,更新初始阈值;重复S5、S6、S7步骤,直到达到预设条件;
S9.根据时间分辨率,获取时间分辨率与动态阈值的关系曲线;
S10.根据关系曲线,获取各通道对应的最佳动态阈值;
S11.根据最佳动态阈值,对对应通道的时间进行最终的标定。
优选的,上述的双阈值比较器包括一个高阈值比较器,用于对每个探测通道标定能量参考电压,高阈值为固定阈值或动态阈值。
优选的,上述的双阈值比较器还包括一个低阈值比较器,用于对每个探测通道标定时间参考电压,低阈值为动态阈值。
优选的,上述的增益值的获取包括以下步骤:
(1)将外部源置于视野中心,采集符合计数;
(2)对所采集到的符合计数进行数据分析,获取各通道的能量分布曲线;
(3)对所述能量分布曲线进行寻峰,获取能量峰值的横坐标值Epeak
(4)通过步骤(3)中获取的Epeak计算各通道的增益值Gain,Gaini,j=511/Epeak,i,j,其中,i为block内通道数,j为PET探测器系统的block数。
优选的,上述的增益值比例关系为各通道的增益值与参考基准的比例,其中,以block内任意一个通道作为参考基准通道,参考基准通道的增益值为参考基准。
优选的,上述的时间标定值
Figure BDA0002513818720000031
其中,Thi,j为初始阈值,meani,j为对应通道的时间分布曲线的峰值的横坐标值,即时间值,Th为经验值,Th=16。
优选的,上述的预设条件为Thi,j=Th。
优选的,上述的时间分辨率T_Resi,j=2.354×sigma;其中sigma为对应时间分布曲线的均方根值。
优选的,上述的更新初始阈值为使用步进迭代的方式更新,即,将T_Vali,j作为对应通道的初始阈值,重复S5、S6、S7步骤,直到达到预设条件。
本发明的有益效果是:基于动态阈值的双阈值比较器实现时间信号的甄别,采用高阈值标定各探测通道的能量参考电压,并获取各通道增益比例关系,基于外部源统计各通道的时间柱状图,采用步进迭代的方式获取各通道的时间分辨率,并将各通道时间分辨率最优的阈值作为该通道的最佳低阈值,最终实现对各探测通道时间的精准标定,从而达到了提升TOF-PET系统时间分辨率及图像质量的目的。不做时间标定的时间分辨率为1097ps;传统时间标定的时间分辨率为512ps;本发明的其中一个实施例的时间分辨率为444ps;较传统的时间标定方法提高了68ps,较传统时间标定方式下的时间分辨率指标提升了13%,提升效果显著,达到了本发明技术方案的目的。
附图说明
下面结合附图和实施例对本发明进一步说明。
图1为PET时间标定流程图;
图2为PET探测器系统模块分布示意图;
图3为模块内block分布示意图;
图4为block内通道数示意图;
图5为能量分布曲线;
图6为能量分布曲线峰值示意图;
图7为时间分布柱状图示意图;
图8为时间标定后的时间柱状图示意图;
图9为动态阈值与时间分辨率关系曲线示意图;
图10为不做时间标定的时间分布曲线;
图11为传统时间校准方式的时间分布曲线;
图12为本发明实施例1的时间分布曲线。
具体实施方式
【实施例1】
为了使本发明的目的、技术方案和优点更加的清晰,结合附图和实施例,对本发明做进一步阐述。
本实施例中的TOF-PET系统时间标定方法是基于动态阈值的双阈值比较器的方法实现。
PET系统共有38个模块,如图2所示;每个模块有24个block,如图3所示,共计912个block;每个block有8×8个SiPM、8×8个晶体以及8×8个电子学通道组成,如图4所示,共计58368个通道。
定义i为block内通道数,j为PET探测器系统的block数。
具体实施步骤如下,如图1所示:
(1)获取TOF-PET系统各探测通道增益值,包括以下步骤:
将外部源置于视野中心,采集符合计数;
对所采集到的符合计数进行数据分析,获取各通道的能量分布曲线;如图5所示;
分别对所有通道的能量分布曲线进行寻峰,获取对应通道的能量峰值的横坐标Epeak,如图6所示;
以block为最小单元,获取各block内所有通道的增益值Gain,即:Gaini,j=511/Epeak,i,j
(2)在block内选取任意一个通道作为参考基准通道,获取block内各通道的增益值比例关系;本实施例中选取通道0为参考基准通道,根据block内所有通道的增益值,获得Gain0,j:Gain1,j:Gain2,j:…Gain63,j=1:1.23:1.06:…:1.36。
(3)获取block内各通道的动态阈值Thi,j比例关系,Th0,j:Th1,j:Th2,j:…Th63,j=Gain0,j:Gain1,j:Gain2,j:…Gain63,j=1:1.23:1.06:…:1.36。
(4)获取各block通道的动态初始阈值,若block内第0通道的动态初始阈值为Th0,j=5,根据上述动态阈值比例关系,则其它通道动态初始阈值分别为:Th1,j=6.15,Th2,j=5.3,…,Th63,j=6.8。
(5)自动加载各通道动态初始阈值,并将外部射线源置于视野中心,采集符合数据。
(6)数据分析,获取各通道的时间分布柱状图,并获取对应通道时间分布曲线的峰值的横坐标meani,j值;如图7所示。
(7)根据所得到的meani,j值,对对应通道做时间标定,时间标定值
Figure BDA0002513818720000051
其中,Th为经验值,本实施例中Th=16。
(8)根据标定值,获取对应通道时间标定后的时间分布柱状图,并获取相应的时间分辨率T_Resi,j=2.354×sigma;其中sigma为对应时间分布曲线的均方根值,如图8所示。
(9)将对应通道的时间标定值作为该通道的动态阈值,并更新该通道的动态初始阈值。
(10)重新加载更新后的动态初始阈值,新动态初始阈值为T_Vali,j,并重复上述(5)-(9)步骤。
(11)当Thi,j=Th时,该通道时间标定结束。
(12)获取各通道的动态阈值与时间分辨率的关系曲线,如图9所示。
(13)获取时间分辨率最优时所对应的动态阈值作为该通道的最终动态阈值,最终完成所有通道的时间标定。
本实施方案的时间标定效果验证,如下:
未做时间标定,TOF-PET系统的时间分布曲线及时间分辨率如图10所示,不做时间标定的时间分辨率为1097ps。
传统时间标定,TOF-PET系统的时间分布曲线及时间分辨率如图11所示,传统时间标定的时间分辨率为512ps。
基于本发明实施例1的TOF-PET系统的时间分布曲线及时间分辨率如图12所示,本发明实施例1的时间分辨率为444ps;较传统的时间标定方法提高了68ps,较传统时间标定方式下的时间分辨率指标提升了13%。提升效果显著,达到了本发明技术方案的目的。

Claims (10)

1.一种PET时间标定方法,其特征在于,采用动态阈值的双阈值比较器实现时间信息的甄别,包括以下步骤:
S1.获取TOF-PET系统各探测通道增益值;
S2.选取block内任意一个通道的增益值作为参考基准,根据所述增益值,获取block内各通道的增益值比例关系;
S3.所述增益值比例关系作为对应block内各通道的动态阈值比例关系;
S4.根据所述动态阈值比例关系,获取各通道的动态阈值的初始阈值;
S5.根据所述初始阈值,获取各通道的时间分布柱状图;
S6.根据所述时间分布柱状图,对各通道进行时间标定;
S7.根据时间标定值获取对应通道时间标定后的时间分布柱状图,并获取相应的时间分辨率;
S8.根据所述动态阈值比例关系及参考基准,更新所述初始阈值;重复S5、S6、S7步骤,直到达到预设条件;
S9.根据所述时间分辨率,获取时间分辨率与动态阈值的关系曲线;
S10.根据所述关系曲线,获取各通道对应的最佳动态阈值;
S11.根据所述最佳动态阈值,对对应通道的时间进行最终的标定。
2.根据权利要求1所述的一种PET时间标定方法,其特征在于:所述双阈值比较器包括一个高阈值比较器和一个低阈值比较器。
3.根据权利要求2所述的一种PET时间标定方法,其特征在于:所述高阈值比较器用于对每个探测通道标定能量参考电压,高阈值为固定阈值或动态阈值。
4.根据权利要求2所述的一种PET时间标定方法,其特征在于:所述低阈值比较器用于对每个探测通道标定时间参考电压,低阈值为动态阈值。
5.根据权利要求1所述的一种PET时间标定方法,其特征在于:所述增益值的获取包括以下步骤:
S1.将外部源置于视野中心,采集符合计数;
S2.对所采集到的符合计数进行数据分析,获取各通道的能量分布曲线;
S3.对所述能量分布曲线进行寻峰,获取能量峰值的横坐标值Epeak
S4.通过步骤S3中获取的Epeak计算各通道的增益值Gain,Gaini,j=511/Epeak,i,j,其中,i为block内通道数,j为PET探测器系统的block数。
6.根据权利要求1所述的一种PET时间标定方法,其特征在于:所述增益值比例关系为各通道的增益值与所述参考基准的比例,其中,以block内任意一个通道作为参考基准通道,参考基准通道的增益值为参考基准。
7.根据权利要求5所述的一种PET时间标定方法,其特征在于:所述时间标定值
Figure FDA0002513818710000021
其中,Thi,j为初始阈值,meani,j为对应通道的时间分布曲线的峰值的横坐标值,即时间值,Th为经验值,Th=16。
8.根据权利要求7所述的一种PET时间标定方法,其特征在于:所述预设条件为Thi,j=Th。
9.根据权利要求7所述的一种PET时间标定方法,其特征在于:所述时间分辨率T_Resi,j=2.354×sigma;其中sigma为对应时间分布曲线的均方根值。
10.根据权利要求7所述的一种PET时间标定方法,其特征在于:S8步骤的所述更新初始阈值为使用步进迭代的方式更新,所述步进迭代的方式为将T_Vali,j作为对应通道的更新后的初始阈值,重复S5、S6、S7步骤,直到达到预设条件。
CN202010469508.1A 2020-05-28 2020-05-28 一种pet时间标定方法 Active CN111736207B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010469508.1A CN111736207B (zh) 2020-05-28 2020-05-28 一种pet时间标定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010469508.1A CN111736207B (zh) 2020-05-28 2020-05-28 一种pet时间标定方法

Publications (2)

Publication Number Publication Date
CN111736207A true CN111736207A (zh) 2020-10-02
CN111736207B CN111736207B (zh) 2022-06-07

Family

ID=72646462

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010469508.1A Active CN111736207B (zh) 2020-05-28 2020-05-28 一种pet时间标定方法

Country Status (1)

Country Link
CN (1) CN111736207B (zh)

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080265167A1 (en) * 2006-01-03 2008-10-30 Koninklijke Philips Electronics N. V. Achieving accurate time-of-flight calibrations with a stationary coincidence point source
US20160209528A1 (en) * 2013-08-30 2016-07-21 Uniwersytet Jagiellonski A method for calibration of tof-pet detectors using cosmic radiation
US20160216386A1 (en) * 2013-08-30 2016-07-28 Uniwersytet Jagiellonski A method for determining parameters of a reaction of a gamma quantum within a scintillator of a pet scanner
CN107843914A (zh) * 2017-10-09 2018-03-27 沈阳东软医疗系统有限公司 一种pet时间标定方法和pet系统
US20190012812A1 (en) * 2017-07-10 2019-01-10 Shenyang Neusoft Medical Systems Co., Ltd. Calibrating Time In PET Device
CN109893154A (zh) * 2019-03-01 2019-06-18 浙江大学 一种基于低秩约束的pet时间校正方法
CN110074806A (zh) * 2019-05-29 2019-08-02 明峰医疗系统股份有限公司 一种基于SiPM探测系统的增益控制方法及装置
CN110179485A (zh) * 2019-05-29 2019-08-30 明峰医疗系统股份有限公司 一种用于pet成像系统增益调整的方法和装置
CN110432922A (zh) * 2019-08-07 2019-11-12 广东明峰医疗科技有限公司 一种提高pet系统时间校准精度的方法
CN110988834A (zh) * 2019-11-22 2020-04-10 航天恒星科技有限公司 一种基于自适应阈值双门限的脉冲到达时间测量方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080265167A1 (en) * 2006-01-03 2008-10-30 Koninklijke Philips Electronics N. V. Achieving accurate time-of-flight calibrations with a stationary coincidence point source
US20160209528A1 (en) * 2013-08-30 2016-07-21 Uniwersytet Jagiellonski A method for calibration of tof-pet detectors using cosmic radiation
US20160216386A1 (en) * 2013-08-30 2016-07-28 Uniwersytet Jagiellonski A method for determining parameters of a reaction of a gamma quantum within a scintillator of a pet scanner
US20190012812A1 (en) * 2017-07-10 2019-01-10 Shenyang Neusoft Medical Systems Co., Ltd. Calibrating Time In PET Device
CN107843914A (zh) * 2017-10-09 2018-03-27 沈阳东软医疗系统有限公司 一种pet时间标定方法和pet系统
CN109893154A (zh) * 2019-03-01 2019-06-18 浙江大学 一种基于低秩约束的pet时间校正方法
CN110074806A (zh) * 2019-05-29 2019-08-02 明峰医疗系统股份有限公司 一种基于SiPM探测系统的增益控制方法及装置
CN110179485A (zh) * 2019-05-29 2019-08-30 明峰医疗系统股份有限公司 一种用于pet成像系统增益调整的方法和装置
CN110432922A (zh) * 2019-08-07 2019-11-12 广东明峰医疗科技有限公司 一种提高pet系统时间校准精度的方法
CN110988834A (zh) * 2019-11-22 2020-04-10 航天恒星科技有限公司 一种基于自适应阈值双门限的脉冲到达时间测量方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
赵然然: "乳腺专用PET 系统PEMi的定量标定方法研究", 《中国医学影像学杂志》 *

Also Published As

Publication number Publication date
CN111736207B (zh) 2022-06-07

Similar Documents

Publication Publication Date Title
US11191510B2 (en) Imaging system and method based on multiple-gamma photon coincidence event
AU2009230876B2 (en) Radiation imaging method with individual signal resolution
US6297506B1 (en) System and method for reducing pile-up errors in multi-crystal gamma ray detector applications
EP2994776B1 (en) Apparatus and method for the evaluation of gamma radiation events
US8716669B2 (en) Line of response estimation for high-resolution PET detector
CN106716179B (zh) 数字正电子发射断层摄影中的飞行时间校准
US10215864B2 (en) System and method to improve image quality of emission tomography when using advanced radionuclides
JP7317586B2 (ja) 医用画像処理装置、方法及びプログラム
US8507842B2 (en) Method for identifying 3-D location of gamma interaction and flat panel gamma imaging head apparatus using the same
Bailey Data acquisition and performance characterization in PET
CN103547942A (zh) 随机事件减少方法、随机事件减少装置以及非暂时性计算机可读存储介质
US11543545B2 (en) Method and apparatus to use a broad-spectrum energy source to correct a nonlinear energy response of a gamma-ray detector
CN111736207B (zh) 一种pet时间标定方法
Ota et al. Direct positron emission imaging: ultra-fast timing enables reconstruction-free imaging
US7304307B2 (en) PMT signal correlation filter
Gu et al. A new pulse pileup rejection method based on position shift identification
CN112932517A (zh) 本底放射性实现时间自校正方法及装置
Kijewski Positron emission tomography (PET) and single-photon emission computed tomography (SPECT) physics
CN111714147A (zh) 能量刻度曲线获取方法、装置、计算机设备和存储介质
CN211834429U (zh) 全模拟核素分布显像系统
CN113253330B (zh) 伽玛射线放射成像装置及能量校准方法
CN111685788B (zh) 一种提高pet信噪比的方法
CN117442224A (zh) 一种用于符合探测系统时间校正的方法
Tayefi Ardebili Evaluation of the NEMA characteristics for the Modular J-PET scanner
Soares Characterizing Low-Activity Imaging for Applications in Stem Cell Tracking with the MiniEXPLORER II PET/CT Scanner

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