CN101357063A - 一种磁共振快速自旋回波成像方法 - Google Patents
一种磁共振快速自旋回波成像方法 Download PDFInfo
- Publication number
- CN101357063A CN101357063A CNA2008100423161A CN200810042316A CN101357063A CN 101357063 A CN101357063 A CN 101357063A CN A2008100423161 A CNA2008100423161 A CN A2008100423161A CN 200810042316 A CN200810042316 A CN 200810042316A CN 101357063 A CN101357063 A CN 101357063A
- Authority
- CN
- China
- Prior art keywords
- echo
- amplitude
- signal peak
- phase
- magnetic resonance
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明涉及磁共振成像技术领域,具体指一种磁共振快速自旋回波(FSE)成像方法。其包括:一、将相位编码梯度幅度置为零,采用FSE脉冲序列进行预扫描,记录回波链上各回波信号峰点的幅度值;二、调节会聚脉冲的相位和倾倒角,使回波链上各回波信号峰点的幅度达到最大值,记录此时各回波信号峰点的幅度值Am和相位值Φm;三、恢复相位编码梯度,采用FSE脉冲序列进行扫描,得到K空间数据;四、利用步骤二得到的各回波信号峰点的幅度Am和相位Φm,对K空间数据进行校正;五、利用经过校正的K空间数据进行图像重建,得到磁共振图像等步骤。解决了现有技术存在“回波幅度相位振荡”和“横向弛豫引起图像模糊”的问题。
Description
技术领域
本发明涉及磁共振成像技术领域,具体指一种磁共振快速自旋回波(FSE)成像方法。
背景技术
磁共振成像(MRI)技术已经成为医学诊断中非常有用的手段。通常,在磁共振成像系统中,当被测样品(如人体组织)处于静磁场B0(B0方向定义为直角坐标系的Z轴方向)中达到平衡时,样品中的原子核(核自旋)因被B0极化而产生一个宏观的磁化矢量M0。M0在射频脉冲的激发下被旋转到水平面(XY平面)里,然后绕Z轴作进动。在被测样品周围放置一个接收线圈,它就会感应出磁化矢量进动信号。接收线圈得到的磁共振信号经过放大以及模数转换(A/D)后,进入计算机进行图像重建。一般而言,为了进行磁共振成像,系统还需产生三路正交的梯度磁场,以便对磁共振信号进行空间编码。
采用常规磁共振成像序列进行扫描,单次扫描时间可达几分钟。如此长的扫描时间不仅会给病人带来不适,而且由于呼吸、心搏、胃肠蠕动以及某些自主运动所造成的伪影,使得图像质量大大恶化。采用快速成像序列可以将单次扫描时间缩短至几秒钟,从而减少图像的运动伪影和病人在扫描过程中的不适。快速自旋回波(FSE)序列(回波链长ETL>=2)是快速成像序列中的一种,已在多数MRI系统中得到应用。
然而,FSE序列在实际应用中存在如下两个主要问题:(1)由于非理想会聚脉冲的存在,回波链上各回波信号的幅度和相位产生振荡,从而导致图像中出现伪影;(2)当FSE序列的回波链较长或者回波时间较长时,由于横向弛豫的影响,回波链上的各回波信号依次衰减,造成图像模糊。
有人提出一种改进方法,它只采集回波信号的两个分量中的一个,从而避免回波信号相位的振荡。但是,该方法只解决了问题(1),而且因为该方法放弃了一半信号,所以图像信噪比大幅下降。
后来,有人提出在未施加相位编码梯度的条件下,调节会聚脉冲的相位和倾倒角度,使回波链上各回波信号的幅度一致,再根据会聚脉冲的相位对回波信号的相位进行校正。该方法消除了信号幅度的振荡,补偿了各回波之间相位的差异,削弱了由信号幅度和相位振荡所引起的图像伪影。并且在一定程度上减小了由于横向弛豫的影响,回波链上的各回波信号依次衰减,而造成的图像模糊。然而该方法仍然存在以下问题:为了使回波信号的幅度一致,需要压制回波链上大部分的回波信号。这些被压制的回波信号通常对应着图像中慢变化的部分,对图像信噪比的贡献非常大。因此,压制信号会造成图像信噪比大幅下降。
发明内容
本发明的目的是解决上述现有技术存在“回波幅度相位振荡”和“横向弛豫引起图像模糊”的问题,而提供一种磁共振快速自旋回波(FSE)成像方法。本发明包括以下步骤:(1)将相位编码梯度幅度置为零,采用FSE脉冲序列进行预扫描,记录回波链上各回波信号峰点的幅度值;(2)调节会聚脉冲的相位和倾倒角,使回波链上各回波信号峰点的幅度达到最大值,记录此时各回波信号峰点的幅度值Am和相位值Φm;(3)恢复相位编码梯度,采用FSE脉冲序列进行扫描,得到K空间数据;(4)利用步骤(2)得到的各回波信号峰点的幅度Am和相位Φm,对K空间数据进行校正;(5)利用经过校正的K空间数据进行图像重建,得到磁共振图像。
本发明还包括以下步骤:(1)将相位编码梯度幅度置为零,采用FSE脉冲序列进行预扫描,记录回波链上各回波信号峰点的幅度值;(2)调节会聚脉冲的相位和倾倒角,使回波链上各回波信号峰点的幅度达到最大值,记录此时各回波信号峰点的幅度值Am和相位值Φm;(3)利用步骤(2)得到的相位值Φm设置接收机的相位列表,恢复相位编码梯度,采用FSE脉冲序列进行扫描,得到K空间数据;(4)利用步骤(2)得到的各回波信号峰点的幅度Am,对K空间数据进行校正;(5)利用经过校正的K空间数据进行图像重建,得到磁共振图像。
本方法与现有技术相比,不会造成图像信噪比大幅下降,特别适用于单激发FSE成像以及回波链较长或者回波时间较长的FSE成像。
附图说明
图1为本发明一种磁共振快速自旋回波成像方法的MRI系统框图;
图2为本发明一种磁共振快速自旋回波成像方法的快速自旋回波序列示意图;
图3为本发明一种磁共振快速自旋回波成像方法的一实施例流程图;
图4为本发明一种磁共振快速自旋回波成像方法的另一实施例流程图。
图例说明:
RF:射频脉冲;90:激发脉冲;180:反转脉冲;Gs:选层梯度;Gp:相位编码梯度;Gr:读梯度;Echo:回波
具体实施方式
以下结合附图对本发明作进一步描述
本发明基于如附图1所示的磁共振成像(MRI)系统结构,并与其结合作如下操作步骤:
(1)将相位编码梯度幅度置为零,采用FSE脉冲序列进行预扫描,记录回波链上各回波信号峰点的幅度值;
(2)调节会聚脉冲的相位和倾倒角,使回波链上各回波信号峰点的幅度达到最大值,记录此时各回波信号峰点的幅度值Am和相位值Φm;
(3)恢复相位编码梯度,采用FSE脉冲序列进行扫描,得到K空间数据;
(4)利用步骤(2)得到的各回波信号峰点的幅度Am和相位Φm,对K空间数据进行校正;
(5)利用经过校正的K空间数据进行图像重建,得到磁共振图像;
或者:
(1)将相位编码梯度幅度置为零,采用FSE脉冲序列进行预扫描,记录回波链上各回波信号峰点的幅度值;
(2)调节会聚脉冲的相位和倾倒角,使回波链上各回波信号峰点的幅度达到最大值,记录此时各回波信号峰点的幅度值Am和相位值Φm;
(3)利用步骤(2)得到的相位值Φm设置接收机的相位列表,恢复相位编码梯度,采用FSE脉冲序列进行扫描,得到K空间数据;
(4)利用步骤(2)得到的各回波信号峰点的幅度Am,对K空间数据进行校正;
(5)利用经过校正的K空间数据进行图像重建,得到磁共振图像。
为对本发明的一种磁共振快速自旋回波成像方法,有进一步的了解,现结合快速自旋回波序列(如附图2所示,为简要起见,附图中只给出了5个回波信号201、202……205。实际扫描中,回波数目(ETL)不限于此5个),进行如下详细说明:
本发明一种实施方式的流程图(如附图3所示)
1.将相位编码梯度幅度置为零,执行附图2所示的脉冲序列,得到ETL个回波。其中,ETL为回波链长。
2.调节会聚脉冲的相位和倾倒角,使回波链上每一个回波信号峰点的幅度达到最大值,记录此时各回波信号峰点的幅度值Am和相位值Φm。其中,m为回波的序数,取1,2……ETL。
3.恢复相位编码梯度,执行附图2所示的脉冲序列,得到K空间数据。
4.1对K空间数据分块分组,并排序,其方法如下:
若相位编码次数为n×ETL(通常相位编码次数是回波链长的整数倍n),则可以将K空间数据分为n块,每一块数据记录的是同一个回波链上的信号。
每一块数据可以分为ETL组,每一组数据记录的是回波链上的一个回波;根据扫描方案,可以对各组在数据块内进行排序;经过排序,每个数据块中,第1组数据表示的回波的幅度与A1相对应,第2组数据表示的回波的幅度与A2相对应,……第ETL组数据表示的回波的幅度与AETL相对应。
4.2对排序后的K空间数据进行幅度校正,其方法如下:
令B1,B2,……BETL满足如下条件:A1×B1=A2×B2=……=AETL×BETL。每一个数据块中,第1组的每一个数据(包括实部和虚部)均与B1相乘,第2组的每一个数据(包括实部和虚部)均与B2相乘,……第ETL组的每一个数据(包括实部和虚部)均与BETL相乘。
4.3对排序后的K空间数据进行相位校正,其方法如下:
每一个数据块中,第1组的每一个数据的相位均与Φ1相减,第2组的每一个数据的相位均与Φ2相减,……第ETL组的每一个数据的相位均与ΦETL相减。
4.4保留经过上一流程4.3和流程4.4校正的K空间数据,将排序恢复至流程3完成时的状态。
5.利用经过校正的K空间数据进行图像重建,得到磁共振图像。
上述流程中,流程4.2与4.3可以调换次序,即流程4.1之后,先执行流程4.3,然后执行流程4.2,其余各流程次序没有变化。
对于多层扫描,对每一层分别执行上述全部流程即可。在图像质量要求不高的情况下,首先选择其中一层执行流程1至流程2,然后执行流程3进行多层扫描,再分别对每一层执行流程4.1至5。
本发明另一种实施方式的流程图(如附图4所示)
1.将相位编码梯度幅度置为零,执行附图2所示的脉冲序列,得到ETL个回波;其中,ETL为回波链长。
2.调节会聚脉冲的相位和倾倒角,使回波链上每一个回波信号峰点的幅度达到最大值,记录此时各回波信号峰点的幅度值Am和相位值Φm;其中,m为回波的序数,取1,2……ETL。
3.1.将接收机的相位列表数组设置为[Φ1,Φ2,……ΦETL]。
3.2.恢复相位编码梯度,执行附图2所示的脉冲序列,得到K空间数据。
4.1.对K空间数据分块分组,并排序,其方法如下:
若相位编码次数为n×ETL(通常相位编码次数是回波链长的整数倍n),则可以将K空间数据分为n块,每一块数据记录的是同一个回波链上的信号。
每一块数据可以分为ETL组,每一组数据记录的是回波链上的一个回波;根据扫描方案,可以对各组在数据块内进行排序;所有数据块经过排序后,每个数据块中,第1组数据表示的回波的幅度与A1相对应,第2组数据表示的回波的幅度与A2相对应,……第ETL组数据表示的回波的幅度与AETL相对应。
4.2.对排序后的K空间数据进行幅度校正,其方法如下:
令B1,B2,……BETL满足如下条件:A1×B1=A2×B2=……=AETL×BETL。每一个数据块中,第1组的每一个数据(包括实部和虚部)均与B1相乘,第2组的每一个数据(包括实部和虚部)均与B2相乘,……第ETL组的每一个数据(包括实部和虚部)均与BETL相乘。
4.3.保留经过上一流程4.2校正的K空间数据,将排序恢复至流程3完成时的状态。
5.利用经过校正的K空间数据进行图像重建,得到磁共振图像。
对于多层扫描,对每一层分别执行上述全部流程即可。在图像质量要求不高的情况下,首先选择其中一层执行流程1至流程3.1,然后执行流程3.2进行多层扫描,再分别对每一层执行流程4.1至流程5。
综上所述,本方法与现有技术相比,解决了因“回波幅度相位振荡”和“横向弛豫引起图像模糊”造成图像信噪比有大幅下降。其特别适用于单激发FSE成像以及回波链较长或者回波时间较长的FSE成像。
Claims (1)
1.一种磁共振快速自旋回波成像方法,其特征,包括以下步骤:(1)将相位编码梯度幅度置为零,采用FSE脉冲序列进行预扫描,记录回波链上各回波信号峰点的幅度值;
(2)调节会聚脉冲的相位和倾倒角,使回波链上各回波信号峰点的幅度达到最大值,记录此时各回波信号峰点的幅度值Am和相位值Фm;
(3)恢复相位编码梯度,采用FSE脉冲序列进行扫描,得到K空间数据;
(4)利用步骤(2)得到的各回波信号峰点的幅度Am和相位Фm,对K空间数据进行校正;
(5)利用经过校正的K空间数据进行图像重建,得到磁共振图像;
或者:
(1)将相位编码梯度幅度置为零,采用FSE脉冲序列进行预扫描,记录回波链上各回波信号峰点的幅度值;
(2)调节会聚脉冲的相位和倾倒角,使回波链上各回波信号峰点的幅度达到最大值,记录此时各回波信号峰点的幅度值Am和相位值Фm;
(3)利用步骤(2)得到的相位值Фm设置接收机的相位列表,恢复相位编码梯度,采用FSE脉冲序列进行扫描,得到K空间数据;
(4)利用步骤(2)得到的各回波信号峰点的幅度Am,对K空间数据进行校正;
(5)利用经过校正的K空间数据进行图像重建,得到磁共振图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008100423161A CN101357063B (zh) | 2008-08-29 | 2008-08-29 | 一种磁共振快速自旋回波成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2008100423161A CN101357063B (zh) | 2008-08-29 | 2008-08-29 | 一种磁共振快速自旋回波成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101357063A true CN101357063A (zh) | 2009-02-04 |
CN101357063B CN101357063B (zh) | 2010-06-16 |
Family
ID=40329758
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2008100423161A Expired - Fee Related CN101357063B (zh) | 2008-08-29 | 2008-08-29 | 一种磁共振快速自旋回波成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101357063B (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102096054A (zh) * | 2010-12-22 | 2011-06-15 | 华东师范大学 | 一种快速自旋回波脉冲序列的梯度系统延时校正方法 |
CN102890255A (zh) * | 2012-10-19 | 2013-01-23 | 浙江大学 | 一种用于缩短回波时间的由内而外平面回波成像方法 |
CN103901376A (zh) * | 2012-12-30 | 2014-07-02 | 上海联影医疗科技有限公司 | 磁共振成像方法与装置 |
CN104459587A (zh) * | 2013-09-17 | 2015-03-25 | 华润万东医疗装备股份有限公司 | 一种用于核磁共振成像系统的流动补偿方法 |
CN104569882A (zh) * | 2014-12-29 | 2015-04-29 | 苏州朗润医疗系统有限公司 | 磁共振成像系统快速自旋回波的相位校正方法及其应用 |
CN105929350A (zh) * | 2016-05-05 | 2016-09-07 | 大连锐谱科技有限责任公司 | 一种单次激发水脂分离成像误差校正系统及方法 |
CN110058185A (zh) * | 2019-04-02 | 2019-07-26 | 佛山瑞加图医疗科技有限公司 | 一种磁共振旋转成像方法及系统 |
CN110794352A (zh) * | 2018-08-02 | 2020-02-14 | 上海康达卡勒幅医疗科技有限公司 | 一种抑制磁共振图像伪影的方法 |
WO2020078131A1 (zh) * | 2018-10-17 | 2020-04-23 | 浙江大学 | 一种基于频率稳定模块的磁共振cest成像序列及装置 |
CN111157932A (zh) * | 2020-01-02 | 2020-05-15 | 华东师范大学 | 一种快速自旋回波脉冲序列中射频脉冲的优化方法 |
CN111352055A (zh) * | 2020-03-17 | 2020-06-30 | 无锡鸣石峻致医疗科技有限公司 | 一种前瞻性相位校正平面回波成像技术 |
WO2021197365A1 (zh) * | 2020-03-31 | 2021-10-07 | 浙江大学 | 磁共振cest成像频率漂移校正方法、装置、介质及成像设备 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100580476C (zh) * | 2007-11-27 | 2010-01-13 | 华东师范大学 | 一种快速自旋回波脉冲序列的校正方法及其应用 |
-
2008
- 2008-08-29 CN CN2008100423161A patent/CN101357063B/zh not_active Expired - Fee Related
Cited By (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102096054A (zh) * | 2010-12-22 | 2011-06-15 | 华东师范大学 | 一种快速自旋回波脉冲序列的梯度系统延时校正方法 |
CN102096054B (zh) * | 2010-12-22 | 2013-03-13 | 华东师范大学 | 一种快速自旋回波脉冲序列的梯度系统延时校正方法 |
CN102890255A (zh) * | 2012-10-19 | 2013-01-23 | 浙江大学 | 一种用于缩短回波时间的由内而外平面回波成像方法 |
CN103901376A (zh) * | 2012-12-30 | 2014-07-02 | 上海联影医疗科技有限公司 | 磁共振成像方法与装置 |
CN103901376B (zh) * | 2012-12-30 | 2017-11-07 | 深圳联影医疗科技有限公司 | 磁共振成像方法与装置 |
CN104459587A (zh) * | 2013-09-17 | 2015-03-25 | 华润万东医疗装备股份有限公司 | 一种用于核磁共振成像系统的流动补偿方法 |
CN104459587B (zh) * | 2013-09-17 | 2018-02-09 | 北京万东医疗科技股份有限公司 | 一种用于核磁共振成像系统的流动补偿方法 |
CN104569882A (zh) * | 2014-12-29 | 2015-04-29 | 苏州朗润医疗系统有限公司 | 磁共振成像系统快速自旋回波的相位校正方法及其应用 |
CN104569882B (zh) * | 2014-12-29 | 2017-08-22 | 苏州朗润医疗系统有限公司 | 磁共振成像系统快速自旋回波的相位校正方法及其应用 |
CN105929350A (zh) * | 2016-05-05 | 2016-09-07 | 大连锐谱科技有限责任公司 | 一种单次激发水脂分离成像误差校正系统及方法 |
CN105929350B (zh) * | 2016-05-05 | 2019-01-25 | 南京拓谱医疗科技有限公司 | 一种单次激发水脂分离成像误差校正系统及方法 |
CN110794352A (zh) * | 2018-08-02 | 2020-02-14 | 上海康达卡勒幅医疗科技有限公司 | 一种抑制磁共振图像伪影的方法 |
CN110794352B (zh) * | 2018-08-02 | 2023-05-26 | 上海康达卡勒幅医疗科技有限公司 | 一种抑制磁共振图像伪影的方法 |
WO2020078131A1 (zh) * | 2018-10-17 | 2020-04-23 | 浙江大学 | 一种基于频率稳定模块的磁共振cest成像序列及装置 |
US11237239B2 (en) | 2018-10-17 | 2022-02-01 | Zhejiang University | Magnetic resonance cest imaging sequence and device based on frequency stabilization module |
CN110058185A (zh) * | 2019-04-02 | 2019-07-26 | 佛山瑞加图医疗科技有限公司 | 一种磁共振旋转成像方法及系统 |
CN110058185B (zh) * | 2019-04-02 | 2021-01-15 | 佛山瑞加图医疗科技有限公司 | 一种磁共振旋转成像方法及系统 |
CN111157932A (zh) * | 2020-01-02 | 2020-05-15 | 华东师范大学 | 一种快速自旋回波脉冲序列中射频脉冲的优化方法 |
CN111157932B (zh) * | 2020-01-02 | 2022-08-30 | 华东师范大学 | 一种快速自旋回波脉冲序列中射频脉冲的优化方法 |
CN111352055A (zh) * | 2020-03-17 | 2020-06-30 | 无锡鸣石峻致医疗科技有限公司 | 一种前瞻性相位校正平面回波成像技术 |
CN111352055B (zh) * | 2020-03-17 | 2022-05-06 | 无锡鸣石峻致医疗科技有限公司 | 一种前瞻性相位校正平面回波成像方法 |
WO2021197365A1 (zh) * | 2020-03-31 | 2021-10-07 | 浙江大学 | 磁共振cest成像频率漂移校正方法、装置、介质及成像设备 |
JP7408184B2 (ja) | 2020-03-31 | 2024-01-05 | 浙江大学 | 磁気共鳴cest結像周波数ドリフト補正方法、装置、媒体及び結像デバイス |
US11927659B2 (en) | 2020-03-31 | 2024-03-12 | Zhejiang University | Method and apparatus for frequency drift correction of magnetic resonance CEST imaging, and medium and imaging device |
Also Published As
Publication number | Publication date |
---|---|
CN101357063B (zh) | 2010-06-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101357063B (zh) | 一种磁共振快速自旋回波成像方法 | |
US7382127B2 (en) | System and method of accelerated MR propeller imaging | |
US9229082B2 (en) | Method for diffusion magnetic resonance imaging | |
US8023705B2 (en) | Method for reconstructing image from echo planar imaging sequence | |
US7688068B2 (en) | System and method for using parallel imaging with compressed sensing | |
CN103238082B (zh) | 使用多点Dixon技术和低分辨率校准的MR成像 | |
JP5443695B2 (ja) | 磁気共鳴イメージング装置 | |
CN101484823B (zh) | 磁共振成像设备和方法 | |
CN108335339A (zh) | 一种基于深度学习和凸集投影的磁共振重建方法 | |
US20120002859A1 (en) | Magnetic resonance partially parallel imaging (ppi) with motion corrected coil sensitivities | |
US20140086469A1 (en) | Mri reconstruction with incoherent sampling and redundant haar wavelets | |
CN103969611A (zh) | 用于运行磁共振系统的方法和控制装置 | |
CN106574954A (zh) | 针对epi的具有奈奎斯特伪影校正的并行mr成像 | |
JP6285669B2 (ja) | 磁気共鳴イメージング装置 | |
CN112368715A (zh) | 用于磁共振成像的运动校正的方法和系统 | |
CN105143906A (zh) | 金属抗mr成像 | |
US10495717B2 (en) | System and method for dual-kernel image reconstruction | |
US11131738B2 (en) | Multiband, multishot magnetic resonance elastography | |
CN110895320B (zh) | 基于深度学习的无变形单次激发平面回波成像方法及装置 | |
CN105738846A (zh) | K空间数据采集方法及其磁共振成像方法 | |
CN102096054A (zh) | 一种快速自旋回波脉冲序列的梯度系统延时校正方法 | |
CN108594146A (zh) | 一种基于深度学习和凸集投影的磁共振弥散加权成像方法 | |
CN108459289A (zh) | 一种基于数据一致性的多次激发弥散加权磁共振成像方法 | |
JP5254633B2 (ja) | Mri装置およびその制御方法 | |
CN115639510A (zh) | 磁共振成像方法、波谱成像方法、装置、设备和存储介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20100616 Termination date: 20170829 |