【技术实现步骤摘要】
一种基于双域傅里叶神经网络的金属伪影去除方法
[0001]本专利技术属于医学影像分析
，具体涉及一种CT金属伪影去除方法。
技术介绍
[0002]计算机断层扫描(CT)是广泛用于临床诊断和筛查的主要方式之一。当存在诸如金属植入物和牙科填充物之类的高密度物体时，CT扫描的原始数据，即弦图，由于射线束硬化和散射等物理因素而被破坏。由此产生的重建图像呈现出严重的射线状伪影，极大地限制了后续诊断。如何有效减少金属伪影仍然具有挑战性。随着深度神经网络的发展，该问题也越来越被关注。
[0003]为了解决这个问题，已有相关工作分别在弦图和图像域中去除金属伪影(MAR)。然而，在单个域中方法无法有效地从损坏的数据中恢复真实的组织。这是因为，一方面，重建图像上的金属伪影是全局存在的，当金属变大时伪影很难在图像域被去除。另一方面，虽然弦图中只有部分区域被破坏，但修复后的原始数据不能保证重建后的理想图像质量，因为弦图中微小的误差个不连续都容易导致图像域的二次伪影。
[0004]目前流行的MAR方法是通过一个可反向传播的Rad...
【技术保护点】
【技术特征摘要】
1.一种基于双域傅里叶神经网络的金属伪影去除方法，其特征在于，在弦图域，采用弦图傅立叶修复网络，借助快速傅里叶卷积在全弦图的感受野，利用全弦图未损坏区域修复损坏区域；同时在图像域，使用傅里叶U
‑
Net，通过图像域的上下文信息从局部到全局恢复图像；此外，引入多窗口损失来帮助网络在不同的动态范围内感知准确的反馈，以在临床应用中取得更好的结果；具体步骤如下：(1)在弦图域，利用基于快速傅里叶卷积的神经网络实现全局信息插值以修复被金属破坏的弦图；具体过程为：设人体二维CT切片的衰减系数分布是：X＝μ(i,j)，其中(i,j)表示二维坐标；设弦图S可以用通过投影求得，即类似地，通过定义反投影过程可以从弦图S重建图像X，即在此处是的逆运算；当金属植入物X
m
出现时，弦图沿金属迹线受损；假设弦图的大小为H
s
×
W
s
，其中H
s
和W
s
分别是检测器的数量和投影角；然后本使用来表示二进制金属轨迹，其中M＝1对应于弦图中被金属的损坏区域；于是，被金属影响的弦图S
mc
和无金属的弦图S之间的关系表示为：S
mc
⊙
(1
‑
M)＝S
⊙
(1
‑
M)
ꢀꢀ
(1)其中，符号
⊙
表示逐元素相乘；弦图域中的金属伪影去除问题是利用未被金属影响区域的信息恢复有金属轨迹标记的损坏区域；对于弦图的金属伪影问题，使用基于快速傅里叶卷积的弦图恢复网络FS
‑
Net，在弦图上利用全局感受野来准确地恢复被金属损坏的弦图；弦图恢复网络首先使用两个下采样降维，之后使用六个包含快速傅里叶卷积的弦图修复模块去除金属伪影，最后通过两个上采样层以得到恢复出的弦图；在每个弦图修复模块中，输入张量首先被分成两个分支，分别称为局部分支和全局分支；对于局部分支和全局分支，首先使用三个卷积模块来捕获多尺度信息；同时，对于全局分支，使用快速傅里叶卷积块在傅里叶域中恢复弦图；FS
‑
Net的输出写成：S
r
＝FS
‑
Net(S
mc
,M)
ꢀꢀ
(2)其中，S
r
是修复后的弦图；最后通过一个可反向传播的Radon层重建已经修复的弦图，以从弦图角度恢复的CT图像X
s
：X
s
＝RL(S
r
⊙
M+S
mc
⊙
(1
‑
M))
ꢀꢀ
(3)使用L1损失函数来衡量弦图域的差异，使用smooth
‑
L1损失函数来衡量图像域的差异，以避免过拟合和梯度爆炸；因此，优化FS
‑
Net的损失定义如下：其中，S
gt
和X
gt
分别是无金属的，真实的弦图和图像；(2)在图像域，首先使用一个局部伪影去除网络IU
‑
Net在图像域去除金属伪影；接着使用包含傅里叶跳连接的U
‑
Net:FR
‑
Net替换传统U
‑
Net，以更好去除二次伪影和全图不一致，具体过程如下：首先使用一个U
‑
Net在图像域初步去除金属伪影，该网络写为：X
u
＝LU
‑
Net(RL(S
mc
))
ꢀꢀ
(5)对于该网络，使用L1 loss作为其损失函数设已有从局部和全局两个角度的不同方法去除伪影的CT图像——X
s
和X
u
；采用一个傅
里叶恢复网络FR
‑
Net使用来自两个粗恢复的图像的信息以得到去伪影后准确的CT图像；所述的FR
‑
Net使用U
‑
Net架构，但使用局部和全局的信息进行重建，以进一步去除伪影，该网络模块关键在于...
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。