CN115542217A - 用于回波平面成像的回波间隔混洗 - Google Patents

用于回波平面成像的回波间隔混洗 Download PDF

Info

Publication number
CN115542217A
CN115542217A CN202210740100.2A CN202210740100A CN115542217A CN 115542217 A CN115542217 A CN 115542217A CN 202210740100 A CN202210740100 A CN 202210740100A CN 115542217 A CN115542217 A CN 115542217A
Authority
CN
China
Prior art keywords
gradient
phase encoding
variable
est
magnetic resonance
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.)
Pending
Application number
CN202210740100.2A
Other languages
English (en)
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 Healthineers AG
University of Glasgow
Original Assignee
Siemens Healthineers AG
University of Glasgow
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 Healthineers AG, University of Glasgow filed Critical Siemens Healthineers AG
Publication of CN115542217A publication Critical patent/CN115542217A/zh
Pending legal-status Critical Current

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]
    • 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/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
    • 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/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56572—Correction of image distortions, e.g. due to magnetic field inhomogeneities caused by a distortion of a gradient magnetic field, e.g. non-linearity of a gradient magnetic field
    • 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/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/56366—Perfusion imaging

Landscapes

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

Abstract

本发明提供了用于回波平面成像的回波间隔混洗。描述了结合可变密度欠采样方案的EPI磁共振成像方法。该方法包括以下步骤：生成HF脉冲(HF)；以可变时间间隔(EST)应用切换的频率编码读出梯度(GR)；以及同时应用具有相位编码梯度(GP)的可变积分值的间歇性标志低幅度相位编码梯度(GP)。执行上述步骤，使得k空间至少部分地被欠采样，一个读出梯度(GR)的时间间隔(EST)根据相位编码梯度(GP)的积分而变化，使得读出梯度(GR)的可变时间间隔(EST)与对应的相位编码梯度(GP)的积分之间的比率被保持高于或处于预定的恒定值，该恒定值与图像质量的预定标准有关。此外，描述了MR控制序列确定装置(40)。此外，描述了MR成像系统(1)。

Description

