深度学习MRI图像加速重建方法及系统.pdf

上传人:奶盖 文档编号:14521842 上传时间:2024-05-19 格式:PDF 页数:18 大小:1.04MB
收藏 版权申诉 举报 下载
深度学习MRI图像加速重建方法及系统.pdf_第1页
第1页 / 共18页
深度学习MRI图像加速重建方法及系统.pdf_第2页
第2页 / 共18页
深度学习MRI图像加速重建方法及系统.pdf_第3页
第3页 / 共18页
文档描述:

《深度学习MRI图像加速重建方法及系统.pdf》由会员分享,可在线阅读,更多相关《深度学习MRI图像加速重建方法及系统.pdf(18页完成版)》请在专利查询网上搜索。

1、(19)国家知识产权局(12)发明专利申请(10)申请公布号 (43)申请公布日 (21)申请号 202410044133.2(22)申请日 2024.01.12(71)申请人 北京航空航天大学杭州创新研究院地址 310051 浙江省杭州市滨江区长河街道炬航弄99号(72)发明人 张磊齐永兴谷双刘建伟关振宇(74)专利代理机构 北京易捷胜知识产权代理有限公司 11613专利代理师 齐胜杰(51)Int.Cl.G06T 11/00(2006.01)(54)发明名称一种深度学习MRI图像加速重建方法及系统(57)摘要本发明涉及MRI重建技术领域,尤其涉及一种深度学习MRI图像加速重建方法及系统,所。

2、述方法,包括:S1、获取MRI扫描图像集,并对所述MRI扫描图像集进行预处理,得到预处理后的MRI扫描图像集;S2、将预处理后的MRI扫描图像集,输入到已训练好的骨干模型中,所述骨干模型输出初步的MRI图像重建结果;所述骨干模型依次包括:灵敏度估计模块SMEB、双域交叉注意力改进模块DCRB以及N1个第一模块组;所述第一模块组依次包括:灵敏度改进模块SMRB和双域交叉注意力改进模块DCRB;S3、采用恢复模型ShiftNet作为细化网络,对初步的MRI图像重建结果进行细化处理,得到最终的MRI图像重建结果。权利要求书4页 说明书11页 附图2页CN 117557675 A2024.02.13C。

3、N 117557675 A1.一种深度学习MRI图像加速重建方法,其特征在于,包括:S1、获取MRI扫描图像集,并对所述MRI扫描图像集进行预处理,得到预处理后的MRI扫描图像集;S2、将预处理后的MRI扫描图像集,输入到已训练好的骨干模型中,所述骨干模型输出初步的MRI图像重建结果;所述骨干模型依次包括:灵敏度估计模块SMEB、双域交叉注意力改进模块DCRB以及N1个第一模块组;所述第一模块组依次包括:灵敏度改进模块SMRB和双域交叉注意力改进模块DCRB;S3、采用恢复模型ShiftNet作为细化网络,对初步的MRI图像重建结果进行细化处理,得到最终的MRI图像重建结果。2.根据权利要求1。

4、所述的深度学习MRI图像加速重建方法,其特征在于,所述预处理包括:图像增强处理、数据填充处理。3.根据权利要求1所述的深度学习MRI图像加速重建方法,其特征在于,所述方法在S1之前还包括:S0、获取初始训练数据集,并基于所述初始训练数据集,对预先获取的骨干模型进行训练,得到训练好的骨干模型。4.根据权利要求3所述的深度学习MRI图像加速重建方法,其特征在于,所述S0具体包括:S01、获取初始训练数据集,所述初始训练数据集为纽约大学fastMRI Initiative数据库的内部24线圈脑MRI数据集;所述24线圈脑MRI数据集包括:多个维度的MRI的图片;其中,多个维度包括:时间维度、切片维度。

5、、视图维度、对比度维度、加速因子维度;S02、使用基线模型CAUnet,采用相邻切片重建法ASR对所述初始训练数据集在指定维度上进行相邻K空间的重建处理,以填充缺失的K空间数据得到重建后的初始训练数据集;S03、采用预先获取的PromptUnet网络对重建后的初始训练数据集中的K空间域数据进行处理,得到相应的提示信息Prompt;所述重建后的初始训练数据集中的K空间域数据为指定的加速因子条件下的MRI的图片;S04、基于所述K空间域数据,以及与该K空间域数据所对应的图像域数据以及提示信息Prompt,对预先获取的骨干模型进行训练,直至满足预先指定的训练结束条件,得到训练好的骨干模型;所述提示信。

