﻿Electrical Impedance Tomography Image Reconstruction Using Iterative Lavrentiev and L-Curve-Based Regularization Algorithm

Journal of Electromagnetic Analysis and Applications
Vol. 2  No. 1 (2010) , Article ID: 1232 , 6 pages DOI:10.4236/jemaa.2010.21007

Electrical Impedance Tomography Image Reconstruction Using Iterative Lavrentiev and L-Curve-Based Regularization Algorithm

Wenqin WANG1,2, Jingye CAI1, Lian YANG1 1School of Communication and Information Engineering, University of Electronic Science and Technology of China, Chengdu, China; 2Henan Key Laboratory of Coal Mine Methane and Fire Prevention, Henan University, Henan, China.

Email: wqwang@uestc.edu.cn

Received June 9th, 2009; revised August 17th, 2009; accepted August 25th, 2009.

Keywords: Electrical Impedance Tomography (EIT), Reconstruction Algorithm, Iterative Lavrentiev, Regularization Parameter, Inverse Problem.

ABSTRACT

Electrical impedance tomography (EIT) is a technique for determining the electrical conductivity and permittivity distribution inside a medium from measurements made on its surface. The impedance distribution reconstruction in EIT is a nonlinear inverse problem that requires the use of a regularization method. The generalized Tikhonov regularization methods are often used in solving inverse problems. However, for EIT image reconstruction, the generalized Tikhonov regularization methods may lose the boundary information due to its smoothing operation. In this paper, we propose an iterative Lavrentiev regularization and L-curve-based algorithm to reconstruct EIT images. The regularization parameter should be carefully chosen, but it is often heuristically selected in the conventional regularization-based reconstruction algorithms. So, an L-curve-based optimization algorithm is used for selecting the Lavrentiev regularization parameter. Numerical analysis and simulation results are performed to illustrate EIT image reconstruction. It is shown that choosing the appropriate regularization parameter plays an important role in reconstructing EIT images.

1. Introduction

Electrical impedance tomography (EIT) is an imaging technique which determines the electrical conductivity and permittivity distribution within a medium using electrical measurement from a series of electrodes on its surface [1–3]. Electrodes are brought into contact with the surface of the object being imaged. A set of voltage (or current) are applied and the corresponding currents (or voltage) are measured . These voltages and currents are then used to estimate the electrical properties of the object using an image reconstruction algorithm . The relatively poor spatial resolution of the reconstructed images in EIT is often quoted as its major disadvantages, compared with already established scanners with good resolution. In this respect, it must be clarified that the motivation of EIT is somewhat different from that of the conventional imaging techniques. Despite its limited resolution, its task is to provide a reliable, real-time, portable and cost efficient imaging tool. Depending on the particular application and the resolution specifications, EIT can sometimes provide an optimum cost-effective imaging solution. Example application areas include geophysical inversion , industrial process monitoring , and medical diagnosis .

However, the process of property estimation in EIT is a highly nonlinear, ill-conditioned, and ill-posed problem. The sensitivity matrix, which relates interior admittivity perturbations to perturbations in the boundary data, is heavily ill-conditioned with respect to inversion. So, it requires special treatment in the form of regularization or a truncation of a singular value expansion . When approaching an ill-posed problem, instead of attempting to solve the original problem one often opts to solve a similar one which is less demanding. Therefore, effective EIT image reconstruction algorithms are required.

Some papers on image reconstruction algorithms have been published [10–14], but little work on EIT image reconstruction is published. The first proposed EIT reconstruction algorithm is the equipotential back-projection . This technique reconstructs images by projecting the change in measurements at each electrode pair across the equipotential region for that current injection pattern, multiplied by an image filter. However, unlike the X-ray used in CAT scanning or general inverse scattering problems . Currents in EIT do not move in a line, but cover the region from the current source to drain. A bias is thus introduced into the image since the equipotential region is an approximation to the region producing the measurement change at each electrode. A finite element-based reconstruction algorithm was proposed in , but limited scalar/vectors can be reconstructed with this algorithm. A reconstruction algorithm for breast tumor imaging based on linearization approach was proposed in , but the number of independent conductivity regions that can be calculated is too small. In  it was assumed that the resistivity distribution could be well approximated as a linear combination of some preselected basis functions. Prior information on the structures and conductivities were used for the construction of these basis functions. The disadvantage of this method is that one may obtain misleading results when prior information is incompatible.

In this paper, we propose an iterative Lavrentiev regularization-based algorithm to reconstruct EIT images using knowledge of the noise variance of the measurements and the covariance of the conductivity distribution. As the regularization parameter should be carefully chosen, an L-curve-based optimization algorithm is applied. Numerical analysis and simulation results are provided. The remaining sections are organized as follows. Section 2 forms the problem of EIT image reconstruction and outlines the motivations of this paper. Section 3 details the iterative Lavrentiev regularization-based EIT reconstruction algorithm, followed by the simulation examples in Section 4. Finally, Section 5 concludes the whole paper.

