第一作者:康瑞瑞(1984—),女,陕西清涧人,博士生.研究方向为图像重建.email:11118381@bjtu.edu.cn.
针对传统二维小波变换不能最优表示MR图像和MR图像数据采集缓慢的问题,基于相似性与联合稀疏变换对纵向压缩感知磁共振成像(LCS-MRI)提出新的正则化模型与相应的重建算法(frsLCS-MRI).为了验证联合稀疏变换与算法frsLCS-MRI的有效性,利用两个医学数据集与基于小波变换和相似性的算法LACS-MRI作了比较.数值实验表明:frsLCS-MRI算法可重建更高精度和更高信噪比的MR图像.使用联合稀疏变换比使用小波变换或剪切波变换可更好地重建MR图像的细节信息,如边界、拐角、轮廓、二维奇异曲线等.此外,利用相似性先验信息可节省数据采集时间.
In order to overcome the deficiencies that the traditional two-dimensional wavelet transform can not provide the optimal representation for MR images and to reduce the scanning time of MR images, a new regularization model and the corresponding reconstruction algorithm (frsLCS-MRI) are proposed based on the joint sparse transform and similarity prior. In order to verify the effectiveness of the joint sparse transform and the proposed algorithm frsLCS-MRI, we compare it with the LACS-MRI algorithm based on single sparse transform and similarity by using two medical datasets. Numerical results indicate that the proposed method frsLCS-MRI can reconstruct MR images with higher accuracy and higher signal to noise ratio. Compared with the method LACS-MRI using only wavelet transform, or shearlet transform, more detail information such as boundaries, corners, contours, can be reconstructed by the proposed method frsLCS-MRI using the joint sparse transform. Moreover, exploiting similarity based on reference image saves acquisition time.
磁共振成像(Magnetic Resonance Imaging, MRI)是一种非侵入性和非电离的成像技术, 它提供了多种对比机制, 很好地实现了解剖结构和生理功能的可视化.但是, 由于MRI中的数据是在时间线上连续地采集
最近发展起来的压缩感知(Compressed Sensing, CS)理论表明, 从高度欠采样的K-空间中可以以压倒性的概率精确重建MR图像[1], 这可以极大地减少MR图像的扫描时间和减轻病人的不适感.此外, 在许多临床MRI中, 利用已经存在的图像信息可以进一步缩短数据采集时间并提高重建图像的信噪比.这里把已经得到的采集图像称为参考图像, 参考图像与被重建的目标图像之间具有相似性[2, 3].Lior Weizman等[2]利用这种相似性与MR图像在小波变换下的稀疏性, 对纵向压缩感知磁共振成像(Longitudinal Compressed Sensing Magnetic Resonance Imaging, LCS-MRI), 提出一种自适应采样和加权重建算法LACS-MRI, 进一步地减少了MR图像扫描时间, 并提高了重建MR图像的信噪比.
虽然文献[2]利用时间域的相似性减少了采样, 但由于其利用了小波变换作为空间域的稀疏变换, 使得重建的图像依然存在类似于曲线状的伪影.这主要是因为二维小波变换是由一维小波变换张量积构成的, 小波变换可以有效地逼近具有点奇异性的一维图像特征但不能高效逼近具有细节信息等高维奇异性的图像特征.与小波变换相比, 剪切波变换[4]能够更加有效地逼近含有如边界、拐角、曲线、尖峰等几何信息的逐段光滑图像特征, 具有多尺度方法提取图像几何信息的强大能力.但是剪切波变换不能有效表示类似于点状的图像特征.很自然的想法是将两者相结合.小波变换与剪切波变换的线性组合构成新的稀疏变换(称为联合稀疏变换), 可以避免重建图像的点状伪影与曲线状的伪影, 表示更丰富的几何图像信息, 既弥补了单一方法的不足又保留了各自的优点.
本文作者利用相似性先验信息, 结合小波变换和剪切波变换, 提出一种新的正则化模型.并为了处理模型中不光滑的正则项和实现快速收敛的目的, 给出一种新的求解算法frsLCS-MRI.数值实验表明, 使用联合变换比只使用小波变换或剪切波变换可以重建更多几何信息和更高信噪比的MR图像, 而且使用相似性减少了MR图像的扫描时间.
CS是一种有效的数学框架, 只要信号是稀疏的或者测量是不相干的, 就可以使用非常少的, 远低于奈奎斯特采样定率所要求的采样重建高维信号.压缩感知磁共振成像(Compressed Sensing Magnetic Resonance Imaging, CS-MRI)重建模型为
式中:
在许多MRI的应用中, 参考图像与被重建的目标图像之间具有相似性.记参考图像为
由于
其中
尽管利用相似性先验信息可以加速采样、减少扫描时间、提高信噪比.但是, 由于有些扫描切片之间的时间间隔较长, 给LCS-MRI带来很大的挑战.这一方面是因为时间间隔较长的扫描切片之间的相似性可能是不存在的, 如有的扫描切片之间发生了大的病理改变或外科技术的改变.另一方面是因为在每个时间点都要求重建高质量的MR图像, 时间域欠采样会导致重建的MR图像不理想.Lior Weizman等[2]根据相似性的相似程度提出一种自适应采样和加权重建算法LACS-MRI.用公式表示如下:
其中
且对
式中:
为了重建MR图像的更多方向特征与高阶奇异性, 本文将小波与剪切波的线性组合作为新的稀疏变换, 并结合相似性先验提出的新模型如下:
式中:
且对
将基于相似性与联合稀疏变换提出快速LCS-MRI成像算法记为frsLCS-MRI, 简称算法1, 详细过程见表1.
注意到算法1中并不假设相似性存在, 因此在第2步加权矩阵初始化为
快速阈值迭代算法(FISTA)[6]原本是求解下面的问题:
最近 Tan等[7]基于FISTA提出一个光滑快速阈值迭代算法(SFISTA), 求解(7)的加权形式问题
其中
我们对算法SFISTA进行扩展, 提出算法rsSFISTA, 简称算法2, 用来求解最小化问题(5).为方便记述, 我们给出下面形式记法
以及
式中:
其中
rsSFISTA算法的具体过程见表2.
算法2中两个稀疏性的平衡性由
本节从不完全测量数据中重建高质量的MR图像来验证本文所提出的模型与算法.我们在两个临床纵向MRI真实数据集上进行了测验, 并与基于相似性与小波变换的算法LACS-MRI作了比较.其中真实数据与部分源码来自于文献[8].两个数据集的部分K-空间数据都是通过对一个全采样的K-空间进行降采样获得的.为了简单起见, 本文使用变密度随机采样机制和快速扫描机制.变密度随机采样机制是根据数据的概率密度函数进行采样, 概率密度函数是根据距离K-空间中心的远近来定义的, 也就是距离K-空间中心越近采样越多, 距离中心越远采样越少[2, 3].
本文使用Daubechies 4小波变换和ShearLab 3D[9, 10]的离散不可分离剪切波变换作为稀疏变换.其中离散不可分离剪切波变换使用默认设置, 也就是对0, 1, 2, 3尺度分别设有4, 4, 8, 8方向滤波的剪切滤波.总共有25个子带.在下面两个实例中, 参数
同一被诊断者的同一个切片在多个对比下都具有相似性, 特别是T2-加权对比与流体衰减反相恢复(Fluid-Attenuated Inversion Recovery, FLAIR)对比的图像在非流体区域具有高度的相似性, 利用这种相似性可以减少采样时间并提高重建MR图像的信噪比.本实例(例1)是基于一个扫描切片在不同对比度下的相似性的MRI, 具体就是基于非流体区域在T2-加权与FLAIR对比下具有高度相似性, 使用FLAIR的15%采样重建FLAIR图像.其中全采样的T2-加权对比扫面切片作为参考图像, 加权矩阵
重建的FLAIR图像如图1所示.从图1的(e)和(f)看出小波能够有效地表示点奇异性图像特征而不能够有效表示曲线奇异性几何特征.从图1的(g)和(h)看出剪切波能够更好地表示边界与曲线, 而对点的奇异性表示则差一些.从图1的(i)和(j)看出, 使用小波与剪切波的线性组合作为稀疏变换比只使用单个稀疏变换(小波变换或剪切波变换)作为稀疏变换, 可以重建更多FLAIR图像细节特征, 包括点奇异性、拐角、边界轮廓曲线等, 并且还抑制了只用小波变换引起的方块状伪影和只用剪切波变换引起的曲线状伪影.
患有脑部肿瘤的病人经常每隔几周或几个月就会进行一次脑部磁共振扫描.同一个病人在时间线上的不同扫描切片之间具有相似性, 利用这种相似性能够加速后续扫描切片的重建.本实例(例2)基于同一个被诊断者在时间(如数周或数月)线上不同扫描切片之间的相似性, 重建后续MR图像.这里把先前扫描的切片作为参考图像.对使用小波变换或剪切波变换作为稀疏变换时, 正则化参数的取值范围设为
重建的后续MR图像如图2所示.从图2看出, 使用小波变换与剪切波变换的线性组合作为稀疏变换, 比使用小波或剪切波单个变换作为稀疏变换可重建更高质量和更高精度的MR图像.只用单个稀疏变换不能重建使用联合变换重建的某些几何图像特征, 这是因为小波与剪切波都有一定局限性, 在抑制某些伪影的同时会引起另一些伪影.两者相结合可以弥补各自的不足, 使得重建的脑部图像有效地抑制了正方形状伪影和曲线状伪影, 并提高了重建图像的峰值信噪比(Peak Signal-to-Noise Ratio, PSNR)与质量.
为了定量说明本文所提出方法的性能, 我们考察了两个数据集的PSNR, 定义如下:
式中:P为峰值信噪比; M是图像中最大像素值;
表3的结果与上述两个例子的结论一致, 即使用联合稀疏变换比使用单个稀疏变换重建的MR图像具有更丰富的图像特征和更高的PSNR.此外, 与第二个例子相比, 第一个例子的PSNR值的变化范围相对比较小, 这解释了T2-加权对比与FLAIR对比之间的相似性只是存在于部分区域而不是整个图像.
本文利用MRI已有的图像与所要重建图像之间的相似性, 结合小波与剪切波提出了新的正则化模型.为了处理模型中的不可微项和达到快速收敛的目的, 对算法SFISTA进行了扩展, 提出frsLCS-MRI算法.在两个临床纵向MRI真实数据集上进行了测验, 与算法LACS-MRI作了比较, 结果表明, 使用联合变换比使用单个变换作为稀疏变换可以更好地重建图像的细节、纹理与边界等特征, 避免了由于只使用单个变换引起的伪影.另外利用相似性减少了MR图像的扫描时间.
The authors have declared that no competing interests exist.
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|