PIER | |

Progress In Electromagnetics Research | ISSN: 1070-4698, E-ISSN: 1559-8985 |

Home > Vol. 155 > pp. 1-18
## ABSOLUTE IMAGING OF LOW CONDUCTIVITY MATERIAL DISTRIBUTIONS USING NONLINEAR RECONSTRUCTION METHODS IN MAGNETIC INDUCTION TOMOGRAPHYBy B. Dekdouk, C. Ktistis, D. W. Armitage, and A. J. Peyton
Abstract:
Magnetic Induction Tomography (MIT) is a newly developing technique of electrical tomography that in principle is able to map the electrical conductivity distribution in the volume of objects. The image reconstruction problem in MIT is similar to electrical impedance tomography (EIT) in the sense that both seek to recover the conductivity map, but differ remarkably in the fact that data being inverted in MIT is derived from induction theory and related sources of noise are different. Progress in MIT image reconstruction is still limited, and so far mainly linear algorithms have been implemented. In difference imaging, step inversion was demonstrated for recovering perturbations within conductive media, but at the cost of producing qualitative images, whilst in absolute imaging, linear iterative algorithms have mostly been employed but mainly offering encouraging results for imaging isolated high conductive targets. In this paper, we investigate the possibility of absolute imaging in 3D MIT within a target for low conductivity application (σ < 5 Sm ^{-1}). For this class of problems, the MIT image reconstruction exhibits non-linearity and ill-posedness that cannot be treated with linear algorithms. We propose to implement for the first time in MIT two effective inversion methods known in non-linear optimization as Levenberg Marquardt (LMM) and Powell's Dog Leg (PDLM) methods. These methods employ damping and trust region techniques for controlling convergence and improving minimization of the objective function. An adaptive version of Gauss Newton is also presented (AGNM), which implements a damping mechanism to the regularization parameter. Here, the level of penalty is varied during the iterative process. As a comparison between the methods, different criteria are examined from image reconstructions using the LMM, PDLM and AGNM. For test examples, volumetric image reconstruction of a perturbation within homogeneous cylindrical background is considered. For inversion, an independent finite element FEM software package Maxwell by Ansys is employed to generate simulated data using a model of a 16 channel MIT system. Numerical results are employed to show different performance characteristics between the methods based on convergence, stability and sensitivity to the choice of the regularization parameter. To demonstrate the effect of scalability of absolute imaging in MIT for more realistic problems, a human head model with an internal anomaly is used to produce reconstructions for different finer resolutions. AGNM is adopted here and employs the Krylov subspace method to replace the computationally demanding direct inversion of the regularized Hessian.
2. Peyton, A. J., Z. Z. Yu, G. Lyon, S. AlZeibak, J. Ferreira, J. Velez, F. Linhares, A. R. Borges, H. L. Xiong, N. H. Saunders, and M. S. Beck, "An overview of electromagnetic inductance tomography: Description of three different systems," 3. Ke, L., X. Lin, and Q. Du, "An improved back-projection algorithm for magnetic induction tomography image reconstruction," 4. Korjenevsky, A., V. Cherepenin, and S. Sapetsky, "Magnetic induction tomography: Experimental realization," 5. Merwa, R., K. Hollaus, P. Brunner, and H. Scharfetter, "Solution of the inverse problem of magnetic induction tomography (MIT)," 6. Scharfetter, H., K. Hollaus, J. Rosell-Ferrer, and R. Merwa, "Single-step 3-D image reconstruction in magnetic induction tomography: Theoretical limits of spatial resolution and contrast to noise ratio," 7. Scharfetter, H., P. Brunner, and R. Merwa, "Magnetic induction tomography: Single-step solution of the 3-D inverse problem for differential image reconstruction," 8. Vauhkonen, M., M. Hamsch, and C. H. Igney, "A measurement system and image reconstruction in magnetic induction tomography," 9. Ma, X., A. J. Peyton, S. R. Higson, A. Lyons, and S. J. Dickinson, "Hardware and software design for an electromagnetic induction tomography (EMT) system for high contrast metal process applications," 10. Wei, H.-Y. and M. Soleimani, "Four dimensional reconstruction using magnetic induction tomography: Experimental study," 11. Wei, H.-Y. and M. Soleimani, "Volumetric magnetic induction tomography," 12. Wei, H.-Y. and M. Soleimani, "Three-dimensional magnetic induction tomography imaging using a matrix free Krylov subspace inversion algorithm," 13. Watson, S., R. J.Williams, W. Gough, and H. Griffiths, "A magnetic induction tomography system for samples with conductivities below 10 Sm 14. Soleimani, M. and W. R. B. Lionheart, "Absolute conductivity reconstruction in magnetic induction tomography using a nonlinear method," 15. Soleimani, M., W. R. B. Lionheart, A. J. Peyton, X. D. Ma, and S. R. Higson, "A three-dimensional inverse finite-element method applied to experimental eddy-current imaging data," 16. Hansen, P. C. and D. P. O'Leary, "The use of the L-curve in the regularization of discrete ill-posed systems," 17. Goharian, M., M. Soleimani, and G. R. Moran, "A trust region subproblem for 3D electrical impedance tomography inverse problem using experimental data," 18. Goharian, M., A. Jegatheesan, and G. R. Moran, "Dogleg trust-region application in electrical impedance tomography," 19. Tan, C., G. Xu, Y. Li, Y. Xuand, and G. Su, "Boundary image reconstruction based on the nonmonotonic and self-adaptive trust region method for electrical impedance tomography," 20. Tan, C., Y. Xuand, and F. Dong, "Determining the boundary of inclusions with known conductivities using a Levenberg-Marquardt algorithm by electrical resistance tomography," 21. Kleefeld, A. and M. Reißel, "The Levenberg-Marquardt method applied to a parameter estimation problem arising from electrical resistivity tomography," 22. Ren, S., F. Dong, Y. Xuand, and C. Tan, "Reconstruction of the three-dimensional inclusion shapes using electrical capacitance tomography," 23. Nielsen, H. B., "Damping parameter in Marquardt's method,", Report IMM-REP-1999-05, Department of Mathematical Modelling, 31 pages, DTU, 1999.
24. Nocedal, J. and S. J. Wright, 25. Frandsen, P. E., K. Jonasson, H. B. Nielsen, and O. Tingleff, 26. Gazzola, S. and J. Nagy, "Generalized Arnoldi-Tikhonov method for sparse reconstruction," 27. Biro, O. and K. Preis, "On the use of the magnetic vector potential in the Finite Element Analysis of 3 dimensional eddy currents," 28. Hollaus, K., C. Magele, R. Merwa, and H. Scharfetter, "Numerical simulation of the eddy current problem in magnetic induction tomography for biomedical applications by edge elements," 29. Morris, A., H. Griffiths, and W. Gough, "A numerical model for magnetic induction tomographic measurements in biological tissues," 30. Pham, M. H. and A. J. Peyton, "A model for the forward problem in magnetic induction tomography using boundary integral equations," 31. Dekdouk, B., W. Yin, C. Ktistis, D.W. Armitage, and A. J. Peyton, "A method to solve the forward problem in magnetic induction tomography based on the weakly coupled field approximation," 32. Dyck, D. N., D. A. Lowther, and E. M. Freeman, "A method of computing the sensitivity of electromagnetic quantities to changes in materials and sources," 33. Gabriel, S., R. W. Lau, and C. Gabriel, "The dielectric properties of biological tissues. 2 Measurements in the frequency range 10 Hz to 20 GHz," 34. Madsen, K., H. B. Nielsen, and O. Tingleff, 35. Horesh, L., "Some novel approaches in modelling and image reconstruction for multi-frequency electrical impedance tomography of the human brain,", PhD Thesis, 137-173, UCL, 2006.
36. Polydorides, N., W. R. B. Lionheart, and H. McCann, "Krylov subspace iterative techniques: On the detection of brain activity with electrical impedance tomography," |

© Copyright 2014 EMW Publishing. All Rights Reserved