一种基于残差网络的单扫描时空编码成像重建方法技术

一种基于残差网络的单扫描时空编码成像重建方法，涉及基于深度学习网络的磁共振图像重建技术。提供在单次扫描中获得二维图像，再使用深度学习方法重建得到一种基于残差网络的单扫描时空编码成像重建方法。将激励脉冲替换为线性扫频脉冲，有效抵抗由于不均匀磁场及化学位移造成的图像扭曲，同时获得和EPI相似的成像速度、分辨率和信噪比。SPEN成像沿相位编码方向都是欠采样的。尽管时空编码成像信号本身无需重建就能够反映成像物的轮廓，但是该轮廓的固有分辨率通常是很低。利用深度学习从低分辨率的信号空间重建SPEN图像，大幅提高图像分辨率，呈现质子密度分布，在获得和传统反卷积重建方法相似的分辨率的同时，获得高的信噪比。

【技术实现步骤摘要】
一种基于残差网络的单扫描时空编码成像重建方法
本专利技术涉及基于深度学习网络的磁共振图像重建技术，具体上说，是涉及一种基于残差网络的单扫描时空编码成像重建技术。
技术介绍
磁共振成像(MRI)是一种无损分析物体内部组织结构信息的成像技术。在临床中，MRI在神经成像、心血管成像以及功能磁共振成像都有着极其重要的作用。磁共振成像是通过对质子施加一定频率的射频脉冲，同时结合外加的三维梯度场来检测所产生的信号。常规的多扫描磁共振序列得到一幅磁共振图像一般需要几分钟，甚至几十分钟，这对于MRI在临床中的应用是不合适的。因此，超快速成像技术不仅能够缩短采样时间，提高采样效率，更能在一定程度上提高图像质量，减弱运动伪影对图像的影响。随着超快速成像技术的快读发展，其中基于Mansfield提出的平面回波(EPI)采集方法可以超快地在单次扫描获取完整的k空间[1]，在灌注成像、动态成像和功能成像中得到了广泛应用。近年来，提出了一种新的基于时空编码(SpatiotemporalEncoding,SPEN)成像技术，可用与EPI相近的成像速度获取二维图像[2]。SPEN的k空间中，原本的相位编码的成像维度被取代为直接位置空间维，也就是只需要经过一维傅里叶变换就可以获得位置空间信息。SPEN成像可以获取与EPI类似的采样时间、分辨率、信噪比等，同时由于SPEN编码维的信号直接来自位置空间的包络，因此提高了抵抗不均匀场的能力，尤其是在全重聚的脉冲设计情形下。由于SPEN扫描方法通过扫频脉冲引入一个与空间位置有关的二次相位分布，结合稳定相位近似定理，可以通过简单的取模获得磁共振图像，但是与EPI通过傅里叶变换的方法相比，模值图像的空间分辨率下降严重。2010年，Ben-Eliezer等人发现，SPEN方法采样得到的信号存在冗余性，即每次采样的相位稳定区域存在重叠[3]。现有的利用信号冗余性进行的超分辨率重建算法，包括共轭梯度重建算法，部分傅里叶重建算法[4]，反卷积重建算法等。虽然已有的方法都可以进行超分辨重建，但是算法运行时间较长，去噪效果不明显。深度学习，一种通过多个隐藏层有效学习输入数据和输出数据之间的非线性映射的算法，随着GPU性能快速提升，在近年取得了越来越多的关注。深度学习网络，特别是卷积神经网络，在各种问题的医学图像分析中广泛应用，包括分类、检测和分割。已有结果表明，基于卷积神经网络的方法在图像超分辨率重建方面远优于传统稀疏表示的方法[5]。作为深度神经网络模型的代表，残差网络能在网络层数很深时解决梯度爆炸/消失的问题。利用训练好的深度学习网络进行SPEN重建，可以在短时间内完成图像重建的过程，有效去除图像噪声，同时获得和传统反卷积重建方法相似的分辨率。参考文献：[1]StehlingMK,TurnerR,MansfieldP.Echo-PlanarImaging:MagneticResonanceImaginginaFractionofaSecond[J].Science,1991,254(5028):43-50.[2]ShrotY,FrydmanL.SpatiallyencodedNMRandtheacquisitionof2Dmagneticresonanceimageswithinasinglescan[J].JournalofMagneticResonance,2005,172(2):179-190.[3]Ben-EliezerN,ShrotY,FrydmanL,etal.Parametricanalysisofthespatialresolutionandsignal-to-noiseratioinsuper-resolvedspatiotemporallyencoded(SPEN)MRI[J].MagneticResonanceinMedicineOfficialJournaloftheSocietyofMagneticResonanceinMedicine,2014,72(2):418-429.[4]ChenY,LiJ,QuX,etal.PartialFouriertransformreconstructionforsingle-shotMRIwithlinearfrequency-sweptexcitation[J].MagneticResonanceinMedicineOfficialJournaloftheSocietyofMagneticResonanceinMedicine,2013,69(5):1326-1336.[5]BaumgartnerCF,OktayO,RueckertD.FullyConvolutionalNetworksinMedicalImaging:ApplicationstoImageEnhancementandRecognition[J].2017.
技术实现思路
本专利技术的目的在于提供在单次扫描中获得完整的二维图像，然后使用深度学习方法重建得到较高分辨的一种基于残差网络的单扫描时空编码成像重建方法。本专利技术包括以下步骤：1)在磁共振成像仪的操作软件中编写时空编码成像序列的代码，调试并编译通过。2)准备好实验样品，将其固定在磁共振成像仪的样品管上，放入磁共振成像仪的腔体中间；3)在磁共振成像仪的操作台上打开磁共振成像仪操作软件，首先用自旋回波序列进行定位，找到合适的成像区域，确定层选信息和感兴趣区域的大小。然后进行常规的调谐、匀场、频率校正和功率校正的操作；4)执行步骤1)所述的时空编码成像序列，设置采样带宽和相关参数，对步骤3)中的感兴趣区域进行采样，得到K空间数据；5)对步骤4)得到的K空间数据进行归一化、充零和一维傅里叶变换，得到图像域的数据；6)在数据模拟软件中生成随机样本，并在软件中对其进行模拟的采样，得到原始的SPEN的K空间数据，对K空间数据进行归一化处理，充零，和对频率编码维进行一维傅里叶变换，得到用于训练深度学习网络的训练数据；7)采用TensorFlow深度学习框架和Python搭建残差网络模型，设置好训练的相关参数；将步骤6)得到的训练数据输入网络进行训练，直至残差网络收敛并达到稳定，得到训练好的网络模型，然后利用训练好的网络模型对步骤5)得到的图像域数据进行重建，得到可靠的基于残差网络的单扫描时空编码成像。在步骤1)中，所述时空编码成像序列的结构依次为：翻转角为90°的线性扫频激发脉冲、翻转角为180°的重聚脉冲、移位梯度和采样回波链。在步骤6)中，所述随机样本可根据实验的待测样品的特征分布使用计算机批量随机生成，随机样本包含实验样本的所有特征，加入不稳定因素，提高网络模型对不理想实验环境的鲁棒性。在步骤7)中，所述训练好的网络可重建一幅图像只需要几十毫秒，远小于传统算法的重建时间。本专利技术将激励脉冲替换为线性扫频脉冲，可以有效抵抗由于不均匀磁场及化学位移造成的图像扭曲，同时获得和EPI相似的成像速度、分辨率和信噪比。SPEN成像沿相位编码方向都是欠采样的。尽管时空编码成像信号本身无需重建就能够反映成像物的轮廓，但是该轮廓的固有分辨率通常是很低的。本专利技术利用深度学习技术从低分辨率的信号空间重建SPEN图像，从而大幅提高图像的分辨率，更准确地呈现本文档来自技高网...

