Biomedical subjects

William R B Lionheart

Publications and source records attributed to William R B Lionheart.

8 recordsLinked to original sources

PubMed2006

Uses and abuses of EIDORS: an extensible software base for EIT.

EIDORS is an open source software suite for image reconstruction in electrical impedance tomography and diffuse optical tomography, designed to facilitate collaboration, testing and new research in these fields. This paper describes recent work to redesign the software structure in order to simplify its use and provide a uniform interface, permitting easier modification and customization. We describe the key features of this software, followed by examples of its use. One general issue with inverse problem software is the difficulty of correctly implementing algorithms and the consequent ease with which subtle numerical bugs can be inadvertently introduced. EIDORS helps with this issue, by allowing sharing and reuse of well-documented and debugged software. On the other hand, since EIDORS is designed to facilitate use by non-specialists, its use may inadvertently result in such numerical errors. In order to address this issue, we develop a list of ways in which such errors with inverse problems (which we refer to as 'cheats') may occur. Our hope is that such an overview may assist authors of software to avoid such implementation issues.

Algorithms↗

PubMed2006

A narrow-band level set method applied to EIT in brain for cryosurgery monitoring.

In this paper, we investigate the feasibility of applying a novel level set reconstruction technique to electrical imaging of the human brain. We focus particularly on the potential application of electrical impedance tomography (EIT) to cryosurgery monitoring. In this application, cancerous tissue is treated by a local freezing technique using a small needle-like cryosurgery probe. The interface between frozen and nonfrozen tissue can be expected to have a relatively high contrast in conductivity and we treat the inverse problem of locating and monitoring this interface during the treatment. A level set method is used as a powerful and flexible tool for tracking the propagating interfaces during the monitoring process. For calculating sensitivities and the Jacobian when deforming the interfaces we employ an adjoint formula rather than a direct differentiation technique. In particular, we are using a narrow-band technique for this procedure. This combination of an adjoint technique and a narrow-band technique for calculating Jacobians results in a computationally efficient and extremely fast method for solving the inverse problem. Moreover, due to the reduced number of unknowns in each step of the narrow-band approach compared to a pixel- or voxel-based technique, our reconstruction scheme tends to be much more stable. We demonstrate that our new method also outperforms its pixel-/voxel-based counterparts in terms of image quality in this application.

Brain Neoplasms↗

PubMed2006

Absolute conductivity reconstruction in magnetic induction tomography using a nonlinear method.

Magnetic induction tomography (MIT) attempts to image the electrical and magnetic characteristics of a target using impedance measurement data from pairs of excitation and detection coils. This inverse eddy current problem is nonlinear and also severely ill posed so regularization is required for a stable solution. A regularized Gauss-Newton algorithm has been implemented as a nonlinear, iterative inverse solver. In this algorithm, one needs to solve the forward problem and recalculate the Jacobian matrix for each iteration. The forward problem has been solved using an edge based finite element method for magnetic vector potential A and electrical scalar potential V, a so called A, A - V formulation. A theoretical study of the general inverse eddy current problem and a derivation, paying special attention to the boundary conditions, of an adjoint field formula for the Jacobian is given. This efficient formula calculates the change in measured induced voltage due to a small perturbation of the conductivity in a region. This has the advantage that it involves only the inner product of the electric fields when two different coils are excited, and these are convenient computationally. This paper also shows that the sensitivity maps change significantly when the conductivity distribution changes, demonstrating the necessity for a nonlinear reconstruction algorithm. The performance of the inverse solver has been examined and results presented from simulated data with added noise.

Algorithms↗

PubMed2005

Reconstruction of spatially inhomogeneous dielectric tensors through optical tomography.

A method to reconstruct weakly anisotropic inhomogeneous dielectric tensors inside a transparent medium is proposed. The mathematical theory of integral geometry is cast into a workable framework that allows the full determination of dielectric tensor fields by scalar Radon inversions of the polarization transformation data obtained from six planar tomographic scanning cycles. Furthermore, a careful derivation of the usual equations of integrated photoelasticity in terms of heuristic length scales of the material inhomogeneity and anisotropy is provided, resulting in a self-contained account about the reconstruction of arbitrary three-dimensional, weakly anisotropic dielectric tensor fields.

