基于ICU-Net的高倍欠采样磁共振图像的重建方法

技术领域

本发明涉及磁共振成像技术领域，主要用于实现对高倍欠采样磁共振图像的有效重建。

背景技术

磁共振成像(Magnetic Resonance Imaging，MRI)是目前临床医学影像中最重要的非侵入式检查方法之一，是继CT之后影像行业的又一大进步，成像质量优于B超和CT，但不会像CT等成像方式那样产生电离辐射；具有温度敏感性、多方位、多参数、多模态、多面成像能力等优点，不仅可显示人体组织的解剖信息，而且可显示功能信息。目前，MRI已经广泛地应用于血管造影、神经病学、心脏病学检查等多种临床应用，在介入治疗中得到广泛认可，被誉为现代医学影像技术皇冠上的明珠。

然而，成像速度过慢和介电伪影较高仍然是横亘在MRI领域多年的两大问题。MRI扫描时间过长、成像较慢会造成以下问题：(1)给病人造成额外的痛苦，尤其是对有“幽闭恐惧症”的病人；(2)器官的非自主运动(如呼吸、眨眼、吞咽等非自主运动)更容易造成图像模糊，增加伪影；(3)无法满足动态实时成像与手术导航的需要。

如何提升MRI的成像速度一直是磁共振人孜孜以求的目标。目前，MRI的加速方法主要有3种：并行MRI(Parallel MRI，PMRI)、压缩感知MRI(Compressed Sensing MRI，CS-MRI)、基于深度神经网络的MRI(Deep Neural Network MRI，DNN-MRI)。PMRI技术(如SENSE、GRAPPA等)利用多个线圈同时成像，但由于使用了欠采样，重建后图像的信噪比会随着加速因子的增加而降低；CS-MRI技术(如BM3D-MRI、CS-SENSE等)受限于数据采集、空间转换、重建算法和处理系统等因素，十几年来一直未能在磁共振领域内广泛应用。

近年来，随着深度学习理论和高性能计算的高速发展，DNN-MRI技术(如DeepCascadeNet、ConvRNN、PC-RNN等)无论是在重建质量还是在重建速度上都取得了突破性的提升。U-Net是现有的DNN-MRI技术所广泛使用的一种深度卷积网络结构，然而，U-Net存在如下问题：(1)参数规模往往比较大(动辄上百兆)，通常都是部署在高性能服务器上的，如果要部署到本地的MRI设备上就需要轻量级的网络结构以及相配套的序列；(2)需要大规模的训练数据，然而，由于受病人隐私保护的限制，获取大规模的训练数据往往并非易事；(3)对于高倍欠采样(如加速因子≥8时)的MRI，基于UNet的DNN-MRI技术通常会产生较为显著的混叠伪影，这限制了其在实际中的广泛使用。

发明内容

为了解决U-Net存在的三个问题，本发明基于稠密连接的思想提出了一种新的网络结构：基于强化连接的U-Net(Intensively Connected U-Net，ICU-Net)，对U-Net进行了改进。ICU-Net能够大幅度节省U-Net的参数规模，更有效地学习欠采样MRI图像中丢失的细节信息，提升MRI图像的重构质量。

本发明提供如下的技术方案：

一种基于ICU-Net的高倍欠采样磁共振图像的重建方法，包括以下步骤：

步骤1输入u倍欠采样的MRI图像

步骤2由ICU-Net f

进一步地，所述步骤1中，x

其中，

再进一步地，所述步骤2中，ICU-Net f

f

