CT图像金属伪影去除方法、装置及计算机可读存储介质

摘要

本发明公开了一种CT图像金属伪影去除方法、装置及计算机可读存储介质，包括步骤：收集包含金属伪影图像和不含金属伪影图像，处理得到对应的投影数据；对生成的两种投影数据分别进行重建，生成包含金属伪影图像和无伪影图像，进一步处理得到初始伪影抑制图像和数据完整性map；设计神经网络和损失函数，神经网络输出去除伪影后的图像等。本发明的CT图像金属伪影去除方法，利用了金属伪影和周围结构的相关性，以及与组织成分的关系，克服了现有的图像域技术不能把金属扫描数据中的数据完整性信息融入到图像域的处理的过程中的问题，提高了伪影校正的准确性和CT图像的质量。

权利要求

1.一种CT图像金属伪影去除方法，其特征在于，包括步骤：S1、创建数据集：收集病人的不同CT图像，包含金属伪影图像和不含金属伪影图像，处理得到对应的包含金属伪影的投影数据和不含金属伪影的投影数据；

S2、数据预处理：对步骤S1中生成的两种投影数据分别进行重建，生成包含金属伪影CT图像 和无伪影的CT图像 ；

对包含金属伪影CT图像 进行分割处理得到仅包含金属的图像，对金属图像进行正投得到金属正投影图像，根据金属伪影的强弱受组织结构以及距离组织结构的距离影响的关系，对金属正投影图像进行处理，得到更正的投影，然后进行重建得到初始伪影抑制图像 ；

对正投得到的金属投影图像进行反投处理，得到数据完整性map；

S3、设计神经网络：神经网络的输入为包含金属伪影CT图像 、初始伪影抑制图像和数据完整性map，神经网络的输出为无伪影图像 ；步骤S2得到的数据用于训练神经网络；

构建的神经网络采用U-net网络或根据U-net网络改进的网络，其损失函数Loss为：其中， 为欧式距离，表示像素点距离最近的金属标记物的距离， 为金属的中心位置； 为组织权重，u表示均值参数， 表示方差；

S4、用于CT图像金属伪影消除：将待处理的CT图像进行步骤S1和S2的处理，得到对应的伪影数据、初始伪影抑制图像 和数据完整性map，输入到步骤S3训练得到的神经网络中，输出经过消除伪影的图像。

2.根据权利要求1所述的一种CT图像金属伪影去除方法，其特征在于，步骤S1顺序执行以下步骤：收集病人的不同CT图像，包括包含金属伪影的第一图像和不包含金属伪影的第二图像；

对第一图像进行分割处理，提取出其中的金属部分作为模拟金属植入物；

采用数值仿真的方法对第二图像进行处理，得到不含金属伪影的第一投影数据；

将得到的模拟金属植入物插入第二图像，通过数值仿真的方法得到包含金属伪影的第二投影数据。

3.根据权利要求1所述的一种CT图像金属伪影去除方法，其特征在于，步骤S2顺序执行以下步骤：对包含金属伪影CT图像 进行分割处理得到仅包含金属的图像；

对金属图像进行正投得到金属正投影图像，按照金属正投影图像的轨迹对被金属植入物影响的投影数据进行线性插值；

用周围组织代替金属部分，得到插值更正的投影，然后进行重建得到初始伪影抑制图像 。

4.一种CT图像金属伪影去除装置，其特征在于，包括处理器、存储器及通信总线；存储器上存储有可被所述处理器执行的计算机可读程序，通信总线用于实现处理器和存储器之间的连接通信；

所述处理器执行所述计算机可读程序时实现如权利要求1-3任意一项所述的CT图像金属伪影去除方法中的步骤。

5.一种计算机可读存储介质，其特征在于：所述计算机可读存储介质存储有至少一个可被处理器执行的指令，其中至少一个指令被处理器执行时，用于执行如权利要求1至3中任一所述CT图像金属伪影去除方法。

说明书

CT图像金属伪影去除方法、装置及计算机可读存储介质

技术领域

[0001] 本发明涉及医学成像技术领域，具体涉及一种CT图像金属伪影去除方法、装置及计算机可读存储介质。

背景技术

[0002] 计算机断层扫描（CT）能够无损检测到人体内部的组织器官结构，因此广泛的应用于临床诊断。受到金属植入物的影响，CT图像中会出现不同程度的金属伪影，表现为条状伪影和亮暗带，严重影响图像质量，阻碍医生精确诊断。因此近40年出现了大量金属伪影校正（metal artifacts reduction，MAR）方法去除金属伪影。如数据采集改进、基于物理的校正、投影数据校正、迭代方法和图像后处理技术，以及新兴的基于深度学习的方法。

