CN110730912B - 并行多切片mr成像 - Google Patents

并行多切片mr成像 Download PDF

Info

Publication number
CN110730912B
CN110730912B CN201880038164.XA CN201880038164A CN110730912B CN 110730912 B CN110730912 B CN 110730912B CN 201880038164 A CN201880038164 A CN 201880038164A CN 110730912 B CN110730912 B CN 110730912B
Authority
CN
China
Prior art keywords
signals
phase
image
slice
encoding
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
CN201880038164.XA
Other languages
English (en)
Other versions
CN110730912A (zh
Inventor
M·P·J·于里森
A·F·N·格里戈莱托博尔戈诺维
T·H·罗济杰恩
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.)
Koninklijke Philips NV
Original Assignee
Koninklijke Philips NV
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 Koninklijke Philips NV filed Critical Koninklijke Philips NV
Publication of CN110730912A publication Critical patent/CN110730912A/zh
Application granted granted Critical
Publication of CN110730912B publication Critical patent/CN110730912B/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/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/483—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy
    • G01R33/4833—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy using spatially selective excitation of the volume of interest, e.g. selecting non-orthogonal or inclined slices
    • G01R33/4835—NMR imaging systems with selection of signals or spectra from particular regions of the volume, e.g. in vivo spectroscopy using spatially selective excitation of the volume of interest, e.g. selecting non-orthogonal or inclined slices of multiple slices
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/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
    • 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/565—Correction of image distortions, e.g. due to magnetic field inhomogeneities
    • G01R33/56509—Correction of image distortions, e.g. due to magnetic field inhomogeneities due to motion, displacement or flow, e.g. gradient moment nulling

Landscapes

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

Abstract

本发明涉及一种对放置在MR设备(1)的检查体积中的对象(10)进行MR成像的方法。所述方法包括以下步骤：通过使所述对象(10)经受包括多切片RF脉冲(21)的多回波成像序列的N次击发来生成MR信号，以同时激发两个或更多个空间上分离的图像切片，其中，切片方向上的相位偏移被赋予给MR信号，其中，所述相位偏移在击发之间变化，采集所述MR信号，其中，所述MR信号是经由在检查体积内具有不同空间灵敏度概况的至少两个RF线圈(11、12、13)的集合并行接收的，并且使用并行重建算法针对每个图像切片从采集的MR信号重建MR图像，其中，基于根据所述RF线圈(11、12、13)的空间灵敏度概况的MR信号的空间编码，并且基于归属于各图像切片和击发的所述相位偏移，分离来自不同图像切片的MR信号贡献。此外，本发明涉及用于执行该方法的MR设备以及要在MR设备上运行的计算机程序。

Description

