基于深度展开神经网络的实时磁共振图像重建方法及系统技术方案

本发明专利技术提供了一种基于深度展开神经网络的实时磁共振图像重建方法及系统，包括如下步骤：步骤S1：搭建基于低秩和稀疏分解的深度展开神经网络；步骤S2：采用仿真数据训练搭建的深度展开神经网络；步骤S3：利用训练好的深度展开神经网络进行图像实时重建。本发明专利技术通过采用将基于低秩和稀疏分解的迭代算法展开到一个深度神经网络中的方式，解决了迭代算法应用于磁共振实时成像中重建时间较长的问题；实现磁共振介入图像的实时重建。磁共振介入图像的实时重建。磁共振介入图像的实时重建。

【技术实现步骤摘要】
基于深度展开神经网络的实时磁共振图像重建方法及系统


[0001]本专利技术涉及磁共振图像重建的
，具体地，涉及一种基于深度展开神经网络的实时磁共振图像重建方法及系统。

技术介绍

[0002]磁共振成像是一种医学成像技术，国际上从一九八二年才正式用于临床。它采用静磁场和射频磁场使人体组织成像，在成像过程中，既不用电子离辐射、也不用造影剂就可获得高对比度的清晰图像。它能够从人体分子内部反映出人体器官失常和早期病变。它在很多地方优于X线CT。虽然X
‑
CT解决了人体影像重叠问题，但由于提供的图像仍是组织对X射线吸收的空间分布图像，不能够提供人体器官的生理状态信息。当病变组织与周围正常组织的吸收系数相同时，就无法提供有价值的信息。只有当病变发展到改变了器官形态、位置和自身增大到给人以异常感觉时才能被发现。磁共振成像装置除了具备X线CT的解剖类型特点即获得无重叠的质子密度体层图像之外，还可借助核磁共振原理精确地测出原子核弛豫时间T1和T2，能将人体组织中有关化学结构的信息反映出来。这些信息通过计算机重建的图像是成分图像(化学结构像)，它有能力将同样密度的不同组织和同一组织的不同化学结构通过影像显示表征出来。这就便于区分脑中的灰质与白质，对组织坏死、恶性疾患和退化性疾病的早期诊断效果有极大的优越性，其软组织的对比度也更为精确。
[0003]公开号为CN112881958A的中国专利技术专利文献公开了一种基于低秩和稀疏分解的磁共振介入成像方法。包括：在介入过程中，采用黄金角径向采样方式连续采集k空间数据，对采集到的k空间数据进行分组，采用基于低秩和稀疏分解和framelet变换的方法重建磁共振介入图像。
[0004]公开号CN111932650A的中国专利技术专利文献公开了一种基于高通量深度展开网络的核磁共振图像重建方法，通过设计了一个高通量深度展开网络来进行训练和重建，该方法包括以下步骤：步骤1：构造训练数据集；步骤2：构造高通量深度展开网络；步骤3：高通量深度展开网络训练过程；步骤4：应用训练好的高通量深度展开网络进行压缩感知磁共振图像重建过程。
[0005]针对上述中的相关技术，专利技术人认为上述方法中的迭代算法，在算法中有多个正则化需要手动调参，非常繁琐耗时，同时需要较多的迭代步骤计算才能获得较好的重建效果，计算时间较长，不适用于实时磁共振介入成像。

技术实现思路

