1. 东北大学 医学与生物信息工程学院，辽宁 沈阳 110169;
2. 东北大学 计算机科学与工程学院，辽宁 沈阳 110169;
3. 东北大学 医学影像智能计算教育部重点实验室，辽宁 沈阳 110169;
4. 深圳技术大学 健康与环境工程学院，广东 深圳 518118
收稿日期：2020-04-07
基金项目：国家重点研发计划项目(2018YFC1314501)。
作者简介：王艳飞(1977-), 男, 内蒙古赤峰人, 东北大学博士研究生;
杨金柱(1977-)，男，内蒙古通辽人，东北大学教授，博士生导师;
康雁(1964-), 男, 辽宁沈阳人, 东北大学教授，深圳技术大学****，博士生导师。
摘要：研究了快速自旋回波产生受激回波的原理，以及受激回波与射频脉冲相位角之间的关系.针对快速自旋回波脉冲序列提出了一种受激回波伪影校正方法.通过两次激发分离出自旋回波与受激回波，然后调整预散相梯度的面积使受激回波与自旋回波中心重合，并修改180°回聚脉冲相位角使受激回波与自旋回波同相位.在1.5 T超导磁共振成像系统对校正方法进行了验证，经过受激回波伪影校正后，快速自旋回波图像的信噪比提高了57.08%，图像伪影降低了76.12%.实验结果表明，所提出的方法可以有效抑制受激回波对快速自旋回波成像的影响，提高了图像质量.
关键词：磁共振成像快速自旋回波受激回波射频脉冲伪影校正
A Method to Correct Stimulated Echo Artifact for Fast Spin Echo
WANG Yan-fei1, TAN Wen-jun2, YANG Jin-zhu2,3, KANG Yan1,4
1. College of Medicine and Biomedical Information Engineering, Northeastern University, Shenyang 110169, China;
2. School of Computer Science & Engineering, Northeastern University, Shenyang 110169, China;
3. Key Laboratory of Intelligent Computing in Medical Image, Ministry of Education, Northeastern University, Shenyang 110169, China;
4. College of Health Science and Environmental Engineering, Shenzhen Technology University, Shenzhen 518118, China
Corresponding author: KANG Yan, E-mail: kangyan@bmie.neu.edu.cn.
Abstract: The principle of fast spin echo (FSE) producing stimulated echo and the relationship between stimulated echo and the phase angle of radio frequency pulse were investigated. For fast spin echo pulse sequence, a stimulated echo artifact correction method was proposed. The spin echo and stimulated echo were separated by two times of excitations, then the area of dispersion gradient was adjusted to make the center of stimulated echo coincide with that of the spin echo, and the phase angle of 180° refocusing pulse was modified to make the phase of stimulated echo and spin echo consistent. The correction method was tested and validated on a 1.5 T superconducting MRI scanner. After stimulated echo artifact was corrected, the signal to noise ratio of the FSE image was improved by 57.08%, while the image artifact was reduced by 76.12%. The experimental results showed that the proposed method can effectively suppress stimulated echo artifacts on FSE image and improve the image quality.
Key words: MRI(magnetic resonance imaging)FSE(fast spin echo)stimulated echoradio frequency pulseartifact correction
快速自旋回波(fast spin echo，FSE)[1]，在90°射频脉冲激励后，通过连续施加多个180°回聚脉冲，重聚散相的横向磁化矢量.每个180°回聚脉冲配合不同的相位编码梯度，这样在一次重复时间可以完成多个回波数据采集，提高了成像速度.目前FSE已经成为最常用的临床诊断脉冲序列之一，广泛用于各种解剖部位的磁共振成像中[2-4].
快速自旋回波在一次重复时间通过施加多个180°回聚脉冲来降低成像时间，但由于射频场不均匀、射频功率校正不准确等因素的影响，180°回聚脉冲的翻转角并不能达到理想的180°，导致接收到的回波信号不但包含了自旋回波，还包含了受激回波[5-7].在编码梯度的作用下，受激回波和自旋回波的空间编码信息不一致，最终导致图像产生伪影.
常用的受激回波伪影校正方法包括：①使用粉碎梯度(crusher gradients)[8-10]；②通过预扫描进行校正.粉碎梯度方法，通常修改180°回聚脉冲所对应的选层梯度，在选层梯度的两侧各加一个极性相同的梯度，通过这两个梯度对横向磁化矢量进行散相和聚相，改变受激回波的相位，抑制受激回波的产生.粉碎梯度方法并不能完全消除受激回波对快速自旋回波成像的影响.预扫描方法，是在脉冲序列正式扫描前，通过采集额外的校正数据分解出受激回波，并对受激回波进行抑制.一种预扫描方法是通过修改快速自旋回波脉冲序列，在每一个读梯度的两侧分别加上一个散相梯度，校正由受激回波引起的相位差[11]；另一种预扫描方法[12-13]，通过多次改变粉碎梯度的强度值以及90°射频脉冲和180°回聚脉冲的相位角，使受激回波或者自旋回波消失，实现对受激回波的单独调节.
本文提出一种基于预扫描法的受激回波伪影校正方法，该方法简单、快速，且不需要更改脉冲序列.校正数据采集过程中，只需要更改一次180°回聚脉冲的相位角，就可以分离出受激回波和自旋回波.根据分离出来的受激回波和自旋回波的回波中心位置差，校正读梯度方向的预散相梯度的面积，根据受激回波和自旋回波的相位差修正180°回聚脉冲的相位角，即可完成受激回波伪影校正.
1 算法与模型1.1 受激回波形成与特点当180°回聚脉冲的翻转角不是理想180°时，回聚脉冲含有90°射频脉冲的等效分量.受激回波可以普遍表示成由3个连续的90°射频脉冲形成[14-15].为了理解受激回波的复杂形成原理，又不失一般性，本文把横向磁化矢量抽象成6个分量，其中a1，a2，a3逆时针散相，b1，b2，b3顺时针散相.经过时间τ，磁化矢量a1，a2，a3，b1，b2，b3的散相角度分别为θ1，θ2，θ3，Φ1，Φ2，Φ3.其中θ1，Φ1∈(0, π/2)，θ2，Φ2=π/2，θ3，Φ3∈(π/2, π).射频脉冲的持续时间非常短，可以认为瞬时完成.
如图 1所示，设z轴与主磁场方向平行，x′，y′轴是与z轴垂直的逻辑坐标系，所施加的3个90°射频脉冲的相位角全为0°.图 1a为t=0时刻沿x′轴方向施加第一个90°射频脉冲，纵向磁化矢量翻转90°落在+y′轴上，形成横向磁化矢量.受主磁场不均匀及T2弛豫物理特性影响，横向磁化矢量开始散相，经过时间τ后，磁化矢量散相如图 1b所示.
图 1(Fig. 1)
 | 图 1 射频脉冲相位角为0°，受激回波形成原理图Fig.1 Schematic diagram of stimulated echo formation with the phase angle of radio frequency pulse equal 0° (a)—90°射频脉冲激发；(b)—横向磁化矢量散相；(c)—横向磁化矢量翻转；(d)—磁化矢量在纵向平面投影；(e)—翻转纵向磁化矢量；(f)—受激回波形成. |
