OSA's Digital Library

Applied Optics

Applied Optics


  • Editor: Joseph N. Mait
  • Vol. 49, Iss. 16 — Jun. 1, 2010
  • pp: 3111–3126

Combined domain-decomposition and matrix-decomposition scheme for large-scale diffuse optical tomography

Fang Yang, Feng Gao, Pingqiao Ruan, and Huijuan Zhao  »View Author Affiliations

Applied Optics, Vol. 49, Issue 16, pp. 3111-3126 (2010)

View Full Text Article

Enhanced HTML    Acrobat PDF (1670 KB)

Browse Journals / Lookup Meetings

Browse by Journal and Year


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools



Image reconstruction in diffuse optical tomography (DOT) is, in general, posed as a model-based, nonlinear optimization problem, which requires repeated use of the three-dimensional (3D) forward and inverse solvers. To cope with the computation and storage problem for some applications, such as breast tumor diagnosis, it is preferable to develop a subdomain-based parallel computation scheme. In this study, we propose a two-level image reconstruction scheme for 3D DOT, which combines the Schwarz-type domain-decomposition (DD)-based forward calculation and the matrix-decomposition (MD)-based inversion. In the forward calculation, the solution to the diffusion equation is initially obtained using a whole-domain finite difference method at a coarse grid, and then updated with a parallel DD scheme at a fine grid. The inversion procedure starts with the wavelet-decomposition-based reconstruction at a coarse grid, and then follows with a Levenberg–Marquardt least-squares solution at a fine grid, where an MD strategy is adopted for the relevant linear inversion. It is demonstrated that the combination of the DD-based forward solver and MD-based inversion allows for coarse-grain parallel implementation of both the forward and inverse issues and effectively reduces computation and storage loads for the large-scale problem. Also, both numerical simulations and phantom experiments show that MD-based linear inversion is superior to the row-fashioned algebraic reconstruction technique.

© 2010 Optical Society of America

OCIS Codes
(170.3010) Medical optics and biotechnology : Image reconstruction techniques
(170.3830) Medical optics and biotechnology : Mammography
(170.3880) Medical optics and biotechnology : Medical and biological imaging

ToC Category:
Medical Optics and Biotechnology

Original Manuscript: December 1, 2009
Revised Manuscript: April 14, 2010
Manuscript Accepted: April 23, 2010
Published: May 28, 2010

Virtual Issues
Vol. 5, Iss. 10 Virtual Journal for Biomedical Optics

Fang Yang, Feng Gao, Pingqiao Ruan, and Huijuan Zhao, "Combined domain-decomposition and matrix-decomposition scheme for large-scale diffuse optical tomography," Appl. Opt. 49, 3111-3126 (2010)

