一种用于电阻抗肺部成像图像的优化方法与流程

1.本发明涉及电阻抗成像，具体涉及一种用于电阻抗肺部成像图像的优化方法。


背景技术：

2.电阻抗层析成像(electrical impedance tomography简称eit)是一种非侵入式的成像技术，它通过场域的边界测量值重建出场域内部的电导率分布，相比其它成像技术如计算机断层扫描(computed tomography,简称ct)和核磁共振成像(magnetic resonance imaging简称mri)，eit有着无辐射、结构简单、速度快等优点，在医学临床监测、多相流流型监测以及地址勘探等领域有着广阔的应用前景。
3.然而，由于eit固有的“软场”特性，其重建图像往往分辨率较低，并且受测量误差和模型误差的影响较大，这限制了eit技术的进一步推广。eit的“软场”特性在数学上表现为逆问题的非线性和病态性，传统重建算法通过线性化的方法进行近似和减小模型误差的影响，并引入正则化项来减弱逆问题的病态性。最近，一些直接非线性的重建算法被提出，比如level-set、d-bar等，而基于神经网络的重建算法也属于此类。
4.近几年来，伴随着深度学习在计算机视觉和图像处理领域的迅速发展，众多深度网络特别是卷积神经网络(convolutional neural network，简称cnn)也纷纷被引入到医学成像领域中。这些数据驱动的新型成像技术被视为医学成像领域继解析法、迭代法后的第三个发展阶段，它们利用大量数据训练好的网络，来补充或代替传统成像算法中人为定义的模型，将大数据内在信息自动融入到成像过程中。其中，一种较为简便易行的办法是在图像域设计并训练一个对粗糙的重建图像加以处理的cnn，使cnn能够学习到医学成像过程中产生的噪声以及伪影等，并在输出图像中移除这些不利因素，这可以被视为一个图像增强或者图像变换过程。
5.然而，在进行图像重建时需要进行大量的计算，并且非常容易受到测量噪声的影响，进而导致图像重建的效率，以及准确度都比较低，无法满足电阻抗成像应用于肺部时对效率、准确度方面的要求。


技术实现要素：

