CN105051563A - 使用相位调制rf脉冲的并行多切片mr成像 - Google Patents

Abstract

本发明涉及一种对被放置在MR设备(1)的检查体积中的目标(10)进行MR成像方法。所述方法包括以下步骤：使目标(10)经受成像序列，所述成像序列包括用于同时激励两个或更多个空间分离的图像切片的相位调制多切片RF脉冲；采集MR信号，其中，MR信号是经由至少两个RF线圈(11、12、13)的集合而被并行接收的，至少两个RF线圈(11、12、13)具有在检查体积内的不同的空间灵敏度分布，以及根据采集到的MR信号来重建针对每个图像切片的MR图像，其中，来自不同的图像切片的MR信号贡献是基于至少两个RF线圈(11、12、13)的空间灵敏度分布并且基于RF脉冲的相位调制方案被分离的。为了优化对MR图像重建的逆问题的调节，从至少两个RF线圈(11、12、13)的空间灵敏度分布导出RF脉冲的相位调制方案。此外，本发明涉及一种用于执行所述方法的MR设备，并且涉及一种在MR设备上运行的计算机程序。

Description

使用相位调制RF脉冲的并行多切片MR成像

技术领域

本发明涉及磁共振(MR)成像领域。本发明关注一种对目标的MR成像的方法。本发明还涉及一种MR设备并且涉及一种在MR设备上运行的计算机程序。

背景技术

采用磁场与核自旋之间的相互作用以便于形成二维或三维图像的图像形成MR方法现今被广泛使用，特别是在医学诊断领域中，这是因为对于对软组织的成像而言所述图像形成方法在许多方面中优于其他成像方法，不要求电离辐射，并且通常是无创的。

发明内容

根据一般的MR方法，所述目标(例如待检查的患者的身体)被布置在强的、均匀的磁场中，所述磁场的方向同时限定测量所基于的坐标系的轴(一般为z轴)。磁场取决于磁场强度而产生针对各个核自旋的不同能量水平，所述磁场强度能够通过施加具有限定的频率(所谓拉莫尔频率，或MR频率)的电磁交变场(RF场)来激励(自旋共振)。从宏观的角度，各个核自旋的分布产生总体磁化，所述总体磁化能够通过施加具有适当频率的电磁脉冲(RF脉冲)而被偏转出平衡状态，使得所述磁化执行关于z轴的进动(precessional)运动。所述进动运动描述了圆锥的表面，所述圆锥的孔径角度被称为翻转角。翻转角的幅度取决于所施加的电磁脉冲的强度和持续时间。在所谓的90°脉冲的情况下，自旋从z轴偏转至横向平面(翻转角90°)。

RF脉冲终止之后，所述磁化弛豫回到原始的平衡状态，其中，在z方向上的磁化再次被建立为具有第一时间常数T₁(自旋晶格或纵向弛豫时间)，并且垂直于z方向的方向上的磁化以第二时间常数T₂(自旋-自旋或横向弛豫时间)来弛豫。能够借助于接收RF线圈来检测磁化的变化，所述接收RF线圈以这样的方式布置并定向在MR设备的检查体积内，所述方式为在垂直于z轴的方向上测量磁化的变化。在施加了例如90°脉冲之后，横向磁化的衰减伴随着(由局部磁场不均匀性诱发的)核自旋从具有相同相位的有序状态到所有相位角度是均匀分布的状态的转变(移相)。能够借助于再聚焦脉冲(例如180°脉冲)来对所述移相进行补偿。这在接收线圈中产生回波信号(自旋回波)。

为了实现在身体中的空间分辨率，沿三个主轴延伸的恒定磁场梯度被重叠在均匀磁场上，引起自旋共振频率的线性空间相关性。在接收线圈中拾取到的信号则含有能够与身体中的不同位置相关联的不同频率的分量。经由接收线圈获取的信号数据与空间频率域相对应，并且被称为k空间数据。k空间数据通常包括利用不同的相位编码采集的多条线。每行通过收集若干样本而被数字化。k空间数据集借助于图像重建算法而被转换成MR图像。