6、息Prompt包括:提示信息的空间域数据以及提示信息的图像域数据;预先指定的训练条件为:所述骨干模型的损失函数最小化。5.根据权利要求4所述的深度学习MRI图像加速重建方法,其特征在于,预先获取的PromptUnet网络是由基线模型CAUnet扩展得到的,具体为:在基线模型CAUnet每个级别的UpBlocks分别一一对应集成用于学习输入类型提示的PromptBlock模块以得到提示信息Prompt。6.根据权利要求5所述的深度学习MRI图像加速重建方法,其特征在于,所述骨干模型依次包括:灵敏度估计模块SMEB、双域交叉注意力改进模块DCRB以及N1权利要求书1/4 页2CN 11755767。

7、5 A2个第一模块组;所述灵敏度估计模块SMEB,用于接收所述K空间域数据,并针对所述K空间域数据,采用以大核卷积为内核的UNet网络对线圈敏感度SM进行粗略估计,得到线圈敏感度SM估计值,以及还用于将所得到的线圈敏感度SM估计值传输至与该灵敏度估计模块SMEB相邻的双域交叉注意力改进模块DCRB中;与该灵敏度估计模块SMEB相邻的双域交叉注意力改进模块DCRB,用于接收图像域数据,并根据第一输入信息以及所接收的由该灵敏度估计模块SMEB所传递的线圈敏感度SM估计值,对所述图像域数据、K空间域数据进行更新,得到更新后的图像域数据和更新后的K空间域数据,并将更新后的图像域数据、更新后的K空间域数。

8、据和所述线圈敏感度SM估计值传输至与该双域交叉注意力改进模块DCRB相邻的第一模块组中;所述第一输入信息包括:提示信息的空间域数据、K空间域数据、提示信息的图像域数据;所述第一模块组中的灵敏度改进模块SMRB,用于基于第二输入信息和所接收到的图像域数据、K空间域数据以及线圈敏感度SM估计值,对所接收到的线圈敏感度SM估计值进行更新,得到更新后的线圈敏感度SM估计值,并将更新后的线圈敏感度SM估计值和所接收到的图像域数据、K空间域数据传输至第一模块组中的双域交叉注意力改进模块DCRB中;所述第二输入信息包括K空间域数据;所述第一模块组中的双域交叉注意力改进模块DCRB,用于基于第一输入信息和所接。

9、收到的图像域数据、K空间域数据以及线圈敏感度SM估计值,对所接收到的图像域数据和K空间域数据进行更新,得到更新后的图像域数据、更新后的K空间域数据,并将更新后的图像域数据、更新后的K空间域数据和所接收到的线圈敏感度SM估计值传输至相邻的下一第一模块组中。7.根据权利要求6所述的深度学习MRI图像加速重建方法,其特征在于,其中,骨干模型中的第t+1个灵敏度改进模块SMRB,采用公式(1)对所接收到的线圈敏感度SM估计值进行更新,得到更新后的线圈敏感度SM估计值;所述公式(1)为:;其中,A*表示A的共轭转置;A=MF,且M是预先设定的二进制掩码;F代表傅里叶变换;是骨干模型中的第t个双域交叉注意。

10、力改进模块DCRB对所接收到的图像域数据进行更新,所得到的更新后的图像域数据;表示的共轭转置;表示K空间域数据;表示展开算子;是在对接收到的线圈敏感度SM估计值进行更新时所使用到的可学习的更新步长;权利要求书2/4 页3CN 117557675 A3表示第三预先设定常数;骨干模型中的第t个灵敏度改进模块SMRB所接收到的线圈敏感度SM估计值的正则化。8.根据权利要求7所述的深度学习MRI图像加速重建方法,其特征在于,所述骨干模型的损失函数为:;表示双域交叉注意力改进模块DCRB上全采样的基准真值;预先设定的第一系数;,其中N为骨干模型中双域交叉注意力改进模块DCRB的数量;是预先设定的固定常数。

11、;其中01;表示骨干模型中第t个双域交叉注意力改进模块DCRB所得到的更新后的K空间域数据;表示骨干模型中第N个双域交叉注意力改进模块DCRB所得到的更新后的K空间域数据;是K空间全采样的基准真值;是骨干模型中的第N个双域交叉注意力改进模块DCRB对所接收到的图像域数据进行更新,所得到的更新后的图像域数据;预先设定的第二系数;预先设定的第三系数;预先设定的第四系数;为基于K空间域数据获取的估计灵敏度;表示在估计灵敏度x的方向的变化率;表示在估计灵敏度y的方向的变化率。9.根据权利要求8所述的深度学习MRI图像加速重建方法,其特征在于,其中,在对所述骨干模型进行训练过程中,采用adam优化器对骨。

