并行多切片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等人(journalofmagneticresonanceimaging,13,313-317,2001)提出了在多切片成像的情况下也在切片方向上应用灵敏度编码以提高扫描效率。breuer等人(magneticresonanceinmedicine,53,684-691,2005)提出的一种称为“controlledaliasinginparallelimagingresultsinhigheracceleration”(caipirinha)的方法改进了此基本想法。该技术通过将相位偏移赋予每个图像切片的mr信号，从而在多切片采集过程中修改了每个个体切片中的混叠伪影的出现，从而改善了后续的并行图像重建过程。因此，caipirinha是一种并行多切片成像技术，与仅使用纯后处理方法的其他多层并行成像概念相比，其效率更高。在caipirinha中，使用相位调制的多切片rf脉冲(类似于已知的hadamard脉冲)同时激发任意厚度和距离的多个切片。采集的mr信号数据是欠采样的，生成叠加的切片图像，这些切片图像表现为彼此相对偏移，对应于不同的相位偏移。根据傅立叶位移定理，通过由rf脉冲的相位调制方案施加的相位偏移来控制混叠切片图像的移位。通过使用此移位，可以改善反向重建问题的数值调节，所述调节将分离的片段的个体信号贡献分开。

us2014/0225612a1公开了一种使用分段回波平面成像(epi)脉冲序列进行mr成像的方法。脉冲序列包括多个多切片rf脉冲，用于同时激发两个或多个空间上分离的图像切片。沿着条带编码方向应用特定的梯度编码方案，以将受控的相位偏移赋予不同的条带。使用并行成像重建方法将采集的mr信号数据重建为mr图像，所述方法将成像数据中的交叠切片分开，以针对跨成像对象的每个切片提供一系列mr图像。

us2016/0018499a1解决了caipirinha使用rf脉冲的固定相位调制方案的问题，使得相邻切片的相对位移为例如视场(fov)尺寸的一半或fov尺寸的其他整数分数。该固定方案的缺点是它不考虑先验信息。因此，没有充分考虑接收rf线圈的阵列的编码能力以及成像问题的基本结构，这可能导致次优的相位调制并因此导致次优的重建性能。作为此问题的解决方案，us2016/0018499a1建议使用可用的线圈灵敏度信息来导出经调整的切片特定的相位偏移，以进一步优化编码过程，并且因此优化展开问题，从而改善最终图像质量。

技术实现要素：

根据本发明，公开了一种对放置在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等人在magneticresonanceinmedicine2013,第69卷,第793-802页)，不同之处在于，不仅考虑单个mr图像，而且考虑了多个图像切片，并且通过将相应的相位信息合并到编码矩阵中，从而将概念扩展到了多切片成像。

本发明的方法的成像序列是多击发序列，其需要以分段方式对k空间进行采样。换句话说，在成像序列的每次击发中仅对k空间的子集进行采样。值得注意的是，其中，采集在由变化的相位偏移赋予的相位编码方向上被欠采样的k空间数据。

在一个可能的实施例中，可变相位偏移可以借助于rf脉冲的相位调制来赋予(像，例如，在已知的pomp技术中，参见glover等人，journalofmagneticresonanceimaging，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来执行的。在每个相位编码梯度blip24的应用过程中，还会在切片编码方向上施加另一个磁场梯度blip25。根据本发明，使用切片编码梯度blip25来向每个图像切片的mr信号赋予相位偏移。尽管可以使用梯度blip25来施加不同的相位偏移，但是也可以通过连续击发的切片选择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方程组设置为：

其中矢量包括需要根据经由单独的rf线圈(“通道”)测得的混叠像素值计算出的图像像素值。混叠像素值由矢量描述。s是线圈灵敏度矩阵，其由所使用的rf线圈的空间灵敏度概况确定。灵敏度矩阵s的大小取决于击发的欠采样、使用的rf线圈数以及由多切片rf脉冲激发的切片数。在根据本发明的成像序列的每一次击发中施加的相位偏移可以通过矩阵φsh来考虑，其是描述在切片方向编码所得到的额外相位的对角矩阵。sense方程可在此基础上写为：

这里，φsh是对角矩阵，其包含每个位置的相位编码：

不仅必须考虑在切片方向(kz)上的不同相位编码，而且还必须考虑在相位编码方向(ky)上的不同相位编码。因此，描述了通过在常规相位编码方向(ky)上的相位编码和在切片方向(kz)上的相位偏移两者赋予给各个像素值的相位。线圈灵敏度编码和相位编码(y和z)可以组合为一个编码矩阵esh：

sφsh＝esh

最后，可以将所有n次击发的公式组合在一个广义的sense重建内核中：

其中，包括最终切片图像的所有像素值，这些最终切片图像是由包括切片方向上的相位编码的多层切片sense重建得到的。的最小二乘解(无噪声去相关和正则化)是：

本发明的多层采集方法可以看作是三维扫描，其中，通过相应的相位编码对ky-kz空间进行采样，其中，kz线的数量等于同时激发的切片的数量。利用本发明的多击发方法，变得有可能以不同的kz编码多次采集给定的ky编码步骤。这使得可以最佳地对三维k空间进行采样(使用欠采样)。图3示出了根据本发明的ky-kz采样的示例，其中，八个切片同时被激发。对epi序列的四个击发a，b，c和d施加了沿切片方向施加的四个不同相位偏移。欠采样是在第ky和kz方向两者上施加，在方向具有交替的采样。在所描绘的实施例中，每个击发应用恒定的kz编码，使得在单次击发的epi序列内不需要在切片方向上施加额外的梯度blip。