Sort:  Author  |  Year  |  Journal  |  Reset  


  1. F. Gao, Y. Xue, and H. Zhao, “Two-dimensional optical tomography of hemodynamic changes in a preterm infant brain,” Chin. Opt. Lett. 5, 472–474 (2007).
  2. G. Strangman, D. Boas, and J. Sutton, “Non-invasive neuroimaging using near-infrared light,” Biol. Psychiatry 52, 679–693 (2002). [CrossRef] [PubMed]
  3. A. Villringer and B. Chance, “Non-invasive optical spectroscopy and imaging of human brain function,” Trends Neurosci. 20, 435–442 (1997). [CrossRef] [PubMed]
  4. Z. Yuan and H. Jiang, “Image reconstruction scheme that combines modified Newton method and efficient initial guess estimation for optical tomography of finger joints,” Appl. Opt. 46, 2757–2768 (2007). [CrossRef] [PubMed]
  5. D. Grosenick, T. Moesta, and H. Wabnitz, “Time-domain optical mammography: initial clinical results on detection and characterization of breast tumor,” Appl. Opt. 42, 3170–3186(2003). [CrossRef] [PubMed]
  6. W. Mo, T. Chan, L. Chen, and N. Chen, “Quantitative characterization of optical and physiological parameters in normal breasts using time-resolved spectroscopy: in vivo results of 19 Singapore women,” J. Biomed. Opt. 14, 064004 (2009). [CrossRef]
  7. X. Intes, J. Ripoll, and Y. Chen, “In vivo continuous-wave optical breast imaging enhanced with Indocyanine Green,” Med. Phys. 30, 1039–1047 (2003). [CrossRef] [PubMed]
  8. B. Kanmani, “Diffuse optical tomography through solving a system of quadratic equations: theory and simulations,” Phys. Med. Biol. 51, 981–998 (2006). [CrossRef] [PubMed]
  9. S. Arridge, “Optical tomography in medical imaging,” Inverse Probl. 15, R41–R93 (1999). [CrossRef]
  10. H. A. Schwarz, “über einen Grenz bergang durch alternirendes Verfahren,” Ges. Math. Abhandlungen 1, 133–143(1870).
  11. M. J. Eppstein, D. E. Dougherty, D. J. Hawrysz, and E. M. Servick-Muraca, “Three-dimensional Bayesian optical image reconstruction with domain decomposition,” IEEE Trans. Med. Imaging 20, 147–163 (2001). [CrossRef] [PubMed]
  12. I.-Y. Son and X. Intes, “A 2-level domain decomposition algorithm for inverse diffuse optical tomography,” in Proceedings of the International Conference on Image Processing (IEEE, 2004), pp. 3315–3318.
  13. K. Kwon, I. Son, and B. Yazici, “Two-level domain decomposition algorithm for a nonlinear inverse DOT problem,” Proc. SPIE 5693, 459–468 (2005). [CrossRef]
  14. K. Kwon, B. Yazici, and M. Guven, “Two-level domain decomposition methods for diffuse tomography,” Inverse Probl. 22, 1533–1559 (2006). [CrossRef]
  15. X. Gu and A. Hielscher, “Parallelization of transport-theory based optical tomography algorithms by domain decomposition,” in Biomedical Optics, OSA Technical Digest (CD) (Optical Society of America, 2008), p. BSuE40.
  16. U. Trottenberg, C. Oosterlee, and A. Schüller, Multigrid(Academic, 2001).
  17. W. L. Briggs, A Multigrid Tutorial (SIAM, 1987).
  18. W. Hackbusch, Multi-Grid Methods and Applications(Springer-Verlag, 1985).
  19. F. Gao, H. Zhao, Y. Yanikawa, and Y. Yamada, “Optical tomographic mapping of cerebral haemodynamics by means of time-domain detection: methodology and phantom validation,” Phys. Med. Biol. 49, 1055–1078 (2004). [CrossRef] [PubMed]
  20. F. Gao, H. Zhao, and Y. Tanikawa, “Time-resolved diffuse optical tomography using a modified generalized spectrum technique,” IEICE Trans. Inf. & Syst. E85-D, 133–142 (2002).
  21. J. Nocedal and S. J. Wright, Numerical Optimization(Springer, 1999). [CrossRef]
  22. S. G. Mallat, “A theory for multiresolution signal decomposition: The wavelet representation,” IEEE Trans. Pattern Anal. Machine Intell. 11, 674–693 (1989). [CrossRef]
  23. W. Zhu and Y. Wang, “Total least-squares reconstruction with wavelets for optical tomography,” J. Opt. Soc. Am. A 15, 2639–2650 (1998). [CrossRef]
  24. W. Zhu and Y. Wang, “A wavelet-based multiresolution regularized least squares reconstruction approach for optical tomography,” IEEE Trans. Med. Imaging 16, 210–217(1997). [CrossRef] [PubMed]
  25. M. Bertero and P. Boccacci, Introduction to Inverse Problems in Imaging (IOP Publishing, 1998). [CrossRef]
  26. A. H. Hielscher and S. Bartel, “Parallel programming of gradient-based iterative image reconstruction schemes for optical tomography,” Comput. Methods Programs Biomed. 73, 101–113 (2004). [CrossRef] [PubMed]
  27. F. Gao, H. Zhao, Y. Tanikawa, K. Homma, and Y. Yamada, “Influences of target size and contrast on near infrared diffuse optical tomography-a comparison between featured-data and full time-resolved schemes,” Opt. Quantum Electron. 37, 1287–1304 (2005). [CrossRef]

Cited By

Alert me when this paper is cited

OSA is able to provide readers links to articles that cite this paper by participating in CrossRef's Cited-By Linking service. CrossRef includes content from more than 3000 publishers and societies. In addition to listing OSA journal articles that cite this paper, citing articles from other participating publishers will also be listed.

« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited