CN109100669B - 基于重叠回波的单扫描同步磁共振扩散及t2成像方法 - Google Patents

基于重叠回波的单扫描同步磁共振扩散及t2成像方法 Download PDF

Info

Publication number
CN109100669B
CN109100669B CN201810765276.7A CN201810765276A CN109100669B CN 109100669 B CN109100669 B CN 109100669B CN 201810765276 A CN201810765276 A CN 201810765276A CN 109100669 B CN109100669 B CN 109100669B
Authority
CN
China
Prior art keywords
echo
gradient
sampling
refocusing
pulse
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
CN201810765276.7A
Other languages
English (en)
Other versions
CN109100669A (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.)
Xiamen University
Original Assignee
Xiamen University
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 Xiamen University filed Critical Xiamen University
Priority to CN201810765276.7A priority Critical patent/CN109100669B/zh
Publication of CN109100669A publication Critical patent/CN109100669A/zh
Application granted granted Critical
Publication of CN109100669B publication Critical patent/CN109100669B/zh
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/50—NMR imaging systems based on the determination of relaxation times, e.g. T1 measurement by IR sequences; T2 measurement by multiple-echo sequences
    • 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
    • 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

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • High Energy & Nuclear Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • Radiology & Medical Imaging (AREA)
  • General Health & Medical Sciences (AREA)
  • Pathology (AREA)
  • Biophysics (AREA)
  • Vascular Medicine (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

基于重叠回波的单扫描同步磁共振扩散及T2成像方法，涉及磁共振成像。用两个相同偏转角的小角度激发脉冲和两个回波链采样产生四个回波，每个回波链采集两个回波信号。在第一个激发脉冲之后有一段演化时间和一对扩散梯度，使得第一次采样中的两个回波信号的横向弛豫时间和扩散加权不同。每个激发脉冲之后都加一个频率维和相位维的移位梯度使得不同的激发脉冲产生的信号在k空间的位置不一样。在第一次采样结束后，用一个重聚脉冲重聚两个回波信号，然后进行第二次采样，在第二次采样中可获得两个具有相同横向弛豫时间和不同扩散加权的回波信号。将采样信号用深度学习进重建得到定量T2和ADC图像，可在单次扫描中获得T2和ADC图像。

Description

基于重叠回波的单扫描同步磁共振扩散及T2成像方法
技术领域
本发明涉及磁共振成像，尤其是涉及基于重叠回波的单扫描同步磁共振扩散(apparent diffusion coefficient，ADC)及T2成像方法。
背景技术
磁共振定量成像可以提供用于以非入侵的方式分析正常和病态活体组织的对比机制，目前已应用于脑、心脏、脊髓微细结构的研究[1,2]。定量成像可以消除独立于组织性质的影响，比如对操作者依赖性、扫描参数差异、磁场空间变化和图像缩放等[3]。多参数的定量成像在临床磁共振成像(magnetic resonance imaging,MRI)中获得越来越多的关注。多参数的定量成像通过采集不同参数的定量图(如T2和扩散图)，并用这些图对感兴趣区域(region of interest，ROI)进行联合分析，从而得到更多更准确的关于人体生理状态的变化或者神经结构的信息[4-6]。T2图需要通过拟合多个不同回波时间(echo time,TE)的自旋回波MRI数据得到，扩散图需要拟合多个不同扩散加权的MRI数据得到，因而一般获取T2和扩散图的时间很长。较长的采样时间也使得T2和扩散图像容易受到运动伪影的影响[7-9]。平面回波成像(echo-planar imaging，EPI)，可以作为单扫描快速成像方法用于T2和扩散成像，但也至少需要两次的EPI 采样才能得到T2或扩散图，病人在扫描间的运动或生理状态的变化仍会导致伪影的出现[10]。现有的用于估计和校正运动伪影的算法无法校正层间运动造成的伪影和运动造成的图像模糊[11,12]。已经提出的通过重叠回波分离(overlapping-echo detachment，OLED)平面成像方法能够在单次扫描中获得高质量的T2图像[13]，其时空分辨率与传统的单扫描EPI图像相当。除了 T2成像，OLED还被用于单扫描扩散成像[14]。OLED的T2和扩散图还显示出对运动伪影和非理想B1场的更强的抵抗力。但传统的OLED图像重建方法在高T2或高扩散值的区域，如脑脊液(cerebrospinal fluid，CSF)，往往会失真，图像的空间分辨率也十分有限。前不久，一种基于深度学习卷积神经网络的OLED T2成像重建方法被提出[15]，该方法重建出的OLED T2图比传统方法得到的OLED T2图与更高的分辨率，提高CSF区域的精度。
参考文献：
1.Zhao B,Lam F,Liang ZP.Model-based MR parameter mapping withsparsity constraints: parameter estimation andperformance bounds.IEEE T MedImaging 2014；33:1832-1844.
2.Garces P,Pereda E,Hernandez-Tamames JA,Del-Pozo F,Maestu F,Pineda-Pardo JA. Multimodal description ofwhole brain connectivity:a comparisonofresting state MEG,fMRI,and DWI.Hum Brain Mapp 2016；37:20-34.
3.Margaret Cheng HL,Stikov N,Ghugre NR,Wright GA.Practical medicalapplications of quantitative MR relaxometry.J Magn Reson Imaging 2012；36:805-824.
4.Calamante F,Jeurissen B,Smith RE,Tournier JD,Connelly A.The role ofwhole-brain diffusion MRI as a tool for studying human in vivo corticalsegregation based on a measure of neurite density.Magn Reson Med 2018；79:2738-2744.
5.Aliotta E,Moulin K,Zhang Z,Ennis DB.Simultaneous measurement of T2and apparent diffusion coefficient(T2+ADC)in the heart with motion-compensated spin echo diffusion-weighted imaging.Magn Reson Med 2018；79:654-662.
6.Staroswiecki E,Granlund KL,Alley MT,Gold GE,HargreavesBA.Simultaneous estimation of T(2)and apparent diffusion coefficient in humanarticular cartilage in vivo with a modified three-dimensional double echosteady state(DESS)sequence at 3 T.Magn Reson Med 2012；67:1086-1096.
7.Le Bihan D,Poupon C,Amadon A,Lethimonnier F.Artifacts and pitfallsin diffusion MRI. J Magn Reson Imaging 2006；24:478-488.
8.Marami B,Scherrer B,Afacan O,Erem B,Warfield SK,Gholipour A.Motion-robust diffusion-weighted brain MRI reconstruction through slice-levelregistration-based motionc tracking.IEEE Trans Med Imaging 2016；35:2258-2269.
9.Arthurs OJ,Edwards A,Austin T,Graves MJ,Lomas DJ.The challenges ofneonatal magnetic resonance imaging.Pediatr Radiol 2012；42:1183-1194.
10.Shi X,Kholmovski EG,Kim SE,Parker DL,Jeong EK.Improvement ofaccuracy of diffusion MRI using real-time self-gated data acquisition.NMRBiomed 2009；22:545-550.
11.Alhamud A,Taylor PA,Laughton B,van der Kouwe AJ,Meintjes EM.Motionartifact reduction in pediatric diffusion tensor imaging using fastprospective correction.J Magn Reson Imaging 2015；41:1353-1364.
12.Brockstedt S,Moore JR,Thomsen C,
Figure BDA0001728884740000021
S,
Figure BDA0001728884740000022
F.High-resolutiondiffusion imaging using phase-corrected segmented echo-planar imaging.MagnReson Imaging 2000；18:649-657.
13.Cai C,Zeng Y,Zhuang Y,Cai S,Chen L,Ding X,Bao L,Zhong J,ChenZ.Single-shot T2 mapping through Overlapping-echo Detachment(OLED)PlanarImaging.IEEE Trans Biomed Eng 2017；64:2450–2461.
14.Ma L,Cai C,Yang H,Cai S,Qian J,Xiao L,Zhong K,Chen Z.Motion-tolerant diffusion mapping based on single-shot overlapping-echo detachment(OLED)planar imaging.Magn Reson Med 2018；80:200–210.
15.Cai C,Wang C,Zeng YQ,Cai SH,Liang D,Wu YW,Chen Z,Ding XH and ZhongJH, Single-Shot T2Mapping using OverLapping-Echo Detachment Planar Imagingand a Deep ConvolutionalNeural Network,Magn Reson in Med,2017.DOI:10.1002/mrm.27205.
发明内容
本发明的目的在于提供将OLED序列进行改进，实现单次扫描可以同时获得四个回波不同扩散和T2加权的信号，然后利用深度学习技术从四个重叠的回波信号重建出可靠T2和扩散图的基于重叠回波的单扫描同步磁共振扩散及T2成像方法。
本发明包括以下步骤：
1)在磁共振成像仪操作台上，打开磁共振成像仪中相应的操作软件，首先对成像物体进行感兴趣区域定位，然后进行调谐、匀场、功率和频率校正；
2)导入事先编译好的DT2M-OLED成像序列，根据具体的实验情况，设置脉冲序列的各个参数；
在步骤2)中，所述DT2M-OLED成像序列的结构可依次为：第1翻转角为α的激发脉冲、脉冲间距δTE、第2翻转角为α的激发脉冲、第2移位梯度、第1翻转角为β的重聚脉冲、第一个采样回波链、第2翻转角为β的重聚脉冲、第二个采样回波链；所述脉冲间距δTE中包含第1移位梯度和扩散梯度，第1移位梯度包含沿频率维(x方向)的梯度Gro1和沿相位维 (y方向)的梯度Gpe1；第2移位梯度包含沿频率维(x方向)的梯度Gro2和沿相位维(y方向)的梯度Gpe2。
所述第1翻转角和第2翻转角为α的激发脉冲结合第1移位梯度和第2移位梯度，使两个回波在k空间的中心产生偏移，重聚脉冲以及两个小角度激发脉冲都与层选方向(z方向) 的层选梯度Gss相结合进行层选；第2翻转角为α的激发脉冲前后分别施加回波延时，其中δTE的长度为1/3回波链长，重聚脉冲前后有x，y，z三个方向的破坏梯度作用。
所述第一个采样回波链和第二个采样回波链是由分别作用在x，y方向的梯度链组成；x 方向的梯度链由一系列正负梯度构成，且每个梯度的面积是所述移位梯度Gro1的三倍；y方向的梯度链是由一系列大小相等的尖峰梯度构成，且尖峰梯度的面积和等于移位梯度Gpe1面积的三倍。
在第一个采样回波链之前，频率和相位方向分别施加了重聚梯度，频率维的重聚梯度的面积是频率维采样梯度Gro梯度面积的一半，方向与Gro相反；相位维的重聚梯度的面积是所述所有尖峰梯度面积的一半，方向与尖峰梯度相反；第一个采样回波链的信号包含第一个回波和第二个回波的信号；
在第一个采样回波链之后，利用重聚脉冲产生两个新的自旋回波，并利用第二个采样回波链采样，第二个采样回波链的信号包含第三个回波和第四个回波的信号。
3)执行步骤2)设置好的DT2M-OLED成像序列，进行数据采样，数据采样完成后，得到两个回波链的信号；
4)对步骤3)得到的两个回波链的信号进行归一化、充零和快速傅里叶变换，将k空间的信号转换到图像域，再将图像域的信号用基于卷积神经网络的方法重建，得实验数据；
5)根据实验样品的特征生成随机的模板，对模板进行模拟采样得到k空间信号，然后对 k空间信号进行归一化、充零和快速傅里叶变换，得到训练数据；
6)使用tensorflow深度学习框架和python搭建卷积神经网络模型，设置好训练的相关参数，得数据训练网络；
7)使用步骤6)得到的数据训练网络，直至数据训练网络收敛并达到稳定，得到训练好的网络模型，然后对步骤4)得到的实验数据进行重建，得到可靠的ADC和T2图像。
在步骤3)中，所述DT2M-OLED成像序列采集到的两个回波链的信号进行分析，并对回波信号磁化矢量演化进行理论推导，在第一个采样回波链采样期，信号表达式如下：
Figure BDA0001728884740000041
式中
Figure BDA0001728884740000042
是扩散因子。
在第二个采样回波链采样期，信号表达式如下：
通过实验发现，当α＝45°，β＝180°时，四个回波信号的强度都相对较高。
在步骤5)中，所述随机的模板是根据实验样本的特征分布使用计算机批量随机生成，同时保证模板的复杂性较高能够包含实验样本的所有特征；在模拟采样过程中，考虑到真实实验环境变化，加入不稳定因素，提高网络模型对不理想实验环境的鲁棒性；所述不稳定因素包括激发脉冲角度偏差、移位梯度偏差和噪声等。
在步骤6)中，所述使用tensorflow深度学习框架和python搭建卷积神经网络模型包括：网络的主体结构以及相关的训练参数，网络模型的目标函数为：
Figure BDA0001728884740000051
其中，N是同一批训练图像的数量，f(·)是训练网络，W和b是网络参数，xk是第k个回波链信号经步骤4)后得到的图像，y是输入图像对应的模板，ychange是对y中小于0.06的值置为0.06之后的矩阵，ymask是对y使用canny算子求得的图像边缘信息，λ表示正则项参数，▽代表梯度算子。
在步骤7)中，所述使用步骤6)得到的数据训练网络由于使用随机的模板进行训练，泛化性较强，能适用于多种样品的重建。
本发明用两个相同偏转角的小角度激发脉冲和两个回波链采样产生四个回波，每个回波链采集两个回波信号。在第一个激发脉冲之后有一段演化时间和一对扩散梯度，使得第一次采样中的两个回波信号的横向弛豫时间和扩散加权不同。同时，在每个激发脉冲之后都加一个频率维和相位维的移位梯度使得不同的激发脉冲产生的信号在k空间的位置不一样。在第一次采样结束后，用一个重聚脉冲重聚两个回波信号，然后进行第二次采样，在第二次采样中可以获得了两个具有相同横向弛豫时间和不同扩散加权的回波信号。然后将采样信号用深度学习的方法进重建得到定量T2图像和ADC图像。该方法可以在单次扫描中获得可靠的T2图像和ADC图像。
附图说明
图1为本发明采用的DT2M-OLED成像序列图。
图2为展示了模型的ADC和T2重建结果；(a)是用EPI序列采模型的结构图，其中的5个用红色圆圈标出的区域是感兴趣区域(region ofinterest，ROI)，分别用数字1～5编号。(b) 是用DT2M-OLED序列采集第一个采样回波链的原始幅值图。(c)是用EPI拟合出来的ADC图像。(d)是DT2M-OLED数据重建得到的ADC图像。(e)是用EPI拟合出来的ADC图像。
图3为展示5个ROI的T2统计结果。
图4为展示5个ROI的ADC统计结果。
具体实施方式
以下实施例将结合附图对本发明作进一步说明。
本发明提供了基于重叠回波的单扫描同步磁共振扩散及T2成像方法，具体实施过程中的各个步骤如下：
(1)在磁共振成像仪操作台上，打开磁共振成像仪中相应的操作软件，首先对成像物体进行感兴趣区域定位，然后进行调谐、匀场、功率和频率校正；
(2)导入事先编译好的DT2M-OLED成像序列：根据具体的实验情况，设置脉冲序列的各个参数。
所述DT2M-OLED成像序列的结构如图1依次为：第1翻转角为α的激发脉冲、脉冲间距δTE、第2翻转角为α的激发脉冲、第2移位梯度、第1翻转角为β的重聚脉冲、第一个采样回波链、第2翻转角为β的重聚脉冲、第二个采样回波链。所述脉冲间距δTE中包含了第1 移位梯度和扩散梯度。第1移位梯度包含了沿频率维(x方向)的梯度Gro1和沿相位维(y方向) 的梯度Gpe1。第2移位梯度包含了沿频率维(x方向)的梯度Gro2和沿相位维(y方向)的梯度Gpe2。
所述第1翻转角和第2翻转角为α的激发脉冲结合第1移位梯度和第2移位梯度，使两个回波在k空间的中心产生偏移，重聚脉冲以及两个小角度激发脉冲都与层选方向(z方向) 的层选梯度Gss相结合进行层选；第2翻转角为α的激发脉冲前后分别施加回波延时，其中δTE的长度为1/3回波链长，重聚脉冲前后有x，y，z三个方向的破坏梯度作用。
采样回波链是由分别作用在x，y方向的梯度链组成；x方向的梯度链由一系列正负梯度构成，且每个梯度的面积是所述移位梯度Gro1的三倍；y方向的梯度链是由一系列大小相等的尖峰梯度构成，且尖峰梯度的面积和等于移位梯度Gpe1面积的三倍。
在第一个采样回波链之前，频率和相位方向分别施加了重聚梯度，频率维的重聚梯度的面积是频率维采样梯度Gro梯度面积的一半，方向与Gro相反；相位维的重聚梯度的面积是所述所有尖峰梯度面积的一半，方向与尖峰梯度相反；第一个采样回波链的信号包含第一个回波和第二个回波的信号。
在第一个采样回波链之后，利用重聚脉冲产生两个新的自旋回波，并利用第二个采样回波链采样，第二个采样回波链的信号包含了第三个回波和第四个回波的信号。
(3)执行步骤(2)设置好的所述DT2M-OLED成像序列，进行数据采样；数据采样完成后，得到两个回波链的信号。
(4)对步骤(3)得到的第一个采样回波链和第二个采样回波链的信号进行归一化、充零和快速傅里叶变换将k空间的信号转换到图像域。对得到的信号进行分析，并对回波信号磁化矢量演化进行理论推导，在第一个采样回波链采样期，信号表达式如下：
式中
Figure BDA0001728884740000071
是扩散因子。
在第二个采样回波链采样期，信号表达式如下：
通过实验发现，当α＝45°，β＝180°时，四个回波信号的强度都相对较高；
(5)步骤(5)～(7)主要是把步骤(4)中图像域的信号用基于卷积神经网络的方法重建。根据实验样品的特征生成随机的模板，对模板进行模拟采样得到k空间信号，然后对信号进行归一化、充零和快速傅里叶变换得到训练数据。
所述随机模板是根据实验样本的特征分布使用计算机批量随机生成，同时保证模板的复杂性较高能够包含实验样本的所有特征；在模拟采样过程中，考虑到真实实验环境变化，加入了一些不稳定因素，比如激发脉冲角度偏差、移位梯度偏差和噪声等来提高网络模型对不理想实验环境的鲁棒性。
(6)使用tensorflow深度学习框架和python搭建卷积神经网络模型，设置好训练的相关参数。
在步骤(6)中，所述卷积神经网络模型主要包括：网络的主体结构以及相关的训练参数。网络模型的目标函数为：
Figure BDA0001728884740000073
其中，N是同一批训练图像的数量,f(·)是训练网络,W和b是网络参数，xk是第k个回波链信号经步骤(4)后得到的图像，y是输入图像对应的模板，ychange是对y中小于0.06的值置为0.06之后的矩阵，ymask是对y使用canny算子求得的图像边缘信息，λ表示正则项参数，▽代表梯度算子。
(7)使用步骤(6)得到的数据训练网络，直至网络收敛并达到稳定得到训练好的网络模型，然后对步骤(4)得到的实验数据进行重建，得到可靠的ADC和T2图像。
所述训练好的网络模型由于使用随机模板进行训练，泛化性较强，能适用于多种样品的重建。
以下给出具体实施例：
用基于重叠回波的单扫描同步磁共振扩散及T2成像方法进行了模型实验，用来验证本发明的可行性。实验前，用水、四水、氯化锰(Mncl2.4H2O)和琼脂糖按照表1的配比配制如图1中的凝胶模型。实验是在小动物核磁共振7T成像仪下进行的。在磁共振成像仪操作台上，打开成像仪中相应的操作软件，首先对成像物体进行感兴趣区域定位，然后进行调谐、匀场、功率和频率校正。为了评价该方法得到图像的有效性，在相同环境下进行平面回波(Echo Planar Imaging，EPI)成像实验作为对比。然后导入已经编译好的DT2M-OLED成像序列(如图1)，根据具体实验情况，设置脉冲序列的各个参数，本实施例的实验参数设置如下：成像视野FOV为55cm×55cm，45°激发脉冲的激发时间为1ms，180°重聚脉冲的激发时间为1.6ms，扩散系数b＝500s/mm2，TE1＝48.61ms，TE2＝25.10ms，TE3＝25.10ms，TE4＝48.61ms，x方向采样点数Nx为96，y方向采样点数Ny为96。将以上实验参数设置好后，直接开始采样。
表1
1 2 3 4 5
总质量(g) 5 5 5 5 50
Mncl<sub>2</sub>.4H<sub>2</sub>O(μg) 12 9 6 3 0
甘油(g) 0 0.35 0.7 1 0
琼脂糖(mg) 50 50 50 50 500
数据采样完成后，按照上述步骤(4)～(7)对数据进行重建，重建出来ADC和T2的图像如图2所示。(a)是用EPI序列采模型的结构图，其中的5个用红色圆圈标出的区域是感兴趣区域(region ofinterest，ROI)，分别用数字1～5编号。(b)是用DT2M-OLED序列采集第一个采样回波链的原始幅值图，图中斜条纹是两个回波信号重叠在一起造成的。(c)是用EPI拟合出来的ADC图像。(d)是DT2M-OLED数据重建得到的ADC图像。(e)是用EPI拟合出来的T2图像。(f)是DT2M-OLED数据重建得到的T2图像。图3给出展示5个ROI的T2统计结果，图4给出展示5个ROI的ADC统计结果，图3和4是对图2中的5个感兴趣区域的ADC (图4)和T2(图3)进行定量统计，即对每个感兴趣区域的ADC和T2值求平均。从图2 中可以看出DT2M-OLED序列得到的ADC和T2图像整体上与EPI一致，图3和4中的统计结果也表明5个ROI区域的ADC和T2值都有较好的准确性。
本发明用两个相同偏转角的小角度激发脉冲和两个回波链采样产生四个回波，每个回波链采集两个回波信号。在第一个激发脉冲之后有一段演化时间和一对扩散梯度，使得第一次采样中的两个回波信号的横向弛豫时间和扩散加权不同。同时，在每个激发脉冲之后都加一个频率维和相位维的移位梯度使得不同的激发脉冲产生的信号在k空间的位置不一样。在第一次采样结束后，用一个重聚脉冲重聚两个回波信号，然后进行第二次采样，在第二次采样中可以获得了两个具有相同横向弛豫时间和不同扩散加权的回波信号。然后将采样信号经过归一化、充零和快速傅里叶变换之后输入到已经训练好的卷积神经网络中重建得到定量T2图像和ADC图像。卷积神经网络的训练数据来源于模拟数据，通过随机生成模板，然后模拟实验环境采样得到网络的输入图像，模板作为标签，通过训练得到输入图像和输出图像之间的映射关系。
符号说明
符号 符号说明
T<sub>1</sub> 纵向弛豫时间
T<sub>2</sub> 横向弛豫时间
T<sub>2</sub>* 受磁场不均匀影响后的横向弛豫时间
b 扩散系数
ADC 表观扩散系数
α 激发脉冲翻转角度
β 重聚脉冲翻转角
G<sub>d</sub> 矢量，扩散梯度
G<sub>ro1</sub> 矢量，频率维第一个移位梯度
G<sub>ro2</sub> 矢量，频率维第二个移位梯度
G<sub>pe1</sub> 矢量，相位维第一个移位梯度，持续时间与G<sub>ro1</sub>相同
G<sub>pe2</sub> 矢量，相位维第二个移位梯度，持续时间与G<sub>ro2</sub>相同
G<sub>1</sub> 矢量，频率维第一个移位梯度和相位维第一个移位梯度的矢量和
G<sub>2</sub> 矢量，频率维第二个移位梯度和相位维第二个移位梯度的矢量和
δ<sub>TE</sub> 两个激发脉冲之间的时间间隔
δ 两个扩散梯度的持续时间间隔
Δ 两个扩散梯度开始的时间间隔
G<sub>ss</sub> 层选脉冲
G<sub>ro</sub> 矢量，频率维采样梯度
G<sub>pe</sub> 矢量，相位维采样梯度，即所有尖峰梯度的和
N<sub>pe</sub> 相位维采样点数
TE<sub>1</sub> 第一个回波第一次重聚的回波时间
TE<sub>2</sub> 第二个回波第一次重聚的回波时间
TE<sub>3</sub> 第一个回波第二次重聚的回波时间
TE<sub>4</sub> 第二个回波第二次重聚的回波时间
Ecoh1 第一个回波信号第一次重聚的中心位置
Echo2 第二个回波信号第一次重聚的中心位置
Ecoh3 第一个回波信号第二次重聚的中心位置
Echo4 第二个回波信号第二次重聚的中心位置
δ<sub>1</sub> G<sub>ro1</sub>的持续时间
δ<sub>2</sub> G<sub>ro2</sub>的持续时间
γ 磁旋比
r 矢量，产生信号的点的空间位移
S<sub>1</sub> 第一个回波第一次重聚产生的信号
S<sub>2</sub> 第二个回波第一次重聚产生的信号
S<sub>3</sub> 第一个回波第二次重聚产生的信号
S<sub>4</sub> 第二个回波第二次重聚产生的信号
本发明提出的方法可以在单次扫描中同时获得可靠的T2图像和ADC图像。