2. Electrical Impedance Tomography (EIT)

Taking EIT imaging in medical applications as an example [18,19], different tissues of the body are shown to have different electrical characteristics. Most tissues can be considered isotropic with the exception of muscles and brain tissue which are anisotropic. It is usually assumed to be homogenous and isotropic, where the constitute parameters such as the conductivity and permittivity are independent of position and direction. Therefore, the underlying relationships that govern the interaction between EIT electricity and magnetism are the Maxwell’s equations. The medium ( ) is modeled as a closed and bounded subset of three-dimensional space with smooth boundary ( ) and uniform conductivity ( ). The electric field ( ) enclosed in is expressed in terms of the scalar potential  , (1)

where (2)

is the vector differential operator, and is the potential in the medium. The current ( ) is given by the multiplication of the conductivity. The electric field can then be computed as . (3)

As there are no interior current sources, EIT simulation in the medium can be described in terms of a scalar voltage potential satisfying Kirchoff’s voltage law . (4)

The boundary current density ( ) can then be represented by . (5)

According to this relationship, the problem of determining the potential inside the medium from boundary measurements can then be carried out .

3. Iterative Lavrentiev Regularization and L-Curve-Based EIT Image Reconstruction

To apply regularization-based image reconstruction methods, EIT problems can be formulated as a system of linear equations (6)

where is the nodal potential distribution, is the potential values on the boundary electrodes, , and are the local matrices. This equation can be represented by (7)

where is the error in the data with and . One standard approach used to solve linear estimation problems is the least square estimation, but this estimate is unsatisfactory because the calculated independent conductivity region is too small.

The inverse of the A cannot be directly computed because the singular values will grow without bound. This will amplify the noise components in the solution associated with the numerical null space of . That is to say, small measurement perturbations may produce large variations in such that the 3-norm residual error is unbounded. Therefore, the matrix is poor conditioned because EIT makes current injection and measurement on the medium surface including higher current densities near the surface where conductivity contrasts will result in more signal than for contrasts in the center. This problem can be resolved by regularization . That is, assume is an invertible and real-valued matrix, it can then be decomposed by singular value decomposition (SVD) (8)

Correspondingly, Equation (7) can be expressed by (9) Figure 1. Tikhonov-based regularization filter at different  Figure 2. The norm of the regularization solution errors as a function of the regularization parameter of As instability may arises due to division by small singular values, to overcome this problem we apply an iterative Lavrentiev-based regularization algorithm.

Lavrentiev regularization replaces Equation (7) by , (10)

where is the regularization parameter. Correspondingly, the iterative Lavrentiev regularization is (11)

This equation can be changed into (12)

Suppose the solution of the Equation (7) is  (13)

where is the normalized function with . Note that refers to the smoothness source conditions. Considering the function (14)

when , it arrives its maximum (15)

Hence, we have (16)

where is a constant parameter for the specific . Therefore, this algorithm is converged with a convergence rate of (17)

In fact, the Lavrentiev regularization just uses the filter function (18)

with a regularization parameter of , as shown in Figure 1. In conventional regularization-based reconstructed algorithms, the regularization parameter is often heuristically selected; however, it determines the cut-off of the filter, as shown in Figure 2. Further analysis results show that when the regularization parameter is very small, filtering of the noise will be inadequate and will be highly oscillatory. In contrast, when the regularization parameter is too large, the reconstructed image will be overly smooth. Therefore, some effective algorithms should be developed to optimize the regularization parameter. To reach this aim, we apply an L-curve - based optimization algorithm.

L-curve is a parametric plot of the squared norm of the regularized solution against the squared norm of the regularized residual for a range of values of regularization parameter. The L-curve criterion for regularization parameter selection is to pick the parameter value corresponding to the “corner” of this curve. Let denote the regularized solution and let denote the regularized residual. Define (19) (20)

We can then select the value of that maximizes the curvature function , (21)

where and are represented, respectively, by , (22) , (23)

The , , and denote the first and second derivatives of , with respect to . As shown in Figure 3, the L-curve has two characteristic parts: the more horizontal where the solution is dominated by the regularization errors, the vertical part where the solution is dominated by the right-hand errors. The solutions are overand under-smoothed, respectively. The corner of the L-curve corresponds to a good balance between minimization of the sizes, and the corresponding regularization parameter is a good one. In this way, an optimum regularization parameter can be optimization determined. Figure 3. The curvature of the L-curve Figure 4. Simulated distribution Figure 5. Iterative Lavrentiev reconstructed image with 4. Simulation Examples