【技术保护点】
1.一种基于残差网络的单扫描时空编码成像重建方法，其特征在于包括以下步骤：1)在磁共振成像仪的操作软件中编写时空编码成像序列的代码，调试并编译通过；2)准备好实验样品，将其固定在磁共振成像仪的样品管上，放入磁共振成像仪的腔体中间；3)在磁共振成像仪的操作台上打开磁共振成像仪操作软件，首先用自旋回波序列进行定位，找到合适的成像区域，确定层选信息和感兴趣区域的大小；然后进行常规的调谐、匀场、频率校正和功率校正的操作；4)执行步骤1)所述的时空编码成像序列，设置采样带宽和相关参数，对步骤3)中的感兴趣区域进行采样，得到K空间数据；5)对步骤4)得到的K空间数据进行归一化、充零和一维傅里叶变换，得到图像域的数据；6)在数据模拟软件中生成随机样本，并在软件中对其进行模拟的采样，得到原始的SPEN的K空间数据，对K空间数据进行归一化处理，充零，和对频率编码维进行一维傅里叶变换，得到用于训练深度学习网络的训练数据；7)采用TensorFlow深度学习框架和Python搭建残差网络模型，设置好训练的相关参数；将步骤6)得到的训练数据输入网络进行训练，直至残差网络收敛并达到稳定，得到训练好的网络模型，然后利用训练好的网络模型对步骤5)得到的图像域数据进行重建，得到可靠的基于残差网络的单扫描时空编码成像。...

【技术特征摘要】
1.一种基于残差网络的单扫描时空编码成像重建方法，其特征在于包括以下步骤：1)在磁共振成像仪的操作软件中编写时空编码成像序列的代码，调试并编译通过；2)准备好实验样品，将其固定在磁共振成像仪的样品管上，放入磁共振成像仪的腔体中间；3)在磁共振成像仪的操作台上打开磁共振成像仪操作软件，首先用自旋回波序列进行定位，找到合适的成像区域，确定层选信息和感兴趣区域的大小；然后进行常规的调谐、匀场、频率校正和功率校正的操作；4)执行步骤1)所述的时空编码成像序列，设置采样带宽和相关参数，对步骤3)中的感兴趣区域进行采样，得到K空间数据；5)对步骤4)得到的K空间数据进行归一化、充零和一维傅里叶变换，得到图像域的数据；6)在数据模拟软件中生成随机样本，并在软件中对其进行模拟的采样，得到原始的SPEN的K空间数据，对K空间数据进行归一化处理，充零，和对频率编码维进行一维傅里叶变换...

【专利技术属性】
技术研发人员：陈忠，周甜甜，蔡聪波，曾坤，
申请(专利权)人：厦门大学，
类型：发明
国别省市：福建,35