Claims (4)

1.基于重叠回波的单扫描同步磁共振扩散及T2成像方法，其特征在于包括以下步骤：
1)在磁共振成像仪操作台上，打开磁共振成像仪中相应的操作软件，首先对成像物体进行感兴趣区域定位，然后进行调谐、匀场、功率和频率校正；
2)导入事先编译好的DT2M-OLED成像序列，根据具体的实验情况，设置脉冲序列的各个参数；
3)执行步骤2)设置好的DT2M-OLED成像序列，进行数据采样，数据采样完成后，得到两个回波链的信号；
所述DT2M-OLED成像序列的结构依次为：第1翻转角为α的激发脉冲、脉冲间距δTE、第2翻转角为α的激发脉冲、第2移位梯度、第1翻转角为β的重聚脉冲、第一个采样回波链、第2翻转角为β的重聚脉冲、第二个采样回波链；所述脉冲间距δTE中包含第1移位梯度和扩散梯度，第1移位梯度包含沿频率维x方向的梯度Gro1和沿相位维y方向的梯度Gpe1；第2移位梯度包含沿频率维x方向的梯度Gro2和沿相位维y方向的梯度Gpe2；
所述第1翻转角和第2翻转角为α的激发脉冲结合第1移位梯度和第2移位梯度，使两个回波在k空间的中心产生偏移，重聚脉冲以及两个小角度激发脉冲都与层选方向z方向的层选梯度Gss相结合进行层选；第2翻转角为α的激发脉冲前后分别施加回波延时，其中δTE的长度为1/3回波链长，重聚脉冲前后有x，y，z三个方向的破坏梯度作用；
所述第一个采样回波链和第二个采样回波链是由分别作用在x，y方向的梯度链组成；x方向的梯度链由一系列正负梯度构成，且每个梯度的面积是所述移位梯度Gro1的三倍；y方向的梯度链是由一系列大小相等的尖峰梯度构成，且尖峰梯度的面积和等于移位梯度Gpe1面积的三倍；
在第一个采样回波链之前，频率和相位方向分别施加了重聚梯度，频率维的重聚梯度的面积是频率维采样梯度Gro梯度面积的一半，方向与Gro相反；相位维的重聚梯度的面积是所述所有尖峰梯度面积的一半，方向与尖峰梯度相反；第一个采样回波链的信号包含第一个回波和第二个回波的信号；
在第一个采样回波链之后，利用重聚脉冲产生两个新的自旋回波，并利用第二个采样回波链采样，第二个采样回波链的信号包含第三个回波和第四个回波的信号；
4)对步骤3)得到的两个回波链的信号进行归一化、充零和快速傅里叶变换，将k空间的信号转换到图像域，再将图像域的信号用基于卷积神经网络的方法重建，得实验数据；
所述DT2M-OLED成像序列采集到的两个回波链的信号进行分析，并对回波信号磁化矢量演化进行理论推导，在第一个采样回波链采样期，信号表达式如下：
Figure FDA0002255715290000021
式中
Figure FDA0002255715290000022
是扩散系数；S1为第一个回波第一次重聚产生的信号，S2为第二个回波第一次重聚产生的信号，α为激发脉冲的角度，β为重聚脉冲角度，T2为横向弛豫时间，TE1为第一个回波第一次重聚的回波时间，TE2为第二个回波第一次重聚的回波时间；
在第二个采样回波链采样期，信号表达式如下：
Figure FDA0002255715290000023
其中，S3为第一个回波第二次重聚产生的信号，S4为第二个回波第二次重聚产生的信号，α为激发脉冲的角度，β为重聚脉冲角度，b为扩散系数，T2为横向弛豫时间，TE1为第一个回波第一次重聚的回波时间，TE2为第二个回波第一次重聚的回波时间，TE3为第一个回波第二次重聚的回波时间，TE4为第二个回波第二次重聚的回波时间；通过实验发现，当α＝45°，β＝180°时，四个回波信号的强度都相对较高；
5)根据实验样品的特征生成随机的模板，对模板进行模拟采样得到k空间信号，然后对k空间信号进行归一化、充零和快速傅里叶变换，得到训练数据；
6)使用tensorflow深度学习框架和python搭建卷积神经网络模型，设置好训练的相关参数，得数据训练网络；
7)使用步骤6)得到的数据训练网络，直至数据训练网络收敛并达到稳定，得到训练好的网络模型，然后对步骤4)得到的实验数据进行重建，得到ADC和T2图像。
2.如权利要求1所述基于重叠回波的单扫描同步磁共振扩散及T2成像方法，其特征在于在步骤5)中，所述随机的模板是根据实验样本的特征分布使用计算机批量随机生成，同时保证模板的复杂性较高能够包含实验样本的所有特征；在模拟采样过程中，加入不稳定因素，提高网络模型对不理想实验环境的鲁棒性；所述不稳定因素包括激发脉冲角度偏差、移位梯度偏差和噪声。
3.如权利要求1所述基于重叠回波的单扫描同步磁共振扩散及T2成像方法，其特征在于在步骤6)中，所述使用tensorflow深度学习框架和python搭建卷积神经网络模型包括：网络的主体结构以及相关的训练参数，网络模型的目标函数为：
Figure FDA0002255715290000024
其中，N是同一批训练图像的数量，f(·)是训练网络，W和b是网络参数，xk是第k个回波链信号经步骤4)后得到的图像，yk是输入图像对应的模板，ychange是对yk中小于0.06的值置为0.06之后的矩阵，ymask是对yk使用canny算子求得的图像边缘信息，λ表示正则项参数，▽代表梯度算子。
4.如权利要求1所述基于重叠回波的单扫描同步磁共振扩散及T2成像方法，其特征在于在步骤7)中，所述使用步骤6)得到的数据训练网络由于使用随机的模板进行训练，适用于多种样品的重建。
CN201810765276.7A 2018-07-12 2018-07-12 基于重叠回波的单扫描同步磁共振扩散及t2成像方法 Active CN109100669B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810765276.7A CN109100669B (zh) 2018-07-12 2018-07-12 基于重叠回波的单扫描同步磁共振扩散及t2成像方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810765276.7A CN109100669B (zh) 2018-07-12 2018-07-12 基于重叠回波的单扫描同步磁共振扩散及t2成像方法

