图像重建模型生成及图像重建方法、装置、设备和介质

1.本发明实施例涉及医学图像处理技术领域，尤其涉及一种图像重建模型生成及图像重建方法、装置、设备和介质。


背景技术：

2.磁共振成像(magnetic resonance imaging，mri)技术是一种先进的成像技术，被广泛应用于临床诊断和医学研究。然而，由于磁共振物理及硬件人体活动、耐受程度的制约，磁共振成像往往空间分辨率受限，无法非常精确评估部分心脑血管疾病，如一些细小的脑血管壁增厚、斑块情况。因此，提高磁共振成像的空间分辨率尤为重要。
3.高分辨率的磁共振图像对应了高密度的k空间采样矩阵，然而采集大量k空间点需要耗费大量时间。现有获取高分辨mri图像的方法主要分为两种方案，第一种是从压缩感知(compressed sensing,cs)重建的角度出发，利用mri图像的稀疏性等先验信息，对k空间进行高度欠采，以节省扫描时间，再利用cs重建算法重建出高分辨图像。另一种方案是基于计算机视觉中的图像超分辨技术，先采集和重建出低分辨的mri图像，然后利用图像超分辨算法生成高分辨率图像。该类方法有两个步骤，第一个步骤是采集低分辨率的k-空间信息，利用重建算法得到低分辨率的图像，第二步为将低分辨率mri图像输入图像超分辨算法模型，得到高分辨率mri图像。
4.但是，现有方法中，基于压缩感知重建的算法在加速倍数较高时，原始采样k空间点过于稀疏，且采样点分散在高和低频各部分，信噪比相对较低，难以重建出高清晰度的高分辨图像。而基于计算机视觉中的图像超分辨技术中，重建和超分辨是两个分开的步骤，没有交互，原始采集的k空间信息不直接参与超分辨步骤，第一步重建算法并不能百分之百还原低分辨图像的信息，会存在一定的误差，该误差会在第二步超分辨步骤中累计、放大，影响高分辨图像的质量。


技术实现要素：

5.本发明实施例提供了一种图像重建模型生成及图像重建方法、装置、设备和介质，以实现在在提高采样加速倍数的同时，能够重建得到高分辨率图像，提升图像质量。
6.第一方面，本发明实施例提供了一种图像重建模型生成方法，该方法包括：
7.获取全采样磁共振k空间数据，并基于预设低分辨率采样模型得到所述全采样磁共振k空间数据对应的欠采样k空间数据；
8.将所述欠采样k空间数据输入至，根据对所述预设低分辨率模型的反问题进行优化求解得到子问题，建立的超分重建生成网络模型，得到欠采样重建图像；
9.将所述欠采样重建图像和所述全采样磁共振k空间数据对应的全采样重建图像同时输入至与所述超分重建生成网络模型关联的判别网络模型，根据所述判别网络模型的输出结果训练所述超分重建生成网络模型，直到所述判别网络模型的损失函数满足预设条件时，完成所述超分重建生成网络模型训练，生成目标图像重建模型。
10.可选的，对所述预设低分辨率模型的反问题进行优化求解得到子问题，建立超分重建生成网络模型的过程包括：
11.对所述预设低分辨率模型的反问题在预设约束下进行优化求解，得到优化问题；
12.设置所述反问题的辅助变量，并确定所述优化问题的罚函数；
13.根据所述罚函数确定所述优化问题的两个子问题；
14.对所述两个子问题分别进行求解，并网络化得到所述超分重建生成网络模型。
15.可选的，所述对所述两个子问题分别进行求解，并网络化得到所述超分重建生成网络模型，包括：
16.采用迫近梯度法分别对所述两个子问题进行求解，得到所述辅助变量的解和所述反问题的解的迭代表示；
17.将所述迭代表示过程网络化，得到所述超分重建生成网络模型。
18.可选的，所述超分重建生成网络模型包括细化模块和数据一致性模块，其中，所述细化模块是与所述辅助变量相关联的神经网络模块。
19.可选的，所述判别网络模型为池化金字塔网络；所述损失函数为对抗生成网络梯度惩罚损失函数。
20.可选的，所述预设低分辨率采样模型由预设低分辨率欠采样算子、傅里叶变换算子及线圈敏感度参数确定。
21.第二方面，本发明实施例还提供了一种图像重建方法，该方法包括：
22.获取基于预设低分辨率采样模型，采样得到的欠采样k空间数据；
23.将所述k空间数据输入至由任一实施例所述的图像重建模型生成方法得到的目标图像重建模型中，得到所述欠采样k空间数据对应的重建图像。
24.第三方面，本发明实施例还提供了一种图像重建模型生成装置，该装置包括：
25.样本处理模块，用于获取全采样磁共振k空间数据，并基于预设低分辨率采样模型得到所述全采样磁共振k空间数据对应的欠采样k空间数据；
26.样本输入模块，用于将所述欠采样k空间数据输入至，根据对所述预设低分辨率模型的反问题进行优化求解得到子问题，建立的超分重建生成网络模型，得到欠采样重建图像；
27.模型生成模块，用于将所述欠采样重建图像和所述全采样磁共振k空间数据对应的全采样重建图像同时输入至与所述超分重建生成网络模型关联的判别网络模型，根据所述判别网络模型的输出结果训练所述超分重建生成网络模型，直到所述判别网络模型的损失函数满足预设条件时，完成所述超分重建生成网络模型训练，生成目标图像重建模型。
28.可选的，图像重建模型生成装置还包括模型构造模块，用于：对所述预设低分辨率模型的反问题进行优化求解得到子问题，建立超分重建生成网络模型。
29.可选的，所述模型构造模块具体用于：
30.对所述预设低分辨率模型的反问题在预设约束下进行优化求解，得到优化问题；
31.设置所述反问题的辅助变量，并确定所述优化问题的罚函数；
32.根据所述罚函数确定所述优化问题的两个子问题；
33.对所述两个子问题分别进行求解，并网络化得到所述超分重建生成网络模型。
34.可选的，所述模型构造模块还用于：
35.采用迫近梯度法分别对所述两个子问题进行求解，得到所述辅助变量的解和所述反问题的解的迭代表示；
36.将所述迭代表示过程网络化，得到所述超分重建生成网络模型。
37.可选的，所述超分重建生成网络模型包括细化模块和数据一致性模块，其中，所述细化模块是与所述辅助变量相关联的神经网络模块。
38.可选的，所述判别网络模型为池化金字塔网络；所述损失函数为对抗生成网络梯度惩罚损失函数。
39.可选的，所述预设低分辨率采样模型由预设低分辨率欠采样算子、傅里叶变换算子及线圈敏感度参数确定。
40.第四方面，本发明实施例还提供了一种图像重建装置，该装置包括：
41.数据获取模块，用于获取基于预设低分辨率采样模型，采样得到的欠采样k空间数据；
42.图像重建模块，用于将所述k空间数据输入至所述，由任一实施例所述的图像重建模型生成方法得到的目标图像重建模型中，得到所述欠采样k空间数据对应的重建图像。
43.第五方面，本发明实施例还提供了一种计算机设备，所述计算机设备包括：
44.一个或多个处理器；
45.存储器，用于存储一个或多个程序；
46.当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现如本发明任意实施例所提供的图像重建模型生成方法或图像重建方法。
47.第六方面，本发明实施例还提供了一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现如本发明任意实施例所提供的图像重建模型生成方法或图像重建方法。
48.上述发明中的实施例具有如下优点或有益效果：
49.本发明实施例，通过采用预设低分辨率采样模型对全采样磁共振k空间数据进行处理，得到对应的欠采样k空间数据，然后将欠采样k空间数据输入至，根据对预设低分辨率模型的反问题进行优化求解得到子问题建立的超分重建生成网络模型，得到欠采样重建图像；将欠采样重建图像和全采样磁共振k空间数据对应的全采样重建图像同时输入至与超分重建生成网络模型关联的判别网络模型，根据判别网络模型的输出结果训练超分重建生成网络模型，直到所述判别网络模型的损失函数满足预设条件时，完成超分重建生成网络模型训练，生成目标图像重建模型，以用于磁共振图像重建。解决了现有技术中磁共振图像重建加速倍数受限及重建图像质量有待提高的问题，实现了在提高采样加速倍数的同时，能够重建得到高分辨率图像，提升图像质量。
附图说明
50.图1是本发明实施例一提供的一种图像重建模型生成方法的流程图；
51.图2是本发明实施例一提供的一种超分重建网络的网络结构示意图；
52.图3是本发明实施例二提供的一种图像重建方法的流程图；
53.图4是本发明实施例三提供的一种图像重建模型生成装置的结构示意图；
54.图5是本发明实施例四提供的一种图像重建装置的结构示意图；
55.图6是本发明实施例五提供的一种计算机设备的结构示意图。
具体实施方式
56.下面结合附图和实施例对本发明作进一步的详细说明。可以理解的是，此处所描述的具体实施例仅仅用于解释本发明，而非对本发明的限定。另外还需要说明的是，为了便于描述，附图中仅示出了与本发明相关的部分而非全部结构。
57.实施例一
58.图1为本发明实施例一提供的一种图像重建模型生成方法的流程图，本实施例可适用于基于全采样磁共振k空间数据进行图像重建模型训练的情况。该方法可以由图像重建模型生成装置执行，该装置可以由软件和/或硬件的方式来实现，集成于具有应用开发功能的电子设备中。
59.如图1所示，图像重建模型生成方法包括以下步骤：
60.s110、获取全采样磁共振k空间数据，并基于预设低分辨率采样模型得到所述全采样磁共振k空间数据对应的欠采样k空间数据。
61.其中，全采样磁共振k空间数据是预先采集的样本数据，可以根据全采样的磁共振k空间数据重建得到高分辨率的磁共振图像。预设低分辨率模型预设低分辨率采样模型可以根据加速倍数的需求由预设低分辨率欠采样算子、傅里叶变换算子及线圈敏感度参数确定。
62.进一步的，加速倍数，是指采集原始数据，即全采样磁共振k空间数据，的加速倍数，也就是欠采倍数。可以通过全采样磁共振k空间数据的总数据点数除以采集的数据点数确定。在本实施例中，加速倍数可以达到12倍以上，也就是说，通过本实施例生成的图像重建模型，加速倍数在12倍以上也能够重建出分辨率足够高的重建图像。
63.具体的，预设低分辨率采样模型可以表示为y＝mhfcx
ꢀꢀ
(1)；其中，m为欠采样算子，f是傅里叶变换算子，c是线圈敏感度矩阵(sensitivity maps)，x∈cn和y∈cm分别是是高分辨率mri图像对应的k空间向量和低分辨率mri图像对应的欠采k空间向量，h:cn→cm
是k空间降采样算子，具体操作可以是截取高分辨率图像对应的k空间的中心部分信息(低频信息)，使之与低分辨图像的k空间对应。分辨率的高低决定于采样矩阵(采样算子)的尺寸，越低分辨率，采集的k空间点越集中在中心区域，也就是低频区域，分辨率提高，采样点扩散到更高频的区域，也就能获得更多高频细节，提升图像分辨率。在该步骤中，可以将全采样磁共振k空间数据作为x，c可以通过espirit或者walsh等算法估计出来，f就是傅里叶变换的算子，m是由0和1构成的向量，决定哪些点采集，哪些点不采集，不采集的点数值为0。
64.s120、将所述欠采样k空间数据输入至，根据对所述预设低分辨率模型的反问题进行优化求解得到子问题，建立的超分重建生成网络模型，得到欠采样重建图像。
65.其中，超分重建生成网络模型是基于预设低分辨率模型的反问题优化求解的结果构造而成的，具体的过程包括如下步骤：
66.首先，对预设低分辨率模型的反问题在预设约束下进行优化求解，得到优化问题。具体的，为了方便表示，可以将公式(1)中的各算子mhfc记作a。公式(1)对应的反问题就可以通过优化问题求解得到公式(2)，
[0067][0068]
其中，r(x)是广义的图像域约束。
[0069]
然后，设置上述反问题的优化求解的辅助变量，并确定优化问题的罚函数，可得到公式(3)。
[0070][0071]
其中，是惩罚项的一个系数，为了公式推导方便将分母设置为2，二范数求导数出来的2可以和它消掉。
[0072]
进一步的，可以利用交替线性化最小化算法，根据罚函数确定优化问题的两个子问题，得到公式(4)。求解算法之所以称之为交替线性化最小算法，是由于在求解过程中对进行了线性化处理，可以记其中，是线性化的二次项系数。这里需要说明的是，η可以用一个其他的字母代替，但是，在本实施例中，了后面推导出的公式(5)式在形式上能够简化，就取线性化的二次项系数取
[0073][0074]
进而，可以采用迫近梯度法分别对两个子问题分别进行求解，得到两个子问题的迭代步骤，如公式(5)所示，
[0075][0076]
其中，proxr是与图像域正则化算子r(即图像域约束)相关的迫近算子；η是更新步长；是图像k空间残差在图像域的表示，具体可以用公式6进行表示：
[0077][0078]
最终，将公式(5)中的迭代化过程网络化，展开即得到超分重建生成网络模型(super resolution reconstruction network,srr-net)，该网络也即一个生成器(generator)，该网络的结构可参考图2中的上面的网络结构，srr-net包含多个网络块，每一个网络块包含两个网络模块，分别是细化模块sk和数据一致性模块xk。具体的，网络化可表示为
[0079][0080]
其中，细化模块sk中，公式(5)中的迫近算子proxr被网络化为残差3d卷积神经网络(3d residual cnn)c。c是一个三层卷积神经网络，α控制残差的保留程度。α和步长γ都是可学习的网络参数，为模型带来了更大的灵活性。进一步的，观察公式(7)可确定，在本模型
中只需要初始化x0即可。在x0初始化后，可根据x0算出s1，进而根据s1算出x1，实现循环迭代过程。在本实施例中，x0的初始化数值取值为欠采样k空间数据对应的补零图数据。
[0081]
s130、将所述欠采样重建图像和所述全采样磁共振k空间数据对应的全采样重建图像同时输入至与所述超分重建生成网络模型关联的判别网络模型，根据所述判别网络模型的输出结果训练所述超分重建生成网络模型，直到所述判别网络模型的损失函数满足预设条件时，完成所述超分重建生成网络模型训练，生成目标图像重建模型。
[0082]
具体的，在该步骤引入判别网络，基于生成对抗网络的思路，对超分重建生成网络模型进行训练。在本实施例中，判别网络模型采用池化金字塔网络，采用对抗生成网络梯度惩罚损失函数(wgan-gp loss)作为池化金字塔网络的损失函数。在将欠采样重建图像和全采样磁共振k空间数据对应的全采样重建图像同时输入至判别网络模型，判别网络模型会判断欠采样重建图像与全采样重建图像之间的差别，当该差别对应的损失函数值满足预设的收敛条件时，即完成模型训练的过程，得到目标图像重建模型。
[0083]
本实施例的技术方案，通过采用预设低分辨率采样模型对全采样磁共振k空间数据进行处理，得到对应的欠采样k空间数据，然后将欠采样k空间数据输入至，根据对预设低分辨率模型的反问题进行优化求解得到子问题建立的超分重建生成网络模型，得到欠采样重建图像；将欠采样重建图像和全采样磁共振k空间数据对应的全采样重建图像同时输入至与超分重建生成网络模型关联的判别网络模型，根据判别网络模型的输出结果训练超分重建生成网络模型，直到所述判别网络模型的损失函数满足预设条件时，完成超分重建生成网络模型训练，生成目标图像重建模型，以用于磁共振图像重建。解决了现有技术中磁共振图像重建加速倍数受限及重建图像质量有待提高的问题，实现了在提高采样加速倍数的同时，能够重建得到高分辨率图像，提升图像质量。
[0084]
实施例二
[0085]
图2为本发明实施例二提供的一种图像重建方法的流程图，本实施例可适用于对欠采样的磁共振k空间数据进行医学图像重建的情况。该方法可以由图像重建装置执行，该装置可以由软件和/或硬件的方式来实现，集成于具有应用开发功能的计算机设备中。
[0086]
如图2所示，图像重建方法包括以下步骤：
[0087]
s210、获取基于预设低分辨率采样模型，采样得到的欠采样k空间数据。
[0088]
具体的，在进行磁共振扫描采样的过程中，可以根据采样速度的需求设置采样的模型，预设低分辨率模型由预设低分辨率欠采样算子、傅里叶变换算子及线圈敏感度参数确定。预设低分辨率欠采样算子决定了采样的加倍速度。
[0089]
预设低分辨率采样模型可以表示为y＝mhfcx；其中，m为欠采样算子，f是傅里叶变换算子，c是线圈敏感度矩阵(sensitivity maps)，x∈cn和y∈cm分别是是高分辨率mri图像对应的k空间向量和低分辨率mri图像对应的欠采k空间向量，h:cn→cm
是k空间降采样算子，具体操作可以是截取高分辨率图像对应的k空间的中心部分信息(低频信息)，使之与低分辨图像的k空间对应。分辨率的高低决定于采样矩阵(采样算子)的尺寸，越低分辨率，采集的k空间点越集中在中心区域，也就是低频区域，分辨率提高，采样点扩散到更高频的区域，也就能获得更多高频细节，提升图像分辨率。在该步骤中，可以将全采样磁共振k空间数据作为x，c可以通过espirit或者walsh等算法估计出来，f就是傅里叶变换的算子，m是由0和1构成的向量，决定哪些点采集，哪些点不采集，不采集的点数值为0。
[0090]
s220、将所述k空间数据输入至由任一实施例所述的图像重建模型生成方法得到的目标图像重建模型中，得到所述欠采样k空间数据对应的重建图像。
[0091]
由上述任一实施例所述的图像重建模型生成方法得到的目标图像重建模型，可以将低分辨k空间数据重建输出得到高分辨率图像。而且，欠采样的低分辨k空间数据加速倍数可以达到12倍以上，也就是说，通过本实施例生成的图像重建模型，加速倍数在12倍以上也能够重建出分辨率足够高的重建图像。
[0092]
本实施例的技术方案，通过预设低分辨率采样模型进行采样，再将采样数据输入至训练好的图像重建模型中，得到高分辨率的重建图像；解决了现有技术中磁共振图像重建加速倍数受限及重建图像质量有待提高的问题，实现了在提高采样加速倍数的同时，能够重建得到高分辨率图像，提升图像质量。
[0093]
实施例三
[0094]
图4为本发明实施例三提供的一种图像重建模型生成装置的结构示意图，本实施例可适用于基于全采样磁共振k空间数据进行图像重建模型训练的情况。
[0095]
如图4所示，图像重建模型生成装置包括样本处理模块310、样本输入模块320和模型生成模块330。
[0096]
其中，样本处理模块310，用于获取全采样磁共振k空间数据，并基于预设低分辨率采样模型得到所述全采样磁共振k空间数据对应的欠采样k空间数据；样本输入模块320，用于将所述欠采样k空间数据输入至，根据对所述预设低分辨率模型的反问题进行优化求解得到子问题，建立的超分重建生成网络模型，得到欠采样重建图像；模型生成模块330，用于将所述欠采样重建图像和所述全采样磁共振k空间数据对应的全采样重建图像同时输入至与所述超分重建生成网络模型关联的判别网络模型，根据所述判别网络模型的输出结果训练所述超分重建生成网络模型，直到所述判别网络模型的损失函数满足预设条件时，完成所述超分重建生成网络模型训练，生成目标图像重建模型。
[0097]
本实施例的技术方案，通过采用预设低分辨率采样模型对全采样磁共振k空间数据进行处理，得到对应的欠采样k空间数据，然后将欠采样k空间数据输入至，根据对预设低分辨率模型的反问题进行优化求解得到子问题建立的超分重建生成网络模型，得到欠采样重建图像；将欠采样重建图像和全采样磁共振k空间数据对应的全采样重建图像同时输入至与超分重建生成网络模型关联的判别网络模型，根据判别网络模型的输出结果训练超分重建生成网络模型，直到所述判别网络模型的损失函数满足预设条件时，完成超分重建生成网络模型训练，生成目标图像重建模型，以用于磁共振图像重建。解决了现有技术中磁共振图像重建加速倍数受限及重建图像质量有待提高的问题，实现了在提高采样加速倍数的同时，能够重建得到高分辨率图像，提升图像质量。
[0098]
可选的，图像重建模型生成装置还包括模型构造模块，用于：对所述预设低分辨率模型的反问题进行优化求解得到子问题，建立超分重建生成网络模型。
[0099]
可选的，所述模型构造模块具体用于：
[0100]
对所述预设低分辨率模型的反问题在预设约束下进行优化求解，得到优化问题；
[0101]
设置所述反问题的辅助变量，并确定所述优化问题的罚函数；
[0102]
根据所述罚函数确定所述优化问题的两个子问题；
[0103]
对所述两个子问题分别进行求解，并网络化得到所述超分重建生成网络模型。
[0104]
可选的，所述模型构造模块还用于：
[0105]
采用迫近梯度法分别对所述两个子问题进行求解，得到所述辅助变量的解和所述反问题的解的迭代表示；
[0106]
将所述迭代表示过程网络化，得到所述超分重建生成网络模型。
[0107]
可选的，所述超分重建生成网络模型包括细化模块和数据一致性模块，其中，所述细化模块是与所述辅助变量相关联的神经网络模块。
[0108]
可选的，所述判别网络模型为池化金字塔网络；所述损失函数为对抗生成网络梯度惩罚损失函数。
[0109]
可选的，所述预设低分辨率采样模型由预设低分辨率欠采样算子、傅里叶变换算子及线圈敏感度参数确定。
[0110]
本发明实施例所提供的图像重建模型生成装置可执行本发明任意实施例所提供的图像重建模型生成方法，具备执行方法相应的功能模块和有益效果。
[0111]
实施例四
[0112]
图5为本发明实施例四提供的一种图像重建装置的结构示意图，本实施例可适用于对欠采样的磁共振k空间数据进行医学图像重建的情况。
[0113]
如图5所示，图像重建装置包括数据获取模块410和图像重建模块420。
[0114]
其中，数据获取模块410，用于获取基于预设低分辨率采样模型，采样得到的欠采样k空间数据；图像重建模块420，用于将所述k空间数据输入至所述，由任一实施例所述的图像重建模型生成方法得到的目标图像重建模型中，得到所述欠采样k空间数据对应的重建图像。
[0115]
本实施例的技术方案，通过预设低分辨率采样模型进行采样，再将采样数据输入至训练好的图像重建模型中，得到高分辨率的重建图像；解决了现有技术中磁共振图像重建加速倍数受限及重建图像质量有待提高的问题，实现了在提高采样加速倍数的同时，能够重建得到高分辨率图像，提升图像质量。
[0116]
本发明实施例所提供的图像重建装置可执行本发明任意实施例所提供的图像重建方法，具备执行方法相应的功能模块和有益效果。
[0117]
实施例五
[0118]
图6为本发明实施例五提供的一种计算机设备的结构示意图。图6示出了适于用来实现本发明实施方式的示例性计算机设备12的框图。图6显示的计算机设备12仅仅是一个示例，不应对本发明实施例的功能和使用范围带来任何限制。计算机设备12可以任意具有计算能力的终端设备，如智能控制器及服务器、手机等终端设备。
[0119]
如图6所示，计算机设备12以通用计算设备的形式表现。计算机设备12的组件可以包括但不限于：一个或者多个处理器或者处理单元16，系统存储器28，连接不同系统组件(包括系统存储器28和处理单元16)的总线18。
[0120]
总线18表示几类总线结构中的一种或多种，包括存储器总线或者存储器控制器，外围总线，图形加速端口，处理器或者使用多种总线结构中的任意总线结构的局域总线。举例来说，这些体系结构包括但不限于工业标准体系结构(isa)总线，微通道体系结构(mac)总线，增强型isa总线、视频电子标准协会(vesa)局域总线以及外围组件互连(pci)总线。
[0121]
计算机设备12典型地包括多种计算机系统可读介质。这些介质可以是任何能够被
计算机设备12访问的可用介质，包括易失性和非易失性介质，可移动的和不可移动的介质。
[0122]
系统存储器28可以包括易失性存储器形式的计算机系统可读介质，例如随机存取存储器(ram)30和/或高速缓存存储器32。计算机设备12可以进一步包括其它可移动/不可移动的、易失性/非易失性计算机系统存储介质。仅作为举例，存储系统34可以用于读写不可移动的、非易失性磁介质(图6未显示，通常称为“硬盘驱动器”)。尽管图6中未示出，可以提供用于对可移动非易失性磁盘(例如“软盘”)读写的磁盘驱动器，以及对可移动非易失性光盘(例如cd-rom,dvd-rom或者其它光介质)读写的光盘驱动器。在这些情况下，每个驱动器可以通过一个或者多个数据介质接口与总线18相连。系统存储器28可以包括至少一个程序产品，该程序产品具有一组(例如至少一个)程序模块，这些程序模块被配置以执行本发明各实施例的功能。
[0123]
具有一组(至少一个)程序模块42的程序/实用工具40，可以存储在例如系统存储器28中，这样的程序模块42包括但不限于操作系统、一个或者多个应用程序、其它程序模块以及程序数据，这些示例中的每一个或某种组合中可能包括网络环境的实现。程序模块42通常执行本发明所描述的实施例中的功能和/或方法。
[0124]
计算机设备12也可以与一个或多个外部设备14(例如键盘、指向设备、显示器24等)通信，还可与一个或者多个使得用户能与该计算机设备12交互的设备通信，和/或与使得该计算机设备12能与一个或多个其它计算设备进行通信的任何设备(例如网卡，调制解调器等等)通信。这种通信可以通过输入/输出(i/o)接口22进行。并且，计算机设备12还可以通过网络适配器20与一个或者多个网络(例如局域网(lan)，广域网(wan)和/或公共网络，例如因特网)通信。如图所示，网络适配器20通过总线18与计算机设备12的其它模块通信。应当明白，尽管图6中未示出，可以结合计算机设备12使用其它硬件和/或软件模块，包括但不限于：微代码、设备驱动器、冗余处理单元、外部磁盘驱动阵列、raid系统、磁带驱动器以及数据备份存储系统等。
[0125]
处理单元16通过运行存储在系统存储器28中的程序，从而执行各种功能应用以及数据处理，例如实现本发实施例所提供的一种图像重建模型生成方法步骤，该方法包括：
[0126]
获取全采样磁共振k空间数据，并基于预设低分辨率采样模型得到所述全采样磁共振k空间数据对应的欠采样k空间数据；
[0127]
将所述欠采样k空间数据输入至，根据对所述预设低分辨率模型的反问题进行优化求解得到子问题，建立的超分重建生成网络模型，得到欠采样重建图像；
[0128]
将所述欠采样重建图像和所述全采样磁共振k空间数据对应的全采样重建图像同时输入至与所述超分重建生成网络模型关联的判别网络模型，根据所述判别网络模型的输出结果训练所述超分重建生成网络模型，直到所述判别网络模型的损失函数满足预设条件时，完成所述超分重建生成网络模型训练，生成目标图像重建模型。
[0129]
或者，还可以实现本发实施例所提供的一种图像重建方法步骤，该方法包括：
[0130]
获取基于预设低分辨率采样模型，采样得到的欠采样k空间数据；
[0131]
将所述k空间数据输入至由任一实施例所述的图像重建模型生成方法得到的目标图像重建模型中，得到所述欠采样k空间数据对应的重建图像。
[0132]
实施例六
[0133]
本实施例六提供了一种计算机可读存储介质，其上存储有计算机程序，该程序被
处理器执行时实现如本发明任意实施例所提供的图像重建模型生成方法，包括：
[0134]
获取全采样磁共振k空间数据，并基于预设低分辨率采样模型得到所述全采样磁共振k空间数据对应的欠采样k空间数据；
[0135]
将所述欠采样k空间数据输入至，根据对所述预设低分辨率模型的反问题进行优化求解得到子问题，建立的超分重建生成网络模型，得到欠采样重建图像；
[0136]
将所述欠采样重建图像和所述全采样磁共振k空间数据对应的全采样重建图像同时输入至与所述超分重建生成网络模型关联的判别网络模型，根据所述判别网络模型的输出结果训练所述超分重建生成网络模型，直到所述判别网络模型的损失函数满足预设条件时，完成所述超分重建生成网络模型训练，生成目标图像重建模型。
[0137]
或者，还可以实现本发实施例所提供的一种图像重建方法步骤，该方法包括：
[0138]
获取基于预设低分辨率采样模型，采样得到的欠采样k空间数据；
[0139]
将所述k空间数据输入至由任一实施例所述的图像重建模型生成方法得到的目标图像重建模型中，得到所述欠采样k空间数据对应的重建图像。
[0140]
本发明实施例的计算机存储介质，可以采用一个或多个计算机可读的介质的任意组合。计算机可读介质可以是计算机可读信号介质或者计算机可读存储介质。计算机可读存储介质例如可以是但不限于：电、磁、光、电磁、红外线、或半导体的系统、装置或器件，或者任意以上的组合。计算机可读存储介质的更具体的例子(非穷举的列表)包括：具有一个或多个导线的电连接、便携式计算机磁盘、硬盘、随机存取存储器(ram)、只读存储器(rom)、可擦式可编程只读存储器(eprom或闪存)、光纤、便携式紧凑磁盘只读存储器(cd-rom)、光存储器件、磁存储器件、或者上述的任意合适的组合。在本文件中，计算机可读存储介质可以是任何包含或存储程序的有形介质，该程序可以被指令执行系统、装置或者器件使用或者与其结合使用。
[0141]
计算机可读的信号介质可以包括在基带中或者作为载波一部分传播的数据信号，其中承载了计算机可读的程序代码。这种传播的数据信号可以采用多种形式，包括但不限于电磁信号、光信号或上述的任意合适的组合。计算机可读的信号介质还可以是计算机可读存储介质以外的任何计算机可读介质，该计算机可读介质可以发送、传播或者传输用于由指令执行系统、装置或者器件使用或者与其结合使用的程序。
[0142]
计算机可读介质上包含的程序代码可以用任何适当的介质传输，包括但不限于：无线、电线、光缆、rf等等，或者上述的任意合适的组合。
[0143]
可以以一种或多种程序设计语言或其组合来编写用于执行本发明操作的计算机程序代码，所述程序设计语言包括面向对象的程序设计语言，诸如java、smalltalk、c++，还包括常规的过程式程序设计语言，诸如“c”语言或类似的程序设计语言。程序代码可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中，远程计算机可以通过任意种类的网络，包括局域网(lan)或广域网(wan)，连接到用户计算机，或者，可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。
[0144]
本领域普通技术人员应该明白，上述的本发明的各模块或各步骤可以用通用的计算装置来实现，它们可以集中在单个计算装置上，或者分布在多个计算装置所组成的网络
上，可选地，他们可以用计算机装置可执行的程序代码来实现，从而可以将它们存储在存储装置中由计算装置来执行，或者将它们分别制作成各个集成电路模块，或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样，本发明不限制于任何特定的硬件和软件的结合。
[0145]
注意，上述仅为本发明的较佳实施例及所运用技术原理。本领域技术人员会理解，本发明不限于这里所述的特定实施例，对本领域技术人员来说能够进行各种明显的变化、重新调整和替代而不会脱离本发明的保护范围。因此，虽然通过以上实施例对本发明进行了较为详细的说明，但是本发明不仅仅限于以上实施例，在不脱离本发明构思的情况下，还可以包括更多其他等效实施例，而本发明的范围由所附的权利要求范围决定。