并行多切片MR成像
技术领域
本发明涉及磁共振(MR)成像的领域。其涉及一种对对象进行MR成像的方法。本发明还涉及MR设备并且涉及要在MR设备上运行的计算机程序。
背景技术
当今广泛地使用图像形成MR方法，其利用磁场与核自旋之间的相互作用以形成二维或三维图像，特别是在医学诊断的领域使用，因为对于对软组织的成像，它们相对于其他方法在许多方面是有优势的，不需要电离辐射并且通常是非侵入性的。
根据一般的MR方法，对象，例如要被检查的患者的身体，被布置于强的均匀的磁场中，所述磁场的方向同时定义的测量所基于的坐标系的轴(通常是z轴)。磁场产生取决于磁场强度的针对个体核自旋不同的能级，所述能级可以通过施加具有限定频率(所谓的拉莫尔频率，或MR频率)的电磁交变场(RF场)而被激发(自旋共振)。从宏观的视角，个体核自旋的分布产生总体磁化，其可以通过施加合适的频率的电磁脉冲(RF脉冲)而被偏离出平衡态，使得自旋执行关于z轴的进动。进动描绘锥形的表面，其孔径角被称为翻转角。翻转角的幅度依赖于所施加的电磁脉冲的强度和持续时间。在所谓的90°脉冲的情况下，自旋被从z轴偏转到横向平面(翻转角90°)。
在RF脉冲结束后，磁化弛豫回初始的平衡态，其中，z方向的磁化以第一时间常数T1(自旋晶格弛豫或纵向弛豫时间)再次建立，并且在垂直于z方向的磁化以第二时间常数T2(自旋-自旋或横向弛豫时间)弛豫。磁化的变化可以借助于接收RF线圈检测到，其以如下的方式在MR设备的检查体积内被布置和取向：使得磁化的变化在垂直于z轴的方向被测量。横向磁化的衰减伴随有，例如，在施加90°脉冲之后，核自旋(由磁场不均匀性引起的)从具有相同相位的有序状态到所有相位角均匀地分布的状态(失相)的转变。所述失相可以借助于重新聚焦脉冲(例如，180°脉冲)来补偿。这在接收线圈中产生回波信号(自旋回波)。
为实现身体中的空间分辨，沿着主轴延伸的恒定磁场梯度被叠加到均匀磁场上，造成自旋共振频率的线性空间依赖性。在所述接收天线中拾取的信号则包括不同频率的分量，所述分量可以与所述身体/对象中的不同位置相关联。
经由所述接收线圈获得的信号数据对应于空间频率域，并且被称作k空间数据。所述k空间数据通常包括用不同的相位编码采集的多条线。每条线都通过收集若干样本进行数字化。k空间数据的集合借助于图像重建算法而被转换成MR图像。
并行采集技术通常用于加速MR采集。该类别中的方法是SENSE(灵敏度编码)、SMASH(空间谐波同时采集)和GRAPPA(广义自动校准部分并行采集)。SENSE、SMASH和GRAPPA以及其他并行采集技术使用从多个并行的RF接收线圈获得的欠采样k空间数据采集。在相应的重建算法中，来自多个线圈的(复数)信号数据与复数权重以抑制最终重建的MR图像中的欠采样伪影(混叠)的方式相组合。这种类型的复杂阵列信号组合有时称为空间滤波，并且包括在k空间域(如SMASH和GRAPPA)或图像域(如SENSE)执行的组合，以及混合的方法。
Larkman等人(Journal of Magnetic Resonance Imaging,13,313-317,2001)提出了在多切片成像的情况下也在切片方向上应用灵敏度编码以提高扫描效率。Breuer等人(Magnetic Resonance in Medicine,53,684-691,2005)提出的一种称为“controlledaliasing in parallel imaging results in higher acceleration”(CAIPIRINHA)的方法改进了此基本想法。该技术通过将相位偏移赋予每个图像切片的MR信号，从而在多切片采集过程中修改了每个个体切片中的混叠伪影的出现，从而改善了后续的并行图像重建过程。因此，CAIPIRINHA是一种并行多切片成像技术，与仅使用纯后处理方法的其他多层并行成像概念相比，其效率更高。在CAIPIRINHA中，使用相位调制的多切片RF脉冲(类似于已知的Hadamard脉冲)同时激发任意厚度和距离的多个切片。采集的MR信号数据是欠采样的，生成叠加的切片图像，这些切片图像表现为彼此相对偏移，对应于不同的相位偏移。根据傅立叶位移定理，通过由RF脉冲的相位调制方案施加的相位偏移来控制混叠切片图像的移位。通过使用此移位，可以改善反向重建问题的数值调节，所述调节将分离的片段的个体信号贡献分开。
US 2014/0225612 A1公开了一种使用分段回波平面成像(EPI)脉冲序列进行MR成像的方法。脉冲序列包括多个多切片RF脉冲，用于同时激发两个或多个空间上分离的图像切片。沿着条带编码方向应用特定的梯度编码方案，以将受控的相位偏移赋予不同的条带。使用并行成像重建方法将采集的MR信号数据重建为MR图像，所述方法将成像数据中的交叠切片分开，以针对跨成像对象的每个切片提供一系列MR图像。
US 2016/0018499 A1解决了CAIPIRINHA使用RF脉冲的固定相位调制方案的问题，使得相邻切片的相对位移为例如视场(FOV)尺寸的一半或FOV尺寸的其他整数分数。该固定方案的缺点是它不考虑先验信息。因此，没有充分考虑接收RF线圈的阵列的编码能力以及成像问题的基本结构，这可能导致次优的相位调制并因此导致次优的重建性能。作为此问题的解决方案，US 2016/0018499 A1建议使用可用的线圈灵敏度信息来导出经调整的切片特定的相位偏移，以进一步优化编码过程，并且因此优化展开问题，从而改善最终图像质量。
发明内容
根据本发明，公开了一种对放置在MR设备的检查体积中的对象进行MR成像方法。所述方法包括以下步骤：
-通过使所述对象经受包括多切片RF脉冲的多回波成像序列的N次击发来生成MR信号，以同时激发两个或更多个空间上分离的图像切片，其中，切片方向上的相位偏移被赋予给每个图像切片的MR信号，其中，所述相位偏移因击发而异，
-采集MR信号，其中，所述MR信号是经由在检查体积内具有不同空间灵敏度概况的至少两个RF线圈的集合并行接收的，并且
-使用并行重建算法针对每个图像切片从采集的MR信号重建MR图像，其中，基于根据RF线圈的空间灵敏度概况的MR信号的空间编码，并且基于归因于各图像切片和击发的相位偏移，分离来自不同图像切片的MR信号贡献。
与(单击发)EPI扫描相比，本发明的方法使用多击发多回波成像序列来改善图像分辨率和/或减少几何畸变和/或增加信噪比(SNR)。根据本发明，在逐击发的基础上改变多切片相位偏移。以这种方式，可以在切片方向上以灵活的方式应用额外的相位编码，这改善了重建步骤中切片图像的分离的稳定性。
本发明使用并行成像重建算法(例如，SENSE或SMASH重建)来从不同图像切片分离MR信号贡献。在多回波成像序列的每个击发中采集的成像信号数据是根据在切片方向上施加的不同相位偏移调制的完整k空间的子集。并行图像重建算法用于从信号数据生成完整(展开)的切片图像。MR信号是经由不同的RF线圈(“通道”)并行采集的，优选地利用在相位编码方向上的欠采样以加快采集速度。在图像重建期间，通过并行重建算法以常规方式考虑了这一点。通过将对应的相位信息额外地并入到并行重建算法所基于的编码矩阵中，可以考虑在击发之间在切片方向上应用的相位偏移。换句话说，根据本发明，使用编码矩阵来解决图像重建的反向问题，所述编码矩阵的矩阵元素不仅由RF线圈的空间灵敏度概况和成像序列的每次击发的k空间采样模式来确定而且也由归属于图像切片和击发的不同相位偏移来确定。这种方法可与已知的IRIS(使用图像空间采样功能的图像重建)重建方案相媲美，所述重建方案可在多击发多回波SENSE成像中进行有效的图像重建(参见Jeoung等人在Magnetic Resonance in Medicine 2013,第69卷,第793-802页)，不同之处在于，不仅考虑单个MR图像，而且考虑了多个图像切片，并且通过将相应的相位信息合并到编码矩阵中，从而将概念扩展到了多切片成像。
本发明的方法的成像序列是多击发序列，其需要以分段方式对k空间进行采样。换句话说，在成像序列的每次击发中仅对k空间的子集进行采样。值得注意的是，其中，采集在由变化的相位偏移赋予的相位编码方向上被欠采样的k空间数据。
在一个可能的实施例中，可变相位偏移可以借助于RF脉冲的相位调制来赋予(像，例如，在已知的POMP技术中，参见Glover等人，Journal of Magnetic Resonance Imaging，1991，第457-461页)。替代地，在切片选择方向上施加的磁场梯度可以用于赋予相位偏移(像在上面引用的US2014/0225612A1中)。
在成像序列的不同击发中采集的MR信号可能包含不同的相位误差，例如由于梯度线圈的加热或由于存在扩散梯度的被成像对象的运动。如果这些相位误差是先验已知的，则通过将相应的相位误差值合并到编码矩阵中，本发明的方法可以将它们考虑在内。类似地，例如由于相反梯度的不同延迟，奇数和偶数回波信号(在存在相反极性的读出磁场梯度的情况下获得)可能包含系统的相位误差。通过将每个完整击发中的回波信号的队列分成两个MR信号数据集，一个包括奇数回波信号，另一个包括偶数回波信号，从重建算法的角度来看，击发数量实际上增加了一倍。通过将相应的相位误差值并入到编码矩阵中，可以校正奇/偶回波相位误差。
为了考虑成像对象的运动，可以从成像对象中收集运动信息，例如在采集期间使用呼吸带，其中，基于运动信息的运动加权被应用于MR图像重建中。与以较少的运动获得的MR信号数据相比，在存在强烈运动的情况下获得的MR信号数据具有减小的加权。
可以通过使对象经受成像序列的击发之间的导航器序列来生成导航器信号。然后可以将从导航器信号导出的运动权重应用于MR图像重建中。替代地，类似于已知的IRIS方法，可以通过将源自导航信号的相应相位信息直接并入到编码矩阵中来考虑运动引起的逐击发相位误差。
在本发明的方法的另一优选实施例中，为了改善SNR，成像序列的M＜N次击发可以被执行两次或更多次。代替经由平均来重复完整采集以提高SNR，也可以仅重复N次击发总数的一个子集。然后可以将对应的MR信号添加到本发明的并行重建方案中。这些额外的击发可以采用与成像序列其他击发不同的相位偏移，以进一步提高重建算法的稳定性。
目前为止描述的本发明的方法可以借助于MR设备来执行，所述MR设备包括：至少一个主磁体线圈，其用于在检查体积内生成均匀的静态磁场；多个梯度线圈，其用于生成检查体积内的不同空间方向上的切换的磁场梯度；RF线圈的集合，其用于并行地接收来自身体的MR信号，所述RF线圈具有不同的空间灵敏度概况；控制单元，其用于控制RF脉冲的时间序列和切换的磁场场梯度；以及重建单元。本发明的方法例如可以通过对重建单元的对应的编程和/或MR设备的控制单元来实现。
本发明的方法可以有利地在临床中当前使用的多数MR设备中实施。为此，仅需要使用控制MR设备的计算机程序，使得其执行本发明的以上解释的方法。所述计算机程序可以存在于数据载体上或者可以存在于数据网络上，使得能够被下载以安装在MR设备的控制单元中。
附图说明
随附附图公开了本发明的优选的实施例。然而，要理解，附图仅被设计用于于图示和说明的目的，并且不作为对本公开的限度的限定。在附图中：
图1示出了用于执行本发明的方法的MR设备；
图2图示了根据本发明所采用的多击发多回波成像序列的一次击发；
图3示意性地示出了根据本发明应用的k空间编码方案。
具体实施方式
参考图1，示出了MR设备1。所述设备包括超导的或常导的主磁体线圈2，使得沿着通过检查体积的z轴创建基本上均匀的、空间上恒定的主磁场。
磁共振生成和操纵系统应用一系列RF脉冲和切换的磁场梯度来反转或激发核磁自旋、诱发磁共振、重新聚焦磁共振、操纵磁共振、空间地或者以其他方式对磁共振进行编码、使自旋饱合，等等，以执行MR成像。
更具体地，梯度放大器3将电流脉冲应用到沿着检查体积的x、y和z轴的全身梯度线圈4、5和6中选定的全身梯度线圈。数字RF频率发射器7经由发送/接收开关8来将RF脉冲或脉冲包发送到全身体积RF线圈9以将RF脉冲发送到检查体积。典型的MR成像序列包括短持续时间的RF脉冲分段的包，其与任何所应用的磁场梯度一起来实现对核磁共振的选定操纵。该RF脉冲被用于饱和、激发共振、反转磁化、重新聚焦共振或者操纵共振并且选择定位在检查体积中的身体10的部分。
为了借助于并行成像来生成身体10的有限区域的MR图像，邻近针对成像所选择的区域来放置一组本地RF线圈11、12、13。
所得的MR信号由RF线圈11、12、13拾取，并由优选地包括前置放大器(未示出)的接收器14解调。接收器14经由发送/接收开关8连接到RF线圈9、11、12和13。
主计算机15控制梯度脉冲放大器3和发射器7来生成多个MR成像序列中的任何一个，例如回波平面成像(EPI)、回波体积成像、梯度和自旋回波成像、快速自旋回波成像等等。针对选定的序列，接收器14在每个RF激发脉冲之后快速地相继接收单个或多个MR数据线。数据采集系统16执行对接收到信号的模数转换并且将每个MR数据线转换为适于进一步处理的数字格式。在现代MR设备中，数据采集系统16是独立的计算机，其专用于采集原始图像数据。
最终，数字原始图像数据通过应用傅立叶变换或其他合适的重建算法的重建处理器17而被重建为图像表示。MR图像可以表示穿过患者的平面切片、平行平面切片的阵列、三维体积等等。图像然后被存储在图像存储器中，其中，它可以被访问以用于例如经由提供得到的MR图像的人类可读的显示的视频监视器18来将切片、投影或者图像表示的其它部分转换为用于可视化的适当格式。
主计算机15和重建处理器17被编程为运行上面和下面描述的本发明的方法。
继续参考图1并且进一步参考图2和3，解释了本发明的成像技术的实施例。
如图2中所示，身体10经受多回波成像序列的多次击发。所述成像序列是经修改的多击发EPI序列。该序列包括用于同时激发两个或更多个空间上分离的图像切片的切片选择性多切片RF脉冲21。在相应的切片选择磁场梯度22的存在下辐射RF脉冲21，以便在若干成像切片中产生横向核磁化。在激发RF脉冲21之后，通过在交替的读出磁场梯度23的存在下对一系列梯度被称为回波的信号进行采样来采集MR信号。MR信号的相位编码是通过一系列相位编码磁场梯度“blip”24来执行的。在每个相位编码梯度blip 24的应用过程中，还会在切片编码方向上施加另一个磁场梯度blip25。根据本发明，使用切片编码梯度blip 25来向每个图像切片的MR信号赋予相位偏移。尽管可以使用梯度blip 25来施加不同的相位偏移，但是也可以通过连续击发的切片选择RF脉冲21的相应相位来赋予相位偏移。通过SENSE编码来加速MR信号的采集。经由具有不同空间灵敏度的RF接收线圈11、12、13并行地采集MR信号，并且通过在相位编码方向上的欠采样来采集MR信号。成像序列的一个击发如图2中所示。多回波序列的多个(N个)击发用于对k空间进行完全采样，以便能够使用SENSE重建从采集的成像回声信号数据重建完整的MR切片图像，其中，根据本发明，赋予不同信号的MR信号的相位偏移随击发而变化。换句话说，kz相位编码在击发之间变化，并且针对每个kz相位编码进行各种ky相位编码。在k空间中的每个(kz，ky)位置，在频率编码方向上对整个kx线进行采样。
使用并行重建算法从采集的MR信号数据重建每个图像切片的MR图像，其中，基于根据RF线圈的空间灵敏度概况的MR信号的空间编码，并且基于归因于各图像切片和击发的相位偏移，分离来自不同图像切片的MR信号贡献。下面将对此进行详细说明：
在具有N次击发的多击发序列中，其中，每个击发表表k空间的常规欠采样，可以将SENSE方程组设置为：
Figure BDA0002308760080000081
其中矢量
Figure BDA0002308760080000082
包括需要根据经由单独的RF线圈(“通道”)测得的混叠像素值计算出的图像像素值。混叠像素值由矢量
Figure BDA0002308760080000083
描述。S是线圈灵敏度矩阵，其由所使用的RF线圈的空间灵敏度概况确定。灵敏度矩阵S的大小取决于击发的欠采样、使用的RF线圈数以及由多切片RF脉冲激发的切片数。在根据本发明的成像序列的每一次击发中施加的相位偏移可以通过矩阵Φsh来考虑，其是描述在切片方向编码所得到的额外相位的对角矩阵。SENSE方程可在此基础上写为：
Figure BDA0002308760080000091
这里，Φsh是对角矩阵，其包含每个位置的相位编码：
Figure BDA0002308760080000092
不仅必须考虑在切片方向(kz)上的不同相位编码，而且还必须考虑在相位编码方向(ky)上的不同相位编码。因此，
Figure BDA0002308760080000097
描述了通过在常规相位编码方向(ky)上的相位编码和在切片方向(kz)上的相位偏移两者赋予给各个像素值的相位。线圈灵敏度编码和相位编码(y和z)可以组合为一个编码矩阵Esh：
SΦsh＝Esh
最后，可以将所有N次击发的公式组合在一个广义的SENSE重建内核中：
Figure BDA0002308760080000093
其中，
Figure BDA0002308760080000094
包括最终切片图像的所有像素值，这些最终切片图像是由包括切片方向上的相位编码的多层切片SENSE重建得到的。
Figure BDA0002308760080000095
的最小二乘解(无噪声去相关和正则化)是：
Figure BDA0002308760080000096
本发明的多层采集方法可以看作是三维扫描，其中，通过相应的相位编码对ky-kz空间进行采样，其中，kz线的数量等于同时激发的切片的数量。利用本发明的多击发方法，变得有可能以不同的kz编码多次采集给定的ky编码步骤。这使得可以最佳地对三维k空间进行采样(使用欠采样)。图3示出了根据本发明的ky-kz采样的示例，其中，八个切片同时被激发。对EPI序列的四个击发A，B，C和D施加了沿切片方向施加的四个不同相位偏移。欠采样是在第ky和kz方向两者上施加，在
Figure BDA0002308760080000101
方向具有交替的采样。在所描绘的实施例中，每个击发应用恒定的kz编码，使得在单次击发的EPI序列内不需要在切片方向上施加额外的梯度blip。

Claims (13)

1.一种对放置在MR设备(1)的检查体积中的对象(10)进行MR成像的方法，所述方法包括以下步骤：
-通过使所述对象(10)经受包括多切片RF脉冲(21)的多回波成像序列的N次击发来生成MR信号，以同时激发两个或更多个空间上分离的图像切片，其中，相位偏移被赋予给每个图像切片的所述MR信号，其中，所述相位偏移在击发之间变化，
-采集所述MR信号，其中，所述MR信号是经由在所述检查体积内具有不同空间灵敏度概况的至少两个RF线圈(11、12、13)的集合并行接收的，并且
-使用并行重建算法针对每个图像切片从采集的MR信号重建MR图像，其中，基于根据所述RF线圈(11、12、13)的空间灵敏度概况的所述MR信号的空间编码，并且基于归属于各图像切片和击发的所述相位偏移，分离来自不同图像切片的MR信号贡献。
2.根据权利要求1所述的方法，其中，在所述成像序列的每次击发中，采集在相位编码方向上欠采样的k空间数据。
3.根据权利要求1所述的方法，其中，在所述成像序列的每次击发中，采集由变化的相位偏移赋予的在相位编码方向上欠采样的k空间数据。
4.根据权利要求1所述的方法，其中，所述相位偏移是借助于所述RF脉冲的相位调制来赋予的。
5.根据权利要求1所述的方法，其中，所述相位偏移是借助于在切片选择方向上施加的磁场梯度(25)来赋予的。
6.根据权利要求1所述的方法，其中，MR图像重建的逆问题通过使用编码矩阵来解决，其中，所述编码矩阵的矩阵元素由所述RF线圈(11、12、13)的空间灵敏度概况、所述成像序列的每次击发的k空间采样、以及归属于所述图像切片和击发的所述相位偏移来确定。
7.根据权利要求6所述的方法，其中，通过将对应的相位误差值并入到所述编码矩阵中来考虑所述MR信号的先验已知相位误差。
8.根据权利要求1-7中的任一项所述的方法，其中，基于运动信息的运动加权在所述MR图像重建中被应用。
9.根据权利要求8所述的方法，其中，通过使所述对象(10)经受所述成像序列的所述击发之间的导航器序列来生成导航器信号，其中，所述运动信息是从所述导航器信号导出的。
10.根据权利要求6或7所述的方法，其中，导航器信号是通过使所述对象(10)经受导航器序列来生成的，其中，通过将从所述导航器信号导出的相位信息并入到所述编码矩阵中，来考虑运动引起的相位误差。
11.根据权利要求1-7中的任一项所述的方法，其中，出于改善信噪比的目的，成像序列的M次击发被执行两次或更多次，其中，M＜N。
12.一种用于执行根据权利要求1-11中的任一项所述的方法的MR设备，所述MR设备(1)包括：至少一个主磁体线圈(2)，其用于在检查体积内生成均匀的静态磁场；若干梯度线圈(4、5、6)，其用于在检查体积内在不同空间方向上产生切换的磁场梯度；至少两个RF线圈(11、12、13)的集合，其具有不同空间灵敏度概况；控制单元(15)，其用于控制RF脉冲和切换的磁场梯度的时间顺序；以及重建单元(17)，其中，所述MR设备(1)被配置为执行以下步骤：
-通过使所述对象(10)经受包括多切片RF脉冲(21)的多回波成像序列的N次击发来生成MR信号，以同时激发两个或更多个空间上分离的图像切片，其中，相位偏移被赋予给每个图像切片的所述MR信号，其中，所述相位偏移在击发之间变化，
-采集所述MR信号，其中，所述MR信号经由RF线圈(11、12、13)的所述集合而被并行地接收，并且
-使用并行重建算法针对每个图像切片从采集的MR信号重建MR图像，其中，基于根据所述RF线圈(11、12、13)的空间灵敏度概况的MR信号的空间编码，并且基于归属于各图像切片和击发的所述相位偏移，分离来自不同图像切片的MR信号贡献。
13.一种存储要在MR设备上运行的计算机程序的数据载体，所述计算机程序包括用于以下的指令：
-生成包括多切片RF脉冲(21)的多回波成像序列的N次击发，以同时激发两个或更多个空间上分离的图像切片，将相位偏移赋予给每个图像切片的MR信号，其中，所述相位偏移在击发之间变化，
-采集所述MR信号，并且
-使用并行重建算法针对每个图像切片从采集的MR信号重建MR图像，其中，基于根据至少两个RF线圈(11、12、13)的集合的空间灵敏度概况的空间编码，并且基于归属于各图像切片和击发的相位偏移，分离来自不同图像切片的MR信号贡献。
CN201880038164.XA 2017-06-09 2018-06-08 并行多切片mr成像 Active CN110730912B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
EP17175156.3 2017-06-09
EP17175156.3A EP3413074A1 (en) 2017-06-09 2017-06-09 Parallel multi-slice mr imaging
PCT/EP2018/065194 WO2018224661A1 (en) 2017-06-09 2018-06-08 Parallel multi-slice mr imaging

