一种基于压缩感知技术的快速有效的动态磁共振成像方法
【专利摘要】现有的成像方法速度比较慢，硬件成本高，本发明公开了一种基于CS压缩感知技术的快速有效的动态磁共振成像方法。本发明在设定的加速因子下，产生k-t空间中不相干采样轨迹，然后按照该轨迹采集数据；接着利用动态MRI图像稀疏和低秩结构的复合信息建立目标函数；最后，通过加入Bregman迭代的kt-SLR方法求解该复合正则化目标函数，能够快速的收敛到最优解，从而重建原始的动态MRI图像序列。本发明突破了经典奈奎斯特采样定理的极限，通过随机采样少量数据点，结合联合稀疏信息，利用加速收敛的快速非线性优化算法重建动态MRI图像序列，从而大幅度的缩短了磁共振成像时间。
【专利说明】一种基于压缩感知技术的快速有效的动态磁共振成像方法
【技术领域】
[0001]本发明涉及医学成像检测【技术领域】，特别涉及医学磁共振成像【技术领域】，具体是指一种基于压缩感知技术的快速有效动态磁共振成像方法。
【背景技术】
[0002]动态磁共振成像(dMRI)是一种利用体外测量的核磁共振信号产生体内器官动态特性的体素图像成像技术，是一种具有较高软组织分辨能力的无损伤诊断方法，具有X-CT等成像方法无法比拟的优点。然而，常规的动态磁共振成像时间比较长，成本过高，而且很难避免被验者身体中生理性运动所产生影像的模糊和对比度的失真，从而限制了它的实际临床应用。
[0003]影响动态磁共振成像速度的因素主要包括两个方面:(I)原始数据采集速度；(2)k-t空间数据采集数量。研究人员通过改进MRI硬件、快速序列设计研究及有效的采集轨迹来提高原始数据采集速度。但是由于硬件和人的生理条件限制，在缩短数据采集时间方面已经接近极限。因此，更多的研究者把精力放在如何在不降低图像质量条件下减少数据采集总量上。
[0004]传统的快速动态磁共振成像方法都受到奈奎斯特定理的限制，例如k-t Blast、k-t SENSE、UNFOLD等。但是近年来提出的压缩感知理论，突破了香农(Shannon)采样定理关于采样速率必须高于两倍信号带宽的极限，因此，一些基于压缩感知理论的动态磁共振成像方法被提了出来，例如k-t SPARSE、k-t FO⑶SS、k-t SLR等。这些基于压缩感知的算法能够更好的利用动态数据k空间和时间的相关性，来减少采集成像所需的数据，不仅减少误差而且进一步加速图像采集，因此成为动态磁共振领域的研究热点。而其中最新的k-tSLR不仅利用了数据的稀疏性，还能够有效的处理低秩结构，该方法可以在不适用门控技术和存在自由呼吸的情况下对心脏等运动部位进行高分辨率扫描，成为一个非常有潜力的研究方向。
【发明内容】
[0005]本发明的目的是为了进一步加快动态磁共振成像速度，提供了一种基于压缩感知技术的快速动态磁共振成像方法。该方法继承了 k-t SLR方法的优点，不需要门控技术仍然能够对存在呼吸运动影响的心脏运动进行高分辨率动态成像，而且采样采用Bregman迭代算法求解该问题，加速了算法的收敛性，保证了重建动态图像的质量。
[0006]本发明基于压缩感知技术的快速有效动态磁共振成像方法具体步骤如下:
(I)采用笛卡尔或者非笛卡尔采样轨迹的动态磁共振成像脉冲序列，从磁共振扫描仪上获取k-t空间的测量样本数据Iii。沁与1€—￡空间信号的采样点相对应，还可以假设测量样本被噪声所干扰Ai =￡；F(x，4)exp(— JAf λ: );? +?j；i= 0,...,5-1,⑴这里，Y(It)表示空时信号，x是空间位置，t代表时间，表示第^个采样位置，:是离散的时间点总数。将上面的表达式写成向量的形式1> = Α（Gamma)+n，其中，A是傅里叶采样运算符，具体的形式由采样轨迹决定。动态磁共振成像的目标就是将信号Yktj从测量的k-t空间采样点中恢复出来。
[0007](2)采用复合优化算法进行动态磁共振成像，即将信号y(x,t)从测量的k -1空间采样点中恢复出来。下面给出算法的推导过程，并给出算法流程。
[0008]为了更方便的利用其相关性，重组空时信号为二维矩阵的形式:
【权利要求】
1.一种基于CS压缩感知技术快速有效的动态磁共振成像方法，其特征在于:在k-tSLR方法中加入了 Bregman迭代运算求解如下复合正则化问题
2.一种基于CS压缩感知技术快速有效的动态磁共振成像方法，其特征在于该方法包括以下步骤:2-1按照设定的加速因子，确定采样轨迹，采集动态磁共振成像扫描仪中的k-t空间数据；2-2对动态磁共振图像序列选择适当的稀疏变换，并结合低秩约束；2-3采用加入Bregman迭代的k-t SLR方法来求解优化目标。
3.根据权利要求1所述的一种基于CS压缩感知技术的快速有效的动态磁共振成像方法，其特征在于:步骤2-2中可以同时考虑了稀疏和低秩结构，能够更充分利用数据冗余性，从而改进了成像质量，最优化正则化参数组合由经验确定。
【文档编号】A61B5/055GK103505207SQ201210200808
【公开日】2014年1月15日 申请日期:2012年6月18日 优先权日:2012年6月18日
【发明者】董恩清, 吕成林, 李贞国, 曹祝楼 申请人:山东大学威海分校, 董恩清