CN112882085A - 温度自适应的pet探测器能量修正方法、系统及计算机可读存储介质 - Google Patents

Abstract

本发明公开了一种温度自适应的PET探测器能量修正方法、系统及计算机可读存储介质，所述方法包括：获取能量和道址的线性函数关系E＝a*x+b和道址和温度的线性函数关系P＝c*T+d；结合所述能量和道址的线性函数关系、所述道址和温度的线性函数关系获取能量和温度、道址的线性函数关系；获取所述PET探测器的实时击中事例的温度信息和能量道址信息，并参照所述能量和温度、道址的线性函数关系获取所述温度信息、能量道址信息所对应的能量修正值，用所述能量修正值替换所述实时击中事例的能量信息，实现对温度和能量的在线实时修正。

Description

温度自适应的PET探测器能量修正方法、系统及计算机可读存 储介质

技术领域

本发明涉及正电子发射计算机断层扫描成像系统技术领域，尤其涉及一种温度自适应的PET探测器能量修正方法、系统及计算机可读存储介质。

背景技术

正电子发射计算机断层扫描成像系统(PET，Positron Emission Tomography)是一种核医学影像设备，主要由探测器系统、电子学系统、数据采集系统和重建系统组成。其中探测器系统由闪烁晶体、光电转换器件和前端电子学构成。目前，基于硅光电倍增管(Silicon Photomultiplier,SiPM)的探测器件由于其良好的能量和时间分辨率以及磁兼容性能越来越多地被应用在PET系统中。其主要工作原理是利用SiPM将闪烁晶体捕获的高能Gamma光子转化成低能可见光，再通过光电效应将其转化为模拟电信号，并利用模拟调理电路将模拟电信号进行放大成形后送入能量、时间和位置测量装置得到该电信号的能量、时间和位置信息，然后利用后端的符合判选等方法筛选出有效信号。

在能量测量装置中，为了排除散射事例，提高信噪比，进而提高图像对比度，通常会选择设置能量窗来挑选能量事例，尽可能将511keV全能峰置于能量窗内而排除掉散射事件。然而对于一台PET系统来说，信号通道多达几万个，各像素之间增益并不是均匀的，因此在设置同一能窗挑选事例前，通常会进行能量修正将所有像素的511keV全能峰修正到相同增益下的同一位置。另外，SiPM是半导体探测器，其性能受温度影响大，温度上升，SiPM增益变小；温度下降，SiPM增益变大。SiPM增益的改变会导致能量峰位漂移，进而提高散射事件在预设能窗内的占比，影响图像信噪比。更为严重的，如果能峰漂移较大，能窗内计算将急剧减少，还会影响探测器灵敏度。

对于能峰漂移，通常会考虑在探测器电路中增加温度反馈电路，通过温度来实时调节SiPM偏压，稳定SiPM增益，但温度反馈电路较为复杂，且会增加成本，因此还是有很多PET系统没有温度调节电路。在这些系统中，现有技术主要通过外部冷却系统，来实现对SiPM的温度控制。但外部冷却系统只能保证SiPM的温度控制在某一个范围内，而不能很精确的将SiPM的温度一直限定在某一个温度，在这个温度范围内，能量也会有一定程度的漂移。尤其是在PET系统开机后，系统从开机到达到稳定工作温度需要很长一段时间，这段时间内温度变化很大，能量漂移严重。另外也可以通过上位机来分析SiPM温度和能量的响应，在每次开机或者每次采集时，根据响应系数重新计算能量查找表并下发至硬件电路，来得到修正的目的，但这种方法实时性差，响应速度较慢，效果不佳。

发明内容

为了克服上述技术缺陷，本发明的目的在于提供一种在不增加额外电路的基础上实现温度反馈、同时修正温度和能量的温度自适应的PET探测器能量修正方法、系统及计算机可读存储介质。

本发明公开了一种温度自适应的PET探测器能量修正方法，包括如下步骤：获取所述PET探测器在工作状态下的各像素的两个全能峰峰位，根据所述两个全能峰峰位的道址获取能量和道址的线性函数关系E＝a*x+b，其中，E为能量值；获取所述PET探测器在不同温度下的各像素的某一全能峰峰位对应的道址，建立道址和温度的线性函数关系P＝c*T+d，其中，P为道址数；结合所述能量和道址的线性函数关系、所述道址和温度的线性函数关系获取能量和温度、道址的线性函数关系；获取所述PET探测器的实时击中事例的温度信息和能量道址信息，并参照所述能量和温度、道址的线性函数关系获取所述温度信息、能量道址信息所对应的能量修正值，用所述能量修正值替换所述实时击中事例的能量信息。

优选地，所述获取所述PET探测器在工作状态下的各像素的两个全能峰峰位，根据所述两个全能峰峰位的道址获取能量和道址的线性函数关系E＝a*x+b包括：在第一工作温度T₀下，设置时间窗和视场，并采集所述PET探测器的所对应的两种能量射线的数据包；获取所述数据包中各像素的能量分布，从所述能量分布中获取两个全能峰峰位所对应的道址；根据所述两个全能峰峰位对应的射线能量值和所述两个全能峰峰位所对应的道址获取系数

其中，s₁、s₂分别为两个全能峰峰位对应的射线真实能量值，u₁、u₂分别为两个全能峰峰位所对应的道址；根据系数a和系数b获取能量和道址的线性函数关系E＝a*x+b。

优选地，所述获取所述PET探测器在不同温度下的各像素的某一全能峰峰位对应的道址，建立道址和温度的线性函数关系P＝c*T+d包括：获取若干个工作温度下各像素能量分布中一预设能量射线的全能峰峰位对应的道址；获取若干个工作温度下所述PET探测器的实际温度；获取所述PET探测器的在不同温度下的同一全能峰峰位道址，建立道址和温度的线性函数关系。

优选地，所述获取若干个工作温度下各像素能量分布中一预设能量射线的全能峰峰位对应的道址包括：在温度T下，将一预设能量射线的放射源置于FOV中心，采集单事例数据包；解析单事例数据包，获取各像素的能量分布；获取各像素能量分布中全能峰峰位对应的道址数P；以预设温度差为间隔取若干个工作温度T₀、T₁、...、T_n，获取若干个工作温度下各像素的同一全能峰峰位对应的道址P₀、P₁、...、P_n。

优选地，所述结合所述能量和道址的线性函数关系、所述道址和温度的线性函数关系获取能量和温度、道址的线性函数关系为：E_i＝a_i*A_measure+a_i*c_i*(T0_β-T_α)+b_i；其中，β为PET探测器单元的编号，i为β号PET探测器单元中的某一像素编号，A_measure为像素i实际测量的能量的道址，T0_β为第一工作温度T₀下β号PET探测器单元的实测温度，T_β为探测器采集实时击中事例时β号PET探测器单元的实际温度。

优选地，所述获取所述PET探测器的实时击中事例的温度信息和能量道址信息，并参照所述能量和温度、道址的线性函数关系获取所述温度信息、能量道址信息所对应的能量修正值包括：实时获取所述PET探测器单元的温度，并取平均值；所述PET探测器实时捕获击中事例，将所述击中事例中的gamma光子的能量转化为模拟电信号，所述模拟电信号经整形滤波并通过模数转换器后进行能量积分获得能量道址值；根据PET探测器被击中的像素编号，查找该像素对应的修正函数参数，利用PET探测器模块的实时温度信息和击中事例的能量道址信息，根据所述能量和温度、道址的线性函数关系计算当前击中事例的能量在当前温度下所对应的能量修正值。

优选地，所述根据PET探测器被击中的像素编号，查找该像素对应的修正函数参数，利用PET探测器模块的实时温度信息和击中事例的能量道址信息，根据所述能量和温度、道址的线性函数关系计算当前击中事例的能量在当前温度下所对应的能量修正值包括：将a_i、b_i、c_i、T0_β作为查找对象，将像素编号作为参考对象建立查找表；将查找到的与像素编号相对应的a_i、b_i、c_i、T0_β代入能量和温度、道址的线性函数关系，获取能量修正值。

本发明还公开了一种温度自适应的PET探测器能量修正系统，位于所述PET探测器的FPGA模块内，包括修正函数查找表单元、位置测量单元、能量测量单元、温度测量单元、修正逻辑单元、模数转化电路和存储单元；所述修正函数查找表单元用于储存探测器各像素的修正函数参数a_i、b_i、c_i、T0_β，每个像素的修正函数参数以像素编号为检索地址存储于所述存储单元中；所述位置测量单元获取伽马光子击中事例中所述PET探测器的像素位置；所述能量测量单元将实时击中探测器的伽马光子的能量转化为模拟电信号，模拟信号经整形滤波并通过模数转换后送入所述FPGA模块内进行积分处理获得事例的能量道址数；所述温度测量单元通过模数转化电路实时读取探测器单元上的温度信息，并取平均值作为输出的温度信息；所述修正逻辑单元根据所述探测器单元的实时温度信息和伽马光子击中像素沉积能量的道址数、并参照所述能量和温度以及道址的响应函数关系计算对应的能量修正值，并用所述能量修正值替换所述实时击中事例的能量信息。

