Optical breast shape capture and finite-element mesh generation for electrical impedance tomography

Advances in electrical impedance tomography-based brain imaging

Novel advances in the field of brain imaging have enabled the unprecedented clinical application of various imaging modalities to facilitate disease diagnosis and treatment. Electrical impedance tomography (EIT) is a functional imaging technique that measures the transfer impedances between electrodes on the body surface to estimate the spatial distribution of electrical properties of tissues. EIT offers many advantages over other neuroimaging technologies, which has led to its potential clinical use. This qualitative review provides an overview of the basic principles, algorithms, and system composition of EIT. Recent advances in the field of EIT are discussed in the context of epilepsy, stroke, brain injuries and edema, and other brain diseases. Further, we summarize factors limiting the development of brain EIT and highlight prospects for the field. In epilepsy imaging, there have been advances in EIT imaging depth, from cortical to subcortical regions. In stroke research, a bedside EIT stroke monitoring system has been developed for clinical practice, and data support the role of EIT in multi-modal imaging for diagnosing stroke. Additionally, EIT has been applied to monitor the changes in brain water content associated with cerebral edema, enabling the early identification of brain edema and the evaluation of mannitol dehydration. However, anatomically realistic geometry, inhomogeneity, cranium completeness, anisotropy and skull type, etc., must be considered to improve the accuracy of EIT modeling. Thus, the further establishment of EIT as a mature and routine diagnostic technique will necessitate the accumulation of more supporting evidence.

Absolute Reconstructions Using Rotational Electrical Impedance Tomography for Breast Cancer Imaging

A rotational Electrical Impedance Tomography (rEIT) methodology is described and shown to produce spatially accurate absolute reconstructions with improved image contrast and an improved ability to distinguish closely spaced inclusions compared to traditional EIT on data recorded from cylindrical and breast-shaped tanks. Rotations of the tank without altering the interior conductivity distribution are used to produce the rEIT data. Quantitatively, rEIT was able to distinguish two inclusions that were 1.5 cm closer together than traditional EIT could achieve for inclusions placed 2 to 3 cm from the center for the cylindrical tank, and rEIT was able to distinguish two tumor-like inclusions where traditional EIT could not reliably do so. Mathematical analysis showed that rEIT improves the number of stable singular vectors by up to 4.2 and 4.7 times than that of traditional EIT for the cylindrical and breast-shaped tanks, respectively, which is an indication of improved resolution. Direct investigations into measurements revealed minimum rotation angles that should yield data uncorrupted by noise. Two inverse approaches (one that inverts then fuses the data (I/DF) and one that fuses the data then inverts (DF/I)) and two mesh modeling approaches were considered. It was found that DF/I produces far better results compared to I/DF and a rotated-mesh approach produces further improvements. The ability to obtain improved absolute reconstructions using rEIT on a practical clinical scenario (breast-shaped tank experiment) is an important step towards using rEIT to improve previous EIT results in medical applications.

Mutual information as a measure of image quality for 3D dynamic lung imaging with EIT

We report on a pilot study of dynamic lung electrical impedance tomography (EIT) at the University of Manchester. Low-noise EIT data at 100 frames per second were obtained from healthy male subjects during controlled breathing, followed by magnetic resonance imaging (MRI) subsequently used for spatial validation of the EIT reconstruction. The torso surface in the MR image and electrode positions obtained using MRI fiducial markers informed the construction of a 3D finite element model extruded along the caudal-distal axis of the subject. Small changes in the boundary that occur during respiration were accounted for by incorporating the sensitivity with respect to boundary shape into a robust temporal difference reconstruction algorithm. EIT and MRI images were co-registered using the open source medical imaging software, 3D Slicer. A quantitative comparison of quality of different EIT reconstructions was achieved through calculation of the mutual information with a lung-segmented MR image. EIT reconstructions using a linear shape correction algorithm reduced boundary image artefacts, yielding better contrast of the lungs, and had 10% greater mutual information compared with a standard linear EIT reconstruction.

Multi-GPU Jacobian accelerated computing for soft-field tomography

Image reconstruction in soft-field tomography is based on an inverse problem formulation, where a forward model is fitted to the data. In medical applications, where the anatomy presents complex shapes, it is common to use finite element models (FEMs) to represent the volume of interest and solve a partial differential equation that models the physics of the system. Over the last decade, there has been a shifting interest from 2D modeling to 3D modeling, as the underlying physics of most problems are 3D. Although the increased computational power of modern computers allows working with much larger FEM models, the computational time required to reconstruct 3D images on a fine 3D FEM model can be significant, on the order of hours. For example, in electrical impedance tomography (EIT) applications using a dense 3D FEM mesh with half a million elements, a single reconstruction iteration takes approximately 15-20 min with optimized routines running on a modern multi-core PC. It is desirable to accelerate image reconstruction to enable researchers to more easily and rapidly explore data and reconstruction parameters. Furthermore, providing high-speed reconstructions is essential for some promising clinical application of EIT. For 3D problems, 70% of the computing time is spent building the Jacobian matrix, and 25% of the time in forward solving. In this work, we focus on accelerating the Jacobian computation by using single and multiple GPUs. First, we discuss an optimized implementation on a modern multi-core PC architecture and show how computing time is bounded by the CPU-to-memory bandwidth; this factor limits the rate at which data can be fetched by the CPU. Gains associated with the use of multiple CPU cores are minimal, since data operands cannot be fetched fast enough to saturate the processing power of even a single CPU core. GPUs have much faster memory bandwidths compared to CPUs and better parallelism. We are able to obtain acceleration factors of 20 times on a single NVIDIA S1070 GPU, and of 50 times on four GPUs, bringing the Jacobian computing time for a fine 3D mesh from 12 min to 14 s. We regard this as an important step toward gaining interactive reconstruction times in 3D imaging, particularly when coupled in the future with acceleration of the forward problem. While we demonstrate results for EIT, these results apply to any soft-field imaging modality where the Jacobian matrix is computed with the adjoint method.