6.(一)解决的技术问题
7.针对现有技术所存在的上述缺点，本发明提供了一种用于电阻抗肺部成像图像的优化方法，能够有效克服现有技术所存在的电阻抗成像应用于肺部时图像重建的效率及准确度都比较低的缺陷。
8.(二)技术方案
9.为实现以上目的，本发明通过以下技术方案予以实现：
10.一种用于电阻抗肺部成像图像的优化方法，包括以下步骤：
11.s1、根据电流信息和电压信息获取电导率信息，并获得电阻抗成像数据，通过迭代得到用于反演新的电阻抗成像数据的下降方向；
12.s2、基于电导率分布初始模型、电阻抗成像正向模型、实测电阻抗成像数据和下降方向进行图像重建；
13.s3、将重建图像输入预设神经网络，输出重建区域及对应置信度；
14.s4、根据置信度获取最优的重建区域，并作为最终电阻抗成像图像。
15.优选地，s1中根据电流信息和电压信息获取电导率信息，并获得电阻抗成像数据，通过迭代得到用于反演新的电阻抗成像数据的下降方向，包括：
16.通过回归方法得到电导率分布、电阻抗成像数据分别与每一步迭代的电导率分布、电阻抗成像数据的差值，从而得到新的电阻抗成像数据的下降方向。
17.优选地，所述电导率信息、电阻抗成像数据为真实数据，或者为所述电阻抗成像正向模型的仿真生成数据。
18.优选地，所述电阻抗成像正向模型为二维或三维，采用包括特定电极数目和电极排布的电极模型。
19.优选地，s2中基于电导率分布初始模型、电阻抗成像正向模型、实测电阻抗成像数据和下降方向进行图像重建，包括：
20.采用下式进行图像重建：
21.σn＝σ
n-1
k(r(σ
n-1
)-d)
22.其中，n为迭代步数，σn为第n步的电导率分布，k为下降方向，r(σ)为电阻抗成像正向模型，d为实测电阻抗成像数据。
23.优选地，s3中将重建图像输入预设神经网络，输出重建区域及对应置信度，包括：
24.获取预设神经网络中的softmax函数，将输出中的最大值作为置信度。
25.优选地，s3中将重建图像输入预设神经网络，输出重建区域及对应置信度，包括：
26.在预设神经网络中设置一条置信度路径，并根据置信度路径设置损失函数l’；
27.将损失函数l’与原始损失函数lo结合，获得总损失函数l；
28.采用总损失函数l对神经网络进行训练，并在训练后，获取损失函数l’的值作为置信度；
29.其中，当神经网络输出的重建区域正确时，损失函数l’的值近于0；当神经网络输出的重建区域错误时，损失函数l’的值近于1。
30.优选地，在预设神经网络中设置一条置信度路径后，神经网络输出的表达式为p’＝c*p (1-c)*y，其中y为目标分类结果，c为神经网络输出的置信度。
31.优选地，所述原始损失函数lo采用交叉熵损失函数计算，所述损失函数l’的计算公式为l’＝sigmoid(-logc)，其中c为神经网络输出的置信度。
32.优选地，s4中根据置信度获取最优的重建区域，并作为最终电阻抗成像图像，包括：
33.对置信度进行降序排列，获取前m个置信度对应的重建区域，将m个重建区域之间的交叉区域作为重叠区域，将包含最多重叠区域的重建区域作为最终电阻抗成像图像。
34.(三)有益效果
35.与现有技术相比，本发明所提供的一种用于电阻抗肺部成像图像的优化方法，具有以下有益效果：
36.1)根据电流信息和电压信息获取电导率信息，并获得电阻抗成像数据，通过迭代
得到用于反演新的电阻抗成像数据的下降方向，并基于电导率分布初始模型、电阻抗成像正向模型、实测电阻抗成像数据和下降方向进行图像重建，省去了在图像重建时大量繁杂的计算过程，有效提升了图像重建的效率；
37.2)将重建图像输入预设神经网络，输出重建区域及对应置信度，根据置信度获取最优的重建区域，并作为最终电阻抗成像图像，能够有效提升电阻抗成像应用于肺部时的准确度，确保电阻抗成像的精准性。
附图说明
38.为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍。显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
39.图1为本发明的流程示意图；
40.图2为本发明中电阻抗肺部成像的优化流程示意图。
具体实施方式
41.为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述。显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
42.一种用于电阻抗肺部成像图像的优化方法，如图1和图2所示，
①
根据电流信息和电压信息获取电导率信息，并获得电阻抗成像数据，通过迭代得到用于反演新的电阻抗成像数据的下降方向，具体包括：
43.通过回归方法得到电导率分布、电阻抗成像数据分别与每一步迭代的电导率分布、电阻抗成像数据的差值，从而得到新的电阻抗成像数据的下降方向。
44.本技术技术方案中，电导率信息、电阻抗成像数据可以为真实数据，也可以为电阻抗成像正向模型的仿真生成数据。其中，电阻抗成像正向模型为二维或三维，采用包括特定电极数目和电极排布的电极模型。
45.②
基于电导率分布初始模型、电阻抗成像正向模型、实测电阻抗成像数据和下降方向进行图像重建，具体包括：
46.采用下式进行图像重建：
47.σn＝σ
n-1
k(r(σ
n-1
)-d)
48.其中，n为迭代步数，σn为第n步的电导率分布，k为下降方向，r(σ)为电阻抗成像正向模型，d为实测电阻抗成像数据。
49.本技术技术方案中，根据电流信息和电压信息获取电导率信息，并获得电阻抗成像数据，通过迭代得到用于反演新的电阻抗成像数据的下降方向，并基于电导率分布初始模型、电阻抗成像正向模型、实测电阻抗成像数据和下降方向进行图像重建，省去了在图像重建时大量繁杂的计算过程(例如采用梯度下降法时在每一步的迭代过程中都需要计算雅克比矩阵)，有效提升了图像重建的效率。
50.③
将重建图像输入预设神经网络，输出重建区域及对应置信度。在本技术技术方案中，可以采用两种方案输出重建区域及对应置信度：
51.⑴
获取预设神经网络中的softmax函数，将输出中的最大值作为置信度；
52.⑵
在预设神经网络中设置一条置信度路径，并根据置信度路径设置损失函数l’，将获取损失函数l’的值作为置信度。
53.上述第二种方法具体包括以下步骤：
54.在预设神经网络中设置一条置信度路径，并根据置信度路径设置损失函数l’；
55.将损失函数l’与原始损失函数lo结合，获得总损失函数l；
56.采用总损失函数l对神经网络进行训练，并在训练后，获取损失函数l’的值作为置信度。
57.其中，当神经网络输出的重建区域正确时，损失函数l’的值近于0；当神经网络输出的重建区域错误时，损失函数l’的值近于1。
58.在预设神经网络中设置一条置信度路径后，神经网络输出的表达式为p’＝c*p (1-c)*y。其中，y为目标分类结果，c为神经网络输出的置信度。
59.原始损失函数lo采用交叉熵损失函数计算，损失函数l’的计算公式为l’＝sigmoid(-logc)。其中，c为神经网络输出的置信度，当置信度趋近于1(神经网络输出的重建区域正确)时，损失函数l’的值近于0；当置信趋近于0(神经网络输出的重建区域错误)时，损失函数l’的值近于1。
60.④
根据置信度获取最优的重建区域，并作为最终电阻抗成像图像，具体包括：
61.对置信度进行降序排列，获取前m个置信度对应的重建区域，将m个重建区域之间的交叉区域作为重叠区域，将包含最多重叠区域的重建区域作为最终电阻抗成像图像。
62.本技术技术方案中，也可以直接将最高置信度对应的重建区域作为最终电阻抗成像图像，但是引入“重叠区域”的概念，使得作为最终电阻抗成像图像的重建区域包含最多“重叠区域”的部分，从而能够有效提升电阻抗成像应用于肺部时的准确度，确保电阻抗成像的精准性。
63.以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不会使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。