[发明专利]一种闪烁脉冲事件的符合处理方法

权利要求书

1.一种闪烁脉冲事件的符合处理方法，其特征在于，所述符合处理方法包括以下步骤：

步骤S1：获取PET系统的实验数据，其中，所述实验数据包括正电子湮灭后产生的一对能量相同方向相反的伽马光子到达对应的两个探测器的时间信息以及所述伽马光子在所述探测器中沉积的能量信息，所述时间信息分别为T_i1和T_i2，所述能量信息分别为E_i1和E_i2；

步骤S2：设定所述PET系统的能量窗和时间窗并获取符合事件，其中，所述能量窗设定为E，所述时间窗设定为T，满足以下条件的所述伽马光子记为一对真符合事件：

∣T_i1－T_i2∣≤T，E_i1≤E，且E_i2≤E；

步骤S3：对所述真符合事件进行加权处理，其中，加强处理采用以下两种公式：

第一种：

第二种：

其中，p_i表示第i条响应线上的投影值，i为自然数，σ和a分别为与PET系统、成像假体以及活度相关的可变参数；

步骤S4：进行图像重建，获取图像结果，所述图像重建的方法为最大似然-期望最大算法，采用的公式为：

其中，表示第j个像素经过n+1次迭代后的值，j为自然数，a_ij表示第i条响应线上第j个像素所占的比重；

步骤S5：对所述图像结果进行分析。

2.根据权利要求1所述的闪烁脉冲事件的符合处理方法，其特征在于，在所述步骤S1中，所述实验数据还包括伽马光子的位置信息，所述位置信息通过对所述PET系统的闪烁晶体条进行编号获得。

3.根据权利要求1所述的闪烁脉冲事件的符合处理方法，其特征在于，在所述步骤S2中，能量窗T的设定范围为150keV-650KeV，能量分辨率设定为20％。

4.根据权利要求1所述的闪烁脉冲事件的符合处理方法，其特征在于，在所述步骤S3中，σ的取值范围为0.1-1000。

5.根据权利要求4所述的闪烁脉冲事件的符合处理方法，其特征在于，σ的取值范围为1.5-100。

6.根据权利要求1所述的闪烁脉冲事件的符合处理方法，其特征在于，在所述步骤S3中，a的取值范围为0.01-100。

7.根据权利要求6所述的闪烁脉冲事件的符合处理方法，其特征在于，a的取值范围为1-50。

8.根据权利要求1所述的闪烁脉冲事件的符合处理方法，其特征在于，σ的取值为2，a的取值为0.5。

9.根据权利要求1所述的闪烁脉冲事件的符合处理方法，其特征在于，在所述步骤S5中，对所述图像结果进行分析采用的对比参数为恢复度对比系数，所述恢复度对比系数包括热区恢复度对比系数以及冷区恢复度对比系数，其中，所述热区恢复度对比系数为：

其中，C_H,j表示重建图像中热区的平均计数，C_B,j表示重建图像中背景的平均计数，a_H表示仿真时热区的活度值，a_B表示仿真时背景的活度值；

冷区恢复度对比系数为：

其中，C_C,j表示重建图像中冷区的平均计数，C_B,j表示重建图像中背景的平均计数。