[0006]针对现有技术中的缺陷，本专利技术的目的是提供一种基于深度展开神经网络的实时磁共振图像重建方法及系统。
[0007]根据本专利技术提供的一种基于深度展开神经网络的实时磁共振图像重建方法，包括如下步骤：
[0008]步骤S1：搭建基于低秩和稀疏分解的深度展开神经网络；
[0009]步骤S2：采用仿真数据训练搭建的深度展开神经网络；
[0010]步骤S3：利用训练好的深度展开神经网络进行图像实时重建。
[0011]优选的，所述步骤S1包括如下步骤：
[0012]步骤S1.1：在磁共振成像中，从介入图像序列x中，将变化的背景和持续动态变化的介入特征分解为低秩矩阵L和稀疏矩阵S，即x＝L+S，通过利用L和S的空间稀疏特性，用于磁共振介入图像重建基于低秩和稀疏分解的模型；
[0013]步骤S1.2：将基于低秩和稀疏分解的模型展开到深度神经网络中；深度展开神经网络由多个迭代块组成，每个迭代块对应基于低秩和稀疏分解的迭代算法中的一次迭代，基于低秩和稀疏分解的迭代算法中的稀疏变换ψ由多个卷积层和激活函数组成的结构替换，卷积层中卷积核的权重参数在深度展开神经网络训练时学习，基于低秩和稀疏分解的迭代算法中的正则化参数在深度展开神经网络训练时学习；
[0014]步骤S1.3：第一个迭代块的第一输入为网络初始化后的参数，第二个及第二个之后的迭代块的第一输入均为上一个迭代块的输出，且每个迭代块的第二输入均包括径向采样得到的k空间数据和线圈敏感度图，最终经过所有迭代块后输出重建图像。
[0015]优选的，在所述步骤S1.1中，用于磁共振介入图像重建的基于低秩和稀疏分解的模型为：
[0016][0017]其中，E＝ΩFC是k空间编码算子，C是线圈敏感度，F是傅里叶变换，Ω为采样方式；d是获取的k空间数据；是沿着S时间方向的TV变换；ψ是稀疏变换，λ
L
、λ
s
、以及是正则化参数；k空间表示傅里叶空间。
[0018]优选的，所述步骤S2包括如下步骤：
[0019]步骤S2.1：通过仿真得到介入图像序列作为训练集，对介入图像序列作非均匀傅里叶变换，得到仿真的径向采样k空间数据；
[0020]步骤S2.2：将仿真的径向采样k空间数据和线圈敏感度图作为待训练的深度展开神经网络的输入，深度展开神经网络输出重建图像；
[0021]步骤S2.3：深度展开神经网络输出的重建图像与介入图像序列计算损失，根据损失通过反向传播的方式更新深度展开神经网络中的参数，深度展开神经网络中的参数包括卷积层中卷积核的权重，以及正则化参数；经过多轮更新训练后得到训练后的网络参数，即得到训练后的深度展开神经网络模型。
[0022]优选的，所述步骤S3包括如下步骤：
[0023]步骤S3.1：在进行实时成像之前，采集全采样的参考图像，用于估计线圈敏感度图；
[0024]步骤S3.2：在进行实时介入成像时，实时采集径向k空间数据，采集到的k空间数据和估计的线圈敏感度图作为已训练好的深度展开神经网络的输入，进而实时重建磁共振介入图像。
[0025]根据本专利技术提供的一种基于深度展开神经网络的实时磁共振图像重建系统，包括
如下模块：
[0026]模块M1：搭建基于低秩和稀疏分解的深度展开神经网络；
[0027]模块M2：采用仿真数据训练搭建的深度展开神经网络；
[0028]模块M3：利用训练好的深度展开神经网络进行图像实时重建。
[0029]优选的，所述模块M1包括如下模块：
[0030]模块M1.1：在磁共振成像中，从介入图像序列x中，将变化的背景和持续动态变化的介入特征分解为低秩矩阵L和稀疏矩阵S，即x＝L+S，通过利用L和S的空间稀疏特性，用于磁共振介入图像重建基于低秩和稀疏分解的模型；
[0031]模块M1.2：将基于低秩和稀疏分解的模型展开到深度神经网络中；深度展开神经网络由多个迭代块组成，每个迭代块对应基于低秩和稀疏分解的迭代算法中的一次迭代，基于低秩和稀疏分解的迭代算法中的稀疏变换ψ由多个卷积层和激活函数组成的结构替换，卷积层中卷积核的权重参数在深度展开神经网络训练时学习，基于低秩和稀疏分解的迭代算法中的正则化参数在深度展开神经网络训练时学习；
[0032]模块M1.3：第一个迭代块的第一输入为网络初始化后的参数，第二个及第二个之后的迭代块的第一输入均为上一个迭代块的输出，且每个迭代块的第二输入均包括径向采样得到的k空间数据和线圈敏感度图，最终经过所有迭代块后输出重建本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于深度展开神经网络的实时磁共振图像重建方法，其特征在于，包括如下步骤：步骤S1：搭建基于低秩和稀疏分解的深度展开神经网络；步骤S2：采用仿真数据训练搭建的深度展开神经网络；步骤S3：利用训练好的深度展开神经网络进行图像实时重建。2.根据权利要求1所述的基于深度展开神经网络的实时磁共振图像重建方法，其特征在于，所述步骤S1包括如下步骤：步骤S1.1：在磁共振成像中，从介入图像序列x中，将变化的背景和持续动态变化的介入特征分解为低秩矩阵L和稀疏矩阵S，即x＝L+S，通过利用L和S的空间稀疏特性，用于磁共振介入图像重建基于低秩和稀疏分解的模型；步骤S1.2：将基于低秩和稀疏分解的模型展开到深度神经网络中；深度展开神经网络由多个迭代块组成，每个迭代块对应基于低秩和稀疏分解的迭代算法中的一次迭代，基于低秩和稀疏分解的迭代算法中的稀疏变换ψ由多个卷积层和激活函数组成的结构替换，卷积层中卷积核的权重参数在深度展开神经网络训练时学习，基于低秩和稀疏分解的迭代算法中的正则化参数在深度展开神经网络训练时学习；步骤S1.3：第一个迭代块的第一输入为网络初始化后的参数，第二个及第二个之后的迭代块的第一输入均为上一个迭代块的输出，且每个迭代块的第二输入均包括径向采样得到的k空间数据和线圈敏感度图，最终经过所有迭代块后输出重建图像。3.根据权利要求2所述的基于深度展开神经网络的实时磁共振图像重建方法，其特征在于，在所述步骤S1.1中，用于磁共振介入图像重建的基于低秩和稀疏分解的模型为：其中，E＝ΩFC是k空间编码算子，C是线圈敏感度，F是傅里叶变换，Ω为采样方式；d是获取的k空间数据；是沿着S时间方向的TV变换；ψ是稀疏变换，λ
L
、λ
s
、以及是正则化参数；k空间表示傅里叶空间。4.根据权利要求1所述的基于深度展开神经网络的实时磁共振图像重建方法，其特征在于，所述步骤S2包括如下步骤：步骤S2.1：通过仿真得到介入图像序列作为训练集，对介入图像序列作非均匀傅里叶变换，得到仿真的径向采样k空间数据；步骤S2.2：将仿真的径向采样k空间数据和线圈敏感度图作为待训练的深度展开神经网络的输入，深度展开神经网络输出重建图像；步骤S2.3：深度展开神经网络输出的重建图像与介入图像序列计算损失，根据损失通过反向传播的方式更新深度展开神经网络中的参数，深度展开神经网络中的参数包括卷积层中卷积核的权重，以及正则化参数；经过多轮更新训练后得到训练后的网络参数，即得到训练后的深度展开神经网络模型。5.根据权利要求1所述的基于深度展开神经网络的实时磁共振图像重建方法，其特征在于，所述步骤S3包括如下步骤：步骤S3.1：在进行实时成像之前，采集全采样的参考图像，用于估计线圈敏感度图；
步骤S3.2：在进行实时介入成像时，实时采集径向k空间数据，采集到的k空间数据和估计的线圈敏感度图作为已训练好的深度展开神经网络的输入，进而实时重建磁共振介入图像。6.一种...

【专利技术属性】
技术研发人员：冯原，何钊，
申请(专利权)人：上海交通大学，
类型：发明
国别省市：