Publications (2)

Publication Number Publication Date
CN110730912A CN110730912A (zh) 2020-01-24
CN110730912B true CN110730912B (zh) 2022-06-14

Family

ID=59034562

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201880038164.XA Active CN110730912B (zh) 2017-06-09 2018-06-08 并行多切片mr成像

Country Status (5)

Country Link
US (1) US11112475B2 (zh)
EP (2) EP3413074A1 (zh)
JP (2) JP7309622B2 (zh)
CN (1) CN110730912B (zh)
WO (1) WO2018224661A1 (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3528002A1 (de) * 2018-02-20 2019-08-21 Siemens Healthcare GmbH Verfahren zur aufnahme eines magnetresonanzdatensatzes, computerprogrammprodukt, datenträger sowie magnetresonanzanlage
US11255941B2 (en) 2020-04-17 2022-02-22 Canon Medical Systems Corporation Apparatus and method using time shifts in magnetic resonance imaging (MRI) to implement controlled aliasing in parallel imaging (CAIPI)
DE102021203269A1 (de) * 2021-03-31 2022-10-06 Siemens Healthcare Gmbh Verfahren zur Reduzierung von Artefakten in mittels Magnetresonanz aufgenommenen Bilddatensätzen

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5245283A (en) * 1991-08-07 1993-09-14 Picker International, Inc. Technique for shifting out-of-slice artifacts to the edge of the field of view
US7977943B2 (en) * 2008-04-10 2011-07-12 General Electric Company Method and system for reconstructing images
CN104685368A (zh) * 2012-10-02 2015-06-03 皇家飞利浦有限公司 耐金属mr成像参考扫描
CN105051563A (zh) * 2013-03-15 2015-11-11 皇家飞利浦有限公司 使用相位调制rf脉冲的并行多切片mr成像
CN105143906A (zh) * 2013-01-25 2015-12-09 皇家飞利浦有限公司 金属抗mr成像
CN106030329A (zh) * 2013-12-12 2016-10-12 皇家飞利浦有限公司 使用多回波分段k空间采集的MR成像
CN106233154A (zh) * 2014-04-22 2016-12-14 皇家飞利浦有限公司 使用预脉冲和导航器的具有运动校正的磁共振成像

Family Cites Families (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE10152734B4 (de) * 2001-10-25 2005-12-29 Siemens Ag Gerät und Verfahren zur Magnet-Resonanz-Bildgebung bei gleichzeitiger Messung zweier benachbarter Schichten
US7408345B2 (en) 2006-02-06 2008-08-05 The Board Of Trustees Of The Leland Stanford Junior University Generalized MRI reconstruction with correction for multiple image distortion
JP5536358B2 (ja) 2009-03-31 2014-07-02 株式会社日立メディコ 磁気共鳴イメージング装置及び感度補正方法
WO2012054768A1 (en) * 2010-10-20 2012-04-26 Advanced Mri Technologies, Llc Multiplexed shifted echo planar imaging
WO2012085796A1 (en) 2010-12-22 2012-06-28 Koninklijke Philips Electronics N.V. Parallel mri method using calibration scan, coil sensitivity maps and navigators for rigid motion compensation
US8981776B2 (en) 2011-04-22 2015-03-17 The General Hospital Corporation Method for magnetic resonance imaging with controlled aliasing
US9886745B2 (en) 2012-06-28 2018-02-06 Duke University Multi-shot scan protocols for high-resolution MRI incorporating multiplexed sensitivity-encoding (MUSE)
CN103767705B (zh) * 2012-10-23 2017-12-22 三星电子株式会社 磁共振成像系统和磁共振成像方法
US9778336B2 (en) 2013-02-13 2017-10-03 The General Hospital Corporation System and method for rapid, multi-shot segmented magnetic resonance imaging
US9989610B2 (en) * 2013-10-24 2018-06-05 Siemens Healthcare Gmbh Multiband slice accelerated imaging with balanced slice-selective gradients
CN104714199B (zh) 2013-12-17 2018-04-24 西门子(深圳)磁共振有限公司 一种磁共振成像方法和装置
US10436866B2 (en) 2014-01-31 2019-10-08 The General Hospital Corporation Simultaneous multislice MRI with random gradient encoding
JP2017508556A (ja) 2014-03-28 2017-03-30 コーニンクレッカ フィリップス エヌ ヴェＫｏｎｉｎｋｌｉｊｋｅ Ｐｈｉｌｉｐｓ Ｎ．Ｖ． Ｓｅｎｓｅに関わるｅｐｉゴースト補正
US9841481B2 (en) 2014-04-16 2017-12-12 The Board Of Trustees Of The Leland Stanford Junior University Multislice acquisition with incoherent aliasing
US9897675B2 (en) * 2014-05-28 2018-02-20 Case Western Reserve University Magnetic resonance fingerprinting (MRF) with simultaneous multivolume acquisition
EP3256870B1 (en) 2015-02-10 2023-09-13 Voxelgrids Innovations Private Limited Method for constructing image from mri data
KR101733800B1 (ko) * 2015-09-18 2017-05-10 삼성전자주식회사 자기 공명 영상 장치가 혈관을 스캔하는 방법 및 그 자기 공명 영상 장치

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5245283A (en) * 1991-08-07 1993-09-14 Picker International, Inc. Technique for shifting out-of-slice artifacts to the edge of the field of view
US7977943B2 (en) * 2008-04-10 2011-07-12 General Electric Company Method and system for reconstructing images
CN104685368A (zh) * 2012-10-02 2015-06-03 皇家飞利浦有限公司 耐金属mr成像参考扫描
CN105143906A (zh) * 2013-01-25 2015-12-09 皇家飞利浦有限公司 金属抗mr成像
CN105051563A (zh) * 2013-03-15 2015-11-11 皇家飞利浦有限公司 使用相位调制rf脉冲的并行多切片mr成像
CN106030329A (zh) * 2013-12-12 2016-10-12 皇家飞利浦有限公司 使用多回波分段k空间采集的MR成像
CN106233154A (zh) * 2014-04-22 2016-12-14 皇家飞利浦有限公司 使用预脉冲和导航器的具有运动校正的磁共振成像

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Sevgi Gökçe Kafalı等.Simultaneous phase-correction and denoising for diffusion-weighted MRI.《2016 24th Signal Processing and Communication Application Conference》.2016,第1-4页. *
同步辐射显微CT的人关节软骨三维成像研究;刘成磊 等;《CT理论与应用研究》;20151130;第24卷(第6期);第793-799页 *

Also Published As

Publication number Publication date
EP3635427B1 (en) 2021-11-17
CN110730912A (zh) 2020-01-24
US20200088822A1 (en) 2020-03-19
EP3413074A1 (en) 2018-12-12
WO2018224661A1 (en) 2018-12-13
EP3635427A1 (en) 2020-04-15
JP7309622B2 (ja) 2023-07-18
US11112475B2 (en) 2021-09-07
JP2023109791A (ja) 2023-08-08
JP2020522344A (ja) 2020-07-30

Similar Documents

Publication Publication Date Title
US10401456B2 (en) Parallel MR imaging with Nyquist ghost correction for EPI
EP2411829B1 (en) Mr imaging using parallel signal acquisition
US9223001B2 (en) MR imaging using navigators
CN106796274B (zh) 具有伪迹抑制的propeller-mr成像
CN111758041A (zh) 使用多梯度回波序列进行的Dixon MR成像
WO2014141055A1 (en) Parallel multi-slice mr imaging using phase-modulated rf pulses
JP2023109791A (ja) 並列マルチスライスｍｒ撮像
US10895619B2 (en) MR imaging with Dixon-type water/fat separation
CN107810425B (zh) 消除非t2加权信号贡献的t2加权mr成像
CN109983358B (zh) Propeller mr成像
CN111656209A (zh) 利用具有固有运动校正的星形堆叠采集的mr成像
US11543482B2 (en) Magnetic resonance imaging using motion-compensated image reconstruction
US11914016B2 (en) Optimized k-space profile ordering for 3D radial MR imaging
EP2581756A1 (en) MR imaging using parallel signal acquisition
EP3185029A1 (en) Mr imaging using propeller acquisition with t2 decay correction
EP3432019A1 (en) Parallel multi-slice mr imaging using signal averaging
CN111164444A (zh) 具有经改进的脂肪位移校正的Dixon型水/脂肪分离MR成像
US11815577B2 (en) Parallel MR imaging using wave-encoding
EP4012434A1 (en) Dixon-type water/fat separation mr imaging
WO2018001759A1 (en) Diffusion weighted mr imaging using multi-shot epi with motion detection and modified sense reconstruction
US20220390538A1 (en) Dixon-type water/fat separation mr imaging
CN113544526A (zh) 使用波编码的并行mr成像

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