最近，已经开发出用于加速MR采集的技术，其被称为并行采集。该类方法是SENSE(灵敏度编码)、SMASH(空间谐波的同时采集)和GRAPPA(概括的自动校准部分并行采集)。SENSE，SMASH和GRAPPA以及其他并行采集技术使用从多个RF接收线圈并行获取的欠采样(undersampled)k-空间数据采集。在这些方法中，来自多个线圈的(合成)信号数据与合成加权以这样的方式，所述方式抑制最终重建的MR图像中的欠采样伪影(混叠)。这种类型的合成阵列信号组合有时被称为空间滤波，并且包括在k空间域中(如在SMASH和GRAPPA中)或在图像域中(如在SENSE中)执行的组合以及混合方法。

Larkman等人(JournalofMagneticResonanceImaging,13,313-317,2001)提出了在多切片成像的情况下也在切片方向中应用灵敏度编码，以提高扫描效率。Breuer等人(MagneticResonanceinMedicine,53,684-691,2005)改进了这种基本想法，提出名为“并行成像中的受控混叠产生较高的加速度”(CAIPIRINHA)的方法。这种技术修改了在多切片采集期间的每个个体切片中的混叠伪影的外观，改进了后续的并行图像重建流程。因此，CAIPIRINHA是与仅使用单纯的后处理方法的其他多切片并行成像概念相比更有效的并行多切片成像技术。在CAIPIRINHA中，任意厚度和距离的多个切片是利用相位调制多切片RF脉冲(类似于已知的Hadamard脉冲)来同时激励的。采集到的MR信号数据是欠采样的，产生相对于彼此移位而出现的重叠的切片图像。混叠切片图像的移位是根据傅立叶移位定理通过对RF脉冲的相位调制方案而被控制的。在相位编码步骤之间，多切片RF脉冲向每个切片的MR信号施加单个相位移位。通过使用这种移位来改进对分离所涉及的切片的个体信号贡献的逆向重建问题的数值调节。CAIPIRINHA具有也在切片彼此相当接近的这种情况下改进对重叠的切片图像的分离的潜力，从而使所使用的RF接收线圈的线圈灵敏度在待成像的个体切片中不会显著地不同。然而，CAIPIRINHA具有局限性。

CAIPIRINHA使用RF脉冲的固定相位调制方案，产生相邻切片的相对移位，例如视场(FOV)的尺寸的一半或FOV的尺寸的一些其他整数分数的移位。这种固定方案的缺点是，其不考虑先验的信息。因此，接收RF线圈的阵列的编码能力以及成像问题的基本结构未被充分考虑。这可能导致次优的相位调制，并因此导致次优的重建性能。

根据前述内容容易意识到，需要一种经改进的并行多切片MR成像技术。

根据本发明，公开了一种对被放置在MR设备的检查体积中的目标进行MR成像的方法。所述方法包括以下步骤：

-使所述目标经受成像序列，所述成像序列包括用于同时激励两个或更多个空间分离的图像切片的相位调制多切片RF脉冲，

-采集MR信号，其中，所述MR信号是经由至少两个RF线圈的集合而被并行接收的，所述至少两个RF线圈具有在所述检查体积内的不同的空间灵敏度分布，以及

–根据采集到的MR信号来重建针对每个图像切片的MR图像，其中，来自不同图像切片的MR信号贡献是基于所述至少两个RF线圈的所述空间灵敏度分布并且基于所述RF脉冲的相位调制方案而被分离的。本发明提出所述RF脉冲的所述相位调制方案是从所述至少两个RF线圈的所述空间灵敏度分布导出的。

