牛 善洲
赣南师范大学数学与计算机科学学院,江西 赣州 341000,
School of Mathematics and Computer Science, Gannan Normal University, Ganzhou 341000, China
赣南师范大学赣州市计算成像重点实验室,江西 赣州 341000,
Ganzhou Key Laboratory of Computational Imaging, Gannan Normal University, Ganzhou 341000, China
刘 宏
赣南师范大学数学与计算机科学学院,江西 赣州 341000,
School of Mathematics and Computer Science, Gannan Normal University, Ganzhou 341000, China
赣南师范大学赣州市计算成像重点实验室,江西 赣州 341000,
Ganzhou Key Laboratory of Computational Imaging, Gannan Normal University, Ganzhou 341000, China
刘 沛沄
赣南师范大学数学与计算机科学学院,江西 赣州 341000,
School of Mathematics and Computer Science, Gannan Normal University, Ganzhou 341000, China
赣南师范大学赣州市计算成像重点实验室,江西 赣州 341000,
Ganzhou Key Laboratory of Computational Imaging, Gannan Normal University, Ganzhou 341000, China
张 梦真
赣南师范大学数学与计算机科学学院,江西 赣州 341000,
School of Mathematics and Computer Science, Gannan Normal University, Ganzhou 341000, China
赣南师范大学赣州市计算成像重点实验室,江西 赣州 341000,
Ganzhou Key Laboratory of Computational Imaging, Gannan Normal University, Ganzhou 341000, China
李 硕
赣南师范大学数学与计算机科学学院,江西 赣州 341000,
School of Mathematics and Computer Science, Gannan Normal University, Ganzhou 341000, China
赣南师范大学赣州市计算成像重点实验室,江西 赣州 341000,
Ganzhou Key Laboratory of Computational Imaging, Gannan Normal University, Ganzhou 341000, China
梁 礼境
赣南师范大学数学与计算机科学学院,江西 赣州 341000,
School of Mathematics and Computer Science, Gannan Normal University, Ganzhou 341000, China
赣南师范大学赣州市计算成像重点实验室,江西 赣州 341000,
Ganzhou Key Laboratory of Computational Imaging, Gannan Normal University, Ganzhou 341000, China
李 楠
赣南师范大学经济管理学院,江西 赣州 341000,
School of Economics and Management, Gannan Normal University, Ganzhou 341000, China
刘 国良
赣南医学院医学信息工程学院,江西 赣州 341000,
School of Medical Information Engineering, Gannan Medical University, Ganzhou 341000, China
赣南师范大学数学与计算机科学学院,江西 赣州 341000,
School of Mathematics and Computer Science, Gannan Normal University, Ganzhou 341000, China
赣南师范大学赣州市计算成像重点实验室,江西 赣州 341000,
Ganzhou Key Laboratory of Computational Imaging, Gannan Normal University, Ganzhou 341000, China
赣南师范大学经济管理学院,江西 赣州 341000,
School of Economics and Management, Gannan Normal University, Ganzhou 341000, China
赣南医学院医学信息工程学院,江西 赣州 341000,
School of Medical Information Engineering, Gannan Medical University, Ganzhou 341000, China
其中,||∙||
F
是Frobenious范数,
β
>0是正则化参数,
R
(
X
) 是正则化函数。
1.2. 低剂量脑灌注CT图像恢复的NLSMD模型
1.2.1. NLSMD模型
首先将图像
x
k
划分为大小为
的图像块,图像块提取算子为{
R
p
x
k
}(
p
=1, 2, …,
P
),
p
是图像块索引,
P
是图像块个数。然后使用相同方法对其他图像进行相同的分块操作,最后将处于相同位置不同帧的图像块组合成矩阵
V
P
:
由于这些图像块具有相似的结构,因此其具有低秩特性。
我们可以将矩阵
V
P
分解为:
V
p
=
L
p
+
S
p
,其中,
L
p
表示低秩矩阵,
S
p
表示稀疏矩阵。矩阵
L
p
和
S
p
可以通过如下问题求出:
其中,||∙||
*
是核范数,||∙||
1
是
范数,
λ
>0正则化参数,
表示梯度算子。进一步,公式(3)可以写成如下的无约束优化问题:
其中,
μ
>0是超参数。
1.2.2. 基于NLSMD的低剂量脑灌注CT图像恢复
基于上述方法,我们提出了低剂量脑灌注CT图像恢复的NLSMD模型:
其中,
,
n
p
是矩阵
V
p
的行数,
δ
>0是超参数。我们使用交替求解算法求解式(5):
综上所述,本文算法的计算步骤如下:
初始化:迭代初始值设置为滤波反投影算法重建的低剂量脑灌注CT图像;
步骤1:根据式(6)对脑灌注CT图像
X
n
进行图像分块处理得到矩阵
V
P
n
,
p
=1, 2, …,
P
;
步骤2:根据式(10),式(13)和式(14)得到矩阵
L
p
n
+1
和矩阵
S
p
n
+1
;
步骤3:根据式(17)更新
X
n
+1
;
步骤4:重复步骤1-3,直到
或者迭代次数超过500次迭代终止。
1.3. 实验数据
根据Niu等
[
24
]
的体膜实验数据仿真得到50帧大小为256×256的低剂量脑灌注CT图像数据。相应的成像几何参数为:(1)投影角度为984,每个投影角度的探测元个数为888,每个探测器单元的大小为2 mm;(2)放射源到探测器的距离为949 mm;(3)射线源到旋转中心的距离为541 mm。X射线的入射光子数目为1×10
5
。
1.4. 对比方法
为了验证NLSMD方法的有效性,NLSMD方法分别与滤波反投影(FBP)算法,低秩与稀疏矩阵分解方法(LSMD)
[
25
]
,低秩与全变分正则化方法(LR-TV)
[
26
]
进行了比较。LSMD方法的目标函数为:
其中,
X
L
和
X
S
分别为脑灌注CT序列图像的低秩和稀疏部分,
μ
>0,
λ
>0为超参数。LR-TV方法的目标函数为:
其中,
λ
1
>0,
λ
2
>0为超参数,Φ
0
=
D
x
,Ψ
0
=
I
,Φ
1
=
D
y
,Ψ
1
=
I
,Φ
2
=
I
,Ψ
2
=
D
t
,
I
是单位矩阵,其中
D
x
,
D
y
和
D
t
分别为沿着
x
,
y
以及
t
方向的差分矩阵。
1.5. NLSMD方法的参数选取
NLSMD方法的参数有图像块
p
的大小
n
p
和参数
δ
。
为实验数据中第25帧图像随参数变化时相应的SSIM值图。从
可以看出,图像块
p
的大小
n
p
取15,正则化参数
δ
取1.6时,结构相似性(SSIM)
[
27
]
值最大为0.9946。
2. 结果
2.1. 低剂量脑灌注CT图像恢复结果
给出了低剂量脑灌注CT图像恢复结果。其中,第一行是体膜图像,第二行是FBP图像,第三行是LR-TV方法恢复的图像,第四行是LSMD方法恢复的图像,第五行是NLSMD方法恢复的图像;从左到右第一列为第5帧图像,第二列为第15帧图像,第三列为第25帧图像,第四列为第35帧图像,第五列为第45帧图像。从
可以看出NLSMD方法的图像与体膜图像最为接近。
给出了
中结果的SSIM值。从
可以看出,FBP方法的平均SSIM值为0.9438,LR-TV方法的平均SSIM值为0.9636,LSMD方法的平均SSIM值为0.9692,NLSMD方法的平均SSIM值为0.9765。
为
中各方法第25帧图像感兴趣区域的局部放大图像。从
可以看出,NLSMD方法的感兴趣区域的局部放大图像与体膜方法最为接近。为了进一步评价本文方法,我们使用了FSIM
[
28
]
值对实验结果进行定量分析。
为
的SSIM值和FSIM值。从
可以看出,FBP方法的SSIM值和FSIM值分别为0.6634和0.9298,LR-TV方法的SSIM值和FSIM值分别为0.813和0.9557,LSMD方法的SSIM值和FSIM值分别为0.8106和0.9556,NLSMD方法的SSIM值和FSIM值分别为0.8727和0.9661。
2.2. 脑血流参数图像
为各方法的CBF参数图像,其中
为体膜图像的CBF参数图像,
为FBP方法的CBF参数图像,
为LR-TV方法的CBF参数图像,
为LSMD方法的CBF参数图像,
为NLSMD方法的CBF参数图像。
为
的SSIM值和FSIM值。从
可以看出,FBP方法的SSIM值和FSIM值分别为0.7005和0.8748,LR-TV方法的SSIM值和FSIM值分别为0.7736和0.8985,LSMD方法的SSIM值和FSIM值分别为0.7557和0.8940,NLSMD方法的SSIM值和FSIM值分别为0.7871和0.9073。从
和
可以看出,NLSMD方法恢复的CBF参数图像与体膜图像的CBF参数图像最为接近。
为
感兴趣区域的局部放大图像。
为
的SSIM值和FSIM值。从
可以看出,FBP方法的SSIM值和FSIM值分别为0.6205和0.8321,LR-TV方法的SSIM值和FSIM值分别为0.7485和0.8784,LSMD方法的SSIM值和FSIM值分别为0.7952和0.8935,NLSMD方法的SSIM值和FSIM值分别为0.8540和0.9260。
为各方法的CBV参数图像,其中
为体膜图像的CBV参数图像,
为FBP方法的CBV参数图像,
为LR-TV方法的CBV参数图像,
为LSMD方法的CBV参数图像,
为NLSMD方法的CBV参数图像。
为
的SSIM值和FSIM值。从
可以看出,FBP方法的SSIM值和FSIM值分别为0.6856和0.8657,LR-TV方法的SSIM值和FSIM值分别为0.778和0.8966,LSMD方法的SSIM值和FSIM值分别为0.7642和0.8949,NLSMD方法的SSIM值和FSIM值分别为0.7972和0.9110。从
和
可以看出,NLSMD方法恢复的CBV参数图像最接近体膜的CBV参数图像。
为
感兴趣区域的局部放大图像。
为
的SSIM值和FSIM值。从
可以看出,FBP方法的SSIM值和FSIM值分别为0.5228和0.7904,LR-TV方法的SSIM值和FSIM值分别为0.6611和0.8381,LSMD方法的SSIM值和FSIM值分别为0.7176和0.8664,NLSMD方法的SSIM值和FSIM值分别为0.7800和0.8938。
3. 讨论
脑灌注CT成像是脑卒中检查的杰出代表,但脑灌注CT检查的辐射剂量较高,优化控制X射线辐射剂量已成为当前研究的热点问题。目前,低剂量脑灌注CT图像恢复方法利用高质量平扫图像和低剂量脑灌注CT图像之间的结构相似性增强低剂量脑灌注CT图像质量
[
19
-
21
]
。但是,先验图像与低剂量图像之间存在差异,会导致恢复结果产生误差。基于脑灌注CT序列图像的结构相关性,本文提出了一种基于非局部低秩稀疏矩阵分解的低剂量脑灌注CT图像恢复方法,该方法利用了低剂量脑灌注CT序列图像自身的结构冗余特性,并对相似的图像块进行非局部低秩稀疏矩阵分解,可以有效地抑制噪声且较好地保持图像细节特征信息。
在数值实验中,NLSMD方法得到的脑灌注CT图像在噪声抑制和细节保持方面有良好的表现。定量实验结果表明,NLSMD方法的平均SSIM值比FBP方法的平均SSIM值提高了3%,与LR-TV方法和LSMD方法相比均提高了1%。进一步,低剂量脑灌注CT图像的感兴趣区域放大图实验结果表明,NLSMD方法的感兴趣区域放大图的SSIM值比FBP方法提高了31.5%,比LR-TV方法提高了7%,比LSMD方法提高了8%。脑血流动力学参数图像能够精准反映脑组织的灌注情况,为临床诊断和治疗提供重要的参考价值
[
29
,
30
]
。CBF参数图像的实验结果表明,FBP方法的CBF参数图像存在大量噪声,LR-TV方法与LSMD方法的CBF参数图像的噪声虽然有所抑制,但其与体膜的CBF参数图像仍然存在差距,而NLSMD方法的CBF参数图像与体膜的CBF参数图像最为接近。CBF参数图像的定量结果进一步表明了NLSMD方法的CBF参数图像恢复效果最好。CBV参数图像实验结果表明,NLSMD方法的CBV参数图像无论从整体图像还是感兴趣区域放大图像都与体膜的CBV参数图像最为接近,并且NLSMD方法的CBV参数图像具有最高的SSIM值和FSIM值。综上所述,NLSMD方法恢复的低剂量脑灌注CT图像相比于其他方法恢复效果更好,并且NLSMD方法的CBF参数图像和CBV参数图像与体膜图像的CBF参数图像和CBV参数图像在结构和细节方面也最为接近。
但是,NLSMD方法也还存在一些局限性。首先,NLSMD方法没有考虑噪声统计特性,降低了图像块匹配的鲁棒性。其次,NLSMD方法在实际应用中所需的计算时间也是一个挑战。我们将在未来的研究中将统计噪声建模纳入NLSMD方法中增强NLSMD方法的恢复效果,并且考虑使用多核CPU和图形处理单元(GPU)硬件来加快NLSMD方法的计算速度。