12、干模型的参数进行调整,从而使得所述骨干模型的损失函数最小化。10.一种深度学习MRI图像加速重建系统,其特征在于,包括:至少一个处理器;以及与所述处理器通信连接的至少一个存储器,其中,所述存储器存储有可被所述处理器权利要求书3/4 页4CN 117557675 A4执行的程序指令,所述处理器调用所述程序指令能够执行如权利要求19中任一所述的深度学习MRI图像加速重建方法。权利要求书4/4 页5CN 117557675 A5一种深度学习MRI图像加速重建方法及系统技术领域0001本发明涉及MRI重建技术领域,尤其涉及一种深度学习MRI图像加速重建方法及系统。背景技术0002磁共振成像(MRI)是。

13、一种无创、无辐射的医学成像技术,广泛用于临床疾病诊断、治疗和筛查。因为MRI需要多次的射频(RF)脉冲以填充K空间,与其他医学成像方式,如X射线或计算机断层扫描(CT)相比,通常需要更长的采集时间。最常用的方案是采用多线圈阵列在相位编码方向进行K空间欠采样。然而,在实际应用中,高加速系数会引起混叠干扰并降低重构图像的信噪比(SNR),显著降低图像质量的同时严重影响诊断。为了提高MRI图像的重建效果,人们提出了多种方法,包括部分傅里叶重建、并行成像(PI)、基于稀疏矩阵的压缩感知(CS)以及深度学习(DL)的模型。PI在加速度因子较低时表现良好,但在更高的扫描速率下出现残余混叠或噪声增强等问题。。

14、此外,基于CS的方法在加速系数变大时,容易出现模糊。发明内容0003鉴于现有技术的上述缺点、不足,本发明提供一种深度学习MRI图像加速重建方法及系统。0004为了达到上述目的,本发明采用的主要技术方案包括:第一方面,本发明实施例提供一种深度学习MRI图像加速重建方法,包括:S1、获取MRI扫描图像集,并对所述MRI扫描图像集进行预处理,得到预处理后的MRI扫描图像集;S2、将预处理后的MRI扫描图像集,输入到已训练好的骨干模型中,所述骨干模型输出初步的MRI图像重建结果;所述骨干模型依次包括:灵敏度估计模块SMEB、双域交叉注意力改进模块DCRB以及N1个第一模块组;所述第一模块组依次包括:灵。

15、敏度改进模块SMRB和双域交叉注意力改进模块DCRB;S3、采用恢复模型ShiftNet作为细化网络,对初步的MRI图像重建结果进行细化处理,得到最终的MRI图像重建结果。0005优选地,所述预处理包括:图像增强处理、数据填充处理。0006优选地,所述方法在S1之前还包括:S0、获取初始训练数据集,并基于所述初始训练数据集,对预先获取的骨干模型进行训练,得到训练好的骨干模型。0007优选地,所述S0具体包括:S01、获取初始训练数据集,所述初始训练数据集为纽约大学fastMRI Initiative数据库的内部24线圈脑MRI数据集;说明书1/11 页6CN 117557675 A6所述24线。

16、圈脑MRI数据集包括:多个维度的MRI的图片;其中,多个维度包括:时间维度、切片维度、视图维度、对比度维度、加速因子维度;S02、使用基线模型CAUnet,采用相邻切片重建法ASR对所述初始训练数据集在指定维度上进行相邻K空间的重建处理,以填充缺失的K空间数据得到重建后的初始训练数据集;S03、采用预先获取的PromptUnet网络对重建后的初始训练数据集中的K空间域数据进行处理,得到相应的提示信息Prompt;所述重建后的初始训练数据集中的K空间域数据为指定的加速因子条件下的MRI的图片;S04、基于所述K空间域数据,以及与该K空间域数据所对应的图像域数据以及提示信息Prompt,对预先获取。

17、的骨干模型进行训练,直至满足预先指定的训练结束条件,得到训练好的骨干模型;所述提示信息Prompt包括:提示信息的空间域数据以及提示信息的图像域数据;预先指定的训练条件为:所述骨干模型的损失函数最小化。0008优选地,预先获取的PromptUnet网络是由基线模型CAUnet扩展得到的,具体为:在基线模型CAUnet每个级别的UpBlocks分别一一对应集成用于学习输入类型提示的PromptBlock模块以得到提示信息Prompt。0009优选地,所述骨干模型依次包括:灵敏度估计模块SMEB、双域交叉注意力改进模块DCRB以及N1个第一模块组;所述灵敏度估计模块SMEB,用于接收所述K空间域数。

