CN111381204B - 基于二维快速自旋回波的磁共振成像方法和装置 - Google Patents
基于二维快速自旋回波的磁共振成像方法和装置 Download PDFInfo
- Publication number
- CN111381204B CN111381204B CN201811618720.9A CN201811618720A CN111381204B CN 111381204 B CN111381204 B CN 111381204B CN 201811618720 A CN201811618720 A CN 201811618720A CN 111381204 B CN111381204 B CN 111381204B
- Authority
- CN
- China
- Prior art keywords
- magnetic resonance
- flip angle
- optimized
- echo signal
- evolution curve
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/543—Control of the operation of the MR system, e.g. setting of acquisition parameters prior to or during MR data acquisition, dynamic shimming, use of one or more scout images for scan plane prescription
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/05—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves
- A61B5/055—Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/483—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/5602—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by filtering or weighting based on different relaxation times within the sample, e.g. T1 weighting using an inversion pulse
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/58—Calibration of imaging systems, e.g. using test probes, Phantoms; Calibration objects or fiducial markers such as active or passive RF coils surrounding an MR active material
- G01R33/583—Calibration of signal excitation or detection systems, e.g. for optimal RF excitation power or frequency
- G01R33/586—Calibration of signal excitation or detection systems, e.g. for optimal RF excitation power or frequency for optimal flip angle of RF pulses
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/4818—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space
- G01R33/4824—MR characterised by data acquisition along a specific k-space trajectory or by the temporal order of k-space coverage, e.g. centric or segmented coverage of k-space using a non-Cartesian trajectory
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/50—NMR imaging systems based on the determination of relaxation times, e.g. T1 measurement by IR sequences; T2 measurement by multiple-echo sequences
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
- G01R33/5615—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
- G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
- G01R33/561—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reduction of the scanning time, i.e. fast acquiring systems, e.g. using echo-planar pulse sequences
- G01R33/5615—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE]
- G01R33/5617—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE] using RF refocusing, e.g. RARE
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
- G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
- G01R33/48—NMR imaging systems
- G01R33/58—Calibration of imaging systems, e.g. using test probes, Phantoms; Calibration objects or fiducial markers such as active or passive RF coils surrounding an MR active material
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- High Energy & Nuclear Physics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- Signal Processing (AREA)
- Biomedical Technology (AREA)
- Pathology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Biophysics (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Optics & Photonics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明实施例中公开了一种基于二维快速自旋回波的磁共振成像方法和装置。其中,方法包括:确定用于预设一磁共振图像的初始对比度的第一调整参数;确定用于得到磁共振图像的优化对比度及该刀锋伪影校正序列的设定数据采集时间的第二调整参数;根据所述第一调整参数和所述第二调整参数来确定一优化的回波信号演化曲线;根据所述优化的回波信号演化曲线计算得到实际可变翻转角链;将所述实际可变翻转角链应用于二维快速自旋回波序列,并采用与所述第二调整参数相应的刀锋伪影校正序列来采集磁共振信号并使磁共振图像满足所述优化对比度。本发明实施例中的技术方案能够在通过增加叶片宽度来加快采集速度以消除伪影时,保证图像的成像质量。
Description
技术领域
本发明涉及磁共振成像技术领域,特别是一种基于二维快速自旋回波的磁共振成像方法和装置。
背景技术
磁共振成像(Magnetic Resonance Imaging,MRI)是利用磁共振现象进行成像的一种技术。磁共振现象的原理主要包括:包含单数质子的原子核,例如人体内广泛存在的氢原子核,其质子具有自旋运动,犹如一个小磁体,并且这些小磁体的自旋轴无一定的规律,如果施加外在磁场,这些小磁体将按外在磁场的磁力线重新排列,具体为在平行于或反平行于外在磁场磁力线的两个方向排列,将上述平行于外在磁场磁力线的方向称为正纵向轴,将上述反平行于外在磁场磁力线的方向称为负纵向轴,原子核只具有纵向磁化分量,该纵向磁化分量既具有方向又具有幅度。用特定频率的射频(RF,Radio Frequency)脉冲激发处于外在磁场中的原子核,使这些原子核的自旋轴偏离正纵向轴或负纵向轴,产生共振,这就是磁共振现象。上述被激发原子核的自旋轴偏离正纵向轴或负纵向轴之后,原子核具有了横向磁化分量。
停止发射射频脉冲后,被激发的原子核发射回波信号,将吸收的能量逐步以电磁波的形式释放出来,其相位和能级都恢复到激发前的状态,将原子核发射的回波信号经过空间编码等进一步处理即可重建图像。上述被激发原子核向激发前状态的恢复过程称为驰豫过程,恢复到平衡状态所需的时间称为驰豫时间。
MRI系统是无创性的,与计算机断层扫描(CT)相比具有较高的组织对比度,并且不因骨组织而产生伪影。此外,MRI系统可在不改变扫描对象的位置的情况下捕获沿期望方向的各种横断面。
k空间是每个横断面的数据空间,通过对k空间执行傅里叶变换,可获得期望的图像。如果相位编码梯度和频率编码梯度的大小在施加RF脉冲序列之后逐渐的改变,则可获得具有各种位置信息的原数据。原数据具有位置信息以及组织对比度信息,k空间表示可形成一个图像的一组原数据。
图1为采用二维快速自旋回波(TSE)的刀锋伪影校正序列(BLADE)采集k空间数据的一个轨迹示意图。如图1所示,在一个重复时间(TR)内获取一个叶片,旋转一定角度(图1为20°)后获取下一个叶片。图1中以叶片宽度为L,对应该叶片宽度L的回波链长度为15的情况为例。其中,每个TR表示在脉冲序列中从90°脉冲到下一个90°脉冲的一个时间段。进行叶片扫描时,图像的基础分辨率m与叶片数量N、叶片宽度L之间的关系满足下述的公式(1),而采集时间Tacq与叶片数量N、重复时间TR之间的关系则满足下面的公式(2):
N*L=π*m/2 (1)
Tacq=N*TR (2)
由于在k空间的中心部分会过采样,因此叶片扫描是一种时间消耗大的数据采集方法。对于部分不配合的患者或进行脑部扫描时由于呼吸运动的存在,使得有些运动伪影仍然不能用BLADE方法进行矫正。甚至这些运动伪影有时会因模糊病理而导致错误的诊断。叶片扫描效率低,阻碍了叶片更多的诊断应用。因此,有必要加快叶片的运动速度来缓解运动伪影。根据公式(1)和(2)可以看出,在维持基础分辨率m不变的情况下,可以通过增加叶片宽度L来减少叶片数量N进而来减少采集时间Tacq。即增加叶片宽度L是提高叶片扫描效率的技术之一。然而叶片越宽,回波链长度越长,横向磁化矢量驰豫时间T2的衰减越多,进而导致图像模糊。
发明内容
有鉴于此,本发明实施例中一方面提出了一种基于二维快速自旋回波的磁共振成像方法,另一方面提出了一种基于二维快速自旋回波的磁共振成像装置,用以在通过增加叶片宽度来缩短采集时间以消除伪影时,保证图像的成像质量。
本发明实施例中提出的一种基于二维快速自旋回波的磁共振成像方法,其利用刀锋伪影校正序列(BLADE)进行磁共振信号的数据采集,包括:确定用于预设一磁共振图像的初始对比度的第一调整参数;确定用于得到磁共振图像的优化对比度及该刀锋伪影校正序列的设定数据采集时间的第二调整参数;根据所述第一调整参数和所述第二调整参数来确定一优化的回波信号演化曲线;根据所述优化的回波信号演化曲线计算得到实际可变翻转角链;将所述实际可变翻转角应用于二维快速自旋回波序列,并采用与所述第二调整参数相应的刀锋伪影校正序列来采集磁共振信号并使磁共振图像满足所述优化对比度。
在一个实施方式中,所述确定用于预设一磁共振图像的初始对比度的第一调整参数之后,还包括:根据所述第一调整参数来预定一初始回波信号演化曲线;根据所述初始回波信号演化曲线计算得到初始可变翻转角链。
在一个实施方式中,所述第二调整参数选自回波间距、回波链长度中的至少一个。
在一个实施方式中,所述根据所述优化的回波信号演化曲线计算得到实际可变翻转角链包括:根据所述优化回波信号演化曲线推算得到优化的可变翻转角;根据所述优化的可变翻转角计算得到一个计算的回波信号演化曲线;将所述计算的回波信号演化曲线和所述优化的回波信号演化曲线进行比较,在二者的差异满足设定要求时,将所述优化的可变翻转角用作所述实际可变翻转角链。
在一个实施方式中,所述第一调整参数为预设的如下三种可变翻转角模式中的一种:纵向磁化矢量驰豫时间T1权重、质子密度PD权重和横向磁化矢量驰豫时间T2权重。
本发明实施例中提出的一种基于二维快速自旋回波的磁共振成像装置,其利用刀锋伪影校正序列进行磁共振信号的数据采集,包括:第一调整参数确定模块,适用于确定用于预设一磁共振图像的初始对比度的第一调整参数;第二调整参数确定模块,适用于确定用于得到磁共振图像的优化对比度及该刀锋伪影校正序列的设定数据采集时间的第二调整参数;第一回波信号演化曲线确定模块,根据所述第一调整参数和所述第二调整参数来确定一优化的回波信号演化曲线;第一可变翻转角链确定模块,根据所述优化的回波信号演化曲线计算得到实际可变翻转角链;磁共振成像模块,将所述实际可变翻转角应用于二维快速自旋回波序列,并采用与所述第二调整参数相应的刀锋伪影校正序列来采集磁共振信号并使磁共振图像满足所述优化对比度。
在一个实施方式中,进一步包括:第二回波信号演化曲线确定模块,根据所述第一调整参数来预定一初始回波信号演化曲线;第二可变翻转角链确定模块,根据所述初始回波信号演化曲线计算得到初始可变翻转角链。
在一个实施方式中,所述第二调整参数确定模块确定用于得到磁共振图像的优化对比度的回波间距和/或回波链长度,以及用于得到该刀锋伪影校正序列的设定数据采集时的回波链长度。
在一个实施方式中,所述第一可变翻转角链确定模块包括:第一单元,用于根据所述预定的回波信号演化曲线推算得到优化的可变翻转角;第二单元,根据所述优化的可变翻转角计算得到一个计算的回波信号演化曲线;第三单元,将所述计算的回波信号演化曲线和所述优化的回波信号演化曲线进行比较,在二者的差异满足设定要求时,将所述优化的可变翻转角用作所述实际可变翻转角链。
在一个实施方式中,所述第一调整参数确定模块确定用于预设一磁共振图像的初始对比度的可变翻转角模式,其为如下三种可变翻转角模式中的一种:纵向磁化矢量驰豫时间T1权重、质子密度PD权重和横向磁化矢量驰豫时间T2权重。
从上述方案中可以看出,由于本发明实施例中首先确定一优化的回波信号演化曲线,基于该回波信号演化曲线能够使得磁共振图像满足预设的对比度,并使得BLADE序列的采集时间满足设定的短时要求,之后基于该优化的回波信号演化曲线计算得到可变翻转角链,并将该可变翻转角链应用于二维快速自旋回波序列,并采用与所述第二调整参数相应的刀锋伪影校正序列来采集磁共振信号。从而可在增加叶片宽度时,通过调整翻转角的值来控制T2的衰减,从而可以在通过增加叶片宽度来加快采集速度以消除伪影时,保证图像的成像质量。
此外,本发明实施例中给出了一种优化可变翻转角的简单可靠的实现方式。
附图说明
下面将通过参照附图详细描述本发明的优选实施例,使本领域的普通技术人员更清楚本发明的上述及其它特征和优点,附图中:
图1为采用TSE BLADE序列采集k空间数据的一个轨迹示意图。
图2为本发明实施例中一种基于二维快速自旋回波的磁共振成像方法的流程示意图。
图3A为本发明一个例子中一种预定的回波信号演化曲线的示意图。
图3B为本发明一个例子中一种优化的可变翻转角的示意图。
图4为本发明实施例中一种基于二维快速自旋回波的磁共振成像装置的结构示意图。
图5A至图5D为本发明一个例子中一组采用传统BLADE序列和采用本发明实施例中的带有翻转角的BLADE序列进行成像时获取的脑部图像的对比图。
图6A至图6D为本发明一个例子中一组采用传统BLADE序列和采用本发明实施例中的带有翻转角的BLADE序列进行成像时获取的膝盖图像的对比图。
图7A和图7B为本发明一个例子中一组采用传统BLADE序列和采用本发明实施例中的带有翻转角的BLADE序列进行成像时获取的手腕图像的对比图。
其中,附图标记如下:
标号 | 含义 |
L | 叶片宽度 |
201-205 | 步骤 |
410 | 第一调整参数确定模块 |
420 | 第二调整参数确定模块 |
430 | 回波信号演化曲线确定模块 |
440 | 可变翻转角链确定模块 |
441 | 第一单元 |
442 | 第二单元 |
443 | 第三单元 |
450 | 磁共振成像模块 |
具体实施方式
为了解决上述通过增加叶片宽度来加快采集速度以便消除伪影,然而又会导致T2衰减过多导致图像模糊的问题,本申请实施例中为BLADE序列引入了可变翻转角,以在增加叶片宽度时,通过调整翻转角的值来控制T2的衰减。即可变翻转角技术是一种在降低T2衰减效果的同时增加回波链长度的技术。
进一步地,针对可变翻转角模式,设置了用于调整组织对比度的成像参数,通过对成像参数进行调整可以实现满足所期望的组织对比度的图像。
为使本发明的目的、技术方案和优点更加清楚,以下举实施例对本发明进一步详细说明。
图2为本发明实施例中一种基于二维快速自旋回波的磁共振成像方法200的流程示意图。如图2所示,该方法200可包括如下步骤:
步骤201,确定用于预设一磁共振图像的初始对比度的第一调整参数。
为了实现所期望的磁共振图像组织对比度,可对已有的或新增的成像参数进行调整。例如,在一个实施方式中,可预先设置不同的可变翻转角模式,如纵向磁化矢量驰豫时间(T1)权重的可变翻转角模式、质子密度(PD)权重的可变翻转角模式和横向磁化矢量驰豫时间(T2)权重的可变翻转角模式;其中,所述T1权重的可变翻转角模式用于使回波信号演化曲线符合T1权重的磁共振成像,所述PD权重的可变翻转角模式用于使回波信号演化曲线符合PD权重的磁共振成像,所述T2权重的可变翻转角模式用于使回波信号演化曲线符合T2权重的磁共振成像。
相应地,本步骤中的第一调整参数可以是上述三个可变翻转角模式中的一个,例如,若期望得到符合T1权重的磁共振成像,则本步骤201可根据所接收的用户或系统选定的对应T1权重的可变翻转角模式的信号确定第一调整参数为T1权重的可变翻转角模式;若期望得到符合T2权重的磁共振成像,则本步骤201可根据所接收的用户或系统选定的对应T2权重的可变翻转角模式的信号确定第一调整参数为T2权重的可变翻转角模式;以此类推。
基于步骤201中确定的第一调整参数,本方法可进一步包括:根据所述第一调整参数来预定一初始回波信号演化曲线;根据所述初始回波信号演化曲线计算得到初始可变翻转角链。
步骤202,确定用于得到磁共振图像的优化对比度及该刀锋伪影校正序列的设定数据采集时间的第二调整参数。
具体实现时,执行步骤201后,可确定磁共振图像的初始对比度,若该初始对比度不能满足设定要求,则可进一步通过执行步骤202来进一步优化对比度,例如可通过调整回波间距和回波链长度中的至少一个来达到满足设定要求的对比度。此外,通过调整回波链长度还可以改变叶片宽度,进而可改变叶片数量和刀锋伪影校正序列的数据采集时间。其中,回波链越长则叶片宽度越大,叶片数量越小,刀锋伪影校正序列的数据采集时间越短。相应地,本步骤中可确定当前调整的回波间距和回波链长度中的至少一个的值。
步骤203,根据所述第一调整参数和所述第二调整参数来确定一优化的回波信号演化曲线。
本实施例中,为了使得回波信号演化曲线符合T2衰减要求,可通过步骤201和步骤202确定相应的参数,如可变翻转角模式,以及回波间距、回波链长度。
其中,优化的回波信号演化曲线即为基于第一调整参数和第二调整参数满足设定要求例如T2衰减要求的回波信号演化曲线。如图3A中的虚线曲线所示,图3A中示出了一个例子中一种预定的回波信号演化曲线的示意图。图3A中横坐标为回波链中的回波数Ne,其中,最大值127对应回波链长度。纵坐标为归一化的信号强度。
步骤204,根据所述优化的回波信号演化曲线计算得到实际可变翻转角链。
在一个例子中,本步骤204可采用如下方法实现:
A、根据所述预定的回波信号演化曲线推算得到优化的可变翻转角。
本步骤中,可基于扩展相位算法(EPG),得到如下计算翻转角的公式(3),基于该公式(3),将预定的回波信号演化曲线中的回波信号强度F1作为已知量,输入公式(3),可迭代计算出优化的翻转角:
其中,在αn为第n个回聚脉冲的翻转角度。E1=exp(-tesp/(2T1))和E2=exp(-tesp/(2T2))是驰豫因子。F1(n-1)可以认为是施加第n-1个射频脉冲后的预定信号(横向),其可以通过协议参数(如回波间距tesp、回波链长度(其等于n的最大取值)、T1、T2)进行调整。Z1(n-1)为施加第n-1个射频脉冲后纵向的磁化强度。F和Z的初始状态为[F1(0),F-1(0),Z1(0)]=[1,0,0],其由90°激发脉冲产生。In是第n个回波的信号强度,其等于E2F-1(n)。
如图3B中的虚线曲线所示,图3B为本发明一个例子中根据图3A所示预定的回波信号演化曲线推算得到的一种优化的可变翻转角的示意图。图3B中横坐标为回波链中的回波数Ne,纵坐标为以弧度(RAD)为单位的可变翻转角。
B、根据所述优化的可变翻转角计算得到一个计算的回波信号演化曲线。
本步骤中,将步骤B中得到的优化翻转角α作为已知量,同样可通过公式(3)推算出对应的回波信号演化曲线I,即计算的回波信号演化曲线。
实际应用中,由于计算误差的存在,经过步骤B反推后得到的回波信号演化曲线与步骤203中得到的优化的回波信号演化曲线之间会存在一些差异。
C、将所述计算的回波信号演化曲线和所述预定的回波信号演化曲线进行比较,在二者的差异满足设定要求时,将所述优化的可变翻转角用作所述实际可变翻转角链。
其中,计算的回波信号演化曲线和预定的回波信号演化曲线之间的差异越小越好。
步骤205,将所述实际可变翻转角链应用于二维快速自旋回波序列,并采用与所述第二调整参数相应的刀锋伪影校正序列来采集磁共振信号并使磁共振图像满足所述优化对比度。
以上对本发明实施例中的基于二维快速自旋回波的磁共振成像方法进行了详细描述,下面再对本发明实施例中的基于二维快速自旋回波的磁共振成像装置进行详细描述。对于本发明装置实施例中未披露的细节可参考本发明方法实施例中的相应描述。
图4为本发明实施例中一种基于二维快速自旋回波的磁共振成像装置400的结构示意图。该实施例中的装置400可用于实施图2所示实施例中的方法。如图4中的实线部分所示,该装置400可包括:第一调整参数确定模块410、第二调整参数确定模块420、第一回波信号演化曲线确定模块430、第一可变翻转角链确定模块440和磁共振成像模块450。
其中,第一调整参数确定模块410适用于确定用于预设一磁共振图像的初始对比度的第一调整参数。
例如,在一个实施方式中,用户可根据所期望的磁共振图像组织对比度,对已有的或新增的成像参数进行调整。例如,可预先设置不同的可变翻转角模式,如T1权重的可变翻转角模式、PD权重的可变翻转角模式和T2权重的可变翻转角模式。其中,所述T1权重的可变翻转角模式用于使回波信号演化曲线符合T1权重的磁共振成像,所述PD权重的可变翻转角模式用于使回波信号演化曲线符合PD权重的磁共振成像,所述T2权重的可变翻转角模式用于使回波信号演化曲线符合T2权重的磁共振成像。相应地,第一调整参数确定模块410可根据用户或系统的选择,确定第一调整参数为相应权重的可变翻转角模式。
第二调整参数确定模块420适用于确定用于得到磁共振图像的优化对比度及该刀锋伪影校正序列的设定数据采集时间的第二调整参数。例如,在一个实施方式中,第二调整参数确定模块420可确定用于得到磁共振图像的优化对比度的回波间距和/或回波链长度,以及用于得到该刀锋伪影校正序列的设定数据采集时的回波链长度。
回波信号演化曲线确定模块430根据所述第一调整参数和所述第二调整参数来确定一优化的回波信号演化曲线。
可变翻转角链确定模块440根据所述优化的回波信号演化曲线计算得到实际可变翻转角链。
磁共振成像模块450将所述实际可变翻转角应用于二维快速自旋回波序列,并采用与所述第二调整参数相应的刀锋伪影校正序列来采集磁共振信号并使磁共振图像满足所述优化对比度。
在一个实施方式中,在第一调整参数确定模块410确定第一调整参数之后,回波信号演化曲线确定模块430可进一步用于根据所述第一调整参数来预定一初始回波信号演化曲线;相应地,可变翻转角链确定模块440可进一步用于根据所述初始回波信号演化曲线计算得到初始可变翻转角链。
在一个实施方式中,所述可变翻转角链确定模块440可包括:第一单元441、第二单元442和第三单元443。
第一单元441用于根据所述优化的回波信号演化曲线推算得到优化的可变翻转角。
第二单元442根据所述优化的可变翻转角计算得到一个计算的回波信号演化曲线。
第三单元443将所述计算的回波信号演化曲线和所述优化的回波信号演化曲线进行比较,在二者的差异满足设定要求时,将所述优化的可变翻转角用作所述实际可变翻转角链。
图5A至图5D示出了本发明一个例子中一组采用传统BLADE序列和采用本发明实施例中的带有翻转角的BLADE序列进行成像时获取的脑部图像的对比图。
其中,图5A为采用传统BLADE序列进行成像时获取的脑部轴向图像,图5B为采用本发明实施例中的带有翻转角的BLADE序列进行成像时获取的脑部轴向图像,可见,图5B所示图像的图像对比度和图像质量与图5A所示图像的图像对比度和图像质量相当,但由于增大了叶片宽度使得扫描速度提高了52%。
图5C为采用传统BLADE黑水序列进行成像时获取的脑部轴向图像,图5D为采用本发明实施例中的带有翻转角的BLADE黑水序列进行成像时获取的脑部轴向图像,可见,图5D所示图像的图像对比度和图像质量与图5C所示图像的图像对比度和图像质量相当,但由于增大了叶片宽度使得扫描速度提高了61%。
图6A至图6D示出了本发明一个例子中一组采用传统BLADE序列和采用本发明实施例中的带有翻转角的BLADE序列进行成像时获取的膝盖图像的对比图。
其中,图6A为采用传统BLADE序列进行成像时获取的膝盖轴向图像,图6B为采用本发明实施例中的带有翻转角的BLADE序列进行成像时获取的膝盖轴向图像,可见,图6B所示图像的图像对比度和图像质量与图6A所示图像的图像对比度和图像质量相当,但由于增大了叶片宽度使得扫描速度提高了42%。
图6C为采用传统BLADE序列进行成像时获取的膝盖冠状面图像,图6D为采用本发明实施例中的带有翻转角的BLADE序列进行成像时获取的膝盖冠状面图像,可见,
图6D所示图像的图像对比度和图像质量与图6C所示图像的图像对比度和图像质量相当,但由于增大了叶片宽度使得扫描速度提高了45%。
图7A和图7B示出了本发明一个例子中一组采用传统BLADE序列和采用本发明实施例中的带有翻转角的BLADE序列进行成像时获取的手腕图像的对比图。
其中,图7A为采用传统BLADE序列进行成像时获取的手腕冠状面图像,图7B为采用本发明实施例中的带有翻转角的BLADE序列进行成像时获取的手腕冠状面图像,可见,图7B所示图像的图像对比度和图像质量与图7A所示图像的图像对比度和图像质量相当,但由于增大了叶片宽度使得扫描速度提高了39%。
从上述方案中可以看出,由于本发明实施例中首先确定一优化的回波信号演化曲线,基于该回波信号演化曲线能够使得磁共振图像满足预设的对比度,并使得BLADE序列的采集时间满足设定的短时要求,之后基于该优化的回波信号演化曲线计算得到可变翻转角链,并将该可变翻转角链应用于二维快速自旋回波序列,并采用与所述第二调整参数相应的刀锋伪影校正序列来采集磁共振信号。从而可在增加叶片宽度时,通过调整翻转角的值来控制T2的衰减,进而可以在保证图像的成像质量和对比度的前提下,可通过增加叶片宽度来减少扫描时间,还可以获得更多的k空间中心的公共数据,从而可更好的进行运动校正,减少伪影。
进一步地,通过设置针对可变翻转角模式的成像参数,从而可以通过对成像参数进行调整实现满足所期望的组织对比度的图像。
此外,本发明实施例中给出了一种优化可变翻转角链的简单可靠的实现方式。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种基于二维快速自旋回波的磁共振成像方法,其利用刀锋伪影校正序列进行磁共振信号的数据采集,其特征在于,包括:
确定用于预设一磁共振图像的初始对比度的第一调整参数;
确定用于得到磁共振图像的优化对比度及该刀锋伪影校正序列的设定数据采集时间的第二调整参数;
根据所述第一调整参数和所述第二调整参数来确定一优化的回波信号演化曲线;
根据所述优化的回波信号演化曲线计算得到实际可变翻转角链;
将所述实际可变翻转角链应用于二维快速自旋回波序列,并采用与所述第二调整参数相应的刀锋伪影校正序列来采集磁共振信号并使磁共振图像满足所述优化对比度。
2.根据权利要求1所述的基于二维快速自旋回波的磁共振成像方法,其特征在于,所述确定用于预设一磁共振图像的初始对比度的第一调整参数之后,还包括:
根据所述第一调整参数来预定一初始回波信号演化曲线;
根据所述初始回波信号演化曲线计算得到初始可变翻转角链。
3.根据权利要求1所述的基于二维快速自旋回波的磁共振成像方法,其特征在于,所述第二调整参数选自回波间距、回波链长度中的至少一个。
4.根据权利要求1所述的基于二维快速自旋回波的磁共振成像方法,其特征在于,所述根据所述优化的回波信号演化曲线计算得到实际可变翻转角链包括:
根据所述优化回波信号演化曲线推算得到优化的可变翻转角;
根据所述优化的可变翻转角计算得到一个计算的回波信号演化曲线;
将所述计算的回波信号演化曲线和所述优化的回波信号演化曲线进行比较,在二者的差异满足设定要求时,将所述优化的可变翻转角用作所述实际可变翻转角链。
5.根据权利要求1至4中任一项所述的基于二维快速自旋回波的磁共振成像方法,其特征在于,所述第一调整参数为预设的如下三种可变翻转角模式中的一种:
纵向磁化矢量驰豫时间T1权重、质子密度PD权重和横向磁化矢量驰豫时间T2权重。
6.一种基于二维快速自旋回波的磁共振成像装置,其利用刀锋伪影校正序列进行磁共振信号的数据采集,其特征在于,包括:
第一调整参数确定模块,适用于确定用于预设一磁共振图像的初始对比度的第一调整参数;
第二调整参数确定模块,适用于确定用于得到磁共振图像的优化对比度及该刀锋伪影校正序列的设定数据采集时间的第二调整参数;
回波信号演化曲线确定模块,根据所述第一调整参数和所述第二调整参数来确定一优化的回波信号演化曲线;
可变翻转角链确定模块,根据所述优化的回波信号演化曲线计算得到实际可变翻转角链;
磁共振成像模块,将所述实际可变翻转角链应用于二维快速自旋回波序列,并采用与所述第二调整参数相应的刀锋伪影校正序列来采集磁共振信号并使磁共振图像满足所述优化对比度。
7.根据权利要求6所述的基于二维快速自旋回波的磁共振成像装置,其特征在于,
所述回波信号演化曲线确定模块进一步用于根据所述第一调整参数来预定一初始回波信号演化曲线;
所述可变翻转角链确定模块进一步用于根据所述初始回波信号演化曲线计算得到初始可变翻转角链。
8.根据权利要求6所述的基于二维快速自旋回波的磁共振成像装置,其特征在于,所述第二调整参数确定模块确定用于得到磁共振图像的优化对比度的回波间距和/或回波链长度,以及用于得到该刀锋伪影校正序列的设定数据采集时间的回波链长度。
9.根据权利要求6所述的基于二维快速自旋回波的磁共振成像装置,其特征在于,所述可变翻转角链确定模块包括:
第一单元,用于根据所述优化的回波信号演化曲线推算得到优化的可变翻转角;
第二单元,根据所述优化的可变翻转角计算得到一个计算的回波信号演化曲线;
第三单元,将所述计算的回波信号演化曲线和所述优化的回波信号演化曲线进行比较,在二者的差异满足设定要求时,将所述优化的可变翻转角用作所述实际可变翻转角链。
10.根据权利要求6至9中任一项所述的基于二维快速自旋回波的磁共振成像装置,其特征在于,所述第一调整参数确定模块确定用于预设一磁共振图像的初始对比度的可变翻转角模式,其为预设的如下三种可变翻转角模式中的一种:
纵向磁化矢量驰豫时间T1权重、质子密度PD权重和横向磁化矢量驰豫时间T2权重。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811618720.9A CN111381204B (zh) | 2018-12-28 | 2018-12-28 | 基于二维快速自旋回波的磁共振成像方法和装置 |
US16/729,799 US20200205693A1 (en) | 2018-12-28 | 2019-12-30 | Magnetic resonance imaging method and apparatus based on two-dimensional fast spin echo |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811618720.9A CN111381204B (zh) | 2018-12-28 | 2018-12-28 | 基于二维快速自旋回波的磁共振成像方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111381204A CN111381204A (zh) | 2020-07-07 |
CN111381204B true CN111381204B (zh) | 2022-07-12 |
Family
ID=71123718
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811618720.9A Active CN111381204B (zh) | 2018-12-28 | 2018-12-28 | 基于二维快速自旋回波的磁共振成像方法和装置 |
Country Status (2)
Country | Link |
---|---|
US (1) | US20200205693A1 (zh) |
CN (1) | CN111381204B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112731235B (zh) * | 2020-12-09 | 2024-04-19 | 中国科学院深圳先进技术研究院 | 一种磁共振化学交换饱和转移成像方法以及相关设备 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102232831A (zh) * | 2010-04-30 | 2011-11-09 | 西门子(深圳)磁共振有限公司 | 一种实现水脂分离的磁共振成像方法 |
CN103118593A (zh) * | 2011-09-22 | 2013-05-22 | 株式会社东芝 | 磁共振成像装置以及磁共振成像方法 |
CN103505210A (zh) * | 2012-06-28 | 2014-01-15 | 西门子(深圳)磁共振有限公司 | 一种实现水脂分离的磁共振成像方法和装置 |
JP2016086991A (ja) * | 2014-10-31 | 2016-05-23 | 株式会社東芝 | 磁気共鳴イメージング装置 |
CN108983133A (zh) * | 2018-07-23 | 2018-12-11 | 中国科学院深圳先进技术研究院 | 磁共振血管壁成像方法、装置、设备及存储介质 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9360545B2 (en) * | 2012-06-26 | 2016-06-07 | Siemens Aktiengesellschaft | Magnetic resonance system and operating method for flow artifact reduction in slab selective space imaging |
WO2017173617A1 (en) * | 2016-04-07 | 2017-10-12 | Shanghai United Imaging Healthcare Co., Ltd. | System and method for fast imaging in magnetic resonance imaging |
EP3607340B1 (en) * | 2017-04-05 | 2023-07-26 | Beth Israel Deaconess Medical Center, Inc. | System and method for improved spin-echo-based magnetic resonance imaging |
CN107638179B (zh) * | 2017-08-28 | 2020-09-04 | 上海联影医疗科技有限公司 | 磁共振翻转角的计算方法、系统和计算机可读存储介质 |
-
2018
- 2018-12-28 CN CN201811618720.9A patent/CN111381204B/zh active Active
-
2019
- 2019-12-30 US US16/729,799 patent/US20200205693A1/en active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102232831A (zh) * | 2010-04-30 | 2011-11-09 | 西门子(深圳)磁共振有限公司 | 一种实现水脂分离的磁共振成像方法 |
CN103118593A (zh) * | 2011-09-22 | 2013-05-22 | 株式会社东芝 | 磁共振成像装置以及磁共振成像方法 |
CN108318843A (zh) * | 2011-09-22 | 2018-07-24 | 东芝医疗系统株式会社 | 磁共振成像装置以及磁共振成像方法 |
CN103505210A (zh) * | 2012-06-28 | 2014-01-15 | 西门子(深圳)磁共振有限公司 | 一种实现水脂分离的磁共振成像方法和装置 |
JP2016086991A (ja) * | 2014-10-31 | 2016-05-23 | 株式会社東芝 | 磁気共鳴イメージング装置 |
CN108983133A (zh) * | 2018-07-23 | 2018-12-11 | 中国科学院深圳先进技术研究院 | 磁共振血管壁成像方法、装置、设备及存储介质 |
Non-Patent Citations (3)
Title |
---|
BLADE刀锋技术详解;翁得河;《磁共振成像》;20100920;第1卷(第05期);第376-379页 * |
BLADE和SSTSE T2WI在头部MR成像中的价值;梁汉欢等;《中国现代医生》;20090818;第47卷(第23期);第105-106页 * |
螺旋桨技术在头部磁共振扫描中应用优势及局限性分析;胡君;《中国民康医学》;20110325;第23卷(第06期);第771-772页 * |
Also Published As
Publication number | Publication date |
---|---|
US20200205693A1 (en) | 2020-07-02 |
CN111381204A (zh) | 2020-07-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10874326B2 (en) | System and method for fast imaging in magnetic resonance imaging | |
JP6513398B2 (ja) | 事前情報に制約される正則化を使用するmr画像再構成 | |
JP4130405B2 (ja) | 磁気共鳴撮影装置 | |
US10262385B2 (en) | Method and apparatus for acquiring magnetic resonance data and entering the data into k-space | |
US8587310B2 (en) | Magnetic resonance imaging device | |
US8648593B2 (en) | Dark blood balanced steady state free precession imaging | |
JPH07323016A (ja) | 磁気共鳴イメージング方法および装置 | |
US10429463B2 (en) | Quiet MRI with spin echo (SE) or fast spin echo (FSE) | |
US10761165B2 (en) | System and method for magnetic resonance image acquisition | |
JP2003225223A (ja) | 磁気共鳴イメージング装置 | |
US20160033608A1 (en) | Magnetic resonance angiography method and apparatus | |
CN111381204B (zh) | 基于二维快速自旋回波的磁共振成像方法和装置 | |
US11022667B2 (en) | System and method for image reconstruction | |
US10156622B2 (en) | Method and apparatus for sectional optimization of radial MR pulse sequences | |
CN116466279A (zh) | 一种用于肺部的磁共振成像方法、系统、电子设备及介质 | |
KR101844514B1 (ko) | 자기 공명 영상 장치 및 자기 공명 영상 획득 방법 | |
US10234529B2 (en) | Method for producing magnet resonance tomography recordings of cyclic movement | |
CN108577841B (zh) | 一种propeller技术中抑制非刚性运动的权重计算方法 | |
EP3451009B1 (en) | System and method for determining flip angles of refocusing radio frequency pulses in magnetic resonance imaging | |
JP6788510B2 (ja) | 磁気共鳴イメージング装置 | |
CN107843862B (zh) | 一种propeller技术中参考位置图像的非迭代生成方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |