OSA's Digital Library

Biomedical Optics Express

Biomedical Optics Express

  • Editor: Joseph A. Izatt
  • Vol. 1, Iss. 1 — Aug. 2, 2010
  • pp: 165–175

Mesh-based Monte Carlo method using fast ray-tracing in Plücker coordinates

Qianqian Fang  »View Author Affiliations

Biomedical Optics Express, Vol. 1, Issue 1, pp. 165-175 (2010)

View Full Text Article

Enhanced HTML    Acrobat PDF (3747 KB)

Browse Journals / Lookup Meetings

Browse by Journal and Year


Lookup Conference Papers

Close Browse Journals / Lookup Meetings

Article Tools



We describe a fast mesh-based Monte Carlo (MC) photon migration algorithm for static and time-resolved imaging in 3D complex media. Compared with previous works using voxel-based media discretization, a mesh-based approach can be more accurate in modeling targets with curved boundaries or locally refined structures. We implement an efficient ray-tracing technique using Plücker Coordinates. The Barycentric coordinates computed from Plücker-formed ray-tracing enables us to use linear Lagrange basis functions to model both media properties and fluence distribution, leading to further improvement in accuracy. The Plücker-coordinate ray-polygon intersection test can be extended to hexahedral or high-order elements. Excellent agreement is found when comparing mesh-based MC with the analytical diffusion model and 3D voxel-based MC code in both homogeneous and heterogeneous cases. Realistic time-resolved imaging results are observed for a complex human brain anatomy using mesh-based MC. We also include multi-threading support in the software and will port it to a graphics processing unit platform in the near future.

© 2010 OSA

OCIS Codes
(170.3660) Medical optics and biotechnology : Light propagation in tissues
(170.5280) Medical optics and biotechnology : Photon migration
(170.7050) Medical optics and biotechnology : Turbid media

ToC Category:
Optics of Tissue and Turbid Media

Original Manuscript: June 8, 2010
Revised Manuscript: July 11, 2010
Manuscript Accepted: July 13, 2010
Published: July 15, 2010

Qianqian Fang, "Mesh-based Monte Carlo method using fast ray-tracing in Plücker coordinates," Biomed. Opt. Express 1, 165-175 (2010)

