2023年11月30日发(作者：)
PET时间分辨率测算、重建系统和方法
技术领域
本发明涉及一种PET时间分辨率测算、重建系统和方法，属于X射线计算
机断层扫描成像技术领域。
背景技术
正电子发射断层设备的原理是测量正电子湮没后的γ光子配对，并利用的γ
光子对的几何方向和时间特性进行成像。其中γ光子对的时间特性是指，两个方
向相反的光子由于飞行时间不同，进入检测器闪烁体晶体的时间产生差异。这个
时间差会受到晶体材质，晶体形态和前、后端电子学性能的影响，在测量时产生
一定的展宽。这个展宽就是时间分辨率的主要特征。只有计算出精准的时间分辨
率，才能够更好地提升系统成像质量。
现有技术中，为了定量地测量时间差展宽(后称为时间分辨率)，一般采用的
技术是利用点源或线源来查看光子计数在不同LOR上的分布。这类方法有着一
定局限性：
(1)一般只能算出系统的总体时间分辨率，无法精确计算到每一个晶体的时间
分辨率。
(2)对于辐射源的形状摆位有严格要求，无法实时计算。
基于此，做出本申请。
发明内容
为了解决现有技术中存在的上述缺陷，本发明提供了一种PET时间分辨率
测算、重建系统和方法。
为了实现上述目的，本发明采取的技术方案如下：
本发明使用任意源对TOF-PET进行分辨率的计算。任意源是指可以是不规
则并具有一定体积的固体或液体，也可以是真实的扫描病人或者其他动物，物体
放置无特定要求，可以放置在检测器环的中心也可以偏心放置。因此本发明具有
很好的实用性，同时也具有很高的准确性。PET诊断成像装置，包括：成像区
域外部的多个辐射探测器环状结构，用来采集到发出的辐射事件的信息。探测器
具有不同的定时延迟。校准体模可以发射时间同时但方向相反的辐射的事件对
(符合事件)的辐射源，辐射事件对与辐射探测器的相互作用路径是LOR，当检测
器环确定后，每一条LOR的几何位置是已知的。一旦每个检测器的定时延迟被
校准，可以认为符合事件的标记时间是在真实符合发生的位置为中心的，并以特
定宽度形成高斯分布。这个特定的宽度(特征展宽)就是本发明要求出的晶体时间
分辨率。
本发明一种PET时间分辨率测算、重建方法，包括如下步骤：
(1)PET原始符合数据采集；
(2)辐射源空间分布获取，对PET原始符合数据进行无飞行时间的重建，得
到光子空间分布图以及光子衰减分布图；目前PET同时借助CT或者MR的成
像结果，产生γ光子在空间中的衰减信息的空间分布，即光子衰减分布图。光子
空间分布图是指PET数据进行重建后，形成的2维或3维空间内强度的分布。
这个强度的分布代表了扫描物体在空间中不同位置的发射γ光子强度。
(3)符合事件累积，对PET原始符合数据进行统计，得到每个晶体的时间直
方图分布H1；进行统计的意思是根据符合事件内的晶体标号，找出同一标号下
的所有符合事件，提取出符合事件的时间差信息，对时间差信息按照特定的时间
间隔dt进行直方图计算。
(4)LOR线追踪累积，通过光子空间分布图计算得到每个晶体的时间直方图
分布H2；线追踪(ray-tracing)方法是指，对所有从特定标号的晶体出发的所有
LOR的方向，按照不同光子的飞行时间，得到多个终点位置，并对这些位置的
像素值求和。
(5)特征宽度计算，根据步骤(3)和步骤(4)得到的H1和H2计算得到第k个晶
体的特征展宽σ 为第k个晶体的时间分辨率；
k
(6)迭代更新，对步骤(5)得到的时间分辨率σ 和上一次时间分辨率σ 进行
k0
更新迭代，得到更新后的σ’ ，比较σ 和σ’ 的差值，判断是否满足设定值，
k0k
若不满足设定值，将σ’ 回传至步骤(4)，进行下一次更新，循环往复直至满足
k
设定值，跳出循环，得到最终的结果σ ；
k
(7)重建，利用重建公式并根据光子衰减分布图对σ 进行重建。
k
本发明的原理和有益技术效果：
1、本发明在对计数统计的同时，还对重建图像进行空间强度计算，通过两
者的相似度来精确计算每一条LOR上的准确时间分辨率。
2、本发明对于任意γ光源，可以精准测试晶体的时间分辨率：本发明方法，
γ光源不管是何种形态，都会在步骤(2)获取到它的空间分布信息即形态信息。这
样，再通过紧密相关的(3)(4)(5)(6)步骤的处理，就可以精准测试晶体的时间分辨
率。
3、对于PET的检测器时间分辨率差异很大时，也可以精准成像：本发明方
法可以精准计算每一个晶体的时间分辨率。同时结合步骤(7)的重建方法，实现
精准重建成像。
因此，相较于传统方法，本发明不仅可以分辨率的计算可以具体到每一个晶
体，而且对于辐射源形态无要求，可以实时测量。
附图说明
图1为本实施例PET时间分辨率测算、重建系统的流程图。
具体实施方式
为了使本发明的技术手段及其所能达到的技术效果，能够更清楚更完善的披
露，兹提供了以下实施例，并结合附图作如下详细说明：
本发明使用任意源对TOF-PET进行分辨率的计算。任意源是指可以是不规
则并具有一定体积的固体或液体，也可以是真实的扫描病人或者其他动物，物体
放置无特定要求，可以放置在检测器环的中心也可以偏心放置。因此本发明具有
很好的实用性，同时也具有很高的准确性。PET诊断成像装置，包括：成像区
域外部的多个辐射探测器环状结构，用来采集到发出的辐射事件的信息。探测器
具有不同的定时延迟。校准体模可以发射时间同时但方向相反的辐射的事件对
(符合事件)的辐射源，辐射事件对与辐射探测器的相互作用路径是LOR，当检测
器环确定后，每一条LOR的几何位置是已知的。一旦每个检测器的定时延迟被
校准，可以认为符合事件的标记时间是在真实符合发生的位置为中心的，并以特
定宽度形成高斯分布。这个特定的宽度(特征展宽)就是本发明要求出的晶体时间
分辨率。
本实施例的一种PET时间分辨率测算、重建方法，包括如下步骤：
(1)PET原始符合数据采集(通过PET原始符合数据采集模块)；PET探测系
统由外罩，机架及探测器组成。探测器包括晶体、光电放大和电子处理电路。探
测器检测从被检体P内部放出的成对湮没γ射线，生成与检测出的成对湮没γ
射线的光量相应的脉冲状电信号。后端信号处理器根据来自光电转换器件的电信
号生成单事件数据(Single)。计数部从重复供给的单事件数据中重复确定容纳在
与预先设定的时间范围内的2个单事件有关的事件数据。时间范围被设定为例如
6ns～18ns左右。该成对的单事件被推测为由来于从同一成对湮没点产生的成对
湮没γ射线。成对的单事件概括地被称为符合事件。符合事件被收集到后，事件
相关的时间空间信息进行组合封装，并在磁盘或者其他数据存储介质上保存，形
成了重建程序可以读取并使用的PET符合数据。有符合事件中的2个单事件，
会各自被被检测器记录的到达检测器时间点，这个2个单事件时间点的差值记为
时间差。符合数据根据PET机器分为两种，一种是含的上述时间差的，一种是
不含有时间差的。本次发明是针对的是有时间差信息的符合事件。
(2)辐射源空间分布获取(通过辐射源空间分布获取模块)，对PET原始符合数
据进行无飞行时间的重建，得到光子空间分布图以及光子衰减分布图；具体来说，
本实施例中它的作用是提取辐射源的实际形态特征，这样可以获得到辐射源的空
间分布。它是利用PET诊断成像装置获取到PET符合数据信号，对信号进行无
飞行时间的重建，重建方法例如，有序子集最大似然期望值法OSEM，最大似
然期望值法(MLEM)，滤波反投影FBP或者其他方法，从而得到了基本除去噪声
和衰减干扰的质量较佳的光子空间分布图。同时PET系统借助CT或者MR，产
生γ光子在空间中的衰减信息的空间分布图。
(3)符合事件累积(通过符合事件累积模块)，对PET原始符合数据进行统计，
得到每个晶体的时间直方图分布H1[c,h]。这里的c代表了晶体的序号，h代表
了出现在该晶体上的事件所记录的时间差。例如，当符合事件i发生时，与它相
关的晶体为c1和c2，时差为ti，则分别在H1[c1,m-ti]和H2[c2,m+ti]上记录，这
里的m是指的H1的直方图长度的1/2。
(4)LOR线追踪累积(通过LOR线追踪累积模块)，从图像的角度(根据光子空
间分布图)计算得到每个晶体的时间直方图分布H2[c,h]，这里的c代表的晶体的
序号，h代表了出现在该晶体上的事件所记录的时间差。H2的记录规则是:统计
所有从晶体c出发的LOR线，其中任意一条LOR线上，记录LOR上每一个位
置的图像强度Intens和距离中心的位置Dis，通过光速与距离的换算因子c，得
到H2上更新的位置是Dis/c。定义射线追踪法获取到晶体c，时间偏移h的强度
为H2，则
公式中，Voxel为步骤(2)计算出的体数据(即光子空间分布)，LORc是与晶
体c相关的所有LOR。x (h),y (h),z (h)是晶体c的第i条LOR上时间差为h的
iii
三维坐标，可以通过晶体的几何位置快速推导。σ′ 是与c的第i个LOR上配对
i
晶体c ’的特征宽度。
i
(5)特征宽度计算(通过特征宽度计算模块)，根据步骤(3)和步骤(4)得到的H1
和H2计算得到第k个晶体的特征展宽σ 为第k个晶体的时间分辨率；
k
本模块通过卷积H2，计算第k个晶体的特征展宽为
其中优化求解最大值算法argmax一般使用线搜索(line search)实现。||代表
计算向量的范数。
(6)迭代更新(通过迭代更新模块)，对步骤(5)得到的时间分辨率σ 和上一次
k
时间分辨率σ 进行更新迭代，得到更新后的σ’ ，比较σ 和σ’ 的差值，判断
0k0k
是否满足设定值，若不满足设定值，将σ’ 回传至步骤(4)，进行下一次更新，
k
循环往复直至满足设定值，跳出循环，得到最终的结果σ ；
k
具体地：
σ’ ＝∈σ +(1-∈)σ
kk0k
σ’ 代表更新后的晶体分辨率宽度，σ 代表上一次分辨率，σ 代表步骤(5)
kk0k
计算的分辨率。∈是经验参数，代表了更新的速度。
更新后的σ’ 会回传给步骤(4)，进行下一次更新。通过比较σ 和σ ，之
kk0k
间差异小到某一预设定值时，跳出循环。得到最终的结果σ 。
k
计算完成后，对于任意一个下标号为k的晶体，都能获得一个σ 与之对应。
k
(7)重建(通过重建模块)，利用重建公式并根据光子衰减分布图对σ 进行重
k
建。
重建方法是含有以下部分或者全部特征的公式：
f代表图像，它的下标代表该图像的像素标号，第一个上标k代表迭代次数，
第二个上标m代表子集的序号。N是图像总的像素点个数。
α为衰减系数，其下标j代表了具体的事件标号。它是通过(1)中所述的衰减
信息的空间分布图计算的。
Lm代表了子集m中的事件序号的集合。
是系统矩阵，它反映了事件j和像素点i之间几何对应关系。
是散射和随机事件形成的校正因子。它是由符合数据通过公知的散射估
计和随机估计计算出来的。
需要说明的是：这里 的计算都是常规方法，有文献支持，不展开说明。
其中本实施例对该公式中的系统矩阵项H的权重分布进行了调整。周知一般
情况下， 含有系统几何信息和飞行时间信息,其中飞行时间信息取决于飞行
时间分辨率。PET系统的时间分辨率为2.236*σ，则 的飞行时间信息的分
布服从[0，σ]的正态分布。本实施例中则是对于任意一个序号为i 符合事件，计
k
算实际的正态分布宽度 其中， 和 是第i 的符合事件所对应的晶
k
体的特征宽度(这里涉及到几个下标：i ，i ，i 和k，步骤6中的k指晶体序
kk1k2
号，i 是符合事件标号。任意一个符合事件i 是2个单事件组成，接收到这2
kk
个单事件的晶体序号为i ,i ，将i ,i 带入(6)的σ 的下标k中，就可以获
k1k2k1k2k
得 最后再计算出 就是符合事件的宽度)。
同时，关于“ 含有系统几何信息和飞行时间信息”，可参考相关文献的描
述(Phys Med Biol.2014 February 7；59(3):541–
:10.1088/0031-9155/59/3/541.)。
其中有和本发明 完全等同变量，原文中记为g ，
ij
发明中指出的几何信息和飞行时间信息分别对应原文中b 和κ ，本文
jσ
与其所指的σ对应，与之唯一不同之处是文献的方法认为所有的事件的都是用相
同的σ，而本发明根据(7)中的 的计算方法，对于会对每一个事件计算出一
个
以上内容是结合本发明的优选实施方式对所提供技术方案所作的进一步详
细说明，不能认定本发明具体实施只局限于上述这些说明，对于本发明所属技术
领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干简单
推演或替换，都应当视为属于本发明的保护范围。
发布评论