磁共振电影成像方法、装置、成像设备及存储介质与流程

本发明实施例涉及磁共振成像领域，尤其涉及一种磁共振电影成像方法、装置、成像设备及存储介质。

背景技术：

磁共振心脏电影成像是一种非侵入式的成像技术，能够用于评估心功能，室壁运动异常等，为心脏临床诊断提供丰富的信息。然而，由于磁共振物理、硬件和心脏运动周期时长的制约，磁共振心脏电影成像往往在时间和空间分辨率方面受限，无法准确评估部分心脏疾病，如心率不齐等。因此，在保证成像质量的前提下，提高磁共振心脏电影成像的速度和空间分辨率尤为重要。

近年来，很多人都在探索深度学习方法在磁共振电影成像领域的使用。比如，基于级联卷积网络的磁共振动态成像(dc-cnn)、卷积递归神经网络(crnn)以及多监督交叉域网络dimension均取得了良好的重建效果。但由于这三种神经网络均是直接学习从欠采图像到全采图像的映射关系，使得它们要么在图像重建过程需要较长的重建时间，要么重建出的磁共振心脏电影图像的质量较低。

综上，现有技术的深度学习方法无法在磁共振电影成像领域同时兼顾图像重建时间和图像质量。

技术实现要素：

本发明实施例提供了一种磁共振电影成像方法、装置、成像设备及存储介质，解决了现有技术的深度学习无法在磁共振电影成像领域同时兼顾图像重建时间和图像质量的问题。

第一方面，本发明实施例提供了一种磁共振电影成像方法，包括：

获取磁共振数据；

将所述磁共振数据输入至已训练的成像模型中以得到磁共振电影图像，其中，所述成像模型是基于admm算法构建的稀疏低秩网络模型，用于控制admm算法根据神经网络模型输出的迭代参数进行相应的迭代求解以得到磁共振电影图像。

进一步，所述神经网络模型用于根据admm算法前一次的迭代求解结果确定admm算法当前迭代求解所需的迭代参数。

进一步，所述成像模型的构建方法包括：

基于admm算法将磁共振信号的欠采重建任务转换为迭代求解数据一致子问题、低秩子问题、稀疏子问题和辅助变量子问题；

控制神经网络模型根据数据一致子问题、低秩子问题、稀疏子问题和辅助变量子问题的前一次的迭代求解结果，确定当前迭代求解所需的迭代参数，以及控制admm算法根据该迭代参数完成当前迭代求解过程，直至当前迭代求解结果符合预设收敛条件。

进一步，所述基于admm算法将磁共振信号的欠采重建任务转换为求解数据一致子问题、低秩子问题、稀疏子问题和辅助变量子问题，包括：

将磁共振信号的欠采重建任务建模为数据一致约束问题、变换域稀疏约束问题和低秩约束问题；

基于admm求解所述数据一致问题、变换域稀疏问题和低秩约束问题，以将磁共振信号的欠采重建任务转换为求解数据一致子问题、低秩子问题、稀疏子问题和辅助变量子问题。

进一步，所述数据一致子问题为：

所述低秩子问题为：

所述稀疏子问题为：

所述辅助变量更新子问题为：

其中，y为磁共振信号，x是图像数据，p是采样矩阵，f是傅里叶变换；z、t是辅助变量，z＝dx，t＝x，d是稀疏变换；α1为稀疏变换中的拉格朗日乘子，α2为低秩变换中的拉格朗日乘子，ρ1为稀疏变换中的惩罚系数，ρ2为低秩变换中的惩罚系数；ist用于表示求取xⁿ的奇异值阈值；s为非线性阈值函数；λ1是用于稀疏变换的正则化系数、λ2是用于低秩变换的正则化系数，η1是稀疏变换中的更新步长，η2是低秩变换中的更新步长。

进一步，迭代参数包括d、λ1和λ2。

进一步，所述磁共振数据为心脏磁共振数据。

第二方面，本发明实施例还提供了一种磁共振电影成像装置，包括：

获取模块，用于获取磁共振数据；

重建模块，用于控制神经网络模型根据admm算法的前一次的迭代求解结果确定当前迭代求解所需的迭代参数，以及控制admm算法根据该迭代参数完成当前迭代求解过程，直至当前迭代求解结果符合预设收敛条件。

第三方面，本发明实施例还提供了一种成像设备，所述成像设备包括：

一个或多个处理器；

存储装置，用于存储一个或多个程序；

当所述一个或多个程序被所述一个或多个处理器执行，使得所述一个或多个处理器实现如任意实施例所述的磁共振电影成像方法。

第四方面，本发明实施例还提供了一种包含计算机可执行指令的存储介质，所述计算机可执行指令在由计算机处理器执行时用于执行如任意实施例所述的磁共振电影成像方法。

本发明实施例提供的磁共振电影成像方法的技术方案，包括：获取磁共振数据；将磁共振数据输入至已训练的成像模型中以得到磁共振电影图像，其中，成像模型是基于admm算法构建的稀疏低秩网络模型，用于控制admm算法根据神经网络模型输出的迭代参数进行相应的迭代求解以得到磁共振电影图像。相较于现有技术。将admm算法与神经网络模型结合，使得神经网络模型可以学习磁共振数据在稀疏和低秩两方面的先验知识，并利用所学的这些先验知识快速准确地确定admm算法每次迭代求解所需的迭代参数，直至admm算法的迭代求解结果符合预设收敛条件，由于迭代参数的确定更加快速准确，因此图像重建过程所需时间大幅减少，重建出的磁共振电影图像的质量大幅提高。

附图说明

为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图做一简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1是本发明实施例一提供的磁共振电影成像方法的流程图；

图2是本发明实施例一提供的admm算法的迭代示意图；

图3a是本发明实施例二提供的磁共振电影成像装置的结构框图；

图3b是本发明实施例二提供的又一磁共振电影成像装置的结构框图；

图4是本发明实施例三提供的成像设备的结构框图。

具体实施方式

为使本发明的目的、技术方案和优点更加清楚，以下将参照本发明实施例中的附图，通过实施方式清楚、完整地描述本发明的技术方案，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

实施例一

图1是本发明实施例一提供的磁共振电影成像方法的流程图。本实施例的技术方案适用于使用基于admm算法构建的稀疏低秩网络模型快速重建出高质量的磁共振电影图像的情况。该方法可以由本发明实施例提供的磁共振电影成像装置来执行，该装置可以采用软件和/或硬件的方式实现，并配置在成像设备的处理器中应用。该方法具体包括如下步骤：

s101、获取磁共振数据。

其中，本实施例的磁共振数据为包含时间信息的磁共振信号。比如，心脏磁共振信号。

s102、将磁共振数据输入至已训练的成像模型中以得到磁共振电影图像，其中，成像模型是基于admm算法构建的稀疏低秩网络模型，用于控制admm算法根据神经网络模型输出的迭代参数进行相应的迭代求解以得到磁共振电影图像。

磁共振数据得到之后，将其输入已训练的成像模型中，由该已训练的成像模型对其分析，以得到磁共振电影图像。

其中，成像模型是基于admm算法构建的稀疏低秩网络模型，具体为将admm算法的迭代求解过程与神经网络模型相结合而成的模型。在该成像模型中，神经网络模型被设置为，根据admm算法前一次的迭代求解结果确定admm算法当前迭代求解所需的迭代参数，直至当前的迭代求解结果符合预设收敛条件，则神经网络模型结束迭代参数的计算，同时该迭代求解结果即为成像模型要输出的磁共振电影图像。

其中，成像模型的构建方法包括：基于admm算法将磁共振信号的欠采重建任务转换为迭代求解数据一致子问题、低秩子问题、稀疏子问题和辅助变量子问题；控制神经网络模型根据数据一致子问题、低秩子问题、稀疏子问题和辅助变量子问题的前一次的迭代求解结果，确定当前迭代求解所需的迭代参数，以及控制admm算法根据该迭代参数完成当前迭代求解过程，直至当前迭代求解结果符合预设收敛条件。

在一些实施例中，基于admm算法将磁共振信号的欠采重建任务转换为求解数据一致子问题、低秩子问题、稀疏子问题和辅助变量子问题的步骤，包括：将磁共振信号的欠采重建任务建模为数据一致约束问题、变换域稀疏约束问题和低秩约束问题；然后基于admm算法求解数据一致问题、变换域稀疏问题和低秩约束问题，以将磁共振信号的欠采重建任务转换为求解数据一致子问题、低秩子问题、稀疏子问题和辅助变量子问题。具体如下：

对于磁共振数据，即k空间数据其对应的欠采心脏电影图像为该心脏电影图像的欠采重建问题可被建模为以下优化问题：

其中，a＝pf为测量矩阵，p为采样矩阵，f为傅里叶变换，d为稀疏变换，g(·)为稀疏约束，||·||*＝∑i(σi,i)是核函数，σ是·的奇异值向量，核范数是对信号的前i个最大奇异值进行求和，反映了信号的低秩特性；λ1和λ2均是正则化系数。

引入辅助变量z＝dx，t＝x，则上述优化问题变为：

公式(2)的增广拉格朗日形式如下：

其中，α1为稀疏变换中的拉格朗日乘子，α2为稀疏变换中的拉格朗日乘子，ρ1为稀疏变换中的惩罚系数，ρ2为低秩变换中的惩罚系数。

利用admm算法对公式(3)进行求解，可以得到：

利用变换a＝fp，t⁽ⁿ⁺¹⁾使用奇异值阈值表示，则得到如下四个子问题。

其中，ist表示对信号进行svd分解得到特征值向量，并对特征值向量进行阈值操作，然后恢复到原信号；s为非线性阈值函数，用于对稀疏矩阵进行阈值过滤；λ1是用于稀疏变换的正则化系数，λ2是用于低秩变换的正则化系数，η1是稀疏变换中的更新步长，η2是低秩变换中的更新步长。

将公式(5)的迭代求解步骤参见图2，将其与神经网络模型相结合。结合之后，神经网络模型可以利用admm算法的前一次迭代求解结果确定其当前迭代求解所需的迭代参数，直至当前迭代求解结果符合预设收敛条件。可以理解的是，将admm算法与神经网络模型相结合，使得神经网络模型可以根据所学习的磁共振数据在稀疏和低秩两方面的先验知识，快速准确地确定出当前迭代求解所需的迭代参数。其中，迭代参数包括正则化系数λ1、λ2和稀疏变换d。

可以理解的是，成像模型构建完成之后，并不能直接使用其进行图像重建，还需要使用一定数量的样本对其进行训练，以生成已训练的成像模型。得到已训练的成像模型之后，即可使用该已训练的成像模型对磁共振数据进行图像重建，以生成磁共振电影图像。

本发明实施例提供的磁共振电影成像方法的技术方案，相较于现有技术。将admm算法与神经网络模型结合，使得神经网络模型可以学习磁共振数据在稀疏和低秩两方面的先验知识，并利用所学的这些先验知识快速准确地确定admm算法每次迭代求解所需的迭代参数，直至admm算法的迭代求解结果符合预设收敛条件，由于迭代参数的确定更加快速准确，因此图像重建过程所需时间大幅减少，重建出的磁共振电影图像的质量大幅提高。

实施例二

图3a是本发明实施例提供的磁共振电影成像装置的结构框图。该装置用于执行上述任意实施例所提供的磁共振电影成像方法，该装置可选为软件或硬件实现。该装置包括：

获取模块11，用于获取磁共振数据；

重建模块12，用于将所述磁共振数据输入至已训练的成像模型中以得到磁共振电影图像，其中，所述成像模型是基于admm算法构建的稀疏低秩网络模型，用于控制admm算法根据神经网络模型输出的迭代参数进行相应的迭代求解以得到磁共振电影图像。

可选地，该装置还包括模型构建模块101(参见图3b)，该模型构建模块包括：

任务转换单元，用于基于admm算法将磁共振信号的欠采重建任务转换为迭代求解数据一致子问题、低秩子问题、稀疏子问题和辅助变量子问题；

结合单元，用于控制神经网络模型根据admm算法的前一次的迭代求解结果确定当前迭代求解所需的迭代参数，以及控制admm算法根据该迭代参数完成当前迭代求解过程，直至当前迭代求解结果符合预设收敛条件。

其中，任务转换单元可选用于将磁共振信号的欠采重建任务建模为数据一致约束问题、变换域稀疏约束问题和低秩约束问题；基于admm求解所述数据一致问题、变换域稀疏问题和低秩约束问题，以将磁共振信号的欠采重建任务转换为求解数据一致子问题、低秩子问题、稀疏子问题和辅助变量子问题。

可选地，该装置还包括训练模块102(参见图3b)，用于接收训练样本数据，并根据所接收的训练样本数据完成成像模型的训练，以生成已训练的成像模型。

本发明实施例提供的磁共振电影成像装置的技术方案，相较于现有技术。将admm算法与神经网络模型结合，使得神经网络模型可以学习磁共振数据在稀疏和低秩两方面的先验知识，并利用所学的这些先验知识快速准确地确定admm算法每次迭代求解所需的迭代参数，直至admm算法的迭代求解结果符合预设收敛条件，由于迭代参数的确定更加快速准确，因此图像重建过程所需时间大幅减少，重建出的磁共振电影图像的质量大幅提高。

本发明实施例所提供的磁共振电影成像装置可执行本发明任意实施例所提供的磁共振电影成像方法，具备执行方法相应的功能模块和有益效果。

实施例三

图4为本发明实施例三提供的成像设备的结构示意图，如图4所示，该设备包括处理器201、存储器202、输入装置203以及输出装置204；设备中处理器201的数量可以是一个或多个，图4中以一个处理器201为例；设备中的处理器201、存储器202、输入装置203以及输出装置204可以通过总线或其他方式连接，图4中以通过总线连接为例。

存储器202作为一种计算机可读存储介质，可用于存储软件程序、计算机可执行程序以及模块，如本发明实施例中的磁共振电影成像方法、装置、成像设备及存储介质对应的程序指令/模块(例如，获取模块11和重建模块12)。处理器201通过运行存储在存储器202中的软件程序、指令以及模块，从而执行设备的各种功能应用以及数据处理，即实现上述的磁共振电影成像方法、装置、成像设备及存储介质。

存储器202可主要包括存储程序区和存储数据区，其中，存储程序区可存储操作系统、至少一个功能所需的应用程序；存储数据区可存储根据终端的使用所创建的数据等。此外，存储器202可以包括高速随机存取存储器，还可以包括非易失性存储器，例如至少一个磁盘存储器件、闪存器件、或其他非易失性固态存储器件。在一些实例中，存储器202可进一步包括相对于处理器201远程设置的存储器，这些远程存储器可以通过网络连接至设备。上述网络的实例包括但不限于互联网、企业内部网、局域网、移动通信网及其组合。

输入装置203可用于接收输入的数字或字符信息，以及产生与设备的用户设置以及功能控制有关的键信号输入。

输出装置204可包括显示屏等显示设备，例如，用户终端的显示屏。

实施例四

本发明实施例四还提供了一种包含计算机可执行指令的存储介质，所述计算机可执行指令在由计算机处理器执行时用于执行一种磁共振电影成像方法、装置、成像设备及存储介质，该方法包括：

获取模块，用于获取磁共振数据；

重建模块，用于控制神经网络模型根据admm算法的前一次的迭代求解结果确定当前迭代求解所需的迭代参数，以及控制admm算法根据该迭代参数完成当前迭代求解过程，直至当前迭代求解结果符合预设收敛条件。

当然，本发明实施例所提供的一种包含计算机可执行指令的存储介质，其计算机可执行指令不限于如上所述的方法操作，还可以执行本发明任意实施例所提供的磁共振电影成像方法中的相关操作。

通过以上关于实施方式的描述，所属领域的技术人员可以清楚地了解到，本发明可借助软件及必需的通用硬件来实现，当然也可以通过硬件实现，但很多情况下前者是更佳的实施方式。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在计算机可读存储介质中，如计算机的软盘、只读存储器(read-onlymemory，简称rom)、随机存取存储器(randomaccessmemory，简称ram)、闪存(flash)、硬盘或光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述的磁共振电影成像方法。

值得注意的是，上述磁共振电影成像装置的实施例中，所包括的各个单元和模块只是按照功能逻辑进行划分的，但并不局限于上述的划分，只要能够实现相应的功能即可；另外，各功能单元的具体名称也只是为了便于相互区分，并不用于限制本发明的保护范围。

注意，上述仅为本发明的较佳实施例及所运用技术原理。本领域技术人员会理解，本发明不限于这里所述的特定实施例，对本领域技术人员来说能够进行各种明显的变化、重新调整和替代而不会脱离本发明的保护范围。因此，虽然通过以上实施例对本发明进行了较为详细的说明，但是本发明不仅仅限于以上实施例，在不脱离本发明构思的情况下，还可以包括更多其他等效实施例，而本发明的范围由所附的权利要求范围决定。