Edge-promoting reconstruction of absorption and diffusivity in optical tomography

In optical tomography a physical body is illuminated with near-infrared light and the resulting outward photon flux is measured at the object boundary. The goal is to reconstruct internal optical properties of the body, such as absorption and diffusivity. In this work, it is assumed that the imaged object is composed of an approximately homogeneous background with clearly distinguishable embedded inhomogeneities. An algorithm for finding the maximum a posteriori estimate for the absorption and diffusion coefficients is introduced assuming an edge-preferring prior and an additive Gaussian measurement noise model. The method is based on iteratively combining a lagged diffusivity step and a linearization of the measurement model of diffuse optical tomography with priorconditioned LSQR. The performance of the reconstruction technique is tested via three-dimensional numerical experiments with simulated measurement data.

[1]  Daniela Calvetti,et al.  Hypermodels in the Bayesian imaging framework , 2008 .

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

[3]  M. Schweiger,et al.  Anisotropic diffusion regularization methods for diffuse optical tomography using edge prior information , 2006 .

[4]  Nuutti Hyvönen,et al.  Edge-Enhancing Reconstruction Algorithm for Three-Dimensional Electrical Impedance Tomography , 2014, SIAM J. Sci. Comput..

[5]  Oliver Dorn,et al.  Fréchet Derivatives for Some Bilinear Inverse Problems , 2002, SIAM J. Appl. Math..

[6]  L. Eldén.  A weighted pseudoinverse, generalized singular values, and constrained least squares problems , 1982 .

[7]  S R Arridge,et al.  3D level set reconstruction of model and experimental data in Diffuse Optical Tomography. , 2010, Optics express.

[8]  Michael Hintermüller,et al.  Nonconvex TVq-Models in Image Restoration: Analysis and a Trust-Region Regularization-Based Superlinearly Convergent Solver , 2013, SIAM J. Imaging Sci..

[9]  M. Schweiger,et al.  Instrumentation and calibration methods for the multichannel measurement of phase and amplitude in optical tomography , 2005 .

[10]  J. Wloka,et al.  Partial differential equations: Strongly elliptic differential operators and the method of variations , 1987 .

[11]  Ali Mohammad-Djafari,et al.  Please Scroll down for Article Journal of Modern Optics Bayesian Inversion for Optical Diffraction Tomography Bayesian Inversion for Optical Diffraction Tomography , 2022 .

[12]  Simon R. Arridge,et al.  Gauss–Newton reconstruction method for optical tomography using the finite element solution of the radiative transfer equation , 2008 .

[13]  Daniela Calvetti,et al.  Left and right preconditioning for electrical impedance tomography with structural information , 2012 .

[14]  Eric L. Miller,et al.  Imaging the body with diffuse optical tomography , 2001, IEEE Signal Process. Mag..

[15]  Simon R. Arridge,et al.  Iterated preconditioned LSQR method for inverse problems on unstructured grids , 2014 .

[16]  S R Arridge,et al.  Optical imaging in medicine: I. Experimental techniques , 1997, Physics in medicine and biology.

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

[18]  Michael A. Saunders,et al.  Algorithm 583: LSQR: Sparse Linear Equations and Least Squares Problems , 1982, TOMS.

[19]  K Paulsen,et al.  Instrumentation and design of a frequency-domain diffuse optical tomography imager for breast cancer detection. , 1997, Optics express.

[20]  S. Arridge,et al.  Optical imaging in medicine: II. Modelling and reconstruction , 1997, Physics in medicine and biology.

[21]  S R Arridge,et al.  Application of temporal filters to time resolved data in optical tomography , 1999, Physics in medicine and biology.

[22]  Tanja Tarvainen,et al.  Modeling boundary measurements of scattered light using the corrected diffusion approximation , 2012, Biomedical optics express.

[23]  D. Delpy,et al.  Optical Imaging in Medicine , 1998, CLEO/Europe Conference on Lasers and Electro-Optics.

[24]  A. Mohammad-Djafari,et al.  Optical diffraction tomography within a variational Bayesian framework , 2012 .

[25]  S. Arridge,et al.  Nonuniqueness in diffusion-based optical tomography. , 1998, Optics letters.

[26]  Britton Chance,et al.  Diffuse optical tomography with physiological and spatial a priori constraints , 2004, Physics in medicine and biology.

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

[28]  Tanja Tarvainen,et al.  MARGINALIZATION OF UNINTERESTING DISTRIBUTED PARAMETERS IN INVERSE PROBLEMS-APPLICATION TO DIFFUSE OPTICAL TOMOGRAPHY , 2011 .

[29]  Michael A. Saunders,et al.  LSQR: An Algorithm for Sparse Linear Equations and Sparse Least Squares , 1982, TOMS.

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

[31]  Curtis R. Vogel,et al.  Iterative Methods for Total Variation Denoising , 1996, SIAM J. Sci. Comput..

[32]  J. Hilgers.  Erratum: On the Equivalence of Regularization and Certain Reproducing Kernel Hilbert Space Approaches for Solving First Kind Problems , 1976 .

[33]  D. Calvetti,et al.  Priorconditioners for linear systems , 2005 .

[34]  Timothy C Zhu,et al.  A review of in‐vivo optical properties of human tissues and its impact on PDT , 2011, Journal of biophotonics.

[35]  L. Rudin,et al.  Nonlinear total variation based noise removal algorithms , 1992 .

[36]  L. Fajardo,et al.  Differentiation of cysts from solid tumors in the breast with diffuse optical tomography. , 2004, Academic radiology.

[37]  M. Schweiger,et al.  Optical tomography of a realistic neonatal head phantom. , 2003, Applied optics.

[38]  Simon R. Arridge,et al.  A Bayesian approach to spectral quantitative photoacoustic tomography , 2014 .

[39]  Nuutti Hyvönen,et al.  Application of a weaker formulation of the factorization method to the characterization of absorbing inclusions in optical tomography , 2005 .

[40]  B. Harrach.  On uniqueness in diffuse optical tomography , 2009 .

[41]  Jitendra Malik,et al.  Scale-Space and Edge Detection Using Anisotropic Diffusion , 1990, IEEE Trans. Pattern Anal. Mach. Intell..

[42]  Daniela Calvetti,et al.  Preconditioned iterative methods for linear discrete ill-posed problems from a Bayesian inversion perspective , 2007 .

[43]  E. Somersalo,et al.  Approximation errors and model reduction with an application in optical diffusion tomography , 2006 .