In many industrial process and biomedical EIT applications, it is often the case that one knows in advance what materials (tissues) are included within the measurement domain. As mixtures of materials are known to have intermediate conductivities, this effectively restricts the admissible solutions, i.e., the pixels of the reconstructed images, to lie within a set of known values. In this sense, the bound-constrained (with bounds on the values of the admittivity distribution) is to locate the detected inhomogeneities.

As an example, one can allow that the admittivity distribution to be reconstructed is mainly homogenous with an unknown number of shaped inclusions. In EIT applications, the current with a frequency of 10-100kHz is widely used. These patterns are similar to those appearing in Figure 4, which shows the simulated admittivity distribution with the three inhomogeneity patterns. Although the geometry (shape and dimension) of these patterns is assumed to be known a priori, their number, admittivity values and exact location are to be recovered from the Figure 6. Regularization selection with L-curve-based optimization algorithm Figure 7. Iterative Lavrentiev reconstructed image with L-curve optimized regularization parameter

image reconstruction. For the reconstructions, simulated measurements contaminated with 2% Gaussian noise have been assumed. Regularization parameter is often used in conventional regularizationbased reconstruction methods, e.g., , but from Figure 5 we notice that it fails to capture the boundary shape and interior gap of the inhomogeneities. To get around this disadvantage, the regularization parameter is optimized with the L-curve algorithm, as shown in Figure 6. Correspondingly, Figure 7 gives the reconstructed image using the optimized regularization parameter (here is ). Comparing the two images, the enhancement is obvious. This result comes in support of the fact that this method is robust to noise.

5. Conclusions

EIT is a technique used to create images of the electrical properties in the interior of a medium from measurements on its boundary, which is particularly important for medical and industrial applications [26,27]. Usually a set of voltage or current measurements is acquired from the boundaries of a conductive volume. In this paper, we presented an iterative Lavrentiev regularization and Lcurve-based algorithm to reconstruct EIT images using knowledge of the noise variance of the measurements and the covariance of the conductivity distribution. The regularization parameter should be carefully selected, but it is often heuristically selected in conventional regularization-based reconstruction algorithms. So, an L-curvebased optimization algorithm is applied to choose the Lavrentiev regularization parameter. The method is validated with numerical analysis and simulation results. Further research efforts are planned to focus on experimental investigations [28,29] and other computational electromagnetic-based algorithms .

REFERENCES