18、据,并针对所述K空间域数据,采用以大核卷积为内核的UNet网络对线圈敏感度SM进行粗略估计,得到线圈敏感度SM估计值,以及还用于将所得到的线圈敏感度SM估计值传输至与该灵敏度估计模块SMEB相邻的双域交叉注意力改进模块DCRB中;与该灵敏度估计模块SMEB相邻的双域交叉注意力改进模块DCRB,用于接收图像域数据,并根据第一输入信息以及所接收的由该灵敏度估计模块SMEB所传递的线圈敏感度SM估计值,对所述图像域数据、K空间域数据进行更新,得到更新后的图像域数据和更新后的K空间域数据,并将更新后的图像域数据、更新后的K空间域数据和所述线圈敏感度SM估计值传输至与该双域交叉注意力改进模块DCRB相邻。

19、的第一模块组中;所述第一输入信息包括:提示信息的空间域数据、K空间域数据、提示信息的图像域数据;所述第一模块组中的灵敏度改进模块SMRB,用于基于第二输入信息和所接收到的图像域数据、K空间域数据以及线圈敏感度SM估计值,对所接收到的线圈敏感度SM估计值进行更新,得到更新后的线圈敏感度SM估计值,并将更新后的线圈敏感度SM估计值和所接收到的图像域数据、K空间域数据传输至第一模块组中的双域交叉注意力改进模块DCRB中;所述第二输入信息包括K空间域数据;所述第一模块组中的双域交叉注意力改进模块DCRB,用于基于第一输入信息和所接收到的图像域数据、K空间域数据以及线圈敏感度SM估计值,对所接收到的图像。

20、域数据和K空间域数据进行更新,得到更新后的图像域数据、更新后的K空间域数据,并将更新后的图说明书2/11 页7CN 117557675 A7像域数据、更新后的K空间域数据和所接收到的线圈敏感度SM估计值传输至相邻的下一第一模块组中。0010优选地,其中,骨干模型中的第t+1个灵敏度改进模块SMRB,采用公式(1)对所接收到的线圈敏感度SM估计值进行更新,得到更新后的线圈敏感度SM估计值;所述公式(1)为:;其中,A*表示A的共轭转置;A=MF,且M是预先设定的二进制掩码;F代表傅里叶变换;是骨干模型中的第t个双域交叉注意力改进模块DCRB对所接收到的图像域数据进行更新,所得到的更新后的图像域数。

21、据;表示的共轭转置;表示K空间域数据;表示展开算子;是在对接收到的线圈敏感度SM估计值进行更新时所使用到的可学习的更新步长;表示第三预先设定常数;骨干模型中的第t个灵敏度改进模块SMRB所接收到的线圈敏感度SM估计值的正则化。0011优选地,所述骨干模型的损失函数为:0012表示双域交叉注意力改进模块DCRB上全采样的基准真值;预先设定的第一系数;,其中N为骨干模型中双域交叉注意力改进模块DCRB的数量;是预先设定的固定常数;其中01;表示骨干模型中第t个双域交叉注意力改进模块DCRB所得到的更新后的K空间域数据;表示骨干模型中第N个双域交叉注意力改进模块DCRB所得到的更新后的K空间域数据;。

22、说明书3/11 页8CN 117557675 A8是K空间全采样的基准真值;是骨干模型中的第N个双域交叉注意力改进模块DCRB对所接收到的图像域数据进行更新,所得到的更新后的图像域数据;预先设定的第二系数;预先设定的第三系数;预先设定的第四系数;为基于K空间域数据获取的估计灵敏度;表示在估计灵敏度x的方向的变化率;表示在估计灵敏度y的方向的变化率。0013优选地,其中,在对所述骨干模型进行训练过程中,采用adam优化器对骨干模型的参数进行调整,从而使得所述骨干模型的损失函数最小化。0014另一方面,本实施例还提供一种深度学习MRI图像加速重建系统,包括:至少一个处理器;以及与所述处理器通信连接。

23、的至少一个存储器,其中,所述存储器存储有可被所述处理器执行的程序指令,所述处理器调用所述程序指令能够执行如上述任一所述的深度学习MRI图像加速重建方法。0015本发明的有益效果是:本发明的一种深度学习MRI图像加速重建方法及系统,由于将预处理后的MRI扫描图像集,输入到已训练好的骨干模型中,所述骨干模型输出初步的MRI图像重建结果;所述骨干模型依次包括:灵敏度估计模块SMEB、双域交叉注意力改进模块DCRB以及N1个第一模块组;所述第一模块组依次包括:灵敏度改进模块SMRB和双域交叉注意力改进模块DCRB,本发明的一种深度学习MRI图像加速重建方法中的训练好的骨干模型是以灵敏度估计模块SMEB。

