CN104063887A - 一种基于Low Rank的动态PET图像重建方法 - Google Patents
一种基于Low Rank的动态PET图像重建方法 Download PDFInfo
- Publication number
- CN104063887A CN104063887A CN201410253071.2A CN201410253071A CN104063887A CN 104063887 A CN104063887 A CN 104063887A CN 201410253071 A CN201410253071 A CN 201410253071A CN 104063887 A CN104063887 A CN 104063887A
- Authority
- CN
- China
- Prior art keywords
- pet
- pet image
- matrix
- low rank
- reestablishment
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Landscapes
- Nuclear Medicine (AREA)
Abstract
本发明公开了一种基于Low Rank的动态PET图像重建方法,该方法通过建立重建问题的数学模型,等价问题模型的转化,并基于Low Rank方法重建PET图像;其中,利用Low Rank模型对PET图像进行重建的过程中,等价模型问题在求解过程中采用交替最小化算法。故本发明有效利用Low-Rank算法,改善了计算机在进行动态PET图像重建的过程中产生的结果低分辨率和噪声干扰的问题;与现有重建方法的实验比较表明,本发明能获得较好的重建效果。
Description
技术领域
本发明属于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+S Y=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浓度分布向量组合而成。
所述的最小二乘正则化模型的表达式如下:
其中:||||*为核范数,||||1为1-范数,||||F为F-范数,λ和μ均为权重系数。
所述的增强型拉格朗日函数的表达式如下:
其中:LA(L,S,γ)为关于L、S和γ的增强型拉格朗日函数,γ为拉格朗日乘子,tr()表示求迹运算;||||*为核范数,||||1为1-范数,||||F为F-范数,λ和μ均为权重系数,β为惩罚系数。
所述的步骤(6)中通过以下迭代方程组对增强型拉格朗日函数进行最小化求解:
Lk+1=Dτ/β(Lk-τD(D(Lk+Sk)-bk)
Sk+1=Sτ/β(Sk-τD(D(Lk+Sk)-bk)
γk+1=γk-β(D(Lk+1+Sk+1)-Y-Rk+1)
其中:Lk和Lk+1分别为第k次和第k+1次迭代的背景,Sk和Sk+1分别为第k次和第k+1次迭代的生物组织,γk和γk+1分别为第k次和第k+1次迭代的拉格朗日乘子,Dτ/β()表示奇异值阈值算子,Sτ/β()表示软收缩算子,μ为权重系数,β为惩罚系数,bk为第k次迭代的中间矩阵,τ为接近变量且ρ(D*D)表示矩阵D*D的谱半径。
所述的中间矩阵bk的表达式如下:
bk=Y+Rk+1+γk/β
通过所述的迭代方程组进行迭代计算,则使达到最大迭代次数或迭代收敛后的背景与生物组织相叠加即作为估计得到的PET浓度分布矩阵X;迭代收敛条件如下:
其中:||||F为F-范数。
本发明通过建立重建问题的数学模型,等价问题模型的转化,并基于Low Rank方法重建PET图像;其中,利用Low Rank模型对PET图像进行重建的过程中,等价模型问题在求解过程中采用交替最小化算法。故本发明有效利用Low-Rank算法,改善了计算机在进行动态PET图像重建的过程中产生的结果低分辨率和噪声干扰的问题;与现有重建方法的实验比较表明,本发明能获得较好的重建效果。
附图说明
图1为本发明PET图像重建方法的流程示意图。
图2(a)为关于肺部体模的真值图像。
图2(b)为采用EM算法重建肺部体模的PET图像。
图2(c)为采用本发明方法重建肺部体模的PET图像。
图2(d)为采用本发明方法分割出的肺部组织的PET图像。
具体实施方式
为了更为具体地描述本发明,下面结合附图及具体实施方式对本发明的技术方案进行详细说明。
如图1所示,一种基于Low Rank的动态PET图像重建方法,包括步骤:
正电子发射断层扫描仪探测人体内发出的放射性信号,经过符合和采集系统处理,形成原始投影线,并以正弦图的方式存放于计算机硬盘中;对原始采集到的sinogram,以sinogram和已知的系统矩阵D为输入项,调用相关模块。
S1.根据PET探测的原理建立重建问题的基本模型;
S2.引入Low-Rank来优化问题模型;
S3.初始化,设置初始点γ0,L0和S0;设置β、λ、μ,
S4.设置的初始点开始,按照交替最小化方法计算增强版本的拉格朗日方程中的噪声rk+1和重建的图像的Low-Rank部分Lk+1以及Sparse部分Sk+1;
S5.更新γk+1;准备将更新后的参数返回步骤S4进行循环;
S6.判断是否满足迭代停止条件,不满足该条件则执行步骤S5,满足则迭代停止;进而使L+S=X得到的PET浓度分布向量X实现PET成像。
为了完成PET图像的重建,PET检测过程的基本模型基于如下方程:
y=Dx+e
其中,y为符合计数向量,这里测量所得的数据一般为sinogram,x为PET浓度分布向量,e为噪声向量,D为系统矩阵(系统矩阵表征了发射光子被探测器接收的概率,其受探测器结构、探测效率、衰减、死时间等因素的决定)。重建的过程即在采集到y后,通过已知的稀疏的D来计算x。
由于D的稀疏性,不能够直接通过求逆来完成重建,因此考虑转化问题模型。故进一步地,通过Low Rank模型来解决上述问题,模型如下:
X=L+S Y=D(L+S)+E
其中,X为PET浓度分布矩阵(由与符合计数向量y对应的多组PET浓度分布向量x组合而成),L为低秩的序列共享的背景,S为分割出的活动的生物组织。重建的过程即在采集到Y(由多组符合计数向量组合成的符合计数矩阵)后,通过已知的稀疏的D来计算L和S,进一步求得X。
进一步地对Low Rank模型进行最小二乘正则化,将问题优化为:
其中,λ为权重系数,||L||*表示核范数,表示L的奇异值之和;||S||1表示S的L1范数,||S||1=Σi,j|Si,j|;μ为随噪声等级调整的权重参数;||||F是弗罗贝尼乌斯范数。
进一步地对上式进行求解,得到对应的增强型拉格朗日函数为:
其中:LA(L,S,γ)为关于L、S和γ的增强型拉格朗日函数,γ为拉格朗日乘子,tr() 表示求迹运算,β为惩罚系数。
通过以下迭代方程组对上述拉格朗日目标函数进行最小化求解:
γk+1=γk-β(D(Lk+1+Sk+1)-Y-Rk+1)
其中:Rk+1的计算表达式如下:
Lk+1的计算表达式如下:
Lk+1=Dτ/β(Lk-τD(D(Lk+Sk)-bk)
Sk+1的计算表达式如下:
Sk+1=Sτ/β(Sk-τD(D(Lk+Sk)-bk)
bk=Y+Rk+1+γk/β
τ为接近变量,ρ(D*D)为矩阵D*D的谱半径;Dτ/β()表示奇异值阈值算子,Sτ/β()表示软收缩算子。
通过上述迭代方程组进行迭代计算,则当达到最大迭代次数或迭代收敛后的背景与生物组织相叠加即作为估计得到的PET浓度分布矩阵X;迭代收敛条件如下:
以下我们采用脑部数字体模模型实验来验证本实施方式的有效性,该模型包含一些高浓度区域。本实验运行环境为:8G内存,3.40GHz,64位操作系统, CPU为intel酷睿双核。
将本实施方式基于Low-Rank的PET重建方法和传统的EM(最大期望算法)方法重建结果做比较,二者使用相同的观测值Y和相同的系统矩阵D以保证结果的可比性,具体参数设置如下:Y为n×n维采集到的sinogram,令m=n×n,D为m×m维事先计算好的系统矩阵;这里n=64,即m=4096。
针对重建图像质量的验证,采用高维度的原始采集数据64个投影角度,每个角度下射束为64条,即m=4096,重建图像大小为64×64,即维度n=4096;初始值设定同上。图2(a)~(c)是真值图像、传统EM方法重建的图像和基于本实施方式重建的图像的比对示意图,可以直观地看出基于本实施方式重建的图像与EM的结果相比,能够恢复出更多的结构,且能够在重建的同时完成图像分割,获得活动组织的图像;图2(d)是分割出来的活动的生物组织图像。
对于相同的数据,分别采用本实施方式和传统的EM方法进行比较,如表1所示;应用本实施方式重建结果在与真值的偏差,方差和均方根误差均小于传统的EM方法,说明对本发明技术方案在提高精确度和降低噪声方面的可行性。
表1
方法 | 偏差 | 方差 | 均方根误差 |
EM | 0.0599 | 0.0136 | 0.1168 |
本实施方法 | 0.0361 | 0.0061 | 0.0781 |
Claims (8)
1.一种基于Low Rank的动态PET图像重建方法,包括如下步骤:
(1)利用探测器对注入有放射性物质的生物组织进行探测,动态采集得到PET的n组符合计数向量,并对所述的符合计数向量进行校正,进而构建PET的符合计数矩阵;n为大于1的自然数;
(2)根据PET成像原理,建立PET的测量方程;
(3)通过对所述的测量方程引入Low Rank,得到PET的Low Rank模型如下:
X=L+S Y=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图像。
2.根据权利要求1所述的动态PET图像重建方法,其特征在于:所述的测量方程的表达式如下:
y=Dx+e
其中:y为校正后的符合计数向量,x为PET浓度分布向量,e为噪声向量。
3.根据权利要求2所述的动态PET图像重建方法,其特征在于:所述的PET浓度分布矩阵X由与符合计数向量对应的n组PET浓度分布向量组合而成。
4.根据权利要求1所述的动态PET图像重建方法,其特征在于:所述的最小二乘正则化模型的表达式如下:
其中:||||*为核范数,||||1为1-范数,||||F为F-范数,λ和μ均为权重系数。
5.根据权利要求1所述的动态PET图像重建方法,其特征在于:所述的增强型拉格朗日函数的表达式如下:
其中:LA(L,S,γ)为关于L、S和γ的增强型拉格朗日函数,γ为拉格朗日乘子,tr()表示求迹运算;||||*为核范数,||||1为1-范数,||||F为F-范数,λ和μ均为权重系数,β为惩罚系数。
6.根据权利要求1所述的动态PET图像重建方法,其特征在于:所述的步骤(6)中通过以下迭代方程组对增强型拉格朗日函数进行最小化求解:
Lk+1=Dτ/β(Lk-τD(D(Lk+Sk)-bk)
Sk+1=Sτ/β(Sk-τD(D(Lk+Sk)-bk)
γk+1=γk-β(D(Lk+1+Sk+1)-Y-Rk+1)
其中:Lk和Lk+1分别为第k次和第k+1次迭代的背景,Sk和Sk+1分别为第k次和第k+1次迭代的生物组织,γk和γk+1分别为第k次和第k+1次迭代的拉格朗日乘子,Dτ/β()表示奇异值阈值算子,Sτ/β()表示软收缩算子,μ为权重系数,β为惩罚系数,bk为第k次迭代的中间矩阵,τ为接近变量且ρ(D*D)表示矩阵D*D的谱半径。
7.根据权利要求6所述的动态PET图像重建方法,其特征在于:所述的中间矩阵bk的表达式如下:
bk=Y+Rk+1+γk/β。
8.根据权利要求6所述的动态PET图像重建方法,其特征在于:通过所述的迭代方程组进行迭代计算,则使达到最大迭代次数或迭代收敛后的背景与生物组织相叠加即作为估计得到的PET浓度分布矩阵X;迭代收敛条件如下:
其中:||||F为F-范数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410253071.2A CN104063887A (zh) | 2014-06-09 | 2014-06-09 | 一种基于Low Rank的动态PET图像重建方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410253071.2A CN104063887A (zh) | 2014-06-09 | 2014-06-09 | 一种基于Low Rank的动态PET图像重建方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104063887A true CN104063887A (zh) | 2014-09-24 |
Family
ID=51551579
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410253071.2A Pending CN104063887A (zh) | 2014-06-09 | 2014-06-09 | 一种基于Low Rank的动态PET图像重建方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104063887A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104657950A (zh) * | 2015-02-16 | 2015-05-27 | 浙江大学 | 一种基于Poisson TV的动态PET图像重建方法 |
CN107991636A (zh) * | 2017-11-24 | 2018-05-04 | 哈尔滨工业大学 | 一种基于适应性结构低秩矩阵的快速磁共振图像重建方法 |
CN108324263A (zh) * | 2018-01-11 | 2018-07-27 | 浙江大学 | 一种基于低秩稀疏约束的无创心脏电生理反演方法 |
CN116503671A (zh) * | 2023-06-25 | 2023-07-28 | 电子科技大学 | 基于有效秩张量近似的残差网络压缩的图像分类方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2007026234A1 (en) * | 2005-08-31 | 2007-03-08 | Ge Healthcare Limited | Method and system of multivariate analysis on normalized volume-wise data in the sinogram domain for improved quality in positron emission tomography studies |
CN101156780A (zh) * | 2007-09-04 | 2008-04-09 | 陈武凡 | 在pet成像中最大后验优化图像重建方法 |
CN103810731A (zh) * | 2014-01-20 | 2014-05-21 | 浙江大学 | 一种基于tv范数的pet图像重建方法 |
-
2014
- 2014-06-09 CN CN201410253071.2A patent/CN104063887A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2007026234A1 (en) * | 2005-08-31 | 2007-03-08 | Ge Healthcare Limited | Method and system of multivariate analysis on normalized volume-wise data in the sinogram domain for improved quality in positron emission tomography studies |
CN101156780A (zh) * | 2007-09-04 | 2008-04-09 | 陈武凡 | 在pet成像中最大后验优化图像重建方法 |
CN103810731A (zh) * | 2014-01-20 | 2014-05-21 | 浙江大学 | 一种基于tv范数的pet图像重建方法 |
Non-Patent Citations (2)
Title |
---|
CHEN S, CHEN Z, HU Z, ET AL: "Simultaneous reconstruction and segmentation for dynamic PET: A low rank framework", 《2014 IEEE 11TH INTERNATIONAL SYMPOSIUM ON BIOMEDICAL IMAGING (ISBI)》 * |
鲍德: "《凸优化》", 31 January 2013 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104657950A (zh) * | 2015-02-16 | 2015-05-27 | 浙江大学 | 一种基于Poisson TV的动态PET图像重建方法 |
CN104657950B (zh) * | 2015-02-16 | 2017-05-03 | 浙江大学 | 一种基于Poisson TV的动态PET图像重建方法 |
CN107991636A (zh) * | 2017-11-24 | 2018-05-04 | 哈尔滨工业大学 | 一种基于适应性结构低秩矩阵的快速磁共振图像重建方法 |
CN107991636B (zh) * | 2017-11-24 | 2020-05-12 | 哈尔滨工业大学 | 一种基于适应性结构低秩矩阵的快速磁共振图像重建方法 |
CN108324263A (zh) * | 2018-01-11 | 2018-07-27 | 浙江大学 | 一种基于低秩稀疏约束的无创心脏电生理反演方法 |
CN116503671A (zh) * | 2023-06-25 | 2023-07-28 | 电子科技大学 | 基于有效秩张量近似的残差网络压缩的图像分类方法 |
CN116503671B (zh) * | 2023-06-25 | 2023-08-29 | 电子科技大学 | 基于有效秩张量近似的残差网络压缩的图像分类方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11501473B2 (en) | Systems and methods for image correction in positron emission tomography | |
CN104657950B (zh) | 一种基于Poisson TV的动态PET图像重建方法 | |
US10765382B2 (en) | Method for mixed tracers dynamic PET concentration image reconstruction based on stacked autoencoder | |
US9155514B2 (en) | Reconstruction with partially known attenuation information in time of flight positron emission tomography | |
CN109009179B (zh) | 基于深度置信网络的相同同位素标记双示踪剂pet分离方法 | |
CN106491151B (zh) | Pet图像获取方法及系统 | |
CN105894550B (zh) | 一种基于tv和稀疏约束的动态pet图像和示踪动力学参数同步重建方法 | |
US20230059132A1 (en) | System and method for deep learning for inverse problems without training data | |
CN105147312A (zh) | Pet图像获取方法及系统 | |
CN108550172B (zh) | 一种基于非局部特性和全变分联合约束的pet图像重建方法 | |
CN103810731A (zh) | 一种基于tv范数的pet图像重建方法 | |
CN102831627A (zh) | 一种基于gpu多核并行处理的pet图像重建方法 | |
CN106204674A (zh) | 基于结构字典和动力学参数字典联合稀疏约束的动态pet图像重建方法 | |
CN102509322B (zh) | 一种基于卡尔曼滤波的pet图像重建方法 | |
CN103279964A (zh) | 一种基于prca的pet图像动态重建方法及系统 | |
CN104063887A (zh) | 一种基于Low Rank的动态PET图像重建方法 | |
CN104851080A (zh) | 一种基于tv的三维pet图像重建方法 | |
CN103230282A (zh) | 一种pet浓度均值与方差的估计方法及系统 | |
Galve et al. | Super-iterative image reconstruction in PET | |
CN103793929A (zh) | 一种基于h2/h∞混合滤波的静态pet图像重建方法 | |
CN103263275A (zh) | 一种基于h无穷滤波的pet生理参数的重构方法 | |
CN105279777A (zh) | 基于改进序贯滤波的静态pet图像重建方法 | |
US20230072175A1 (en) | Attenuation distribution image creating device, image processing device, radiation computed tomography system, attenuation distribution image creating method, and image processing method | |
CN112927132B (zh) | 提升pet系统空间分辨率均匀性的pet图像重建方法 | |
Shen et al. | Limited view PET reconstruction of tissue radioactivity maps |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20140924 |
|
RJ01 | Rejection of invention patent application after publication |