用于回波平面成像的回波间隔混洗
技术领域
本发明涉及一种结合可变密度欠采样方案的EPI磁共振成像方法。此外，本发明涉及一种MR控制序列确定装置。此外，本发明涉及一种MR成像系统。
背景技术
在现代成像方法的帮助下，通常会生成二维或三维图像数据，这些二维或三维图像数据可以用于对要成像的患者例如人或动物进行可视化，并且也可以用于其他应用。
使用磁共振成像(缩写为MRI)实现一种特殊类型的医学成像。
磁共振成像是一种提供生物组织的诊断相关信息的医学成像技术。其诊断相关性源于这种成像方式的信息能力和多功能性，它以非侵入性方式提供差异化和精确的结构信息。磁共振成像基于对患者体内核自旋的受控操纵以及随后对核自旋响应的检测。自旋响应的空间相关编码允许重建患者的结构组成和功能构成，以供放射科医师执行诊断读取。
在磁共振系统中，待检查的身体通常借助基本场磁体系统暴露于较高的基本场磁场，例如1.5特斯拉、3特斯拉或7特斯拉。在施加基本场之后，检查对象中的原子核将自身与非零核磁偶极矩(通常也被称为自旋)沿基本场对齐。自旋系统的这种集体行为用宏观的“磁化”来描述。宏观磁化是对象在某一位置的所有微观磁矩的矢量和。除了基本场之外，还使用梯度系统施加磁场梯度，借助于磁场梯度，相应位置处的磁共振频率(拉莫尔(Larmor)频率)被确定。然后经由使用合适的天线装置的高频传输系统来传输高频激励信号(RF/HF脉冲)，这将导致被限定区域周围的该高频场(即，在相应位置处存在的拉莫尔频率下)激励的某些原子核的核自旋的翻转角相对于基本磁场的磁场线倾斜。如果这样的RF/HF脉冲作用于已经被激发的自旋，则它们可以被翻转到不同的角位置，甚至被折叠回平行于基本磁场的初始状态。在激发的核自旋的弛豫期间，所谓的磁共振信号的高频信号被共振发射，借助于合适的接收天线(也被称为磁共振线圈或接收线圈)被接收，然后被解调和数字化，并且然后作为所谓的“原始数据”被进一步处理。磁共振信号的获取发生在空间频率空间，即所谓的“k空间”中，其中在例如切片的测量期间，随着时间沿由梯度脉冲的切换限定的“梯度轨迹”(也被称为“k空间轨迹”)遍历k空间。此外，RF/HF脉冲必须在合适的时间被发出。从以这种方式获取的原始数据，在通常也取决于获取方法的另外的处理步骤之后，最终可以借助于二维傅里叶变换来重建期望的图像数据。
替选地，现在也可以以限定的方式激发和读取三维体积，在进一步的处理步骤之后将原始数据依次分类到三维k空间中。然后可以使用三维傅里叶变换重建三维图像数据体。
通常在测量期间使用某些预定的脉冲序列(即，以不同方向和读出窗口限定的RF/HF脉冲和梯度脉冲的序列)来控制磁共振断层扫描系统，在此期间，接收天线被切换以接收并且磁共振信号被接收和处理。在所谓的测量协议的帮助下，这些序列被预先参数化以用于期望的检查，例如计算出的图像的特定对比度。测量协议还可以包含用于测量的另外的控制数据。有许多磁共振序列技术，根据这些磁共振序列技术可以构建脉冲序列。磁共振成像(MR成像)未来发展面临的巨大挑战之一是加速磁共振序列技术，而不会在分辨率、对比度和对伪影(artifact)的敏感性方面做出广泛妥协。
当前的临床MR成像几乎完全基于所谓的笛卡尔或直角(直线)成像，其中扫描的k空间点(即，记录原始数据的k空间中的扫描点)位于矩形网格或图案的网格点上。
在这样做的过程中，所谓的并行成像方法已经成功地显著加速了临床MR成像。在并行MR成像的特殊变体(被称为压缩感知)的情况下，数据获取被缩短，因为在k空间中没有获取实际上对于重建无卷积图像所必需的一些网格线。这些缺失的线稍后会在k空间中的图像重建期间被替代，或者在图像空间中去除由欠采样产生的卷积伪影。能够使用并行成像方法的先决条件是利用数个接收线圈(天线)接收高频信号，其中必须知道各个接收线圈的空间灵敏度。在所谓的线圈校准数据的帮助下计算接收线圈的空间灵敏度。通常必须对线圈校准数据进行充分采样。由于灵敏度通常在空间中缓慢地变化，因此如果线圈校准数据具有低空间分辨率，那么这通常就足够了。通常，必须针对每位患者重新测量线圈校准数据。
最重要的并行成像方法之一是所谓的GRAPPA方法，它例如在以下文献中被描述：Marc Griswold等人，文章“Generalized Autocalibrating Partially ParallelAcquisitions(GRAPPA)”，医学磁共振(Magnetic Resonance in Medicine)，2002年第47期，1202页至1210页。
对于一般的并行成像或者特别是压缩感知，可变密度欠采样方案有利于加速的MRI中的图像质量。使用这种方法，k空间的中心部分被更密集地采样，而外围的更多线或点被忽略。将压缩感知与基于回波平面成像(被缩写为EPI)的快速获取技术相结合有潜在的好处。然而，标准EPI与可变密度欠采样不兼容，因为它会破坏非共振信号的相位演变，从而导致图像伪影。
利用EPI进行空间编码是通过在快速梯度切换期间对信号进行采样来实现的，这会在k空间中生成锯齿形轨迹。该轨迹在读出方向(kx)(被缩写为RO方向)与相位编码方向(ky)(被缩写为PE方向)之间存在差异。沿着RO方向，使用高梯度强度对k空间进行快速采样，而沿着PE方向，以与k空间中低采样率对应的规则间隔应用小的梯度标志(blip)。在组织界面处，例如在人脑的额叶皮层中的水和空气之间的边界处，可能会出现强烈的敏感性差异。这具有有助于背景磁场梯度的效果，其在幅度上可与标志PE梯度相比。这会扰乱从该组织区域发出的信号的k空间轨迹，从而导致最终图像中的空间失真。
为了减少这些失真，可以应用并行成像，其使用常规欠采样来缩短回波链(echotrain)，从而缩短在PE方向上遍历k空间所需的时间。这有效地增加了相对于背景梯度的标志PE梯度的幅度并降低了它们的影响。类似地，也可以使用多次激发(multi-shot)、分段k空间技术来减少回波链的长度，尽管以扫描时间的整体增加为代价。
存在以下问题：上述EPI方法依赖于k空间中的定期采样来保持图像质量，并且这与CS对可变密度欠采样的要求不兼容。这是因为可变采样会导致与非共振信号的恒定相位演化不匹配的k空间中的不均匀步长，这在具有背景磁场梯度的区域中是特别的问题。实际上，这会产生作为空间频率的函数的空间失真，从而导致最终图像中出现拖尾和细节丢失。
一种可能的方法是：始终使标志PE梯度的最小幅度高于某个阈值，使得在获取期间根本不会出现失真。高于该阈值，则可以实现可变密度欠采样模式。然而，由于需要较低的阈值，回波链的长度不能长于大约11个回波。这降低了整体获取速度，因此不会使该方法适用于成像速度至关重要的应用，比如功能MRI(fMRI)或扩散加权成像(DWI)。
因此存在开发具有足够图像质量的快速MR成像方法的问题。
发明内容
该目的通过根据本发明的技术方案的结合可变密度欠采样方案的EPI磁共振成像方法、根据本发明的技术方案的MR控制序列确定装置和根据本发明的技术方案的MR成像系统来实现。
根据本发明的结合可变密度欠采样方案的EPI磁共振成像方法包括以下步骤：生成HF脉冲；以可变时间间隔应用切换的频率编码读出梯度；以及同时应用具有相位编码梯度的可变积分值的间歇性标志低幅度相位编码梯度。可变时间间隔意味着所放出(playout)的读出梯度中的至少两个的时间间隔包括不同的值。相位编码梯度的可变积分值是指所放出的相位编码梯度中的至少两个相位编码梯度的积分值彼此不同。执行上述步骤，使得k空间至少部分地被欠采样；以及一个读出梯度的时间间隔根据相位编码梯度的积分值而变化，使得读出梯度的可变时间间隔与对应的相位编码梯度的积分值之间的比率被保持高于或处于预定的恒定值，该预定的恒定值与图像质量的预定标准有关。一个读出梯度的时间间隔也被称为回波间隔。相位编码梯度的积分值与k空间中的局部欠采样因子成正比。读出梯度的可变时间间隔与对应的相位编码梯度的积分值之间的比率被称为有效回波间隔。有效回波间隔也可以被定义为读出梯度的回波间隔与欠采样因子之间的比率。与欠采样因子对应的回波间隔的变化也被称为回波间隔混洗(shuffle)。
有利地，根据本发明的方法允许将欠采样方案与EPI获取方案相结合，EPI获取方案使用长回波链在每个激励处获取k空间的大区段。这样的结合进一步减少了扫描时间。由于有效回波间隔被保持高于预定值，因此图像质量被保持在预定水平，但是欠采样方案在整个k空间中发生变化。因此，可以使用与EPI脉冲序列结合的任意可变密度欠采样方案，而不会引入比预先预定的量更高的图像伪影。
根据本发明的MR控制序列确定装置包括：HF脉冲生成单元，其用于放出HF脉冲；梯度生成单元，其用于以可变时间间隔生成切换的频率编码读出梯度并同时生成具有相位编码梯度的可变幅度或可变积分值的间歇性标志低幅度相位编码梯度。梯度生成单元被设计成放出梯度，使得k空间至少部分地被欠采样，并且一个读出梯度的时间间隔根据相位编码梯度的积分值而变化，使得读出梯度的可变时间间隔与对应的相位编码梯度的积分值之间的比率被保持高于或处于预定的恒定值，该预定的恒定值与图像质量的预定标准有关。梯度生成单元共享根据本发明的EPI磁共振成像方法的优点。
MR成像系统包括根据本发明的扫描单元和MR控制序列确定装置。MR成像系统共享MR控制序列确定装置的优点。
根据本发明的MR控制序列确定装置的一些部件可以大部分以软件组件的形式被设计。这尤其适用于梯度生成单元和HF脉冲生成单元的部分。然而，原则上，这些部件中的一些也可以以软件支持的硬件的形式例如FPGA等来实现，尤其是当涉及特别快速的计算时。同样，所需的接口——例如如果仅仅是从其他软件组件传输数据的问题——可以被设计为软件接口。然而，它们也可以被设计为由合适的软件控制的基于硬件的接口。
很大程度上基于软件的实现方式具有以下优点：已经被使用的MR成像系统的计算机单元或控制单元可以通过软件更新容易地被改装，以便以根据本发明的方式工作。在这方面，该目的还通过具有计算机程序的相应计算机程序产品来实现，该计算机程序可以直接被加载到MR成像系统的计算机单元或控制单元的存储器装置中，该计算机程序包括程序段，以便当该计算机程序在MR成像系统的计算机单元或控制单元中被执行时，执行所有步骤以执行根据本发明的方法。
除了计算机程序之外，这样的计算机程序产品可以可选地包括附加组件例如文档和/或包括用于使用软件的硬件部件例如硬件密钥(加密狗等)的附加组件。
一种计算机可读介质例如记忆棒、硬盘或一些其他可移动或永久安装的数据载体，该计算机可读介质上存储有计算机程序的程序段，所述计算机程序可以被读取和执行，所述计算机程序可以用于传输到MR成像系统的计算机单元的存储装置和/或用于存储在MR成像系统的计算机单元上。例如，计算机单元可以为此目的而具有一个或更多个协作的微处理器等。
此外，本发明的特别有利的配置和发展由本发明的技术方案以及以下描述和附图而产生，由此一个权利要求类别中的独立权利要求也可以类似于另一权利要求类别中的从属权利要求及其描述部分来发展。
在根据本发明的EPI磁共振成像方法的优选变体中，相位编码梯度的幅度是可变的，并且一个读出梯度的时间间隔根据对应的相位编码梯度的幅度而变化，使得读出梯度的可变时间间隔与对应的相位编码梯度的幅度之间的比率被保持高于或处于预定的恒定值，该预定的恒定值与图像质量的预定标准有关。在前述变体中，相位编码梯度包括形状，其中，幅度与相位编码梯度的积分值成比例。例如，矩形或三角形形状的相位编码梯度就是这种情况。
在根据本发明的EPI磁共振成像方法的变体中，读出梯度的可变时间间隔与对应的相位编码梯度的积分值和/或幅度之间的比率在k空间的整个采样过程中被保持高于或处于预定的恒定值。有利地，失真的水平在整个k空间中被保持恒定以维持图像质量。
在根据本发明的EPI磁共振成像方法的优选变体中，欠采样方案包括可变密度方案，该可变密度方案被用来获取用于并行图像重建的原始数据。并行成像与EPI的结合使得能够进一步加速EPI获取。
在根据本发明的EPI磁共振成像方法的非常优选的变体中，欠采样方案包括可变密度方案，该可变密度方案被用来获取用于压缩感知重建的原始数据。压缩感知与EPI的结合还使得能够加速EPI获取过程。此外，这样的结合使得能够在短时间内对k空间的扩展区域进行采样，这对于诸如功能性磁共振成像或扩散加权磁共振成像方法的时间相关应用是非常重要的。
附图说明
下面基于示例性实施方式参考附图更详细地说明了本发明。在各个附图中，相同的部件具有相同的附图标记。
图1示出了不同EPI脉冲序列的一些脉冲序列图；
图2示出了与欠采样方案结合的标志相位编码梯度EPI方案的脉冲序列的一些k空间轨迹；
图3示出了根据本发明的实施方式的与可变密度欠采样方案结合的EPI磁共振成像方法的流程图；
图4示出了关于根据本发明的实施方式的MR控制序列确定装置的示意图；
图5示出了根据本发明的实施方式的MR成像系统的示意图；
图6示出了与压缩感知结合的EPI脉冲序列的脉冲序列图；以及
图7示出了与图6所示的脉冲序列对应的k空间轨迹。
具体实施方式
图1示出了根据现有技术的三个不同的脉冲序列图以及三个不同的EPI脉冲序列的对应k空间轨迹。第一脉冲序列表示具有连续相位编码梯度GP的EPI序列，其在图1的上部中被示出。脉冲序列图的右侧处的对应k空间轨迹示出频率和经编码的相位的线性变化。与此相比，图1中的中间位置所示的第二脉冲序列图包括具有所谓的标志梯度GP的脉冲序列，所谓的标志梯度GP具有频率的变化和相位的变化被去耦的效果，这在第二脉冲序列图的右侧上的k空间图中被示出。第三脉冲序列图示出了一些交替的相位编码梯度GP和读出梯度GR，它们导致k空间中的螺旋形轨迹，如对应脉冲序列图的右侧上所示。
图2示出了与欠采样方案结合的标志相位编码梯度EPI方案的脉冲序列的一些k空间轨迹。在左侧上，常规欠采样EPI方案被表示，其中，虚线表示忽略的读出轨迹线，实线表示EPI方案的采样的轨迹线。具有常规欠采样的压缩感知方案的问题是出现重复伪影，这会降低图像质量。
在图2的右侧上示出了具有可变欠采样的欠采样EPI方案，其中，忽略的读出轨迹线的数目从数目2变化到数目3。这样的可变欠采样的问题是：可变采样将导致与非共振信号的恒定相位演化不匹配的k空间中的不均匀步长，这对于具有背景磁场梯度的区域很重要。这种效果会产生作为空间频率的函数的空间失真，从而导致最终图像中出现拖尾和细节丢失。
图3示出了根据本发明的实施方式的与可变密度欠采样方案结合的EPI磁共振成像方法的流程图300。
在步骤3.I中，HF脉冲HF被放出。在步骤3.II中，以可变时间间隔放出经切换的频率编码读出梯度GR，并且在步骤3.III中，同时以可变幅度放出间歇性标志低幅度相位编码梯度GP，使得k空间至少部分地被欠采样，并且一个读出梯度GR的时间间隔根据相位编码梯度GP的积分值而变化。一个读出梯度GR的时间间隔被改变，使得读出梯度GR的可变时间间隔与对应相位编码梯度GP的积分值之间的比率被保持高于或处于预定的恒定值，该预定的恒定值与图像质量的预定标准有关。在步骤3.IV中，梯度回波被读出。
在图4中，示出了关于根据本发明的实施方式的MR控制序列确定装置40的示意图。
MR控制序列确定装置40包括用于生成HF脉冲的HF脉冲生成单元41。此外，MR控制序列确定装置40还包括用于生成相位编码脉冲序列GP的相位编码梯度生成单元42a和用于生成读出脉冲序列GR的读出梯度生成单元42b。
读出梯度生成单元42b在可变时间间隔内生成经切换的频率编码读出梯度，相位编码梯度生成单元42a同时生成具有可变幅度A的间歇性标志低幅度相位编码梯度GP，使得k空间至少部分地被欠采样，并且一个读出梯度GR的时间间隔EST根据相位编码梯度GP的积分值而变化，使得读出梯度GR的可变时间间隔EST与对应的相位编码梯度GP的积分值之间的比率被保持高于或处于预定的恒定值，该预定的恒定值与图像质量的预定标准有关。
在图5中，示意性地示出了根据本发明的实施方式的磁共振系统50(以下被简称为“MR系统”)。一方面，磁共振系统50包括具有检查室53或患者通道的实际磁共振扫描仪52，患者O或在此患者或测试主体可以在床58上进入检查室53或患者通道，在患者O或在此患者或测试主体的体内存在例如某个要被成像的器官。
磁共振扫描仪52以通常方式装配有基本场磁体系统54、梯度系统56以及RF发射天线系统55和RF接收天线系统57。在所示的示例性实施方式中，HF发射天线系统55是永久地安装在磁共振扫描仪52中的全身线圈，而HF接收天线系统57由布置在患者或测试者身上的局部线圈构成(在图5中仅由单个局部线圈表示)。然而，原则上，全身线圈55也可以被用作RF接收天线系统，并且局部线圈57也可以被用作RF发射天线系统，只要这些线圈可以各自被切换到不同的操作模式即可。
MR系统50还具有用于控制MR系统50的中央控制装置63。该中央控制装置63包括用于脉冲序列控制的序列控制单元64。这被用于根据选择的成像序列PS来控制高频脉冲(RF脉冲/HF脉冲)和梯度脉冲GR、GP的时间序列。例如，可以在测量或控制协议P内指定这样的成像序列。用于不同测量的不同控制协议P通常被存储在存储器69中并且可以由操作员选择(并且如果需要，则可以被改变)，以及然后被用于执行测量。序列控制单元64还包括根据本发明的控制序列确定装置40，该控制序列确定装置40具有图4所示的结构。控制序列确定装置40生成启用激励模块和读出模块的序列的控制数据SD，序列控制单元64输出该控制数据SD以用于脉冲序列控制。
为了输出单独的RF脉冲，中央控制装置63具有高频发射装置65，该高频发射装置65生成RF脉冲、对RF脉冲进行放大并经由合适的接口(未详细示出)将RF脉冲馈送到RF发射天线系统55中。为了控制梯度系统56的梯度线圈，中央控制装置63具有梯度系统接口66。序列控制单元64以合适的方式例如通过发出序列控制数据SD来进行通信，其中高频发射装置65和梯度系统接口66用于按照由控制序列确定装置40产生的顺序发出脉冲序列PS。控制装置63还具有高频接收装置67(同样以适当的方式与序列控制单元64进行通信)，以协调地获取由RF发射天线系统57接收到的磁共振信号。重建单元68将获取的数据在解调和数字化之后接受为原始数据RD并且从其重建MR图像数据。然后，这些图像数据BD例如可以被存储在存储器69中。
可以经由具有输入单元60和显示单元59的终端61来操作中央控制装置63，因此也可以由操作者经由终端61操作整个MR系统50。MR图像也可以被显示在显示单元59上，并且可以借助于输入单元60、可能地与显示单元59结合来计划和开始测量，并且特别地，可以如上所述选择具有合适测量序列的合适控制协议，并且必要时对这些合适控制协议进行修改。
根据本发明的MR系统50以及特别是中央控制装置63也可以具有大量的其他部件，这些部件在此没有详细地被示出但通常存在于诸如网络接口这样的装置上，以将整个系统连接至网络并能够交换原始数据RD和/或图像数据BD或参数卡以及其他数据，例如患者相关数据或控制协议。
对于本领域技术人员而言，如何通过照射RF/HF脉冲和产生梯度场来获取合适的原始数据RD以及可以如何从其重建MR图像BD原则上是已知的，因此在此不再对其进行更详细地说明。同样，对于本领域技术人员来说，用于生成动态或静态图像的多种多样的测量序列例如不同的EPI序列基本上是已知的。
在图6中，示出了结合压缩感知的EPI脉冲序列的脉冲序列图。在上面的行中，示出了取决于时间的HF脉冲HF，其与第三行中示出的读出脉冲GR结合地被放出。通过放出第二行中示出的具有第一低幅度A的标志相位编码梯度GP来生成第一回波，之后生成具有与相位编码梯度GP的低幅度成比例的短回波间隔时间的读出梯度GR。接下来，放出具有作为标志相位编码梯度GP的第一幅度的两倍高的幅度的标志相位编码梯度GP。然后以与所分配的相位编码梯度GP的幅度或积分值成比例的较长的回波间隔时间放出读出梯度GR。接下来，以低幅度放出标志相位编码梯度GP，之后以短回波间隔时间EST放出读出梯度GR。接下来，以作为先前标志相位编码梯度GP的低幅度的三倍高的高幅度放出标志相位编码梯度GP，之后，以作为先前读出梯度GR的回波间隔时间EST的三倍长的长回波间隔时间EST放出对应的读出梯度GR。最后，放出作为低幅度相位编码梯度GP的两倍高的标志相位编码梯度GP，并且放出作为短读出梯度GR的两倍长的对应读出梯度GR。因此，作为相位编码梯度GP的幅度或积分值与随后的读出梯度GR的回波间隔时间EST之间的比率的有效回波间隔是恒定的。恒定的有效回波间隔产生预定的图像质量。
在图7中，示出了与图6所示的脉冲序列对应的k空间轨迹。如从图7可以看出，通过忽略k空间的虚线来实现欠采样。此外，k空间轨迹的在ky方向即在相位编码方向上的一部分的长度与k空间轨迹在kx方向上的后续部分的长度成比例。由于欠采样变化，轨迹的在kx方向上的长度也变化，使得轨迹的在kx方向上的一部分的长度与轨迹的在ky方向上的前面部分的长度之间的比率是恒定的。这导致恒定的有效回波间隔。必须提及的是，有效回波间隔也可以高于预定恒定值，使得图像质量高于最低要求。实际上，较高的有效回波间隔会导致较长的回波间隔时间或降低的欠采样因子，这意味着成像过程需要更多时间和/或获取更多的原始数据。
最后再次指出，以上详细描述的方法和装置仅是示例性实施方式，本领域技术人员可以在不脱离本发明的范围的情况下以各种方式对其进行修改。此外，不定冠词“一”或“一个”的使用不排除相关特征出现不止一次的可能性。同样，术语“单元”不排除相关部件由数个相互作用的子部件构成的事实，这些子部件也可以空间地被分布。

