一种基于成像模型与机器学习的磁共振超分辨率成像方法与流程

本发明涉及一种磁共振成像方法，尤其是涉及一种基于成像模型与机器学习的磁共振超分辨率成像方法。

背景技术：

近年来，为了能节省磁共振成像时间，更好地发挥磁共振成像技术的在医学中的临床价值，越来越多的方法被研究而出，用于加速磁共振成像。在这其中，从单一的快速采集的低分辨率图像中重建出其对应的高分辨率图像的超分辨率成像技术是一项极具潜能的技术。而要想实现超分辨率成像，可靠的先验信息与实际采集的低频信息均需被充分利用。

目前主流的两大类超分辨率成像算法为，基于成像模型的超分辨率成像算法与基于机器学习的超分辨率算法。两者皆存在着客观的不足。基于成像模型的超分辨率成像算法，多采用最大后验概率模型进行求解。其中似然项融入了成像模型，实现了成像结果与采集的低频信息的一致性，先验项为人为构造的包含了先验信息的正则项，这类正则项大多形式单一，所能融入的先验信息十分有限，导致最终所成超分辨率图像高频细节不理想。基于机器学习的超分辨率算法，通过训练神经网络，实现了基于输入的低分辨率图像生成对应高分辨率图像的工作。该类算法能使神经网络从训练数据集中自主地学习并利用大量的先验信息，生成较为理想的高频细节，但却因为没有成像模型等限制，低频信息难以保证与采集信息一致，出现伪影等问题。另外，如何更好地学习先验信息，仍是基于机器学习的超分辨率算法所需考虑的问题。

技术实现要素：

本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种更好地融合了先验与采集的信息的基于成像模型与机器学习的磁共振超分辨率成像方法。

本发明的目的可以通过以下技术方案来实现：

一种基于成像模型与机器学习的磁共振超分辨率成像方法，包括以下步骤：

s1，数据集准备：准备低分辨率数据集和高分辨率数据集，其中，高分辨率数据集分为不同等级分辨率的多个；

s2，神经网络训练：按等级由低到高依次训练多个神经网络实现层进式学习，其中，第一个神经网络基于低分辨率图像生成最低等级的高分辨率图像，第二个神经网络基于低分辨率图像和第一个神经网络的结果，以此类推，最后一个神经网络基于低分辨率图像和上一个神经网络的结果生成最高等级的高分辨率数据集；

s3，成像阶段：通过磁共振成像采集得到相应的低分辨率图像，将该低分辨率图像依次通过已经训练好的神经网络生成高分辨图像，将生成的高分辨图像作为先验信息，将采集到的低分辨率图像代入广义级数模型，求解广义级数模型，重建出最终的超分辨率图像。

所述的低分辨率数据集获取方式为：将高分辨率数据集变换回k-空间，仅保留k-空间中心数据，再变换回图像域，得到低分辨率数据集。

所述的广义级数模型通过admm算法求解。

所述的数据集中的图像为大脑图像，高分辨率数据集中的不同等级依据图像纹理清晰程度划分。

所述的高分辨率数据集数量为四个，分别包含以下图像：高分辨率大脑白质和灰质的几何分割图像，高分辨率带白质纹理的大脑图像，高分辨率带白质和灰质纹理的大脑图像，高分辨率带白质、灰质及脑脊液纹理纹理的大脑图像。

所述的广义级数模型中，待求解的超分辨率图像x(s)表示如下：

其中，p(s)为神经网络所生成的高分辨率先验图像，n为基函数个数，c(n)为相应基函数权重，s为空间域坐标，△k为频域平移量。

所述的神经网络为unet神经网络。

所述的广义级数模型的求解过程为利用admm算法将整体优化求解问题分解为多个子优化问题进行求解，具体为：

其中y为采集低频信号，f为傅里叶变换算子，ω为k-空间采样算子，p为机器学习先验高分辨率图像，n为基函数个数，c(n)为相应基函数权重，δ为冲激函数，k为频域坐标，△k为频域平移量，||||1表示一范数，||||2表示二范数。

与现有技术相比，本发明具有以下优点：

(1)提升机器学习性能，生成更优的高分辨大脑先验图像：本技术中提出的“层进式”训练方式，相对于传统的直接训练神经网络基于低分辨率图像生成高分辨率图像而言，更加精细化，生成了更加准确的高分辨率边缘，纹理等信息，在峰值信噪比，均方误差，结构相似度这三类常用的衡量生成高分辨率图像质量的定量指标上均取得了显著优于传统机器学习方式的结果。

(2)更好地融合了先验与采集的信息，生成了更为准确的高分辨率大脑图像：传统基于成像模型的超分辨率成像技术，仅有好的数据一致性，而没有可靠的先验，高频细节不理想。传统基于机器学习的超分辨率成像技术仅有好的先验，但没有满足数据一致性需求，低频信息被破坏，出现伪影。本发明将神经网络生成的高分辨率大脑图像作为先验信息，与采集到的低频信息一同代入广义级数模型，实现了既包含可靠先验信息，又实现了成像结果与所采集信息一致的需求，最终生成了具有更为准确高低频信息，去除了伪影的超分辨率大脑图像。

附图说明

图1为本实施例基于成像模型与机器学习的磁共振超分辨率成像方法流程图；

图2为本实施整体成像步骤示意图；

图3(a)～3(d)分别为本实施例第一个～第四个神将网络学习过程示意图。

具体实施方式

下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案为前提进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于下述的实施例。

实施例

如图1所示，一种基于成像模型与机器学习的磁共振超分辨率成像方法，包括以下步骤：

s1，数据集准备：准备低分辨率数据集和高分辨率数据集，其中，高分辨率数据集分为不同等级分辨率的多个；

s2，神经网络训练：按等级由低到高依次训练多个神经网络实现层进式学习，其中，第一个神经网络基于低分辨率图像生成最低等级的高分辨率图像，第二个神经网络基于低分辨率图像和第一个神经网络的结果，以此类推，最后一个神经网络基于低分辨率图像和上一个神经网络的结果生成最高等级的高分辨率数据集；

s3，成像阶段：通过磁共振成像采集得到相应的低分辨率图像，将该低分辨率图像依次通过已经训练好的神经网络生成高分辨图像，将生成的高分辨图像作为先验信息，将采集到的低分辨率图像代入广义级数模型，通过admm算法求解广义级数模型，重建出最终的超分辨率图像。

低分辨率数据集获取方式为：将高分辨率数据集变换回k-空间，仅保留k-空间中心数据，再变换回图像域，得到低分辨率数据集。

本实施例中，数据集中的图像为大脑图像，高分辨率数据集中的不同等级依据图像纹理清晰程度划分。高分辨率数据集数量为四个，分别包含以下图像：高分辨率大脑白质和灰质的几何分割图像，高分辨率带白质纹理的大脑图像，高分辨率带白质和灰质纹理的大脑图像，高分辨率带白质、灰质及脑脊液纹理纹理的大脑图像。

相应的，第一个神经网络，被训练用于基于低分辨率大脑图像生成高分辨率大脑几何图像(大脑白质和灰质的分割图像，无纹理信息)；

第二个神经网络，被训练用于基于第一步结果与低分辨率图像生成带高分辨率白质纹理的大脑几何图像；

第三个神经网络，被训练用于基于第二步结果与低分辨率图像生成含高分辨率灰质与白质纹理的大脑图像；

第四个神经网络，被训练用于基于第三步结果与低分辨率图像生成含有高分辨率灰质、白质、脑脊液等纹理信息的最终高分辨率大脑图像。

整体成像步骤如图2所示，主要分为两步：

(1)将采集所得低分辨率图像经由一组经“层进式”训练而得的神经网络，生成高分辨率先验信息。

如图3(a)～3(d)所示，将采集到的低分辨率图片依次经过四个已经训练完成的神经网络，即可得到包含了高分辨率组织纹理，边界信息的先验高分辨率大脑图像。

(2)将先验高分辨率图像与采集所得低分辨率图像代入基于广义级数模型的成像模型用于最终高分辨率图像求解。

在广义级数模型中，待求解的高分辨率图像x(s)表示如下：

其中p(s)为网络所生成的高分辨率先验图像，n为基函数个数，c(n)为相应基函数权重，s为空间域坐标，△k为频域平移量。

基于广义级数模型，该技术引入成像模型来保证成像结果与采集的信息一致，求解过程即为如下优化过程：

其中y为采集低频信号，f为傅里叶变换算子，ω为k-空间采样算子，p为机器学习先验高分辨率图像，n为基函数个数，c(n)为相应基函数权重，δ为冲激函数，k为频域坐标，△k为频域平移量，||||1表示一范数，||||2表示二范数。

通过求解上式，便能得到满足了数据一致性需求，同时融入了充足先验信息的成像结果。