本发明的技术总体上与已知的CAIPIRINHA方法相对应。借助于相位调制多切片RF脉冲的如以上所描述的图像切片之间的相对移位的应用增大了在重叠的切片图像中折叠在一起的体素的位置处的线圈灵敏度之间的实际差异，并因此改进了对重建问题的调节。然而，根据本发明，不应用如CAIPIRINHA中的固定相位调制方案。作为替代，可用的线圈灵敏度信息被用于导出经调节的切片特异相位调制，以优化编码过程，并因此优化对展开问题的调节，以便改进最终图像质量。

在本发明的优选实施例中，相位编码的MR信号是借助于所述成像序列而被采集的，所述相位调制方案包括变化的相移，使得个体相位循环被应用到每个图像切片的所述MR信号。优选地，所述相移是在相位编码步骤之间线性递增的，其中，个体相位增量被应用到每个图像切片。这样，每个切片图像的个体移位是根据傅立叶移位定理通过所述RF脉冲的所述相位调制方案来控制的。

根据本发明的优选实施例，所述MR信号是利用在所述图像切片的平面内方向上的欠采样来采集的。所述图像切片的所述MR图像能够在这种情况下通过使用如SENSE、SMASH或GRAPPA的本身已知的并行图像重建算法来重建。

以下更详细地对本发明的技术进行解释：

首先，在没有所述多切片RF脉冲的相位调制的情况下考虑在所有N个不同图像切片中的相同平面内位置r处的一个体素。M个接收线圈中的每个中的对应的采集到的MR信号能够以矩阵向量记法被写为

Sx＝y

其中，向量y表示M个RF线圈中的每个中的采集到的MR信号，如N个不同图像切片的灵敏度加权的信号贡献x_i的线性组合，而矩阵S表示(N*M)的灵敏度矩阵。作为范例，针对使用M＝4个线圈和N＝3个切片的示范性实施例，以上方程能够被写为

$\left(\begin{array}{ccc}{S}_{11}& S_{12}& S_{13}\\ S_{21}& S_{22}& S_{23}\\ S_{31}& S_{32}& S_{33}\\ S_{41}& S_{42}& S_{43}\end{array}\right)\left(\begin{array}{c}{x}_1\\ x_2\\ x_3\end{array}\right)=\left(\begin{array}{c}{y}_1\\ y_2\\ y_3\\ y_4\end{array}\right)$

其中，S_ij是针对线圈j和切片i的线圈灵敏度。对包括编码矩阵的逆矩阵的该方程组求解产生向量x，向量x含有对应的N个切片特异MR信号：

(S^HS)^-1S^Hy＝x

矩阵(S^HS)^-1S^H为S的伪逆，其范数描述了从MR信号采集到最终图像中的误差传播。该范数在良好调节的情况下是小的。

根据本发明，通过在针对不同图像切片的MR信号采集过程期间应用适当的相位调制来影响得到的条件数量，所述相位调制最终改变S的各个元素。实际上，只有(N*M-M)个元素必须被考虑，这是因为针对第一图像切片的相位调制不是必须被执行的，这意味着只考虑那些具有i>l的矩阵元素S_ij就够了。

(在该示范性实施例的三切片激励试验中的)切片相关的平面内移位Δ₁和Δ₂被适当地选择为在相位编码方向上的FOV的分数，以修改MR信号采集期间的实际线圈灵敏度：

$\left(\begin{array}{ccc}{S}_{11}\left(r\right)& S_{12}(r+{\mathrm{\Δ}}_1)& S_{13}(r+{\mathrm{\Δ}}_2)\\ S_{21}\left(r\right)& S_{22}(r+{\mathrm{\Δ}}_1)& S_{23}(r+{\mathrm{\Δ}}_2)\\ S_{31}\left(r\right)& S_{32}(r+{\mathrm{\Δ}}_1)& S_{33}(r+{\mathrm{\Δ}}_2)\\ S_{41}\left(r\right)& S_{42}(r+{\mathrm{\Δ}}_1)& S_{43}(r+{\mathrm{\Δ}}_2)\end{array}\right)\left(\begin{array}{c}{x}_1\left(r\right)\\ x_2(r+{\mathrm{\Δ}}_1)\\ x_3(r+{\mathrm{\Δ}}_2)\end{array}\right)=\left(\begin{array}{c}{y}_1\\ y_2\\ y_3\\ y_4\end{array}\right)$