Publications (2)

Publication Number Publication Date
CN109100669A CN109100669A (zh) 2018-12-28
CN109100669B true CN109100669B (zh) 2020-01-14

Family

ID=64846217

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810765276.7A Active CN109100669B (zh) 2018-07-12 2018-07-12 基于重叠回波的单扫描同步磁共振扩散及t2成像方法

Country Status (1)

Country Link
CN (1) CN109100669B (zh)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109814056B (zh) * 2019-02-01 2020-11-10 沈阳工业大学 一种获取磁共振精准定量图像的方法
CN110174632B (zh) * 2019-06-10 2021-06-01 上海东软医疗科技有限公司 磁共振成像方法、装置、成像设备及磁共振成像系统
US11009578B2 (en) * 2019-07-17 2021-05-18 GE Precision Healthcare LLC Systems and methods for predicting B1+ maps from magnetic resonance calibration images
CN110346743B (zh) * 2019-07-22 2021-09-14 上海东软医疗科技有限公司 一种磁共振弥散加权成像方法和装置
CN110782452B (zh) * 2019-11-05 2022-08-12 厦门大学 一种t2定量图像成像方法及系统
CN112304811B (zh) * 2020-10-22 2023-02-03 无锡鸣石峻致医疗科技有限公司 一种对运动不敏感的低场核磁共振表观扩散系数测量方法
CN112924913B (zh) * 2021-02-02 2022-08-12 厦门大学 一种时空编码磁共振成像超分辨率重建方法及系统
CN113805133B (zh) * 2021-09-13 2022-05-20 上海交通大学 一种减少图像失真的磁共振平面回波成像方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016207035A1 (en) * 2015-06-26 2016-12-29 Koninklijke Philips N.V. Phase corrected dixon magnetic resonance imaging
CN105232045B (zh) * 2015-11-26 2017-12-29 厦门大学 基于双回波的单扫描定量磁共振扩散成像方法
CN105232046B (zh) * 2015-11-26 2017-12-29 厦门大学 一种基于重叠回波的单扫描定量磁共振t2成像方法
CN107045115A (zh) * 2017-05-04 2017-08-15 厦门大学 基于双回波的单扫描定量磁共振t2*成像方法
CN108010100B (zh) * 2017-12-07 2020-02-18 厦门大学 一种基于残差网络的单扫描磁共振定量t2成像重建方法

Also Published As

Publication number Publication date
CN109100669A (zh) 2018-12-28

Similar Documents

Publication Publication Date Title
CN109100669B (zh) 基于重叠回波的单扫描同步磁共振扩散及t2成像方法
CN108010100B (zh) 一种基于残差网络的单扫描磁共振定量t2成像重建方法
Renvall et al. Automatic cortical surface reconstruction of high-resolution T1 echo planar imaging data
Shi et al. Parallel imaging and compressed sensing combined framework for accelerating high‐resolution diffusion tensor imaging using inter‐image correlation
US10552955B2 (en) Imaging acceleration methods for MRI parameter mapping
CN108132274B (zh) 不均匀磁场下回波平面成像无参考扫描图像畸变矫正方法
CN108663644B (zh) 基于残差网络重建的单扫描定量磁共振t2*成像方法
US10338178B2 (en) System and method for high-resolution spectroscopic imaging
US20100329528A1 (en) Imaging of moving objects
CN110889897B (zh) 一种体素内不相干运动磁共振成像参数重建方法及系统
Liu et al. Improved low-rank filtering of magnetic resonance spectroscopic imaging data corrupted by noise and $ B_0 $ field inhomogeneity
US20240095889A1 (en) Systems and methods for magnetic resonance image reconstruction with denoising
US11009577B2 (en) System and method for Nyquist ghost correction in medical imaging
Fang et al. Deepsti: towards tensor reconstruction using fewer orientations in susceptibility tensor imaging
Gan et al. SS-JIRCS: Self-supervised joint image reconstruction and coil sensitivity calibration in parallel MRI without ground truth
CN115137347B (zh) 三维超短回波时间磁共振指纹成像的髓鞘定量成像方法
Dai et al. Distortion-free diffusion imaging using self-navigated Cartesian Echo-planar time resolved acquisition and joint magnitude and phase constrained reconstruction
Wang et al. Simultaneous pure T2 and varying T2’-weighted BOLD fMRI using Echo Planar Time-resolved Imaging (EPTI) for mapping laminar fMRI responses
WO2022212244A1 (en) Distortion-free diffusion and quantitative magnetic resonance imaging with blip up-down acquisition of spin- and gradient-echoes
Cheng et al. deep MR parametric mapping with unsupervised multi-tasking framework
CN113808176B (zh) 基于装置成像坐标的mr图像配准方法、分析方法及其装置
US20240183924A1 (en) Distortion-free diffusion and quantitative magnetic resonance imaging with blip up-down acquisition of spin- and gradient-echoes
Baqee A Novel Approach of Obtaining Optimal Solution for Iterative Self-Consistent Parallel Imaging Reconstruction
Odille Motion-Corrected Reconstruction
CN117331010A (zh) 基于多重重叠回波的多参数同时定量成像方法和系统

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