1. M. Cheney, D. Isaacson and J. C. Newell, “Cooling water system design,” SIAM Review, Vol. 41, pp. 85–101, 1999.
2. S. Jaruwatanadilok and A. Ishimaru, “Electromagetic coherent tomography array imaging in random scattering media,” IEEE Antennas and Wireless Propagation Letters, Vol. 7, pp. 524–527, 2008.
3. B. H. Brown, “Medical impedance tomography and process impedance tomography: a brief review,” Measurement Science and Technology, Vol. 12, pp. 991–996, 2001.
4. E. Martini, G. Carli and S. Maci, “An equivalance theorem based on the use of electric currents radiating in free space,” IEEE Antennas and Wireless Propagation Letters, Vol. 7, pp. 421–424, 2008.
5. U. A. Ahsan, S. Toda, T. Takemae, Y. Kosugi and M. Hongo, “A new method for electric impedance imaging using an eddy current with a tetrapolar circuit,” IEEE Transaction on Biomedical Engineering, Vol. 56, No. 2, pp. 400–406.
6. P. Metherall, D. C. Barber, R. H. Smallwood and B. H. Brown, “Three dimensional electrical impedance tomography,” Nature, Vol. 380, pp. 509–512, 1996.
7. N. K. Soni, K. D. Paulsen, H. Dehghani and A. Hartov, “Finite element implementation of Maxwell’s equations for image reconstruction in electrical impedance tomography,” IEEE Transaction on Medical Imaging, Vol. 25, No. 1, pp. 55–61, 2006.
8. M. H. Choi, T. J. Kao, D. Isaacson and G. J. Saulnier, “A reconstruction algorithm for breast cancer imaging with electrical impedance tomography in mammography geometry,” IEEE Transaction on Biomedical Engineering, Vol. 54, No. 4, pp. 700–710, 2007.
9. W. R. B. Lionheart, “EIT reconstruction algorithms pitfalls, challenges and recent developments,” Physiological Measurement, Vol. 25, pp. 125–142, 2004.
10. M. R. Hajihashemi and M. EI-Shenawee, “Shape reconstruction using the level set method for microwave applications,” IEEE Antennas and Wireless Propagation Letters, Vol. 7, pp. 92–96, 2008.
11. G. Wang, L. J. Schultz and J. Qi, “Bayesian image reconstruction for improving detection performance of Muon tomography,” IEEE Transaction on Image Processing, Vol. 18, No. 5, pp. 1080–1089, 2009.
12. S. K. Padhi, A. Fhager, M. Persson and J. Howard, “Measured antenna response of a proposed microwave tomography system using an efficient 3-D FTD model,” IEEE Antennas and Wireless Propagation Letters, Vol. 7, pp. 689–692, 2008.
13. J. Zhou, L. Senhadji, J. L. Coatrieux and L. Luo, “Iterative PET image reconstruction using translation invariant wavelet transform,” IEEE Transaction on Nuclear Science, Vol. 56, No. 1, part 1, pp. 116–128, 2009.
14. J. Trzasko and A. Manduca, “Highly undersampling magnetic resonance image reconstruction via homotopic -minization,” IEEE Transaction on Medical Imaging, Vol. 28, No. 1, part 1, pp. 106–121, 2009.
15. Adler and R. Guardo, “Electrical impedance tomography regularized imaging and contrast detection,” IEEE Transaction on Medical Imaging, Vol. 15, No. 2, pp. 170–179, 1999.
16. M. Benedetti, M. Donelli, D. Lesselier and A. Massa, “A two-step inverse scattering procedure for the qualitative imaging of homogenous cracks in known host media-preliminary results,” IEEE Antennas and Wireless Propagation Letters, Vol. 6, pp. 592–595, 2007.
17. M. Vauhkonen, J. P. Kaipio, E. Somersalo and P. A. Karjalainen, “Electrical impedance tomography with basis constraints,” Inverse Problems, Vol. 13, pp. 523–530, 1997.
18. P. M. Edic, G. J. Saulnier, J. C. Newell and D. Isaacson, “A real-time electrical impedance tomography,” IEEE Transaction on Biomedical Engineering, Vol. 42, No. 9, pp. 849–859, 1995.
19. R. D. Cook, G. J. Saulnier, D. G. Gisser, J. C. Goble, J. C. Newell, and D. Isaacson, “ACT3: a high-speed, highprecision electrical impedance tomography,” IEEE Transaction on Biomedical Engineering, Vol. 41, No. 9, pp. 713–722, 1994.
20. L. Di Rienzo, S. Yuferev, N. Ida, “Computation of the impedance mxtrix of multiconductor transmisson lines using high-order surface impedance boundary conditions,” IEEE Transaction on Electromagnetic Compatibility, Vol. 50, No. 4, pp. 974–984, 2008.
21. N. Polydorides, Image reconstruction algorithms for soft-field tomography, Ph. D thesis, University of Manchester Institute of Science and Technology, 2002.
22. Y. Z. O’Connor and J. A. Fessler, “Fast predictions of variance images for fan-beam transmission tomography with quadratic regularization,” IEEE Transaction on Medical Imaging, Vol. 26, No. 3, pp. 335–346, 2007.
23. S. Morigi, L. Reichel and F. Sgallari, “An iterative Lavrentiev regularization method,” BIT, Vol. 43, No. 2, pp. 1–14, 2002.
24. P. C. Hansen and D. P. O’Leary, “The use of l-curve in the regularization of discrete ill-posed problems,” SIAM Journal on Scientific Computing, Vol. 14, No. 6, pp.1487–1503, 1993.
25. P. C. Hansen, “Rank-deficient and discrete ill-posed problems: Numerical aspects of linear inversion,” SIAM, Philadelphia, 1998.
26. V. Kolehmainen, M. Lassas and P. Ola, “Electrical impedance tomography problem with inaccuracy known boundary and contact impedance,” IEEE Transaction on Medical Imaging, Vol. 27, No. 10, pp. 1404–1414, 2008.
27. F. Ferraioli, A. Formisano, and R. Martone. Ramos, “Effective exploitation of prior information in electrical impedance tomography for thermal monitoring of hyperthermia treatments,” IEEE Transaction on Magnetics, Vol. 45, No. 3, pp. 1554–1557, 2009.
28. E. Hartinger, R. Guardo, A. Alder, and H. Gagnon, “Real-time management of faulty electrodes in electrical impedance tomography,” IEEE Transaction on Biomedical Engineering, Vol. 56, No. 2, pp. 369–377, 2009.
29. P. M. Ramos, F. M. Janeiro, M. Tlemcani and A. C. Seira, “Recent developments on impedance measurements with DSP-based ellipse-fitting algorithms,” IEEE Transaction on Instrumentation and Measurement, Vol. 58, No. 5, pp. 1680–1689, 2009.
30. Ergul, I. van den Bosch and L. Gurel, “Two-step lagrange interpolation method form multilevel fast multipole algorithm,” IEEE Antennas and Wireless Propagation Letters, Vol. 8, pp. 69–71, 2009.

NOTES

This work was supported in part by the Open Fund of the Henan Key Laboratory of Coal Mine Methane and Fire Prevention, Henan University under the contract number of HKLGF200803.