本发明还公开了一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现上述任一所述的方法的步骤。

采用了上述技术方案后，与现有技术相比，具有以下有益效果：

1.在不额外增加反馈电路的基础上，实现了探测器的温度的反馈，成本较低；

2.直接在所述PET探测器的FPGA模块内进行，无需通过所述PET探测器的计算机系统操作，反应速度快，效果佳，可在线完成温度修正和能量修正。

附图说明

图1为本发明提供的温度自适应的PET探测器能量修正方法的整体流程图；

图2为本发明提供的温度自适应的PET探测器能量修正方法的步骤S4的流程图；

图3为本发明提供的温度自适应的PET探测器能量修正方法的PET探测器模块像素i的LYSO本底符合能量分布图；

图4为本发明提供的温度自适应的PET探测器能量修正系统的PET探测器模块像素i的peak-温度响应曲线。

具体实施方式

以下结合附图与具体实施例进一步阐述本发明的优点。

这里将详细地对示例性实施例进行说明，其示例表示在附图中。下面的描述涉及附图时，除非另有表示，不同附图中的相同数字表示相同或相似的要素。以下示例性实施例中所描述的实施方式并不代表与本公开相一致的所有实施方式。相反，它们仅是与如所附权利要求书中所详述的、本公开的一些方面相一致的装置和方法的例子。

在本公开使用的术语是仅仅出于描述特定实施例的目的，而非旨在限制本公开。在本公开和所附权利要求书中所使用的单数形式的“一种”、“所述”和“该”也旨在包括多数形式，除非上下文清楚地表示其他含义。还应当理解，本文中使用的术语“和/或”是指并包含一个或多个相关联的列出项目的任何或所有可能组合。

应当理解，尽管在本公开可能采用术语第一、第二、第三等来描述各种信息，但这些信息不应限于这些术语。这些术语仅用来将同一类型的信息彼此区分开。例如，在不脱离本公开范围的情况下，第一信息也可以被称为第二信息，类似地，第二信息也可以被称为第一信息。取决于语境，如在此所使用的词语“如果”可以被解释成为“在……时”或“当……时”或“响应于确定”。

在本发明的描述中，需要理解的是，术语“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。

在本发明的描述中，除非另有规定和限定，需要说明的是，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是机械连接或电连接，也可以是两个元件内部的连通，可以是直接相连，也可以通过中间媒介间接相连，对于本领域的普通技术人员而言，可以根据具体情况理解上述术语的具体含义。

在后续的描述中，使用用于表示元件的诸如“模块”、“部件”或“单元”的后缀仅为了有利于本发明的说明，其本身并没有特定的意义。因此，“模块”与“部件”可以混合地使用。

参见附图1，本发明公开了一种温度自适应的PET探测器能量修正方法，包括如下步骤：

S1、获取PET探测器在工作状态下的各像素的两个全能峰峰位，根据两个全能峰峰位的道址获取能量和道址的线性函数关系E＝a*x+b，其中，E为能量值；

S2、获取PET探测器在不同温度下的各像素的某一全能峰峰位对应的道址，建立道址和温度的线性函数关系P＝c*T+d，其中，P为道址数；

S3、结合能量和道址的线性函数关系、道址和温度的线性函数关系获取能量和温度、道址的线性函数关系；

S4、获取PET探测器的实时击中事例的温度信息和能量道址信息，并参照能量和温度、道址的线性函数关系获取温度信息、能量道址信息所对应的能量修正值，用能量修正值替换实时击中事例的能量信息。

步骤S1具体包括：在第一工作温度T₀下，设置时间窗和视场角，并采集PET探测器所对应的两种能量射线的数据包；获取数据包中各像素的能量分布，从能量分布中获取两个全能峰峰位所对应的道址；根据两个全能峰峰位对应的射线能量值和两个全能峰峰位所对应的道址获取系数

其中，s₁、s₂分别为两个全能峰峰位对应的射线能量值，u₁、u₂分别为两个全能峰峰位所对应的道址；根据系数a和系数b获取能量和道址的线性函数关系E＝a*x+b。

步骤S2具体包括：在温度T下，将某一特定能量射线放射源置于FOV中心，采集单事例数据包；解析单事例数据包，获取各像素的能量分布；获取各像素能量分布中全能峰峰位对应的道址数P；以预设温度差为间隔取若干个工作温度T₀、T₁、...、T_n，获取若干个工作温度下各像素在该特定能量射线下的全能峰峰位对应的道址P₀、P₁、...、P_n；获取若干个工作温度T₀、T₁、...、T_n下的PET探测器单元的实时温度T_0β、T_1β、...、T_nβ，β为探测器单元编号；根据PET探测器各像素在不同温度T_0β、T_1β、...、T_nβ下的道址建立曲线关系，根据曲线关系获取道址和温度的线性函数关系P＝c*T+d，其中，c、d为曲线关系中获取的函数系数。

根据线性函数关系E＝a*x+b和线性函数关系P＝c*T+d获取步骤S3中的能量、温度、道址的线性函数关系为：E_i＝a_i*A_measure+a_i*c_i*(T0_β-T_α)+b_i；其中，β为PET探测器单元的编号，i为β号PET探测器中的某一像素编号，A_measure为像素i实际测量的能量的道址，T0_β为第一工作温度T₀下β号PET探测器单元的实测温度，T_β为探测器采集实时击中事例时β号PET探测器单元的实际温度。

参见附图2，步骤S4具体包括：实时获取PET探测器单元的温度，并在预设时间段内取平均值以获得稳定的PET探测器单元温度；PET探测器实时捕获击中事例，将击中事例中的gamma光子的能量转化为模拟电信号，模拟电信号经整形滤波并通过模数转换器后进行能量积分获得能量道址值；将a_i、b_i、c_i、T0_β作为查找对象，将像素编号作为参考对象建立查找表，并存储于FPGA内的RAM资源中；通过伽马事例实时击中探测器像素的位置编号从RAM中查询对应的函数修正参数a_i、b_i、c_i、T0_β代入能量和温度、道址的线性函数关系，通过实时捕获的击中事例的能量道址以及模块单元温度计算能量修正值，用能量修正值替换实时击中事例的能量信息，实现对能量的实时在线修正。

本发明还公开了一种温度自适应的PET探测器能量修正系统，包括：

-修正函数查找表单元，用于储存探测器各像素的修正函数参数a_i、b_i、c_i、T0_β，每个像素的修正函数参数以像素编号为检索地址存储于存储单元中；

-位置测量单元，获取伽马光子击中事例中PET探测器的像素位置；

-能量测量单元，将实时击中探测器的伽马光子的能量转化为模拟电信号，模拟信号经整形滤波并通过模数转换后送入FPGA模块内进行积分处理获得事例的能量道址数；

-温度测量单元，通过模数转化电路实时读取探测器单元上的温度信息，并取平均值作为输出的温度信息；

-修正逻辑单元，根据探测器单元的实时温度信息和伽马光子击中像素沉积能量的道址数、并参照能量和温度以及道址的响应函数关系计算对应的能量修正值，并用能量修正值替换实时击中事例的能量信息。

本发明提供一种优选实施例，首先获取能量修正系数：将环境温度设置为T₀，待PET系统的温度稳定后，控制系统通过接收指令读取PET探测器整环的各探测器单元的温度T_0β，并保存至存储器，β为探测器单元的编号。本优选实施例中的探测器采用LYSO闪烁体，对于采用其他无本底放射性的闪烁晶体(如BGO)的PET探测器，可使用放射源标定能量响应曲线。在T₀温度下，设置时间窗和FOV的参数，采集LYSO本底符合数据包；解析符合数据包，提取各符合事件中两个gamma信号的能量和击中位置信息，以像素位置为单位，统计各像素的能量分布；通过202keV和307keV两个本底峰获取各像素能量的响应曲线E＝a*x+b，如附图3所示，获取各像素能量分布中202keV和307keV全能峰峰位(peak)对应的模拟数字转换器ADC的道址数u₁和u₂，得到的能量-ADC道址对应关系，并使用各像素能量的响应曲线，获得各像素的能量修正参数a和b，将各像素能量修正参数a和b以像素位置为单位保存至存储器。

