706 research outputs found
Expectation Propagation for Nonlinear Inverse Problems -- with an Application to Electrical Impedance Tomography
In this paper, we study a fast approximate inference method based on
expectation propagation for exploring the posterior probability distribution
arising from the Bayesian formulation of nonlinear inverse problems. It is
capable of efficiently delivering reliable estimates of the posterior mean and
covariance, thereby providing an inverse solution together with quantified
uncertainties. Some theoretical properties of the iterative algorithm are
discussed, and the efficient implementation for an important class of problems
of projection type is described. The method is illustrated with one typical
nonlinear inverse problem, electrical impedance tomography with complete
electrode model, under sparsity constraints. Numerical results for real
experimental data are presented, and compared with that by Markov chain Monte
Carlo. The results indicate that the method is accurate and computationally
very efficient.Comment: Journal of Computational Physics, to appea
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
Quadrature formula for computed tomography
AbstractWe give a bivariate analog of the Micchelli–Rivlin quadrature for computing the integral of a function over the unit disk using its Radon projections
Quantitative photoacoustic imaging in radiative transport regime
The objective of quantitative photoacoustic tomography (QPAT) is to
reconstruct optical and thermodynamic properties of heterogeneous media from
data of absorbed energy distribution inside the media. There have been
extensive theoretical and computational studies on the inverse problem in QPAT,
however, mostly in the diffusive regime. We present in this work some numerical
reconstruction algorithms for multi-source QPAT in the radiative transport
regime with energy data collected at either single or multiple wavelengths. We
show that when the medium to be probed is non-scattering, explicit
reconstruction schemes can be derived to reconstruct the absorption and the
Gruneisen coefficients. When data at multiple wavelengths are utilized, we can
reconstruct simultaneously the absorption, scattering and Gruneisen
coefficients. We show by numerical simulations that the reconstructions are
stable.Comment: 40 pages, 13 figure
- …