24、、双域交叉注意力改进模块DCRB以及灵敏度改进模块SMRB为主干,以Prompt为辅助的双流网络结构,充分利用了MRI图像和K空间特征信息,能够捕捉输入数据的局部特征和全局关联,提高了骨干模型对输入数据的理解和建模能力,从而能进一步优化MRI图像重建的质量。附图说明0016图1为本发明的一种深度学习MRI图像加速重建方法流程图;图2为现有的基线模型CAUnet结构示意图;图3为本发明实施例中的PromptUnet网络结构示意图;图4为本发明实施例中的骨干模型的结构示意图。说明书4/11 页9CN 117557675 A9具体实施方式0017为了更好的解释本发明,以便于理解,下面结合附图,通过具。

25、体实施方式,对本发明作详细描述。0018为了更好的理解上述技术方案,下面将参照附图更详细地描述本发明的示例性实施例。虽然附图中显示了本发明的示例性实施例,然而应当理解,可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了能够更清楚、透彻地理解本发明,并且能够将本发明的范围完整的传达给本领域的技术人员。0019实施例一0020参见图1,本实施例提供一种深度学习MRI图像加速重建方法,包括:S1、获取MRI扫描图像集,并对所述MRI扫描图像集进行预处理,得到预处理后的MRI扫描图像集;所述预处理包括:图像增强处理、数据填充处理。0021S2、将预处理后的MRI扫描图像。

26、集,输入到已训练好的骨干模型中,所述骨干模型输出初步的MRI图像重建结果;所述骨干模型依次包括:灵敏度估计模块SMEB、双域交叉注意力改进模块DCRB以及N1个第一模块组;所述第一模块组依次包括:灵敏度改进模块SMRB和双域交叉注意力改进模块DCRB;S3、采用恢复模型ShiftNet作为细化网络,对初步的MRI图像重建结果进行细化处理,得到最终的MRI图像重建结果。0022本实施例中通过使用强大的恢复模型ShiftNet作为细化网络,进一步探索图像域中的帧间/对比度一致性,以进行多帧/对比度特征聚合,该细化网络采用堆叠Unet和分组时空移位操作来扩展有效感受野。0023在本实施例的实际应用中。

27、,所述方法在S1之前还包括:S0、获取初始训练数据集,并基于所述初始训练数据集,对预先获取的骨干模型进行训练,得到训练好的骨干模型。0024所述S0具体包括:S01、获取初始训练数据集,所述初始训练数据集为纽约大学fastMRI Initiative数据库的内部24线圈脑MRI数据集;所述24线圈脑MRI数据集包括:多个维度的MRI的图片;其中,多个维度包括:时间维度、切片维度、视图维度、对比度维度、加速因子维度;S02、使用基线模型CAUnet,采用相邻切片重建法ASR对所述初始训练数据集在指定维度上进行相邻K空间的重建处理,以填充缺失的K空间数据得到重建后的初始训练数据集;本实施例中,使用。

28、CAUnet网络填充缺失的K空间数据,很大程度上消除图像混叠伪影。0025S03、采用预先获取的PromptUnet网络对重建后的初始训练数据集中的K空间域数据进行处理,得到相应的提示信息Prompt;所述重建后的初始训练数据集中的K空间域数据为指定的加速因子条件下的MRI的图片;说明书5/11 页10CN 117557675 A10S04、基于所述K空间域数据,以及与该K空间域数据所对应的图像域数据以及提示信息Prompt,对预先获取的骨干模型进行训练,直至满足预先指定的训练结束条件,得到训练好的骨干模型;所述提示信息Prompt包括:提示信息的空间域数据以及提示信息的图像域数据;预先指定的。

29、训练条件为:所述骨干模型的损失函数最小化。0026本实施例中的提示信息的空间域数据参见图4中的KR,本实施例中的提示信息的图像域数据参见图4中的XR。0027其中,预先获取的PromptUnet网络是由基线模型CAUnet扩展得到的,具体为:在基线模型CAUnet每个级别的UpBlocks分别一一对应集成用于学习输入类型提示的PromptBlock模块以得到提示信息Prompt。0028需要说明的是,本实施例中的基线模型CAUnet是一个具有通道注意力的Unet式网络,参见图2,基线模型CAUnet具有3级编码器解码器结构,每个级别都包含一个编码器、解码器以及对应的skip连接。同时,Bott。

30、leneckBlock也是神经网络中常见的结构,以降低计算复杂性,同时保持网络的表达能力,得到了本实施例中的预先获取的PromptUnet网络结构如图3所示。0029参见图4,所述骨干模型依次包括:灵敏度估计模块SMEB、双域交叉注意力改进模块DCRB以及N1个第一模块组;所述灵敏度估计模块SMEB,用于接收所述K空间域数据,并针对所述K空间域数据,采用以大核卷积为内核的UNet网络对线圈敏感度SM进行粗略估计,得到线圈敏感度SM估计值,以及还用于将所得到的线圈敏感度SM估计值传输至与该灵敏度估计模块SMEB相邻的双域交叉注意力改进模块DCRB中;与该灵敏度估计模块SMEB相邻的双域交叉注意力。

