Importance sampling-accelerated simulation of full-spectrum backscattered diffuse reflectance.

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.