在t=τ时刻，沿x′轴方向施加第二个90°射频脉冲.受射频脉冲激发，所有横向磁化矢量翻转到x′z平面，如图 1c所示.其中，a1和b1位于-z平面，a2和b2位于x′轴，a3和b3位于+z平面.此后，所有磁化矢量在x′y′平面方向，继续以恒定的速度绕z轴散相，在t=2τ时刻，所有的磁化矢量的水平分量全部转到-y′平面，产生自旋回波(也称为哈恩回波).
在第二个90°射频脉冲施加完成后，仅考虑图 1c中磁化矢量在纵向平面的变化情况.此时a2和b2位于x′轴，没有纵向分量，a3和b3的纵向分量位于+z轴，a1和b1的纵向分量位于-z轴，如图 1d所示.此后，所有纵向磁化矢量以T1弛豫时间沿z轴方向恢复，磁化矢量在纵向的大小发生变化，但方向并没有变化.
如图 1e所示，在t=τ+T时刻，沿x′轴方向施加第三个90°射频脉冲.其中，T为混合时间.此时a3和b3的纵向分量翻转到+y′轴，a1和b1的纵向分量翻转到-y′轴.a2和b2不受射频脉冲的影响，方向不发生变化.此后，所有磁化矢量继续以恒定的速度在x′y′平面绕z轴散相.如图 1f所示，经过τ时间后，a1逆时针转θ1角位于第三象限，b3顺时针旋转Φ3角位于第三象限，a2和b2重合位于-y′轴，b1和a3旋转到第二象限.此时，所有横向磁化矢量全部位于-y′平面，在-y′轴产生净磁化矢量，形成受激回波.
上述分析中所有的90°射频脉冲相位角全等于0°，接下来的分析，改变第二个和第三个90°射频脉冲的相位角为90°.如图 2a所示，在t=0时刻施加第一个翻转角等于90°、相位角等于0°的射频脉冲.如图 2b所示，在t=τ时刻，散相的横向磁化矢量受第二个90°射频脉冲作用，由于第二个射频脉冲相位角等于90°，这时横向磁化矢量会绕y′轴旋转，全部横向磁化矢量都在y′z平面，见图 2c.此时，仅考虑磁化矢量在纵向的变化情况，从y′轴观察，a1，a2，a3位于-z轴，b1，b2，b3位于+z轴，如图 2d所示.此后，纵向磁化矢量沿z向以T1弛豫时间恢复.在t=τ+T时刻，施加相位角、翻转角均为90°的第三个射频脉冲，见图 2e.之后各横向磁化矢量以恒定的速度绕z轴旋转，经过τ时间后所有的纵向磁化矢量全部位于-y′轴，形成受激回波，见图 2f.
图 2(Fig. 2)
 | 图 2 第二和第三个射频脉冲相位角为90°，受激回波形成原理图Fig.2 Schematic diagram of stimulated echo formation with phase angles of the second and the third radio frequency pulse equal 90° (a)—90°射频脉冲激发；(b)—横向磁化矢量散相；(c)—横向磁化矢量翻转；(d)—沿y′轴观察磁化矢量；(e)—翻转磁化矢量；(f)—受激回波形成. |
