Anisotropic TV Regularization in Electrical Impedance Tomography: An Experimental Study ()
1. Introduction
Electrical Impedance Tomography (EIT) is a method of reconstructing interior conductivity distribution of the imaging target. In EIT, a series of low frequency current is injected through the electrodes attached around the boundary and we measure the induced voltages. Compared with the widely used imaging methods, such as computed tomography (CT), magnetic resonance imaging (MRI) and ultrasound imaging, EIT has a promising clinical application due to its advantages of non-invasive, no radicalization, high time resolution, etc. EIT also has very important applications in other fields, for example, to locate high flow drilling [1], to apply in Electrical Cell-substrate Impedance Spectroscopy (ECIS) [2].
However, EIT is a typical ill-posed problem. To be precise, the voltage response due to the conductivity changes decreases rapidly with the distance between the sources of conductivity changes increases [3]. Moreover, EIT measurements often suffer from noises and artifacts especially in clinical environments. Regularization is a widely used technique to deal with the ill-posedness. The main idea of regularization is to approximate the ill-posed problem by a well-posed one.
Depending on the prior information of the target conductivity distribution, the researchers proposed Tikhonov regularization [4], sparsity based regularization [5], total variation (TV) based regularization [6] and so on. Among all these regularizers, except TV based method, the remaining will blur the edges of the internal structure. TV regularizer is widely used in EIT, because it has the ability of preserving the discontinuities structure of the imaging target.
There exists difficulty in using the TV regularization due to its non-differentiability structure [6]. Many methods have been proposed to overcome this difficulty, such as Newton’s method [7] and the Primal Dual-Interior Point Methods (PDIPM) [8]. Both methods, however, are either unstable or too time consuming [9]. Split Bregman method is also used to deal with the non-differentiability, however, the ability of preserving the edges is decreased [10].
There are two common variants of TV, namely the isotropic TV and the anisotropic TV [11]. However, the anisotropic TV may distort the internal inhomogeneities along coordinate axis. Two possible reasons may cause this distortion, one is the characteristics of the model itself, and another is the algorithm to solve the inverse problem. Gonzalez [11] overcame the non-differentiability by introducing an auxiliary smoothing parameter
and using the Guass-Newton method to solve the TV regularization problem, and verified the distortion characteristics by numerical simulation and tank experiments. However, in EIT, there is still a lack of real tank experiments to solve anisotropic TV regularization problems with ADMM [12] and verify its possible distortion.
In this article, we use an iteration scheme based on the alternating direction method of multipliers (ADMM) to overcome the non-differentiability of anisotropic TV regularization method and do the tank experiment to verify its property of distortion.
2. Forward and Inverse Problems in EIT
To simplify EIT imaging process, we assume that we use 16-channel EIT system and adopt adjacent current injection and measurement pattern. To be specific, low-frequency current with magnitude I is injected through j-th neighboring electrodes pair
for
around the boundary of imaging target
, where we denote
. Then the potential distribution of j-th injection
is governed by the following equations [13]:
Here,
is the outward unit vector to
, ds is the surface element. The voltage between
subject to the j-th injection can be measured:
Using EIT scanner we can measure the following datum:
where we neglect the measurement near the driving electrode to minimize the measurement error.
The relation between
and
can be expressed approximately by the reciprocity principle
(1)
where
is a position inside
. The EIT problem is to reconstruct the conductivity distribution
using the measured voltage datum
and the relation (1). However, the above equation is nonlinear. We consider to linearize the above equation. To be specific, assuming
, we can approximate
by replacing
by
:
(2)
where
the potential computed under the reference conductivity
. For the purpose of computerized reconstruction, we discretize the imaging target into finite element elements (
). Assume that on each element
the conductivity is a constant. Let
denote the value of
on the kth element. Then
can be approximated by
. Thus, (2) can be written as
where
is the sensitivity matrix (or jacobian matrix) given by
with
,
,
.
In practice, the number of measured data
is fewer than N (the total number of elements for
). We can find an estimate of
by minimizing the following least square problem:
Since the above problem is ill-posed, the regularization technique is widely used to deal with this difficulty. The most common used regularization method is Tikhonov regularization (TR), which is to solve
where
is the regularization matrix,
is the regularization parameter. However, it has an excessively smooth effect on the solution, which will blur the edge of reconstructed image.
One technique to preserve the discontinuous boundary is the Total Variation (TV) regularization. It was first introduced in [6] for image denoising and applied in EIT inverse problem in [8]. The two basic variants of TV are isotropic TV
(3)
and anisotropic TV
(4)
where
are the first order discrete partial derivative operators in the horizontal direction and the vertical direction respectively.
To solve the non-differentiability problem (3) and (4), an auxillary smoothing parameter
was introduced in the time marching method for the corresponding Euler-Lagrange equation [6] and the primal-dual algorithms [8]. However, both of the methods are too time consuming to be used in the online mode. We intend to solve non-differentiability problem (4) using ADMM algorithm. It bypasses the difficulty by using splitting scheme and soft thresholding. To be precise, we first denote
. By introducing the a new variable
, (4) can be written into a constrained minimization problem:
(5)
The augmented Lagrangian functional for (5) is defined as:
where the variable
is the Lagrangian multiplier,
is a penalty parameter.
By updating
and
in an alternate order, we can get its ADMM iteration scheme:
The detail of solving the above equations can refer to [9].
3. Experimental Studies
In this section, we compare the performance of isotropic TV using PDIPM and anisotropic TV using ADMM by tank experiments. The results in this article are carried out by EIDORS [14].
In order to quantitatively evaluate the regularization method, two indicators are introduced, respectively
and
. Since we don’t know the true conductivity in the tank experiments, we use the conductivity reconstructed by TV regularization as a reference value. We also compare the time spent in the reconstruction process.
We did two groups of tank experiments using the Sciospec 16-channel EIT system [15] shown in Figure 1(a). The cylindrical tank with a radius of 10 cm is filled with the tap water. Sixteen electrodes are attached to the perimeter of the tank. We use the adjacent current injection and measurement pattern. The current with the maximum amplitude 1 mA and the frequency 10 kHz is injected using EIT device. The speed of the data acquisition is about 20 frames/s. The sensitivity matrix
is computed under the calculated conductivity of the water using the measured voltage data (
). We set the imaging plane to be the electrode plane. The finite element model of the imaging plane is shown in Figure 1(b).
In the first group of experiments, carrot, cucumber and bean curd are placed inside the tank, separately. In the second group of tests, carrot, cucumber and bean curd are placed inside the tank in pairs. The parameters for the anisotropic TV method are set to
,
. The parameters of TR and isotropic TV are set to be optimal empirically. The reconstructed images of single inclusion using Tikhonov regularization (TR), isotropic TV and anisotropic TV are shown in Figure 2. Figure 3 shows the reconstructed images of two inclusions cases. Table 1 and Table 2 illustrate the behavior of RE, PSNR of the single and two inclusions tank experiments, respectively. Table 3 compares the computational time of the tank experiments.
There are several observations from the results. All the three regularization methods can well capture the main feature of the inner object, including position and shape. Numerical indicators show that using ADMM to solve the anisotropic
(a) (b)
Figure 1. (a) Tank measurement setup. (b) Finite element model of imaging plane.
Figure 2. Reconstructed images of single inclusion using different regularization methods. The three rows show respectively the reconstructed images of carrot, bean curd and cucumber. The first column is the images of the tank. The last three columns are the reconstructed images using Tikhonov regularization (TR), isotropic TV with PDIPM and anisotropic TV with ADMM.
Figure 3. Reconstructed images of two different inclusions using different regularization methods. From top to bottom: 1) carrot and bean curd, 2) carrot and cucumber, 3) bean curd and cucumber. The first column is the images of the tank. The last three columns are the reconstructed images using Tikhonov regularization (TR), isotropic TV with PDIPM and anisotropic TV with ADMM.
Table 1. The behaviors of RE and PSNR for single inclusion tank model.
Table 2. The behaviors of RE and PSNR for two inclusions tank model.
Table 3. Comparison of computational time for tank experiments.
TV regularization problem can obtain an accurate image and take less time than using PDIPM. However, the images reconstructed by isotropic TV have obvious ladder effect. As a result, there exists pseudo edge. As expected, the edges of the reconstructed images using the anisotropic TV distort along the coordinate axes, whether it has one or two inclusions.
4. Conclusion and Future Work
In this article, we use experimental data to verify the difference between isotropic TV and anisotropic TV. Experimental results demonstrate that the reconstructed image with anisotropic TV regularization will cause geometric distortions along the coordinate axis. This proves that it is the model itself that causes the distortion. In the future work, we will focus on a method that can avoid the distortions along the coordinate axis and do not depend on the selection of regularization parameters.
Acknowledgements
The author would like to thank editor and referees for their valuable advice for the improvement of this article.