Magnetic induction tomography (MIT) of biological tissue is used to reconstruct the changes in the complex conductivity distribution inside an object under investigation. The measurement principle is based on determining the perturbation DeltaB of a primary alternating magnetic field B0, which is coupled from an array of excitation coils to the object under investigation. The corresponding voltages DeltaV and V0 induced in a receiver coil carry the information about the passive electrical properties (i.e. conductivity, permittivity and permeability). The reconstruction of the conductivity distribution requires the solution of a 3D inverse eddy current problem. As in EIT the inverse problem is ill-posed and on this account some regularization scheme has to be applied. We developed an inverse solver based on the Gauss-Newton-one-step method for differential imaging, and we implemented and tested four different regularization schemes: the first and second approaches employ a classical smoothness criterion using the unit matrix and a differential matrix of first order as the regularization matrix. The third method is based on variance uniformization, and the fourth method is based on the truncated singular value decomposition. Reconstructions were carried out with synthetic measurement data generated with a spherical perturbation at different locations within a conducting cylinder. Data were generated on a different mesh and 1% random noise was added. The model contained 16 excitation coils and 32 receiver coils which could be combined pairwise to give 16 planar gradiometers. With 32 receiver coils all regularization methods yield fairly good 3D-images of the modelled changes of the conductivity distribution, and prove the feasibility of difference imaging with MIT. The reconstructed perturbations appear at the right location, and their size is in the expected range. With 16 planar gradiometers an additional spurious feature appears mirrored with respect to the median plane with negative sign. This demonstrates that a symmetrical arrangement with one ring of planar gradiometers cannot distinguish between a positive conductivity change at the true location and a negative conductivity change at the mirrored location.

The detection and continuous monitoring of brain oedema is of particular interest in clinical applications because existing methods (invasive measurement of the intracranial pressure) may cause considerable distress for the patients. A new non-invasive method for continuous monitoring of an oedema promises the use of multi-frequency magnetic induction tomography (MIT). MIT is an imaging method for reconstructing the changes of the conductivity deltakappa in a target object. The sensitivity of a single MIT-channel to a spherical oedematous region was analysed with a realistic model of the human brain. The model considers the cerebrospinal fluid around the brain, the grey matter, the white matter, the ventricle system and an oedema (spherical perturbation). Sensitivity maps were generated for different sizes and positions of the oedema when using a coaxial coil system. The maps show minimum sensitivity along the coil axis, and increasing values when moving the perturbation towards the brain surface. Parallel to the coil axis, however, the sensitivity does not vary significantly. When assuming a standard deviation of 10(-7) for the relative voltage change due to the system's noise, a centrally placed oedema with a conductivity contrast of 2 with respect to the background and a radius of 20 mm can be detected at 100 kHz. At higher frequencies the sensitivity increases considerably, thus suggesting the capability of multi-frequency MIT to detect cerebral oedema.

Magnetic induction tomography (MIT) is a low-resolution imaging modality which aims at the three-dimensional (3D) reconstruction of the electrical conductivity in objects from alternating magnetic fields. In MIT systems the magnetic field perturbations to be detected are very small when compared to the excitation field (ppm range). The voltage which is induced by the excitation field in the receiver coils must be suppressed for providing sufficient dynamic range. In the past, two very efficient strategies were proposed: adjusted planar gradiometers (PGRAD) and the orientation of a receiver coil with respect to the excitation coil such that the net magnetic flow is zero (zero flow coil, ZFC). In contrast to the PGRAD no voltage is induced in the ZFC by the main field. This is advantageous because two comparatively high voltages in the two gradiometer coils can never be subtracted perfectly, thus leaving a residual voltage which is prone to drift. However, a disadvantage of the ZFC is the higher susceptibility to interferences from far RF sources. In contrast, in the gradiometer such interferences are cancelled to a high degree. We developed a new type of gradiometer (zero flow gradiometer, ZFGRAD) which combines the advantages of ZFC and PGRAD. All three systems were compared with respect to sensitivity and perturbation to signal ratio (PSR) defined as the ratio of the signal change due to a magnetic perturbation field at the carrier frequency and the signal change due to shifting a metallic sphere between two test points. The spatial sensitivity of the three systems was found to be very similar. The PSR of the ZFGRAD was more than 12 times lower than that of the ZFC. Finally, the feasibility of image reconstruction with two arrays of eight excitation coils and eight ZFGRAD, respectively, was shown with a single-step Gauss-Newton reconstructor and simulated measurement data generated for a cylindrical tank with a spherical perturbation. The resulting images show a clear, bright feature at the correct position of the sphere and are comparable to those with PGRAD arrays.

Magnetic induction tomography (MIT) is used for reconstructing the changes of the conductivity in a target object using alternating magnetic fields. Applications include, for example, the non-invasive monitoring of oedema in the human brain. A powerful software package has been developed which makes it possible to generate a finite element (FE) model of complex structures and to calculate the eddy currents in the object under investigation. To validate our software a model of a previously published experimental arrangement was generated. The model consists of a coaxial coil system and a conducting sphere which is moved perpendicular to the coil axis (a) in an empty space and (b) in a saline-filled cylindrical tank. The agreement of the measured and simulated data is very good when taking into consideration the systematic measurement errors in case (b). Thus the applicability of the simulation algorithm for two-compartment systems has been demonstrated even in the case of low conductivities and weak contrast. This can be considered an important step towards the solution of the inverse problem of MIT.

We developed a 14-channel multifrequency magnetic induction tomography system (MF-MIT) for biomedical applications. The excitation field is produced by a single coil and 14 planar gradiometers are used for signal detection. The object under measurement was rotated (16 steps per turn) to obtain a full data set for image reconstruction. We make measurements at frequencies from 50 kHz to 1 MHz using a single frequency excitation signal or a multifrequency signal containing several frequencies in this range. We used two acquisition boards giving a total of eight synchronous channels at a sample rate of 5 MS s(-1) per channel. The real and imaginary parts of DeltaB/B(0) were calculated using coherent demodulation at all injected frequencies. Calibration, averaging and drift cancellation techniques were used before image reconstruction. A plastic tank filled with saline (D = 19 cm) and with conductive and/or paramagnetic perturbations was measured for calibration and test purposes. We used a FEM model and an eddy current solver to evaluate the experimental results and to reconstruct the images. Measured equivalent input noise voltage for each channel was 2 nV Hz(-1/2). Using coherent demodulation, with an integration time of 20 ms, the measured STD for the magnitude was 7 nV(rms) (close to the theoretical value only taking into account the amplifier's thermal noise). For long acquisition times the drift in the signal produced a bigger effect than the input noise (typical STD was 10 nV with a maximum of 35 nV at one channel) but this effect was reduced using a drift cancellation technique based on averaging. We were able to image a 2 S m(-1) agar sphere (D = 4 cm) inside the tank filled with saline of 1 S m(-1).