Validation of a finite-element solution for electrical impedance tomography in an anisotropic medium

Electrical impedance tomography is an imaging method, with which volumetric images of conductivity are produced by injecting electrical current and measuring boundary voltages. It has the potential to become a portable non-invasive medical imaging technique. Until now, implementations have neglected anisotropy even though human tissues such as bone, muscle and brain white matter are markedly anisotropic. We present a numerical solution using the finite-element method that has been modified for modelling anisotropic conductive media. It was validated in an anisotropic domain against an analytical solution in an isotropic medium after the isotropic domain was diffeomorphically transformed into an anisotropic one. Convergence of the finite element to the analytical solution was verified by showing that the finite-element error norm decreased linearly related to the finite-element size, as the mesh density increased, for the simplified case of Laplace's equation in a cubic domain with a Dirichlet boundary condition

A direct reconstruction algorithm for the anisotropic inverse conductivity problem based on Calderon's method in the plane

A direct reconstruction algorithm based on Calderon's linearization method for the reconstruction of isotropic conductivities is proposed for anisotropic conductivities in two-dimensions. To overcome the non-uniqueness of the anisotropic inverse conductivity problem, the entries of the unperturbed anisotropic tensors are assumed known a priori, and it remains to reconstruct the multiplicative scalar field. The quasi-conformal map in the plane facilitates the Calderon-based approach for anisotropic conductivities. The method is demonstrated on discontinuous radially symmetric conductivities of high and low contrast.Peer reviewe

A Direct Reconstruction Method for Anisotropic Electrical Impedance Tomography

A novel computational, non-iterative and noise-robust reconstruction method is introduced for the planar anisotropic inverse conductivity problem. The method is based on bypassing the unstable step of the reconstruction of the values of the isothermal coordinates on the boundary of the domain. Non-uniqueness of the inverse problem is dealt with by recovering the unique isotropic conductivity that can be achieved as a deformation of the measured anisotropic conductivity by \emph{isothermal coordinates}. The method shows how isotropic D-bar reconstruction methods have produced reasonable and informative reconstructions even when used on EIT data known to come from anisotropic media, and when the boundary shape is not known precisely. Furthermore, the results pave the way for regularized anisotropic EIT. Key aspects of the approach involve D-bar methods and inverse scattering theory, complex geometrical optics solutions, and quasi-conformal mapping techniques.Comment: 30 pages, 8 figure

Regularisation methods for imaging from electrical measurements

In Electrical Impedance Tomography the conductivity of an object is estimated from boundary measurements. An array of electrodes is attached to the surface of the object and current stimuli are applied via these electrodes. The resulting voltages are measured. The process of estimating the conductivity as a function of space inside the object from voltage measurements at the surface is called reconstruction. Mathematically the ElT reconstruction is a non linear inverse problem, the stable solution of which requires regularisation methods. Most common regularisation methods impose that the reconstructed image should be smooth. Such methods confer stability to the reconstruction process, but limit the capability of describing sharp variations in the sought parameter. In this thesis two new methods of regularisation are proposed. The first method, Gallssian anisotropic regularisation, enhances the reconstruction of sharp conductivity changes occurring at the interface between a contrasting object and the background. As such changes are step changes, reconstruction with traditional smoothing regularisation techniques is unsatisfactory. The Gaussian anisotropic filtering works by incorporating prior structural information. The approximate knowledge of the shapes of contrasts allows us to relax the smoothness in the direction normal to the expected boundary. The construction of Gaussian regularisation filters that express such directional properties on the basis of the structural information is discussed, and the results of numerical experiments are analysed. The method gives good results when the actual conductivity distribution is in accordance with the prior information. When the conductivity distribution violates the prior information the method is still capable of properly locating the regions of contrast. The second part of the thesis is concerned with regularisation via the total variation functional. This functional allows the reconstruction of discontinuous parameters. The properties of the functional are briefly introduced, and an application in inverse problems in image denoising is shown. As the functional is non-differentiable, numerical difficulties are encountered in its use. The aim is therefore to propose an efficient numerical implementation for application in ElT. Several well known optimisation methods arc analysed, as possible candidates, by theoretical considerations and by numerical experiments. Such methods are shown to be inefficient. The application of recent optimisation methods called primal- dual interior point methods is analysed be theoretical considerations and by numerical experiments, and an efficient and stable algorithm is developed. Numerical experiments demonstrate the capability of the algorithm in reconstructing sharp conductivity profiles

Development of the VHP-Female Full-Body Computational Model and Its Applications for Biomedical Electromagnetic Modeling