31、改进模块DCRB,用于接收图像域数据,并根据第一输入信息以及所接收的由该灵敏度估计模块SMEB所传递的线圈敏感度SM估计值,对所述图像域数据、K空间域数据进行更新,得到更新后的图像域数据和更新后的K空间域数据,并将更新后的图像域数据、更新后的K空间域数据和所述线圈敏感度SM估计值传输至与该双域交叉注意力改进模块DCRB相邻的第一模块组中;所述第一输入信息包括:提示信息的空间域数据、K空间域数据、提示信息的图像域数据;所述第一模块组中的灵敏度改进模块SMRB,用于基于第二输入信息和所接收到的图像域数据、K空间域数据以及线圈敏感度SM估计值,对所接收到的线圈敏感度SM估计值进行更新,得到更新后的线。

32、圈敏感度SM估计值,并将更新后的线圈敏感度SM估计值和所接收到的图像域数据、K空间域数据传输至第一模块组中的双域交叉注意力改进模块DCRB中;所述第二输入信息包括K空间域数据;所述第一模块组中的双域交叉注意力改进模块DCRB,用于基于第一输入信息和所接收到的图像域数据、K空间域数据以及线圈敏感度SM估计值,对所接收到的图像域数据和K空间域数据进行更新,得到更新后的图像域数据、更新后的K空间域数据,并将更新后的图像域数据、更新后的K空间域数据和所接收到的线圈敏感度SM估计值传输至相邻的下一第一模块组中。0030本实施例中训练好的骨干模型中最后一个第一模块组中的双域交叉注意力改进说明书6/11 页。

33、11CN 117557675 A11模块DCRB所输出的更新后的图像域数据就作为所述骨干模型输出初步的MRI图像重建结果。0031本实施例中整个骨干模型由1个灵敏度估计模块SMEB、双域交叉注意力改进模块DCRB以及N1个第一模块组组成。所述第一模块组依次包括:灵敏度改进模块SMRB和双域交叉注意力改进模块DCRB;这样,线圈敏感度SM就可以根据最新更新的图像来调整其“焦点”,同时重建的图像也可以受益于“重新聚焦”的线圈敏感度SM。0032在本实施例中,骨干模型中的第t+1个灵敏度改进模块SMRB,采用公式(1)对所接收到的线圈敏感度SM估计值进行更新,得到更新后的线圈敏感度SM估计值;其中,。

34、所述公式(1)为:;其中,A*表示A的共轭转置;A=MF,且M是预先设定的二进制掩码;F代表傅里叶变换;是骨干模型中的第t个双域交叉注意力改进模块DCRB对所接收到的图像域数据进行更新,所得到的更新后的图像域数据;表示的共轭转置;表示K空间域数据;表示展开算子;是在对接收到的线圈敏感度SM估计值进行更新时所使用到的可学习的更新步长;是骨干模型中的第t+1个灵敏度改进模块SMRB,对所接收到的线圈敏感度SM估计值进行更新,得到的更新后的线圈敏感度SM估计值;表示第三预先设定常数;骨干模型中的第t个灵敏度改进模块SMRB所接收到的线圈敏感度SM估计值的正则化。0033在本实施例的一种实施方式中,骨。

35、干模型中的第t+1个双域交叉注意力改进模块DCRB,是采用公式(2)对所接收到的图像域数据Xt进行更新,得到更新后的图像域数据xt+1;其中,所述公式(2)为:;其中,A*表示A的共轭转置;A=MF,且M是二进制掩码,其中非零元素表示数据采集;F代表傅里叶变换;是在对接收到的图像域数据进行更新时所使用到的可学习的更新步长;表示展开算子;表示缩减算子;表示第一预先设定常数;表示通过Swin Transformer得到的图像域数据的正则化;表示第二预先设定常数;表示骨干模型中第t个双域交叉注意力改进模块DCRB所得到的更新后的K空说明书7/11 页12CN 117557675 A12间域数据的正则。

36、化。0034本实施例中,将SwinT(Swin Transformer)与大核卷积分别作用于图像域数据和K空间域数据的双域交叉注意力改进模块DCRB。由于混叠通常会导致图像间的全局伪影,而SwinT实现跨层连接和多尺度特征融合,因此在图像域数据中采用SwinT执行操作。同时,由于K空间域数据中的每个元素编码全局感受野的信息,因此使用Unet形状的大核卷积来执行操作,用于对SwinT得到的结果进行有效改进和补充。本实施例中将Swin Transformer和大核卷积分别应用于图像域数据和K空间域数据的双域交叉注意力改进模块DCRB,能更好地利用参考模态中的先验信息,以改进目标模态的重构效果。通过。

