679 research outputs found
Deep D-Bar: Real-Time Electrical Impedance Tomography Imaging With Deep Neural Networks
The mathematical problem for electrical impedance tomography (EIT) is a highly nonlinear ill-posed inverse problem requiring carefully designed reconstruction procedures to ensure reliable image generation. D-bar methods are based on a rigorous mathematical analysis and provide robust direct reconstructions by using a low-pass filtering of the associated nonlinear Fourier data. Similarly to low-pass filtering of linear Fourier data, only using low frequencies in the image recovery process results in blurred images lacking sharp features, such as clear organ boundaries. Convolutional neural networks provide a powerful framework for post-processing such convolved direct reconstructions. In this paper, we demonstrate that these CNN techniques lead to sharp and reliable reconstructions even for the highly nonlinear inverse problem of EIT. The network is trained on data sets of simulated examples and then applied to experimental data without the need to perform an additional transfer training. Results for absolute EIT images are presented using experimental EIT data from the ACT4 and KIT4 EIT systems
Discrete Geometric Structures in Homogenization and Inverse Homogenization with application to EIT
We introduce a new geometric approach for the homogenization and inverse
homogenization of the divergence form elliptic operator with rough conductivity
coefficients in dimension two. We show that conductivity
coefficients are in one-to-one correspondence with divergence-free matrices and
convex functions over the domain . Although homogenization is a
non-linear and non-injective operator when applied directly to conductivity
coefficients, homogenization becomes a linear interpolation operator over
triangulations of when re-expressed using convex functions, and is a
volume averaging operator when re-expressed with divergence-free matrices.
Using optimal weighted Delaunay triangulations for linearly interpolating
convex functions, we obtain an optimally robust homogenization algorithm for
arbitrary rough coefficients. Next, we consider inverse homogenization and show
how to decompose it into a linear ill-posed problem and a well-posed non-linear
problem. We apply this new geometric approach to Electrical Impedance
Tomography (EIT). It is known that the EIT problem admits at most one isotropic
solution. If an isotropic solution exists, we show how to compute it from any
conductivity having the same boundary Dirichlet-to-Neumann map. It is known
that the EIT problem admits a unique (stable with respect to -convergence)
solution in the space of divergence-free matrices. As such we suggest that the
space of convex functions is the natural space in which to parameterize
solutions of the EIT problem
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
The Linearized Inverse Problem in Multifrequency Electrical Impedance Tomography
This paper provides an analysis of the linearized inverse problem in
multifrequency electrical impedance tomography. We consider an isotropic
conductivity distribution with a finite number of unknown inclusions with
different frequency dependence, as is often seen in biological tissues. We
discuss reconstruction methods for both fully known and partially known
spectral profiles, and demonstrate in the latter case the successful employment
of difference imaging. We also study the reconstruction with an imperfectly
known boundary, and show that the multifrequency approach can eliminate
modeling errors and recover almost all inclusions. In addition, we develop an
efficient group sparse recovery algorithm for the robust solution of related
linear inverse problems. Several numerical simulations are presented to
illustrate and validate the approach.Comment: 25 pp, 11 figure
Graph- and finite element-based total variation models for the inverse problem in diffuse optical tomography
Total variation (TV) is a powerful regularization method that has been widely
applied in different imaging applications, but is difficult to apply to diffuse
optical tomography (DOT) image reconstruction (inverse problem) due to complex
and unstructured geometries, non-linearity of the data fitting and
regularization terms, and non-differentiability of the regularization term. We
develop several approaches to overcome these difficulties by: i) defining
discrete differential operators for unstructured geometries using both finite
element and graph representations; ii) developing an optimization algorithm
based on the alternating direction method of multipliers (ADMM) for the
non-differentiable and non-linear minimization problem; iii) investigating
isotropic and anisotropic variants of TV regularization, and comparing their
finite element- and graph-based implementations. These approaches are evaluated
on experiments on simulated data and real data acquired from a tissue phantom.
Our results show that both FEM and graph-based TV regularization is able to
accurately reconstruct both sparse and non-sparse distributions without the
over-smoothing effect of Tikhonov regularization and the over-sparsifying
effect of L regularization. The graph representation was found to
out-perform the FEM method for low-resolution meshes, and the FEM method was
found to be more accurate for high-resolution meshes.Comment: 24 pages, 11 figures. Reviced version includes revised figures and
improved clarit
A Data-Driven Edge-Preserving D-bar Method for Electrical Impedance Tomography
In Electrical Impedance Tomography (EIT), the internal conductivity of a body
is recovered via current and voltage measurements taken at its surface. The
reconstruction task is a highly ill-posed nonlinear inverse problem, which is
very sensitive to noise, and requires the use of regularized solution methods,
of which D-bar is the only proven method. The resulting EIT images have low
spatial resolution due to smoothing caused by low-pass filtered regularization.
In many applications, such as medical imaging, it is known \emph{a priori} that
the target contains sharp features such as organ boundaries, as well as
approximate ranges for realistic conductivity values. In this paper, we use
this information in a new edge-preserving EIT algorithm, based on the original
D-bar method coupled with a deblurring flow stopped at a minimal data
discrepancy. The method makes heavy use of a novel data fidelity term based on
the so-called {\em CGO sinogram}. This nonlinear data step provides superior
robustness over traditional EIT data formats such as current-to-voltage
matrices or Dirichlet-to-Neumann operators, for commonly used current patterns.Comment: 24 pages, 11 figure
High-order regularized regression in Electrical Impedance Tomography
We present a novel approach for the inverse problem in electrical impedance
tomography based on regularized quadratic regression. Our contribution
introduces a new formulation for the forward model in the form of a nonlinear
integral transform, that maps changes in the electrical properties of a domain
to their respective variations in boundary data. Using perturbation theory the
transform is approximated to yield a high-order misfit unction which is then
used to derive a regularized inverse problem. In particular, we consider the
nonlinear problem to second-order accuracy, hence our approximation method
improves upon the local linearization of the forward mapping. The inverse
problem is approached using Newton's iterative algorithm and results from
simulated experiments are presented. With a moderate increase in computational
complexity, the method yields superior results compared to those of regularized
linear regression and can be implemented to address the nonlinear inverse
problem
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
- …