电阻抗成像OpenEIT 番外篇简单算法

FakeOccupational

于 2023-01-24 23:00:00 发布

阅读量1.5k

点赞数

分类专栏：算法文章标签：算法 EIT 牛顿法

本文链接：https://blog.csdn.net/ResumeProject/article/details/128753066

版权

算法专栏收录该内容

54 篇文章 4 订阅

订阅专栏

算法意义

在这里插入图片描述
$C\rightarrow Y \rightarrow V \rightarrow e$

符号

符号	含义	–
$\Omega$	研究图像区域
$\partial\Omega$	研究图像区域的边界
$\phi(x,y)$	代求电导率
$\sigma(x,y)$	节点电压 e 电导率
$j_n(x,y)$	注入表面电流密度分布
N	同态分割网格的个数
$\Delta x$	x轴向的网格边长
$\Delta y$	y轴向的网格边长，令 $p_i=\Delta x\Delta y$
C ={ $c_{x,y}$ }	连续电导率分布在实数中离散化, $c_{x,y}$

在该方法中，离散化域被建模为导纳网络，其中基尔霍夫电流压定律的应用为每个单元产生以下线性方程：

符号	含义	–
$V_{x,y}$	（x，y）节点处的离散电势
$Y_{y^+}$	导纳项

电压 = 电压/导纳率
$V_{x,y}=\frac{1}{Y_{x^{-}}+ Y_{x^{+}}+ Y_{y^{-}}+ Y_{y^{+}} }( Y_{x^{+}}V_{x+\Delta x,y}+ Y_{x^{-}}V_{x-\Delta x,y}+ Y_{y^{+}}V_{x,y+\Delta y}+ Y_{y^{-}}V_{x,y-\Delta y} - I_{ie}) \tag{1}$

符号	含义	–
$V_{x,y}$	（x，y）节点处的离散电势
$I_{ie}$	方程的已知电源项（源）

导纳项计算如下（由C推断Y（C是连续的初始假设的，Y是离散的）,以 $Y_{y^+}$ 为例）：
$Y_{y^+}=c_{x^+}*\frac{\Delta y}{\Delta x}=\frac{2c_{x,y} c_{x+\Delta x,y} }{c_{x,y}+ c_{x+\Delta x,y}} *\frac{\Delta y}{\Delta x} \tag{2}$

为了在N个未知数中求解N个方程的平方线性系统，可以使用迭代或直接方法。

Iterative solution

令 $V^{n+1}$ 表示第n+1次

过松弛参数1 < $\alpha$ < 2, $V (x, y)$ 的计算使用公式1.

Non-iterative solution（线性系统）

$Y_{(N-1,n-1)}*V_{(N-1,1)} =I_{(N-1,1)}\tag{4}$

符号	含义	–
$V=[V_1… V_{N-1} ]^T$	未知离散电压分布
$I=[I_1… I_{N-1} ]$	已知离散电流分布
$Y=[y_{ij}]_{i=1…N-1}^{j=1…N-1}$	电导率矩阵

，其中每行都有所有零元素，但通常有五个元素除外，如等式（1）所示
通过求Y逆，可以获得线性系统的解向量V

Computed voltages

为了收集计算的电压，通过放置在身体表面上的电极施加多个P电流注入。
所有实现的重构算法都需要为每个注入模式计算定义良好的电极对之间的电压差的数量G。

所计算的电压被存储到M = P*G个向量 $g^{comp}$ 中，并且可以通过求解系统（4）的P倍并考虑V向量的对应元素之间的差异来评估它们：

$g^{comp}(c)= TV'=TY'^{-1}I' \tag 5$

其中列向量c是通过堆叠先前引入的c矩阵的行而获得的，其通用元素ci表示第i个像素的电导率；T是一个变换矩阵，它只从与电极相关的单元中提取电压，并且其行具有除两个值为1和-1之外的所有零元素；
在这里插入图片描述
是两个向量，其中Vj是 $I=I_j$ 时的系统（4）的解；Y’是块对角矩阵
其中每个块仅由Y给出。因此， $g^{comp}(c)$ 的元素 $g_{ij}$ 是第i个和第j个电极之间存在的电压差。