一种减少图像失真的磁共振平面回波成像方法
阅读说明：本技术 一种减少图像失真的磁共振平面回波成像方法 (Magnetic resonance plane echo imaging method for reducing image distortion ) 是由 张志勇 陈浩 于 2021-09-13 设计创作，主要内容包括：本申请提供一种减少图像失真的磁共振平面回波成像方法。使用两个射频脉冲α和β在不同时刻分别将纵向磁化矢量翻转到横向平面,得到可观测的横向磁化矢量。通过梯度脉冲在两个不同的回波链中分别重聚出通过α和β翻转纵向磁化矢量得到的横向磁化矢量。在两个梯度回波链之前分别采集导航回波信号,通过两个导航回波信号的幅值差异校正两个梯度回波链信号的幅值差异。使用低秩矩阵重建方法矫正两个梯度回波链采集的信号,得到两幅图像数据,对两幅图像数据取平均值,得到最终图像。本申请的技术方案相比于现有技术中的单次扫描方法,显著提高了成像精度,有效减少图像的畸变和失真。同时,与现有技术中的多次扫描方法相比,成像效率显著提升。(The application provides a magnetic resonance planar echo imaging method for reducing image distortion. And (3) respectively overturning the longitudinal magnetization vector to a transverse plane at different moments by using two radio frequency pulses alpha and beta to obtain an observable transverse magnetization vector. Transverse magnetization vectors, which are obtained by inverting the longitudinal magnetization vectors by α and β, are each refocused in two different echo trains by means of gradient pulses. And respectively acquiring navigation echo signals before the two gradient echo chains, and correcting the amplitude difference of the two gradient echo chain signals through the amplitude difference of the two navigation echo signals. And correcting signals acquired by the two gradient echo chains by using a low-rank matrix reconstruction method to obtain two image data, and averaging the two image data to obtain a final image. Compared with a single scanning method in the prior art, the technical scheme of the application obviously improves the imaging precision and effectively reduces the distortion and the distortion of the image. Meanwhile, compared with a multi-scanning method in the prior art, the imaging efficiency is obviously improved.)
技术领域
本申请涉及磁共振成像(MRI,Magnetic Resonance Imaging)领域，尤其涉及一种减少图像失真的磁共振平面回波成像方法。
背景技术
磁共振成像技术(Magnetic Resonance Imaging,MRI)能够无损伤地检测活体内组织清晰的解剖构造、反映器质性病变的图像，提供满足不同诊断需求的生理学信息，是现今最常用的医学检测手段之一。磁共振平面回波成像(Echo Planar Imaging,EPI)仅需要一次射频激发就能完成整个二维k空间数据的扫描采集，在几十毫秒内就能获得完整的二维图像，是目前成像速度最快的MRI技术之一。EPI突破了扫描速度对MRI技术应用的限制，进一步拓展了MRI技术的应用范围。EPI已经被广泛应用于扩散成像、灌注成像及功能成像等磁共振影像扫描中，是临床疾病诊断和科学研究的重要工具。
但EPI对磁场的不均匀性以及化学位移伪影都极为敏感，在非均匀磁场中容易出现图像畸变和信号丢失，尤其是在空气和组织交界等部位，例如大脑扫描中的前额、眼眶等部。尽管螺旋采样等非笛卡尔采样轨迹也能有效缩短MRI数据的采样时间，然而这些非笛卡尔采样方式对磁共振梯度性能的不完善性也很敏感，稳定性还需要进一步研究。另一方面，一些通过预扫描获取校正信息或者多次扫描可以校正或减少磁场不均匀引起的图像失真，但会降低扫描效率并面临多次扫描间可能的不一致问题。
在现有技术中，多采用单次扫描EPI或多次激发扫描EPI方法。
现有的单次扫描EPI通过连续梯度回波采样和相位编码blip梯度，实现k空间的连续快速编码采样，从而在一次激发扫描中就能完成整个k空间的数据填充，进而得到二维图像。然而梯度回波的采样方式无法像自旋回波那样重聚磁场不均匀产生的演化相位，并且随着空间编码的进行，相位差逐步累积，导致最终的图像出现扭曲和失真。现有技术中的单次扫描EPI序列如图1所示。
现有的多次激发扫描EPI方法将k空间分多段进行相位编码，能够缩短相位编码维的时间间隔，降低了磁场不均匀引起的相位差的累积，可以有效减少EPI图像失真并提高分辨率，但k空间被分成多段需要多次激发扫描的数据组合而成，无法在单次激发扫描中完成填充，容易受到检测部位运动的干扰而出现运动伪影。虽然依靠PROPELLER等技术可以改善多扫描EPI的运动伪影问题，但是多次激发扫描的方式使得扫描时间变长，扫描效率降低，在一些对时间分辨率要求较高的场景无法适用，这也限制了EPI图像质量的提升和进一步应用。现有技术中的k空间填充轨迹如图2所示。
因此，本领域的技术人员致力于开发一种减少图像失真的磁共振平面回波成像方法，不仅能够缩短空间编码的时间，降低磁场不均匀相位的累积，减少图像的畸变失真，而且不增加扫描次数或扫描时间，保证扫描效率。
发明内容
为实现上述目的，本申请提供了一种减少图像失真的磁共振平面回波成像方法，基于一个翻转角为α的第一射频脉冲、一个翻转角为β的第二射频脉冲、一个自旋回波模块和用于采集平面回波信号的第一梯度回波链、第二梯度回波链，其特征在于，具体包括以下步骤：
步骤一、使用所述第一射频脉冲激发层面，将第一部分磁化矢量激发到横向平面，将第二部分磁化矢量保持在纵向平面；
步骤二、经过一段延迟，使用所述第二射频脉冲激发层面，将所述第二部分磁化矢量激发到横向平面；
步骤三、施加所述自旋回波模块；
步骤四、使用梯度脉冲重聚出所述第一部分磁化矢量，在所述第一梯度回波链中检测横向的所述第一部分磁化矢量的信号，得到第一信号；
步骤五、使用梯度脉冲重聚出所述第二部分磁化矢量，在所述第二梯度回波链中检测横向的所述第二部分磁化矢量的信号，得到第二信号；
步骤六、将所述第一信号和所述第二信号沿相位编码方向交替存放，重组成k空间，进而得到最终图像。
可选地，所述第一射频脉冲的翻转角为47°；所述第二射频脉冲的翻转角为122°。
可选地，步骤二中，所述延迟的长度为所述第一梯度回波链与所述第二梯度回波链中心的时间间隔。
可选地，步骤三中，施加包含180°射频脉冲的所述自旋回波模块，以获得失真减少的T2加权图像。
可选地，步骤三中，在所述180°射频脉冲两侧施加扩散加权梯度，以获得失真减少的扩散加权图像。
可选地，步骤三中，不施加所述自旋回波模块，以获得失真减少的T2*加权图像。
可选地，在步骤四、步骤五中，选取相位编码梯度的动量为2/(γH*FoVPE)，其中γH为氢原子核的旋磁比，FoVPE为所述相位编码方向成像视野大小。
可选地，在步骤四、步骤五中，选取相位编码梯度的动量为4/(γH*FoVPE)，其中γH为氢原子核的旋磁比，FoVPE为所述相位编码方向成像视野大小。
可选地，步骤六中，在所述第一梯度回波链和所述第二梯度回波链之前分别采集导航回波信号，通过两个所述导航回波信号的幅值差异校正所述第一梯度回波链和所述第二梯度回波链信号的幅值差异，所述导航回波为未进行相位编码的梯度回波。
可选地，步骤六中，所述第一梯度回波链和所述第二梯度回波链的信号采用MUSSELS方法重建，得到两幅图像数据，对所述两幅图像数据取平均值，得到所述最终图像。
可选地，步骤四、步骤五中，所述第一梯度回波链和所述第二梯度回波链采用不同极性的相位编码梯度，并将所述第一梯度回波链和所述第二梯度回波链的信号分别重建成图像，通过两幅图像估算磁场分布，进而校正图像畸变。
可选地，还包括第三射频脉冲以及第三梯度回波链，采用所述第一梯度回波链、所述第二梯度回波链、所述第三梯度回波链采集信号，在步骤四、步骤五中，选取相位编码梯度的动量为3/(γH*FoVPE)，其中γH为氢原子核的旋磁比，FoVPE为所述相位编码方向成像视野大小。
相比于现有技术中的单次扫描方法，本申请提供的技术方案通过基于多射频激发的平面回波成像方法，显著提高了成像精度，有效减少图像的畸变和失真。同时，本申请的技术方案仍然属于单次扫描，因此与现有技术中的多次扫描方法相比，本申请的技术方案的成像效率显著提升。
以下将结合附图对本申请的构思、具体结构及产生的技术效果作进一步说明，以充分地了解本申请的目的、特征和效果。
附图说明
图1是现有技术中单次扫描EPI序列示意图；
图2是现有技术中k空间填充轨迹示意图；
图3是本申请的一个实施例的脉冲序列示意图；
图4是本申请的一个实施例的k空间填充轨迹示意图；
图5是采用本申请的一个实施例的第9层人脑图像；
图6是采用现有技术的第9层人脑图像；
图7是采用本申请的一个实施例的第15层人脑图像；
图8是采用现有技术的第15层人脑图像。
具体实施方式
以下参考说明书附图介绍本申请的多个优选实施例，使其技术内容更加清楚和便于理解。本申请可以通过许多不同形式的实施例来得以体现，本申请的保护范围并非仅限于文中提到的实施例。
本申请的总体设计思路如下：
1、使用两个射频脉冲α和β在不同时刻分别将纵向磁化矢量翻转到横向平面，得到可观测的横向磁化矢量。
2、通过梯度脉冲在两个不同的回波链中分别重聚出通过α和β翻转纵向磁化矢量得到的横向磁化矢量。
3、在两个梯度回波链之前分别采集导航回波信号，通过两个导航回波信号的幅值差异校正两个梯度回波链信号的幅值差异。
4、使用MUSSELS方法重建两个梯度回波链采集的信号，得到两幅图像数据，对两幅图像数据取平均值，得到最终图像。
其中，相位编码blip梯度的动量可以选取为2/(γH*FoVPE)或4/(γH*FoVPE)。其中γH为氢原子核的旋磁比，FoVPE为所述相位编码方向成像视野大小。从而进一步降低图像失真，得到更高分辨率图像。
其中，在第一梯度回波链和第二梯度回波链中分别使用不同极性的相位编码blip梯度，并将第一梯度回波链和第二梯度回波链的信号分别重建成图像，通过两幅图像估算磁场分布，进而校正图像畸变。
使用三个脉冲激发信号，并使用三个梯度回波链采集信号，重建得到图像。
其中，具体地，施加包含180°射频脉冲的自旋回波模块，以获得失真减少的T2加权图像；在180°射频脉冲两侧施加扩散加权梯度，以获得失真减少的扩散加权图像；不施加自旋回波模块，以获得失真减少的T2*加权图像。
具体地，本实施例是一种减少图像失真的磁共振平面回波成像方法，基于一个翻转角为α的第一射频脉冲、一个翻转角为β的第二射频脉冲、一个自旋回波模块和用于采集平面回波信号的第一梯度回波链、第二梯度回波链。其中，第一射频脉冲的翻转角为47°；第二射频脉冲的翻转角为122°具体包括以下步骤：
步骤一、使用第一射频脉冲激发层面，将第一部分磁化矢量激发到横向平面，将第二部分磁化矢量保持在纵向平面；
步骤二、经过一段延迟，使用第二射频脉冲激发层面，将第二部分磁化矢量激发到横向平面；延迟的长度为第一梯度回波链与第二梯度回波链中心的时间间隔；
步骤三、施加自旋回波模块；其中，施加包含180°射频脉冲的自旋回波模块，以获得失真减少的T2加权图像；在180°射频脉冲两侧施加扩散加权梯度，以获得失真减少的扩散加权图像；不施加自旋回波模块，以获得失真减少的T2*加权图像。
步骤四、使用梯度脉冲重聚出第一部分磁化矢量，在第一梯度回波链中检测横向的第一部分磁化矢量的信号，得到第一信号；
步骤五、使用梯度脉冲重聚出第二部分磁化矢量，在第二梯度回波链中检测横向的第二部分磁化矢量的信号，得到第二信号；
其中，在步骤四与步骤五中，相位编码梯度的动量可以选取为2/(γH*FoVPE)或4/(γH*FoVPE)，其中γH为氢原子核的旋磁比，FoVPE为相位编码方向成像视野大小。在一些实施例中，采用三个脉冲激发信号，并使用三个梯度回波链采集信号，此时相位编码梯度的动量可以选取为3/(γH*FoVPE)。
步骤六、将第一信号和第二信号沿相位编码方向交替存放，重组成k空间，进而得到最终图像。具体地，在第一梯度回波链和第二梯度回波链之前分别采集导航回波信号，通过两个导航回波信号的幅值差异校正第一梯度回波链和第二梯度回波链信号的幅值差异，导航回波为未进行相位编码的梯度回波。同时，第一梯度回波链和第二梯度回波链的信号采用MUSSELS方法重建，得到两幅图像数据，对两幅图像数据取平均值，得到最终图像。
如图3所示为本实施例提供的基于多射频激发的平面回波成像序列，该平面回波成像序列包括一个翻转角为α的第一射频脉冲、一个翻转角为β的第二射频脉冲、一个包含180°射频脉冲的自旋回波模块以及两个用于采集平面回波信号的梯度回波链(第一梯度回波链和第二梯度回波链)，RF对应射频脉冲，GRO对应读出方向所施加梯度脉冲，GPE对应相位编码方向所施加梯度脉冲，Gss对应层选方向所施加梯度脉冲。如图4所示为本实施例提供的基于多射频激发的平面回波成像序列对应的k空间填充轨迹，kx对应读出方向k空间，ky对应相位编码方向k空间，两条梯度回波链采集的信号按照图4中交替组合的方式填充到k空间中。
将现有技术的单平面回波成像序列和根据本申请的实施例的平面回波成像序列，在配有32个通道头部阵列接收线圈的3T磁共振成像系统上对一位健康志愿者进行了头部扫描的对照测量。测试中采用了两个序列：基于现有技术的平面回波成像序列和根据本实施例的平面回波成像序列。两个平面回波成像序列的对照实验中，成像视野为240×240mm2，激发层厚为5mm，TR为4000ms，TE为92ms，回波间隔为0.57ms，平面内分辨率为2.5×2.5mm2，采集层数为19。对于基于现有技术的平面回波成像序列，使用一个翻转角为90°的射频脉冲激发信号，然后施加包含180°射频脉冲的自旋回波模块，使用一条梯度回波链采集信号，梯度回波链长度为96。对于根据本实施例的平面回波成像序列，使用一个翻转角为47°的射频脉冲和一个翻转角为122°的射频脉冲激发信号，然后施加包含180°射频脉冲的自旋回波模块，使用两条梯度回波链采集信号，每条梯度回波链长度为48，相位编码梯度幅度为基于现有技术的平面回波序列的两倍。对于根据本实施例的平面回波成像序列，使用MUSSELS方法重建两个梯度回波链采集的信号，得到两幅图像数据，对两幅图像数据取平均值，得到最终图像。
图5～图8是采用了上述两种序列的志愿者实验结果。其中图5是采用根据本实施例的平面回波成像序列的第9层图像，图7是采用根据本实施例的平面回波成像序列的第15层图像，图6是采用基于现有技术的平面回波成像序列的第9层图像，图8是采用基于现有技术的平面回波成像序列的第15层图像。图5与图6相比，图5中前额部分的形状畸变大为减少(箭头所指位置)。图7与图8相比，很明显图7中眼部的形状畸变大为减少(箭头b所指位置)，信号堆积情况也显著改善(箭头a所指位置)。
以上详细描述了本申请的较佳具体实施例。应当理解，本领域的普通技术无需创造性劳动就可以根据本申请的构思作出诸多修改和变化。因此，凡本技术领域中技术人员依本申请的构思在现有技术的基础上通过逻辑分析、推理或者有限的实验可以得到的技术方案，皆应在由权利要求书所确定的保护范围内。