[0003] 数据采集改进方法，需要保证扫描区域内没有金属，如拿掉耳环，项链等，但在大多数情况下不能满足，如金属植入物在人体内，并调节扫描参数等抑制金属伪影但是该方法很难获得足够好的图像质量。基于物理的校正方法，校正射束硬化和光子饥饿等，由于X光成像系统的复杂性，很难对成像过程做出精确的模拟，通常的方法都是基于理想的近似，因此方法误差较大，无法取得令人满意的结果。基于投影数据校正的方法容易引入新的图像伪影，迭代方法虽然能提高图像质量，但是需要较大的计算量，和复杂的物理模型。另外迭代算法，通常的模型针对的是数据中的统计噪声，因此更加适合处理数据中的噪声，但是，对于存在金属伪影的数据中，更需要的是对于确实信号的恢复。因此，针对金属伪影的问题，引入了图像域后处理方法，可以抑制轻微的伪影，不能处理严重的情况。但是，这些图像域的方法给予的多是图像形态学、甚至是深度学习的方法，由于金属伪影的复杂性，很难仅仅从图像上达到满意的效果。

[0004] 正是基于上述几个原因，考虑到金属伪影的严重程度，不仅与植入物的材料和形状相关，而且与植入位置以及组织背景密切相关，不同的位置和组织产生的伪影严重程度不一样。如图1所示为金属伪影不同表现的示意图。

[0005] 文献【1】（即Y. Zhang and H. Yu, “Convolutional neural network based metal artifact reduction in x-ray computed tomography,” arXiv preprint arXiv:1709.01581, 2017.）中提出的金属伪影的抑制技术如图2所示，该方法是基于图像域的处理，虽然增加了输入信息，比如输入图像包括三种传统处理方法得到的CT图像，三个图像输入到网络以后，还需要经过后处理和正投，重建等复杂的过程才能得到临床上可以接受的图像。

[0006] 现有图像域的方法，大多是从原始图像或者插值投影的图像来进行图像处理的。但是金属伪影和周围结构的相关性，还有组织的成分关系很大，比如软组织周围的金属伪影就比较少，而骨组织和金属周围的伪影就会更加严重。其次现有的图像域的技术不能把金属扫描数据中的数据完整性信息融入到图像域的处理的过程中，这个数据的完整性也是对金属伪影处理的重要意义。比如数据完整的区域不存在金属伪影，而数据缺失的区域金属伪影会比较强。而且金属伪影的方向和这个数据的完整性也有着重要的关系。如图3所示金属导致的图像数据缺失，从左到右依次为金属植入物，植入物所在投影域的图像，以及数据完整性map。在map中数值越高表示完整性越差，伪影越严重。

发明内容

[0007] 技术目的：针对现有技术中的不足，本发明提出一种基于人体组织结构和数据完整性的深度学习CT图像金属伪影去除方法、装置及计算机可读存储介质，用以提高CT图像质量。

[0008] 技术方案：为实现上述技术目的，本发明采用了如下技术方案：

[0009] 一种CT图像金属伪影去除方法，其特征在于，包括步骤：

[0010] S1、创建数据集：收集病人的不同CT图像，包含金属伪影图像和不含金属伪影图像，处理得到对应的包含金属伪影的投影数据和不含金属伪影的投影数据；

[0011] S2、数据预处理：对步骤S1中生成的两种投影数据分别进行重建，生成包含金属伪影CT图像 和无伪影的CT图像 ；

[0012] 对包含金属伪影CT图像 进行分割处理得到仅包含金属的图像，对金属图像进行正投得到金属正投影图像，根据金属伪影的强弱受组织结构以及距离组织结构的距离影响的关系，对金属正投影图像进行处理，得到更正的投影，然后进行重建得到初始伪影抑制图像 ；

[0013] 对正投得到的金属投影图像进行反投处理，得到数据完整性map；

[0014] S3、设计神经网络：神经网络的输入为包含金属伪影CT图像 、初始伪影抑制图像 和数据完整性map，神经网络的输出为无伪影图像 ；步骤S2得到的数据用于训练神经网络；

[0015] S4、用于CT图像金属伪影消除：将待处理的CT图像进行步骤S1和S2的处理，得到对应的伪影数据、初始伪影抑制图像 和数据完整性map，输入到步骤S3训练得到的神经网络中，输出经过消除伪影的图像。

[0016] 优选地，步骤S1顺序执行以下步骤：

[0017] 收集病人的不同CT图像，包括包含金属伪影的第一图像和不包含金属伪影的第二图像；

[0018] 对第一图像进行分割处理，提取出其中的金属部分作为模拟金属植入物；

[0019] 采用数值仿真的方法对第二图像进行处理，得到不含金属伪影的第一投影数据；

[0020] 将得到的模拟金属植入物插入第二图像，通过数值仿真的方法得到包含金属伪影的第二投影数据。

[0021] 优选地，步骤S2顺序执行以下步骤：

[0022] 对包含金属伪影CT图像 进行分割处理得到仅包含金属的图像；

[0023] 对金属图像进行正投得到金属正投影图像，按照金属正投影图像的轨迹对被金属植入物影响的投影数据进行线性插值；

[0024] 用周围组织代替金属部分，得到插值更正的投影，然后进行重建得到初始伪影抑制图像 。

[0025] 优选地，步骤S3中构建的神经网络采用U-net网络或根据U-net网络改进的网络，其损失函数Loss为：