37、引入注意力机制和大核卷积,骨干模型可以更好地捕获全局和局部的上下文信息,而双域交叉注意力改进模块DCRB可以有助于在两个模态之间传递和整合信息。0035在本实施例中,所述骨干模型的损失函数为:;表示双域交叉注意力改进模块DCRB上全采样的基准真值;预先设定的第一系数;,其中N为骨干模型中双域交叉注意力改进模块DCRB的数量;是预先设定的固定常数;其中01;表示骨干模型中第t个双域交叉注意力改进模块DCRB所得到的更新后的K空间域数据;表示骨干模型中第N个双域交叉注意力改进模块DCRB所得到的更新后的K空间域数据;是K空间全采样的基准真值;是骨干模型中的第N个双域交叉注意力改进模块DCRB对所接。

38、收到的图像域数据进行更新,所得到的更新后的图像域数据;预先设定的第二系数;预先设定的第三系数;预先设定的第四系数;为基于K空间域数据获取的估计灵敏度;表示在估计灵敏度x的方向的变化率;表示在估计灵敏度y的方向的变化率。说明书8/11 页13CN 117557675 A130036其中,在对所述骨干模型进行训练过程中,采用adam优化器对骨干模型的参数进行调整,从而使得所述骨干模型的损失函数最小化。0037另外,本实施例还提供一种深度学习MRI图像加速重建系统,包括:至少一个处理器;以及与所述处理器通信连接的至少一个存储器,其中,所述存储器存储有可被所述处理器执行的程序指令,所述处理器调用所述程。

39、序指令能够执行如实施例中所述的深度学习MRI图像加速重建方法。本实施例中的处理器为图形处理单元(GPU ,Graphics Processing Unit)。0038实施例二0039本实施例提供一种深度学习MRI图像加速重建方法,以缩短扫描时间并提高图像质量。本实施例二提供的一种深度学习MRI图像加速重建方法按照以下步骤实施:第1步骤、使用NYU fastMRI Initiative数据库的内部24线圈脑MRI数据集。其中,NYU fastMRI Initiative数据库是一个合作项目,旨在通过开发新的技术来加速磁共振成像(MRI),提高MRI扫描的效率。该数据库包含使用不同成像协议和扫描仪。

40、获得的MRI数据集。NYU fastMRI Initiative数据库的内部24线圈脑MRI数据集是由3T扫描仪(uMR 890)获取,包括62对T1FLAIR和T2w成对切片(T1FLAIR和T2w分别指代两种不同的MRI序列,T1FLAIR表示T1加权流体衰减反演恢复,T2w表示T2加权,T1代表一种MRI序列,T2也是一种MRI序列),维度为24320320。0040将24线圈脑MRI数据集分割,其中,25名被试(503个成对切片)用于训练,12名被试(241个成对切片)和25名被试(504个成对切片)分别用于验证和测试。最后采用最小最大归一化。0041第2步骤、将24线圈脑MRI数据集。

41、中每名被试(503个成对切片)进行随机刚性变换,包括在1.8,1.8之间的角度旋转和16,16之间的体素平移,其中为尺度因子 =1。然后对训练数据进行水平翻转。0042第3步骤、使用PyTorch实现提出的深度学习MRI图像加速重建方法。其中,优化器为Adam优化器,初始学习率为5e4,每三分之一的总epoch衰减一半。Minibatch大小为1。大脑数据集的训练epoch数为300。双域交叉注意力改进模块DCRB的数量为N =25。0至3分别设置为10、10、1和100。SwinT初始窗口大小:初始窗口大小为4。DCRB的patch大小:设置为4、2和1,对于N个DCRB,patch大小交替。

42、变化。0043第4步骤、获取初始训练数据集,使用NYU fastMRI Initiative数据库的内部24线圈脑MRI数据集。0044数据集包括多个维度的MRI图片,如时间维度、切片维度、视图维度、对比度维度、加速因子维度。0045使用基线模型CAUnet,采用相邻切片重建法ASR对初始训练数据集在指定维度上进行相邻K空间的重建处理。0046使用预先获取的PromptUnet网络对重建后的初始训练数据集中的K空间域数据进行处理,得到相应的提示信息Prompt。0047基于K空间域数据、图像域数据以及提示信息Prompt对预先获取的骨干模型进行训练,直至满足预先指定的训练结束条件,得到训练好的。

43、骨干模型。0048其中,骨干模型包括:灵敏度估计模块SMEB、双域交叉注意力改进模块DCRB以及N说明书9/11 页14CN 117557675 A141个第一模块组。第一模块组包括:灵敏度改进模块SMRB和双域交叉注意力改进模块DCRB。可以使用已有的深度学习网络架构,如UNet或ResNet等,进行骨干模型的训练。0049灵敏度估计模块SMEB接收K空间域数据,使用UNet网络对线圈敏感度SM进行粗略估计,得到线圈敏感度SM估计值,并传递至相邻的DCRB模块。0050双域交叉注意力改进模块DCRB接收图像域数据和来自SMEB的线圈敏感度SM估计值,对图像域数据和K空间域数据进行更新,传递至。

44、第一模块组。0051第5步骤、使用恢复模型ShiftNet对初步的MRI图像重建结果进行细化处理,得到最终的MRI图像重建结果。0052第6步骤、采用所述骨干模型的损失函数,结合SSIM等指标对模型性能进行评估。根据评估结果,进行骨干模型参数调整和优化,以提高重建质量和加速效果。0053第7步骤、完成骨干模型的训练和调优后,可以将其部署到实际的MRI扫描设备上,用于加速脑部MRI图像的重建。通过减少扫描时间,可以提高患者的舒适度,并在临床上更快地获取高质量的MRI图像。0054其中,所使用的处理器为12th Gen Intel(R)Core(TM)i512400F CPU,主频为2.50 GH。

45、z,内存32GB,显卡为GeForce RTX 3060。Ubuntu操作系统和PyCharm 2022。0055本实施例中深度学习MRI图像加速重建方法在脑部MRI数据上得到了详细的应用和调优,以满足医学影像重建的需求。0056本实施例中的一种深度学习MRI图像加速重建方法及系统,由于将预处理后的MRI扫描图像集,输入到已训练好的骨干模型中,所述骨干模型输出初步的MRI图像重建结果;所述骨干模型依次包括:灵敏度估计模块SMEB、双域交叉注意力改进模块DCRB以及N1个第一模块组;所述第一模块组依次包括:灵敏度改进模块SMRB和双域交叉注意力改进模块DCRB,本实施例中的一种深度学习MRI图像。

46、加速重建方法中的训练好的骨干模型是以灵敏度估计模块SMEB、双域交叉注意力改进模块DCRB以及灵敏度改进模块SMRB为主干,以Prompt为辅助的双流网络结构,充分利用了MRI图像和K空间特征信息,能够捕捉输入数据的局部特征和全局关联,提高了骨干模型对输入数据的理解和建模能力,从而能进一步优化MRI图像重建的质量。0057在本发明的描述中,需要理解的是,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”的特征可以明示或者隐含地包括一个或者更多个该特征。在本发明的描述中,“多个”的含义是两个或两个以上,除非另有。

47、明确具体的限定。0058在本发明中,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”、“固定”等术语应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或成一体;可以是机械连接,也可以是电连接;可以是直接相连,也可以通过中间媒介间接相连;可以是两个元件内部的连通或两个元件的相互作用关系。对于本领域的普通技术人员而言,可以根据具体情况理解上述术语在本发明中的具体含义。0059在本发明中,除非另有明确的规定和限定,第一特征在第二特征“上”或“下”,可以是第一和第二特征直接接触,或第一和第二特征通过中间媒介间接接触。而且,第一特征在第二特征“之上”、“上方”和“上面”,可以是第一特征在。

48、第二特征正上方或斜上方,或仅仅表示第一特征水平高度高于第二特征。第一特征在第二特征“之下”、“下方”和“下面”,可以说明书10/11 页15CN 117557675 A15是第一特征在第二特征正下方或斜下方,或仅仅表示第一特征水平高度低于第二特征。0060在本说明书的描述中,术语“一个实施例”、“一些实施例”、“实施例”、“示例”、“具体示例”或“一些示例”等的描述,是指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外,在不相互矛盾的情况下,本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。0061尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在本发明的范围内可以对上述实施例进行改动、修改、替换和变型。说明书11/11 页16CN 117557675 A16图 1图 2说明书附图1/2 页17CN 117557675 A17图 3图 4说明书附图2/2 页18CN 117557675 A18。

展开阅读全文
内容关键字: 深度 学习 MRI 图像 加速 重建 方法 系统
关于本文
本文标题:深度学习MRI图像加速重建方法及系统.pdf
链接地址:https://www.zhuanlichaxun.net/pdf/14521842.html
关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服 - 联系我们

copyright@ 2017-2018 zhuanlichaxun.net网站版权所有
经营许可证编号:粤ICP备2021068784号-1