文章来源:郑翊宸,陈硕,孙爱琦,等. 无需心电门控与呼吸控制的心脏磁共振实时成像.中国医学影像学杂志,2020,28(07):543-549 |
摘要 |
目的 为解决传统心脏电影成像易受屏气不完全和心律不齐等影响的缺陷,并获得高时间分辨率的心脏图像,本文实现并验证一种无需心电门控与呼吸控制的实时心脏磁共振成像和重建方法。 资料与方法 使用低秩模型与稀疏约束实现实时成像方法。前瞻性纳入6名健康志愿者及1名心律不齐患者进行心脏磁共振扫描。通过仿真优化成像重建参数;通过健康志愿者实验比较实时方法与传统电影成像方法的一致性;进行心律不齐患者实验以验证实时方法对异常心动周期的捕捉能力。 结果 模型在参数为阶数L=25、约束系数λ=0.5时表现最佳;实时方法与电影方法对于射血分数、心输出量、每搏输出量指数、左心室质量的皮尔逊相关系数分别为0.9934、0.9997、0.9936、0.9916;Bland-Altman与相关系数分析结果表明,实时方法与电影方法在评估健康志愿者的左心室面积与心功能参数时高度一致,实时方法能够捕捉心律不齐患者的发病状态。 结论 心脏实时成像方法与传统电影方法具有高度一致性,并能够实现对心律不齐患者的异常心动周期的捕捉,具有进一步研究及临床应用价值。 【关键词】心脏磁共振成像;实时成像;低秩模型;稀疏约束;磁共振成像 【中图分类号】R445.2;R541 【DOI】10.3969/j.issn.1005-5185.2020.07.016 |
心脏病是一种严重危害人类健康的常见疾病,具有高患病率、高病死率、高治疗费用等特点
[1] 。心脏磁共振电影成像( cinecardiovascular magnetic resonance , CMR )已初步成为心功能评估的“金标准”,并在心脏疾病的精准诊断中发挥重要作用
[2] 。然而由于 MRI 扫描速度较慢,且受心脏自身节律性运动及呼吸运动的影响,其位置与形态在成像期间不断变化, CMR 检查时往往需要患者屏气配合并使用心电门控,通过采集多个心动周期的数据,最终得到心脏在一个心动周期内随时间变化的运动情况。然而,该方法对于心律不齐患者存在扫描效率较低,甚至无法正常扫描等问题;同时该方法依赖于患者屏气状况,故不完全屏气会造成严重的运动伪影;此外,电影成像结果为经过平均的心动周期,无法反映不同心动周期间,尤其是心率异常周期内的心脏形态及功能。 为克服传统电影方法的问题,研究者们提出了一系列新的心脏成像方法,如并行成像、径向采集、压缩感知等,以实现心脏磁共振快速成像
[3] 。然而,这些方法存在时间分辨率不足、适用范围受限或依旧需要门控配合等问题,难以获得多心动周期、高时空分辨率的心脏图像,也难以用于心律不齐患者。 本研究拟通过基于低秩模型
[4-6] 和稀疏约束
[7] 联合并行成像的方法实现实时心脏成像,实现在不使用心电门控及屏气配合条件下的心脏磁共振采集,从而获得高时间分辨率多个心动周期内心脏动态图像,并且能够捕捉异常心动周期的心脏搏动情况。
1
资料与方法
1.1理论分析 为实现心脏实时成像,本研究使用低秩模型与稀疏约束对心脏信号进行重建恢复。通常,所采集到关于心脏的磁共振信号可表达为公式( 1 )。
![a49bc4a2e3ad2f61116a49a98386e766.png](https://img-blog.csdnimg.cn/img_convert/a49bc4a2e3ad2f61116a49a98386e766.png)
其中, 代表第
i 个接收线圈在
t 时刻采集的
k 空间数据(
i 的取值范围是
1-Nc ,
Nc 代表接收线圈的个数), 表示
t 时刻
r 位置的心脏图像数据, 表示第
i 个接收线圈的敏感度图。 表示第
i 个接收线圈在
t 时刻采集的
k 空间的测量噪声。图像重建是从采集获得的数据 求解心脏动态图像 的逆问题。 考虑离散的图像模型,每个图像序列可通过基于时间空间分布的 Casorati 矩阵表示,如公式( 2 )。
![2471383997fc94c391426452d039cb3f.png](https://img-blog.csdnimg.cn/img_convert/2471383997fc94c391426452d039cb3f.png)
其中,
M 代表时间帧的数目,
N 代表一个时间帧内对应获取的空间域中像素个数。 由于心脏动态图像中时间与空间高度相关,故 Casorati 矩阵 C 具有低秩特性,即公式( 3 )。 rank(C)≤L