English (United Kingdom)

https://curated-unify.zendy.io/wp-json/zendy-region/v1/featured_content/oa?rat=en

https://curated-unify.zendy.io/wp-json/zendy-region/v1/highlighted_journal/

Global: Zendy Plus annual

Presents the access of premium content as premium feature

Premium Content

Presents the keyphrase highlighting as premium feature

Keyphrase Highlighting

Presents the summarisation as premium feature

Summarisation

Insights

Presents the pdf analysis as premium feature

PDF Analysis

Presents the zaia usage as premium feature

ZAIA

Global: Zendy Tools annual

Global: Zendy Free Plan

Global: Zendy Tools monthly

Global: Zendy Plus monthly

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. The algorithm for the inverse problem involves three components: a profile log-likelihood evaluation, importance sampling, and an optimization procedure. In this article, we present two C++/OpenMP codes: MCMLpar, the parallel version of MCML, and MCSLinv, a program to solve the inverse problem for a single layer. In MCMLpar, [3] a particle class is introduced to represent photons. The particles all enter an optical medium at the same point, travelling parallel to the vector normal to the surface. The particle weights are initially set to the probability of transmitting into the material, calculated by the Fresnel equation using the air and first layer indices of refraction. Each particle’s path length to an interaction is calculated by sampling the exponential distribution with mean μt = μa + μs. If the interaction is within the layer, the particle’s direction is updated by sampling the Henyey-Greenstein distribution for the polar angle θ relative to the direction of travel. The relative azimuthal angle φ is sampled uniformly on [0, 2π). The particle’s

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