并行磁共振的图像重建方法、装置及并行磁共振成像系统的制作方法
【专利摘要】本发明公开了一种并行磁共振的图像重建方法，属于磁共振成像技术领域。该方法将用于单线圈MRI重建的二维k空间数据的相位条件扩展到适用于并行成像重建的三维k空间数据集；通过构造基于多线圈k空间数据的低秩数据矩阵，并利用矩阵填充的方法来重建这个数据矩阵从而实现并行磁共振图像的重建。从而可用较少的采样数据和更短的图像重建时间获得更好的重建图像质量，且重建过程中不需要自校正环节，可以更好地抑制噪声。本发明还公开了一种并行磁共振的图像重建装置及一种并行磁共振成像系统。
【专利说明】
并行磁共振的图像重建方法、装置及并行磁共振成像系统
技术领域
[0001 ] 本发明设及磁共振成像(Magnetic Resonance Imaging,简称MRI)技术领域，尤其 设及一种并行磁共振的图像重建方法、装置及并行磁共振成像系统。
【背景技术】
[0002] 磁共振成像已经被广泛应用于医学诊断之中。它也在后续的单一放射疗法或多种 放射疗法结合的治疗中起到了重要的作用。减少扫描时间是MRI检查的一个核屯、问题。对于 图像重建而言，减少扫描时间通常是通过减少采样点的数量来实现的。总体而言，受约束的 重建和并行成像是两个实现运一目的的两种有效方法。
[0003] 由于在减少数据获取和伪影方面的优越性能，受约束的重建在现代MRI中受到广 泛的重视。受约束的方法运用数学工具和先验知识在重建过程中对未采样的实验数据进行 估计和补偿。在运些方法中，基于先验知识的相关关系被用来描述可W减少空间采样量的 数据相关性。相关关系表明数据之间存在着冗余信息。在已知数据的基础上，利用运些冗余 信息可W推测出未采样的数据。相位约束重建是一种常用的受约束的重建方法。它利用了 傅里叶变换的对称性。
[0004] 并行成像是另外一种用于缩短MRI扫描时间的常用手段，它利用多个空间分布的 线圈同时接受信号。多线圈接受的信号提供了数据冗余，运些数据冗余能够更好地重建未 采样数据。常用的并行成像方法包括：（1)基于明确的敏感度测量的重建方法，例如SMASH和 SENSE;(2)基于自校正信号(ACSs)获取的重建方法，例如GRAPPA和SPIRiT;(3)无需自校正 信号(ACSs)的重建方法，例如SAKE。
[0005] 在一个MRI重建方法中将相位约束和并行成像相结合会有很多益处。两个最主要 的优点是更少的数据采集量和更好的重建结果。在现有的此类方法中，相位约束条件和并 行成像的结合是简单地遵循一定顺序的，即先进行并行成像再进行相位约束重建或者先进 行相位约束重建再进行并行成像，且往往需要自校正环节，图像重建过程比较繁琐，其所需 的采样数据及计算量均较大，且成像质量也需进一步提高。
【发明内容】
[0006] 本发明所要解决的技术问题在于克服现有技术不足，提供一种并行磁共振的图像 重建方法，用较少的采样数据和更短的图像重建时间即可获得更好的重建图像质量。
[0007] 本发明具体采用W下技术方案解决上述技术问题：
[000引一种并行磁共振的图像重建方法，首先利用多个并行磁共振线圈的欠采样数据重 建完整的=维k空间数据集，然后根据所述完整的=维k空间数据集生成二维磁共振图像； 所述完整的=维k空间数据集的重建方法具体如下：
[0009] 步骤1、利用W下方法将=维k空间数据集映射为一个低秩的二维数据矩阵D:
[0010] 假设并行磁共振k空间数据集尺寸为LXMXN，其中N为并行磁共振线圈数目，每个 并行磁共振线圈的k空间数据集的尺寸为LXM;用A[m，v|表示第i个并行磁共振线圈上k空 间[U，v]处的采样点;则所述二维数据矩阵D的构造如下：
[00" ] D = [01，02，…，01，...，D(l-2R)X(M-2R) ]，1 二U-R+ (V-R-I) X (L-2 X R)
[0012]其中，R为预设的k空间中W[u，v]位置为中屯、的邻域半径;〇1是D中第1列的列向 量，它的形式如下：
[0013]
[0014] 分和并行重建部分，具体结构如下：
[0015]
[0016]
[0017]其中，上标*表示复数的共辆;片':表示吗中的第t个元素，它的定义如下：
[001 引
[0019]其中，[11可，乂-9]，[化皆]￡罕圳中圳表示^[11，乂]位置为中屯、，1^为半径的邻域 中的每一个数据点；
[0020] 步骤2、W所述欠采样数据作为所构造二维数据矩阵D中的已知元素，对二维数据 矩阵D进行填充；
[0021] 步骤3、将填充完毕的二维数据矩阵D恢复为完整的=维k空间数据集。
[0022] 根据相同的发明思路还可W得到W下技术方案：
[0023] -种并行磁共振的图像重建装置，包括用于利用多个并行磁共振线圈的欠采样数 据重建完整的=维k空间数据集的k空间数据集重建单元，W及用于根据k空间数据集重建 单元所重建的完整的=维k空间数据集生成二维磁共振图像的图像生成单元;所述k空间数 据集重建单元包括：
[0024] 低秩矩阵构造模块，其功能是利用W下方法将=维k空间数据集映射为一个低秩 的二维数据矩阵D:
[0025] 假设并行磁共振k空间数据集尺寸为LXMXN，其中N为并行磁共振线圈数目，每个 并行磁共振线圈的k空间数据集的尺寸为LXM;用表示第i个并行磁共振线圈上k空 间[u，v]处的采样点;则所述二维数据矩阵D的构造如下：
[0026] D = [01，，…，01，...，D(l-2R)X(M-2R) ]，1 二U-R+ (V-R-I) X a-2 X R)
[0027] 其中，R为预设的k空间中W[u，v]位置为中屯、的邻域半径;〇1是D中第1列的列向 量，它的形式如下：
[0028
[0029 分和并行重建部分，具体结构如下：
[0030
[0031
[0032] 其中，上标*表示复数的共辆;片'表示吗中的第t个元素，它的定义如下：
[0033]
[0034] 其中，[u-p，v-q]，[p，y]E卒化外[-化刊表示W[u，v]位置为中屯、，Wr为半径的邻域 中的每一个数据点；
[0035] 矩阵填充模块，其功能是W所述欠采样数据作为低秩矩阵构造模块所构造二维数 据矩阵D中的已知元素，对二维数据矩阵D进行填充；
[0036] k空间数据恢复模块，其功能是将矩阵填充模块所填充完毕的二维数据矩阵D恢复 为完整的=维k空间数据集。
[0037] -种并行磁共振成像系统，包括用于利用多个并行磁共振线圈对检查对象进行磁 共振数据采样的采样单元，W及如上所述并行磁共振的图像重建装置。
[0038] 相比现有技术，本发明具有W下有益效果：
[0039] 1、相比现有技术，本发明所需要的采样数量更少；
[0040] 2、本发明通过低秩矩阵填充，可一次性同时完成相位约束部分和并行部分的成 像，因此可用更短的重建时间获得更好的磁共振图像重建效果；
[0041] 3、本发明图像重建过程不需要自校正环节，可W更好地抑制噪声。
【附图说明】
[0042] 图1为仿真数据集示例；
[0043] 图2为仿真数据集重建示例和重建误差；
[0044] 图3为真实数据集示例；
[0045] 图4为真实数据集添加噪声重建示例和重建误差。
【具体实施方式】
[0046] 下面结合附图对本发明的技术方案进行详细说明：
[0047] 本发明的思路是将用于单线圈MRI重建的二维k空间数据的相位条件扩展到适用 于并行成像重建的=维k空间数据集;通过构造基于多线圈k空间数据的低秩数据矩阵，并 利用矩阵填充的方法来重建运个数据矩阵从而实现并行磁共振图像的重建。从而可用较少 的采样数据和更短的图像重建时间获得更好的重建图像质量，且重建过程中不需要自校正 环节，可W更好地抑制噪声。
[0048] 为了便于公众理解，下面W-个优选实施例来对本发明技术方案进行详细说明。
[0049] 本实施例中的磁共振图像重建过程具体包括W下步骤：
[0050] 步骤1、将=维k空间数据集构造成二维数据矩阵：
[0051] 在构造数据矩阵时，k空间采样点和该点k空间对称位置处各个线圈相同邻域构成 的数据块满足如下线性相关关系：
[00 对
（1)
[0053]式（1)表示并行MRI的k空间中的一个采样点可W由各个线圈上W它的k空间对称 位置为中屯、的一个相同大小的区域内全部采样点的线性组合来表示。运一线性组合关系取 决于并行MRI系统各个线圈的敏感度W及各个线圈上磁共振图像在截止邻域内的相位平滑 程度。依据运一关系可W进一步扩展构造一个低秩矩阵D。
[0054]具体地，假设并行磁共振k空间数据集尺寸为LXMXN，其中N为并行磁共振线圈数 目，每个并行磁共振线圈的k空间数据集的尺寸为LXM;用A[w，v|表示第i个并行磁共振线 圈上k空间[u，v]处的采样点;则所述二维数据矩阵D的构造如下：
[005引 D = [01，，…，01，...，D(l-2R)X(M-2R) ]，1 二U-R+ (V-R-I) X a-2 X R) (2)
[0056]其中，R为预设的k空间中W[u，v]位置为中屯、的邻域半径;〇1是D中第1列的列向 量，它的形式如下：
[0化7]
[005引 蛋构如下：
[0化9]
[0060]
[0061] ：
[0062] (6)
[0063]
[0064] 位置为中屯、，WR为半径的邻域中的每一个数据点可W证明，数据矩阵D有块 化nkel结构，局部副对角线元素相同。而有块化nkel结构的矩阵有低秩的性质。因此，并行 磁共振=维k空间数据集的重建问题就转化成了一个具有结构特征的低秩矩阵的填充问 题。
[0065] 步骤2、W所述欠采样数据作为所构造二维数据矩阵D中的已知元素，对二维数据 矩阵D进行填充：
[0066] 矩阵填充是继压缩感知之后的又一引人注目的新研究领域，它应用在很多科学领 域和工程领域，如协同过滤、图像修复、机器学习、控制、计算机视觉和在传感网络工程中预 测丢失的数据等等，运个问题也就是从一些观测到的矩阵元素来填充运个矩阵。通过对许 多实际情况的总结，研究人员建立了矩阵填充的理论模型:在适当的条件下(样本数目的控 制和矩阵相干性），可W在较高的概率下通过求解一个凸优化问题(最小核范数问题或核范 数正则线性最小二乘问题)来精确或近似精确的恢复低秩或近似低秩矩阵。
[0067] 对二维数据矩阵D进行填充可W采用拉格朗日乘子法、h范数法、SVT法等，本发明 优选采用SVT法，下面对该方法进行简要介绍。
[0068] 数据矩阵D中全部的已知元素构成一个集看运样， 问题可W写成如下公式：
minimize rank (X) subject !;〇 JtT 托 _y)二 Z)(x, _y), (X, ) e Q
[0070] 运里X是D的估计。
[0071] 然而，最小化秩是一个非确定性多项式时间问题。当预先确定数据矩阵D的秩时， 可W用一个凸的严格近似方法来简化运一复杂的NP难问题。该方法提出公式(8)来重新解 释运一凸优化问题
[007^
佩
[0073] 运里Mx II*表示矩阵X的核范数，它定义为
[0074]
(9)
[0075] 公式(8)中的问题可W用化dzow算法来解决。Cadzow算法是一个迭代算法，每一次 循环包含了=个步骤：
[0076] 第一步:用奇异值分解的方法计算当前要估计的数据矩阵Xn的奇异值矩阵En
[0077] Xn = UnEnVn (10)
[0078] 第二步:取En的r个最大值，其余置零来构造新的奇异值矩阵En+1。新估计的数据 矩阵Xn+1计算如下：
[0079] Xwl = UnEwlVn (10)
[0080] 第=步:对新估计的Xn+1要保持数据的一致性，表示如下：
[0081] Xn+I(x，y) =D(x，y)，（x，y) G Q (11)
[0082] 如果达到预先设定的迭代次数或者两次连续的估计之间的误差小于一个阔值，那 么迭代过程终止。
[0083] 步骤3、将填充完毕的二维数据矩阵D恢复为完整的=维k空间数据集：
[0084] =维k空间数据集的恢复实际上是步骤1中映射的逆过程。本发明优选采用W下方 法进行=维k空间数据集的恢复：
[0085] 假设第i个通道k空间[u，v]处的数据点为Ki[u，v]，则可根据W下重建方程恢复Ki [U，v]:
[0086]
[0087] 式中，片(0是巧的第t项;A是相位约束条件的权重，表示相位约束条件在整个重建 中所占比例。
[0088] 步骤4、根据完整的=维k空间数据集生成二维磁共振图像：
[0089] 得到完整的S维k空间数据集后，即可利用SENSE、GRAPPA等方法生成二维磁共振 图像，本发明优选采用ESPIRiT方法。
[0090] 为了验证本发明的效果，将其与现有的SAKE方法(简称SAKE)进行仿真数据和真实 数据的对比实验。实验一使用的仿真数据集为:利用Guerquin-Kern的方法将一个相位平滑 的化epp-Logan幻影的二维k空间数据集扩展成S维k空间数据集，它的完全采集的k空间数 据的重建图像如图1所示。实验二使用的真实数据集为:Tl脑部数据集，它的完全采集的k空 间数据的重建图像如图3所示。
[0091] 对于同样的实验对象，本发明方法(简称PCLR)的数据使用量分别为:100%、80%。
[0092] 1)不加噪声，100次迭代重建，5巧截止邻域，采样率分别为40%，重建结果和重建 误差分别对应图2中第一行、第二行所示，其中第一列为采用SAKE算法得到的重建结果，第 二列为本发明方法使用全部采样点的重建结果，第=列为本发明方法使用80%采样点的重 建结果。
[00W] 2伽高斯噪声，100次迭代重建，5巧截止邻域，40 %采样率，信噪比分别为30地，重 建结果和重建误差分别对应图4中第一行、第二行所示，其中第一列为采用SAKE算法得到的 重建结果，第二列为本发明方法使用全部采样点的重建结果，第=列为本发明方法使用 80%采样点的重建结果。
[0094]根据图2和图4的对比实验结果可见本发明方法相比于现有的SAKE算法具有如下 优势：
[00%] 1)可W正确地得到重建后的磁共振图像，更好地保持磁共振图像的结构特征；
[0096] 2)可W需要更少的采样点
[0097] 3)抑制高斯噪声的影响；
[0098] 4)更快速地实现重建。
【主权项】
1. 一种并行磁共振的图像重建方法，首先利用多个并行磁共振线圈的欠采样数据重建 完整的三维k空间数据集，然后根据所述完整的三维k空间数据集生成二维磁共振图像;其 特征在于，所述完整的三维k空间数据集的重建方法具体如下：步骤1、利用以下方法将三维 k空间数据集映射为一个低秩的二维数据矩阵D: 假设并行磁共振k空间数据集尺寸为LXMXN，其中N为并行磁共振线圈数目，每个并行 磁共振线圈的k空间数据集的尺寸为LXM;用表示第i个并行磁共振线圈上k空间[u， v]处的采样点;则所述二维数据矩阵D的构造如下： D=[D1,D2,. . . ,?1,. . . ,D(l_2R)x(m_2R)],1=u-R+(v-R-1)X(L-2XR) 其中，R为预设的k空间中以[u，v]位置为中心的邻域半径;D1是D中第1列的列向量，它的 形式如下：其中，<和W分别表示k空间相位约束重建部分和并行重建部分，具体结构如下：其中，上标*表示复数的共辄;纪表示碑中的第t个元素，它的定义如下：其中，[u-p，v-q]，[/λ?/] e 4表示以[u，v]位置为中心，以R为半径的邻域中的每 一个数据点； 步骤2、以所述欠采样数据作为所构造二维数据矩阵D中的已知元素，对二维数据矩阵D 进行填充； 步骤3、将填充完毕的二维数据矩阵D恢复为完整的三维k空间数据集。2. 如权利要求1所述并行磁共振的图像重建方法，其特征在于，利用SVT方法对二维数 据矩阵D进行填充。3. 如权利要求1所述并行磁共振的图像重建方法，其特征在于，所述将填充完毕的二维 数据矩阵D恢复为完整的三维k空间数据集的具体方法如下： 根据以下重建方程恢复第i个通道k空间[u，v]处的数据点式中，W⑴是旮的第t项;λ是相位约束条件在重建中的权重。4. 如权利要求1所述并行磁共振的图像重建方法，其特征在于，根据所述完整的三维k 空间数据集，使用ESPIRiT方法生成二维磁共振图像。5. -种并行磁共振的图像重建装置，包括用于利用多个并行磁共振线圈的欠采样数据 重建完整的三维k空间数据集的k空间数据集重建单元，以及用于根据k空间数据集重建单 元所重建的完整的三维k空间数据集生成二维磁共振图像的图像生成单元;其特征在于，所 述k空间数据集重建单元包括： 低秩矩阵构造模块，其功能是利用以下方法将三维k空间数据集映射为一个低秩的二 维数据矩阵D: 假设并行磁共振k空间数据集尺寸为LXMXN，其中N为并行磁共振线圈数目，每个并行 磁共振线圈的k空间数据集的尺寸为LXM;用/3, hv]表示第i个并行磁共振线圈上k空间[u， v]处的采样点;则所述二维数据矩阵D的构造如下： D=[D1,D2,. . . ,?1,. . . ,D(l_2R)x(m_2R)],1=u-R+(v-R-1)X(L-2XR) 其中，R为预设的k空间中以[u，v]位置为中心的邻域半径;D1是D中第1列的列向量，它的 形式如下：其中，名和软分别表示k空间相位约束重建部分和并行重建部分，具体结构如下：其中，上标*表示复数的共辄;允表示对中的第t个元素，它的定义如下：其中，[u-p，v_q]，[?，以表示以[u，v]位置为中心，以R为半径的邻域中的每 一个数据点； 矩阵填充模块，其功能是以所述欠采样数据作为低秩矩阵构造模块所构造二维数据矩 阵D中的已知元素，对二维数据矩阵D进行填充； k空间数据恢复模块，其功能是将矩阵填充模块所填充完毕的二维数据矩阵D恢复为完 整的三维k空间数据集。6. 如权利要求5所述并行磁共振的图像重建装置，其特征在于，所述矩阵填充模块利用 SVT方法对二维数据矩阵D进行填充。7. 如权利要求5所述并行磁共振的图像重建装置，其特征在于，所述k空间数据恢复模 块将填充完毕的二维数据矩阵D恢复为完整的三维k空间数据集的具体方法如下： 根据以下重建方程恢复第i个通道k空间[u，v]处的数据点式中，故(〇是旮的第t项;λ是相位约束条件在重建中的权重。8. 如权利要求5所述并行磁共振的图像重建装置，其特征在于，所述图像生成单元使用 ESPIRiT方法生成二维磁共振图像。9. 一种并行磁共振成像系统，包括用于利用多个并行磁共振线圈对检查对象进行磁共 振数据采样的采样单元，以及如权利要求5~8任一项所述并行磁共振的图像重建装置。
【文档编号】G06T5/00GK105957117SQ201610264588
【公开日】2016年9月21日
【申请日】2016年4月26日
【发明人】贺润国, 姜龙玉, 宋文博, 洪亚萍, 伍家松, 舒华忠
【申请人】东南大学