[发明专利]一种基于IST的次晶体级PET系统时间修正方法

权利要求书

1.一种基于IST的次晶体级PET系统时间修正方法，包括如下步骤，IST表示迭代收缩阈值；

(1)对PET系统中的每个探测器进行次晶体级分割，每个探测器通过分割对应得到一个由n×n个晶体单元组成的晶体单元阵列，n为大于1小于m的自然数，m为探测器原晶体阵列的维度；

(2)利用次晶体级分割后的探测器对注入放射性示踪剂的生物组织进行扫描探测，得到多组LOR数据；

(3)对每一组LOR数据进行预处理，以剔除每组LOR数据中时间信息有极大偏差的Single model数据记录，并确定每组LOR数据的探测时延；

(4)将PET系统时间修正过程转化为线性方程，通过L1范数对该线性方程进行约束得到以下目标函数L；进而根据由各组LOR数据探测时延组成的探测时延序列ΔT，对目标函数L进行最小化求解得到PET系统的时间修正序列T_cal；

L＝||AT_cal-ΔT||²+λ||T_cal||₁

其中：A为系统矩阵，|| ||为L2范数，|| ||₁为L1范数，λ为预设的权重系数；

(5)最后根据所述的时间修正序列T_cal对预处理后的各组LOR数据Single model数据记录中晶体单元的探测时间进行修正。

2.根据权利要求1所述的PET系统时间修正方法，其特征在于：所述的步骤(1)中根据计算机内存空间、系统中探测器的数量及探测器所含的晶体数对探测器进行次晶体级分割，n为除去1和m以外m的任一约数。

3.根据权利要求1所述的PET系统时间修正方法，其特征在于：所述的步骤(2)中的每组LOR数据对应一对探测到同一偶合事件且分属于不同探测器内的晶体单元J_A和J_B，每组LOR数据包含多条Singlemodel数据记录，其中每条Singlemodel数据记录对应一个探测时间差即晶体单元J_A和J_B对于同一偶合事件的探测时间之差，每条Singlemodel数据记录包含晶体单元J_A和J_B的编号、对应探测时间差精度范围内晶体单元J_A和J_B对于同一偶合事件的探测计数以及对应每次探测计数晶体单元J_A和J_B对于同一偶合事件各自的探测时间。

4.根据权利要求3所述的PET系统时间修正方法，其特征在于：所述的步骤(3)中对每组LOR数据进行预处理的具体方法为：对于任一组LOR数据，将该组LOR数据中所有Single model数据记录按探测计数做成直方图，并计算所有Single model数据记录探测计数的均值，进而剔除探测计数小于均值的Single model数据记录。

5.根据权利要求1所述的PET系统时间修正方法，其特征在于：所述的步骤(3)中确定每组LOR数据探测时延的具体方法为：对于预处理后的任一组LOR数据，将其中所有Single model数据记录对应探测时间差的最大值作为该组LOR数据的探测时延。

6.根据权利要求1所述的PET系统时间修正方法，其特征在于：所述的步骤(4)中线性方程的表达式为A·T_cal＝ΔT。