一种电容耦合电阻层析成像图像重建方法技术

本发明专利技术公开了一种电容耦合电阻层析成像图像重建方法，具体包含下列步骤：(1)根据实验获得边界测量值，并计算出灵敏度矩阵；(2)提出本发明专利技术的目标函数；(3)为了便于计算，引入1个辅助变量w，写出目标函数的最小化增广拉格朗日函数；(4)求解目标函数的算法，包括设置初始化参数，更新正则化参数λ，更新辅助变量w，更新电导率分布判断迭代是否符合迭代终止条件；(5)根据求解出的最优灰度值完成图像的重建。本发明专利技术有效提高了电容耦合电阻层析成像重建图像的质量和图像分辨率，一定程度上减少了伪影，简化了算法的复杂度，增强了算法的实时性。时性。时性。

【技术实现步骤摘要】
一种电容耦合电阻层析成像图像重建方法


[0001]本专利技术属于电容耦合电阻层析成像
，具体涉及一种电容耦合电阻层析成像图像重建方法，即通过新型的全变分正则化方法实现电容耦合电阻层析成像图像的重建。

技术介绍

[0002]电学层析成像技术(Electrical Tomography，ET)出现于二十世纪八十年代后期，是一种基于电学特性(电导率/介电系数/磁导率/复导纳)敏感机理的过程层析成像技术。它通过边界测量值得出被测区域内介质的分布信息，运用相应的成像算法，进而对电学特性的分布信息进行成像。电学层析成像常用于工业检测领域，在生物医学成像领域也有很大的应用前景。与其他成像技术相比，电阻层析成像具有产生高分辨率图像的优点，同时具有无辐射、无侵入、相对较低的成本、快速响应和结构简单等优点。然而，在传统的电阻层析成像中，电极和导电介质之间的直接接触会引发问题。电阻层析成像图像对电极属性敏感，如接触阻抗。在医疗应用中，接触阻抗的值会随着身体运动和研究区域而变化。另外，它对接触层的特性很敏感，因此在临床实验中缺乏边界特性可能导致不准确。在工程领域，与导电液体长时间接触后，电极的电化学腐蚀效应和极化效应可能会导致严重的误差，电极的污染可能会带来严重的测量偏差。2010年，王等人提出了一种称为电容耦合电阻层析成像(CCERT)的非接触方法。基于电容耦合非接触电导率检测(C4D)技术，CCERT通过在电极和导电液体之间插入绝缘层来避免接触误差。与其它电学层析成像技术一样，CCERT也存在高度非线性和不适定的逆问题，阻碍了电学层析成像技术的发展。用于解决逆问题的传统算法包括非迭代方法和迭代方法，面临重建速度和精度的挑战。到目前为止，已经提出了许多方法用于CCERT中的图像重建。如B L Wang等人提出了一种结合Tikhonov正则化和同步迭代重建技术的混合图像重建方法，引入Tikhonov正则化获得初始重建图像，SIRT用于获得最终的重建图像，用于恢复CCERT中的电导率分布。W H Tan等人是通过将Levenberg
‑
Marquardt(L
‑
M)方法与同步代数重建技术相结合来实现CCERT的图像重建。以L
‑
M方法得到的图像灰度分布为初始迭代值，采用同步代数重建技术实现重建。Y X Wang等人提出了一种将线性反投影算法和K均值聚类算法相结合的图像重建方法。利用线性反投影算法获得原始重建图像，利用K均值聚类算法获得灰度阈值，进而实现图像重建。值得注意的是，在CCERT中对图像重建的研究很少。虽然已经获得了一些令人满意的结果，但是重建图像的空间分辨率仍然需要提高。目前解决不适定问题比较常用的方法是正则化方法，通过在目标函数中加入正则项的方法来稳定解。为了处理典型的不适定反问题，Tikhonov正则化是一种经典的方法。通过添加平滑正则化项，不适定问题被正则化，然而，Tikhonov正则化方法不能很好地保留图像的边缘，从而降低了重建图像的分辨率。与Tikhonov正则化相比，全变分正则项的使用使得图像质量有了很大的改善，全变分的优势是在抑制重建图像噪声的同时，仍保留图像的边缘细节信息。然而，全变分正则化方法在重建图像的平滑区域时经常出现不想要的阶梯伪影，称为“阶梯效应”。为了抑制阶梯效应，提高图像的分辨率，本专利技术提
出了一种新型的电容耦合电阻层析成像图像重建方法，背景清晰，目标物重建准确，抗噪性能好。

技术实现思路