Algorithms↗

PubMed2004

EIT reconstruction algorithms: pitfalls, challenges and recent developments.

We review developments, issues and challenges in electrical impedance tomography (EIT) for the 4th Conference on Biomedical Applications of Electrical Impedance Tomography, held at Manchester in 2003. We focus on the necessity for three-dimensional data collection and reconstruction, efficient solution of the forward problem, and both present and future reconstruction algorithms. We also suggest common pitfalls or 'inverse crimes' to avoid.

Algorithms↗

PubMed2003

Uniqueness and reconstruction in magnetic resonance-electrical impedance tomography (MR-EIT).

Magnetic resonance-electrical impedance tomography (MR-EIT) was first proposed in 1992. Since then various reconstruction algorithms have been suggested and applied. These algorithms use peripheral voltage measurements and internal current density measurements in different combinations. In this study the problem of MR-EIT is treated as a hyperbolic system of first-order partial differential equations, and three numerical methods are proposed for its solution. This approach is not utilized in any of the algorithms proposed earlier. The numerical solution methods are integration along equipotential surfaces (method of characteristics), integration on a Cartesian grid, and inversion of a system matrix derived by a finite difference formulation. It is shown that if some uniqueness conditions are satisfied, then using at least two injected current patterns, resistivity can be reconstructed apart from a multiplicative constant. This constant can then be identified using a single voltage measurement. The methods proposed are direct, non-iterative, and valid and feasible for 3D reconstructions. They can also be used to easily obtain slice and field-of-view images from a 3D object. 2D simulations are made to illustrate the performance of the algorithms.

Algorithms↗

PubMed2002

Krylov subspace iterative techniques: on the detection of brain activity with electrical impedance tomography.

In this paper, we review some numerical techniques based on the linear Krylov subspace iteration that can be used for the efficient calculation of the forward and the inverse electrical impedance tomography problems. Exploring their computational advantages in solving large-scale systems of equations, we specifically address their implementation in reconstructing localized impedance changes occurring within the human brain. If the conductivity of the head tissues is assumed to be real, the pre-conditioned conjugate gradients (PCGs) algorithm can be used to calculate efficiently the approximate forward solution to a given error tolerance. The performance and the regularizing properties of the PCG iteration for solving ill-conditioned systems of equations (PCGNs) is then explored, and a suitable preconditioning matrix is suggested in order to enhance its convergence rate. For image reconstruction, the nonlinear inverse problem is considered. Based on the Gauss-Newton method for solving nonlinear problems we have developed two algorithms that implement the PCGN iteration to calculate the linear step solution. Using an anatomically detailed model of the human head and a specific scalp electrode arrangement, images of a simulated impedance change inside brain's white matter have been reconstructed.

Algorithms↗

PubMed2002

Generation of anisotropic-smoothness regularization filters for EIT.

In the inverse conductivity problem, as in any ill-posed inverse problem, regularization techniques are necessary in order to stabilize inversion. A common way to implement regularization in electrical impedance tomography is to use Tikhonov regularization. The inverse problem is formulated as a minimization of two terms: the mismatch of the measurements against the model, and the regularization functional. Most commonly, differential operators are used as regularization functionals, leading to smooth solutions. Whenever the imaged region presents discontinuities in the conductivity distribution, such as interorgan boundaries, the smoothness prior is not consistent with the actual situation. In these cases, the reconstruction is enhanced by relaxing the smoothness constraints in the direction normal to the discontinuity. In this paper, we derive a method for generating Gaussian anisotropic regularization filters. The filters are generated on the basis of the prior structural information, allowing a better reconstruction of conductivity profiles matching these priors. When incorporating prior information into a reconstruction algorithm, the risk is of biasing the inverse solutions toward the assumed distributions. Simulations show that, with a careful selection of the regularization parameters, the reconstruction algorithm is still able to detect conductivities patterns that violate the prior information. A generalized singular-value decomposition analysis of the effects of the anisotropic filters on regularization is presented in the last sections of the paper.

Algorithms↗