CN111352054B - 一种振荡梯度准备的3d梯度自旋回波成像方法及设备 - Google Patents

一种振荡梯度准备的3d梯度自旋回波成像方法及设备 Download PDF

Info

Publication number
CN111352054B
CN111352054B CN202010243339.XA CN202010243339A CN111352054B CN 111352054 B CN111352054 B CN 111352054B CN 202010243339 A CN202010243339 A CN 202010243339A CN 111352054 B CN111352054 B CN 111352054B
Authority
CN
China
Prior art keywords
gradient
module
spin echo
diffusion
epi
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
Application number
CN202010243339.XA
Other languages
English (en)
Other versions
CN111352054A (zh
Inventor
吴丹
徐义程
李浩天
张祎
孙毅
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN202010243339.XA priority Critical patent/CN111352054B/zh
Publication of CN111352054A publication Critical patent/CN111352054A/zh
Application granted granted Critical
Publication of CN111352054B publication Critical patent/CN111352054B/zh
Priority to PCT/CN2021/074841 priority patent/WO2021196865A1/zh
Priority to JP2022558224A priority patent/JP7304106B2/ja
Priority to US17/952,341 priority patent/US20230125982A1/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/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/563—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution of moving material, e.g. flow contrast angiography
    • G01R33/56341—Diffusion 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/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
    • G01R33/4826—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 in three dimensions
    • 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
    • G01R33/4833—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy using spatially selective excitation of the volume of interest, e.g. selecting non-orthogonal or inclined slices
    • G01R33/4835—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy using spatially selective excitation of the volume of interest, e.g. selecting non-orthogonal or inclined slices of multiple slices
    • 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/5611—Parallel magnetic resonance imaging, e.g. sensitivity encoding [SENSE], simultaneous acquisition of spatial harmonics [SMASH], unaliasing by Fourier encoding of the overlaps using the temporal dimension [UNFOLD], k-t-broad-use linear acquisition speed-up technique [k-t-BLAST], k-t-SENSE
    • 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/5616—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 gradient refocusing, e.g. EPI
    • 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/5618—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 both RF and gradient refocusing, e.g. 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/5619—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 by temporal sharing of data, e.g. keyhole, block regional interpolation scheme for k-Space [BRISK]

Landscapes

  • Physics & Mathematics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • High Energy & Nuclear Physics (AREA)
  • General Health & Medical Sciences (AREA)
  • Signal Processing (AREA)
  • Engineering & Computer Science (AREA)
  • Radiology & Medical Imaging (AREA)
  • Vascular Medicine (AREA)
  • Optics & Photonics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

本发明公开一种振荡梯度准备的3D梯度自旋回波成像方法及设备。该成像方法包括以下步骤：首先，利用全局饱和模块中破坏先前残余的横向磁化；其次，通过扩散编码模块，将一对梯形余弦振荡梯度嵌入到90°x‑180°y‑90°‑x射频脉冲中，将扩散编码与信号采集分开；然后，使用脂肪饱和模块对脂肪信号进行抑制；最后，采用梯度自旋回波的读出方式采集信号，并使用多路复用灵敏度编码重建校正了多次激发之间的相位误差。与在3T临床系统上使用的基于2D平面回波的振荡梯度弥散序列相比，3D振荡梯度准备的梯度自旋回波序列有效提高了成像时间和信噪比，有助于时间依赖性的弥散MRI技术的临床转化。

Description

一种振荡梯度准备的3D梯度自旋回波成像方法及设备
技术领域
本申请涉及磁共振技术领域，尤其涉及弥散磁共振成像领域。
背景技术
扩散MRI(dMRI)可以根据生物组织中的水分子在微结构环境下受限扩散的情况对组织的微结构进行探测。通过改变扩散编码方案，可以采集到多维的 dMRI信号，该信号强度不仅取决于扩散的方向和b值，还取决于扩散时间和弛豫时间。具体来说，扩散时间决定了dMRI对在不同空间尺度下的微结构特征的敏感性。在不同的扩散时间测量dMRI信号，即可获取水分子扩散的时间依赖性，而时间依赖性曲线反映了重要的微结构特性，例如细胞大小，细胞内分数，表面积与体积比等。相比于采用传统脉冲梯度dMRI(PG-dMRI)序列，基于振荡梯度dMRI(OG-dMRI)的采集技术可能实现更短的扩散时间(<5ms)，更有利于扩散参数的短时测量。
振荡编码梯度由多个周期的余弦或梯形的余弦波形组成，目前已经被广泛用于在动物成像系统中对动物和水模的时间依赖性dMRI的研究。然而，该技术转化到临床应用面临的主要挑战是需要很强的振荡梯度。由于OG波形的b值与其振荡频率的立方成反比，所以即使使用最大梯度为80mT/m的高端临床扫描仪，最高振荡频率也被限制在40-60Hz之间，对应的b值约为300-500s/mm2。这导致了临床扫描仪上可以实现的短扩散时间有限且信号的对比度-噪声比较低。同时，常规的OG-dMRI序列仍然存在其他的问题会降低图像质量并阻碍OG-dMRI 的临床转换。一方面，通常需要具有多个周期的长振荡梯度才能在临床扫描仪上达到合理的b值，因此图像采集时需要采用较长的回波时间(TE)，进而导致信噪比(SNR)低的问题。另一方面，对于2D多层面(multislice)采集，由于需要对每个层面重复应用强振荡梯度，因此会导致较高的占空比和系统产热量。综上所述，与PG-dMRI相比，OG-dMRI需要使用较长的重复时间(TRs)以降低占空比，使得总采集时间变长。
发明内容
为了克服现有序列中存在的不足，本发明提出了一种振荡梯度准备的3D梯度自旋回波序列(OGprep-GRASE)，以解决上述提及的两个OG-dMRI序列的问题，提高dMRI成像性能。
为了达到上述目的，本发明采用以下技术方案予以实现：
第一方面，本发明提供了一种振荡梯度准备的3D梯度自旋回波成像方法，其包括：
S1：在序列开始的位置添加全局饱和模块，以破坏之前残余的横向磁化；然后等待一段饱和后延迟(PSD)时间，使纵向磁化矢量恢复；
S2：在饱和后延迟(PSD)之后添加扩散准备模块，将一对梯形余弦振荡梯度或脉冲梯度嵌入到90°x-180°y-90°-x射频脉冲中，以实现扩散编码与信号采集的分离；
S3：在扩散准备模块后添加脂肪饱和模块，抑制脂肪信号；
S4：在脂肪饱和模块后，使用梯度自旋回波序列模块对3D K空间进行信号采集；
S5：使用多路复用灵敏度编码重建(MUSE)模块校正多次激发之间存在的相位误差。
基于上述第一方面的方案，各步骤还可以进一步提供以下优选的实现方式。需要注意的是，各优选方式中的技术特征在没有冲突的情况下均可进行相互组合。当然这些优选方式也可以通过其他能够实现相同技术效果的方式实现，不构成限制。
作为优选，所述步骤S1中的全局饱和模块建立方法如下：
每施加一次射频脉冲，同时在X、Y、Z三个方向施加梯度，将上次信号采集时残留的横向磁化完全消除；施加3次射频脉冲后，等待一段饱和后延迟(PSD) 时间，所述饱和后延迟(PSD)为饱和度模块和扩散编码模块之间的时间间隔。
作为优选，所述步骤S2中的扩散准备模块建立方法如下：
首先，沿X轴方向进行90°翻转角的非层选性硬脉冲激励；然后，沿着设置弥散的方向施加梯形余弦振荡梯度；接着，沿Y轴方向进行180°翻转角的非层选性绝热正切脉冲，然后施加完全相同的梯形余弦振荡梯度；之后沿Z轴方向施加稳定器梯度(Gstb)；最后，沿-X轴方向进行90°翻转角的非层选性方硬脉冲激励，将横向磁化矢量转化为纵向磁化矢量。
作为优选，所述步骤S3中的脂肪饱和模块建立方法如下：
施加110°翻转角的频率选择性的高斯脉冲，将脂肪纵向磁化矢量的cos(110°) 部分进行翻转并等待T1*ln(2)时间后将其置0，并同时施加Z轴方向的梯度，将脂肪的横向磁化矢量完全消除，其中T1为脂肪的纵向弛豫时间。
作为优选，所述步骤S4中的梯度自旋回波序列模块建立方法如下：
通过在Y方向上进行回波平面成像(EPI)编码以及在Z方向上进行涡轮自旋回波(TSE)编码来实现3D GRASE读出；设定自旋回波个数为NSE，梯度回波个数为NEPI；沿EPI方向可执行分段读取，而沿TSE方向上强制执行中心编码；在每段EPI读出的前后沿Z轴方向施加稳定器梯度。
作为优选，所述步骤S5中的多路复用灵敏度编码重建模块建立方法如下：
首先使用常规的灵敏度编码(SENSE)技术来估计运动引起的多个EPI段之间的相位变化，然后同时在所有交错EPI中联合计算混叠体素的幅度信号。
第二方面，本发明提供了一种磁共振成像设备，其包括磁共振扫描器以及控制单元，所述控制单元中存储有成像序列，所述成像序列被执行时，用于实现如上述第一方面中任一方案所述的振荡梯度准备的3D梯度自旋回波成像方法。
相对于现有技术，本发明具有以下特点：提出了一种3D振荡梯度制备梯度自旋回波序列，首次将振荡梯度扩散编码序列与梯度自旋回波序列结合起来。与传统的临床扫描仪上OG-dMRI序列进行的2D-EPI采集相比,3D OGprep-GRASE 序列不仅缩短了图像的扫描时间而且提高了信噪比。通过比较PG和OG-dMRI 对成人大脑进行成像的结果，表明无论使用GRASE或者EPI读出序列，人脑的灰质和白质都表现出明显的扩散时间依赖性。
附图说明
图1是3D OGprep-GRASE序列的时序图。
图2是对分别使用两次激发的OGprep-GRASE和PGprep-GRASE序列进行 MUSE重建的前后结果对比。
图3是对OG-dMRI分别使用2D EPI和3D GRASE读出序列的扫描时间和信噪比的对比结果。
图4是OG-dMRI分别使用2D EPI和3D GRASE读出序列对DTI量化指标：表观扩散系数图、各向异性分数图、以及扩散方向编码彩图的对比结果。
图5是OG-dMRI分别使用2D EPI和3D GRASE读出序列得到的DTI量化指标的扩散时间依赖性对比结果。
具体实施方式
如附图1下面基于本发明提出的方法结合实施例展示其具体的技术效果，以便本领域技术人员更好地理解本发明的实质。
参见图1所示，展示了本发明的一种较优实现方式中提供的振荡梯度准备的 3D梯度自旋回波序列(3D OGprep-GRASE)，该序列是一种新型的振荡梯度场弥散磁共振三维成像序列。下面详细描述该序列具体方法包括以下步骤：
步骤一：在序列开始的位置添加全局饱和模块，以破坏之前残余的横向磁化。全局饱和模块的具体做法为：在序列时序开始时，首先施加B1射频脉冲；每施加一次射频脉冲，同时在X、Y、Z三个方向施加梯度，三个梯度的幅值完全相等。在施加3次射频脉冲后，可以认为已经将上次信号采集时残留的横向磁化完全消除。之后等待一段被称为饱和后延迟(PSD)的时间，等待纵向磁化矢量的恢复。本发明中，将饱和度模块和扩散编码模块之间的时间间隔记为PSD。
步骤二：在PSD之后紧接着扩散准备模块，将一对梯形余弦振荡梯度或脉冲梯度嵌入到90°x-180°y-90°-x射频脉冲中，以实现扩散编码与信号采集的分离。扩散准备模块的具体做法为：首先，沿X轴方向发射90°翻转角的非层选性的硬脉冲激励；然后，沿着设置弥散的方向施加梯形余弦振荡梯度，梯度的幅值大小与设置的b值有关；接着，由于绝热正切脉冲对B1场的不敏感性，沿Y轴方向进行180°翻转角的绝热正切脉冲激励，将散向的横向磁化矢量进行重聚；然后，施加与之前施加的幅值、形状、持续时间等完全相同的梯形余弦振荡梯度；之后，沿Z轴方向施加稳定器梯度(Gstb)，以减轻与相位误差有关的信号调制，并将梯度幅值设置为与GRASE读出模块中的破碎机梯度(crusher)相同的幅度。最后，沿-X轴方向发射90°翻转角的非片选性的硬脉冲激励，将重聚的横向磁化矢量转化为纵向磁化矢量。扩散准备模块中的回波时间表示为TE1。
步骤三：在扩散准备模块后添加脂肪饱和模块，抑制脂肪信号。脂肪饱和模块的具体做法为：施加110°翻转角的频率选择性的高斯脉冲，将脂肪纵向磁化矢量的cos(110°)部分进行翻转，等待T1*ln(2)时间，当翻转后恢复的磁化矢量为0时施加步骤四的成像模型，并同时施加Z轴方向的梯度，将脂肪的横向磁化矢量完全消除。
步骤四：在脂肪饱和模块后，使用梯度自旋回波序列模块对3D K空间进行信号采集。梯度自旋回波序列模块的具体做法为：通过在Y方向上进行回波平面成像(EPI)编码以及在Z方向上进行涡轮自旋回波(TSE)编码来实现3D GRASE读出。在EPI方向上的相位编码的数目被表示为NEPI，而在TSE方向上的相位编码的数目被表示为NSE，即自旋回波个数为NSE，梯度回波个数为NEPI。沿EPI方向可执行分段读取，此处分段读取是沿EPI方向对NEPI进行整段或分段采集，而沿TSE方向上强制执行中心编码，即第一次聚焦脉冲后的EPI读取信号的方式在kz方向填充k-space中心。稳定器梯度沿Z轴方向在每段EPI读出的前后与重聚焦脉冲旁边的破碎机梯度同时施加。GRASE模块中的回波时间表示为TE2，它等于之后以TSE方式进行信号采集的回波间隔(ESP)。因此，所获取信号的总回波时间为TE1+TE2。
当把EPI方向的相位编码分成N段进行信号采集时，对每一段的信号采集都需要重复一遍步骤一到步骤四，这样相对于常规采集方式，分段采集方式的 NEPI减少，每次读出时间较短，信号衰减较小，可提高图像分辨率和减小EPI的变形效应，但采集时间变为原来的N倍，同时对运动较为敏感，需要矫正每段之间的运动误差。但是极大地提高了图像的信噪比。
步骤五：使用多路复用灵敏度编码(MUSE)重建模块校正多次激发之间存在的相位误差。也就是说，步骤四中沿EPI方向如果使用分段信号进行采集，则需要MUSE重建来矫正每段采集到的信号之间的运动误差。MUSE重建时，首先使用常规的灵敏度编码(SENSE)技术来估计运动引起的多个EPI段之间的相位变化，然后同时在所有交错EPI中联合计算混叠体素的幅度信号。与常规的 SENSE方法相比，MUSE方法大大改善了矩阵求逆条件，因此可以得到更高信噪比的图像。
由此可见，通过上述步骤一到五，本发明中的3D OGprep-GRASE序列，不仅缩短了图像的扫描时间而且提高了信噪比。为了进一步展示该做法的效果，本发明还进行了下列步骤。
步骤六：经过步骤一到五处理后，再进行后处理得到DTI的量化指标：表观扩散系数(ADC)、各向异性分数(FA)、以及扩散方向编码彩图(DEC)。其中后处理方法如下：将获取到的沿X、Y、Z三个方向施加弥散梯度的不同b 值的图像取出，分别计算图像每个像素点在X、Y、Z三个方向的扩散值Dxx、 Dyy、Dzz，三者取均值即可得到ADC值；其中任意一个方向的扩散值D的计算通式为：
D＝-ln(S2/S1)/(b2-b1)
式中S1和S2分别为当前计算方向中两个不同b值下的图像像素值。
另外，根据扩散的计算公式求得6个方向的D值(D的下标表示方向)，6 个方向的D值可以构建成矩阵D：
Figure BDA0002433282700000061
通过求取矩阵D的3个特征向量和3个特征值λ1、λ2、λ3，再通过以下公式求得图像中每个体素的FA值：
Figure BDA0002433282700000062
将3个特征值中最大特征值λ对应的特征向量当作扩散的主方向，根据三维空间方向的不同绘制了扩散方向编码彩图。
下面基于上述方法的步骤一到六，结合实施例对其技术效果进行展示，以便本领域技术人员更好地理解本发明的实质。
实施例
将上述的3D振荡梯度准备的梯度自旋回波序列方法在7名健康的年轻男性志愿者中进行了测试，下面介绍此处的具体参数：核磁共振扫描使用了西门子 Prisma 3.0T扫描仪，所有扫描均使用64通道头部线圈在最大梯度为80mT/m上进行。本实施例中进行了两组实验。
实验一：为了比较3D OGprep-GRASE和2D OG-EPI序列的扫描时间和图像信噪比，以50Hz的振荡梯度，b＝500s/mm2、12个方向，2次重复，FOV＝ 220×220mm，分辨率为2.75×2.75mm2，切片厚度为3mm(3D GRASE为 2.75×2.75×3mm3)执行OG-dMRI扫描，采用以下方案：
(1)单次激发的GRASE读出序列，NEPI＝80，NSE＝10，turbo因子＝12，EPI因子＝63，层面过采样＝20％，部分傅立叶因子＝61％，带宽＝2,500Hz， TE1/TE2/TR＝124/33.6/3000ms，扫描时间＝1.5min；与之对应的2D-EPI读出序列，扫描了10个层面，部分傅立叶因子＝6/8，带宽＝1,690Hz，TE/TR＝158/4200ms (最小TR)，扫描时间＝2.0min。
(2)两次激发的GRASE读出序列，NEPI＝80，NSE＝20，turbo因子＝24，EPI 因子＝27，TE1/TE2/TR＝124/22.9/3000ms，扫描时间＝2.9min，其他参数与单次激发的GRASE读出序列的参数一致；与之对应的2D-EPI读出序列，扫描了20个层面，部分傅里叶因子＝6/8，TE/TR＝147/8400ms(最小TR)，扫描时间＝4.0min，其他参数与单次激发对应的2D-EPI相同。
(3)三次激发的GRASE读出序列，NEPI＝80，NSE＝30，turbo因子＝36，EPI 因子＝21，TE1/TE2/TR＝124/20.2/3000ms，扫描时间为4.3min，其他参数与单次激发的GRASE读出序列的参数一致；与之对应的2D-EPI读出序列，扫描了20 个层面，部分傅里叶因子＝6/8，TE/TR＝144/12600ms(最小TR)，扫描时间＝6.0min，其他参数与单次激发对应的2D-EPI相同。
实验二：为了测试本发明中的序列对不同扩散时间下扩散系数测量的可靠性，以25Hz和50Hz的振荡频率，b＝600s/mm2，6个方向，TE1/TE2/TR＝84/32/3000ms 执行单次激发的OGprep-GRASE序列，矩阵尺寸＝80×80×10，成像分辨率＝2.75×2.75×5mm；PGprep-GRASE序列，其δ＝20ms，Δ＝30和60ms，其他参数与OGprep序列一致。作为比较，还以25Hz，50Hz的振荡梯度和δ＝20ms，Δ＝30 /60ms的脉冲梯度以及其他与prep-GRASE扫描一致的参数采集了2D EPI数据。最后采集了一组受试者的两次激发PG-dMRI和OG-dMRI，以测试多路复用灵敏度编码(MUSE)重建算法对运动伪影校正的有效性。
附图2显示了MUSE重建前后使用OG-和PG-prepared GRASE序列采集得到的图像。从图中可以看出，对于两次激发的GRASE采集方式，运动伪影在 OG-dMRI图像中不明显，但在PG-dMRI图像中明显，并且MUSE可以在很大程度上纠正运动引起的相位误差。
附图3显示了在实验1中对于10层面和20层面的采集协议，与2D EPI序列相比，3DOGprep-GRASE的扫描速度分别提高了1.34和1.38倍。对于10层面协议，GRASE和EPI序列显示出相似的图像信噪比，但是对于3D GRASE序列，当层面数量增加到20层面时，图像信噪比增加了一倍，而2D EPI数据的图像信噪比保持不变。
附图4是OG-dMRI分别使用2D EPI和3D GRASE的读出数据经过弥散张量重建后得到的指标表观弥散系统(ADC)，各向异性(FA)和方向编码彩图 (DEC)的结果。从图像中可以看到在10、20和30层面协议的2D EPI数据中可以观察到噪声点。在10层面的3D GRASE数据中，噪声类似于EPI数据，但是随着层面数量增加到20和30，图像逐渐变得清晰。说明3DGRASE读出序列采集得到了较高信噪比的图像。
在OG-dMRI(50Hz和25Hz)和PG-dMRI(30ms和60ms)的3D GRASE 和2D EPI序列成像得到的图像中，分别求得在人脑的皮层下白质和深部灰质中的ADC数据。如附图5所示，来自两个序列的ADC值在白质和灰质区域均表现出显著的时间依赖性(p<0.0001)，并且通过双向ANOVA测试并未发现显著的序列差异，说明3D OGprep-GRASE序列可有效测量ADC的时间依赖信，结果与普通2D EPI序列一致。
由此可见，与在3T临床系统上使用的基于2D平面回波的振荡梯度弥散序列相比，本发明的3D OGprep-GRASE序列有效提高了成像时间和信噪比，有助于时间依赖性的弥散MRI技术的临床转化。
基于上述3D OGprep-GRASE序列，在其他实施例中还可以进一步提供一种磁共振成像设备，其包括磁共振扫描器以及控制单元，控制单元中存储有成像序列，即3D OGprep-GRASE序列。该成像序列被执行时，用于实现前述步骤一到五中的振荡梯度准备的3D梯度自旋回波成像方法。
需要注意的是，此处的磁共振扫描器和控制单元的硬件结构可采用常规磁共振成像系统实现。磁共振扫描器应当包括磁体部分和磁共振波谱仪部分，控制单元中应当含有实现成像必要的数据处理和图像重建软件，3D OGprep-GRASE序列可由控制单元读取并用于控制磁共振扫描器。也就是说，3D OGprep-GRASE 序列可以直接应用至相应的磁共振成像设备中。
需要指出的是，以上所述的实施例只是本发明的一种较佳的方案，然其并非用以限制本发明。有关技术领域的普通技术人员，在不脱离本发明的精神和范围的情况下，还可以做出各种变化和变型。因此凡采取等同替换或等效变换的方式所获得的技术方案，均落在本发明的保护范围内。

Claims (7)

1.一种振荡梯度准备的3D梯度自旋回波成像方法，其特征在于，包括：
S1：在序列开始的位置添加全局饱和模块，以破坏之前残余的横向磁化；然后等待一段饱和后延迟(PSD)时间，使纵向磁化矢量恢复；
S2：在饱和后延迟(PSD)之后添加扩散准备模块，将一对梯形余弦振荡梯度或脉冲梯度嵌入到90°x-180°y-90°-x射频脉冲中，以实现扩散编码与信号采集的分离；
S3：在扩散准备模块后添加脂肪饱和模块，抑制脂肪信号；
S4：在脂肪饱和模块后，使用梯度自旋回波序列模块对3D K空间进行信号采集；
S5：使用多路复用灵敏度编码重建模块校正多次激发之间存在的相位误差。
2.如权利要求1所述的振荡梯度准备的3D梯度自旋回波成像方法，其特征在于，所述步骤S1中的全局饱和模块建立方法如下：
每施加一次射频脉冲，同时在X、Y、Z三个方向施加梯度，将上次信号采集时残留的横向磁化完全消除；施加3次射频脉冲后，等待一段饱和后延迟(PSD)时间，所述饱和后延迟(PSD)为饱和度模块和扩散编码模块之间的时间间隔。
3.如权利要求1所述的振荡梯度准备的3D梯度自旋回波成像方法，其特征在于，所述步骤S2中的扩散准备模块建立方法如下：
首先，沿X轴方向进行90°翻转角的非层选性硬脉冲激励；然后，沿着设置弥散的方向施加梯形余弦振荡梯度；接着，沿Y轴方向进行180°翻转角的非层选性绝热正切脉冲，然后施加完全相同的梯形余弦振荡梯度；之后沿Z轴方向施加稳定器梯度(Gstb)；最后，沿-X轴方向进行90°翻转角的非层选性方硬脉冲激励，将横向磁化矢量转化为纵向磁化矢量。
4.如权利要求1所述的振荡梯度准备的3D梯度自旋回波成像方法，其特征在于，所述步骤S3中的脂肪饱和模块建立方法如下：
施加110°翻转角的频率选择性的高斯脉冲，将脂肪纵向磁化矢量的cos(110°)部分进行翻转并等待T1*ln(2)时间后将其置0，并同时施加Z轴方向的梯度，将脂肪的横向磁化矢量完全消除，其中T1为脂肪的纵向弛豫时间。
5.如权利要求1所述的振荡梯度准备的3D梯度自旋回波成像方法，其特征在于，所述步骤S4中的梯度自旋回波序列模块建立方法如下：
通过在Y方向上进行回波平面成像(EPI)编码以及在Z方向上进行涡轮自旋回波(TSE)编码来实现3D GRASE读出；设定自旋回波个数为NSE，梯度回波个数为NEPI；沿EPI方向可执行分段读取，而沿TSE方向上强制执行中心编码；在每段EPI读出的前后沿Z轴方向施加稳定器梯度。
6.如权利要求1所述的振荡梯度准备的3D梯度自旋回波成像方法，其特征在于，所述步骤S5中的多路复用灵敏度编码重建模块建立方法如下：
首先使用常规的灵敏度编码(SENSE)技术来估计运动引起的多个EPI段之间的相位变化，然后同时在所有交错EPI中联合计算混叠体素的幅度信号。
7.一种磁共振成像设备，其特征在于，包括磁共振扫描器以及控制单元，所述控制单元中存储有成像序列，所述成像序列被执行时，用于实现如权利要求1～6任意所述的振荡梯度准备的3D梯度自旋回波成像方法。
CN202010243339.XA 2020-03-31 2020-03-31 一种振荡梯度准备的3d梯度自旋回波成像方法及设备 Active CN111352054B (zh)

Priority Applications (4)

Application Number Priority Date Filing Date Title
CN202010243339.XA CN111352054B (zh) 2020-03-31 2020-03-31 一种振荡梯度准备的3d梯度自旋回波成像方法及设备
PCT/CN2021/074841 WO2021196865A1 (zh) 2020-03-31 2021-02-02 一种振荡梯度准备的3d梯度自旋回波成像方法及设备
JP2022558224A JP7304106B2 (ja) 2020-03-31 2021-02-02 振動勾配準備の３ｄ勾配スピンエコーイメージング方法及び機器
US17/952,341 US20230125982A1 (en) 2020-03-31 2022-09-26 3d oscillating gradient-prepared gradient spin-echo imaging method, and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010243339.XA CN111352054B (zh) 2020-03-31 2020-03-31 一种振荡梯度准备的3d梯度自旋回波成像方法及设备

Publications (2)

Publication Number Publication Date
CN111352054A CN111352054A (zh) 2020-06-30
CN111352054B true CN111352054B (zh) 2020-10-13

Family

ID=71197640

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010243339.XA Active CN111352054B (zh) 2020-03-31 2020-03-31 一种振荡梯度准备的3d梯度自旋回波成像方法及设备

Country Status (4)

Country Link
US (1) US20230125982A1 (zh)
JP (1) JP7304106B2 (zh)
CN (1) CN111352054B (zh)
WO (1) WO2021196865A1 (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111352054B (zh) * 2020-03-31 2020-10-13 浙江大学 一种振荡梯度准备的3d梯度自旋回波成像方法及设备
CN112763958B (zh) * 2020-12-10 2022-06-21 复旦大学 一种基于神经网络的多次激发平面回波磁共振成像方法
CN113476031B (zh) * 2021-08-04 2024-03-29 浙江大学 反转恢复准备的3d梯度自旋回波扩散成像方法、介质及设备
CN113970716B (zh) * 2021-10-18 2024-04-23 上海联影医疗科技股份有限公司 梯度回波序列压脂方法及磁共振成像方法及设备

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6320380B1 (en) * 2000-10-03 2001-11-20 Marconi Medical Systems, Inc. MRI method and apparatus for increasing the efficiency of echo lanar imaging and other late echo techniques
CN107548464A (zh) * 2014-12-15 2018-01-05 皇家飞利浦有限公司 自旋回波mr成像
CN108720834A (zh) * 2018-02-06 2018-11-02 苏州朗润医疗系统有限公司 一种梯度回波多回波水脂分离方法及应用该方法的磁共振成像系统
CN109115820A (zh) * 2018-07-02 2019-01-01 清华大学 基于平面回波成像的磁共振水脂分离和定量方法及装置

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5245282A (en) * 1991-06-28 1993-09-14 University Of Virginia Alumni Patents Foundation Three-dimensional magnetic resonance imaging
US7898254B2 (en) * 2008-02-19 2011-03-01 Advanced Mri Technologies, Llc Arterial spin labeled, segmented, interleaved 3D GRASE MRI
US9030201B2 (en) * 2011-01-27 2015-05-12 Siemens Medical Solutions Usa, Inc. System and method for independent manipulation of a fat and a water component in magnetic resonance imaging
US10531812B2 (en) * 2013-01-16 2020-01-14 Beth Israel Deaconess Medical Center, Inc. System and method for improved cardiac imaging of subjects with adverse cardiac conditions
US10092199B2 (en) * 2014-04-01 2018-10-09 Siemens Healthcare Gmbh MR imaging apparatus and method for generating a perfusion image with motion correction
CN106344015A (zh) * 2015-07-15 2017-01-25 四川大学华西医院 一种异常扩散程度加权的弥散磁共振成像方法
US10330762B2 (en) * 2016-04-12 2019-06-25 The Johns Hopkins University Measurement of blood volume using velocity-selective pulse trains on MRI
EP3235424A1 (en) * 2016-04-19 2017-10-25 Université d'Aix Marseille Method for enhancing the ihmt sensitivity of steady-state gradient echo acquisitions in an mri system
US20210228096A1 (en) * 2020-01-28 2021-07-29 The Johns Hopkins University Measurement of blood volume using fourier-transform based velocity-selective pulse trains on mri
CN111352054B (zh) * 2020-03-31 2020-10-13 浙江大学 一种振荡梯度准备的3d梯度自旋回波成像方法及设备

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6320380B1 (en) * 2000-10-03 2001-11-20 Marconi Medical Systems, Inc. MRI method and apparatus for increasing the efficiency of echo lanar imaging and other late echo techniques
CN107548464A (zh) * 2014-12-15 2018-01-05 皇家飞利浦有限公司 自旋回波mr成像
CN108720834A (zh) * 2018-02-06 2018-11-02 苏州朗润医疗系统有限公司 一种梯度回波多回波水脂分离方法及应用该方法的磁共振成像系统
CN109115820A (zh) * 2018-07-02 2019-01-01 清华大学 基于平面回波成像的磁共振水脂分离和定量方法及装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
A Simple Approach for Three-Dimensional ,apping of Baseline Cerebrospinal Fluid Volume Fraction;Qin Qin;《MAGNETIC RESONANCE IN MEDICINE》;20101130;第385-391页 *
基于扩散磁共振成像的大脑白质微结构检测研究进展;张力新等;《纳米技术与精密工程》;20150731;第13卷(第4期);第276-286页 *

Also Published As

Publication number Publication date
WO2021196865A1 (zh) 2021-10-07
US20230125982A1 (en) 2023-04-27
CN111352054A (zh) 2020-06-30
JP2023510031A (ja) 2023-03-10
JP7304106B2 (ja) 2023-07-06

Similar Documents

Publication Publication Date Title
CN111352054B (zh) 一种振荡梯度准备的3d梯度自旋回波成像方法及设备
EP2414861B1 (en) Magnetic resonance imaging with improved imaging contrast
US8406849B2 (en) Systems and methods for magnetic resonance imaging
US9778336B2 (en) System and method for rapid, multi-shot segmented magnetic resonance imaging
KR101625733B1 (ko) 다수의 자기 공명 이미지를 결정하기 위한 방법 및 자기 공명 시스템
US9645211B2 (en) Magnetic resonance imaging using steering-propeller
EP1846777A1 (en) Mr multi-slice steady state free precession imaging
EP3673281B1 (en) Dixon-type water/fat separation mr imaging
CN109219757A (zh) Dixon型水/脂肪分离MR成像
CN109983358B (zh) Propeller mr成像
CN108431625B (zh) 具有对运动引起的扩散梯度不一致性的修正的dti
WO2018114554A1 (en) Dixon-type water/fat separation mr imaging
WO2012140543A1 (en) Mri of chemical species having different resonance frequencies using an ultra-short echo time sequence
CN113805133B (zh) 一种减少图像失真的磁共振平面回波成像方法
CN111164444B (zh) 具有经改进的脂肪位移校正的Dixon型水/脂肪分离MR成像
EP3185029A1 (en) Mr imaging using propeller acquisition with t2 decay correction
US20050116710A1 (en) Magnetic resonance imaging
Chang et al. Navigator-free submillimeter diffusion imaging using multishot-encoded simultaneous multi-slice (MUSIUM)
CN113476031B (zh) 反转恢复准备的3d梯度自旋回波扩散成像方法、介质及设备
RU2773486C2 (ru) Магнитно-резонансная визуализация с разделением вода/жир методом диксона
US20230366962A1 (en) Dixon-type water/fat separation mr imaging
Dietrich Single-shot pulse sequences
von Mengershausen et al. 3D diffusion tensor imaging with 2D navigated turbo spin echo
Pipe et al. SAR reduction for PROPELLER DWI at 3T
Pipe Pulse Sequences forDiffusion-Weighted MRI

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