其中，Θ＝{Θ

具体的，公式(2)中的编码模块E(·；Θ

其中，K为编码模块的深度，

这里，

其中，z

其中，[·，·，…，·]表示拼接运算，即：将参与运算的各特征图在通道方向上拼接起来；

公式(2)中的解码模块D(·；Θ

其中，

这里，

公式(2)中的数据一致性函数f

其中，

进一步地，所述步骤2中，ICU-Net f

本发明的技术构思为：U-Net是现有的基于深度学习的快速磁共振成像方法所广泛使用的一种网络结构，但U-Net模型通常需要占用较大的存储空间，且难以有效重建高倍欠采样的MRI图像。对于高倍欠采样的MRI图像，网络所生成的浅层特征尤为重要。稠密连接的思想已成功应用于自然图像的识别、超分辨率重建等领域，不仅可以有效压缩网络参数的规模，而且可以充分利用网络的底层特征提升高层特征的非线性重构能力。因此，本发明基于稠密连接的思想，提出了一种新的网络结构：基于强化连接的U-Net(IntensivelyConnected U-Net，ICU-Net)。ICU-Net对U-Net进行了3点改进：(1)基于稠密连接和残差连接构建残差稠密模块(Residual Dense Block，RDB)，用RDB代替U-Net中编码模块和解码模块内部的卷积模块，强化编码阶段和解码阶段的特征学习能力；(2)将残差连接和特征图拼接相结合，代替U-Net中原有的特征图拼接，强化从编码模块到解码模块的连接；(3)将残差连接和数据一致性连接相结合，连接U-Net的输入端和输出端，强化U-Net的端到端的连接。ICU-Net能够大幅度节省U-Net的参数规模，更有效地学习欠采样MRI图像中丢失的细节信息，提升MRI图像的重构质量。

相对于现有的基于深度学习的快速磁共振成像算法，本发明的有益效果主要体现在：(1)ICU-Net由于使用了强化稠密连接，可以有效地进行特征重用，能够大幅度缩减网络参数的规模，从而能够较好地适用于较小规模训练样本的MRI图像的快速重建；(2)ICU-Net可以有效地利用网络低层模块所生成的特征，有效地提升高倍欠采样的磁共振图像的重建效果。

附图说明

图1是本发明所提出的ICU-Net的网络结构图，其中，Conv表示卷积核大小大于1×1的卷积运算，1×1Conv表示卷积核大小为1×1的卷积运算，Sum表示求和运算，Cat表示特征图拼接运算(详见图6)，ConvDn表示步长(stride)为2的卷积运算，Deconv为反卷积运算。

图2(a)为卷积模块(Convolutional Block，CB)的结构图，图2(b)为ICU-Net的核心模块之一“残差稠密模块(RDB)”的结构图。

图3(a)为U-Net的网络结构图，图3(b)为RDU-Net的网络结构图。

图4(a)和(b)分别为k空间中的8倍欠采样的MRI图像和全采样的MRI图像的示例，(c)和(d)分别为像素空间中的由零填充方法得到的MRI图像和全采样的MRI图像的示例。

图5(a)是中心欠采样掩码矩阵的示例，图5(b)是高斯欠采样掩码矩阵的示例。

图6是特征图拼接运算的示意图，即：将参与运算的各特征图在通道方向拼接起来。

图7是对数据一致性算子f

图8给出了MR训练图像集

图9为各种方法对8倍欠采样MRI图像的重构效果图：(a)全采样的MRI图像；(b)由0填充方法重构的MRI图像；(c)由U-Net重构的MRI图像；(d)由RDU-Net重构的MRI图像；(e)由ICU-Net重构的MRI图像。

图10为各种方法对16倍欠采样MRI图像的重构效果图：(a)全采样的MRI图像；(b)由0填充方法重构的MRI图像；(c)由U-Net重构的MRI图像；(d)由RDU-Net重构的MRI图像；(e)由ICU-Net重构的MRI图像。

具体实施方式

下面结合附图进一步说明本发明的实施步骤。

参照图1～图10，一种基于ICU-Net的快速磁共振成像方法，包括如下步骤：

步骤1输入u倍欠采样的MRI图像

步骤2由ICU-Netf

所述步骤1中，

其中，

所述步骤2中，ICU-Net f

f

其中，Θ＝{Θ

具体的，公式(4)中的编码模块E(·；Θ

其中，K为编码模块的深度，

这里，

其中，z

其中，[·，·，…，·]表示拼接运算，即：将参与运算的各特征图在通道方向上拼接起来，图6给出了拼接运算的示意图；

公式(4)中的解码模块D(·；Θ

其中，

这里，

公式(4)中的数据一致性函数f

其中，

所述步骤2中，ICU-Net f

进一步地，通过实验来验证本发明所提出的ICU-Net的有益效果。使用飞利浦Ingenia3T扫描仪对5名受试者进行扫描，获取T2WI模态的MRI图像，对其进行预处理，去除没有任何脑组织的切片，将每张MR切片裁剪为336×336×261大小。图8给出了所使用的的训练集

表1

表1对比了0填充、U-Net、RDU-Net和ICU-Net等四种方法在4倍、8倍、16倍等3种欠采样MRI图像上的定量评价结果。为了定量评价网络的重建性能，使用峰值信噪比(PSNR)和结构相似度(SSIM)对重构的MRI图像进行质量评估。可以看出，ICU-Net表现出了优于其它方法的重构精度。

图9和图10分别针对8倍和16倍欠采样MRI图像，对0填充、U-Net、RDU-Net和ICU-Net的重构效果进行了视觉对比，可以看出，ICU-Net在视觉上也表现出了优于其它方法的重构效果。