Monte Carlo simulation of optical coherence tomography for turbid media with arbitrary spatial distributions

We developed a Monte Carlo-based simulator of optical coherence tomography (OCT) imaging for turbid media with arbitrary spatial distributions. This simulator allows computation of both Class I diffusive reflectance due to ballistic and quasiballistic scattered photons and Class II diffusive reflectance due to multiple scattered photons. It was implemented using a tetrahedron-based mesh and importance sampling to significantly reduce computational time. Our simulation results were verified by comparing them with results from two previously validated OCT simulators for multilayered media. We present simulation results for OCT imaging of a sphere inside a background slab, which would not have been possible with earlier simulators. We also discuss three important aspects of our simulator: (1) resolution, (2) accuracy, and (3) computation time. Our simulator could be used to study important OCT phenomena and to design OCT systems with improved performance.