1,459 research outputs found
An Analysis of Finite Element Approximation in Electrical Impedance Tomography
We present a finite element analysis of electrical impedance tomography for
reconstructing the conductivity distribution from electrode voltage
measurements by means of Tikhonov regularization. Two popular choices of the
penalty term, i.e., -norm smoothness penalty and total variation
seminorm penalty, are considered. A piecewise linear finite element method is
employed for discretizing the forward model, i.e., the complete electrode
model, the conductivity, and the penalty functional. The convergence of the
finite element approximations for the Tikhonov model on both polyhedral and
smooth curved domains is established. This provides rigorous justifications for
the ad hoc discretization procedures in the literature.Comment: 20 page
EIT Reconstruction Algorithms: Pitfalls, Challenges and Recent Developments
We review developments, issues and challenges in Electrical Impedance
Tomography (EIT), for the 4th Workshop on Biomedical Applications of EIT,
Manchester 2003. We focus on the necessity for three dimensional data
collection and reconstruction, efficient solution of the forward problem and
present and future reconstruction algorithms. We also suggest common pitfalls
or ``inverse crimes'' to avoid.Comment: A review paper for the 4th Workshop on Biomedical Applications of
EIT, Manchester, UK, 200
Lipschitz stability for the electrostatic inverse boundary value problem with piecewise linear conductivities
We consider the electrostatic inverse boundary value problem also known as
electrical impedance tomography (EIT) for the case where the conductivity is a
piecewise linear function on a domain and we show
that a Lipschitz stability estimate for the conductivity in terms of the local
Dirichlet-to-Neumann map holds true.Comment: 28 pages. arXiv admin note: text overlap with arXiv:1405.047
A Partially Reflecting Random Walk on Spheres Algorithm for Electrical Impedance Tomography
In this work, we develop a probabilistic estimator for the voltage-to-current
map arising in electrical impedance tomography. This novel so-called partially
reflecting random walk on spheres estimator enables Monte Carlo methods to
compute the voltage-to-current map in an embarrassingly parallel manner, which
is an important issue with regard to the corresponding inverse problem. Our
method uses the well-known random walk on spheres algorithm inside subdomains
where the diffusion coefficient is constant and employs replacement techniques
motivated by finite difference discretization to deal with both mixed boundary
conditions and interface transmission conditions. We analyze the global bias
and the variance of the new estimator both theoretically and experimentally. In
a second step, the variance is considerably reduced via a novel control variate
conditional sampling technique
A Direct D-Bar Method for Partial Boundary Data Electrical Impedance Tomography With a Priori Information
Electrical Impedance Tomography (EIT) is a non-invasive imaging modality that uses surface electrical measurements to determine the internal conductivity of a body. The mathematical formulation of the EIT problem is a nonlinear and severely ill-posed inverse problem for which direct D-bar methods have proved useful in providing noise-robust conductivity reconstructions. Recent advances in D-bar methods allow for conductivity reconstructions using EIT measurement data from only part of the domain (e.g., a patient lying on their back could be imaged using only data gathered on the accessible part of the body). However, D-bar reconstructions suffer from a loss of sharp edges due to a nonlinear low-pass filtering of the measured data, and this problem becomes especially marked in the case of partial boundary data. Including a priori data directly into the D-bar solution method greatly enhances the spatial resolution, allowing for detection of underlying pathologies or defects, even with no assumption of their presence in the prior. This work combines partial data D-bar with a priori data, allowing for noise-robust conductivity reconstructions with greatly improved spatial resolution. The method is demonstrated to be effective on noisy simulated EIT measurement data simulating both medical and industrial imaging scenarios
On moduli of rings and quadrilaterals: algorithms and experiments
Moduli of rings and quadrilaterals are frequently applied in geometric
function theory, see e.g. the Handbook by K\"uhnau. Yet their exact values are
known only in a few special cases. Previously, the class of planar domains with
polygonal boundary has been studied by many authors from the point of view of
numerical computation. We present here a new -FEM algorithm for the
computation of moduli of rings and quadrilaterals and compare its accuracy and
performance with previously known methods such as the Schwarz-Christoffel
Toolbox of Driscoll and Trefethen. We also demonstrate that the -FEM
algorithm applies to the case of non-polygonal boundary and report results with
concrete error bounds
Spatially structured oscillations in a two-dimensional excitatory neuronal network with synaptic depression
We study the spatiotemporal dynamics of a two-dimensional excitatory neuronal network with synaptic depression. Coupling between populations of neurons is taken to be nonlocal, while depression is taken to be local and presynaptic. We show that the network supports a wide range of spatially structured oscillations, which are suggestive of phenomena seen in cortical slice experiments and in vivo. The particular form of the oscillations depends on initial conditions and the level of background noise. Given an initial, spatially localized stimulus, activity evolves to a spatially localized oscillating core that periodically emits target waves. Low levels of noise can spontaneously generate several pockets of oscillatory activity that interact via their target patterns. Periodic activity in space can also organize into spiral waves, provided that there is some source of rotational symmetry breaking due to external stimuli or noise. In the high gain limit, no oscillatory behavior exists, but a transient stimulus can lead to a single, outward propagating target wave
- …