换言之，目的是使S的伪逆中的误差传播最小化。为了实现该目的，需要在所有可能的设定中确定使(S^HS)^～1S^H的范数最小化那些Δ₁和Δ₂。然而，以上的方程仅描述了所考虑的多个图像切片的集合中的一个体素位置。最优移位参数Δ₁、Δ₂必须被确定为优选地针对图像切片中的所有体素的折衷。

出于该目的，将以上针对所有体素的方程中给出的单个体素的问题连接起来，能够写出线性方程组。得到的单个灵敏度矩阵S'包括根据移位参数A_d的针对所有体素的较小的矩阵S_K，S_K也能够被解释为向量。对应的向量元素x"和y'是由以上限定的向量x和y形成的：

$\left(\begin{array}{cccc}{S}_1(r_1,{\mathrm{\Δ}}_1,{\mathrm{\Δ}}_2)& & & \\ & S_2(r_2,{\mathrm{\Δ}}_1,{\mathrm{\Δ}}_2)& & \\ & & S_3(r_3,{\mathrm{\Δ}}_1,{\mathrm{\Δ}}_2)& \\ & & & ...\end{array}\right)\left(\begin{array}{c}{x^{\′}}_1\\ {x^{\′}}_2\\ {x^{\′}}_3\\ ...\end{array}\right)=S^{\′}{x}^{\′}=\left(\begin{array}{c}{y^{\′}}_1\\ {y^{\′}}_2\\ {y^{\′}}_3\\ ...\end{array}\right)$

计算出的伪逆(S'^HS')¹S'^H的适当的范数能够用作用来选择Δ_d的最优解的量度：

Δ_d＝argmin||(S'^HS')^-1S'^H||

得到的移位Δ_d能够根据本发明被转译为对应的多切片RF脉冲的相位调制方案。

矩阵S'是稀疏带或对角矩阵，针对所述稀疏带或对角矩阵存在有效的数值逆例程。然而，能够采取若干措施来降低数值努力。

在一个实施例中，如果不是图像切片中的所有体素都被用于优化所述相位调制方案，则可以减小矩阵S'的大小。出于该目的，潜在切片内位置密度能够被降低(低分辨率方法)，由此减少S'中的元素的数量。这也暗示了，所述线圈灵敏度是在空间中相当平滑的函数。这样，额外地，性能能够进一步被集中在高度地临床感兴趣(ROI)的图像的重要区域上。

此外，在从所述RF脉冲的所述空间灵敏度分布导出所述RF脉冲的所述相位调制方案的步骤中能够考虑关于所述目标的几何结构的现有知识。该现有知识例如能够是从在实际图像采集之前执行的SENSE参考扫描获得的。

用来减小优化问题的大小的另外的方法是使用例如奇异值分解(SVD)或类似的技术来将所有涉及的RF接收线圈变换为虚拟线圈系统。因此，只有N个(图像切片的数量)最重要的虚拟线圈元件必须在评估中被考虑。

然而，作为对所述伪逆的基本计算的元素，针对每个位置直接计算S^HS并使用所有RF线圈也是高效的，并且避免忽略RF线圈信息的部分。

根据本发明的对所述相位调制方案的优化也能够通过选择用于关于所选择的Δ_d来评估所述伪逆的高效矩阵范数来加速。可以在对矩阵(S'^HS')^-1S'^H的个体对角元素的评估期间应用统计手段。另外，能够以分层次的方式来执行对最优Δ_d的搜索，开始形成Δ_d空间中的粗分辨率，所述粗分辨率能够在优化进展期间被调整。为了避免在重建中需要重组(regridding)，可以有利地只考虑是体素大小的整数倍的移位。