接着获取温度修正系数：在T₀温度下，将Na22，G68等正电子源的放射源放置于FOV中心，采集单事例数据；解析单事例数据包，统计各像素的能量分布；获取各像素能量分布中511keV峰对应的ADC道址数P_0i，i为像素索引编号，将P_0i以像素位置为单位保存至存储器；将环境温度设置为T₀+2℃，记为T₁，获取各像素511keV全能峰对应的ADC道址数P_1i，并保存至存储器；以2℃为温差间隔重复获取T₂、T₃、T₄、T₅、T₆温度下各像素511keV全能峰对应的ADC道址数P_2i、P_3i、P_4i、P_5i、P_6i，并保存至存储器。在环境温度T₀、T₁、T₂、T₃、T₄、T₅、T₆下待PET系统温度稳定后，读取探测器的整环上的各探测器单元模块温度T_0β、T_1β、T_2β、T_3β、T_4β、T_5β、T_6β，并保存至存储器。参见附图4，以像素位置为单位，描绘511keV峰的(P_0i～P_6i)对温度(T_0β～T_6β)的响应曲线，从图中可见此响应曲线随温度变高，峰值呈线性减小，根据该响应曲线，获得各像素的温度反馈修正系数，以附图4中的像素为例，c为-0.229±0.01202、d为136.8±0.3146，最后将各像素的温度反馈修正系数以像素编号为查找的参考对象保存至存储器。

接着获取能量-温度联合修正查找表：根据能量修正系数和温度修正系数a、b、c、T0_β，可得各像素的能量和温度的联合修正函数为：E_i＝a_i*A_measure+a_i*c_i*(T0_β-T_α)+b_i；将a_i、b_i、c_i、T0_β以像素编号为查找的参考对象做成查找表LUT，发送保存至PET探测器前端的FPGA模块内的RAM资源中。

最后根据能量修正查找表修正正常工作状态下的能量值：系统启动数据采集后，FPGA通过模数转化电路实时读取单位模块上温度传感器信息，并通过取平均的方式获得稳定的温度信息；PET探测器实时捕获伽马光子事例，将gamma光子能量转化为模拟电信号，模拟信号经整形滤波并通过模数转换器后传输到FPGA中积分获得原始能量道址值。将gamma光子击中PET探测器的像素位置作为RAM地址索引号，从LUT表中查寻得到该像素的能量-温度修正系数，结合温度信息和信号能量道址、利用能量和温度的联合修正函数，在FPGA模块内对原始能量值进行实时修正，并输出修正后的能量值。

本发明还公开了一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现上述任一的方法的步骤。

应当注意的是，本发明的实施例有较佳的实施性，且并非对本发明作任何形式的限制，任何熟悉该领域的技术人员可能利用上述揭示的技术内容变更或修饰为等同的有效实施例，但凡未脱离本发明技术方案的内容，依据本发明的技术实质对以上实施例所作的任何修改或等同变化及修饰，均仍属于本发明技术方案的范围内。

Claims (9)

1.一种温度自适应的PET探测器能量修正方法，其特征在于，包括如下步骤：

获取所述PET探测器在工作状态下的各像素的两个全能峰峰位，根据所述两个全能峰峰位的道址获取能量和道址的线性函数关系E＝a*x+b，其中，E为能量值；

获取所述PET探测器在不同温度下的各像素的某一全能峰峰位对应的道址，建立道址和温度的线性函数关系P＝c*T+d，其中，P为道址数；

结合所述能量和道址的线性函数关系、所述道址和温度的线性函数关系获取能量和温度、道址的线性函数关系；

获取所述PET探测器的实时击中事例的温度信息和能量道址信息，并参照所述能量和温度、道址的线性函数关系获取所述温度信息、能量道址信息所对应的能量修正值，用所述能量修正值替换所述实时击中事例的能量信息。

2.根据权利要求1所述的PET探测器能量修正方法，其特征在于，所述获取所述PET探测器在工作状态下的各像素的两个全能峰峰位，根据所述两个全能峰峰位的道址获取能量和道址的线性函数关系E＝a*x+b包括：

在第一工作温度T₀下，设置时间窗和视场，并采集所述PET探测器的所对应的两种能量射线的数据包；

获取所述数据包中各像素的能量分布，从所述能量分布中获取两个全能峰峰位所对应的道址；

根据所述两个全能峰峰位对应的射线能量值和所述两个全能峰峰位所对应的道址获取系数

其中，s₁、s₂分别为两个全能峰峰位对应的射线真实能量值，u₁、u₂分别为两个全能峰峰位所对应的道址；

