An Investigation of Modeling Error of EIT Reconstruction

Electrical Impedance Tomography is an imaging method which attempts to reveal the conductivity distribution of a domain based on the electrical boundary condition. For time difference EIT, the voltage difference at two time steps is employed for reconstruction. This is an ill-posed inverse problem, especially, it is non-linear. The currently available EIT devices are all based on linearized reconstruction algorithms. The linearized reconstruction employs a reconstruction matrix which is essentially a regularized pseudo inverse of the Jacobian matrix. This reconstruction matrix multiplying the voltage differences will provide a distribution of conductivity changes. However, the linearized reconstruction contains modelling error. In this paper, we study the modelling error caused by linearization based on the shunt model through simulations. Specifying a current injection pattern in simulation, at each time step a simulated voltage measurement can be calculated from Maxwell’s equations. The voltage difference between two time steps can be obtained. On the other hand, according to the assumption of linearized reconstruction, the voltage difference is assumed to be the Jacobian matrix multiplying the conductivity distribution changes. The discrepancy between these two voltage differences will be studied.


Introduction
Electrical Impedance Tomography (EIT) is a radiation-free imaging method.It attempts to reveal the conductivity distribution changes inside the human body of two time instants through electrical data obtained via the electrodes attached to the boundary.In lung EIT, commonly 16 electrodes are placed equidistantly on the boundary of a horizontal chest plane.
We denote the conductivity of the domain changes between two time steps by σ ∆ and the measured voltage changes on the electrodes by a vector V ∆ .Under the FEM framework with M elements, the conductivity change σ ∆ is represented by a M 1 × vector.Approximately, there exists the following relation: where J denotes the Jacobian matrix calculated at the constant conductivity 1: Jacobian matrix is commonly calculated by studying the first order perturbation of conductivity on each element [1].Briefly, under FEM framework, the potential distribution can be solved by forward model.The perturbations on each element can be determined by the simulated potential information [1].
However, such linearized forward and inverse model contains modeling error.EIT reconstruction is a non-linear inverse problem.There are many modelling error affects the imaging quality of EIT.For example, the deformations of the boundary shape the inhomogeneous background conductivity as well as the shift of contact impedance on the electrodes.In this study we focus only on the modelling error caused by linearization.To this end, we simplify the simulation models by specifying the domain shape to be unit disk and the known background conductivity.

EIT Current Injection and Measurement
On the currently available EIT devices, all electrical boundary conditions are collected through the electrode-band located around a plane of the human body.The electrical current injection and voltage measurement are always performed with a period scheme.The most commonly employed of the current injection and voltage measurement pattern is called adjacent pattern.At each time step, alternative currents with fixed amplitude are sequentially injected through two neighboring electrodes (Figure 1).For each current injection, the induced voltages are measured through adjacent electrodes with a fixed order.By fixing the injection and measurement pattern, the voltage measurement can be saved as a vector V .According to Ohm's law this voltage measurement is only depending on the conductivity distribution within the domain.In simulations, the voltage measurements are calculated through forward model.We will outline this forward model in Section 2.2.

Voltage Measurement Calculated from Maxwell's Equation
In this subsection we outline the forward model calculation based on Maxwell's equation.In order to simplify the formulations, we employ the shunt model boundary value problem.This model neglects the effect of contact impedance.In the absence of interior current sources, there is Ohm's law: .This voltage distribution is calculated by using EIDORS toolbox [2].Hence the voltage measurement differences along the boundary between two time steps can also be exactly calculated.

Voltage Calculated from the Linearized Model
With a given conductivity distribution, the Jacobian matrix J , also called the sensitivity matrix, can be calculated.Intuitively, a row of J represents the response of the voltage changes with respect to small impedance changes of each FEM element.In real applications, the actual conductivity distribution is unknown.In EIT society, the Jacobian matrix is commonly calculate with respect to the homogenous conductivity distribution: We adopt this convention in this study.According to the linearized model, expressed as Equation ( 1), the voltage difference between two time points can be approximately calculated by simply multiply the Jacobian matrix J and the conductivity changes.On the other hand, the exact voltage differences between two time steps can be calculated according to Section 2.2.Under the specified adjacent simulation settings, the discrepancy between these two voltage differences is because of the modeling error induced by linearization.Explicitly, the following discrepancy is identified as the modeling error induced by linearization: Linearization error: = J σ V ⋅ ∆ − ∆ (4)

Simulation Phantom
The discrepancy between the voltage differences calculated from Maxwell's equ-

Results and Discussion
In this section, we present the first simulation results.We employ the 2D FEM model with disk shape for simulation as displayed in section 2.4.The linearization modeling errors, pronounced by the voltage discrepancy descript in the last section will be studied.In below the modeling errors induced by linearization are studied with respect real to the voltage changes.The relation between modeling errors and voltage differences will be investigated according to each of the four given phantoms.For the phantom A with contrast conductivity

Conclusion
In this study, the modeling errors induced by linearization are investigated according to simulation.The first simulation results indicated that such modeling errors may have linear correspondence with the voltage differences.An exceptional case appears when the contrast's shape has large asymmetry.For this situation, the linear correspondence is not obvious.

Figure 1 .
Figure 1.Adjacent current injection and voltage measurement pattern.In this pattern, 32 electrodes are attached around the boundary of the domain.Alternative current with 5mA are injected from the neighboring electrodes.The induced voltages are measured through the adjacent electrodes.

Figure 2 .
Figure 2. Conductivity distribution and the voltage distribution changes.For a given conductivity distribution (shown in left image), the voltage distribution changes with respect to a current injection can be calculated through Maxwell's equation.This voltage distribution changes is exact.The simulated voltage measurment changes arround the boundary can be derived from this voltage distribution changes.
ation and the linearized model will be studied through simulation.A 2D disk FEM model was constructed with the help of MATLAB toolbox EIDORS.This model has 1215 triangular elements.Around the boundary of the disk phantom, 16 electrodes are attached equidistantly.The contact impedance between each electrode and the disk domain is fixed to be 0.01 Ω•m.The simulated current with the amplitude of 10 mA was injected through these electrodes with the adjacent pattern.Different conductivity contrasts are embedded into this disk phantom.In this first study, we increase the constant conductivity of the contrast region from 1.01 Sm −1 to 2 Sm −1 .These simulated phantoms reflex the different aspects of the contrasts.Such as different location, size as well as asymmetry.
Figure 4.The corresponding modeling error, defined by Equation (4), is plotted in the same figure.For three different conductivities:

Figure 4 .
Figure 4. Voltage differences and the modeling error induced by linearization.For the phantom A with contrast conductivity 1 1 Sm σ

Figure 5 .
Figure 5. Voltage differences and the modeling error induced by linearization.This is the same plot as Figure 4 with respect to the phantom B.

Figure 6 .
Figure 6.Voltage differences and the modeling error induced by linearization.This is the same plot as Figure 4 with respect to the phantom C.

Figure 7 .
Figure 7. Voltage differences and the modeling error induced by linearization.This is the same plot as Figure 4 with respect to the phantom D.