在优选实施例中，典型地应用在SENSE重建方法中的额外的信息(例如噪声相关矩阵、正则化MR图像、或者在最简单的情况下的将所述线圈灵敏度分布限制到MR信号被期望来自于其中的区域的二进制掩模)也可以被包括在对所述相位调制方案的优化中的对伪逆矩阵的计算中。

在本发明的实际实施例中，所述相位调制方案必须针对每个位置，即每组同时被激励的图像切片而被分别确定。所述线圈灵敏度分布和所述目标几何结构能够针对不同的位置而显著地变化，导致最优相位调制和对应的移位的不同集合。

现在为止所描述的本发明的方法能够借助于MR设备来执行，所述MR设备包括：用于在检查体积内生成均匀静磁场的至少一个主磁体线圈、用于在所述检查体积内在不同空间方向上生成切换的磁场梯度的多个梯度线圈、用于并行地从所述身体接收MR信号的RF线圈的集合、具有不同空间灵敏度分布的所述RF线圈、用于控制RF脉冲和切换的磁场梯度的时间次序的控制单元、以及重建单元。例如，本发明的方法能够通过对所述MR设备的所述重建单元和/或所述控制单元的对应的编程来实现。

本发明的方法能够有利地在目前临床使用的大部分MR设备中执行。为此，只需要采用计算机程序，所述MR设备是通过所述计算机程序来控制的，使得所述MR设备执行本发明的以上解释的方法步骤。所述计算机程序可以被呈现在数据载体上，或者被呈现在数据网络中，从而被下载以用于在所述MR设备的所述控制单元中安装。

附图说明

附图公开了本发明的优选实施例。然而，应当理解，附图是仅被设计用于图示的目的，而不是对本发明的限制的限定。在附图中：

图1示出了用于执行本发明的方法的MR设备；

图2示意性地图示了本发明的MR图像采集方案。

具体实施方式

参考图1，示出了MR设备1。设备包括超导的或常导的主磁体线圈2，是的沿通过检查体积的z轴创建基本均匀的、暂时恒定的主磁场。

磁共振生成和操纵系统应用一系列RF脉冲和切换的磁场梯度来反转或激励核磁自旋、诱发磁共振、聚焦磁共振、操纵磁共振、在空间上并以其他方式对磁共振进行编码、使自旋饱和等来执行MR成像。

更具体地，梯度脉冲放大器3沿检查体积的x、y和z轴向全身梯度线圈4、5和6中的选定的一个施加电流脉冲。数字RF频率发射器7经由发送/接收开关8向全身体积RF线圈9发射RF脉冲或脉冲包，以将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成像序列中的任何序列，例如回波平面成像(EPF)、回波体积成像、梯度和自旋回波成像、快速自旋回波成像等。对于选定的序列，接收器14在每个RF激励脉冲之后紧接着接收单条或多条MR数据线。数据采集系统16执行对接收到的信号的模数转换，并将每条MR数据线转换为适合于进一步处理的数字格式。在现代MR设备中，数据采集系统16是专门采集原始图像数据的单独的计算机。

最终，数字原始图像数据由应用傅立叶变换或其他适当的重建算法的重建处理器17重建为图像表示。MR图像可以表示通过患者的平面切片、平行平面切片的阵列、三维体积等。图像接着被存储在图像存储器中，其中，所述图像存储器可以被访问，以例如经由提供得到的MR图像的人类可读显示的视频监视器18来将切片、投影或图像表示的其他部分转换为用于可视化的适当的格式。

继续参考图1并进一步参考图2，解释了本发明的成像方法的实施例。

