Abstract:Significance: Perturbation and differential Monte Carlo (pMC/dMC) methods, used in conjunction with nonlinear optimization methods, have been successfully applied to solve inverse problems in diffuse optics. Application of pMC to systems over a large range of optical properties requires optimal "placement" of baseline conventional Monte Carlo (cMC) simulations to minimize the pMC variance. The inability to predict the growth in pMC solution uncertainty with perturbation size limits the application of pMC, especially for multispectral datasets where the variation of optical properties can be substantial. Aim: We aim to predict the variation of pMC variance with perturbation size without explicit computation of perturbed photon weights. Our proposed method can be used to determine the range of optical properties over which pMC predictions provide sufficient accuracy. This method can be used to specify the optical properties for the reference cMC simulations that pMC utilizes to provide accurate predictions over a desired optical property range. Approach: We utilize a conventional error propagation methodology to calculate changes in pMC relative error for Monte Carlo simulations. We demonstrate this methodology for spatially resolved diffuse reflectance measurements with ±20% scattering perturbations. We examine the performance of our method for reference simulations spanning a broad range of optical properties relevant for diffuse optical imaging of biological tissues. Our predictions are computed using the variance, covariance, and skewness of the photon weight, path length, and collision distributions generated by the reference simulation. Results: We find that our methodology performs best when used in conjunction with reference cMC simulations that utilize Russian Roulette (RR) method. Specifically, we demonstrate that for a proximal detector placed immediately adjacent to the source, we can estimate the pMC relative error within 5% of the true value for scattering perturbations in the range of [-15%,+20%]. For a distal detector placed at ∼3 transport mean free paths relative to the source, our method provides relative error estimates within 20% for scattering perturbations in the range of [-8%,+15%]. Moreover, reference simulations performed at lower (μs'/μa) values showed better performance for both proximal and distal detectors. Conclusions: These findings indicate that reference simulations utilizing continuous absorption weighting (CAW) with the Russian Roulette method and executed using optical properties with a low (μs'/μa) ratio spanning the desired range of μs values, are highly advantageous for the deployment of pMC to obtain radiative transport estimates over a wide range of optical properties.

Application of Monte Carlo Simulation-Based Photon Migration for Enhanced Understanding of Near-Infrared (nir) Diffuse Reflectance. Part Ii: Photon Radial Diffusion in Nir Chemical Images

Application of Monte Carlo Simulation-Based Photon Migration for Enhanced Understanding of Near-Infrared (nir) Diffuse Reflectance. Part I: Depth of Penetration in Pharmaceutical Materials

A Sequential Method for Measuring the Optical Properties of Two-Layer Media with Spatially-Resolved Diffuse Reflectance: Simulation Study

Pharmaceutical applications of separation of absorption and scattering in near-infrared spectroscopy (NIRS).

Numerical study of reflectance imaging using a parallel Monte Carlo method.

A Pervasive Approach for Determining the Optimal Tissue Windows for Near-Infrared Fluorescence Imaging

2-D Image Localization in Hyperspectral Image Analysis of Pharmaceutical Materials

Analysis of relative error in perturbation Monte Carlo simulations of radiative transport

Near-Infrared-II Bioimaging for in Vivo Quantitative Analysis

Monte Carlo-based data generation for efficient deep learning reconstruction of macroscopic diffuse optical tomography and topography applications

Estimation of Photon Distribution within Biological Tissue Using Monte Carlo Simulation

Efficient inversion strategies for estimating optical properties with Monte Carlo radiative transport models

Near-Infrared-II (NIR-II) Bioimaging via Off-Peak NIR-I Fluorescence Emission

Near-infrared-II Deep Tissue Fluorescence Microscopy and Application

Direct reconstruction of pharmacokinetic-rate images of optical fluorophores from NIR measurements

A NIR Optical Imaging System with Diffuse Photons

Drug Quantification in Turbid Media by Fluorescence Imaging Combined with Light-Absorption Correction Using White Monte Carlo Simulations.

Exploring NIR Spectroscopy Data: A Practical Chemometric Tutorial for Analyzing Freeze-Dried Pharmaceutical Formulations

Photon Migration in Turbid Media: a Finite-Element Solution for the Most Favorable Path

Quantitative interpretations of Visible-NIR reflectance spectra of blood

Finding The Most Probable Diffuse Path Of Photons In Tissue By Monte Carlo Method