A quantitative algorithm for parameter estimation in magnetic induction tomography

Magnetic induction tomography aims at producing images of the spatial distribution of electric conductivity and/or magnetic permeability in a given region of space. The associated inverse problem is ill posed and strongly nonlinear, making it very hard to deal with. However, if 'a priori' knowledge is incorporated into the solution, the problem may be simplified and an inversion procedure that produces quantitative results for some situations may be obtained. This paper describes an approach that follows this rationale. The spatial distribution in the region of interest is restricted to a single object of cylindrical shape and the physical property to be imaged to the magnetic permeability. Solving the problem will mean in this sense estimating only four parameters: two spatial coordinates for the position of the cylinder's centre, its radius and its magnetic permeability. The geometry of data collection is simple enough for the direct problem to be solved analytically and allows one to concentrate on the main characteristics of the inversion. It is shown that the problem is invertible, although in the limit, when the radius tends to zero, this is no longer true. A Newton–Raphson- type algorithm was developed for inversion purposes. It converges in a small number of iterations and produces exact values of the parameters for a predefined range of magnetic permeability and cylinder radius, if no noise is present. The inversion procedure is studied and its performance evaluated in several situations, related to the amount of S/N in the measurement vector, the degree of overdetermination in the system of equations and the choice of the sampling points.

[1]  Tat Soon Yeo,et al.  Full-wave computed tomography. Part 2: Resolution limits , 1984 .

[2]  R. Bates,et al.  Full-wave computed tomography. Part 4: Low-frequency electric current CT , 1985 .

[3]  D. Isaacson Distinguishability of Conductivities by Electric Current Computed Tomography , 1986, IEEE Transactions on Medical Imaging.

[4]  J. Newell,et al.  Current topics in impedance imaging. , 1987, Clinical physics and physiological measurement : an official journal of the Hospital Physicists' Association, Deutsche Gesellschaft fur Medizinische Physik and the European Federation of Organisations for Medical Physics.

[5]  C. Trowbridge,et al.  The Analytical and Numerical Solution of Electric and Magnetic Fields , 1992 .

[6]  Petri Ola,et al.  Electromagnetic Inverse Problems and Generalized Sommerfeld Potentials , 1996, SIAM J. Appl. Math..

[7]  Joaquim Ferreira,et al.  An overview of electromagnetic inductance tomography: Description of three different systems , 1996 .

[8]  Anthony J. Peyton,et al.  Development of Electromagnetic Tomography (EMT) for Industrial Applications. Part 1: Sensor Design and Instrumentation , 1999 .

[9]  H. Griffiths Magnetic induction tomography , 2001 .

[10]  Anthony J. Peyton,et al.  A combined digital camera – EMT system to measure human body composition. , 2003 .

[11]  Manuchehr Soleimani,et al.  Forward problem in 3D magnetic induction tomography (MIT) , 2003 .

[12]  Hermann Scharfetter,et al.  Detection of brain oedema using magnetic induction tomography: a feasibility study of the likely sensitivity and detectability. , 2004, Physiological measurement.