快速多对比度磁共振扩散成像和重建方法及系统

1.本发明涉及磁共振扩散成像领域，具体地，涉及一种快速多对比度磁共振扩散成像和重建方法及系统。


背景技术：

2.人体组织内的水分子的扩散能力的变化可以反映出生理组织的构成以及病变情况。磁共振扩散成像(diffusion mri,dmri)是一项能够无创地检查活体组织扩散能力的成像技术。被广泛地应用在微血管灌注和扩散效应活体定量研究中，在急性脑卒中、肝脏疾病、癌症等疾病的早期诊断以及神经科学研究中具有重要的作用。扩散效应的定量图像的计算依赖于多个具有不同b值，方向，扩散时间等参数的磁共振扩散加权图像(diffusion weighted imaging,dwi)。由于需要获得多幅不同对比度的图像，磁共振扩散成像需要较长的扫描时间。与此同时，因为扩散梯度的加入，信号衰减速度快，采集的分辨率因此受限，同时图像的畸变大。
3.现有的磁共振扩散成像主要有两类方法。
4.第一类是单次激发平面回波序列。这类方法在一次激发之后就完成整个k空间的数据采集，扫描时间短，对运动鲁棒。但是因为较长的回波链，导致分辨率低，信噪比低，图像畸变大。
5.第二类是多次激发平面回波序列。这类方法采用多次激发策略，每次激发采集k空间的一部分数据，多次激发获取的数据合并为完整的k空间，可以实现更高的空间分辨率和图像质量。然而，扩散编码梯度导致不同激发之间存在着相位的变化，此外多次激发间还可能存在着运动；直接合并这些数据将导致图像伪影。此外，对于adc、dti、dki、ivim等扩散定量成像应用，需要获得多幅多次激发的图像以实现高分辨率，扫描时间变长。
6.专利文献(申请号201210388038.1)公开了一种快速扩散磁共振成像和重建方法，包括以下步骤：(s1)通过多通道线圈在n个扩散加权方向上对被测目标进行信号采集；(s2)将所获得的互补的k空间数据进行合并，以获得满采样k空间数据kc；(s3)基于与所述k空间数据对应的不同扩散加权方向的图像和所述满采样k空间数据kc所对应的图像进行初步重建，以获得初步图像以及(s4)基于所述像初步图像进行正则化重建，并迭代直至收敛，以获得所需的最终图像(i1，
…
，in)。该发明使用普通而非keyhole回波平面成像序列进行数据采集，需要对采集的扩散图像进行k空间的合并，同时还需要求解线圈灵敏度参数，难以避免合并时带来的重建误差，也不能避免线圈灵敏度参数估计不准确给重建带来的不准确性。
7.专利文献cn103675737a(申请号：cn201310659202.2)公开了一种扩散磁共振成像和重建方法。该方法包括：s1、使用多个通道线圈，采用多次激发扩散成像方式对被测目标进行信号采集，以获取k空间数据；s2、计算线圈敏感度图，进行迭代初始化；s3、根据采集到的k空间数据、计算得到的线圈敏感度图和初始化参数，基于pocs算法，对所需的扩散图像进行迭代重建。该发明采用多次激发策略，需要先对采集的多线圈扩散图像进行k空间的合并，不能通过变化对比度在多对比度的扩散成像中减少成像时间。


技术实现要素：

