CN103654779B - 磁共振成像方法和装置 - Google Patents
磁共振成像方法和装置 Download PDFInfo
- Publication number
- CN103654779B CN103654779B CN201210362793.2A CN201210362793A CN103654779B CN 103654779 B CN103654779 B CN 103654779B CN 201210362793 A CN201210362793 A CN 201210362793A CN 103654779 B CN103654779 B CN 103654779B
- Authority
- CN
- China
- Prior art keywords
- pulse
- spectrum
- reversing time
- flip angle
- suppressor
- 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
Classifications
-
- 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/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/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
- G01R33/56527—Correction of image distortions, e.g. due to magnetic field inhomogeneities due to chemical shift effects
-
- 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/5607—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by reducing the NMR signal of a particular spin species, e.g. of a chemical species for fat suppression, or of a moving spin species for black-blood 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/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/5613—Generating steady state signals, e.g. low flip angle sequences [FLASH]
-
- 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/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
- G01R33/56509—Correction of image distortions, e.g. due to magnetic field inhomogeneities due to motion, displacement or flow, e.g. gradient moment nulling
Landscapes
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- General Health & Medical Sciences (AREA)
- Radiology & Medical Imaging (AREA)
- Engineering & Computer Science (AREA)
- High Energy & Nuclear Physics (AREA)
- Signal Processing (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
公开了磁共振成像方法和装置。所述磁共振成像方法包括步骤:根据谱抑制脉冲所抑制的谱成分的纵向磁化分量的稳定状态条件和所述纵向磁化分量的过零点条件,计算出该谱抑制脉冲的翻转角和/或反转时间;应用包括设有所述翻转角和/或所述反转时间的所述谱抑制脉冲的磁共振成像序列进行磁共振成像。
Description
技术领域
本发明涉及磁共振成像方法和装置,特别涉及包括谱抑制脉冲的磁共振成像方法和装置。
背景技术
磁共振成像是随着计算机技术、电子电路技术、超导体技术的发展而迅速发展起来的一种生物磁学核自旋成像技术。在磁共振成像中,人体组织被置于静磁场B0中,随后用频率与氢原子核的进动频率相同的射频脉冲激发人体组织内的氢原子核,引起氢原子核共振,并吸收能量;在停止射频脉冲后,氢原子核按特定频率发出射电信号,并将吸收的能量释放出来,由体外的接受器收录,经计算机处理后获得图像。
在磁共振成像中,为了获得较好的成像质量,往往需要对例如脂肪信号、水信号、硅胶信号(乳房植入物)等特定谱成分的信号进行抑制。例如,在用于腹部和胸部等的磁共振成像检查中,普遍需要抑制脂肪信号,以突出感兴趣组织或病灶的显示。至今为止,已经提出了许多脂肪抑制技术,例如FatSat(fatsaturation,脂肪饱和)技术,SPAIR(SpectralPresaturationAttenuatedInversionRecovery,频谱预饱和衰减反转恢复)技术,STIR(shortinversiontimeinversionrecovery,短反转时间反转恢复)技术,DIXON水脂分离技术等。在各种脂肪抑制技术中,除了DIXON技术,其他技术的效能均取决于向人体组织施加的脂肪抑制脉冲的翻转角(例如对于FatSat)或反转时间(例如对于SPAIR和STIR等反转恢复技术)。众所周知,翻转角是指在射频脉冲的激发下,宏观磁化矢量偏离静磁场B0的角度。该翻转角的大小是由施加的射频脉冲的强度和作用时间决定的。反转时间是反转恢复脉冲序列中180度反转脉冲与激发脉冲之间的时间间隔。
如前所述,通过优化脂肪抑制脉冲的翻转角或反转时间可以实现脂肪抑制的优化。针对这一情况,目前已经提出了一种用于计算FatSat脉冲的翻转角的方法。但是该方法仅仅考虑了采用自旋回波序列进行磁共振成像时FatSat脉冲的翻转角的计算方法,而没有考虑在应用短射频脉冲间隔(TR)的损毁梯度回波序列(SpoiledGRE,例如VIBE序列)的腹部磁共振成像中,在每个脂肪抑制脉冲之后施加的多个激发脉冲。另外,该方法也无法用于优化脂肪抑制脉冲的反转时间(例如对于SPAIR和STIR等反转恢复脉冲)。
发明内容
有鉴于此,本发明提出了一种磁共振成像方法,包括步骤:根据谱抑制脉冲所抑制的谱成分的纵向磁化分量的稳定状态条件和所述纵向磁化分量的过零点条件,计算出该谱抑制脉冲的翻转角和/或反转时间;应用包括设有所述翻转角和/或所述反转时间的所述谱抑制脉冲的磁共振成像序列进行磁共振成像。
本发明还提出了一种磁共振成像装置,包括:一计算部件,根据谱抑制脉冲所抑制的谱成分的纵向磁化分量的稳定状态条件和优化的谱抑制条件,计算出该谱抑制脉冲的翻转角和/或反转时间;一成像部件,应用包括设有所述翻转角和/或所述反转时间的所述谱抑制脉冲的磁共振成像序列进行磁共振成像。
所述谱抑制脉冲可以是例如包括脂肪饱和脉冲和水饱和脉冲的谱饱和脉冲,包括谱预饱和衰减反转恢复脉冲、短反转时间反转恢复脉冲和体衰减反转恢复脉冲的反转恢复脉冲等各种谱抑制脉冲。
从上述方案中可以看出,由于本发明的磁共振成像方法和装置中使用具有优化的翻转角和/或反转时间的谱抑制脉冲以尽可能地抑制相应的谱成分信号,因此可以获得高质量的图像。
附图说明
下面将通过参照附图详细描述本发明的优选实施例,使本领域的普通技术人员更清楚本发明的上述及其它特征和优点,附图中:
图1为设置有脂肪抑制脉冲的损毁梯度回波序列的示意图。
图2示出了在连续的脂肪抑制脉冲之间施加多个激发脉冲的序列图。
图3示出了根据本发明实施例的使用设置有谱抑制脉冲的MRI序列进行磁共振成像的方法的流程图。
图4示出了使用设置有谱抑制脉冲的MRI序列进行磁共振成像的装置的框图。
图5(a)和6(b)示出了当谱抑制脉冲为FatSat脉冲时的仿真结果。
图6(a)和7(b)示出了当谱抑制脉冲为SPAIR脉冲时的仿真结果。
图7(a)和7(b)示出了根据本发明的、利用设置有优化的谱抑制脉冲的MRI序列采集的图像;图7(c)示出了使用未配置谱抑制脉冲的MRI序列采集的图像。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,以下举实施例对本发明进一步详细说明。
(实施例1)
本实施例中将脂肪抑制脉冲作为具体的谱抑制脉冲来进行描述。
首先对本发明中涉及的MRI序列进行简要描述。
图1为设置有脂肪抑制脉冲的损毁梯度回波序列的示意图。由于损毁梯度回波序列的射频脉冲间隔TR较短,并且在连续的脂肪抑制脉冲之间可以施加多个激发脉冲,因此减少了测量时间。因而,该序列适合于在需要短采样时间以避免运动伪影的腹部和胸部等的MRI检查中使用。
如图1所示,首先施加脂肪抑制脉冲,并且由于施加该脂肪抑制脉冲而在XY平面中产生的横向磁化矢量通过损毁梯度而被消除。随后施加构成梯度回波序列的射频激发脉冲,选层梯度,相位编码梯度等等(为了清楚起见,图中仅示出一个激发脉冲),并在施加激发脉冲之后,产生如图1中所示的梯度回波。
图2图示了在连续的脂肪抑制脉冲之间施加多个激发脉冲的序列图。为了简化附图,图2中未示出梯度脉冲(选择性梯度、相编码梯度、损毁梯度等)和射频损毁机制。
从图2中可以看出,序列中的一个脉冲串由一个脂肪抑制脉冲和多个射频激发脉冲构成。所述脂肪抑制脉冲是例如前述的FatSat脉冲、SPAIR脉冲、STIR脉冲等。图中所示的β是该脂肪抑制脉冲的翻转角,α是激发脉冲的翻转角。
如图2所示,TI1表示脂肪抑制脉冲的中心和脂肪抑制脉冲之后的第一激发脉冲的中心之间的时间间隔,Tr表示磁共振成像序列的一个脉冲串中最后一个激发脉冲的中心与下一脉冲串中的脂肪抑制脉冲的中心之间的时间间隔,Trep表示连续的激发脉冲之间的时间。此外,在本发明中,用TI1Min表示TI1的最小值,TIFill表示第一反转时间。容易理解,TI1Min+TIFill=TI1。类似的,在本发明中,用TrMin表示Tr的最小值,TrFill表示第二反转时间,并且TrMin+TrFill=Tr。
下面,对本发明中涉及的K空间和K空间中心线(KSpaceCenterLine)进行简要描述。
K空间也称傅里叶空间,是磁共振信号原始数据的填充空间,而这些磁共振信号是在脉冲序列运行时采集的,在进行傅立叶变换后,就能变成图像。在二维图像的信号采集过程中,每个磁共振信号的频率编码梯度场的大小和方向保持不变,而相位编码梯度的强度则以一定的步阶发生变化,每个磁共振信号的相位编码变化一次,采集到的磁共振信号填充K空间Ky方向的一条线。在本发明中,多个激发脉冲中的每一个激发脉冲对应于K中间中的一条线。K空间的不同位置的数据对最终图像的贡献是不同的,K空间中心的数据主要决定图像的信噪比和对比度信息,K空间的边缘部分的信号主要贡献图像的分辨能力方面的信息。
以下将描述根据本实施例的磁共振成像方法中采用的优化的脂肪抑制脉冲。如前所述,脂肪抑制的能效很大程度上取决于脂肪抑制脉冲的翻转角和反转时间,因此通过优化脂肪抑制脉冲的翻转角或反转时间可以实现脂肪抑制的优化。因此,如何确定优化的翻转角和反转时间成为一个关键。对此,本实施例提出了一种通用的解决方案,其既可以用于计算脂肪抑制脉冲的翻转角(对于例如FatSat)也可以用于计算脂肪抑制脉冲的反转时间(对于包括SPAIR和STIR在内的反转脉冲),同时考虑了在损毁梯度回波序列的每个脂肪抑制脉冲后跟随的多个激发脉冲的情况。下面首先对该通用解决方案进行描述。
1.从布洛赫方程推导出脂肪成分的磁化矢量表达式
假设损毁梯度回波序列的损毁机制理想地工作,从而在XY平面的横向磁化分量在每个激发脉冲(α脉冲)之前完全损毁并且对于随后的激发脉冲没有贡献。还假设脂肪抑制脉冲(β脉冲)仅激发脂肪成分。在这样的条件下,根据经典布洛赫方程,每个脉冲之前的脂肪成分的纵向磁化矢量表示可以如下所示:
Mβ(1)=M0(1)
Mα(1)=M0(1-ETI1)+Mβ(m)·cos(β)·ETI1(2)
Mα(n)=M0(1-ETrep)+Mα(n-1)·cos(α)·ETrep(3)
Mβ(m+1)=M0(1-ETr)+Mα(N)·cos(α)·ETr(4)
其中:
ETI1=exp(-TI1/T1)
ETrep=exp(-Trep/T1)
ETr=exp(-Tr/T1)
T1:脂肪组织的纵向弛豫时间,
M0:脂肪组织的初始磁化矢量,
Mα(n):第n个α脉冲之前的脂肪的纵向磁化分量,
Mβ(m):第m个β脉冲之前的脂肪的纵向磁化分量,
N:两个β脉冲之间的α脉冲的数目,即每个脉冲串填充K空间的线数。
2.建立通式
由于在计算期间,M0将被抵消掉,因此可以设置为1,即M0=1。
因此,等式(1)到(4)可以写成
Mβ(1)=1(5)
Mα(1)=A+B·Mβ(m)(6)
Mα(n)=C+D·Mα(n-1)(7)
Mβ(m+1)=E+F·Mα(N)(8)
其中
A=1-ETI1B=cos(β)·ETI1
C=1-ETrepD=cos(α)·ETrep
E=1-ETrF=cos(α)·ETr
根据等式(6)和(7),可以得到
Mα(1)=A+B·Mβ(m)
Mα(2)=C+D·Mα(1)
=C+A·D+B·D·Mβ(m)
Mα(3)=C+D·Mα(2)
=C+C·D+A·D2+B·D2·Mβ(m)
Mα(4)=C+D·Mα(3)
=C+C·D+C·D2+A·D3+B·D3·Mβ(m)
……
Mα(n)=C+C·D+C·D2+C·D3…+C·Dn-2+A·Dn-1+B·Dn-1·Mβ(m)
=C(Dn-1-1)/(D-1)+A·Dn-1+B·Dn-1·Mβ(m)(9)
(其中,n=1,2,3,…,并且D≠1因为cos(α)·ETrep≠1)
根据等式(8)和(9),可以得到
Mβ(m+1)=E+F·Mα(N)
=E+F·{C·(DN-1-1)/(D-1)+A·DN-1+B·DN-1·Mβ(m)}
=E+F·C·(DN-1-1)/(D-1)+F·A·DN-1+F·B·DN-1·Mβ(m)(10)
针对上面的等式,在本发明中通过引入稳定状态条件和纵向磁化分量的过零点条件来进行推导。
稳定状态条件是本领域中公知的。具体的,在稳定状态时,每个β脉冲之前的脂肪的纵向磁化分量应当是相同的,因此稳定状态条件为:
Mβ(m+1)=Mβ(m)(11)。
另外,基于前述K空间和K空间填充线的特性,认识到当对应于K空间中一条或多条填充线的激发脉冲之前的脂肪的纵向磁化分量为零时,能够获得优化的脂肪抑制。特别的,当对应于K空间中心线的激发脉冲之前的脂肪的纵向磁化分量为零时,具有最优的脂肪抑制。因此可以设定纵向磁化分量的过零点条件为:
Ma(KSpaceCenterLine)=0(12),
其中,KSpaceCenterLine是K空间中心线的索引编号。
根据稳定状态条件(11)和等式(10),可以得到
Mβ(m)=E+F·C·(DN-1-1)/(D-1)+F·A·DN-1+F·B·DN-1·Mβ(m),
从而
Mβ(m)={E+F·C·(DN-1-1)/(D-1)+F·A·DN-1}/(1-F·B·DN-1)(13)
另外,结合纵向磁化分量的过零点条件(12)和等式(9),可以得到
0=C·(DKSpaceCenterLine-1-1)/(D-1)+A·DKSpaceCenterLine-1+B·DKSpaceCenterLine-1·Mβ(m)
从而
Mβ(m)={C·(DKSpaceCenterLine-1-1)/(D-1)+A·DKSpaceCenterLine-1}/(-B·DKSpaceCenterLine-1)
(14)
根据等式(13)和(14),可以得到
{E+F·C·(DN-1-1)/(D-1)+F·A·DN-1}/(1-F·B·DN-1)
={C·(DKSpaceCenterLine-1-1)/(D-1)+A·DKSpaceCenterLine-1}/(-B·DKSpaceCenterLine-1)
(15)
由此得到最终的通式(15)。通过该通式,本领域技术人员可以容易地计算出翻转角β以及反转时间TIFill和TrFill。尽管如此,为了清楚起见,下文中以FatSat脉冲的翻转角β和SPAIR脉冲的反转时间TIFill和TrFill为例,给出了计算翻转角和反转时间的示例。
FatSat脉冲的翻转角β:
对于设置有FatSat脉冲的损毁梯度回波序列,以下参数是可以根据进行磁共振成像的硬件的条件以及用户的要求等预先确定的:Trep,α,KSpaceCenterLine,N,Tr=TrMin,TI1=TI1Min。
根据等式(15),可以得到
β=arccos(X/ETI1)
其中
X=G1/(G1·G4–G2·G3)
G1=C·(DKSpaceCenterLine-1-1)/(D-1)+A·DKSpaceCenterLine-1
G2=DKSpaceCenterLine-1
G3=E+F·C·(DN-1-1)/(D-1)+F·A·DN-1
G4=F·DN-1
因此,可以容易的计算出β。
SPAIR脉冲的第一反转时间TIFill和第二反转时间TrFill:
对于设置有SPAIR脉冲的损毁梯度回波序列,以下参数是可以根据进行磁共振成像的硬件的条件以及用户的要求等预先确定的:Trep,α,β=180°,KspaceCenterLine,N,TrMin,TI1Min。并且如前所述,TI1=TIFill+TI1Min,Tr=TrFill+TrMin。
TrFill和TIFill对脂肪的纵向磁化分量为零的位置的影响是相反的,即较长的TrFill将增大脂肪信号为零的线索引,而较长的TIFill将减小脂肪信号为零的线索引。并且对于任何参数集,反转时间的结果将是(TIFill>=0,TrFill=0)或(TIFill=0,TrFill>=0)。
所以对于SPAIR脉冲的第一反转时间TIFill和第二反转时间TrFill,通过以下步骤进行计算:
步骤A:保持TrFill=0,即Tr=TrMin,计算TIFill。
由于根据等式(15),可以得到
TIFill=TI1-TI1Min=-T1·ln(X/cos(β))-TI1Min
其中
X=(G1+G2)/{G1·G4+G2/cos(β)–G2·G3}
G1=C·(DKSpaceCenterLine-1-1)/(D-1)
G2=DKSpaceCenterLine-1
G3=E+F·C·(DN-1-1)/(D-1)
G4=F·DN-1
因此,可以容易的计算出TIFill。如果TIFill>=0,则意味着对于目前的参数集已经计算出了反转时间,其中TIFill>=0,TrFill=0;而如果TIFill<0,则意味着当TrFill=0时计算得到的结果集(TIFill,TrFill=0)不符合要求,因此需要执行步骤B。
步骤B:保持TIFill=0,即TI1=TIMin,计算TrFill。
由于根据等式(15),可以得到
TrFill=Tr-TrMin=-T1·ln(Y/cos(β))–TrMin
其中
Y=(G1+G2)/(G1·G4-G2·G3)
G1=C·(DKSpaceCenterLine-1-1)/(D-1)+A·DKSpaceCenterLine-1
G2=B·DKSpaceCenterLine-1
G3=C·(DN-1-1)/(D-1)+A·DN-1-1/cos(α)
G4=B·DN-1
因此,可以容易的计算出TrFill。如果TrFill>=0,则意味着对于目前的参数集已经计算出了反转时间,其中TIFill=0,TrFill>=0;而如果TIFill<0,则意味着对于目前的参数集没有对应的反转时间,结果将是TrFill=0,TIFill=0。
需要说明的是,虽然在上面的示例中,在步骤A中保持TrFill=0,计算出对应的TIFill,随后在步骤B中保持TIFill,计算出对应的TrFill,但这仅仅是一种示例。在步骤A中保持TIFill=0,计算出对应的TrFill,随后如果需要则在步骤B中保持TrFill,计算出对应的TrFill也是可以的。
以上已经详细描述了用于计算优化的翻转角和反转时间的通用解决方案。下面,将参考图3对根据本实施例的磁共振成像方法进行描述。
如图3所示,在步骤301中,根据脂肪抑制脉冲所抑制的脂肪成分的纵向磁化分量的稳定状态条件和所述纵向磁化分量的过零点条件,计算出该谱抑制脉冲的翻转角或反转时间。在该步骤中,预先确定脂肪抑制脉冲的各个参数。所述参数包括例如两个脂肪抑制脉冲之间的激发脉冲的数目N、脂肪抑制脉冲的中心和该脂肪抑制脉冲之后的第一激发脉冲的中心之间的最小时间间隔TI1Min、MRI序列的一个脉冲串中的最后一个激发脉冲的中心与下一个脉冲串的脂肪抑制脉冲的中心之间的时间间隔TrMin、连续的激发脉冲之间的时间Trep、激发脉冲的翻转角α、K空间中心线的索引编号KSpaceCenterLine等。随后,根据脂肪的纵向磁化分量的稳定状态条件和纵向磁化分量的过零点条件,基于所述各个参数计算出该脂肪抑制脉冲的翻转角或反转时间。具体的,根据布洛赫方程得出脂肪的磁化矢量表示,并应用脂肪的纵向磁化分量的稳定状态条件和纵向磁化分量的过零点条件得出用于计算脂肪抑制脉冲的翻转角或反转时间的通式,随后根据已确定的各个参数计算出所述翻转角或反转时间。应用脂肪的纵向磁化分量的稳定状态条件Mβ(m+1)=Mβ(m)和纵向磁化分量的过零点条件Ma(KSpaceCenterLine)=0是得到优化的翻转角和反转时间的一个关键所在。随后,在步骤302中,应用包括设有所述翻转角或所述反转时间的脂肪抑制脉冲的MRI序列进行磁共振成像。在MRI序列中应用脂肪抑制脉冲的一种示例方式是使得在每一个脂肪抑制脉冲后跟随一个或多个激发脉冲。MRI序列可以是任意序列,只要其能够在施加下一个激发脉冲之前使得残留的横向磁化分量被损毁即可。例如,该MRI序列可以是损毁快速自旋回波(TSE)序列、损毁梯度回波(GRE)序列等。
(实施例2)
上述实施例1中描述了针对脂肪抑制脉冲的优化方法,该优化方法也可以应用于水抑制脉冲。
在磁共振成像中,除了脂肪抑制之外,可能需要通过施加例如水饱和脉冲、FLAIR(体衰减反转恢复)脉冲来对人体组织的水信号进行抑制。针对该水饱和脉冲和FLAIR脉冲,均可以适用本发明实施例1中给出的通式(15)。更明确地说,对于水饱和脉冲,可以通过与实施例1中描述的用于FatSat的解决方案相同的方法计算出优化翻转角,从而得到优化的水抑制脉冲,并应用实施例1中参照图3描述的磁共振成像方法来成像;同样,对于FLAIR,可以通过与实施例1中描述的用于SPAIR的解决方案相同的方法计算出优化的TIfill和Trfill值,从而得到优化的水抑制脉冲,并应用实施例1中参照图3描述的磁共振成像方法来优化水抑制脉冲。
事实上,除了脂肪抑制脉冲和水抑制脉冲之外,本发明实施例1中描述的优化方法可以应用于任何特定谱抑制脉冲,仅仅需要改变例如T1或者抑制脉冲的中心频率等参数即可。也就是说,本发明实施例1中描述的磁共振成像方法实际上可以应用任意谱抑制脉冲。
以上已经对根据本发明的磁共振成像方法进行了详细描述。在该磁共振成像方法中,通过利用设置有优化的谱抑制脉冲的MRI序列进行成像,可以很好地抑制特定谱成分的信号,从而可以获得高质量的图像。
图4示出了使用设置有谱抑制脉冲的MRI序列进行磁共振成像的装置400。如图4所示,该磁共振成像设备包括:计算部件401,根据谱抑制脉冲所抑制的谱成分的纵向磁化分量的稳定状态条件和优化的谱抑制条件,计算出该谱抑制脉冲的翻转角或反转时间;成像部件402,应用包括设有所述翻转角或所述反转时间的所述谱抑制脉冲的磁共振成像序列进行磁共振成像。
如前所述,根据本发明的一方面的谱抑制脉冲的优化方法具有良好的通用性,可以广泛使用。并且,通过优化谱抑制脉冲,提高了谱信号抑制的能效。此外,利用设置有优化的谱抑制脉冲的MRI序列进行成像可以很好地抑制特定谱成分的信号,从而可以获得高质量的图像。下面将结合图5到图7对应用本发明上述方法的MATLAB仿真结果和活体内测试结果进行说明。
MATLAB仿真结果
图5(a)和5(b)示出了当谱抑制脉冲为FatSat脉冲时的仿真结果。其中各参数设置如下:采用损毁3D梯度回波序列进行成像,α=9°,每个脉冲串对应的K空间填充线数量=32,KSpaceCenterLine=#12,Trep=3.92ms,TrMin=12.298ms,TI1min=10.482ms,T1=290.0ms(3T)。根据本发明的上述通式(15)计算出的FatSat翻转角为:β=153.113°,其中在线#12处脂肪信号为0%。
图5(a)为具有上述给定参数和计算的翻转角的脂肪信号的仿真结果。图5(b)是图5(a)中的矩形框内的曲线部分的放大示意图。
图6(a)和6(b)示出了当谱抑制脉冲为SPAIR脉冲时的仿真结果。其中各参数设置如下:采用损毁3D梯度回波序列进行成像。α=9°,β=180°,每个脉冲串对应的K空间填充线数量=32,KSpaceCenterLine=#12,Trep=3.92ms,TrMin=14.918ms,TI1min=17.342ms,T1=290.0ms(3T)。根据本发明的上述通式(15)计算出的反转时间为:TrFill=0,TIFill=0.91ms,其中在线#12处脂肪信号为0%。
图6(a)为具有上述给定参数和计算的TrFill和TIFill的脂肪信号的仿真结果。图6(b)是图6(a)中的矩形框内的曲线部分的放大示意图。
活体测试结果
图7(a)和7(b)示出了根据本发明的、利用设置有优化的谱抑制脉冲的MRI序列形成的图像。具体的,图7(a)示出了谱抑制脉冲为FatSat脂肪抑制脉冲、每个FatSat脉冲串对应的K空间填充线数量为32线、计算出FatSat翻转角=153°时采集的图像;图7(b)示出了谱抑制脉冲为SPAIR脂肪抑制脉冲、每个SPAIR脉冲串对应的K空间填充线数量为32线、计算出TrFill=0,TIFill=62.3ms时采集的图像。图7(c)示出了使用未配置谱抑制脉冲的MRI序列采集的图像。在图7(a)-7(c)中显示的图像均是利用短T1损毁3D梯度回波序列获得的,并且在测试时均要求志愿者屏住呼吸以消除运动伪影。
图7(a)-7(c)下方的箭头均指示皮下区域。如图7(c)所示,当未进行脂肪抑制时,皮下区域显示高亮度的脂肪信号。而相应区域在图7(a)-7(b)则显示低暗信号。因此,根据本发明的谱抑制脉冲优化方法对于FatSat和SPAIR技术能够很好地工作。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (7)
1.一种磁共振成像方法,包括步骤:
根据谱抑制脉冲所抑制的谱成分的纵向磁化分量的稳定状态条件和所述纵向磁化分量的过零点条件,计算出该谱抑制脉冲的翻转角和/或反转时间;
应用包括设有所述翻转角和/或所述反转时间的所述谱抑制脉冲的磁共振成像序列进行磁共振成像,
其中在所述磁共振成像序列中所述谱抑制脉冲设置在一个或多个激发脉冲之前,所述纵向磁化分量的过零点条件是对应于K空间中心线的激发脉冲之前所述谱成分的纵向磁化分量为零,
其中,反转时间包括第一反转时间和第二反转时间,
所述计算出该谱抑制脉冲的翻转角或反转时间的步骤包括:利用通式计算谱抑制脉冲的翻转角或反转时间,所述通式为:
{E+F·C·(DN-1-1)/(D-1)+F·A·DN-1}/(1-F·B·DN-1)
={C·(DKSpaceCenterLine-1-1)/(D-1)+A·DKSpaceCenterLine-1}/(-B·DKSpaceCenterLine-1),
其中,
A=1-ETI1
B=cos(β)·ETI1
C=1-ETrep
D=cos(α)·ETrep
E=1-ETr
F=cos(α)·ETr
ETI1=exp(-TI1/T1)
ETrep=exp(-Trep/T1)
ETr=exp(-Tr/T1)
TI1=TI1Min+TIFill
Tr=TrMin+TrFill
α是激发脉冲的翻转角,
β是谱抑制脉冲的翻转角,
T1是谱抑制脉冲所抑制的谱成分的纵向弛豫时间,
N是两个谱抑制脉冲之间的激发脉冲的数目,
KSpaceCenterLine是K空间中心线的索引编号,
TI1Min是谱抑制脉冲的中心和该谱抑制脉冲之后的第一个激发脉冲的中心之间的最小时间间隔,
Trep是连续的激发脉冲之间的时间,
TrMin是磁共振成像序列的一个脉冲串中的最后一个激发脉冲的中心与下一个脉冲串的谱抑制脉冲的中心之间的最小时间间隔,
TIFill是所述第一反转时间,
TrFill是所述第二反转时间。
2.如权利要求1所述的磁共振成像方法,还包括:将所述第一反转时间TIFill和所述第二反转时间TrFill设置为0,
在所述计算出谱抑制脉冲的翻转角或反转时间的步骤中,计算所述翻转角β。
3.如权利要求1所述的磁共振成像方法,还包括步骤:将所述翻转角β设定为180度,
所述计算出所述谱抑制脉冲的翻转角或反转时间的步骤包括:
将所述第二反转时间TrFill设定为0,
计算所述第一反转时间TIFill。
4.如权利要求3所述的磁共振成像方法,其中当计算出的所述第一反转时间TIFill小于0时,所述计算出所述谱抑制脉冲的翻转角或反转时间的步骤进一步包括:
将所述第一反转时间TIFill设定为0,
计算所述第二反转时间TrFill。
5.如权利要求1所述的磁共振成像方法,还包括步骤:将所述翻转角β设定为180度,
所述计算出所述谱抑制脉冲的翻转角或反转时间的步骤包括:
将所述第一反转时间TIFill设定为0,
计算所述第二反转时间TrFill。
6.如权利要求5所述的磁共振成像方法,其中当计算出的所述第二反转时间TrFill小于0时,所述计算出所述谱抑制脉冲的翻转角或反转时间的步骤进一步包括:
将所述第二反转时间TrFill设定为0,
计算所述第一反转时间TIFill。
7.一种磁共振成像装置,包括:
一计算部件,根据谱抑制脉冲所抑制的谱成分的纵向磁化分量的稳定状态条件和所述纵向磁化分量的过零点条件,计算出该谱抑制脉冲的翻转角和/或反转时间;
一成像部件,应用包括设有所述翻转角和/或所述反转时间的所述谱抑制脉冲的磁共振成像序列进行磁共振成像,
其中在所述磁共振成像序列中所述谱抑制脉冲设置在一个或多个激发脉冲之前,所述纵向磁化分量的过零点条件是对应于K空间中心线的激发脉冲之前所述谱成分的纵向磁化分量为零,
其中,反转时间包括第一反转时间和第二反转时间,
所述计算出该谱抑制脉冲的翻转角或反转时间的步骤包括:利用通式计算谱抑制脉冲的翻转角或反转时间,所述通式为:
{E+F·C·(DN-1-1)/(D-1)+F·A·DN-1}/(1-F·B·DN-1)
={C·(DKSpaceCenterLine-1-1)/(D-1)+A·DKSpaceCenterLine-1}/(-B·DKSpaceCenterLine-1),
其中,
A=1-ETI1
B=cos(β)·ETI1
C=1-ETrep
D=cos(α)·ETrep
E=1-ETr
F=cos(α)·ETr
ETI1=exp(-TI1/T1)
ETrep=exp(-Trep/T1)
ETr=exp(-Tr/T1)
TI1=TI1Min+TIFill
Tr=TrMin+TrFill
α是激发脉冲的翻转角,
β是谱抑制脉冲的翻转角,
T1是谱抑制脉冲所抑制的谱成分的纵向弛豫时间,
N是两个谱抑制脉冲之间的激发脉冲的数目,
KSpaceCenterLine是K空间中心线的索引编号,
TI1Min是谱抑制脉冲的中心和该谱抑制脉冲之后的第一个激发脉冲的中心之间的最小时间间隔,
Trep是连续的激发脉冲之间的时间,
TrMin是磁共振成像序列的一个脉冲串中的最后一个激发脉冲的中心与下一个脉冲串的谱抑制脉冲的中心之间的最小时间间隔,
TIFill是所述第一反转时间,
TrFill是所述第二反转时间。
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210362793.2A CN103654779B (zh) | 2012-09-26 | 2012-09-26 | 磁共振成像方法和装置 |
KR1020130112894A KR20140040654A (ko) | 2012-09-26 | 2013-09-23 | 자기 공명 영상 촬영을 위한 방법 및 장치 |
US14/037,619 US9541624B2 (en) | 2012-09-26 | 2013-09-26 | Method and apparatus for magnetic resonance imaging |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210362793.2A CN103654779B (zh) | 2012-09-26 | 2012-09-26 | 磁共振成像方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103654779A CN103654779A (zh) | 2014-03-26 |
CN103654779B true CN103654779B (zh) | 2016-03-30 |
Family
ID=50294422
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210362793.2A Active CN103654779B (zh) | 2012-09-26 | 2012-09-26 | 磁共振成像方法和装置 |
Country Status (3)
Country | Link |
---|---|
US (1) | US9541624B2 (zh) |
KR (1) | KR20140040654A (zh) |
CN (1) | CN103654779B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109085523A (zh) * | 2018-07-23 | 2018-12-25 | 中国科学院深圳先进技术研究院 | 血管壁成像中脑脊液信号的抑制方法、装置、设备及介质 |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102014204995B4 (de) * | 2014-03-18 | 2015-10-01 | Siemens Aktiengesellschaft | Verfahren und Magnetresonanzanlage zur Fettsättigung |
US10335049B2 (en) | 2014-04-22 | 2019-07-02 | Case Western Reserve University | Magnetic resonance fingerprinting with steady state precession (MRF-FISP) |
CN105395200B (zh) * | 2015-12-02 | 2018-08-03 | 沈阳东软医疗系统有限公司 | 磁共振成像系统扫描方法、射频功率校准方法及装置 |
CN107219484B (zh) * | 2016-03-22 | 2023-09-22 | 西门子(深圳)磁共振有限公司 | 多层并行激发弥散成像的激发翻转角度确定方法和装置 |
CN109073720B (zh) * | 2016-04-26 | 2022-02-11 | 皇家飞利浦有限公司 | 静音磁共振指纹识别 |
CN106943143A (zh) * | 2017-03-28 | 2017-07-14 | 上海东软医疗科技有限公司 | 磁共振成像方法和装置 |
US10684343B2 (en) * | 2017-05-31 | 2020-06-16 | Canon Medical Systems Corporation | Magnetic resonance imaging apparatus and magnetic resonance imaging method |
CN108209918B (zh) * | 2017-12-30 | 2021-09-07 | 上海联影医疗科技股份有限公司 | 磁共振成像方法及磁共振系统 |
CN109620228B (zh) * | 2018-12-20 | 2023-01-03 | 上海联影医疗科技股份有限公司 | 磁共振成像中脂肪零点偏差校正方法及磁共振成像方法 |
CN109738841B (zh) * | 2019-02-21 | 2021-05-28 | 奥泰医疗系统有限责任公司 | 化学位移选择饱和压脂脉冲翻转角的优化方法 |
CN113805129B (zh) * | 2020-06-11 | 2024-01-02 | 西门子(深圳)磁共振有限公司 | 数据采集装置、方法及磁共振成像装置 |
CN111990997B (zh) * | 2020-08-25 | 2022-10-28 | 李友 | 一种基于反转恢复实图重建的水脂分离方法及系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7292039B1 (en) * | 2006-08-21 | 2007-11-06 | Siemens Medical Solutions Usa, Inc. | Contrast enhancement in magnetic resonance imaging |
US7706855B1 (en) * | 2004-07-20 | 2010-04-27 | General Electric Company | System and method for MR data acquisition with uniform fat suppression |
CN102772206A (zh) * | 2011-05-12 | 2012-11-14 | 上海联影医疗科技有限公司 | 磁共振成像方法 |
Family Cites Families (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6272369B1 (en) | 1999-01-22 | 2001-08-07 | Ge Medical Systems Global Technology Company Llc | Method for optimizing fat suppression using the chemical shift selective MR imaging technique |
US8010179B2 (en) * | 2003-09-25 | 2011-08-30 | General Electric Company | Method and apparatus of gradient echo imaging with on-the-fly optimization of tissue suppression |
WO2008044501A1 (fr) * | 2006-10-06 | 2008-04-17 | Hitachi Medical Corporation | Dispositif d'imagerie par résonance magnétique |
DE102006061177B4 (de) * | 2006-12-22 | 2009-04-02 | Siemens Ag | 3D-MR-Bildgebung mit Fettunterdrückung |
DE102007011807B3 (de) * | 2007-03-12 | 2008-10-23 | Siemens Ag | Sequenz für die Magnet-Resonanz-Bildgebung und Magnet-Resonanz-Gerät hierzu |
JP5398149B2 (ja) * | 2007-03-27 | 2014-01-29 | 株式会社東芝 | 磁気共鳴イメージング装置 |
WO2009093517A1 (ja) * | 2008-01-23 | 2009-07-30 | Hitachi Medical Corporation | 磁気共鳴イメージング装置及びマルチコントラスト画像取得方法 |
JP5236356B2 (ja) * | 2008-05-22 | 2013-07-17 | 株式会社日立メディコ | 核磁気共鳴イメージング装置 |
US9211081B2 (en) * | 2009-04-16 | 2015-12-15 | Ghost Medical Technologies, Ltd | Method for non-contrast enhanced magnetic resonance angiography |
WO2012060252A1 (ja) * | 2010-11-04 | 2012-05-10 | 株式会社 東芝 | 磁気共鳴イメージング装置及び磁気共鳴イメージング方法 |
US9081073B2 (en) * | 2011-01-27 | 2015-07-14 | Siemens Medical Solutions Usa, Inc. | System for suppression of artifacts in MR imaging |
US20130274592A1 (en) * | 2012-02-29 | 2013-10-17 | Taehoon SHIN | Time-resolved early-to-late gadolinium enhancement magnetic resonance imaging |
US9700230B2 (en) * | 2012-08-29 | 2017-07-11 | Toshiba Medical Systems Corporation | Enhanced fat saturation in myocardial infarction MRI |
DE102012220449B4 (de) * | 2012-11-09 | 2014-08-21 | Siemens Aktiengesellschaft | Verfahren und Magnetresonanzgerät zur automatisierten Klassifizierung einer Bildeigenschaft eines Magnetresonanzbildes |
DE102013215883B3 (de) * | 2013-08-12 | 2015-01-08 | Siemens Aktiengesellschaft | MR-Bildgebung mit Signalunterdrückung einer Spinspezies |
JP2015116474A (ja) * | 2013-11-13 | 2015-06-25 | 株式会社東芝 | 磁気共鳴イメージング装置 |
US9498139B2 (en) * | 2014-02-03 | 2016-11-22 | Toshiba Medical Systems Corporation | Background suppression by time dependent flip angle of saturation pulses |
US10335049B2 (en) * | 2014-04-22 | 2019-07-02 | Case Western Reserve University | Magnetic resonance fingerprinting with steady state precession (MRF-FISP) |
-
2012
- 2012-09-26 CN CN201210362793.2A patent/CN103654779B/zh active Active
-
2013
- 2013-09-23 KR KR1020130112894A patent/KR20140040654A/ko not_active Application Discontinuation
- 2013-09-26 US US14/037,619 patent/US9541624B2/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7706855B1 (en) * | 2004-07-20 | 2010-04-27 | General Electric Company | System and method for MR data acquisition with uniform fat suppression |
US7292039B1 (en) * | 2006-08-21 | 2007-11-06 | Siemens Medical Solutions Usa, Inc. | Contrast enhancement in magnetic resonance imaging |
CN102772206A (zh) * | 2011-05-12 | 2012-11-14 | 上海联影医疗科技有限公司 | 磁共振成像方法 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109085523A (zh) * | 2018-07-23 | 2018-12-25 | 中国科学院深圳先进技术研究院 | 血管壁成像中脑脊液信号的抑制方法、装置、设备及介质 |
CN109085523B (zh) * | 2018-07-23 | 2020-05-22 | 中国科学院深圳先进技术研究院 | 血管壁成像中脑脊液信号的抑制方法、装置、设备及介质 |
Also Published As
Publication number | Publication date |
---|---|
KR20140040654A (ko) | 2014-04-03 |
US20140084922A1 (en) | 2014-03-27 |
US9541624B2 (en) | 2017-01-10 |
CN103654779A (zh) | 2014-03-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103654779B (zh) | 磁共振成像方法和装置 | |
CN105005011B (zh) | 在三维快速自旋回波中的混合采集的磁共振成像方法 | |
JP5591493B2 (ja) | 磁気共鳴イメージング装置 | |
US7786729B2 (en) | Method for magnetic resonance imaging | |
US20080278159A1 (en) | Turbospin echo imaging sequence with long echo trains and optimized t1 contrast | |
CN107949325A (zh) | 磁共振成像装置、扩散加权图像的生成方法以及图像处理装置 | |
EP2610632A1 (en) | MRI with Dixon-type water/fat separation and prior knowledge about inhomogeneity of the main magnetic field | |
JPS6047945A (ja) | Nmr方法 | |
JP6285669B2 (ja) | 磁気共鳴イメージング装置 | |
CN103505214A (zh) | 用于mr图像的同时双平板获取的系统 | |
CN110095742A (zh) | 一种基于神经网络的平面回波成像方法和装置 | |
US6586933B1 (en) | Method and system for MRI with lipid suppression | |
JPH05285116A (ja) | 磁気共鳴イメージング方法 | |
US9535148B2 (en) | Dynamic contrast enhanced magnetic resonance imaging with high spatial-temporal resolution | |
US20150301145A1 (en) | Multislice acquisition with incoherent aliasing | |
JP7023954B2 (ja) | プロペラmrイメージング | |
US7167003B2 (en) | Measuring method of spatially resolved frequency selective images | |
WO2018213358A1 (en) | Methods, systems, and computer readable media for partition-encoded simultaneous multislice (prism) imaging | |
CN101571579B (zh) | 基于反转恢复的组织分离成像方法 | |
US10114095B2 (en) | Multi-spectral MRI scan with magnetization recovery | |
JP4076655B2 (ja) | 磁気共鳴イメージング装置 | |
US7821263B2 (en) | Self-refocused spatial-spectral pulse | |
US11609293B2 (en) | Fat saturation technique with linear flexible reordering and bloch equation calculation techniques | |
US11255937B2 (en) | K-space data acquisition device and method, and magnetic resonance imaging device and method | |
US20230366965A1 (en) | Highly Accelerated Sub-Millimeter Resolution 3D GRASE with Controlled T2 Blurring in T2-Weighted Functional MRI at 7 Tesla |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | 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 |