CN112085809A - 神经网络的训练方法、消除金属伪影的方法及其装置 - Google Patents
神经网络的训练方法、消除金属伪影的方法及其装置 Download PDFInfo
- Publication number
- CN112085809A CN112085809A CN202010971211.5A CN202010971211A CN112085809A CN 112085809 A CN112085809 A CN 112085809A CN 202010971211 A CN202010971211 A CN 202010971211A CN 112085809 A CN112085809 A CN 112085809A
- Authority
- CN
- China
- Prior art keywords
- neural network
- image
- metal
- projection
- projection data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 239000002184 metal Substances 0.000 title claims abstract description 123
- 229910052751 metal Inorganic materials 0.000 title claims abstract description 116
- 238000013528 artificial neural network Methods 0.000 title claims abstract description 75
- 238000012549 training Methods 0.000 title claims abstract description 47
- 238000000034 method Methods 0.000 title claims abstract description 45
- 238000012545 processing Methods 0.000 claims abstract description 54
- 238000003062 neural network model Methods 0.000 claims abstract description 35
- 238000001914 filtration Methods 0.000 claims abstract description 17
- 238000011176 pooling Methods 0.000 claims description 27
- 210000000988 bone and bone Anatomy 0.000 claims description 13
- 210000004872 soft tissue Anatomy 0.000 claims description 13
- 230000004913 activation Effects 0.000 claims description 12
- 238000002156 mixing Methods 0.000 claims description 9
- 210000001519 tissue Anatomy 0.000 claims description 7
- 230000008859 change Effects 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 6
- 230000011218 segmentation Effects 0.000 claims description 6
- XOOUIPVCVHRTMJ-UHFFFAOYSA-L zinc stearate Chemical compound [Zn+2].CCCCCCCCCCCCCCCCCC([O-])=O.CCCCCCCCCCCCCCCCCC([O-])=O XOOUIPVCVHRTMJ-UHFFFAOYSA-L 0.000 claims description 6
- 238000012937 correction Methods 0.000 claims description 5
- 230000001902 propagating effect Effects 0.000 claims description 4
- 238000001228 spectrum Methods 0.000 claims description 4
- 230000006835 compression Effects 0.000 claims description 3
- 238000007906 compression Methods 0.000 claims description 3
- 238000010606 normalization Methods 0.000 claims description 3
- 230000003313 weakening effect Effects 0.000 abstract description 2
- 238000002591 computed tomography Methods 0.000 description 25
- 230000000694 effects Effects 0.000 description 12
- 238000010586 diagram Methods 0.000 description 4
- 230000006872 improvement Effects 0.000 description 4
- 230000003902 lesion Effects 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000002238 attenuated effect Effects 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- ORILYTVJVMAKLC-UHFFFAOYSA-N Adamantane Natural products C1C(C2)CC3CC1CC2C3 ORILYTVJVMAKLC-UHFFFAOYSA-N 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 210000003484 anatomy Anatomy 0.000 description 1
- 238000013329 compounding Methods 0.000 description 1
- 238000013170 computed tomography imaging Methods 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 238000013135 deep learning Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000007943 implant Substances 0.000 description 1
- 238000002372 labelling Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 150000002739 metals Chemical class 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/008—Specific post-processing after tomographic reconstruction, e.g. voxelisation, metal artifact correction
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/136—Segmentation; Edge detection involving thresholding
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10072—Tomographic images
- G06T2207/10081—Computed x-ray tomography [CT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20081—Training; Learning
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20084—Artificial neural networks [ANN]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2210/00—Indexing scheme for image generation or computer graphics
- G06T2210/41—Medical
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Health & Medical Sciences (AREA)
- Evolutionary Computation (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明提供一种神经网络的训练方法、消除金属伪影的方法及其装置，所述神经网络模型包括第一神经网络、快速滤波反投影算子和第二神经网络；第一神经网络用于接收含有伪影的投影数据Sma和Sma对应的金属投影数据Mp，减弱Sma上的金属伪影并得到投影Sse；所述快速滤波反投影算子用于对投影Sse进行重建处理并得到图像Xse，以及对投影Sma进行重建处理并得到图像Xma；第二神经网络用于在图像域上、对图像Xse和Xma进行减弱二次金属伪影的处理，最终输出消除伪影的图像。基于若干投影数据Sma、以及每个Sma对应的金属投影数据Mp，以及相对应的不含伪影的投影标签Sgt和重建图像标签Xgt，对所述神经网络模型进行训练。该神经网络模型能够很好的去除金属伪影。
Description
技术领域
本发明涉及电子计算机断层扫描技术领域，尤其涉及一种神经网络的训练方法、消除金属伪影的方法及其装置。
背景技术
CT(Computed Tomography，电子计算机断层扫描)技术利用精确准直的射线(例如：X射线、γ射线或超声波等)、以及高灵敏度的探测器一同围绕人体的某一部位作连续的断面扫描，并生成该部位的影像，具有扫描时间快且图像清晰等优点。其过程通常可以归纳为：以不同的角度向人体射出射线，在经过人体时，这些射线会被人体吸收并衰减，之后，探测器会接受到衰减后的射线，并得到原始投影数据，然后通过图像重建(例如：滤波反投影等)算法生成最终的影像学图像。
在发明人的长期实践中，发现如果在CT成像视野范围内存在金属物体时，该金属物体会导致射线束硬化，同时噪声和散射效应被放大；甚至会导致射线完全衰减，没有光子达到接收器，导致没有有效信息用于重建。这些因素会导致或者加剧金属伪影。因此，设计一种去除金属伪影的方法，就成为一个亟待解决的问题。
发明内容
本发明的目的在于提供一种神经网络的训练方法、消除金属伪影的方法及其装置。
为了实现上述发明目的之一，本发明一实施方式提供了一种神经网络的训练方法，包括以下步骤：获取投影数据Sma、二值图Xm和Sma对应的金属投影数据Mp，所述投影数据Sma包含金属伪影，在所述二值图Xm中，值为1的像素值用于表征金属的形状；创建神经网络模型，所述神经网络模型包括第一神经网络、快速滤波反投影算子和第二神经网络；第一神经网络用于接收Sma和Mp，减弱Sma上的金属伪影并得到投影Sse；所述快速滤波反投影算子用于对投影Sse进行重建处理并得到图像Xse，以及对投影Sma进行重建处理并得到图像Xma，并且滤波反投影算子是可微的，能够将第二神经网络的训练中的梯度前向传播至第一神经网络；第二神经网络用于在图像域上，结合金属二值图Xm对图像Xse和Xma进行减弱二次金属伪影的处理；基于若干投影数据Sma、以及每个Sma对应的金属投影数据Mp，对所述神经网络模型进行有监督的训练。
作为本发明一实施方式的进一步改进，所述第一神经网络包括：第一编码器、第二编码器和第一解码器；第一编码器接收投影数据Sma，第一编码器包含有M个排成一列的编码层，其中，M为自然数；每个编码层均包含有卷积模块、激活模块和池化模块，编码层的输入经过卷积模块提取特征信息，再经过激活模块进行非线性变化，最后通过池化模块压缩高维的特征并输出所述特征；第二编码器接收金属投影数据Mp，第二编码器包含有N个排成一列的处理层，其中，N为自然数且N≤M，每个处理层都唯一对应到编码层；每个处理层均包含有池化模块，处理层的输入经过池化模块的平均池化处理、并得到池化结果；所述池化结果与对应的编码层的输出进行连接处理、并得到连接结果，所述连接结果作为对应的编码层的下一个编码层的输入；第一解码器包含有M个排成一列的解码层，解码层包含有卷积模块、激活模块和池化模块，编码层的输入经过卷积模块提取特征信息，再经过激活模块进行非线性变化，最后通过池化模块得到更高维的特征；第1个解码层接收第M个编码层的输出，第i个解码器接收第M-i+1编码层的输出和第i-1个解码层的输出的连接处理的结果；第一神经网络最后输出消除金属伪影的投影数据Sse。
作为本发明一实施方式的进一步改进，所述“获取投影数据Sma”具体包括：获取不包含伪影的CT图Xgt；将Xgt分割为软组织图Xw以及骨组织图Xb，对Xw和Xb进行处理并分别得到X’w和X’b，所述处理具体包括：在确定所述二值图Xm中的第一像素的值为1时，则将Xw中的与第一像素对应的第二像素的值置为0，将Xb中的与第一像素对应的第三像素的值置为0；基于正投影算子分别对X’w、X’b和Xm分别进行处理，并分别得到若干投影数据Pw、若干投影数据Pb和若干投影数据Pm，其中，不同的Pw表征软组织在不同能级下的投影数据，不同的Pb表征骨组织在不同能级下的投影数据，不同的Pm表征金属在不同能级下的投影数据；基于不同组织在不同能级下的衰减密度曲线、对Pw、Pb和Pm进行能谱密度归一化处理，从而得到不含其它噪声的投影数据S* ma；对S* ma叠加泊松噪声和散射噪声，之后，进行射束硬化校正处理，从而得到投影数据Sma。
作为本发明一实施方式的进一步改进，所述“将Xgt分割为软组织图Xw以及骨组织图Xb”具体包括：基于软阈值法，将Xgt分割为软组织图Xw以及骨组织图Xb。
作为本发明一实施方式的进一步改进，第二神经网络包括：第一卷积器、第二卷积器、第三卷积器和U-Net网络，首先对图像Xma和二值图Xm进行堆叠处理并得到图像X’ma，对图像Xse和二值图Xm进行堆叠处理从而分别得到图像X’se；之后，第一卷积器获取图像X’ma的融合有金属位置的第一特征层，第二卷积器获取图像X’se的融合有金属位置的第二特征层；之后，对第一、第二特征层进行堆叠处理，并将堆叠处理输入到U-Net网络，U-Net网络得到包含有多个通道第三特征层；第三卷积器得到通道数据为1的减弱伪影的结果图Xout。
本发明实施例提供了一种神经网络的训练装置，包括以下模块：数据获取模块，用于获取投影数据Sma、二值图Xm和Sma对应的金属投影数据Mp，所述投影数据Sma包含金属伪影，在所述二值图Xm中，值为1的像素值用于表征金属的形状；神经网络创建模块，用于创建神经网络模型，所述神经网络模型包括第一神经网络、快速滤波反投影算子和第二神经网络；第一神经网络用于接收Sma和Mp，减弱Sma上的金属伪影并得到投影Sse；所述快速滤波反投影算子用于对投影Sse进行重建处理并得到图像Xse，以及对投影Sma进行重建处理并得到图像Xma，并且滤波反投影算子是可微的，能够将第二神经网络的训练中的梯度前向传播至第一神经网络；第二神经网络用于在图像域上，结合金属二值图Xm对图像Xse和Xma进行减弱二次金属伪影的处理；训练模块，用于基于若干投影数据Sma、以及每个Sma对应的金属投影数据Mp，对所述神经网络模型进行有监督的训练。
本发明实施例提供了一种电子设备，包括：存储器，用于存储可执行指令；
处理器，用于执行所述存储器中存储的可执行指令时，实现上述的训练方法。
本发明实施例提供了一种计算机可读存储介质，存储有可执行指令，用于引起处理器执行时，实现上述的训练方法。
本发明实施例提供了一种消除金属伪影的方法，包括以下步骤：获取包含有金属伪影的CT图像Xma，通过阈值分割从Xma得到金属的二值图像Xm，将Xma和Xm通过正投影算子得到含金属伪影的投影数据Sma和金属投影数据Mp，并将Sma、Mp和Xm输入上述的神经网络模型，并得到减弱金属伪影的CT图像Xout。
本发明实施例提供了一种消除金属伪影的装置，包括以下模块：处理模块，用于获取包含有金属伪影的CT图像Xma，通过阈值分割从Xma得到金属的二值图像Xm，将Xma和Xm通过正投影算子得到含金属伪影的投影数据Sma和金属投影数据Mp，并将Sma、Mp和Xm输入上述的神经网络模型，并得到减弱金属伪影的CT图像Xout。
相对于现有技术，本发明的技术效果在于：本发明实施例提供一种神经网络的训练方法、消除金属伪影的方法及其装置，所述神经网络模型包括第一神经网络、快速滤波反投影算子和第二神经网络；第一神经网络用于接收含有伪影的投影数据Sma和Sma对应的金属投影数据Mp，减弱Sma上的金属伪影并得到投影Sse；所述快速滤波反投影算子用于对投影Sse进行重建处理并得到图像Xse，以及对投影Sma进行重建处理并得到图像Xma；第二神经网络用于在图像域上、对图像Xse和Xma进行减弱二次金属伪影的处理，最终输出消除伪影的图像。基于若干投影数据Sma、以及每个Sma对应的金属投影数据Mp，以及相对应的不含伪影的投影标签Sgt和重建图像标签Xgt，对所述神经网络模型进行训练。该神经网络模型能够很好的去除金属伪影。
附图说明
图1是本发明实施例中的生成投影数据的原理图；
图2是本发明实施例中的神经网络模型的结构示意图；
图3是本发明实施例中的神经网络的训练方法的流程示意图；
图4是本发明实施例中的神经网络模型的去除金属伪影的效果图；
图5是本发明实施例中的展示使用投影数据左右方进行循环padding的一个示例图。
具体实施方式
以下将结合附图所示的各实施方式对本发明进行详细描述。但这些实施方式并不限制本发明，本领域的普通技术人员根据这些实施方式所做出的结构、方法、或功能上的变换均包含在本发明的保护范围内。
本文使用的例如“上”、“上方”、“下”、“下方”等表示空间相对位置的术语是出于便于说明的目的来描述如附图中所示的一个单元或特征相对于另一个单元或特征的关系。空间相对位置的术语可以旨在包括设备在使用或工作中除了图中所示方位以外的不同方位。例如，如果将图中的设备翻转，则被描述为位于其他单元或特征“下方”或“之下”的单元将位于其他单元或特征“上方”。因此，示例性术语“下方”可以囊括上方和下方这两种方位。设备可以以其他方式被定向(旋转90度或其他朝向)，并相应地解释本文使用的与空间相关的描述语。
本发明实施例一提供了一种神经网络的训练方法，如图3所示，包括以下步骤：
步骤301：获取投影数据Sma、二值图Xm和Sma对应的金属投影数据Mp，所述投影数据Sma包含金属伪影，在所述二值图Xm中，值为1的像素值用于表征金属的形状；这里，Mp为金属的二值图Xm在当前几何下的投影，例如，使用正投影算子对二值图Xm进行处理，从而得到金属投影数据Mp。
步骤302：创建神经网络模型，所述神经网络模型包括第一神经网络、快速滤波反投影算子和第二神经网络；第一神经网络用于接收Sma和Mp，减弱Sma上的金属伪影并得到投影Sse；所述快速滤波反投影算子用于对投影Sse进行重建处理并得到图像Xse，以及对投影Sma进行重建处理并得到图像Xma，并且滤波反投影算子是可微的，能够将第二神经网络的训练中的梯度前向传播至第一神经网络；第二神经网络用于在图像域上，结合金属二值图Xm对图像Xse和Xma进行减弱二次金属伪影的处理；这里，第一神经网络1可以具有多个尺度，第一神经网络1能够在投影域上直接纠正原始数据中的异常值，第二神经网络2模型能够进一步减弱二次伪影，可以理解的是，该神经网络模型结合了投影法和影像法的优势。
这里，快速滤波反投影算子是一个使用并行计算实现的快速反投影重建算子，其几何参数与训练数据的几何参数一致，同时，快速滤波反投影算子是可微的，即在训练过程中，快速滤波反投影算子后面的第二网络结构的误差梯度可以通过快速滤波反投影算子传递至前面的第一网络结构，这一特性使得投影域和图像域的网络可以共同训练。
步骤303：基于若干投影数据Sma、以及每个Sma对应的金属投影数据Mp，对所述神经网络模型进行有监督的训练。
本实施例中，所述第一神经网络1包括：
第一编码器、第二编码器和第一解码器；
第一编码器接收投影数据Sma，第一编码器包含有M个排成一列的编码层，其中，M为自然数；每个编码层均包含有卷积模块、激活模块和池化模块，编码层的输入经过卷积模块提取特征信息，再经过激活模块进行非线性变化，最后通过池化模块压缩高维的特征并输出所述特征；
第二编码器接收金属投影数据Mp，第二编码器包含有N个排成一列的处理层，其中，N为自然数且N≤M，每个处理层都唯一对应到编码层；每个处理层均包含有池化模块，处理层的输入经过池化模块的平均池化处理、并得到池化结果；所述池化结果与对应的编码层的输出进行连接处理、并得到连接结果，所述连接结果作为对应的编码层的下一个编码层的输入；
第一解码器包含有M个排成一列的解码层，解码层包含有卷积模块、激活模块和池化模块，编码层的输入经过卷积模块提取特征信息，再经过激活模块进行非线性变化，最后通过池化模块得到更高维的特征；第1个解码层接收第M个编码层的输出，第i个解码器接收第M-i+1编码层的输出和第i-1个解码层的输出的连接处理的结果；
第一神经网络最后输出消除金属伪影的投影数据Sse。
这里，可以将特征信息融合在一起，在图像恢复的过程中起到保持特征细节的作用。
这里，可以根据投影数据的特性(投影几何为在一个圆周内均匀采样)，第一编码器11的padding操作在探测器方向上采用循环padding的操作，即，使用投影数据一侧的边缘数据部分对另一侧进行padding(见图5)。由于只需要恢复受到金属影响的投影数据，最后减弱金属伪影的投影数据定义为第一编码器11的输出和Sma的复合：Sse＝φSE(Sma,Mp)eMt+Smae(1-Mt)，其中，Mt为二值化的Mp，像素值为1的区域表示受金属伪影影响的区域，φSE代表第一神经网络。
本实施例中，所述“获取投影数据Sma”具体包括：
步骤1：获取不包含伪影的CT图Xgt；这里，如图1所示，该CT图Xgt和二值图Xm为在同一次实验中得到的数据，其中，CT图Xgt为一个不包含伪影的CT影像图，而二值图Xm描述了高密度金属的形状。可选的，在该二值图中，值为1的像素值构成一片区域，该区域表征了金属的形状。将不含伪影的图像标签Xgt经过正投影(FP，Forward Projection)算子得到不含伪影的投影数据标签Sgt。
步骤2：将Xgt分割为软组织图Xw以及骨组织图Xb，对Xw和Xb进行处理并分别得到X’w和X’b，所述处理具体包括：在确定所述二值图Xm中的第一像素的值为1时，则将Xw中的与第一像素对应的第二像素的值置为0，将Xb中的与第一像素对应的第三像素的值置为0；这里，可以理解的是，软组织的密度通常比较低，而骨组织的密度通常比较高；如图1所示，Xm、Xw和Xb是相对应的，因此，对于二值图Xm中的第一像素来说，其在Xw和Xb都会存在着对应的像素，分别为第二像素和第三像素。
步骤3：基于正投影算子分别对X’w、X’b和Xm分别进行处理，并分别得到若干投影数据Pw、若干投影数据Pb和若干投影数据Pm，其中，不同的Pw表征软组织在不同能级下的投影数据，不同的Pb表征骨组织在不同能级下的投影数据，不同的Pm表征金属在不同能级下的投影数据；
可选的，考虑到金属投影的局部容积效应(在CT扫描中，凡小于层厚的病变，其CT值受层厚内其它组织的影响，所测出的CT值不能代表病变的真正的CT值：如在高密度组织中较小的低密度病灶，其CT值偏高；反之，在低密度组织中的较小的高密度病灶，其CT值偏低，这种现象称为部分容积效应)，可以将Pm中的带状投影数据边缘处的像素值变为原来的1/4。
步骤4：基于不同组织在不同能级下的衰减密度曲线、对Pw、Pb和Pm进行能谱密度归一化处理，从而得到不含其它噪声的投影数据S* ma；可选的，对于不同类型的金属，可以使用不同的衰减密度曲线，来提高算法的鲁棒性。
步骤5：对S* ma叠加泊松噪声和散射噪声，之后，进行射束硬化校正(BHC，BeamHardening Correction)处理，从而得到投影数据Sma。
这里，泊松噪声是一种信号依赖噪声，对于一张图像而言，每个像素点的值都满足泊松分布，且每个像素点的泊松分布的均值是无噪图像在该像素点对应的值。可选的，可以使用MATLAB软件中的imnoise(I，'poisson')来叠加泊松噪声。
这里，散射效应也被称作康普顿效应。康普顿散射可以在任何物质中发生。当光子从光子源发出，射入散射物质(一般指金属)时，主要是与电子发生作用。如果光子的能量比较低，与电子束缚能同数量级时，则主要产生光电效应，原子吸收光子而产生电离。一般情况下，如果光子的能量相当大，远超过电子的束缚能时，电子与光子发生弹性碰撞，电子获得光子的一部分能量而反弹，失去部分能量的光子则改变方向飞出，而产生康普顿效应，即散射噪声。
射束硬化是由X射线束的多能性和衰减系数与能量的相关性造成的。多能X射线穿过物体时，能量较低的X光子由于光电效应优先被吸收，使得穿透物质的X射线高能成分比例增加，表现为平均能量升高，从而使射线随贯穿长度增加，变得更易穿透，频谱分布的峰值向较高的能量方向移动，这就是射束硬化效应，它会导致密度均匀的被测物体切片在重建CT图像上表现为亮度不一，图像上的像素值分布呈中间黑边缘亮的“茶杯”状，即“杯状”伪影。BHC校正处理可以为：多项式校正法、迭代法、单能法和双能法等。
本实施例中，所述“将Xgt分割为软组织图Xw以及骨组织图Xb”具体包括：基于软阈值法，将Xgt分割为软组织图Xw以及骨组织图Xb。
本实施例中，第二神经网络包括：第一卷积器、第二卷积器、第三卷积器和U-Net网络，首先对图像Xma和二值图Xm进行堆叠处理并得到图像X’ma，对图像Xse和二值图Xm进行堆叠处理从而分别得到图像X’se；之后，第一卷积器获取图像X’ma的融合有金属位置的第一特征层，第二卷积器获取图像X’se的融合有金属位置的第二特征层；之后，对第一、第二特征层进行堆叠处理，并将堆叠处理输入到U-Net网络，U-Net网络得到包含有多个通道第三特征层；第三卷积器得到通道数据为1的减弱伪影的结果图Xout。
这里，可以理解的是，第一神经网络1的每层的结构为“卷积-激活-卷积-激活-降采样层”，第二神经网络2的层的结构为“反卷积-卷积-激活-卷积-激活”，这里反卷积用于不断恢复深层特征的空间信息。
这里，用来训练的loss定义如下，其中αse＝αrc＝αie＝1，可以使用Adam优化器训练模型，相应的参数设置为(β1,β2)＝(0.5,0.999)，训练中学习率设置为从0.0002开始每训练30轮后衰减原来的1/2，整个模型一共训练了201轮，每次迭代的数据量为2。
Ltotal＝αse||Sse-Sgt||1+(αrc||Xse-Xgt||1+αie||Xout-Xgt||1)e(1-M)
这里，该神经网络模型保留投影数据中受金属影响的部分投影数据，在成像几何下的投影先验知识Mp的前提下，通过深度学习网络在投影域恢复投影数据中的细节信息，同时在图像域结合金属伪影的CT图像Mp和金属二值图Xm进一步抑制二次金属伪影。
这里，CT图像表示线性衰减系数的空间分布，能体现人体的解剖结构信息。X(E)代表依赖于能量E的衰减系数图。根据Lambert-Beer法则，探测器接收到的理想投影数据S可表示为(记为公式1)：
S＝-ln∫η(E)e-P(X(E))dE，其中，η(E)代表E能量分数，Ρ代表投影算子。
当金属出现在成像视野范围内时，X(E)随着能量E的不同变化剧烈，可以分解成两部分(记为公式2)：
X(E)＝Xr+Xm(E)＝Xr+λm(E)ρmM，
其中，第一部分为随着能量变化缓慢的人体组织Xr，第二部分为金属植入物的衰减系数图Xm(E)，Xm(E)可以进一步分解为金属二值图像Xm，金属密度ρm和随着能量剧烈变化的金属衰减系数λm(E)，根据投影算子Ρ的线性特性，金属的投影可以写成(记为公式3)：
P(X(E))＝λm(E)ρmP(M)＝λm(E)ρmMp
将公式3式带入公式1式后，可以得到：
其中，等式右边第一项是来自于Xr的投影数据P(Xr)，而第二项则给图像中引入金属伪影，在已知金属类型和能谱密度的情况下，第二项可以写为f(Mp)。因此，利用Mp这个先验知识，可以使投影域的伪影去除更加有效，同时最大程度的保留来自于有效区域的投影信息的细节。
本发明实施例二提供了一种神经网络的训练装置，包括以下模块：
数据获取模块，用于获取投影数据Sma、二值图Xm和Sma对应的金属投影数据Mp，所述投影数据Sma包含金属伪影，在所述二值图Xm中，值为1的像素值用于表征金属的形状；
神经网络创建模块，用于创建神经网络模型，所述神经网络模型包括第一神经网络、快速滤波反投影算子和第二神经网络；第一神经网络用于接收Sma和Mp，减弱Sma上的金属伪影并得到投影Sse；所述快速滤波反投影算子用于对投影Sse进行重建处理并得到图像Xse，以及对投影Sma进行重建处理并得到图像Xma，并且滤波反投影算子是可微的，能够将第二神经网络的训练中的梯度前向传播至第一神经网络；第二神经网络用于在图像域上，结合金属二值图Xm对图像Xse和Xma进行减弱二次金属伪影的处理；
训练模块，用于基于若干投影数据Sma、以及每个Sma对应的金属投影数据Mp，对所述神经网络模型进行有监督的训练。
本实施例还提供了一种电子设备，包括：
存储器，用于存储可执行指令；
处理器，用于执行所述存储器中存储的可执行指令时，实现实施例一中的训练方法。
本实施例还提供了一种计算机可读存储介质，存储有可执行指令，用于引起处理器执行时，实现实施例一中的训练方法。
本发明实施例三提供了一种消除金属伪影的方法，包括以下步骤：
获取包含有金属伪影的CT图像Xma，通过阈值分割从Xma得到金属的二值图像Xm，将Xma和Xm通过正投影算子得到含金属伪影的投影数据Sma和金属投影数据Mp，并将Sma、Mp和Xm输入权利要求1至5任一项所得到的神经网络模型，并得到减弱金属伪影的CT图像Xout。
本发明实施例四提供了一种消除金属伪影的装置，包括以下模块：
处理模块，用于获取包含有金属伪影的CT图像Xma，通过阈值分割从Xma得到金属的二值图像Xm，将Xma和Xm通过正投影算子得到含金属伪影的投影数据Sma和金属投影数据Mp，并将Sma、Mp和Xm输入权利要求1至5任一项所得到的神经网络模型，并得到减弱金属伪影的CT图像Xout。
这里，图4展示该神经网络模型的去除金属伪影的效果。
应当理解，虽然本说明书按照实施方式加以描述，但并非每个实施方式仅包含一个独立的技术方案，说明书的这种叙述方式仅仅是为清楚起见，本领域技术人员应当将说明书作为一个整体，各实施方式中的技术方案也可以经适当组合，形成本领域技术人员可以理解的其他实施方式。
上文所列出的一系列的详细说明仅仅是针对本发明的可行性实施方式的具体说明，它们并非用以限制本发明的保护范围，凡未脱离本发明技艺精神所作的等效实施方式或变更均应包含在本发明的保护范围之内。
Claims (10)
1.一种神经网络的训练方法，其特征在于，包括以下步骤：
获取投影数据Sma、二值图Xm和Sma对应的金属投影数据Mp，所述投影数据Sma包含金属伪影，在所述二值图Xm中，值为1的像素值用于表征金属的形状；
创建神经网络模型，所述神经网络模型包括第一神经网络、快速滤波反投影算子和第二神经网络；第一神经网络用于接收Sma和Mp，减弱Sma上的金属伪影并得到投影Sse；所述快速滤波反投影算子用于对投影Sse进行重建处理并得到图像Xse，以及对投影Sma进行重建处理并得到图像Xma，并且滤波反投影算子是可微的，能够将第二神经网络的训练中的梯度前向传播至第一神经网络；第二神经网络用于在图像域上，结合金属二值图Xm对图像Xse和Xma进行减弱二次金属伪影的处理；
基于若干投影数据Sma、以及每个Sma对应的金属投影数据Mp，对所述神经网络模型进行有监督的训练。
2.根据权利要求1所述的训练方法，其特征在于，所述第一神经网络包括：
第一编码器、第二编码器和第一解码器；
第一编码器接收投影数据Sma，第一编码器包含有M个排成一列的编码层，其中，M为自然数；每个编码层均包含有卷积模块、激活模块和池化模块，编码层的输入经过卷积模块提取特征信息，再经过激活模块进行非线性变化，最后通过池化模块压缩高维的特征并输出所述特征；
第二编码器接收金属投影数据Mp，第二编码器包含有N个排成一列的处理层，其中，N为自然数且N≤M，每个处理层都唯一对应到编码层；每个处理层均包含有池化模块，处理层的输入经过池化模块的平均池化处理、并得到池化结果；所述池化结果与对应的编码层的输出进行连接处理、并得到连接结果，所述连接结果作为对应的编码层的下一个编码层的输入；
第一解码器包含有M个排成一列的解码层，解码层包含有卷积模块、激活模块和池化模块，编码层的输入经过卷积模块提取特征信息，再经过激活模块进行非线性变化，最后通过池化模块得到更高维的特征；第1个解码层接收第M个编码层的输出，第i个解码器接收第M-i+1编码层的输出和第i-1个解码层的输出的连接处理的结果；
第一神经网络最后输出消除金属伪影的投影数据Sse。
3.根据权利要求1所述的训练方法，其特征在于，所述“获取投影数据Sma”具体包括：
获取不包含伪影的CT图Xgt；
将Xgt分割为软组织图Xw以及骨组织图Xb，对Xw和Xb进行处理并分别得到X’w和X’b，所述处理具体包括：在确定所述二值图Xm中的第一像素的值为1时，则将Xw中的与第一像素对应的第二像素的值置为0，将Xb中的与第一像素对应的第三像素的值置为0；
基于正投影算子分别对X’w、X’b和Xm分别进行处理，并分别得到若干投影数据Pw、若干投影数据Pb和若干投影数据Pm，其中，不同的Pw表征软组织在不同能级下的投影数据，不同的Pb表征骨组织在不同能级下的投影数据，不同的Pm表征金属在不同能级下的投影数据；
基于不同组织在不同能级下的衰减密度曲线、对Pw、Pb和Pm进行能谱密度归一化处理，从而得到不含其它噪声的投影数据S* ma；
对S* ma叠加泊松噪声和散射噪声，之后，进行射束硬化校正处理，从而得到投影数据Sma。
4.根据权利要求3所述的训练方法，其特征在于，所述“将Xgt分割为软组织图Xw以及骨组织图Xb”具体包括：基于软阈值法，将Xgt分割为软组织图Xw以及骨组织图Xb。
5.根据权利要求3所述的训练方法，其特征在于，
第二神经网络包括：第一卷积器、第二卷积器、第三卷积器和U-Net网络，首先对图像Xma和二值图Xm进行堆叠处理并得到图像X’ma，对图像Xse和二值图Xm进行堆叠处理从而分别得到图像X’se；之后，第一卷积器获取图像X’ma的融合有金属位置的第一特征层，第二卷积器获取图像X’se的融合有金属位置的第二特征层；之后，对第一、第二特征层进行堆叠处理，并将堆叠处理输入到U-Net网络，U-Net网络得到包含有多个通道第三特征层；第三卷积器得到通道数据为1的减弱伪影的结果图Xout。
6.一种神经网络的训练装置，其特征在于，包括以下模块：
数据获取模块，用于获取投影数据Sma、二值图Xm和Sma对应的金属投影数据Mp，所述投影数据Sma包含金属伪影，在所述二值图Xm中，值为1的像素值用于表征金属的形状；
神经网络创建模块，用于创建神经网络模型，所述神经网络模型包括第一神经网络、快速滤波反投影算子和第二神经网络；第一神经网络用于接收Sma和Mp，减弱Sma上的金属伪影并得到投影Sse；所述快速滤波反投影算子用于对投影Sse进行重建处理并得到图像Xse，以及对投影Sma进行重建处理并得到图像Xma，并且滤波反投影算子是可微的，能够将第二神经网络的训练中的梯度前向传播至第一神经网络；第二神经网络用于在图像域上，结合金属二值图Xm对图像Xse和Xma进行减弱二次金属伪影的处理；
训练模块，用于基于若干投影数据Sma、以及每个Sma对应的金属投影数据Mp，对所述神经网络模型进行有监督的训练。
7.一种电子设备，其特征在于，包括：
存储器，用于存储可执行指令；
处理器，用于执行所述存储器中存储的可执行指令时，实现权利要求1至5任一项所述的训练方法。
8.一种计算机可读存储介质，其特征在于，存储有可执行指令，用于引起处理器执行时，实现权利要求1至5任一项所述的训练方法。
9.一种消除金属伪影的方法，其特征在于，包括以下步骤：
获取包含有金属伪影的CT图像Xma，通过阈值分割从Xma得到金属的二值图像Xm，将Xma和Xm通过正投影算子得到含金属伪影的投影数据Sma和金属投影数据Mp，并将Sma、Mp和Xm输入权利要求1至5任一项所得到的神经网络模型，并得到减弱金属伪影的CT图像Xout。
10.一种消除金属伪影的装置，其特征在于，包括以下模块：
处理模块，用于获取包含有金属伪影的CT图像Xma，通过阈值分割从Xma得到金属的二值图像Xm，将Xma和Xm通过正投影算子得到含金属伪影的投影数据Sma和金属投影数据Mp，并将Sma、Mp和Xm输入权利要求1至5任一项所得到的神经网络模型，并得到减弱金属伪影的CT图像Xout。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010971211.5A CN112085809B (zh) | 2020-09-16 | 2020-09-16 | 神经网络的训练方法、消除金属伪影的方法及其装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010971211.5A CN112085809B (zh) | 2020-09-16 | 2020-09-16 | 神经网络的训练方法、消除金属伪影的方法及其装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112085809A true CN112085809A (zh) | 2020-12-15 |
| CN112085809B CN112085809B (zh) | 2024-04-16 |
Family
ID=73736440
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010971211.5A Active CN112085809B (zh) | 2020-09-16 | 2020-09-16 | 神经网络的训练方法、消除金属伪影的方法及其装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112085809B (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113554563A (zh) * | 2021-07-23 | 2021-10-26 | 上海友脉科技有限责任公司 | 一种医学图像处理方法、介质及电子设备 |
| CN113570586A (zh) * | 2021-08-02 | 2021-10-29 | 苏州工业园区智在天下科技有限公司 | 神经网络系统的创建、处理ct图像的方法及其装置 |
| CN114241070A (zh) * | 2021-12-01 | 2022-03-25 | 北京长木谷医疗科技有限公司 | Ct图像去金属伪影以及模型训练的方法和装置 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107871332A (zh) * | 2017-11-09 | 2018-04-03 | 南京邮电大学 | 一种基于残差学习的ct稀疏重建伪影校正方法及系统 |
| CN110458762A (zh) * | 2019-07-02 | 2019-11-15 | 山东科技大学 | 一种基于可调双因子的ct图像射束硬化伪影校正系统 |
| CN110916708A (zh) * | 2019-12-26 | 2020-03-27 | 南京安科医疗科技有限公司 | 一种ct扫描投影数据伪影校正方法、ct图像重建方法 |
-
2020
- 2020-09-16 CN CN202010971211.5A patent/CN112085809B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107871332A (zh) * | 2017-11-09 | 2018-04-03 | 南京邮电大学 | 一种基于残差学习的ct稀疏重建伪影校正方法及系统 |
| CN110458762A (zh) * | 2019-07-02 | 2019-11-15 | 山东科技大学 | 一种基于可调双因子的ct图像射束硬化伪影校正系统 |
| CN110916708A (zh) * | 2019-12-26 | 2020-03-27 | 南京安科医疗科技有限公司 | 一种ct扫描投影数据伪影校正方法、ct图像重建方法 |
Non-Patent Citations (3)
| Title |
|---|
| CHENGTAO PENG 等: "A Cross-domain Metal Trace Restoring Network for Reducing X-ray CT Metal Artifacts", IEEE TRANSACTIONS ON MEDICAL IMAGING, 29 June 2020 (2020-06-29) * |
| CHENGTAO PENG 等: "A Cross-domain Metal Trace Restoring Network for Reducing X-ray CT Metal Artifacts", IEEE TRANSACTIONS ON MEDICAL IMAGING, 31 December 2020 (2020-12-31) * |
| 马建华;陈武凡;黄静;杨迪;毕一鸣;: "基于最大互信息量熵差分割的CT金属伪影消除", 电子学报, no. 08, 15 August 2009 (2009-08-15) * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113554563A (zh) * | 2021-07-23 | 2021-10-26 | 上海友脉科技有限责任公司 | 一种医学图像处理方法、介质及电子设备 |
| CN113554563B (zh) * | 2021-07-23 | 2024-05-14 | 上海友脉科技有限责任公司 | 一种医学图像处理方法、介质及电子设备 |
| CN113570586A (zh) * | 2021-08-02 | 2021-10-29 | 苏州工业园区智在天下科技有限公司 | 神经网络系统的创建、处理ct图像的方法及其装置 |
| CN114241070A (zh) * | 2021-12-01 | 2022-03-25 | 北京长木谷医疗科技有限公司 | Ct图像去金属伪影以及模型训练的方法和装置 |
| CN114241070B (zh) * | 2021-12-01 | 2022-09-16 | 北京长木谷医疗科技有限公司 | Ct图像去金属伪影以及模型训练的方法和装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112085809B (zh) | 2024-04-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Zhang et al. | Convolutional neural network based metal artifact reduction in x-ray computed tomography | |
| CN112085809B (zh) | 神经网络的训练方法、消除金属伪影的方法及其装置 | |
| Karimi et al. | Segmentation of artifacts and anatomy in CT metal artifact reduction | |
| US7362845B2 (en) | Method and apparatus of global de-noising for cone beam and fan beam CT imaging | |
| US7831097B2 (en) | System and method for image reconstruction | |
| Wang et al. | Metal artifact reduction in CT using fusion based prior image | |
| Chang et al. | Prior-guided metal artifact reduction for iterative X-ray computed tomography | |
| Park et al. | Metal artifact reduction in CT by identifying missing data hidden in metals | |
| Peng et al. | A cross-domain metal trace restoring network for reducing X-ray CT metal artifacts | |
| US10657679B2 (en) | Multi-energy (spectral) image data processing | |
| CN114387359A (zh) | 一种三维x射线低剂量成像方法及装置 | |
| Karimi et al. | Metal artifact reduction for CT-based luggage screening | |
| CN111223156A (zh) | 一种用于牙科锥束ct系统的金属伪影消除方法 | |
| CN106233327A (zh) | 去噪图像数据中对低对比度结构的恢复 | |
| US11049295B2 (en) | Detection and/or correction of residual iodine artifacts in spectral computed tomography (CT) imaging | |
| Chen et al. | Novel method for reducing high-attenuation object artifacts in CT reconstructions | |
| Wang et al. | Hybrid pre-log and post-log image reconstruction for computed tomography | |
| Pua et al. | A pseudo-discrete algebraic reconstruction technique (PDART) prior image-based suppression of high density artifacts in computed tomography | |
| CN111080740A (zh) | 一种图像校正方法、装置、设备及介质 | |
| Mehranian et al. | Metal artifact reduction in CT-based attenuation correction of PET using sobolev sinogram restoration | |
| Wang et al. | A new approach for reducing beam hardening artifacts in polychromatic X-ray computed tomography using more accurate prior image | |
| Kinahan et al. | Quantitative attenuation correction for PET/CT using iterative reconstruction of low-dose dual-energy CT | |
| KR102590461B1 (ko) | 이중 에너지 기법 기반 다중 물질 분리 방법 및 장치 | |
| CN111091516B (zh) | 一种基于人工智能的抗散射光栅方法及装置 | |
| Safdari et al. | A new method for metal artifact reduction in CT scan images |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |