Approximation errors and model reduction with an application in optical diffusion tomography

Model reduction is often required in several applications, typically due to limited available time, computer memory or other restrictions. In problems that are related to partial differential equations, this often means that we are bound to use sparse meshes in the model for the forward problem. Conversely, if we are given more and more accurate measurements, we have to employ increasingly accurate forward problem solvers in order to exploit the information in the measurements. Optical diffusion tomography (ODT) is an example in which the typical required accuracy for the forward problem solver leads to computational times that may be unacceptable both in biomedical and industrial end applications. In this paper we review the approximation error theory and investigate the interplay between the mesh density and measurement accuracy in the case of optical diffusion tomography. We show that if the approximation errors are estimated and employed, it is possible to use mesh densities that would be unacceptable with a conventional measurement model.

[1]  J. L. Melsa,et al.  Decision and Estimation Theory , 1981, IEEE Transactions on Systems, Man, and Cybernetics.

[2]  A. Mandelbaum,et al.  Linear estimators and measurable linear transformations on a Hilbert space , 1984 .

[3]  Erkki Somersalo,et al.  Linear inverse problems for generalised random variables , 1989 .

[4]  Sabine Van Huffel,et al.  Total least squares problem - computational aspects and analysis , 1991, Frontiers in applied mathematics.

[5]  Sabine Van Huffel,et al.  The total least squares problem , 1993 .

[6]  Yves Nievergelt.  Total Least Squares: State-of-the-Art Regression in Numerical Analysis , 1994, SIAM Rev..

[7]  Richard A Williams,et al.  Process tomography : principles, techniques and applications , 1995 .

[8]  M. Schweiger,et al.  The finite element method for the propagation of light in scattering media: boundary and source conditions. , 1995, Medical physics.

[9]  Charles A. Bouman,et al.  Optical diffusion tomography by iterative- coordinate-descent optimization in a Bayesian framework , 1999 .

[10]  S. Arridge.  Optical tomography in medical imaging , 1999 .

[11]  Margaret J. Eppstein,et al.  3-D Bayesian Optical Image Reconstruction with Domain Decomposition , 2001, IEEE Trans. Medical Imaging.

[12]  Charles A. Bouman,et al.  Nonlinear multigrid algorithms for Bayesian optical diffusion tomography , 2001, IEEE Trans. Image Process..

[13]  E. Sevick-Muraca,et al.  Three-dimensional Bayesian optical image reconstruction with domain decomposition , 2001, IEEE Transactions on Medical Imaging.

[14]  Kevin J Webb,et al.  Three-dimensional Bayesian optical diffusion tomography with experimental data. , 2002, Optics letters.

[15]  Erkki Somersalo,et al.  Estimation of optical absorption in anisotropic background , 2002 .

[16]  Ville Kolehmainen,et al.  Time-series estimation of biological factors in optical diffusion tomography. , 2003, Physics in medicine and biology.

[17]  Albert Tarantola,et al.  Inverse problem theory - and methods for model parameter estimation , 2004 .

[18]  E. Somersalo,et al.  Statistical and computational inverse problems , 2004 .

[19]  S R Arridge,et al.  Recent advances in diffuse optical imaging , 2005, Physics in medicine and biology.

[20]  M. Schweiger,et al.  Gauss–Newton method for image reconstruction in diffuse optical tomography , 2005, Physics in medicine and biology.

[21]  Ville Kolehmainen,et al.  Hybrid radiative-transfer-diffusion model for optical tomography. , 2005, Applied optics.

[22]  Jari P. Kaipio,et al.  Aristotelian prior boundary conditions , 2006 .

[23]  E. Somersalo,et al.  Statistical inverse problems: discretization, model reduction and inverse crimes , 2007 .