Computational modeling offers better insight into a wide range of bioelectrical and biomechanical problems with improved tools for the design of medical devices and the diagnosis of pathologies. Electromagnetic modeling at low and high frequencies is particularly necessary. Modeling electromagnetic, structural, thermal, and acoustic response of the human body to different internal and external stimuli is limited by the availability of numerically efficient computational human models. This study describes the development to date of a computational full-body human model - Visible Human Project (VHP) - Female Model. Its unique feature is full compatibility both with MATLAB and specialized FEM computational software packages such as ANSYS HFSS/Maxwell 3D. This study also describes progress made to date in using the newly developed tools for segmentation. A visualization tool is implemented within MATLAB and is based on customized version of the constrained 2D Delaunay triangulation method for intersecting objects. This thesis applies a VHP - Female Model to a specific application, transcranial Direct Current Stimulation (tDCS). Transcranial Direct Current Stimulation has been beneficial in the stimulation of cortical activity and treatment of neurological disorders in humans. The placement of electrodes, which is cephalic versus extracephalic montages, is studied for optimal targeting of currents for a given functional area. Given the difficulty of obtaining in vivo measurements of current density, modeling of conventional and alternative electrode montages via the FEM has been utilized to provide insight into the tDCS montage performance. An insight into future work and potential areas of research, such as study of bone quality have been presented too

The use of charge -charge correlation in impedance measurements: A test of the EPET method

It is well known that biological tissues possess impedance properties that might be useful in medical diagnostics and treatment. Electrical Impedance Tomography (EIT) images internal electrical properties by using numerical methods to solve Laplace\u27s differential equation. The indirect reconstruction method (IRM), a common method in application, predicts internal electrical property distribution by iteratively computing a forward and inverse solution. This approach reduces the non-linear Laplace\u27s equation into a poorly conditioned series of linear equations, which are solved simultaneously. This method suffers from high computational effort and is susceptible to prediction errors that stem from measurement noise.;As an alternative to Laplace\u27s differential equation, this research applies the quasi-static approximation, Dirichlet boundary conditions and a rectangular shaped domain (with corresponding Green\u27s function for Cartesian coordinates) to solve the integral form of Poisson\u27s equation (Green\u27s 2nd identity). The result is the charge-charge correlation method (CCCM), a well-conditioned relationship between static charge build-up at internal structures and induced domain boundary charge build-up (which corresponds to measured boundary current). The CCCM is applied in a reconstruction technique called Electrical Property Enhanced Tomography (EPET). While related to the existing impedance imaging methods, EPET does not attempt to create the image with the electrical data but rather adds electrical property information to an existing conventional imaging modality (CT or MI) and, in fact, requires the data from the other modality to locate the position of internal structures in the object. Predicted electrical properties are then superimposed over the a priori structural image to yield the electrical property distribution.;To test the feasibility of the CCCM, experiments using agar media placed in a saline bath were performed. The position, size and conductivity of the agar were varied and the CCCM applied to predict the conductivities from external boundary current measurements. Predicted conductivities yielded relative errors less than 10%, results that are equal to or better than the IRM. Additionally, CCCM was able to compute these results with a 104 improvement in speed over the IRM

Variation in Reported Human Head Tissue Electrical Conductivity Values

Electromagnetic source characterisation requires accurate volume conductor models representing head geometry and the electrical conductivity field. Head tissue conductivity is often assumed from previous literature, however, despite extensive research, measurements are inconsistent. A meta-analysis of reported human head electrical conductivity values was therefore conducted to determine significant variation and subsequent influential factors. Of 3121 identified publications spanning three databases, 56 papers were included in data extraction. Conductivity values were categorised according to tissue type, and recorded alongside methodology, measurement condition, current frequency, tissue temperature, participant pathology and age. We found variation in electrical conductivity of the whole-skull, the spongiform layer of the skull, isotropic, perpendicularly- and parallelly-oriented white matter (WM) and the brain-to-skull-conductivity ratio (BSCR) could be significantly attributed to a combination of differences in methodology and demographics. This large variation should be acknowledged, and care should be taken when creating volume conductor models, ideally constructing them on an individual basis, rather than assuming them from the literature. When personalised models are unavailable, it is suggested weighted average means from the current meta-analysis are used. Assigning conductivity as: 0.41 S/m for the scalp, 0.02 S/m for the whole skull, or when better modelled as a three-layer skull 0.048 S/m for the spongiform layer, 0.007 S/m for the inner compact and 0.005 S/m for the outer compact, as well as 1.71 S/m for the CSF, 0.47 S/m for the grey matter, 0.22 S/m for WM and 50.4 for the BSCR

DICOM for EIT

With EIT starting to be used in routine clinical practice [1], it important that the clinically relevant information is portable between hospital data management systems. DICOM formats are widely used clinically and cover many imaging modalities, though not specifically EIT. We describe how existing DICOM specifications, can be repurposed as an interim solution, and basis from which a consensus EIT DICOM ‘Supplement’ (an extension to the standard) can be writte