一种基于模糊聚类的电阻抗层析成像优化算法的制作方法

本发明属于电学层析成像图像重建算法领域，具体涉及一种基于模糊聚类的电阻抗层析成像优化算法。

背景技术：

电学层析成像(electricaltomography,简称et)技术是一种重建被测敏感场域内部电学特性分布的技术，以其无辐射、非侵入性、响应迅速、成本低廉等优势在医学临床监护和工业测量等领域有着广阔的应用前景。在电阻层析成像(ert)中，利用阵列电极对被测物场施加电流激励，测量边界电压数据，进而重建被测场域内的电导率分布。以16电极的电阻层析成像传感器为例，采用相邻电流激励、相邻电压测量的方式。在每次激励下，除去激励电极，共有13个测量值。如：以1、2号电极作为激励电极，测量3、4，4、5，5、6......15、16电极对上的电压信号。共有16次激励，即有13×16＝208个测量值。为了兼顾时间分辨率及空间分辨率，一般将被测敏感场域划分为812个像素点。电阻层析成像图像重建问题即为求解这812个像素点的灰度值。传统的线性反投影算法(lbp)即为根据边界电压值和被测场域投影区的关系进行电压数据的反投，但由于硬件设备的限制，导致反投数据远远少于需要求解的灰度值，使得重建图像有较大的伪影，目标及背景的边界模糊，空间分辨率较差。

技术实现要素：

本发明的目的是提供一种可以减少伪影，具有更高的空间分辨率的电阻抗层析成像优化算法。本发明的技术方案如下：

一种基于模糊聚类的电阻抗层析成像优化算法，从传感器采集得到的边界电压值中提取有用信息，形成特征向量，进而进行聚类分析，方法如下：

(1)得到电阻抗层析成像的图像矩阵；

(2)从图像矩阵中提取代表每个像素信息的特征，特征一：均值向量，即为每个像素点所对应的各个边界电压值的均值；特征二：方差矩阵，在均值矩阵的基础上，计算每个像素点的方差，每个像素点的2维特征向量形成特征矩阵；

(3)将特征矩阵作为模糊聚类算法的输入数据集进行模糊聚类，进而得到每个像素点的隶属度值，表征每个像素点对每一类的隶属度大小；

(4)根据隶属度矩阵进行类别的归属判断。

本发明由于采取以上技术方案，其具有以下优点：

1、与传统的lbp算法相比，本发明算法减少了伪影，使目标与背景之间的界限更加清晰，提高了成像的空间分辨率。

2、聚类算法采用模糊聚类算法而非传统的硬划分算法，使得目标边缘更加平滑，贴近实际边缘。

附图说明

图1为本发明计算像素点方差时的示意图；

图2为本发明的算法流程图；

图3为本发明的实施例中原始被测物场模型图；

图4为对图3仿真模型用lbp算法进行图像重建的结果；

图5为对图3仿真模型用本发明的算法进行图像重建的结果。

具体实施方式

本发明的目的是克服现有图像重建算法空间分辨率较低的不足，提出了一种基于模糊聚类的电阻抗层析成像优化算法。模糊聚类的对象是一组数据的集合，通过提取数据的特征向量来对数据进行划分，使类内差别较小，类间差异较大。基于模糊聚类的电阻抗层析成像优化算法从传感器采集得到的边界电压值中提取有用信息，形成特征向量，进而进行聚类分析。由于ert成像独特的背景，可将812个像素点分成2类，其中一类代表目标，另外一类代表背景。

对于每一个像素点而言，每次激励时有且仅有一个投影区覆盖该像素点，16次循环激励后每个像素点有16个边界电压值与之相对应，形成一个16维的向量，整幅图像形成812×16的矩阵。为了有效地对数据点进行聚类，从812×16的矩阵中提取代表每个像素信息的特征。本发明为每个像素点选取两个特征。特征一：均值向量。即为每个像素点所对应的16个边界电压值的均值。特征二：方差矩阵。在均值矩阵的基础上，计算每个像素点的方差，如图1所示。图1表示以点(i，j)为中心的9个像素点，ave(i,j)表示像素点(i，j)的均值，该点的方差如公式(1)所示：

每个像素点的2维特征向量形成812×2的特征矩阵，将特征矩阵作为模糊聚类算法的输入数据集进行模糊聚类，进而得到每个像素点的隶属度值，表征每个像素点对每一类的隶属度大小。最后根据隶属度矩阵进行类别的归属判断。图2为算法流程图。

下面结合附图和实施例对本发明的进行详细的描述。

在comsol和matlab计算软件中对本发明所提供的一种基于模糊聚类的电阻抗层析成像优化算法进行仿真实验。实验方法如下：

1、在comsol中设计仿真模型，如图3所示。红色部分表示目标，蓝色部分表示背景。其中，被测场域圆心坐标为(0,0)，半径为0.5cm。目标为圆形，圆心坐标为(0,0.65)，半径为0.2cm。

2、用matlab对求解后的模型进行计算，得到208个边界电压数据。

3、分配反投电压数据。将边界电压值按照投影区进行反投，每个像素点得到16个反投的电压值。

4、计算特征向量矩阵。计算每个像素点的ave及var值，进而形成812×2的特征向量矩阵。

5、执行模糊聚类算法。对上述4中的特征向量矩阵执行模糊聚类算法，得到隶属度矩阵。

6、重建图像。根据隶属度矩阵对每个像素点进行类别的划分，进而对被测场域进行图像重建。

针对图3的仿真模型，分别用lbp及本发明的算法进行图像重建，如图4及图5所示。其中，图4为用lbp算法进行图像重建的结果，图5为用本发明算法进行图像重建的结果。

由仿真实验结果可以看出，本发明所提供基于模糊聚类的电阻抗层析成像优化算法与传统的lbp算法相比减少了伪影，具有更高的空间分辨率。

技术特征：

技术总结
本发明涉及一种基于模糊聚类的电阻抗层析成像优化算法，步骤如下：(1)得到电阻抗层析成像的图像矩阵；(2)从图像矩阵中提取代表每个像素信息的特征，特征一：均值向量，即为每个像素点所对应的各个边界电压值的均值；特征二：方差矩阵，在均值矩阵的基础上，计算每个像素点的方差，每个像素点的2维特征向量形成特征矩阵；(3)将特征矩阵作为模糊聚类算法的输入数据集进行模糊聚类，进而得到每个像素点的隶属度值，表征每个像素点对每一类的隶属度大小；根据隶属度矩阵进行类别的归属判断。

技术研发人员：李奇;岳士弘;高晓峰
受保护的技术使用者：天津大学
技术研发日：2018.07.19
技术公布日：2018.12.25