8.针对现有技术中的缺陷，本发明提供一种快速多对比度磁共振扩散成像和重建方法及系统，能有效提升分辨率、降低形变伪影、提升扫描效率。
9.根据本发明提供的一种快速多对比度磁共振扩散成像和重建方法，包括：
10.步骤s1：通过多个接收线圈，使用keyhole平面回波序列欠采多幅可变对比度扩散加权图像的k空间数据；
11.步骤s2：根据采集结果，构建包括多线圈图像保真项和基于线圈、空间及对比度域张量低秩性的图像正则项的目标函数；
12.步骤s3：通过最优化算法最小化目标函数，重建所有接收线圈的扩散图像并根据对比度变化的具体方式获取相关定量参数图像。
13.优选地，在所述步骤s1中：
14.所述接收线圈数量为线圈通道数量大于等于1个；
15.所述可变对比度包括：在平面回波序列的不同tr间的对比度参数能够任意变化，包括但不限于b值、扩散编码方向、和扩散时间的任意变化；
16.针对特定一个对比度参数，所述任意变化为不同tr间的对比度参数保持不变或存在不同；
17.所述keyhole-epi序列使用以下变密度k空间采集：
18.在中心k空间区域以欠采样率rc密集采集；
19.在周围k空间区域以欠采样率r
p
均匀移位欠采，以使相邻的tr采集循环移位的相位编码行；
20.所述keyhole-epi序列的变密度k空间采集通过如下方式实施：
21.确定跨越单个相位编码行的梯度动量；并根据所计算结果，计算其他所需的跨越不同相位编码行的blip梯度动量；将所需的对应动量的blip梯度，根据预设采样的轨迹进行时序排列；根据以上计算结果，采用bipolar epi梯度对k空间数据进行采集。
22.优选地，在所述步骤s2中：
23.待求解图像为多线圈下的可变对比度扩散加权图像，而非多线圈合并后的扩散加权图像；
24.目标函数中不包括线圈灵敏度，求解该目标函数无需先求解线圈灵敏度或进行并行成像的自校正；
25.基于线圈、空间及对比度域张量低秩性的图像保真项利用多线圈可变对比度扩散加权图像的张量低秩性对待求解图像的变化范围做出约束，包括基于全局空间内张量低秩性的约束和基于局部空间张量低秩性的约束；
26.目标函数为：
[0027][0028]
其中，y表示预处理后得到多维k空间数据，d表示欠采样算子，f表示傅里叶变换算子，x表示待重建的多线圈多对比度图像，α为正则化系数，bix为以第i个像素为中心的线圈、空间及对比度三维张量，r为该三维张量基于tucker分解或parafac分解下的求秩算子，nv
为预先设定的值，如使用全局低秩约束，则nv＝1，如使用局部低秩约束，则nv可等于像素数量。
[0029]
优选地，在所述步骤s3中：
[0030]
其目标函数能够通过任意一种最优化算法求解，包括梯度下降法、共轭梯度法、拟牛顿法、交替乘子法、近端梯度法、投影梯度法、坐标下降法，以重建所有接收线圈的扩散图像并根据具体应用场景获取相关定量参数图像。
[0031]
优选地，该成像方法能够用于扩散成像应用场景包括：多次激发下的扩散加权成像、表观扩散系数定量成像、扩散张量成像、扩散峰度成像、体素内不相干运动成像。
[0032]
根据本发明提供的一种快速多对比度磁共振扩散成像和重建系统，包括：
[0033]
模块m1：通过多个接收线圈，使用keyhole平面回波序列欠采多幅可变对比度扩散加权图像的k空间数据；
[0034]
模块m2：根据采集结果，构建包括多线圈图像保真项和基于线圈、空间及对比度域张量低秩性的图像正则项的目标函数；
[0035]
模块m3：通过最优化算法最小化目标函数，重建所有接收线圈的扩散图像并根据对比度变化的具体方式获取相关定量参数图像。
[0036]
优选地，在所述模块m1中：
[0037]
所述接收线圈数量为线圈通道数量大于等于1个；
[0038]
所述可变对比度包括：在平面回波序列的不同tr间的对比度参数能够任意变化，包括但不限于b值、扩散编码方向、和扩散时间的任意变化；
[0039]
针对特定一个对比度参数，所述任意变化为不同tr间的对比度参数保持不变或存在不同；
[0040]
所述keyhole-epi序列使用以下变密度k空间采集：
[0041]
在中心k空间区域以欠采样率rc密集采集；
[0042]
在周围k空间区域以欠采样率r
p
均匀移位欠采，以使相邻的tr采集循环移位的相位编码行；
[0043]
所述keyhole-epi序列的变密度k空间采集通过如下方式实施：
[0044]
确定跨越单个相位编码行的梯度动量；并根据所计算结果，计算其他所需的跨越不同相位编码行的blip梯度动量；将所需的对应动量的blip梯度，根据预设采样的轨迹进行时序排列；根据以上计算结果，采用bipolar epi梯度对k空间数据进行采集。
[0045]
优选地，在所述模块m2中：
[0046]
待求解图像为多线圈下的可变对比度扩散加权图像，而非多线圈合并后的扩散加权图像；
[0047]
目标函数中不包括线圈灵敏度，求解该目标函数无需先求解线圈灵敏度或进行并行成像的自校正；
[0048]
基于线圈、空间及对比度域张量低秩性的图像保真项利用多线圈可变对比度扩散加权图像的张量低秩性对待求解图像的变化范围做出约束，包括基于全局空间内张量低秩性的约束和基于局部空间张量低秩性的约束。
[0049]
目标函数为：
[0050][0051]
其中，y表示预处理后得到多维k空间数据，d表示欠采样算子，f表示傅里叶变换算子，x表示待重建的多线圈多对比度图像，α为正则化系数，bix为以第i个像素为中心的线圈、空间及对比度三维张量，r为该三维张量基于tucker分解或parafac分解下的求秩算子，nv为预先设定的值，如使用全局低秩约束，则nv＝1，如使用局部低秩约束，则nv可等于像素数量。
[0052]
优选地，在所述模块m3中：
[0053]
其目标函数能够通过任意一种最优化算法求解，包括梯度下降法、共轭梯度法、拟牛顿法、交替乘子法、近端梯度法、投影梯度法、坐标下降法，以重建所有接收线圈的扩散图像并根据具体应用场景获取相关定量参数图像。
[0054]
优选地，该成像方法能够用于扩散成像应用场景包括：多次激发下的扩散加权成像、表观扩散系数定量成像、扩散张量成像、扩散峰度成像、体素内不相干运动成像。
[0055]
与现有技术相比，本发明具有如下的有益效果：
[0056]
1、本发明相对于单次激发平面回波成像类技术，通过keyhole平面回波序列成像方法减少了epi回波链的采集时间，可以实现更高的分辨率，更高的信噪比，更小的图像畸变；
[0057]
2、本发明相对于多次激发平面回波成像类技术，通过在多次激发过程中引入对比度变化，在较短的成像时间内得到多张不同对比度扩散加权图像以及相应的定量图像，通过使用线圈-空间-对比度域的低秩张量约束，避免了求解线圈灵敏度图的步骤，减小了变密度epi成像所导致的图像伪影，提高了重建结果的鲁棒性和精确性。
附图说明
[0058]
通过阅读参照以下附图对非限制性实施例所作的详细描述，本发明的其它特征、目的和优点将会变得更明显：
[0059]
图1为本技术的一个实施例的扩散磁共振成像和重建方法的流程图；
[0060]
图2为本技术的一个实施例的keyhole平面回波序列示意图；
[0061]
图3为本技术的一个具体实施例的keyhole平面回波序列轨迹示意图；
[0062]
图4为本技术的一个具体实施例的扩散磁共振成像和重建方法中的迭代重建流程图；
[0063]
图5为本技术的一个具体实施例与现有技术的回顾性实验重建结果对比示意图；
[0064]
图6为本技术的一个具体实施例与现有技术的前瞻性实验重建结果对比示意图。
具体实施方式
[0065]
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明，但不以任何形式限制本发明。应当指出的是，对本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变化和改进。这些都属于本发明的保护范围。
[0066]
实施例1：
[0067]
根据本发明提供的一种快速多对比度磁共振扩散成像和重建方法，如图1-图6，包括：
[0068]
步骤s1：通过多个接收线圈，使用keyhole平面回波序列欠采多幅可变对比度扩散加权图像的k空间数据；
[0069]
具体地，在所述步骤s1中：
[0070]
所述接收线圈数量为线圈通道数量大于等于1个；
[0071]
所述可变对比度包括：在平面回波序列的不同tr间的对比度参数能够任意变化，包括但不限于b值、扩散编码方向、和扩散时间的任意变化；
[0072]
针对特定一个对比度参数，所述任意变化为不同tr间的对比度参数保持不变或存在不同；
[0073]
所述keyhole-epi序列使用以下变密度k空间采集：
[0074]
在中心k空间区域以欠采样率rc密集采集；
[0075]
在周围k空间区域以欠采样率r
p
均匀移位欠采，以使相邻的tr采集循环移位的相位编码行；
[0076]
所述keyhole-epi序列的变密度k空间采集通过如下方式实施：
[0077]
确定跨越单个相位编码行的梯度动量；并根据所计算结果，计算其他所需的跨越不同相位编码行的blip梯度动量；将所需的对应动量的blip梯度，根据预设采样的轨迹进行时序排列；根据以上计算结果，采用bipolar epi梯度对k空间数据进行采集。
[0078]
步骤s2：根据采集结果，构建包括多线圈图像保真项和基于线圈、空间及对比度域张量低秩性的图像正则项的目标函数；
[0079]
具体地，在所述步骤s2中：
[0080]
待求解图像为多线圈下的可变对比度扩散加权图像，而非多线圈合并后的扩散加权图像；
[0081]
目标函数中不包括线圈灵敏度，求解该目标函数无需先求解线圈灵敏度或进行并行成像的自校正；
[0082]
基于线圈、空间及对比度域张量低秩性的图像保真项利用多线圈可变对比度扩散加权图像的张量低秩性对待求解图像的变化范围做出约束，包括基于全局空间内张量低秩性的约束和基于局部空间张量低秩性的约束；
[0083]
目标函数为：
[0084][0085]
其中，y表示预处理后得到多维k空间数据，d表示欠采样算子，f表示傅里叶变换算子，x表示待重建的多线圈多对比度图像，α为正则化系数，bix为以第i个像素为中心的线圈、空间及对比度三维张量，r为该三维张量基于tucker分解或parafac分解下的求秩算子，nv为预先设定的值，如使用全局低秩约束，则nv＝1，如使用局部低秩约束，则nv可等于像素数量。
[0086]
步骤s3：通过最优化算法最小化目标函数，重建所有接收线圈的扩散图像并根据
对比度变化的具体方式获取相关定量参数图像。
[0087]
具体地，在所述步骤s3中：
[0088]
其目标函数能够通过任意一种最优化算法求解，包括梯度下降法、共轭梯度法、拟牛顿法、交替乘子法、近端梯度法、投影梯度法、坐标下降法，以重建所有接收线圈的扩散图像并根据具体应用场景获取相关定量参数图像。
[0089]
具体地，该成像方法能够用于扩散成像应用场景包括：多次激发下的扩散加权成像、表观扩散系数定量成像、扩散张量成像、扩散峰度成像、体素内不相干运动成像。
[0090]
根据本发明提供的一种快速多对比度磁共振扩散成像和重建系统，包括：
[0091]
模块m1：通过多个接收线圈，使用keyhole平面回波序列欠采多幅可变对比度扩散加权图像的k空间数据；
[0092]
具体地，在所述模块m1中：
[0093]
所述接收线圈数量为线圈通道数量大于等于1个；
[0094]
所述可变对比度包括：在平面回波序列的不同tr间的对比度参数能够任意变化，包括但不限于b值、扩散编码方向、和扩散时间的任意变化；
[0095]
针对特定一个对比度参数，所述任意变化为不同tr间的对比度参数保持不变或存在不同；
[0096]
所述keyhole-epi序列使用以下变密度k空间采集：
[0097]
在中心k空间区域以欠采样率rc密集采集；
[0098]
在周围k空间区域以欠采样率r
p
均匀移位欠采，以使相邻的tr采集循环移位的相位编码行；
[0099]
所述keyhole-epi序列的变密度k空间采集通过如下方式实施：
[0100]
确定跨越单个相位编码行的梯度动量；并根据所计算结果，计算其他所需的跨越不同相位编码行的blip梯度动量；将所需的对应动量的blip梯度，根据预设采样的轨迹进行时序排列；根据以上计算结果，采用bipolar epi梯度对k空间数据进行采集。
[0101]
模块m2：根据采集结果，构建包括多线圈图像保真项和基于线圈、空间及对比度域张量低秩性的图像正则项的目标函数；
[0102]
具体地，在所述模块m2中：
[0103]
待求解图像为多线圈下的可变对比度扩散加权图像，而非多线圈合并后的扩散加权图像；
[0104]
目标函数中不包括线圈灵敏度，求解该目标函数无需先求解线圈灵敏度或进行并行成像的自校正；
[0105]
基于线圈、空间及对比度域张量低秩性的图像保真项利用多线圈可变对比度扩散加权图像的张量低秩性对待求解图像的变化范围做出约束，包括基于全局空间内张量低秩性的约束和基于局部空间张量低秩性的约束；
[0106]
目标函数为：
[0107][0108]
其中，y表示预处理后得到多维k空间数据，d表示欠采样算子，f表示傅里叶变换算
子，x表示待重建的多线圈多对比度图像，α为正则化系数，bix为以第i个像素为中心的线圈、空间及对比度三维张量，r为该三维张量基于tucker分解或parafac分解下的求秩算子，nv为预先设定的值，如使用全局低秩约束，则nv＝1，如使用局部低秩约束，则nv可等于像素数量。
[0109]
模块m3：通过最优化算法最小化目标函数，重建所有接收线圈的扩散图像并根据对比度变化的具体方式获取相关定量参数图像。
[0110]
具体地，在所述模块m3中：
[0111]
其目标函数能够通过任意一种最优化算法求解，包括梯度下降法、共轭梯度法、拟牛顿法、交替乘子法、近端梯度法、投影梯度法、坐标下降法，以重建所有接收线圈的扩散图像并根据具体应用场景获取相关定量参数图像。
[0112]
具体地，该成像方法能够用于扩散成像应用场景包括：多次激发下的扩散加权成像、表观扩散系数定量成像、扩散张量成像、扩散峰度成像、体素内不相干运动成像。
[0113]
实施例2：
[0114]
实施例2为实施例1的优选例，以更为具体地对本发明进行说明。
[0115]
根据本发明提供的一种快速多对比度磁共振扩散成像和重建方法，包括：
[0116]
步骤s1：通过多个接收线圈，使用keyhole平面回波(keyhole-epi)序列欠采多幅可变对比度扩散加权图像的k空间数据；
[0117]
步骤s2：构建包括多线圈图像保真项和基于线圈-空间-对比度域张量低秩性的图像正则项的目标函数；
[0118]
步骤s3：通过最优化算法最小化目标函数，重建所有接收线圈的扩散图像并根据对比度变化的具体方式获取相关定量参数图像。
[0119]
优选地，在所述步骤s1中：
[0120]
所述接收线圈数量为线圈数量大于等于1个；
[0121]
所述可变对比度包括：在平面回波序列的不同tr间的对比度参数能够任意变化，包括但不限于b值、扩散编码方向、和扩散时间的任意变化；
[0122]
针对特定一个对比度参数，所述任意变化为不同tr间的对比度参数保持不变或存在不同。
[0123]
优选地，所述keyhole-epi序列使用以下变密度k空间采集：
[0124]
在中心k空间区域以欠采样率rc密集采集；
[0125]
在周围k空间区域以欠采样率r
p
均匀移位欠采，以使相邻的tr采集循环移位的相位编码行；
[0126]
所述keyhole-epi序列的变密度k空间采集通过如下方式实施：
[0127]
确定跨越单个相位编码行的梯度动量；并根据所计算结果，计算其他所需的跨越不同相位编码行的blip梯度动量；将所需的对应动量的blip梯度，根据预设采样的轨迹进行时序排列；根据以上计算结果，采用bipolarepi梯度对k空间数据进行采集。
[0128]
优选地，在所述步骤s2中：
[0129]
待求解图像为多线圈下的可变对比度扩散加权图像，而非多线圈合并后的扩散加权图像。
[0130]
目标函数中不包括线圈灵敏度，求解该目标函数无需先求解线圈灵敏度或进行并
行成像的自校正。
[0131]
基于线圈-空间-对比度域张量低秩性的图像保真项利用多线圈可变对比度扩散加权图像的张量低秩性对待求解图像的变化范围做出约束，包括基于全局空间内张量低秩性的约束和基于局部空间张量低秩性的约束。
[0132]
具体地，根据本发明的一个实施例，使用局部张量低秩作为低秩约束时，构建的目标函数可以描述为：
[0133][0134]
其中，y表示预处理后得到多维k空间数据，d表示欠采样算子，f表示傅里叶变换算子，x表示待重建的多线圈多对比度图像，α为正则化系数，bix为以第i个像素为中心的线圈-空间-对比度三维张量，r为该三维张量基于tucker分解或parafac分解下的求秩算子，nv为预先设定的值，如使用全局低秩约束则nv＝1，如使用局部低秩约束则nv可等于像素数量。
[0135]
根据本发明的上述快速多对比度磁共振扩散成像和重建方法，可以避免求解线圈灵敏图参数同时避开直接合并具有不同相位的k空间，提高了重建的图像质量，允许了更高的分辨率和更高的采集效率。
[0136]
根据本发明的一个实施例，目标函数可以通过构造辅助变量，改写成如下述公式：
[0137][0138]
s.t bix＝mi[0139]
mi＝gi×
lai
×2bi×3ci，j＝1，2，
…
,nv[0140]
其中，y表示预处理后得到多维k空间数据，d表示欠采样算子，f表示傅里叶变换算子，x表示待重建的多线圈多对比度图像。α为正则化系数，bix为以第i个像素为中心的线圈-空间-对比度三维张量，r为该三维张量基于tucker分解或parafac分解下的求秩算子，nv为预先设定的值，如使用全局低秩约束则nv＝1，如使用局部低秩约束则nv可等于像素数量。mi为一个tensor类型的辅助变量，通过parafac或tucker分解得到，mi＝gi×1ai×2bi×3ci，i＝1，2，
…
，nv[0141]
根据本发明的一个实施例，目标函数可以改写为增广拉格朗日形式，如下述公式：
[0142][0143]
其中，y表示预处理后得到多维k空间数据，d表示欠采样算子，f表示傅里叶变换算子，x表示待重建的多线圈多对比度图像。α为正则化系数，bix为以第i个像素为中心的线圈-空间-对比度三维张量，r为该三维张量基于tucker分解或parafac分解下的求秩算子，nv为预先设定的值，如使用全局低秩约束则nv＝1，如使用局部低秩约束则nv可等于像素数量。re为取实数算子。mi为一个tensor类型的辅助变量，通过parafac或tucker分解得到，mi＝gi×1ai×2bi×3ci，i＝1，2，
…
，nv，μ为增广拉格朗日系数，λi为拉格朗日乘子。
[0144]
优选地，在所述步骤s3中：
[0145]
其目标函数可通过任意一种最优化算法求解，包括梯度下降法(gradient descent)，共轭梯度法(conjugategradients)，拟牛顿法(quasi-newton)，交替乘子法(alternatingdirectionmethod ofmultipliers)，近端梯度法(proximalgradient)，投影梯度法(projectedgradient)，坐标下降法(coordinatedescent)等算法，以重建所有接收线圈的扩散图像并根据具体应用场景获取相关定量参数图像。
[0146]
该成像方法可用于多次激发下的扩散加权成像、apparent diffusion coefficient(adc)定量成像、扩散张量成像(diffusion tensor imaging)、扩散峰度成像(diffusion kurtosis imaging)、体素内不相干运动成像(ivim)等多种扩散成像相关应用场景。
[0147]
所述根据应用场景获取相关定量参数图像表示：在需要获取定量参数图像时，参数图像和扩散图像为最终结果；在不需要获取定量参数图像时，重建的扩散图像即为最终结果。
[0148]
具体地，在本发明的一个实施例中，通过交替乘子法迭代最小化目标函数，交替乘子法的计算步骤如下：
[0149]
s31初始化：在正式迭代重建之前，需要初始化待求解扩散图像，即确定图像初始值，本发明的实施例中，初始值可以设为0，为方便描述，在公式中描述为，x选择初始点x0。
[0150]
s32分别求解子问题：在admm中，需要分别求解只有待求解扩散图像变量，辅助变量以及拉格朗日乘子变量的子问题，步骤s321-323分别描述了上文提到的对应的各个子问题的求解细节。
[0151]
步骤s321：求解xk，存在解析解
[0152][0153]
其中h表示共轭转置。
[0154]
步骤s322:求解子目标如下
[0155][0156]
优选地，首先进行张量的parafac或tucker分解
[0157][0158]
随后使用奇异值软阈值法：
[0159][0160][0161]
[0162]
步骤s323:求解
[0163][0164]
s33完成以上步骤后，自动地对每个通道的线圈图像进行更新，并进行循环迭代，达到下列条件之一则迭代停止，：(1)迭代次数大于50次；(2)梯度的无穷范数小于1*10-10
；(3)步长小于1*10-12
。判断是否收敛或者达到最大迭代次数。如果没有满足迭代停止条件，重新进入步骤s32继续迭代。如果满足迭代停止条件，迭代终止，得到磁共振扩散加权图像。
[0165]
s34根据采集的扩散加权参数类型，以及实际需求，求解adc、dti、dki、ivim等扩散定量成像，如果不需要求解定量图像，则磁共振扩散加权图像即为最终结果。
[0166]
实施例3：
[0167]
实施例3为实施例1的优选例，以更为具体地对本发明进行说明。
[0168]
本技术的目的在于提出一种快速多对比度磁共振扩散成像和重建方法。
[0169]
本技术提供的快速多对比度磁共振扩散成像和重建方法的总体步骤如图1所示，包括：步骤s1：通过多个接收线圈，使用keyhole平面回波(keyhole-epi)序列欠采多幅可变对比度扩散加权图像的k空间数据；步骤s2：构建包括多线圈图像保真项和基于线圈-空间-对比度域张量低秩性的图像正则项的目标函数；步骤s3：通过最优化算法最小化目标函数，重建所有接收线圈的扩散图像并根据对比度变化的具体方式获取相关定量参数图像。
[0170]
在下述的对本技术实施例的快速高分辨率扩散磁共振成像和重建方法的说明中，采用了线圈-空间-对比度域局部张量低秩约束(locally low rank tensor,llrt)结合数据保真项形成目标函数，并优选交替方向乘子法(alternating direction method of multipliers，admm)作为迭代重建算法，重建多对比度扩散加权图像并获得adc图像。
[0171]
下面结合图2，3，4对上述各步骤进行详细说明。
[0172]
(s101)在任意接收线圈数量下，通过keyhole平面回波序列采集k空间；图3为所述可变对比度keyhole平面回波序列的示意图。图3中，单位相位编码k空间的动量为mu，中心k空间blip梯度的动量为mu，周围k空间blip梯度的动量为r
pmu
,r
p
为周围k空间欠采率。相邻的tr在相位编码方向上，prephase梯度的动量相差rcmu/2。实现在中心k空间区域密集采集，在周围k空间区域以欠采样率r
p
均匀移位采集。不同的tr间具有的对比度参数如b值、扩散编码方向、扩散时间等可以任意变化。在本实施例中，在图2中的tr1与tr2具有不同的扩散对比度参数。在图4中，给出了一个更为详细的实施例。不同的tr间具有连续变化的b值。
[0173]
(s102)对获取到的k空间数据进行数据预处理，进行奈奎斯特鬼影及梯形采样校正。奈奎斯特鬼影校正及梯形采样校正可以使用任意一种已知方法。本实施例中，使用采集机器umr790集成的方法进行处理。
[0174]
(s103)通过结合线圈-空间-对比度域张量低秩正则化项与数据保真项，构建目标方程。本实例优选待求图像在线圈-空间-对比度域的张量低秩性(locally low rank tensor,llrt)结合数据保真项形成目标函数，并选用交替方向乘子法作为迭代重建算法，获取所需的扩散图像。
[0175]
需要特别说明的是，本发明在重建时不需要获取线圈灵敏度参数；
[0176]
具体地，该目标函数可以描述为：
[0177][0178]
其中，y表示预处理后得到多维k空间数据，d表示欠采样算子，f表示傅里叶变换算子，x表示待重建的多线圈多对比度图像，α为正则化系数，bix为以第i个像素为中心的线圈-空间-对比度三维张量，r为该三维张量基于tucker分解或parafac分解下的求秩算子，nv为预先设定的值，如使用全局低秩约束则nv＝1，如使用局部低秩约束则nv可等于像素数量。根据本发明的上述快速多对比度磁共振扩散成像和重建方法，可以避免求解线圈灵敏图参数同时避开直接合并具有不同相位的k空间，提高了重建的图像质量，允许了更高的分辨率和更高的采集效率。
[0179]
(s104)本技术的实施例中，基于交替方向乘子法作为迭代重建算法，获取所需的扩散图像。使用上述方法前，需构造辅助变量，改写目标函数为下式。
[0180][0181]
s.t bix＝mi[0182]
mi＝gi×1ai×2bi×3ci，i＝1，2，
…
，nv[0183]
其中，y表示预处理后得到多维k空间数据，d表示欠采样算子，f表示傅里叶变换算子，x表示待重建的多线圈多对比度图像。α为正则化系数，bix为以第i个像素为中心的线圈-空间-对比度三维张量，r为该三维张量基于tucker分解或parafac分解下的求秩算子，nv为预先设定的值，如使用全局低秩约束则nv＝1，如使用局部低秩约束则nv可等于像素数量。mi为一个tensor类型的辅助变量，通过parafac或tucker分解得到，mi＝gi×1ai×2bi×3ci，i＝1，2，
…
，nv[0184]
根据本发明的一个实施例，目标函数可以改写为增广拉格朗日形式如下述公式：
[0185][0186]
其中，y表示预处理后得到多维k空间数据，d表示欠采样算子，f表示傅里叶变换算子，x表示待重建的多线圈多对比度图像。α为正则化系数，bix为以第i个像素为中心的线圈-空间-对比度三维张量，r为该三维张量基于tucker分解或parafac分解下的求秩算子，nv为预先设定的值，如使用全局低秩约束则nv＝1，如使用局部低秩约束则nv可等于像素数量。re为取实数算子。mi为一个tensor类型的辅助变量，通过parafac或tucker分解得到，mi＝gi×1ai×2bi×3ci，i＝1，2，
…
，nv，μ为增广拉格朗日系数，λi为拉格朗日乘子。
[0187]
(s105)在正式迭代重建之前，需要初始化待求解扩散图像，即确定图像初始值，本发明的实施例中，初始值可以设为0，为方便描述，在公式中描述为，x选择初始点x0。
[0188]
(s106)分别求解子问题。在admm中，需要分别求解只有待求解扩散图像变量，辅助变量以及拉格朗日乘子变量的子问题，下述子步骤分别描述了上文提到的对应的各个子问题的求解细节。
[0189]
求解xk，存在解析解
[0190][0191]
其中h表示共轭转置。
[0192]
求解子目标如下
[0193][0194]
优选地，首先进行张量的parafac或tucker分解
[0195][0196]
随后使用奇异值软阈值法：
[0197][0198][0199][0200]
求解
[0201][0202]
(s107)完成以上步骤后，自动地对每个通道的线圈图像进行更新，并进行循环迭代，达到下列条件之一则迭代停止：(1)迭代次数大于50次；(2)梯度的无穷范数小于1*10-10
；(3)步长小于1*10-12
。判断是否收敛或者达到最大迭代次数。
[0203]
(s108)如果没有满足迭代停止条件，重新进入步骤s106继续迭代。
[0204]
(s109)根据采集的扩散加权参数类型，以及实际需求，求解adc、dti、dki、ivim等扩散定量成像，如果不需要求解定量图像，则磁共振扩散加权图像即为最终结果。
[0205]
上文公开了一种实施例来实现本技术的不同结构；所描述的实施例仅仅是本技术的一部分实施例，而不是全部的实施例。可以通过其他的方法与组合实现本发明的不同结构。
[0206]
使用本发明进行以下实施例，在此实施例中可以更加体现本发明的优点。
[0207]
使用可变对比度多次激发平面回波序列和根据本技术实施例的可变对比度keyhole平面回波序列成像方法对一位健康受试者进行头部多b值磁共振扩散成像对照实验。实验前进行了口头和书面的知情同意。扫描设备为联影医疗umr790 3t磁共振。使用32通道头线圈从b值0-1000s/mm2以50s/mm2为间隔采集21个不同b值的dwi全采样的k空间数据，每个b值下采集3个扩散方向。其他序列扫描参数为：
[0208]
视野/图像大小/tr/te/带宽/重复次数/翻转角＝240
×
240mm2/128
×
128/3000ms/132.3ms/1790hz/1/90
°
。随后分别以中心k空间区域密集采样，周围k空间区域以欠采样率r
p
＝4；中心k空间区域密集采样，周围k空间区域以欠采样率r
p
＝8，扫描前瞻性欠采样数据。其他序列扫描参数为：视野/图像大小/tr/te/带宽/重复次数/翻转角＝240
×
240mm2/240
×
240/3000ms/100.4ms/1790hz/1/90
°
。
[0209]
图5展示了低分辨率下周围k空间欠采样率为r
p
＝4和8时，对比的回顾性实验重建结果。对比方法命名为epi spa-llr和kepi llr,根据实施例重建的结果为kepi llrt(可变对比度keyhole-epi)。由图可知，kepi llrt显示出相比kepi llr和epi spa-llr更好的图像(更接近金标准)，特别epi spa-llr在部分区域存在较大误差。
[0210]
图6展示了高分辨率下周围k空间欠采样率为r
p
＝4和8时，对比的前瞻性实验重建结果以及对应的扩散参数图像。对比方法命名为epi spa-llr和kepi llr；根据实施例重建的结果为kepi llrt。epi spa-llr在图像的全局存在最高强度的残留伪影。根据实施构建的kepi llrt准确地反映了图像的细节，没有残留伪影，具有更高的信噪比。
[0211]
由此，本发明实现了快速高分辨磁共振扩散成像和重建，通过使用任意变化对比度keyhole平面回波序列进行k空间采集，并利用扩散加权图像在图像维度，对比度维度以及线圈维度的低秩性进行重建，避免了直接合并k空间造成的相位和运动伪影，避免了求解线圈灵敏度参数。提高了鲁棒性，实现了更高的分辨率，更高的信噪比，可以在扩散定量成像等应用中显著减少成像时间。在分辨率与采集效率上达成了平衡，在临床应用中具有重要的意义。
[0212]
本领域技术人员知道，除了以纯计算机可读程序代码方式实现本发明提供的系统、装置及其各个模块以外，完全可以通过将方法步骤进行逻辑编程来使得本发明提供的系统、装置及其各个模块以逻辑门、开关、专用集成电路、可编程逻辑控制器以及嵌入式微控制器等的形式来实现相同程序。所以，本发明提供的系统、装置及其各个模块可以被认为是一种硬件部件，而对其内包括的用于实现各种程序的模块也可以视为硬件部件内的结构；也可以将用于实现各种功能的模块视为既可以是实现方法的软件程序又可以是硬件部件内的结构。
[0213]
以上对本发明的具体实施例进行了描述。需要理解的是，本发明并不局限于上述特定实施方式，本领域技术人员可以在权利要求的范围内做出各种变化或修改，这并不影响本发明的实质内容。在不冲突的情况下，本技术的实施例和实施例中的特征可以任意相互组合。