Claims (11)

1.一种EPI磁共振成像方法，所述方法与可变密度欠采样方案结合并且包括以下步骤：
-生成HF脉冲(HF)，
-以可变时间间隔(EST)应用切换的频率编码读出梯度(GR)，以及
-同时应用具有相位编码梯度(GP)的可变积分值的间歇性标志低幅度相位编码梯度(GP)，使得
-k空间至少部分地被欠采样，以及
-一个读出梯度(GR)的时间间隔(EST)根据所述相位编码梯度(GP)的积分值而变化，使得所述读出梯度(GR)的可变时间间隔(EST)与对应的相位编码梯度(GP)的积分值之间的比率被保持高于或处于预定的恒定值，所述预定的恒定值与图像质量的预定标准有关。
2.根据权利要求1所述的EPI磁共振成像方法，其中，
-所述相位编码梯度(GP)的幅度是可变的，以及
-一个读出梯度(GR)的时间间隔(EST)根据所述相位编码梯度(GP)的幅度(A)而变化，使得所述读出梯度(GR)的可变时间间隔(EST)与所述对应的相位编码梯度(GP)的幅度(A)之间的比率被保持高于或处于预定的恒定值，所述预定的恒定值与图像质量的预定标准有关。
3.根据权利要求1或2所述的EPI磁共振成像方法，其中，所述读出梯度(GR)的可变时间间隔(EST)与所述对应的相位编码梯度(GP)的积分或幅度(A)之间的比率在k空间的整个采样过程中被保持高于或处于所述预定的恒定值。
4.根据前述权利要求中任一项所述的EPI磁共振成像方法，其中，所述欠采样方案包括沿相位编码方向的可变密度方案，所述可变密度方案被用来基于并行成像方案获取用于图像重建的原始数据。
5.根据权利要求1至3中任一项所述的EPI磁共振成像方法，其中，所述欠采样方案包括沿所述相位编码方向的可变密度方案，所述可变密度方案被用来基于压缩感知方案获取用于图像重建的原始数据。
6.根据权利要求1至5中任一项所述的EPI磁共振成像方法，包括功能性磁共振成像方法。
7.根据权利要求6所述的方法，包括扩散加权磁共振成像方法。
8.一种MR控制序列确定装置(40)，包括：
-HF脉冲生成单元(41)，所述HF脉冲生成单元用于生成HF脉冲，
-梯度生成单元(42a，42b)，所述梯度生成单元用于：
-以可变时间间隔(EST)生成切换的频率编码读出梯度(GR)以及
-同时生成具有相位编码梯度(GP)的可变积分值的间歇性标志低幅度相位编码梯度(GP)，使得
-k空间至少部分地被欠采样，以及
-一个读出梯度(GR)的时间间隔(EST)根据所述相位编码梯度(GP)的积分值而变化，使得所述读出梯度(GR)的可变时间间隔(EST)与对应的相位编码梯度(GP)的积分值之间的比率被保持高于或处于预定的恒定值，所述预定的恒定值与图像质量的预定标准有关。
9.一种MR成像系统(50)，包括：
-扫描单元(52)，
-根据权利要求8所述的MR控制序列确定装置(40)。
10.一种具有计算机程序的计算机程序产品，所述计算机程序能够直接被加载到MR成像系统(50)的控制单元(63)的存储单元中，所述计算机程序具有程序段，以便当所述计算机程序在所述控制单元(63)中被执行时执行根据权利要求1至7中任一项所述的方法的所有步骤。
11.一种计算机可读介质，所述计算机可读介质上存储有能够由计算机单元执行的程序段，以便当所述程序段由所述计算机单元执行时执行根据权利要求1至7中任一项所述的方法的所有步骤。
CN202210740100.2A 2021-06-30 2022-06-28 用于回波平面成像的回波间隔混洗 Pending CN115542217A (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
EP21182848.8 2021-06-30
EP21182848.8A EP4113149A1 (en) 2021-06-30 2021-06-30 Echo-spacing shuffling for echo-planar-imaging

Publications (1)

Publication Number Publication Date
CN115542217A true CN115542217A (zh) 2022-12-30

Family

ID=76730406

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210740100.2A Pending CN115542217A (zh) 2021-06-30 2022-06-28 用于回波平面成像的回波间隔混洗

Country Status (3)

Country Link
US (1) US11815583B2 (zh)
EP (1) EP4113149A1 (zh)
CN (1) CN115542217A (zh)

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10605883B2 (en) * 2016-04-22 2020-03-31 Sunnybrook Research Institute System and method for producing distortion free magnetic resonance images using dual-echo echo-planar imaging
KR20180085976A (ko) * 2017-01-20 2018-07-30 삼성전자주식회사 자기 공명 영상 장치 및 자기 공명 영상 장치의 제어 방법
EP3683595A1 (en) 2019-01-15 2020-07-22 Siemens Healthcare GmbH Segment-specific undersampling in readout-segmented magnetic resonance imaging
US11231476B2 (en) * 2019-04-26 2022-01-25 Regents Of The University Of Minnesota Accelerated magnetic resonance imaging acquisition using two-dimensional pulse segments as virtual receivers

Also Published As

Publication number Publication date
EP4113149A1 (en) 2023-01-04
US20230003820A1 (en) 2023-01-05
US11815583B2 (en) 2023-11-14

Similar Documents

Publication Publication Date Title
US9733328B2 (en) Compressed sensing MR image reconstruction using constraint from prior acquisition
US9575153B2 (en) MR imaging using a multi-point dixon technique
EP2193385B1 (en) Magnetic resonance using multi-dimensional rf excitation pulses
JP6270709B2 (ja) 磁気共鳴生データの再構成方法および装置
US8274286B2 (en) System and method for multi-spectral MR imaging near metal
US9977108B2 (en) Metal resistant MR imaging reference scan
US10338176B2 (en) Method and apparatus for actuation of a magnetic resonance scanner for the simultaneous acquisition of multiple slices
EP2496954B1 (en) Mr imaging using navigators
US7821264B2 (en) System and method for MR imaging in inhomogeneous magnetic fields
US9664762B2 (en) System and method for reduced field of view magnetic resonance imaging
US10203387B2 (en) MR imaging with enhanced susceptibility contrast
US20120046539A1 (en) Dual-contrast mr imaging using fluid-attenuation inversion recovery (flair)
EP2656095B1 (en) Fast dual contrast mr imaging
CN107076818A (zh) 零回波时间mr成像
US10732242B2 (en) T2-weighted MR imaging with elimination of non-T2-weighted signal contributions
US12019134B2 (en) MR electric properties tomography without contrast agent
US11815583B2 (en) Echo-spacing shuffling for echo-planar-imaging
WO2020173688A1 (en) Epi mr imaging with distortion correction
EP4137834A1 (en) Multi-region mr fingerprinting
WO2016124397A1 (en) Mr imaging with b1 mapping

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