CN113476031B - 反转恢复准备的3d梯度自旋回波扩散成像方法、介质及设备 - Google Patents
反转恢复准备的3d梯度自旋回波扩散成像方法、介质及设备 Download PDFInfo
- Publication number
- CN113476031B CN113476031B CN202110892182.8A CN202110892182A CN113476031B CN 113476031 B CN113476031 B CN 113476031B CN 202110892182 A CN202110892182 A CN 202110892182A CN 113476031 B CN113476031 B CN 113476031B
- Authority
- CN
- China
- Prior art keywords
- diffusion
- gradient
- module
- spin echo
- imaging
- 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
- 238000009792 diffusion process Methods 0.000 title claims abstract description 77
- 238000003384 imaging method Methods 0.000 title claims abstract description 43
- 238000011084 recovery Methods 0.000 title claims abstract description 11
- 238000000034 method Methods 0.000 claims abstract description 21
- 238000002360 preparation method Methods 0.000 claims abstract description 15
- 210000001175 cerebrospinal fluid Anatomy 0.000 claims abstract description 14
- 230000010355 oscillation Effects 0.000 claims abstract description 7
- 238000005516 engineering process Methods 0.000 claims abstract 3
- 230000005415 magnetization Effects 0.000 claims description 27
- 238000002075 inversion recovery Methods 0.000 claims description 21
- 238000012545 processing Methods 0.000 claims description 14
- 230000005284 excitation Effects 0.000 claims description 11
- 238000002595 magnetic resonance imaging Methods 0.000 claims description 8
- 238000004590 computer program Methods 0.000 claims description 6
- 239000003381 stabilizer Substances 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 5
- 238000000264 spin echo pulse sequence Methods 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 3
- 238000000926 separation method Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 9
- 238000005259 measurement Methods 0.000 abstract description 6
- 210000005013 brain tissue Anatomy 0.000 abstract description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 abstract description 4
- 230000009286 beneficial effect Effects 0.000 abstract description 2
- 210000004556 brain Anatomy 0.000 abstract description 2
- 230000009466 transformation Effects 0.000 abstract description 2
- 239000011159 matrix material Substances 0.000 abstract 1
- 230000036962 time dependent Effects 0.000 abstract 1
- 230000001419 dependent effect Effects 0.000 description 4
- 210000001519 tissue Anatomy 0.000 description 4
- 230000001054 cortical effect Effects 0.000 description 3
- 230000000971 hippocampal effect Effects 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 210000004884 grey matter Anatomy 0.000 description 2
- 210000001320 hippocampus Anatomy 0.000 description 2
- 210000001103 thalamus Anatomy 0.000 description 2
- 206010020649 Hyperkeratosis Diseases 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 210000000877 corpus callosum Anatomy 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000002597 diffusion-weighted imaging Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 230000002739 subcortical effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 210000004885 white matter Anatomy 0.000 description 1
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
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/0033—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room
- A61B5/004—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part
- A61B5/0042—Features or image-related aspects of imaging apparatus classified in A61B5/00, e.g. for MRI, optical tomography or impedance tomography apparatus; arrangements of imaging apparatus in a room adapted for image acquisition of a particular organ or body part for the brain
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2576/00—Medical imaging apparatus involving image processing or analysis
- A61B2576/02—Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part
- A61B2576/026—Medical imaging apparatus involving image processing or analysis specially adapted for a particular organ or body part for the brain
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Animal Behavior & Ethology (AREA)
- Veterinary Medicine (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Medical Informatics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Biophysics (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Neurology (AREA)
- High Energy & Nuclear Physics (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physiology (AREA)
- Psychiatry (AREA)
- Signal Processing (AREA)
- Computer Graphics (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Magnetic Resonance Imaging Apparatus (AREA)
Abstract
本发明公开一种反转恢复准备的3D梯度自旋回波扩散成像方法、介质及设备。该方法包括以下步骤：首先，在反转恢复准备模块中，施加180°反转射频脉冲，设置相应反转恢复时间以抑制脑脊液和自由水信号；其次，通过扩散编码模块，将一对梯形余弦振荡梯度或脉冲梯度嵌入到90°x‑180°y‑90°‑x射频脉冲中，将扩散编码与信号采集分开；然后，使用脂肪饱和模块对脂肪信号进行抑制；最后，采用3D梯度自旋回波的读出方式采集信号。本发明不仅改善了磁共振扩散成像的时间和信噪比，而且有利于脑组织扩散磁共振信号时间依赖性的精准测量，特别对于脑室附近由于部分容积效应受脑脊液影响的脑区效果明显，将有效推进时间依赖性扩散磁共振技术的临床转化。
Description
技术领域
本发明涉及磁共振技术领域，尤其涉及扩散磁共振成像领域。
背景技术
扩散磁共振(dMRI)通过检测被微观结构环境限制的水的扩散来揭示组织微观结构的各种细节。因此，dMRI测量的表观扩散系数(ADC)通常取决于水分子在周围环境中的扩散时间(td)。通过振荡梯度(OG)和脉冲梯度的结合可以测量一系列td范围(几毫秒到几十毫秒)的水分子扩散，并推断组织微结构，目前已经被广泛用于td依赖性dMRI的研究。
然而，由于临床MRI系统的梯度强度较低，因此需要较长的回波时间和重复时间才能达到合理的振荡频率和b值。这导致了访问高频振荡梯度的能力有限、成像时间较长、图像信噪比较弱等问题，所以该技术在人体临床应用上面临挑战。另一方面，我们观察到脑脊液(CSF)中自由水的ADC值大约是正常脑组织的3倍，因此在脑沟和脑室附近扫描时的部分容积效应和点扩展函数效应可能会导致脑组织ADC值受到脑脊液的影响，被高估约15-30％。因此，对于靠近脑室和脑沟的组织，如海马和皮质灰质，td依赖性测量可能不准确，由此估算出的组织微结构参数也受到影响。
发明内容
为了克服现有序列中存在的不足和时域弥散信号测量的问题，本发明提出了一种反转恢复准备的3D梯度自旋回波扩散成像序列(3D IR-OGprep-GRASE)，以解决上述提及的dMRI序列的问题，提高时域扩散成像的临床实用性。
为了达到上述目的，本发明采用以下技术方案予以实现：
第一方面，本发明提供了一种反转恢复准备的3D梯度自旋回波扩散成像方法，其包括：
S1：在成像序列的开始位置执行一个反转恢复模块，以抑制脑脊液的纵向磁化信号；
S2：在反转恢复模块之后执行扩散准备模块，将一对扩散梯度嵌入到90°x-180°y-90°-x射频脉冲中，以实现扩散编码与信号采集的分离；所述扩散梯度为梯形余弦振荡扩散梯度或脉冲扩散梯度；
S3：在扩散准备模块后执行脂肪饱和模块，抑制脂肪信号；
S4：在脂肪饱和模块后执行梯度自旋回波模块，使用3D梯度回波-自旋回波序列进行信号采集，得到3D K-空间数据；
S5：在梯度自旋回波模块后执行通用自动校准并行采集重建模块，对3D K-空间数据进行重建得到完整图像。
基于该方案，各步骤还可以进一步提供以下优选的实现方式。需要注意的是，各优选方式中的技术特征在没有冲突的情况下均可进行相互组合。当然这些优选方式也可以通过其他能够实现相同技术效果的方式实现，不构成限制。
作为优选，所述反转恢复模块中的处理流程如下：
先对整个3D成像空间施加非空间选择性的180°翻转角的射频脉冲激励，再等待反转时间(TI)，使得脑脊液的纵向磁化矢量完全消除。
作为优选，所述扩散准备模块中的处理流程如下：
首先，在TI时刻沿X轴方向进行90°翻转角的射频激发脉冲；然后，沿着预定的扩散方向施加所述扩散梯度；接着，沿Y轴方向进行180°翻转角的射频回聚脉冲，然后再施加完全相同的所述扩散梯度；再后，沿Z轴方向施加稳定器梯度(Gstb)；最后，沿-X轴方向进行90°翻转角的射频激发脉冲，将横向磁化矢量转化为纵向磁化矢量。
作为优选，所述脂肪饱和模块中的处理流程如下：
施加三个95°翻转角的频率选择性的高斯脉冲，将脂肪纵向磁化矢量的sin(95°)部分进行翻转并等待T1*ln(2)时间后将其置0，使脂肪的大部分纵向磁化矢量转化为横向磁化矢量，同时施加梯度，将脂肪的横向磁化矢量完全消除；其中T1为脂肪的纵向弛豫时间。
作为优选，所述梯度自旋回波序列(GRASE)模块中的处理流程如下：
通过在Y方向上进行回波平面成像(EPI)编码以及在Z方向上进行涡轮自旋回波(TSE)编码，从而实现3D GRASE读出；其中沿EPI方向执行并行成像采集，而沿TSE方向上执行中心编码；在每段EPI读出时刻的前后沿Z轴方向施加稳定器梯度。
作为优选，所述通用自动校准并行采集(GRAPPA)重建模块中的处理流程如下：
首先使用3-navigator-echo技术来估计EPI奇偶行之间的相位变化，然后使用GRAPPA技术将K空间中沿EPI方向执行的并行成像采集的空缺行填充完整。
第二方面，本发明提供了一种计算机可读存储介质，所述存储介质上存储有计算机程序，当所述计算机程序被处理器执行时，用于实现如第一方面中任一所述的反转恢复准备的3D梯度自旋回波扩散成像方法。
第三方面，本发明提供了一种磁共振成像设备，其特征在于，包括磁共振扫描器以及控制单元，所述控制单元中存储有成像序列，所述成像序列被执行时，用于实现如第一方面中任一所述的反转恢复准备的3D梯度自旋回波扩散成像方法。
相对于现有技术，本发明具有以下有益效果：
本发明提出了一种反转恢复准备的3D梯度自旋回波扩散成像方法，首次将反转恢复序列、扩散编码、3D成像序列结合起来。与传统的2D-EPI采集相比，本发明的序列不仅缩短了振荡梯度扩散成像的扫描时间而且提高了信噪比，而且反转恢复模块抑制了脑脊液对脑室、脑沟附近的ADC值。实验表明只有使用具有反转恢复准备的3D IR-OGprep-GRASE序列，才能准确测量这些区域的扩散信号的td依赖性。
附图说明
图1是3D IR-OGprep-GRASE序列的时序图。
图2是对分别使用3D OGprep-GRASE和3D IR-OGprep-GRASE序列获取得到的b0和ADC图，并显示了7个脑区的感兴趣区域。
图3是对分别使用3D OGprep-GRASE和3D IR-OGprep-GRASE序列在6名被试上获取得到的ADC值在各个感兴趣区域的对比结果。
具体实施方式
下面基于本发明提出的方法结合实施例展示其具体的技术效果，以便本领域技术人员更好地理解本发明的实质。
在本发明的一种较优实现方式中，上述反转恢复准备的3D梯度自旋回波扩散成像方法，其对应于一种反转恢复准备的3D梯度自旋回波扩散成像(3DIR-OGprep-GRASE)序列，如图1是3D IR-OGprep-GRASE的序列时序图。该成像方法具体包括以下步骤：
步骤一：在成像序列时序开始时先执行一个反转恢复模块，以抑制脑脊液的纵向磁化信号。在本实施例中，反转恢复模块中的处理流程如下：首先对整个3D成像空间施加非空间选择性的180°翻转角的射频脉冲激励，在施加完射频脉冲后等待一段被称为反转时间(TI)的时间使脑脊液信号的纵向磁化矢量完全消失。此时进行信号采集或再次激发，可以认为已经将脑脊液信号的纵向磁化完全消除而脑组织的纵向磁化依然存在。
步骤二：在反转恢复模块之后紧接着执行扩散准备模块，将一对扩散梯度嵌入到90°x-180°y-90°-x射频脉冲中，以实现扩散编码与信号采集的分离，其中扩散梯度可以是梯形余弦振荡扩散梯度或脉冲扩散梯度。在本实施例中，该扩散准备模块中的具体处理流程如下：首先，在TI时刻即前述反转时间结束的时刻，沿X轴方向施加90°翻转角的非层选性的硬脉冲激励(作为射频激发脉冲)；然后，沿着设置扩散的方向施加一个扩散梯度，该扩散梯度可以是一定频率梯形余弦振荡扩散梯度或一定扩散间隔的脉冲扩散梯度，梯度的幅值大小取决于预先设定的b值；接着，沿Y轴方向进行180°翻转角的绝热正切脉冲激励(作为射频回聚脉冲)，使得散向的横向磁化矢量回聚；然后，再施加与之前施加的幅值、形状、持续时间等参数完全相同的扩散梯度；之后，沿Z轴方向施加稳定器梯度(Gstb)，以减轻与相位误差有关的信号调制，并将梯度幅值设置为与之后的GRASE读出模块中的破碎机梯度(crusher)相同的幅度；最后，沿-X轴方向发射90°翻转角的非层选性的硬脉冲激励，将横向磁化矢量转化为纵向磁化矢量。扩散准备模块中的回波时间表示为TE1。
步骤三：在扩散准备模块后执行脂肪饱和模块，抑制脂肪信号。在本实施例中，脂肪饱和模块中的处理流程如下：施加三个95°翻转角的频率选择性的高斯脉冲，将脂肪纵向磁化矢量的sin(95°)部分进行翻转并等待T1*ln(2)时间后将其置0，使脂肪的大部分纵向磁化矢量转化为横向磁化矢量，同时施加梯度，将脂肪的横向磁化矢量完全消除；其中T1为脂肪的纵向弛豫时间。
步骤四：在脂肪饱和模块后执行梯度自旋回波(GRASE)模块，使用3D梯度回波-自旋回波序列进行信号采集，得到3D K-空间数据。在本实施例中，梯度自旋回波模块的处理流程如下：通过在Y方向上进行回波平面成像(EPI)编码以及在Z方向上进行涡轮自旋回波(TSE)编码来实现3D-GRASE读出。在EPI方向上的相位编码的数目被表示为NEPI，而在TSE方向上的相位编码的数目被表示为NSE。并行采集是沿EPI方向对NEPI进行间隔加速因子值(Na)的行数进行采集，而沿TSE方向上执行中心编码，即第一次回聚脉冲后的EPI读取信号的方式在kz方向填充k-space中心。稳定器梯度沿Z轴方向在每段EPI读出时刻的前后与回聚脉冲旁边的破碎机梯度同时施加。GRASE模块中的回波时间表示为TE2，它等于之后以TSE方式进行信号采集的回波间隔(ESP)。因此，所获取信号的总回波时间为TE1+TE2。
当沿EPI方向对NEPI进行并行成像采集时，相对于常规采集方式，不同之处在于NEPI减少为原来的1/Na，采集时间降低从而导致TE减小，因此有效地提高了图像的信噪比。
步骤五：如果使用Na>1，则需要在梯度自旋回波模块后执行通用自动校准并行采集(GRAPPA)重建模块，来填充未采集行的K空间数据，对3D K-空间数据进行重建得到完整图像。在本实施例中，通用自动校准并行采集重建模块中的处理流程如下：重建首先使用3-navigator-echo技术来估计EPI奇偶行之间的相位变化，对所有采集到的EPI奇偶行信号进行相位进行校正。然后使用GRAPPA技术填充未采集行的K空间数据，由此将K空间中沿EPI方向执行的并行成像采集的空缺行填充完整。需要说明的是，如果使用的Na＝1，则无需执行GRAPPA重建模块。与常规的采集方法相比，通用自动校准部分并行采集方法大大提高了采集效率，不仅缩短了采集时间，而且降低了TE，因此可以得到更高信噪比的图像。
由此可见，上述上述步骤一到五中的各模块本质上相当于3DIR-OGprep-GRASE序列中的不同部分。为了进一步展示该做法的效果，本发明还进行了下列效果验证步骤：
步骤六：通过步骤五处理得到对应的图像后，再进行计算得到扩散加权成像的量化指标，如表观扩散系数(ADC)等。
下面基于上述方法的步骤一到步骤六，结合实施例对其技术效果进行展示，以便本领域技术人员更好地理解本发明的实质。
实施例
将上述的反转恢复准备的3D梯度自旋回波扩散成像方法即对应的3DIR-OGprep-GRASE序列以及没有反转恢复模块的3D梯度自旋回波扩散成像方法即对应的3D OGprep-GRASE序列在6名健康的年轻男性志愿者中进行了测试，下面介绍此处的具体参数：磁共振扫描使用了西门子Prisma 3T扫描仪(最大梯度为80mT/m、最大切换率200mT/m)，所有扫描均使用64通道头部线圈。
实验：为了比较3D OGprep-GRASE和3D IR-OGprep-GRASE序列对脑脊液附近区域扩散信息的td依赖性测量的准确性，使用脉冲扩散梯度(0Hz)以及20Hz、40Hz、60Hz的振荡扩散梯执行扫描，其他成像参数一致：扩散加权＝420s/mm2、6个方向，2次重复，扫描视野＝220×200mm，分辨率为2×2mm2，切片厚度为4mm。3D-GRASE读出部分采用以下方案：
单次激发的GRASE读出，NEPI＝39，NSE＝10，Na＝2，turbo因子＝12，层面方向过采样＝20％，部分傅立叶因子＝76％，带宽＝2094Hz，TE1/TE2/TR/TI＝103.46/31.94/9000/2500ms，扫描时间＝3min。3D OGprep-GRASE序列除去除了IR模块,其余参数与3D IR-OGprep-GRASE序列完全一致。
图2显示了使用3D IR-OGprep-GRASE和OGprep-GRASE序列分别获得的b0和ADC图。在b0图像上勾画出了感兴趣区域(ROI)，包括海马头部(HPH)，海马体部(HPB)，海马尾部(HPT)，皮层灰质(CGM)，丘脑(THA)，皮层下白质(PWM)和胼胝体压部(SCC)。
图3显示了使用3D IR-OGprep-GRASE和OGprep-GRASE序列在不同频率(60Hz、40Hz、20Hz、0Hz)下测量扩散时间(td)依赖的ADC。在3DIR-OGprep-GRASE序列中，所有ROI均显示出ADC的td依赖性变化。而在无IR模块的3D OGprep-GRASE序列中，靠近脑室和脑沟的区域，例如：海马、皮质灰质和胼胝体压部等区域，未见td依赖效应。
另外，在其他实施例中，还可以提供一种计算机可读存储介质，所述存储介质上存储有计算机程序，当所述计算机程序被处理器执行时，用于实现如上述步骤一到步骤五所述的反转恢复准备的3D梯度自旋回波扩散成像方法。
需要注意的是，上述的存储器可以包括随机存取存储器(Random Access Memory，RAM)，也可以包括非易失性存储器(Non-Volatile Memory，NVM)，例如至少一个磁盘存储器。上述的处理器可以是通用处理器，包括中央处理器(Central Processing Unit，CPU)、网络处理器(Network Processor，NP)等；还可以是数字信号处理器(Digital SignalProcessing，DSP)、专用集成电路(Application Specific Integrated Circuit，ASIC)、现场可编程门阵列(Field-Programmable Gate Array，FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。当然，还装置中还应当具有实现程序运行的必要组件，例如电源、通信总线等等。
同样的，在其他实施例中，还可以提供一种磁共振成像设备，其包括磁共振扫描器以及控制单元，所述控制单元中存储有成像序列即前述的3DIR-OGprep-GRASE序列，该成像序列被执行时，用于实现如上述步骤一到步骤五所述的反转恢复准备的3D梯度自旋回波扩散成像方法。
上述的磁共振扫描器和控制单元的硬件结构可采用常规磁共振成像系统实现。磁共振扫描器应当包括磁体部分和磁共振波谱仪部分，控制单元中应当含有实现磁共振成像必要的数据处理和图像重建软件，3D IR-OGprep-GRASE序列可由控制单元读取并用于控制磁共振扫描器。也就是说，3D IR-OGprep-GRASE序列可以直接应用至相应的磁共振成像设备中。
需要指出的是，以上所述的实施例只是本发明的一种较佳的方案，然其并非用以限制本发明。有关技术领域的普通技术人员，在不脱离本发明的精神和范围的情况下，还可以做出各种变化和变型。因此凡采取等同替换或等效变换的方式所获得的技术方案，均落在本发明的保护范围内。
Claims (4)
1.一种反转恢复准备的3D梯度自旋回波扩散成像方法，其特征在于，包括：
S1：在成像序列的开始位置执行一个反转恢复模块，以抑制脑脊液的纵向磁化信号；
S2：在反转恢复模块之后执行扩散准备模块，将一对扩散梯度嵌入到90°x-180°y-90°-x射频脉冲中，以实现扩散编码与信号采集的分离；所述扩散梯度为梯形余弦振荡扩散梯度或脉冲扩散梯度；
S3：在扩散准备模块后执行脂肪饱和模块，抑制脂肪信号；
S4：在脂肪饱和模块后执行梯度自旋回波模块，使用3D梯度回波-自旋回波序列进行信号采集，得到3D K-空间数据；
S5：在梯度自旋回波模块后执行通用自动校准并行采集重建模块，对3D K-空间数据进行重建得到完整图像；
所述脂肪饱和模块中的处理流程如下：
施加三个95°翻转角的频率选择性的高斯脉冲，将脂肪纵向磁化矢量的sin(95°)部分进行翻转并等待T1*ln(2)时间后将其置0，使脂肪的大部分纵向磁化矢量转化为横向磁化矢量，同时施加梯度，将脂肪的横向磁化矢量完全消除；其中T1为脂肪的纵向弛豫时间；
所述通用自动校准并行采集GRAPPA重建模块中的处理流程如下：
首先使用3-navigator-echo技术来估计EPI奇偶行之间的相位变化，然后使用GRAPPA技术将K空间中沿EPI方向执行的并行成像采集的空缺行填充完整；
所述反转恢复模块中的处理流程如下：
先对整个3D成像空间施加非空间选择性的180°翻转角的射频脉冲激励，再等待反转时间TI，使得脑脊液的纵向磁化矢量完全消除；
所述梯度自旋回波GRASE模块中的处理流程如下：
通过在Y方向上进行回波平面成像EPI编码以及在Z方向上进行涡轮自旋回波TSE编码，从而实现3D GRASE读出；其中沿EPI方向执行并行成像采集，而沿TSE方向上执行中心编码；在每段EPI读出时刻的前后沿Z轴方向施加稳定器梯度。
2.如权利要求1所述的反转恢复准备的3D梯度自旋回波扩散成像方法，其特征在于，所述扩散准备模块中的处理流程如下：
首先，在TI时刻沿X轴方向进行90°翻转角的射频激发脉冲；然后，沿着预定的扩散方向施加所述扩散梯度；接着，沿Y轴方向进行180°翻转角的射频回聚脉冲，然后再施加完全相同的所述扩散梯度；再后，沿Z轴方向施加稳定器梯度Gstb；最后，沿-X轴方向进行90°翻转角的射频激发脉冲，将横向磁化矢量转化为纵向磁化矢量。
3.一种计算机可读存储介质，其特征在于，所述存储介质上存储有计算机程序，当所述计算机程序被处理器执行时，用于实现如权利要求1或2所述的反转恢复准备的3D梯度自旋回波扩散成像方法。
4.一种磁共振成像设备，其特征在于，包括磁共振扫描器以及控制单元，所述控制单元中存储有成像序列，所述成像序列被执行时，用于实现如权利要求1或2所述的反转恢复准备的3D梯度自旋回波扩散成像方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110892182.8A CN113476031B (zh) | 2021-08-04 | 2021-08-04 | 反转恢复准备的3d梯度自旋回波扩散成像方法、介质及设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110892182.8A CN113476031B (zh) | 2021-08-04 | 2021-08-04 | 反转恢复准备的3d梯度自旋回波扩散成像方法、介质及设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113476031A CN113476031A (zh) | 2021-10-08 |
| CN113476031B true CN113476031B (zh) | 2024-03-29 |
Family
ID=77945497
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110892182.8A Active CN113476031B (zh) | 2021-08-04 | 2021-08-04 | 反转恢复准备的3d梯度自旋回波扩散成像方法、介质及设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113476031B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115184851A (zh) * | 2022-07-07 | 2022-10-14 | 北京万东医疗科技股份有限公司 | 磁共振成像方法、装置、设备及存储介质 |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000023937A (ja) * | 1998-07-13 | 2000-01-25 | Toshiba Corp | Ｍｒｉ装置およびｍｒイメージング方法 |
| CN101259018A (zh) * | 2006-12-07 | 2008-09-10 | Ge医疗系统环球技术有限公司 | 磁共振成像设备和磁共振成像方法 |
| CN102772206A (zh) * | 2011-05-12 | 2012-11-14 | 上海联影医疗科技有限公司 | 磁共振成像方法 |
| CN103961097A (zh) * | 2013-02-04 | 2014-08-06 | 上海联影医疗科技有限公司 | 一种磁共振扫描短te成像方法及磁共振扫描系统 |
| CN104856677A (zh) * | 2015-06-14 | 2015-08-26 | 西安交通大学医学院第一附属医院 | 磁化转移联合层面选择反转恢复预脉冲的磁共振成像方法 |
| CN105829906A (zh) * | 2013-12-19 | 2016-08-03 | 皇家飞利浦有限公司 | 具有水/脂肪分离的相位敏感的反转恢复mri |
| CN109085523A (zh) * | 2018-07-23 | 2018-12-25 | 中国科学院深圳先进技术研究院 | 血管壁成像中脑脊液信号的抑制方法、装置、设备及介质 |
| CN111352054A (zh) * | 2020-03-31 | 2020-06-30 | 浙江大学 | 一种振荡梯度准备的3d梯度自旋回波成像方法及设备 |
-
2021
- 2021-08-04 CN CN202110892182.8A patent/CN113476031B/zh active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2000023937A (ja) * | 1998-07-13 | 2000-01-25 | Toshiba Corp | Ｍｒｉ装置およびｍｒイメージング方法 |
| CN101259018A (zh) * | 2006-12-07 | 2008-09-10 | Ge医疗系统环球技术有限公司 | 磁共振成像设备和磁共振成像方法 |
| CN102772206A (zh) * | 2011-05-12 | 2012-11-14 | 上海联影医疗科技有限公司 | 磁共振成像方法 |
| CN103961097A (zh) * | 2013-02-04 | 2014-08-06 | 上海联影医疗科技有限公司 | 一种磁共振扫描短te成像方法及磁共振扫描系统 |
| CN105829906A (zh) * | 2013-12-19 | 2016-08-03 | 皇家飞利浦有限公司 | 具有水/脂肪分离的相位敏感的反转恢复mri |
| CN104856677A (zh) * | 2015-06-14 | 2015-08-26 | 西安交通大学医学院第一附属医院 | 磁化转移联合层面选择反转恢复预脉冲的磁共振成像方法 |
| CN109085523A (zh) * | 2018-07-23 | 2018-12-25 | 中国科学院深圳先进技术研究院 | 血管壁成像中脑脊液信号的抑制方法、装置、设备及介质 |
| CN111352054A (zh) * | 2020-03-31 | 2020-06-30 | 浙江大学 | 一种振荡梯度准备的3d梯度自旋回波成像方法及设备 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113476031A (zh) | 2021-10-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Hargreaves et al. | Variable‐rate selective excitation for rapid MRI sequences | |
| US7639010B2 (en) | Turbospin echo imaging sequence with long echo trains and optimized T1 contrast | |
| US9329252B2 (en) | Apparatus for real-time phase correction for diffusion-weighted magnetic resonance imaging using adaptive RF pulses | |
| EP2702423B1 (en) | Method for r2* quantification with mri with correction for macroscopic magnetic field inhomogeneities | |
| US10288703B2 (en) | MRI method of hybrid acquisition in 3D TSE | |
| US7015696B2 (en) | Magnetic resonance imaging apparatus and magnetic resonance imaging method | |
| US8072214B2 (en) | Method for magnetic resonance imaging using sequential excitation and refocusing | |
| US20090267604A1 (en) | Magnetic resonance imaging apparatus and magnetic resonance imaging method | |
| US20170350951A1 (en) | System and method for chemical shift magnetic resonance imaging using magnetization transfer | |
| EP1506417A1 (en) | Reduction of susceptibility artifacts in subencoded single-shot magnetic resonance imaging | |
| CN111352054B (zh) | 一种振荡梯度准备的3d梯度自旋回波成像方法及设备 | |
| KR101694610B1 (ko) | 자기 공명 장치를 이용하여 파라미터 맵을 획득하기 위한 방법 및 자기 공명 장치 | |
| US9645211B2 (en) | Magnetic resonance imaging using steering-propeller | |
| US20160170000A1 (en) | Method and apparatus for acquiring magnetic resonance data and entering the data into k-space | |
| US20090129648A1 (en) | Method of reducing imaging time in propeller-MRI by under-sampling and iterative image reconstruction | |
| CN113476031B (zh) | 反转恢复准备的3d梯度自旋回波扩散成像方法、介质及设备 | |
| US6586933B1 (en) | Method and system for MRI with lipid suppression | |
| US5459401A (en) | MRI method for producing images having weak through medium T2 weighing employing a turbo-spin echo sequence | |
| US11016162B2 (en) | Method and magnetic resonance apparatus for generating a combination image dataset | |
| Finsterbusch et al. | Diffusion‐weighted single‐shot line scan imaging of the human brain | |
| EP3102103B1 (en) | Magnetic resonance imaging method with t2* mapping based on echo planar imaging | |
| Maier | Slab scan diffusion imaging | |
| US11333733B2 (en) | Method and apparatus for recording a magnetic resonance data record | |
| CN110873856B (zh) | 确定最佳磁共振成像扫描嵌套方式的方法及装置 | |
| US20060074291A1 (en) | Magnetic resonance imaging system and method |
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 |