[0003]本专利技术解决的技术问题是提供了一种电容耦合电阻层析成像图像重建方法，相比于全变分正则化方法，本专利技术所提出的新型的全变分正则化方法分辨率更高，背景更清晰，在去阶梯效应方面展现出了优越的性能，明显的抑制了重建图像的阶梯效应，而且在保证抑制阶梯效应效果的基础上，能够进一步增强对边缘信息的保留。
[0004]本专利技术的技术方案为：本专利技术提出了一种新型的全变分正则化算法来求解电容耦合电阻层析成像的逆问题，其最小化的目标函数为：其中，||
·
||为欧几里得范数，S代表灵敏度矩阵，它反映了在被测区域中由电导率变化引起的电阻变化，g表示电导率变化，H表示等效电阻的变化，D
p
表示p阶有限差分，D是离散梯度，是预测的最优电导率，λ是正则化参数，用于平衡保真度项和正则化项。
[0005]在目标函数中引入1个辅助变量w，目标函数还可以表示为：
[0006][0007]由于上式的不可微性和非线性，仍然难以直接有效地求解，将上式表示为最小化增广拉格朗日函数：
[0008][0009]式中，T表示矩阵的转置算子，v,β和μ是增广拉格朗日乘子。
[0010]重建的算法包含以下步骤：(1)针对被测对象，获取重建所需的边界测量值和灵敏度矩阵：边界测量值是将被测对象置于电容耦合电阻层析成像测量系统中，被测场域外均匀分布12个电极，采用电压激励、电流测量的模式，采集单电极激励、单电极测量下的边界测量值，共获得132个边界测量值。(2)提出本专利技术的目标函数。(3)为了便于计算，引入1个辅助变量w，写出目标函数的最小化增广拉格朗日函数。(4)求解目标函数的算法，包括设置初始化参数，更新正则化参数λ，更新辅助变量w，更新电导率分布判断迭代是否符合迭代终止条件。(5)根据求解出的最优灰度值完成图像的重建。
[0011]本专利技术与现有技术相比具有以下有益效果：本专利技术提出了一种新型的全变分正则化方法，用于实现电容耦合电阻层析成像图像的重建。本专利技术对这种新型的全变分正则化方法的性能进行了定性和定量分析。结果表明，所提出的新型的全变分正则化方法成像的质量更高，背景更清晰，在保证有效抑制阶梯效应的基础上，可以进一步地增强对重建图像边缘信息的保留。
附图说明
[0012]图1为本专利技术的新型的全变分正则化电容耦合电阻层析成像方法图像重建的流程框图；
[0013]图2为本专利技术的电容耦合电阻层析成像传感器电极阵列图；
[0014]图3为选取的三种经典模型在相同电导率下，分别通过Landweber方法，Newton
‑
Raphson方法，Tikhonov方法和新型的全变分正则化电容耦合电阻层析成像方法图像重建的结果图；
[0015]图4为四种方法在不同电导率下对模型重建的结果图；
[0016]图5为三种经典模型在相同电导率下，通过四种方法进行重建时计算的模糊半径对比表；
[0017]图6为三种经典模型在相同电导率下，通过四种方法进行重建时计算的时间对比表；
[0018]图7为在相同电导率下，四种方法在1％噪声水平下对三种模型重建的结果图；
[0019]图8为在不同电导率下，四种方法在1％噪声水平下对模型重建的结果图；
[0020]图9为在相同电导率下，四种方法在1％噪声水平下，对三种模型重建时计算的模糊半径对比表。
具体实施方式
[0021]结合附图详细对本专利技术的新型的全变分正则化电容耦合电阻层析成像方法加以说明。
[0022本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种电容耦合电阻层析成像图像重建方法，其特征在于具体步骤为：步骤一：针对被测对象，获取重建所需的边界测量值和灵敏度矩阵：边界测量值是将被测对象置于电容耦合电阻层析成像测量系统中，被测场域外均匀分布12个电极，采用电压激励、电流测量的模式，采集单电极激励、单电极测量下的边界测量值，共获得132个边界测量值；灵敏度矩阵是根据不含内含物的空场的边界测量值，结合灵敏度理论计算得到的，计算公式为：其中S
ij
是第j个电极对对第i个电极对的灵敏度系数，是所有元素电导率为σ0的第i个元素的电阻测量值，是第j个元素电导率从σ0变为σ1时第i个元素的电阻测量值；步骤二：为了便于计算，引入1个辅助变量w，目标函数表示为：由于上式的不可微性和非线性，仍然难以直接有效地求解，将上式表示为最小化增广拉格朗日函数：式中，T表示矩阵的转置算子，v,β和μ是增广拉格朗日乘子；步骤三：求解目标函数算法如下：1)初始化：0＜δ,μ,η＜1,C0＝L
A
(w0,g0)；2)更新辅助变量w：式中，m表示第m次迭代；3)更新电导率分布g
m+1
＝(β(D
p
)
T
D
p
+μS
T
S)
+
((D
...

【专利技术属性】
技术研发人员：施艳艳，廖娟娟，王萌，杨坷，李亚婷，刘镇琨，
申请(专利权)人：河南师范大学，
类型：发明
国别省市：