CN109143134A - 一种分段读出扩散加权成像方法、装置及可存储介质 - Google Patents

一种分段读出扩散加权成像方法、装置及可存储介质 Download PDF

Info

Publication number
CN109143134A
CN109143134A CN201710442188.9A CN201710442188A CN109143134A CN 109143134 A CN109143134 A CN 109143134A CN 201710442188 A CN201710442188 A CN 201710442188A CN 109143134 A CN109143134 A CN 109143134A
Authority
CN
China
Prior art keywords
pulse
diffusion
weighted
imaging
gradient
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.)
Granted
Application number
CN201710442188.9A
Other languages
English (en)
Other versions
CN109143134B (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.)
Siemens Shenzhen Magnetic Resonance Ltd
Original Assignee
Siemens Shenzhen Magnetic Resonance Ltd
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 Siemens Shenzhen Magnetic Resonance Ltd filed Critical Siemens Shenzhen Magnetic Resonance Ltd
Priority to CN201710442188.9A priority Critical patent/CN109143134B/zh
Priority to EP18177302.9A priority patent/EP3415941A1/en
Priority to US16/007,332 priority patent/US10928476B2/en
Publication of CN109143134A publication Critical patent/CN109143134A/zh
Application granted granted Critical
Publication of CN109143134B publication Critical patent/CN109143134B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56545—Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by finite or discrete sampling, e.g. Gibbs ringing, truncation artefacts, phase aliasing artefacts
    • 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/46—NMR spectroscopy
    • G01R33/4616—NMR spectroscopy using specific RF pulses or specific modulation schemes, e.g. stochastic excitation, adiabatic RF pulses, composite pulses, binomial pulses, Shinnar-le-Roux pulses, spectrally selective pulses not being used for spatial selection
    • 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
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/543—Control of the operation of the MR system, e.g. setting of acquisition parameters prior to or during MR data acquisition, dynamic shimming, use of one or more scout images for scan plane prescription
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/5602—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution by filtering or weighting based on different relaxation times within the sample, e.g. T1 weighting using an inversion pulse
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/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/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/56554—Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by acquiring plural, differently encoded echo signals after one RF excitation, e.g. correction for readout gradients of alternating polarity in 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/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/5617—Echo train techniques involving acquiring plural, differently encoded, echo signals after one RF excitation, e.g. using gradient refocusing in echo planar imaging [EPI], RF refocusing in rapid acquisition with relaxation enhancement [RARE] or using both RF and gradient refocusing in gradient and spin echo imaging [GRASE] using RF refocusing, e.g. RARE

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • General Physics & Mathematics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • General Health & Medical Sciences (AREA)
  • Radiology & Medical Imaging (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Vascular Medicine (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Optics & Photonics (AREA)
  • Medical Informatics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Biophysics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Magnetic Resonance Imaging Apparatus (AREA)

Abstract

本发明实施方式公开了一种分段读出扩散加权成像方法、装置及可存储介质。包括：针对磁化强度矢量施加非线性相位射频激发脉冲，并在选层方向上施加持续时间与所述非线性相位射频激发脉冲相对应的选层梯度脉冲，以将磁化强度矢量翻转到X‑Y平面；针对翻转到X‑Y平面上的磁化强度矢量执行扩散加权处理；利用读出方向分段采样序列从扩散加权处理后的磁化强度矢量中分段读出成像数据，其中在选层方向上施加视角倾斜梯度脉冲。本发明实施方式可以减少图像畸变，还可以降低图像模糊，从而提高成像质量。

Description

一种分段读出扩散加权成像方法、装置及可存储介质
技术领域
本发明涉及扩散成像技术领域，特别是涉及一种分段读出扩散加权成像方法、装置及可存储介质。
背景技术
扩散加权成像与传统的磁共振(Magnetic Resonance，MR)技术不同，它主要依赖于水分子的运动而非组织的自旋质子密度、T1值或T2值，为组织成像对比提供了一种崭新的技术。组织中的水分子扩散是一种随机的热运动，扩散的方向与幅度受局部组织微环境的影响，例如在脑梗、肿瘤等病灶区水分子扩散受到限制，使得扩散度降低。MR通过H+的磁化来标记分子而不干扰它的扩散过程。在任一常规MR成像序列中加入扩散梯度突出扩散效应，即可执行扩散加权成像，可以对组织中水分子的扩散行为直接进行检测。
常规的扩散成像序列通常使用单次激发平面回波(single-shot EPI)序列进行图像采集，这种序列在磁场不均匀区域(例如空气组织交界面和骨骼肌肉交界面)容易产生相位误差，导致图像发生几何畸变。目前已经出现了读出方向分段采集的平面回波(ReadoutSegmentation Of Long Variable Echo-trains，RESOLVE)序列。与单次激发平面回波相比，RESOLVE序列可以降低磁敏感性伪影(Susceptibility Artifact)，提高成像质量。
然而，对于具有复杂几何图形和严重磁敏感性伪影的区域(比如，颈部)，RESOLVE序列所采集到的图像失真仍然可能较为严重，导致成像质量不理想。
发明内容
本发明实施方式提出一种分段读出扩散加权成像方法、装置及可存储介质，从而提高成像质量。
本发明实施方式的技术方案如下：
一种分段读出扩散加权成像方法，包括：
针对磁化强度矢量施加非线性相位射频激发脉冲，并在选层方向上施加持续时间与所述非线性相位射频激发脉冲相对应的选层梯度脉冲，以将磁化强度矢量翻转到X-Y平面；
针对翻转到X-Y平面上的磁化强度矢量执行扩散加权处理；
利用读出方向分段采样序列从扩散加权处理后的磁化强度矢量中分段读出成像数据，其中在选层方向上施加视角倾斜梯度脉冲。
在一个实施方式中，所述非线性相位射频激发脉冲在时间域上为方波或近似方波。
在一个实施方式中，所述非线性相位射频激发脉冲的非线性相位大小为K；
K≤3.6/(BW2*FTW)，其中BW为非线性相位射频激发脉冲的带宽，FTW为非线性相位射频激发脉冲的部分过渡带宽，FTW＝(ws-wp)/BW；ws是非线性相位射频激发脉冲的止带频率，wp是非线性相位射频激发脉冲的通带频率。
在一个实施方式中，所述针对翻转到X-Y平面上的磁化强度矢量执行扩散加权处理包括：针对翻转到X-Y平面上的磁化强度矢量施加扩散加权梯度脉冲，在射频方向上施加回聚脉冲，在选层方向上施加选层梯度脉冲和损毁梯度脉冲。
在一个实施方式中，所述利用读出方向分段采样序列从扩散加权处理后的磁化强度矢量中分段读出成像数据包括：在读出方向上施加成像回波读出梯度脉冲以读出成像数据，在读出方向上施加导航回波梯度读出脉冲以读出导航数据，用以校正分段采集的成像回波之间的相位误差，在相位编码方向上施加成像回波相位编码梯度脉冲和导航回波相位编码梯度脉冲。
在一个实施方式中，在选层方向上施加视角倾斜梯度脉冲包括：在选层方向上导航回波的回聚脉冲的两边施加视角倾斜梯度脉冲，所述视角倾斜梯度脉冲包括成像回波的视角倾斜梯度脉冲和导航回波的视角倾斜梯度脉冲；该方法还包括：
在选层方向上施加成像回波的视角倾斜梯度脉冲的预散相梯度脉冲、成像回波的视角倾斜梯度脉冲的聚相梯度脉冲和导航回波的视角倾斜梯度脉冲的预散相梯度脉冲。
一种分段读出扩散加权成像装置，包括：
翻转模块，用于针对磁化强度矢量施加非线性相位射频激发脉冲，并在选层方向上施加持续时间与所述非线性相位射频激发脉冲相对应的选层梯度脉冲，以将磁化强度矢量翻转到X-Y平面；
扩散加权模块，用于针对翻转到X-Y平面上的磁化强度矢量执行扩散加权处理；
分段读出模块，用于利用读出方向分段采样序列从扩散加权处理后的磁化强度矢量中分段读出成像数据，其中在选层方向上施加视角倾斜梯度脉冲。
在一个实施方式中，所述非线性相位射频激发脉冲在时间域上为方波或近似方波。
在一个实施方式中，所述非线性相位射频激发脉冲的非线性相位大小为K；
K≤3.6/(BW2*FTW)，其中BW为非线性相位射频激发脉冲的带宽，FTW为非线性相位射频激发脉冲的部分过渡带宽，FTW＝(ws-wp)/BW；ws是非线性相位射频激发脉冲的止带频率，wp是非线性相位射频激发脉冲的通带频率。
在一个实施方式中：
扩散加权模块，用于针对翻转到X-Y平面上的磁化强度矢量施加扩散加权梯度脉冲，在射频方向上施加回聚脉冲，在选层方向上施加选层梯度脉冲和损毁梯度脉冲；和/或
分段读出模块，用于在读出方向上施加成像回波读出梯度脉冲以读出成像数据，在读出方向上施加导航回波读出梯度脉冲以读出导航数据，用以校正分段采集的成像回波之间的相位误差，在相位编码方向上施加成像回波相位编码梯度脉冲和导航回波相位编码梯度脉冲。
在一个实施方式中，分段读出模块，用于在选层方向上导航回波的回聚脉冲的两边施加视角倾斜梯度脉冲，所述视角倾斜梯度脉冲包括成像回波的视角倾斜梯度脉冲和导航回波的视角倾斜梯度脉冲；在选层方向上施加成像回波的视角倾斜梯度脉冲的预散相梯度脉冲、成像回波的视角倾斜梯度脉冲的聚相梯度脉冲和导航回波的视角倾斜梯度脉冲的预散相梯度脉冲。
一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现如上任一项所述的分段读出扩散加权成像方法的步骤。
从上述技术方案可以看出，在本发明实施方式中，针对磁化强度矢量施加非线性相位射频激发脉冲，并在选层方向上施加持续时间与非线性相位射频激发脉冲相对应的选层梯度脉冲，以将磁化强度矢量翻转到X-Y平面；针对翻转到X-Y平面上的磁化强度矢量执行扩散加权处理；利用读出方向分段采样序列从扩散加权处理后的磁化强度矢量中分段读出成像数据，其中在选层方向上施加视角倾斜梯度脉冲。由此可见，本发明实施方式将视角倾斜梯度脉冲应用到读出方向分段采样序列中，可以补偿由于场不均匀性引起的相位误差，改善图像畸变。另外，本发明实施方式还利用非线性相位射频激发脉冲克服视角倾斜梯度脉冲导致的图像模糊。因此，本发明实施方式可以提高成像质量。
附图说明
图1为根据本发明实施方式分段读出扩散加权成像方法的流程图。
图2为根据本发明实施方式非线性相位激发脉冲的示范性幅度图。
图3为根据本发明实施方式非线性相位激发脉冲的示范性相位图。
图4为根据本发明实施方式的采用VAT和非线性相位激发脉冲的RESOLVE序列图。
图5A为根据快速自旋回波(TSE)序列获取的采样图像示范图。
图5B为根据现有技术的RESOLVE序列获取的采样图像示范图。
图5C为根据本发明实施方式的RESOLVE序列获取的采样图像示范图。
图6为根据本发明实施方式的分段读出扩散加权成像装置结构图。
具体实施方式
为了使本发明的技术方案及优点更加清楚明白，以下结合附图及实施方式，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施方式仅仅用以阐述性说明本发明，并不用于限定本发明的保护范围。
为了描述上的简洁和直观，下文通过描述若干代表性的实施方式来对本发明的方案进行阐述。实施方式中大量的细节仅用于帮助理解本发明的方案。但是很明显，本发明的技术方案实现时可以不局限于这些细节。为了避免不必要地模糊了本发明的方案，一些实施方式没有进行细致地描述，而是仅给出了框架。下文中，“包括”是指“包括但不限于”，“根据……”是指“至少根据……，但不限于仅根据……”。由于汉语的语言习惯，下文中没有特别指出一个成分的数量时，意味着该成分可以是一个也可以是多个，或可理解为至少一个。
申请人经过大量试验发现：现有RESOLVE序列的图像失真问题主要归因于场不均匀性引起的相位误差。申请人还发现，将视角倾斜(View Angle Tilting，VAT)技术应用到RESOLVE序列中以补偿因场不均匀性所导致的相位误差，可以减少图像畸变。
另外，本领域技术人员可以意识到目前的扩散加权成像序列通常采用线性相位射频激发脉冲翻转磁化强度矢量。申请人发现：当采用线性相位射频激发脉冲时，VAT技术与RESOLVE序列的结合可能导致不期望的图像模糊。申请人经过大量试验还发现：利用非线性相位射频激发脉冲代替线性相位射频激发脉冲，可以显著降低甚至克服VAT技术导致的图像模糊。
基于上述分析，图1为根据本发明实施方式分段读出扩散加权成像方法的流程图。
如图1所示，该方法包括：
步骤101：针对磁化强度矢量施加非线性相位射频激发脉冲，并在选层方向上施加持续时间与非线性相位射频激发脉冲相对应的选层梯度脉冲，以将磁化强度矢量翻转到X-Y平面。
在磁共振扩散加权成像中，通常将主磁场B0的方向定义为三维笛卡尔坐标系的Z向，将与主磁场B0正交的其他两个方向分别定义为X向和Y向。X向和Y向共同构成X-Y平面。
在本发明实施方式中，利用非线性相位射频激发脉冲代替线性相位射频激发脉冲。
非线性相位射频激发脉冲是一种具有非线性相位的射频脉冲。传统的线性相位射频脉冲可以看作将选频范围内的磁化矢量在同一时刻翻转到X-Y平面上，因此其能量集中在一个持续时间短的主峰中，使得脉冲幅度分布不均。通过改变射频脉冲的相位，可以将选频范围内的磁化矢量在脉冲持续时间内依次翻转到X-Y平面上，因此其能量不用集中在持续时间短的主峰中，这种相位改变后的射频脉冲即为具有非线性相位的射频脉冲。利用非线性相位的脉冲设计方法，可以在保证脉冲选层轮廓的同时得到幅度分布较为平均的激励脉冲。
优选地，本发明实施方式的非线性相位射频激发脉冲在时间域上呈现方波或近似方波形状。在时间域上呈现方波或近似方波形状的非线性相位射频激发脉冲，幅度分布均匀，可以减少K空间的滤波作用，并降低图像模糊。
非线性相位大小K是非线性相位射频激发脉冲的重要参数。当K增大后，射频脉冲的幅度在其持续时间内的分布更加平均，但是会影响脉冲的选层轮廓。
在本发明实施方式中，K优选满足下列标准：K≤3.6/(BW2*FTW)，其中BW是射频脉冲的带宽，FTW是部分过度带宽，FTW＝(ws-wp)/BW，其中Ws和WP分别是非线性相位射频激发脉冲的止带频率和通带频率。
K的增大会导致射频脉冲长度的增加，因此在实际应用中适当截断低于一定幅度的射频脉冲从而得到合适的射频脉冲长度。但是，过度的截断也会影响射频脉冲的选层轮廓。
在本发明实施方式中，优选截断幅度低于5％最大幅度的射频脉冲部分，即将幅度值小于最大幅度值5％的射频脉冲幅度置零，从而减小脉冲持续时间。
具体地，本发明实施方式可以通过多种方式设计非线性相位射频激发脉冲。
比如，可以采用希纳-勒鲁(Shinnar-Le Roux，SLR)方法设计本发明实施方式的非线性相位射频激发脉冲。在SLR方法中，射频脉冲可以看作是两个复数多项式A和B的组合，通过重新设计复数多项式B，可以改变射频脉冲的形状及其选层轮廓。当复数多项式B确定后，可以通过最小化脉冲能量计算得到复数多项式A。复数多项式B多被设计成线性相位，最大相位和最小相位，此时脉冲幅度分布不均，应用在视角倾斜技术上会引发严重的滤波作用。通过改变复数多项式B的相位，可以在保证脉冲选层轮廓的同时改变脉冲的幅度分布。在本发明实施方式中，期望得到幅度分布较为均匀的激励脉冲，可以利用SLR方法设计得到复数多项式B，然后通过给其施加不同的非线性相位，优化得到幅度分布平均的激励脉冲。
图2为根据本发明实施方式非线性相位激发脉冲的示范性幅度图；图3为根据本发明实施方式非线性相位激发脉冲的示范性相位图。在图2中，横轴为激发脉冲的取样点，纵轴为激发脉冲的归一化幅度；在图3中，横轴为激发脉冲的取样点，纵轴为激发脉冲的归一化相位。
由图2和图3可见，相比具有单主峰的线性相位激发脉冲，非线性相位激发脉冲具有多个主峰，因此主峰宽度更长，幅度分布平均。
以上示范性描述了本发明实施方式所采用的非线性相位射频激发脉冲。本领域技术人员可以意识到，这种描述仅是示范性的，并不用于限定本发明实施方式的保护范围。
步骤102：针对翻转到X-Y平面上的磁化强度矢量执行扩散加权处理。
在这里，可以针对翻转到X-Y平面上的磁化强度矢量施加扩散加权梯度脉冲，在射频方向上施加回聚脉冲，在选层方向上施加选层梯度脉冲和损毁梯度脉冲，从而执行扩散加权处理。
步骤103：利用RESOLVE序列从扩散加权处理后的磁化强度矢量中分段读出成像数据，其中在选层方向上施加VAT梯度。
VAT技术最初是用来解决传统自旋回波序列中由于场不均匀导致的图像畸变以及读出方向上的化学位移问题。VAT技术还被应用到平面回波序列：即在施加相位编码梯度的同时，额外在选层方向上施加梯度。此时读出的MR信号在这两个梯度的同时作用下，会形成一个在相位编码和选层方向上的视角倾斜，从而减少重建出来的信号由于场不均匀导致的图像畸变。由于RESOLVE序列的回波间距远小于单次激发平面回波成像，因此更加适合引入视角倾斜技术以进一步改善场不均匀导致的图像畸变。
在相位编码方向上施加相位编码梯度时同时在选层方向上施加VAT梯度，磁共振信号S可以用下面的公式表示：
s(tm,tn)＝∫xx∫y∫ΔSl+Δz(x，y)ρ(x,y,z)exp(-jγnGvattbz)×exp(-jγΔB(x,y)nTesp)×exp[-jγ(mGxΔtxx+nGytby)]dxdydz；
其中：ΔB(x,y)是在(x，y)位置处的B0偏置；Gvat是在施加相位编码梯度的同时施加在选层方向上的梯度；ρ(x,y,z)是自旋密度；γ是旋磁比；Δtx是读出方向上的采样率；Tesp是回波间距；tb是相位编码梯度的时间；Δz(x,y)是由于场不均匀性引起的在选层方向上的位移，Δz(x,y)＝-ΔB(x,y)/Gz，其中Gz是选层梯度；ΔSl是层厚；m,n分别是读出方向和频率编码方向上的K空间坐标。
在一个实施方式中，在选层方向上施加VAT梯度脉冲包括：在选层方向上导航回波的回聚脉冲的两边施加VAT梯度脉冲，该VAT梯度脉冲包括成像回波的VAT梯度脉冲和导航回波的VAT梯度脉冲；该方法还包括：在选层方向上施加成像回波的VAT梯度脉冲的预散相梯度脉冲、成像回波的VAT梯度脉冲的聚相梯度脉冲和导航回波的VAT梯度脉冲的预散相梯度脉冲。
在一个实施方式中，利用RESOLVE序列从扩散加权处理后的磁化强度矢量中分段读出成像数据包括：在读出方向上施加成像回波读出梯度脉冲以读出成像数据，在读出方向上施加导航回波读出梯度脉冲以读出导航数据，用以校正分段采集的成像回波之间的相位误差，在相位编码方向上施加成像回波相位编码梯度脉冲和导航回波相位编码梯度脉冲。
下面示范性描述采用本发明实施方式的RESOLVE序列的执行过程。
图4为根据本发明实施方式的采用VAT和非线性相位激发脉冲的RESOLVE序列图。
在图4中，从上到下分别示出射频(R)脉冲、扩散梯度(GD)、选层梯度(GS)、读出方向(GR)和相位编码(GP)方向。
首先，施加90度的非线性相位射频激发脉冲30，而且在选层方向上施加持续时间与该非线性相位射频激发脉冲30相对应的选层梯度脉冲31，将磁化强度矢量翻转到X-Y平面上。
接着，针对翻转到X-Y平面上的磁化强度矢量执行扩散加权处理。具体包括：施加成像回波的回聚脉冲33(比如，180度的射频脉冲)，以及在扩散梯度方向上施加成像回波的回聚脉冲33两边的扩散加权梯度脉冲32，还在选层方向上施加持续时间与成像回波的回聚脉冲33相对应的选层梯度脉冲和损毁梯度脉冲34，其中损毁梯度脉冲用于损耗由180度射频脉冲33引起的自由感应延迟因子(FID)，从而针对翻转到X-Y平面上的磁化强度矢量执行扩散加权处理。
再接着，施加导航回波的回聚脉冲35(比如，180度的射频脉冲)，在读出方向上导航回波的回聚脉冲35的两边分别施加成像回波读出梯度脉冲37和导航回波读出梯度脉冲38，其中成像回波读出梯度脉冲37用于分段读出成像数据；导航回波读出梯度脉冲38用于读出导航数据，用以校正分段采集的成像回波之间的相位误差。在相位编码方向上，在导航回波的回聚脉冲35的两边分别施加成像回波相位编码梯度脉冲39和导航回波相位编码梯度脉冲40。在选层方向上，施加下列梯度脉冲：
(1)、位于导航回波的回聚脉冲35的两边的VAT梯度脉冲36。VAT梯度脉冲36包括：成像回波的VAT梯度脉冲(位于导航回波的回聚脉冲35的左边)和导航回波的VAT梯度脉冲(位于导航回波的回聚脉冲35的右边)。
(2)、成像回波的VAT梯度脉冲的预散相梯度脉冲41。
(3)、持续时间与导航回波的回聚脉冲35相对应的梯度组合体42。梯度组合体42包括：与导航回波的回聚脉冲35相对应的选层梯度和损毁梯度、成像回波的VAT梯度脉冲的聚相梯度和导航回波的VAT梯度脉冲的预散相梯度。
可见，VAT梯度脉冲36补偿由于场不均匀性导致的相位误差，并改善图像畸变。而且，非线性相位射频激发脉冲30由于幅度分布均匀，在时间域上近似矩形，可以改善VAT梯度脉冲36可能导致的图像模糊。
图5A为根据快速自旋回波(TSE)序列获取的采样图像示范图；图5B为根据现有技术的RESOLVE序列获取的采样图像示范图；图5C为根据本发明实施方式的RESOLVE序列获取的采样图像示范图。
由图5A-图5C可见，图5B的成像质量优于图5A的成像质量，而且图5C的成像质量明显优于图5B的成像质量。
基于上述描述，本发明实施方式还提出了一种分段读出扩散加权成像装置。
图6为根据本发明实施方式的分段读出扩散加权成像装置结构图。
如图6所示，该分段读出扩散加权成像装置60包括：
翻转模块61，用于针对磁化强度矢量施加非线性相位射频激发脉冲，并在选层方向上施加持续时间与非线性相位射频激发脉冲相对应的选层梯度脉冲，以将磁化强度矢量翻转到X-Y平面；
扩散加权模块62，用于针对翻转到X-Y平面上的磁化强度矢量执行扩散加权处理；
分段读出模块63，用于利用RESOLVE序列从扩散加权处理后的磁化强度矢量中分段读出成像数据，其中在选层方向上施加VAT梯度脉冲。
在一个实施方式中：非线性相位射频激发脉冲在时间域上为方波或近似方波。
在一个实施方式中：
所述非线性相位射频激发脉冲的非线性相位大小为K；
K≤3.6/(BW2*FTW)，其中BW为非线性相位射频激发脉冲的带宽，FTW为非线性相位射频激发脉冲的部分过渡带宽，FTW＝(ws-wp)/BW；ws是非线性相位射频激发脉冲的止带频率，wp是非线性相位射频激发脉冲的通带频率。
在一个实施方式中：
扩散加权模块62，用于针对翻转到X-Y平面上的磁化强度矢量施加扩散加权梯度脉冲，在射频方向上施加回聚脉冲，在选层方向上施加选层梯度脉冲和损毁梯度脉冲。
在一个实施方式中：
分段读出模块63，用于在读出方向上施加成像回波读出梯度脉冲以读出成像数据，在读出方向上施加导航回波读出梯度脉冲以读出导航数据，用以校正分段采集的成像回波之间的相位误差，在相位编码方向上施加成像回波相位编码梯度脉冲和导航回波相位编码梯度脉冲。
在一个实施方式中：
分段读出模块63，用于在选层方向上导航回波的回聚脉冲的两边施加VAT梯度脉冲，所述VAT梯度脉冲包括成像回波的VAT梯度脉冲和导航回波的VAT梯度脉冲；在选层方向上施加成像回波的VAT梯度脉冲的预散相梯度脉冲、成像回波的VAT梯度脉冲的聚相梯度脉冲和导航回波的VAT梯度脉冲的预散相梯度脉冲。
可以遵循一定规范的应用程序接口，将本发明实施方式所提出的分段读出扩散加权成像方法编写为安装到磁共振扩散成像控制主机、个人电脑、移动终端等中的插件程序，也可以将其封装为应用程序以供用户自行下载使用。
可以通过指令或指令集存储的储存方式将本发明实施方式所提出的分段读出扩散加权成像方法存储在各种存储介质上。这些存储介质包括但是不局限于：软盘、光盘、DVD、硬盘、闪存等。另外，还可以将本发明实施方式所提出的分段读出扩散加权成像方法应用到基于闪存(Nand flash)的存储介质中，比如U盘、CF卡、SD卡、SDHC卡、MMC卡、SM卡、记忆棒、xD卡等。上述计算机可读存储介质上存储有计算机程序，该计算机程序被处理器执行时实现如上所述的分段读出扩散加权成像方法的步骤。
综上所述，在本发明实施方式中，针对磁化强度矢量施加非线性相位射频激发脉冲，并在选层方向上施加持续时间与非线性相位射频激发脉冲相对应的选层梯度脉冲，以将磁化强度矢量翻转到X-Y平面；针对翻转到X-Y平面上的磁化强度矢量执行扩散加权处理；利用RESOLVE序列从扩散加权处理后的磁化强度矢量中分段读出成像数据，其中在选层方向上施加VAT梯度脉冲。本发明实施方式通过将VAT梯度脉冲应用到RESOLVE序列，可以补偿由于场不均匀性引起的相位误差，并改善图像畸变，还可以利用非线性相位射频激发脉冲克服VAT梯度脉冲导致的图像模糊，从而提高成像质量。
以上所述，仅为本发明的较佳实施方式而已，并非用于限定本发明的保护范围。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (12)

1.一种分段读出扩散加权成像方法，其特征在于，包括：
针对磁化强度矢量施加非线性相位射频激发脉冲，并在选层方向上施加持续时间与所述非线性相位射频激发脉冲相对应的选层梯度脉冲，以将磁化强度矢量翻转到X-Y平面；
针对翻转到X-Y平面上的磁化强度矢量执行扩散加权处理；
利用读出方向分段采样序列从扩散加权处理后的磁化强度矢量中分段读出成像数据，其中在选层方向上施加视角倾斜梯度脉冲。
2.根据权利要求1所述的分段读出扩散加权成像方法，其特征在于，所述非线性相位射频激发脉冲在时间域上为方波或近似方波。
3.根据权利要求1所述的分段读出扩散加权成像方法，其特征在于，所述非线性相位射频激发脉冲的非线性相位大小为K；
K≤3.6/(BW2*FTW)，其中BW为非线性相位射频激发脉冲的带宽，FTW为非线性相位射频激发脉冲的部分过渡带宽，FTW＝(ws-wp)/BW；ws是非线性相位射频激发脉冲的止带频率，wp是非线性相位射频激发脉冲的通带频率。
4.根据权利要求1所述的分段读出扩散加权成像方法，其特征在于，
所述针对翻转到X-Y平面上的磁化强度矢量执行扩散加权处理包括：针对翻转到X-Y平面上的磁化强度矢量施加扩散加权梯度脉冲，在射频方向上施加回聚脉冲，在选层方向上施加选层梯度脉冲和损毁梯度脉冲。
5.根据权利要求1所述的分段读出扩散加权成像方法，其特征在于，
所述利用读出方向分段采样序列从扩散加权处理后的磁化强度矢量中分段读出成像数据包括：在读出方向上施加成像回波读出梯度脉冲以读出成像数据，在读出方向上施加导航回波读出梯度脉冲以读出导航数据，用以校正分段采集的成像回波之间的相位误差，在相位编码方向上施加成像回波相位编码梯度脉冲和导航回波相位编码梯度脉冲。
6.根据权利要求1所述的分段读出扩散加权成像方法，其特征在于，
所述在选层方向上施加视角倾斜梯度脉冲包括：在选层方向上导航回波的回聚脉冲的两边施加视角倾斜梯度脉冲，所述视角倾斜梯度脉冲包括成像回波的视角倾斜梯度脉冲和导航回波的视角倾斜梯度脉冲；该方法还包括：
在选层方向上施加成像回波的视角倾斜梯度脉冲的预散相梯度脉冲、成像回波的视角倾斜梯度脉冲的聚相梯度脉冲和导航回波的视角倾斜梯度脉冲的预散相梯度脉冲。
7.一种分段读出扩散加权成像装置，其特征在于，包括：
翻转模块，用于针对磁化强度矢量施加非线性相位射频激发脉冲，并在选层方向上施加持续时间与所述非线性相位射频激发脉冲相对应的选层梯度脉冲，以将磁化强度矢量翻转到X-Y平面；
扩散加权模块，用于针对翻转到X-Y平面上的磁化强度矢量执行扩散加权处理；
分段读出模块，用于利用读出方向分段采样序列从扩散加权处理后的磁化强度矢量中分段读出成像数据，其中在选层方向上施加视角倾斜梯度脉冲。
8.根据权利要求7所述的分段读出扩散加权成像装置，其特征在于，所述非线性相位射频激发脉冲在时间域上为方波或近似方波。
9.根据权利要求7所述的分段读出扩散加权成像装置，其特征在于，
所述非线性相位射频激发脉冲的非线性相位大小为K；
K≤3.6/(BW2*FTW)，其中BW为非线性相位射频激发脉冲的带宽，FTW为非线性相位射频激发脉冲的部分过渡带宽，FTW＝(ws-wp)/BW；ws是非线性相位射频激发脉冲的止带频率，wp是非线性相位射频激发脉冲的通带频率。
10.根据权利要求7所述的分段读出扩散加权成像装置，其特征在于，
扩散加权模块，用于针对翻转到X-Y平面上的磁化强度矢量施加扩散加权梯度脉冲，在射频方向上施加回聚脉冲，在选层方向上施加选层梯度脉冲和损毁梯度脉冲；和/或
分段读出模块，用于在读出方向上施加成像回波读出梯度脉冲以读出成像数据，在读出方向上施加导航回波读出梯度脉冲以读出导航数据，用以校正分段采集的成像回波之间的相位误差，在相位编码方向上施加成像回波相位编码梯度脉冲和导航回波相位编码梯度脉冲。
11.根据权利要求7所述的分段读出扩散加权成像装置，其特征在于，
分段读出模块，用于在选层方向上导航回波的回聚脉冲的两边施加视角倾斜梯度脉冲，所述视角倾斜梯度脉冲包括成像回波的视角倾斜梯度脉冲和导航回波的视角倾斜梯度脉冲；在选层方向上施加成像回波的视角倾斜梯度脉冲的预散相梯度脉冲、成像回波的视角倾斜梯度脉冲的聚相梯度脉冲和导航回波的视角倾斜梯度脉冲的预散相梯度脉冲。
12.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现如权利要求1至6中任一项所述的分段读出扩散加权成像方法的步骤。
CN201710442188.9A 2017-06-13 2017-06-13 一种分段读出扩散加权成像方法、装置及可存储介质 Active CN109143134B (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN201710442188.9A CN109143134B (zh) 2017-06-13 2017-06-13 一种分段读出扩散加权成像方法、装置及可存储介质
EP18177302.9A EP3415941A1 (en) 2017-06-13 2018-06-12 Readout-segmented diffusion-weighted mri
US16/007,332 US10928476B2 (en) 2017-06-13 2018-06-13 Readout-segmented diffusion-weighted imaging method, apparatus, and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710442188.9A CN109143134B (zh) 2017-06-13 2017-06-13 一种分段读出扩散加权成像方法、装置及可存储介质

Publications (2)

Publication Number Publication Date
CN109143134A true CN109143134A (zh) 2019-01-04
CN109143134B CN109143134B (zh) 2021-01-29

Family

ID=62620762

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710442188.9A Active CN109143134B (zh) 2017-06-13 2017-06-13 一种分段读出扩散加权成像方法、装置及可存储介质

Country Status (3)

Country Link
US (1) US10928476B2 (zh)
EP (1) EP3415941A1 (zh)
CN (1) CN109143134B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113805133A (zh) * 2021-09-13 2021-12-17 上海交通大学 一种减少图像失真的磁共振平面回波成像方法

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102020209382A1 (de) * 2020-07-24 2022-01-27 Siemens Healthcare Gmbh Verfahren zur Aufnahme von Messdaten mittels einer Magnetresonanzanlage mit einer Korrektur der verwendeten k-Raumtrajektorien

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1499218A (zh) * 2002-10-31 2004-05-26 西门子公司 扩散加权磁共振成像中用稳态序列确定adc系数的方法
US20040257078A1 (en) * 2003-04-25 2004-12-23 Porter David Andrew Magnetic resonance imaging method and apparatus with spatial coding using readout segmentation
CN1683939A (zh) * 2004-04-13 2005-10-19 西门子公司 在磁共振断层造影中融合加权成像的运动校正的多点方法
US20070249929A1 (en) * 2006-03-31 2007-10-25 University Of Utah Systems and methods for magnetic resonance imaging
CN101143093A (zh) * 2006-09-11 2008-03-19 西门子（中国）有限公司 磁共振扩散成像方法
US20100237865A1 (en) * 2009-03-23 2010-09-23 Alto Stemmer Magnetic resonance method and apparatus to generate an image using a parallel acquisition technique
US20100237864A1 (en) * 2009-03-23 2010-09-23 Alto Stemmer Magnetic resonance method and apparatus to generate an image using a parallel acquisition technique
CN101879063A (zh) * 2009-05-04 2010-11-10 西门子公司 用于扩散加权地拍摄磁共振信号的方法和磁共振设备
CN102542534A (zh) * 2010-12-31 2012-07-04 北京海思威科技有限公司 基于图像轮廓的图像畸变校正方法和装置
CN102525469A (zh) * 2011-12-12 2012-07-04 中国科学院深圳先进技术研究院 针对小目标的磁共振成像方法和系统
CN102772210A (zh) * 2011-05-10 2012-11-14 西门子公司 弥散加权磁共振成像
US20150084629A1 (en) * 2013-09-24 2015-03-26 Siemens Aktiengesellschaft Method and magnetic resonance system for acquiring magnetic resonance data
CN104569881A (zh) * 2013-10-29 2015-04-29 西门子公司 利用扩散加权的磁共振序列记录磁共振数据的方法和设备

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108020796B (zh) * 2016-10-31 2021-10-29 西门子(深圳)磁共振有限公司 一种磁共振扩散加权成像方法和装置

Patent Citations (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1499218A (zh) * 2002-10-31 2004-05-26 西门子公司 扩散加权磁共振成像中用稳态序列确定adc系数的方法
US20040257078A1 (en) * 2003-04-25 2004-12-23 Porter David Andrew Magnetic resonance imaging method and apparatus with spatial coding using readout segmentation
CN1683939A (zh) * 2004-04-13 2005-10-19 西门子公司 在磁共振断层造影中融合加权成像的运动校正的多点方法
US20050237057A1 (en) * 2004-04-13 2005-10-27 Porter David A Movement-corrected multi-shot method for diffusion-weighted imaging in magnetic resonance tomography
US20070249929A1 (en) * 2006-03-31 2007-10-25 University Of Utah Systems and methods for magnetic resonance imaging
CN101143093A (zh) * 2006-09-11 2008-03-19 西门子（中国）有限公司 磁共振扩散成像方法
US20100237865A1 (en) * 2009-03-23 2010-09-23 Alto Stemmer Magnetic resonance method and apparatus to generate an image using a parallel acquisition technique
US20100237864A1 (en) * 2009-03-23 2010-09-23 Alto Stemmer Magnetic resonance method and apparatus to generate an image using a parallel acquisition technique
CN101879063A (zh) * 2009-05-04 2010-11-10 西门子公司 用于扩散加权地拍摄磁共振信号的方法和磁共振设备
CN102542534A (zh) * 2010-12-31 2012-07-04 北京海思威科技有限公司 基于图像轮廓的图像畸变校正方法和装置
CN102772210A (zh) * 2011-05-10 2012-11-14 西门子公司 弥散加权磁共振成像
CN102525469A (zh) * 2011-12-12 2012-07-04 中国科学院深圳先进技术研究院 针对小目标的磁共振成像方法和系统
US20150084629A1 (en) * 2013-09-24 2015-03-26 Siemens Aktiengesellschaft Method and magnetic resonance system for acquiring magnetic resonance data
CN104569881A (zh) * 2013-10-29 2015-04-29 西门子公司 利用扩散加权的磁共振序列记录磁共振数据的方法和设备

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
DAVID A. PORTER 等: "High Resolution Diffusion-Weighted Imaging Using Readout-Segmented Echo-Planar Imaging，Parallel Imaging and a Two-Dimensional Navigator-Based Reacquisition", 《MAGNETIC RESONANCE IN MEDICINE》 *
ROLF F. SCHULTE 等: "Equi-ripple design of quadratic-phase RF pulses", 《MAGNETIC RESONANCE》 *
SINYEOB AHN 等: "View Angle Tilting Echo Planar Imaging for Distortion Correction", 《MAGNETIC RESONANCE IN MEDICINE》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113805133A (zh) * 2021-09-13 2021-12-17 上海交通大学 一种减少图像失真的磁共振平面回波成像方法

Also Published As

Publication number Publication date
EP3415941A1 (en) 2018-12-19
US10928476B2 (en) 2021-02-23
CN109143134B (zh) 2021-01-29
US20180356487A1 (en) 2018-12-13

Similar Documents

Publication Publication Date Title
CN108020796B (zh) 一种磁共振扩散加权成像方法和装置
CN107126211B (zh) 磁共振成像方法和系统
CN104569881B (zh) 利用扩散加权的磁共振序列记录磁共振数据的方法和设备
US9529068B2 (en) System and method for image warp correction for magnetic resonance imaging
EP2588878A1 (en) Parallel magnetic resonance imaging using undersampled coil data for coil sensitivity estimation
US9103896B2 (en) Method and device for establishing excitation parameters for MR imaging
CN1102320A (zh) 梯度和自旋回波磁共振成像中的重影抑制
US9841481B2 (en) Multislice acquisition with incoherent aliasing
CN113313778B (zh) 磁共振图像的重建方法、计算机设备及存储介质
DE10318990A1 (de) Bildgebungsverfahren für die Magnetresonanz-Tomographie
CN106597337A (zh) 一种磁共振t2*加权快速成像方法及装置
CN109143134A (zh) 一种分段读出扩散加权成像方法、装置及可存储介质
US7508211B2 (en) Regularized species separation
Li et al. Imaging with referenceless distortion correction and flexible regions of interest using single-shot biaxial spatiotemporally encoded MRI
US8208986B2 (en) Steady state dark blood magnetic resonance imaging
CN107430179B (zh) 使用阵列空间伪灵敏度编码技术的多次激发扩散加权成像的设备和方法
Noll et al. A spectral approach to analyzing slice selection in planar imaging: optimization for through‐plane interpolation
US8625871B2 (en) Method for processing magnetic resonance images and apparatus for processing magnetic resonance images
CN107219484B (zh) 多层并行激发弥散成像的激发翻转角度确定方法和装置
CN108802649B (zh) 一种快速空间饱和的射频扰相方法和装置
KR101549184B1 (ko) 압축 센싱 지원 mmme 시퀀스를 이용한 초고속 3d 스핀 에코 영상법 및 그 영상장치
US10928473B2 (en) Methods and systems for reduced shading and blurring in magnetic resonance imaging
CN114062988B (zh) 磁共振波谱成像方法、装置、计算机设备和存储介质
US20130271131A1 (en) Varying Blipped Trajectory
US10139459B2 (en) Passive MRI-guided tracking of metallic objects using artifact-reduced imaging with selective excitation

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