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.

A Stepwise Method for Estimating Optical Properties of Two-Layer Turbid Media from Spatial-Frequency Domain Reflectance

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

Measuring the Optical Properties of Two-Layer Agricultural Tissue with a Sequential Method from Spatially-Resolved Diffuse Reflectance, Part I: Analysis and Simplification

Improved Algorithm for Estimating the Optical Properties of Food Products Using Spatially-Resolved Diffuse Reflectance

A Two-Step Parameter Optimization Algorithm for Improving Estimation of Optical Properties Using Spatial Frequency Domain Imaging

Optimization of Spatial Frequency Domain Imaging Technique for Estimating Optical Properties of Food and Biological Materials

Finite element modeling of light propagation in turbid media under illumination of a continuous-wave beam

Reconstruction of Optical Properties in Turbid Media: Omitting the Need of the Collimated Transmission for an Integrating Sphere Setup

Spatial-frequency domain imaging coupled with frequency optimization for estimating optical properties of two-layered food and agricultural products

Effects of Optical Variables in a Single Integrating Sphere System on Estimation of Scattering Properties of Turbid Media

An Integrated Fiber-Optic Probe Combined with Support Vector Regression for Fast Estimation of Optical Properties of Turbid Media

Finite Element Simulation of Light Transfer in Turbid Media under Structured Illumination.

Estimation of optical properties of turbid media using spatially resolved diffuse reflectance combined with LSTM-attention network

Multi-parameter Estimation in Semitransparent Graded-Index Media Based on Coupled Optical and Thermal Information

Methods Comparision of Estimating Optical Turbulence Profile Using Conventional Meteorology Parameters

Tutorial on methods for estimation of optical absorption and scattering properties of tissue

A Fast and Accurate Method for the Simulation of the Diffusing Temporal Light Correlation in Multi-Layered Turbid Media

Quantitative estimation of optical properties in bilayer media within the subdiffusive regime using tilted fiber-optic probe diffuse reflectance spectroscopy, part 2: probe design, realization, and experimental validation

Numerical Study of Turbid Slab Optical Properties Reconstruction from Multiple Scattering Signals Using Time-Based Markov Chain Model

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

Estimation of optical constants from multiple-scattered light using approximations for single particle scattering characteristics.