[0026]

[0027]

[0028]

[0029] 其中， 为欧式距离，表示像素点距离最近的金属标记物的距离，为金属的中心位置； 为组织权重，u表示均值参数， 表示方差。

[0030] 本发明还公开了一种CT图像金属伪影去除装置，其特征在于，包括处理器、存储器及通信总线；存储器上存储有可被所述处理器执行的计算机可读程序，通信总线用于实现处理器和存储器之间的连接通信；

[0031] 所述处理器执行所述计算机可读程序时实现所述的CT图像金属伪影去除方法中的步骤。

[0032] 本发明还公开了一种计算机可读存储介质，其特征在于：所述计算机可读存储介质存储有至少一个可被处理器执行的指令，其中至少一个指令被处理器执行时，用于执行所述CT图像金属伪影去除方法。

[0033] 技术效果：由于采用了上述技术方案，本发明具有如下技术效果：

[0034] 本发明的CT图像金属伪影去除方法，利用了金属伪影和周围结构的相关性，以及与组织成分的关系，克服了现有的图像域技术不能把金属扫描数据中的数据完整性信息融入到图像域的处理的过程中的问题，提高了伪影校正的准确性和CT图像的质量。

附图说明

[0035] 图1为金属伪影不同表现的示意图；

[0036] 图2为参考文献1的处理过程的示意图；

[0037] 图3为金属导致的图像数据缺失的示意图，从左到右依次为金属植入物的图像、植入物所在投影域的图像、数据完整性map；

[0038] 图4为本发明的金属伪影去除方法的示意图一；

[0039] 图5为本发明的金属伪影去除方法的示意图二；

[0040] 图6为本发明的金属伪影去除方法的流程图；

[0041] 图7为更正结果的示意图，从左到右依次为伪影图像、NMAR更正结果和采用本发明方法得到的图像。

具体实施方式

[0042] 本发明提出一种基于人体组织结构和数据完整性的深度学习金属伪影抑制方法，用于提高CT图像质量。本发明的金属伪影去除方法的示意图如图4至图6所示。本发明的目的在于通过深度学习的方法校正金属伪影，帮助医生快速精准诊断。具体步骤如下：

[0043] 1、创建数据集：包含金属伪影图像和不含金属伪影图像，以及其对应的投影图像，参考文献1。具体方法为收集病人不同CT图像，分为包含金属和不包含金属两部分；从包含金属植入物的CT图像中分割提取出其中的金属部分作为模拟植入物；通过数值仿真的方法对不包含金属伪影的图像处理，得到不含金属的投影数据；把模拟金属植入物插入不含金属伪影的CT图像，通过数值仿真的方法得到包含金属的投影数据。本发明用到的数值仿真方法可参考文献1中的方法，并对该方法进行修改。文献1中模拟了射束硬化和泊松分布的特点，没有考虑物体散射情况。为了确保与真实情况相似，影域增加一个低频信号（例如原始信号的10%），模拟散射的贡献。更精确的散射模拟，可以通过Monte Carlo仿真来实现。在本发明中，对图像细化，分为空气、水、软组织、皮肤、骨组织和金属，采用多能量多物质线性积分的方法得到投影数据。

[0044] 2、数据预处理：对步骤1中生成的投影数据重建生成无伪影的CT图像 和包含金属伪影CT图像 ；对 进行分割处理得到仅包含金属的图像，对金属图像进行正投得到金属投影图像，对金属图像进行正投操作，探测器上数值非0的点即为金属正投影图像的轨迹。

[0045] 按照金属在投影图像的轨迹对被金属植入物影响的投影数据进行线性插值，用周围组织代替金属部分，得到插值更正的投影，具体插值方法可以通过行、列和角度方向上，没有被金属污染的投影数据来恢复被污染的像素。如：线性插值、双线性插值和三线性插值以及区域填充的方法等。

[0046] 然后进行重建得到初始伪影抑制图像 ，重建方法用FBP传统解析法。

[0047] 对正投得到的金属投影图像进行反投处理，得到数据完整性map，map中的值越大，表示该位置数据缺失越明显。

[0048] 3、设计网络结构和损失函数：网络结构如U-net及其变种等，输入为 、和map输出为无伪影图像 。金属伪影的强弱受组织结构以及距离组织结构的距离影响，损失函数如下：

[0049]

[0050]

[0051]

[0052] 其中， 为欧式距离，表示像素点距离最近的金属标记物的距离，为金属的中心位置； 为组织权重。例如：5000为金属HU阈值，大于该值
权重为0，例如高斯形状的权重，参数均值u=2000，方差 =200；把步骤2得到的数据输入到网络训练调整参数，直到网络收敛。

[0053] 4、将伪影数据和初始伪影抑制图像 和完整性map输入到网络，得到伪影校正结果。如图7所示为校正结果的示意图，从左到右依次为伪影图像、NMAR更正结果和采用本发明方法得到的图像。

[0054] 以上所述仅是本发明的优选实施方式，应当指出：对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。