Related Experiment Video
Updated: Jul 15, 2025

Measuring Spatially- and Directionally-varying Light Scattering from Biological Material
Published on: May 20, 2013
Importance sampling-accelerated simulation of full-spectrum backscattered diffuse reflectance
Jianing Mao1, Yuye Ling1, Ping Xue2
1Department of Electronic Engineering, Shanghai Jiao Tong University, Shanghai 200240, China.
Abstract:
The Monte Carlo (MC) method is one of the most widely used numerical tools to model the light interaction with tissue. However, due to the low photon collection efficiency and the need to simulate the entire emission spectrum, it is computationally expensive to simulate the full-spectrum backscattered diffuse reflectance (F-BDR). Here, we propose an acceleration scheme based on importance sampling (IS). We derive the biasing sampling function tailored for simulating BDR based on the two-term scattering phase function (TT). The parameters of the TT function at different wavelengths are directly obtained by fitting the Mie scattering phase function. Subsequently, we incorporate the TT function and its corresponding biased function into the redefined IS process and realize the accelerated simulation of F-BDR. Phantom simulations based on the Fourier-domain optical coherence tomography (FD-OCT) are conducted to demonstrate the efficiency of the proposed method. Compared to the original simulator without IS, our proposed method achieves a 373× acceleration in simulating the F-BDR of the multi-layer phantom with a relative mean square error (rMSE) of less than 2%. Besides, by parallelly computing A-lines, our method enables the simulation of an entire B-scan in less than 0.4 hours. To our best knowledge, it is the first time that a volumetric OCT image of a complex phantom is simulated. We believe that the proposed acceleration method can be readily applied to fast simulations of various F-BDR-dependent applications. The source codes of this manuscript are also publicly available online.
Related Concept Videos
Attenuated Total Reflectance (ATR) Infrared Spectroscopy: Overview
The ATR process begins by directing a beam...
Aliasing
If the sampling frequency is below the Nyquist rate, these replicas overlap, preventing the original...