图2示出了三种不同的MR信号采集方法的方案。目标10经受成像序列(例如自旋-回波序列，如已知的TSE序列)，所述成像序列包括用于同时激励具有任意切片厚度的三个空间独立的、共平面图像切片(切片1、切片2、切片3)的多切片RF脉冲。图2d示意性示出了多切片RF脉冲在切片选择方向z上的激励分布。Φι和Φ₂分别指示切片2和切片3中的RF激励的相移增量。切片1可以总体上被激励，而不需要相移(0°)。MR信号经由具有在切片选择的方向上的不同空间灵敏度的RF线圈11、12、13而被并行地接收。如图2中右侧所示，采集到的MR信号产生在“折叠图像”中的重叠的切片图像。针对每个图像切片的MR图像是根据采集到的MR信号来重建的，其中，来自不同图像切片的MR信号贡献是基于RF线圈11、12、13的空间灵敏度分布而被分离的。

图2a示出了标准的并行多切片方法(见Larkman等人的JournalofMagneticResonance,13,313-317,2001)。在这种情况下，相移增量为Φ₁＝Φ₂＝0°，即没有应用相位调制。右侧的折叠的MR图像示出了所有的切片图像显著地交叠在重叠的MR图像中，使得分离非常困难(即数值病态)

在图2b中，图示了常规CAIPIRINHA方法。在相位编码步骤之间，多切片RF脉冲向切片1、2和3中的每个的MR信号应用个体相移。对逆向重建问题的数值调节是通过使用这种移位来改进的，即所述移位由于“折叠”图像中的更分散的折叠表现而变得明显。CAIPIRINHA使用RF脉冲的固定相位调制方案，得到所描绘的实施例中的相邻切片的0、FOV的1/3、以及FOV的2/3的相对移位。对应的固定RF相移增量对于每个相位编码步骤为Φ₁＝120°、Φ₂＝240°。这种固定方案的缺点是其不考虑先验信息。因此，接收RF线圈的阵列的编码能力以及成像问题的基本结构没有被充分地考虑。这可能导致次优的相位调制，并因此导致次优的重建性能。

图2C示出了本发明的方法，在所述方法中没有应用如CAIPIRINHA中的固定相位调制方案。根据本发明，可用线圈灵敏度信息被用于导出经调节的切片特异相位调制，以优化编码过程，并因此优化对展开问题的调节，以改进最终的图像质量。在所描绘的实施例中，优化流程得到相邻切片的0、FOV的1.15/3、以及FOV的2.3/3的相对移位。对应的RF相移增量对于每个相位编码步骤为Φ₁＝140°、Φ₂＝280°，其也满足整数体素移位条件。能够利用由这两个相移的最小公倍数给出的周期性来计算对应的多切片RF脉冲的集合。接着在MR信号采集期间以如下的相继顺序来应用该集合(括号中的数字分别指示应用到切片1、2和3的相移)：

1.RF脉冲(0°、140°、280°)，

2.RF脉冲((0°、280°、200°)，

3.RF脉冲(0°、60°、120°)，

4.RF脉冲(0°、200°、40°)，

5.RF脉冲(0°、340°、320°)，

6.……

对于线性相位编码方案中的每个相继的相位编码步骤，对应的RF脉冲以轮转的次序循环。因此，对分离问题的调节被改进，这在图2c中被示意性地示出。同常规CAIPIRINHA方法(图2b)相比，减小了“折叠”图像中的交叠MR信号强度的量。

在本发明的另外的实施例中，可以执行如以上范例中的多切片梯度回波扫描，但具有使用适当的缩减因子R(例如R＝2)的、在平面内相位编码方向上对SENSE的额外应用。为了优化对应的相移增量，也可以考虑对MR信号的平面内SENSE-诱发折叠。

Claims (9)

1.一种对被放置在MR设备(1)的检查体积中的目标(10)进行MR成像的方法，所述方法包括以下步骤：

-使所述目标(10)经受成像序列，所述成像序列包括用于同时激励两个或更多个空间分离的图像切片的相位调制多切片RF脉冲，

-采集MR信号，其中，所述MR信号是经由至少两个RF线圈(11、12、13)的集合而被并行接收的，所述至少两个RF线圈(11、12、13)具有在所述检查体积内的不同的空间灵敏度分布，以及

-根据采集到的MR信号来重建针对每个图像切片的MR图像，其中，来自不同图像切片的MR信号贡献是基于所述至少两个RF线圈(11、12、13)的所述空间灵敏度分布并且基于所述RF脉冲的所述相位调制方案而被分离的，

其中，所述RF脉冲的所述相位调制方案是从所述至少两个RF线圈(11、12、13)的所述空间灵敏度分布导出的。

2.根据权利要求1所述的方法，其中，相位编码的MR信号是借助于所述成像序列而被采集的，所述相位调制方案包括变化的相移，使得个体相位循环被应用到每个图像切片的所述MR信号。

3.根据权利要求2所述的方法，其中，所述相移是在相位编码步骤之间线性递增的，其中，个体相位增量被应用到每个图像切片。

4.根据权利要求1-3中的任一项所述的方法，其中，所述RF脉冲的所述相位调制方案被确定为使得对所述MR图像的重建的逆问题的调节被优化。

5.根据权利要求1-4中的任一项所述的方法，其中，在从所述RF脉冲的所述空间灵敏度分布导出所述RF脉冲的所述相位调制方案的步骤中考虑所述目标(10)的几何结构的现有知识。

6.根据权利要求1-5中的任一项所述的方法，其中，所述相位调制方案是通过使范数(S^HS)^-1S^H最小化而被导出的，其中，S是所述至少两个RF线圈(11、12、13)的灵敏度矩阵，所述S考虑了针对由所述RF脉冲的相位调制诱发的每个图像切片的在所述相位编码方向上的个体空间移位。

7.根据权利要求1-6中的任一项所述的方法，其中，所述MR信号是利用在所述图像切片的平面内方向上的欠采样来采集的，其中，所述MR图像是使用如SENSE、SMASH或GRAPPA的并行图像重建算法来重建的。

8.一种用于执行根据权利要求1-5所述的方法的MR设备，所述MR设备(1)包括：用于在检查体积内生成均匀静磁场的至少一个主磁体线圈(2)、用于在所述检查体积内在不同空间方向上生成切换的磁场梯度的多个梯度线圈(4、5、6)、具有不同空间灵敏度分布的至少两个RF线圈(11、12、13)、用于控制RF脉冲和切换的磁场梯度的时间次序的控制单元(15)、以及重建单元(17)，其中，所述MR设备(1)被布置为执行以下步骤：

-使目标(10)经受成像序列，所述成像序列包括用于同时激励两个或更多个空间分离的图像切片的相位调制多切片RF脉冲，其中，所述RF脉冲的相位调制方案是从所述至少两个RF线圈(11、12、13)的所述空间灵敏度分布导出的，

-采集MR信号，以及

-根据采集到的MR信号来重建针对每个图像切片的MR图像，

其中，来自不同图像切片的MR信号贡献是基于所述至少两个RF线圈(11、12、13)的所述空间灵敏度分布并且基于所述RF脉冲的所述相位调制方案而被分离的。

9.一种在MR设备上运行的计算机程序，所述计算机程序包括指令，所述指令用于：

-生成成像序列，所述成像序列包括用于同时激励两个或更多个空间分离的图像切片的相位调制多切片RF脉冲，其中，所述RF脉冲的相位调制方案是从至少两个RF线圈(11、12、13)的空间灵敏度分布导出的，

-采集MR信号，以及

-根据采集到的MR信号来重建针对每个图像切片的MR图像，其中，来自不同图像切片的MR信号贡献是基于所述至少两个RF线圈(11、12、13)的所述空间灵敏度分布并且基于所述RF脉冲的所述相位调制方案而被分离的。