根据系数a和系数b获取能量和道址的线性函数关系E＝a*x+b。

3.根据权利要求1所述的PET探测器能量修正方法，其特征在于，所述获取所述PET探测器在不同温度下的各像素的某一全能峰峰位对应的道址，建立道址和温度的线性函数关系P＝c*T+d包括：

获取若干个工作温度下各像素能量分布中一预设能量射线的全能峰峰位对应的道址；

获取若干个工作温度下所述PET探测器的实际温度；

获取所述PET探测器的在不同温度下的同一全能峰峰位道址，建立道址和温度的线性函数关系。

4.根据权利要求3所述的PET探测器能量修正方法，其特征在于，所述获取若干个工作温度下各像素能量分布中一预设能量射线的全能峰峰位对应的道址包括：

在温度T下，将一预设能量射线的放射源置于FOV中心，采集单事例数据包；

解析单事例数据包，获取各像素的能量分布；

获取各像素能量分布中全能峰峰位对应的道址数P；

以预设温度差为间隔取若干个工作温度T₀、T₁、...、T_n，获取若干个工作温度下各像素的同一全能峰峰位对应的道址P₀、P₁、...、P_n。

5.根据权利要求2所述的PET探测器能量修正方法，其特征在于，所述结合所述能量和道址的线性函数关系、所述道址和温度的线性函数关系获取能量和温度、道址的线性函数关系为：E_i＝a_i*A_measure+a_i*c_i*(T0_β-T_α)+b_i；

其中，β为PET探测器单元的编号，i为β号PET探测器单元中的某一像素编号，A_measure为像素i实际测量的能量的道址，T0_β为第一工作温度T₀下β号PET探测器单元的实测温度，T_β为探测器采集实时击中事例时β号PET探测器单元的实际温度。

6.根据权利要求1所述的PET探测器能量修正方法，其特征在于，所述获取所述PET探测器的实时击中事例的温度信息和能量道址信息，并参照所述能量和温度、道址的线性函数关系获取所述温度信息、能量道址信息所对应的能量修正值包括：

实时获取所述PET探测器单元的温度，并取平均值；

所述PET探测器实时捕获击中事例，将所述击中事例中的gamma光子的能量转化为模拟电信号，所述模拟电信号经整形滤波并通过模数转换器后进行能量积分获得能量道址值；

根据PET探测器被击中的像素编号，查找该像素对应的修正函数参数，利用PET探测器模块的实时温度信息和击中事例的能量道址信息，根据所述能量和温度、道址的线性函数关系计算当前击中事例的能量在当前温度下所对应的能量修正值。

7.根据权利要求6所述的PET探测器能量修正方法，其特征在于，所述根据PET探测器被击中的像素编号，查找该像素对应的修正函数参数，利用PET探测器模块的实时温度信息和击中事例的能量道址信息，根据所述能量和温度、道址的线性函数关系计算当前击中事例的能量在当前温度下所对应的能量修正值包括：

将a_i、b_i、c_i、T0_β作为查找对象，将像素编号作为参考对象建立查找表；

将查找到的与像素编号相对应的a_i、b_i、c_i、T0_β代入能量和温度、道址的线性函数关系，获取能量修正值。

8.一种温度自适应的PET探测器能量修正系统，其特征在于，位于所述PET探测器的FPGA模块内，包括修正函数查找表单元、位置测量单元、能量测量单元、温度测量单元、修正逻辑单元、模数转化电路和存储单元；

所述修正函数查找表用于储存探测器各像素的修正函数参数a_i、b_i、c_i、T0_β，每个像素的修正函数参数以像素编号为检索地址存储于所述存储单元中；

所述位置测量单元获取伽马光子击中事例中所述PET探测器的像素位置；

所述能量测量单元将实时击中探测器的伽马光子的能量转化为模拟电信号，模拟信号经整形滤波并通过模数转换后送入所述FPGA模块内进行积分处理获得事例的能量道址数；

所述温度测量单元通过模数转化电路实时读取探测器单元上的温度信息，并取平均值作为输出的温度信息；

所述修正逻辑单元根据所述探测器单元的实时温度信息和伽马光子击中像素沉积能量的道址数、并参照所述能量和温度以及道址的响应函数关系计算对应的能量修正值，并用所述能量修正值替换所述实时击中事例的能量信息。

9.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现权利要求1-7中任一所述的方法的步骤。

Images