Sort:  Author  |  Year  |  Journal  |  Reset  


  1. B. C. Wilson and G. Adam, “A Monte Carlo model for the absorption and flux distributions of light in tissue,” Med. Phys. 10(6), 824–830 (1983). [CrossRef] [PubMed]
  2. P. van der Zee, and D. T. Delpy, “Simulation of the point spread function for light in tissue by a Monte Carlo method,” in Oxygen Transport to Tissue I. A. Silver and A. Silver, eds. (Plenum, New York,), Vol. IX. (1987)
  3. Y. Hasegawa, Y. Yamada, M. Tamura, and Y. Nomura, “Monte Carlo simulation of light transmission through living tissues,” Appl. Opt. 30(31), 4515–4520 (1991). [CrossRef]
  4. M. Hiraoka, M. Firbank, M. Essenpreis, M. Cope, S. R. Arridge, P. van der Zee, and D. T. Delpy, “A Monte Carlo investigation of optical pathlength in inhomogeneous tissue and its application to near-infrared spectroscopy,” Phys. Med. Biol. 38(12), 1859–1876 (1993). [CrossRef] [PubMed]
  5. L. H. Wang, S. L. Jacques, and L. Q. Zheng, “MCML - Monte Carlo modeling of light transport in multilayered tissues,” Comput. Meth. Prog. Bio. 47(2), 131–146 (1995). [CrossRef]
  6. J. Ripoll, M. Nieto-Vesperinas, S. R. Arridge, and H. Dehghani, “Boundary conditions for light propagation in diffusive media with nonscattering regions,” J. Opt. Soc. Am. A 17(9), 1671–1681 (2000). [CrossRef] [PubMed]
  7. S. R. Arridge, H. Dehghani, M. Schweiger, and E. Okada, “The finite element model for the propagation of light in scattering media: a direct method for domains with nonscattering regions,” Med. Phys. 27(1), 252–264 (2000). [CrossRef] [PubMed]
  8. F. Martelli, D. Contini, A. Taddeucci, and G. Zaccanti, “Photon migration through a turbid slab described by a model based on diffusion approximation. II. Comparison with Monte Carlo results,” Appl. Opt. 36(19), 4600–4612 (1997). [CrossRef] [PubMed]
  9. Y. Xu, Q. Zhang, and H. Jiang, “Optical image reconstruction of non-scattering and low scattering heterogeneities in turbid media based on the diffusion approximation model,” J. Opt. A, Pure Appl. Opt. 6(1), 29–35 (2004). [CrossRef]
  10. A. H. Hielscher, R. E. Alcouffe, and R. L. Barbour, “Comparison of finite-difference transport and diffusion calculations for photon migration in homogeneous and heterogeneous tissues,” Phys. Med. Biol. 43(5), 1285–1302 (1998). [CrossRef] [PubMed]
  11. A. Custo, W. M. Wells, A. H. Barnett, E. M. Hillman, and D. A. Boas, “Effective scattering coefficient of the cerebral spinal fluid in adult head models for diffuse optical imaging,” Appl. Opt. 45(19), 4747–4755 (2006). [CrossRef] [PubMed]
  12. G. Ma, J. F. Delorme, P. Gallant, and D. A. Boas, “Comparison of simplified Monte Carlo simulation and diffusion approximation for the fluorescence signal from phantoms with typical mouse tissue optical properties,” Appl. Opt. 46(10), 1686–1692 (2007). [CrossRef] [PubMed]
  13. D. R. Kirkby and D. T. Delpy, “Parallel operation of Monte Carlo simulations on a diverse network of computers,” Phys. Med. Biol. 42(6), 1203–1208 (1997). [CrossRef] [PubMed]
  14. A. Colasanti, G. Guida, A. Kisslinger, R. Liuzzi, M. Quarto, P. Riccio, G. Roberti, and F. Villani, “Multiple processor version of a Monte Carlo code for photon transport in turbid media,” Comput. Phys. Commun. 132(1–2), 84–93 (2000). [CrossRef]
  15. S. R. Arridge, M. Schweiger, M. Hiraoka, and D. T. Delpy, “A finite element approach for modeling photon transport in tissue,” Med. Phys. 20(2 Pt 1), 299–309 (1993). [CrossRef] [PubMed]
  16. K. Ren, G. S. Abdoulaev, G. Bal, and A. H. Hielscher, “Algorithm for solving the equation of radiative transfer in the frequency domain,” Opt. Lett. 29(6), 578–580 (2004). [CrossRef] [PubMed]
  17. CGAL, Computational Geometry Algorithms Library, URL: http://www.cgal.org
  18. Q. Fang and D. Boas, “Tetrahedral mesh generation from volumetric binary and gray-scale images,” Proceedings of IEEE International Symposium on Biomedical Imaging 2009, 1142–1145 (2009).
  19. T. Binzoni, T. S. Leung, R. Giust, D. Rüfenacht, and A. H. Gandjbakhche, “Light transport in tissue by 3D Monte Carlo: influence of boundary voxelization,” Comput. Methods Programs Biomed. 89(1), 14–23 (2008). [CrossRef] [PubMed]
  20. D. A. Boas, J. P. Culver, J. J. Stott, and A. K. Dunn, “Three dimensional Monte Carlo code for photon migration through complex heterogeneous media including the adult human head,” Opt. Express 10(3), 159–170 (2002). [PubMed]
  21. E. Alerstam, T. Svensson, and S. Andersson-Engels, “Parallel computing with graphics processing units for high-speed Monte Carlo simulation of photon migration,” J. Biomed. Opt. 13(6), 060504 (2008). [CrossRef] [PubMed]
  22. Q. Fang and D. A. Boas, “Monte Carlo simulation of photon migration in 3D turbid media accelerated by graphics processing units,” Opt. Express 17(22), 20178–20190 (2009). [CrossRef] [PubMed]
  23. Y. Fukui, Y. Ajichi, and E. Okada, “Monte Carlo prediction of near-infrared light propagation in realistic adult and neonatal head models,” Appl. Opt. 42(16), 2881–2887 (2003). [CrossRef] [PubMed]
  24. E. Margallo-Balbás and P. J. French, “Shape based Monte Carlo code for light transport in complex heterogeneous Tissues,” Opt. Express 15(21), 14086–14098 (2007). [CrossRef] [PubMed]
  25. N. Ren, J. Liang, X. Qu, J. Li, B. Lu, and J. Tian, “GPU-based Monte Carlo simulation for light propagation in complex heterogeneous tissues,” Opt. Express 18(7), 6811–6823 (2010). [CrossRef] [PubMed]
  26. N. Platis and T. Theoharis, “Fast ray-tetrahedron intersection using Plücker coordinates,” Journal of Graphics Tools 8(4), 37–48 (2003).
  27. H.S.M. Coxeter, Non-Euclidean geometry, Univ. Toronto Press pp. 88–90 (1965)
  28. Q. Fang, “Mesh-based Monte Carlo – the software”, URL: http://mcx.sourceforge.net/mmc/
  29. H. Shen and G. Wang, “A tetrahedron-based inhomogeneous Monte Carlo optical simulator,” Phys. Med. Biol. 55(4), 947–962 (2010). [CrossRef] [PubMed]
  30. M. Meyer, H. Lee, A. Barr, and M. Desbrun, “Generalized Barycentric Coordinates on Irregular Polygons,” Journal of Graphics Tools 7(1), 13–22 (2002).
  31. C. A. Felippa, Advanced Finite Element Methods, URL: http://www.colorado.edu/engineering/cas/courses.d/AFEM.d/
  32. OpenMP Architecture Review Board, OpenMP Application Program Interface, Version 3.0, URL: http://www.openmp.org/mp-documents/spec30.pdf (2008)
  33. NVIDIA, CUDA Compute Unified Device Architecture - Programming Guide, Version 3.0 (2010)
  34. Khronos Group, OpenCL 1.0 Specification, URL: http://www.khronos.org/opencl/ (2008)
  35. H. Si, and K. Gaertner, “Meshing piecewise linear complexes by constrained Delaunay tetrahedralizations,” In Proc. of the 14th Int. Meshing Roundtable, pp. 147–163, (2005)
  36. P. Fleischmann, B. Haindl, R. Kosik, and S. Selberherr, “Investigation of a mesh criterion for three-dimensional finite element diffusion simulation,” Proc. Int. Conf. Simulation Semiconductor Processes Devices (SISPAD '99). 71–74 (1999)
  37. D. A. Boas, M. A. O’Leary, B. Chance, and A. G. Yodh, “Scattering of diffuse photon density waves by spherical inhomogeneities within turbid media: analytic solution and applications,” Proc. Natl. Acad. Sci. U.S.A. 91(11), 4887–4891 (1994). [CrossRef] [PubMed]
  38. D. L. Collins, A. P. Zijdenbos, V. Kollokian, J. G. Sled, N. J. Kabani, C. J. Holmes, and A. C. Evans, “Design and construction of a realistic digital brain phantom,” IEEE Trans. Med. Imaging 17(3), 463–468 (1998). [CrossRef] [PubMed]
  39. J. Vollmer, R. Mencel, and H. Müller, ““Improved Laplacian smoothing of noisy surface meshes,” In Proc,” EuroGraphics 99, 131–138 (1999).
  40. A. N. Yaroslavsky, P. C. Schulze, I. V. Yaroslavsky, R. Schober, F. Ulrich, and H.-J. Schwarzmaier, “Optical properties of selected native and coagulated human brain tissues in vitro in the visible and near infrared spectral range,” Phys. Med. Biol. 47(12), 2059–2073 (2002). [CrossRef] [PubMed]

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.


Fig. 1 Fig. 2 Fig. 3
Fig. 4

Supplementary Material

» Media 1: AVI (1689 KB)     

« Previous Article  |  Next Article »

OSA is a member of CrossRef.

CrossCheck Deposited