由此可以得出：①连续3个翻转角为90°的射频脉冲会产生受激回波.实际上，只要射频脉冲的翻转角不等于180°，就会有90°翻转角等效分量，就会形成受激回波；② 3个射频脉冲的相位角分别等于0°产生的受激回波，与3个相位角分别等于0°，90°，90°的射频脉冲产生的受激回波方向相反.
1.2 受激回波校正原理基于前文分析的受激回波形成原理及射频脉冲相位角对受激回波的影响，本文设计的受激回波和自旋回波分离算法如图 3所示.图 3a为施加相位角等于0°、翻转角等于90°的射频脉冲对被成像物体进行激发，其中tet代表回波时间.在tet/2时间施加翻转角等于180°、相位角等于0°的回聚射频脉冲，在3tet/2时间再一次施加翻转角等于180°、相位角等于0°的回聚射频脉冲，在2tet时间采集信号，这时采集到的信号S1既包含自旋回波信号SSE1，又包含受激回波信号SSTE1，可表示为
 | (1) |
 | 图 3 自旋回波、受激回波分离原理图Fig.3 Schematic diagram of spin echo and stimulated echo separation (a)—射频脉冲0°相位角激发，0°相位角回聚；(b)—射频脉冲0°相位角激发，90°相位角回聚. |
图 3b为施加翻转角等于90°、相位角等于0°的射频脉冲对被成像物体进行激发，之后的180°回聚脉冲相位角均为90°，这时采集到的信号同样包含自旋回波信号SSE2和受激回波信号SSTE2，可表示为
 | (2) |
 | (3) |
 | (4) |
 | (5) |
由于回聚脉冲的翻转角不能达到理想180°，使自旋回波的相位角φSE和受激回波的相位角φSTE不一致，因此造成图像伪影.对180°回聚脉冲的相位修正为
 | (6) |
FSE脉冲序列的时序如图 4所示，Gx，Gy和Gz分别为读梯度、相位编码梯度和选层梯度，Tx/Rx为射频发射和信号接收，ACQ为采集窗.脉冲序列的扫描参数：重复时间为3 600 ms，回波时间为105 ms，扫描矩阵为256×256，回波链长度为16.
图 4(Fig. 4)
 | 图 4 快速自旋回波脉冲序列时序图Fig.4 The pulse sequence diagram of fast spin echo |
2.2 数据分析为了定量分析受激回波伪影校正对快速自旋回波图像的伪影及信噪比影响，本文以国家食品药品监督管理局发布的“YY/T 0482—2010医用成像磁共振设备主要图像质量参数的测定”为依据，采用如下方法进行定量分析.
对于伪影，计算伪影与信号的比值.采用3个测量值：①在伪影处画感兴趣区域(region of interest，ROI)，找到最大伪影值Gmax；②在测量模体内画ROI，计算信号平均值Savg；③在背景处画ROI，计算背景噪声平均值Navg，伪影与信号的比值g为
 | (7) |
 | (8) |
图 5(Fig. 5)
 | 图 5 FSE水模图像Fig.5 The phantom images of fast spin echo (a)—没有校正；(b)—校正后. |
为了定量分析本文所提出的方法对图像的影响，对圆柱状水模进行了成像，并在图像上画了3个ROI，如图 6所示.其中ROI 1用于标记水模图像的统计信息，ROI 2用于标记背景噪声的统计信息，ROI 3用于标记伪影的统计信息，ROI的统计值见表 1.
图 6(Fig. 6)
 | 图 6 快速自旋回波图像伪影及信噪比测量Fig.6 Ghost and signal to noise ratio measurement of fast spin echo images (a)—校正后；(b)—没有校正. |
表 1(Table 1)
 表 1 ROI统计信息Table 1 The statistics values of ROI
| 表 1 ROI统计信息 Table 1 The statistics values of ROI |
伪影和信噪比的计算结果见表 2.从表 2可以看出，经过受激回波伪影校正后，快速自旋回波图像的信噪比约提高了57.08%，而伪影约降低了76.12%.
表 2(Table 2)
表 2 受激回波伪影校正前后的伪影和信噪比定量对比Table 2 Quantitative comparison of ghost and signal to noise ratio with or without stimulated echo artifact correction
| 表 2 受激回波伪影校正前后的伪影和信噪比定量对比 Table 2 Quantitative comparison of ghost and signal to noise ratio with or without stimulated echo artifact correction |
如图 7所示，FSE脉冲序列共激发5个成像层面，图中显示了每一个层面校正前、后自旋回波与受激回波的回波中心与相位的变化情况.从图 7a可以看出，校正前每一个层面的回波中心位置与相位均有较大差别，而经过校正后，图 7b中所有层面的回波中心位置与相位基本一致.
图 7(Fig. 7)
 | 图 7 校正前后相位与回波中心变化情况Fig.7 The changes of phase and echo center before and after correction (a)—校正前；(b)—校正后. |
由于回聚脉冲的翻转角不是理想的180°，受激回波和自旋回波的相位不一致，接收到的信号是二者的合成信号，所以会造成图像信号降低.而经过校正后，受激回波和自旋回波的相位基本一致，接收到的信号不仅有自旋回波信号，还在相同的相位方向增加了受激回波信号，信号强度增大，所以图像的信噪比得到了提高.
从第二个180°回聚脉冲开始，根据回波链的长度，每次回聚，散相的磁化矢量都会产生受激回波，受激回波被施加的空间编码梯度所编码，采样到k空间，扰乱了k空间的编码轨迹.在通过傅里叶变换进行图像重建时，把像素在图像域映射到了错误的位置，因此形成伪影.由于形成过程复杂，受激回波伪影很难通过图像后处理进行校正.
4 结论1) 本文根据受激回波的形成过程，分析出了射频脉冲的相位角对磁化矢量的影响，并以此为依据，分解出受激回波和自旋回波.
2) 通过校正读梯度方向的预散相梯度面积和180°回聚脉冲的相位角，可以实现受激回波伪影校正.
3) 从快速自旋回波的模体图像可以看出，本文提出的方法能有效降低图像伪影，提高图像信噪比.
参考文献
| [1] | Hennig J, Nauerth A, Friedburg H. RARE imaging:a fast imaging method for clinical MR[J]. Magnetic Resonance in Medicine, 1986, 3(6): 823-833. DOI:10.1002/mrm.1910030602 |
| [2] | Olivieri I, Salvarani C, Cantini F, et al. Fast spin echo-T2-weighted sequences with fat saturation in dactylitis of spondylarthritis[J]. Clinical Rheumatology, 2002, 46(11): 2964-2967. |
| [3] | Narumi S, Sasaki M, Miyazawa H, et al. T1-weighted magnetic resonance carotid plaque imaging:a comparison between conventional and fast spin-echo techniques[J]. Journal of Stroke & Cerebrovascular Diseases, 2017, 26(2): 273-279. |
| [4] | Sonin A H, Pensy R A, Mulligan M E, et al. Grading articular cartilage of the knee using fast spin-echo proton density-weighted MR imaging without fat suppression[J]. American Journal of Roentgenol, 2002, 179(5): 1159-1166. DOI:10.2214/ajr.179.5.1791159 |
| [5] | Chan R W, Lau A Z, Detzler G, et al. Evaluating the accuracy of multicomponent T2 parameters for luminal water imaging of the prostate with acceleration using inner-volume 3DGRASE[J]. Magnetic Resonance in Medicine, 2019, 81(1): 466-476. DOI:10.1002/mrm.27372 |
| [6] | Hennig J. Multiecho imaging sequences with low refocusing flip angles[J]. Journal of Magnetic Resonance, 1988, 78(3): 397-407. |
| [7] | Larson P E Z, Kerr A B, Reed G D, et al. Generating super stimulated-echoes in MRI and their application to hyperpolarized C-13 diffusion metabolic imaging[J]. IEEE Transactions on Medical Imaging, 2012, 31(2): 265-275. DOI:10.1109/TMI.2011.2168235 |
| [8] | Bernstein M A, King K F, Zhou X J. Handbook of MRI pulse sequences[M]. London: Academic Press, 2004: 305-316. |
| [9] | Nagy Z, Thomas D L, Weiskopf N. Orthogonalizing crusher and diffusion-encoding gradients to suppress undesired echo pathways in the twice-refocused spin echo diffusion sequence[J]. Magnetic Resonance in Medicine, 2014, 71(2): 506-515. DOI:10.1002/mrm.24676 |
| [10] | Kleinnijenhuis M, Mollink J, Lam W W, et al. Choice of reference measurements affects quantification of long diffusion time behaviour using stimulated echoes[J]. Magnetic Resonance in Medicine, 2017, 79(2): 952-959. |
| [11] | 唐昕, 姜忠德, 陈铭明, 等.磁共振成像系统快速自旋回波的相位校正方法及其应用: 201410834223.8[P].2015-04-29. (Tang Xin, Jiang Zhong-de, Chen Ming-ming, et al.Phase correction method of fast spin echo in magnetic resonance imaging scanner and its application: 201410834223.8[P].2015-04-29. ) |
| [12] | 宁瑞鹏, 李钰, 胡彦, 等.一种快速自旋回波脉冲序列中射频脉冲相位的优化方法: 201310207957.9[P].2013-05-30. (Ning Rui-peng, Li Yu, Hu Yan, et al.The phase optimization method of radio frequency pulse for fast spin echo pulse sequence: 201310207957.9[P].2013-05-30. ) |
| [13] | 宁瑞鹏, 杨光, 谢海滨, 等.一种快速自旋回波脉冲序列中射频脉冲的优化方法: 201910042695.2[P].2019-04-30. (Ning Rui-peng, Yang Guang, Xie Hai-bin, et al.A radio frequency pulse optimization method for fast spin echo pulse sequence: 201910042695.2[P].2019-04-30. ) |
| [14] | Frahm J, Merboldt K D, H?nicke W, et al. Stimulated echo imaging[J]. Journal of Magnetic Resonance, 1985, 64(1): 81-93. |
| [15] | Nessaiver M. All you really need to know about MRI physics[M]. Baltimore: Simply Physics, 1996: 71-74. |
