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

摘要

本发明公开了一种基于IST的次晶体级PET系统时间修正方法，包括(1)扫描生物组织获取single model数据；(2)对探测器单元进行次晶体级分割；(3)对数据进行预处理；(4)将L1范数约束加到线性方程式中；(5)使用IST方法求得时间修正序列；本发明方法通过将L1范数引入到时间修正线性过程中，并将时间修正水平从探测器级提升到次晶体级，来获取更精确的时间修正序列。这种方法有效地加强了修正算法对噪声的过滤，提高了PET时间修正的准确率，提升了PET系统的空间分辨率，使PET系统能在医学领域上具有更高的使用价值。

主权项

一种基于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系统的时间修正序列Tcal；L＝||ATcal‑ΔT||2+λ||Tcal||1其中：A为系统矩阵，|| ||为L2范数，|| ||1为L1范数，λ为预设的权重系数；(5)最后根据所述的时间修正序列Tcal对预处理后的各组LOR数据Single model数据记录中晶体单元的探测时间进行修正。