Algorithm for rapid determination of optical scattering parameters

Abstract

Preliminary experiments at the NIST Spectral Tri-function Automated Reference Reflectometer (STARR) facility have been conducted with the goal of providing the diffuse optical properties of a solid reference standard with optical properties similar to human skin. Here, we describe an algorithm for determining the best-fit parameters and the statistical uncertainty associated with the measurement. The objective function is determined from the profile log likelihood, including both experimental and Monte Carlo uncertainties. Initially, the log likelihood is determined over a large parameter search box using a relatively small number of Monte Carlo samples such as 2·104. The search area is iteratively reduced to include the 99.9999% confidence region, while doubling the number of samples at each iteration until the experimental uncertainty dominates over the Monte Carlo uncertainty. Typically this occurs by 1.28·106 samples. The log likelihood is then fit to determine a 95% confidence ellipse. The inverse problem requires the values of the log likelihood on many points. Our implementation uses importance sampling to calculate these points on a grid in an efficient manner. Ultimately, the time-to-solution is approximately six times the cost of a Monte Carlo simulation of the radiation transport problem for a single set of parameters with the largest number of photons required. The results are found to be 64 times faster than our implementation of Particle Swarm Optimization.

General Method for Determining Light Scattering and Absorption of Nanoparticle Composites

“…The two‐background method depends on the accurate characterization of a white and a black reference material. Other approaches to the inverse problem are computationally complex and involve Monte‐Carlo simulations and numerical solutions of the general radiative transfer equation . A simpler alternative is offered by inversion schemes using the two‐flux theory, but it has been found that the T and R values obtained from the solutions for S and K sometimes do not converge to the experimental data.…”

MCMLpar and MCSLinv: A Parallel Version of MCML and an Inverse Monte Carlo Algorithm to Calculate Optical Scattering Parameters

“…Our recent research on Monte Carlo simulation of light propagation in biomedical phantoms has involved the implementation of an algorithm to solve for the optical attenuation and scattering coefficients of a single-layer material, or µa and µs [1]. It also involved the development of a parallel version of MCML [2], a popular light propagation model with over 2100 citations, according to ScienceDirect.…”

Correction of an adding-doubling inversion algorithm for the measurement of the optical parameters of turbid media

“…The index of refraction of PDMS was estimated from the mean value of the Fresnel reflectance of the three samples (uncertainty 0.015 at k = 1). A set of three specular phantoms from INO was measured to estimate the index of refraction of the polyurethane material [26]. We assume that the B0430 samples are made of an identical base material (uncertainty 0.015 at k = 1).…”