[发明专利]一种基于Low Rank的动态PET图像重建方法

说明书

技术领域

本发明属于PET成像技术领域，具体涉及一种基于Low Rank(低秩)的动态PET图像重建方法。

背景技术

正电子发射断层成像(Positron emission tomography，PET)是一种基于核物理学、分子生物学的医学影像技术，它能够从分子水平上观察细胞的代谢活动，为早期疾病的检测和预防提供了有效依据。在进行PET扫描时首先需要将由放射性同位核素标记的药物注入人体内，通过血液循环系统，这些物质在人体内各组织器官中将形成一定的分布。由于放射性同位核素不稳定，将发生衰变，衰变过程中所产生的正电子与组织中的电子发生湮灭反应，产生一对方向几乎相反的伽玛光子，经由符合采集系统对这些带有放射性药物分布信息的成对光子进行处理生成投影数据。通过相应的数学方法对投影数据进行反演求解，可重建出放射性物质的空间浓度分布。

动态正电子发射断层扫描测量放射性在生物体组织的空间分布，通过连续数据采集，生理状态成像组织可以从动态PET图像序列来获得，这使得动态PET在临床诊断的强大工具和生物医学研究。活动地图进行重构从正弦图，通过PET扫描仪收集到的原始数据。随后，图像通常被分割为多个ROI(感兴趣区域)代谢参数估计，并进一步病理分析。

传统上，放射性浓度分布重建和感兴趣区域边界分割问题被视为两个顺序步骤。放射性浓度分布重建往往采用统计迭代方法，其中包括著名的MLEM(最大似然期望最大化)、MAP(最大后验)和SAGE(惩罚似然)算法，然后再应用以聚类，形变模型或图割(Graph Cut)为基础的方法来分割PET图像。然而，由于测量数据的复杂性，所得到的放射性浓度分布存在低分辨率和噪声干扰的问题，造成后续分割的困难。

发明内容

针对现有技术所存在的上述技术问题，本发明提供了一种基于Low Rank的动态PET图像重建方法，解决了计算机在进行图像重建的过程中产生的结果低分辨率和噪声干扰的问题。

一种基于Low Rank的动态PET图像重建方法，包括如下步骤：

(1)利用探测器对注入有放射性物质的生物组织进行探测，动态采集得到PET的n组符合计数向量，并对所述的符合计数向量进行校正，进而构建PET的符合计数矩阵；n为大于1的自然数；

(2)根据PET成像原理，建立PET的测量方程；

(3)通过对所述的测量方程引入Low Rank，得到PET的Low Rank模型如下：

X＝L+SY＝D(L+S)+E

其中：D为系统矩阵，Y为符合计数矩阵，X为PET浓度分布矩阵，E为噪声矩阵，L为低秩的序列共享的背景，S为分割出的活动的生物组织；

(4)对所述的Low Rank模型进行最小二乘正则化，得到对应的最小二乘正则化模型；

(5)对所述的最小二乘正则化模型进行求解，得到对应的增强型拉格朗日函数；

(6)对所述的增强型拉格朗日函数进行最小化求解，同时获得背景L和生物组织S，通过使背景L与生物组织S相叠加得到PET浓度分布矩阵X，并对其进行PET成像从而获得n帧连续的PET图像。

所述的测量方程的表达式如下：

y＝Dx+e

其中：y为校正后的符合计数向量，x为PET浓度分布向量，e为噪声向量。

所述的PET浓度分布矩阵X由与符合计数向量对应的n组PET浓度分布向量组合而成。

所述的最小二乘正则化模型的表达式如下：