Accelerating nonlinear reconstruction in laminar optical tomography by use of recursive SVD inversion.

Image reconstruction in the most model-based biophotonic imaging modalities essentially poses an ill-posed nonlinear inverse problem, which has been effectively tackled in the diffusion-approximation-satisfied scenarios such as diffuse optical tomography. Nevertheless, a nonlinear implementation in high-resolution laminar optical tomography (LOT) is normally computationally-costly due to its strong dependency on a dense source-detector configuration and a physically-rigorous photon-transport model. To circumvent the adversity, we herein propose a practical nonlinear LOT approach to the absorption reconstruction. The scheme takes advantage of the numerical stability of the singular value decomposition (SVD) for the ill-posed linear inversion, and is accelerated by adopting an explicitly recursive strategy for the time-consuming repeated SVD inversion, which is based on a scaled expression of the sensitivity matrix. Experiments demonstrate that the proposed methodology can perform as well as the traditional nonlinear one, while the computation time of the former is merely 26.27% of the later on average.

[1]  Feng Gao,et al.  Full domain-decomposition scheme for diffuse optical tomography of large-sized tissues with a combined CPU and GPU parallelization. , 2014, Applied optics.

[2]  Cameron Musco,et al.  Randomized Block Krylov Methods for Stronger and Faster Approximate Singular Value Decomposition , 2015, NIPS.

[3]  Saurabh Gupta,et al.  Singular value decomposition based computationally efficient algorithm for rapid dynamic near-infrared diffuse optical tomography. , 2009, Medical physics.

[4]  Baohong Yuan,et al.  A system for high-resolution depth-resolved optical imaging of fluorescence and absorption contrast. , 2009, The Review of scientific instruments.

[5]  S Arridge,et al.  A gradient-based optimisation scheme foroptical tomography. , 1998, Optics express.

[6]  Hamid Dehghani,et al.  An efficient Jacobian reduction method for diffuse optical image reconstruction. , 2007, Optics express.

[7]  Huabei Jiang,et al.  Mesh-based enhancement schemes in diffuse optical tomography. , 2003, Medical physics.

[8]  Yukio Yamada,et al.  Diffuse optical tomography: Present status and its future , 2014 .

[9]  Simon R. Arridge,et al.  Differential imaging in heterogeneous media: limitations of linearization assumptions in optical tomography , 2001, SPIE BiOS.

[10]  A. Joshi,et al.  Adaptive finite element based tomography for fluorescence optical imaging in tissue. , 2004, Optics express.

[12]  D. Hanahan,et al.  Neovascularization is associated with a switch to the export of bFGF in the multistep development of fibrosarcoma , 1991, Cell.

[13]  M. Schweiger,et al.  Linear and nonlinear reconstruction for optical tomography of phantoms with nonscattering regions. , 2005, Applied optics.

[14]  David A Boas,et al.  Laminar optical tomography: demonstration of millimeter-scale depth-resolved imaging in turbid media. , 2004, Optics letters.

[15]  Huabei Jiang,et al.  Quantitative photoacoustic tomography based on the radiative transfer equation. , 2009, Optics letters.

[16]  Haowei Huang,et al.  GPU accelerated simulations of 3D deterministic particle transport using discrete ordinates method , 2011, J. Comput. Phys..

[17]  Hanli Liu,et al.  Comprehensive investigation of three-dimensional diffuse optical tomography with depth compensation algorithm. , 2010, Journal of biomedical optics.

[18]  Luis Mateus Rocha,et al.  Singular value decomposition and principal component analysis , 2003 .

[19]  David A Boas,et al.  Monte Carlo simulation of photon migration in 3D turbid media accelerated by graphics processing units. , 2009, Optics express.

[20]  Nirmala Ramanujam,et al.  Scaling method for fast Monte Carlo simulation of diffuse reflectance spectra from multilayered turbid media. , 2007, Journal of the Optical Society of America. A, Optics, image science, and vision.

[21]  J. Spanier,et al.  Perturbation Monte Carlo methods to solve inverse photon migration problems in heterogeneous tissues. , 2001, Optics letters.

[22]  S. Jacques Optical properties of biological tissues: a review , 2013, Physics in medicine and biology.

[23]  Alexander D. Klose,et al.  Gradient-based iterative image reconstruction scheme for time-resolved optical tomography , 1999, IEEE Transactions on Medical Imaging.

[24]  Timothy J. Muldoon,et al.  Analysis of skin lesions using laminar optical tomography , 2012, Biomedical optics express.

[25]  Yukio Yamada,et al.  Hybrid model of light propagation in random media based on the time-dependent radiative transfer and diffusion equations , 2014 .

[26]  Vasilis Ntziachristos,et al.  Three‐dimensional multispectral optoacoustic mesoscopy reveals melanin and blood oxygenation in human skin in vivo , 2016, Journal of biophotonics.

[27]  Mengyu Jia,et al.  Image reconstruction method for laminar optical tomography with only a single Monte-Carlo simulation , 2014 .

[28]  Feng Gao,et al.  A modified laminar optical tomography system with small dip-angle and the initial validation , 2016, SPIE BiOS.

[29]  D Boas,et al.  Transport-based image reconstruction in turbid media with small source-detector separations. , 2000, Optics letters.

[30]  S. Arridge,et al.  Optical tomography: forward and inverse problems , 2009, 0907.2586.

[31]  Vasilis Ntziachristos,et al.  Advances in real-time multispectral optoacoustic imaging and its applications , 2015, Nature Photonics.

[32]  V. Venugopalan,et al.  Radiative transport in the delta-P1 approximation for semi-infinite turbid media. , 2008, Medical physics.

[33]  Manolis I. A. Lourakis,et al.  Estimating the Jacobian of the Singular Value Decomposition: Theory and Applications , 2000, ECCV.

[34]  Joseph P Culver,et al.  Fast and efficient image reconstruction for high density diffuse optical imaging of the human brain. , 2015, Biomedical optics express.