Fast method for inverse determination of optical parameters from two measured signals.

Solving inverse problems requires multiple forward calculations of measured signals. We present a fast method combining graphic processing unit-accelerated Monte Carlo simulations of individual photons and a new perturbation scheme for a 300-fold speedup in comparison to conventional CPU-based approaches. The method allows rapid calculations of the diffuse reflectance and transmittance signals from a turbid sample of absorption coefficient μ(a), scattering coefficient μ(s), and anisotropy factor g based on the principle of correlated sampling. To demonstrate its strong utility, we have applied the method for determining the optical parameters of diluted intralipid samples with satisfactory results.

[1]  J. Spanier,et al.  Perturbation Monte Carlo methods to solve inverse photon migration problems in heterogeneous tissues. , 2001, Optics letters.

[2]  Angelo Sassaroli.  Fast perturbation Monte Carlo method for photon migration in heterogeneous turbid media. , 2011, Optics letters.

[3]  Stefan Andersson-Engels,et al.  Next-generation acceleration and code optimization for light transport in turbid media using GPUs , 2010, Biomedical optics express.

[4]  Huafeng Ding,et al.  A primary method for determination of optical parameters of turbid samples and application to intralipid between 550 and 1630 nm. , 2006, Optics express.

[5]  E. Gelbard,et al.  Monte Carlo Principles and Neutron Transport Problems , 2008 .

[6]  Jun Q. Lu,et al.  Optical properties of porcine skin dermis between 900 nm and 1500 nm , 2001, Physics in medicine and biology.

[7]  X. H. Hu,et al.  Monte carlo simulation of converging laser beams propagating in biological materials. , 1999, Applied optics.

[8]  Xin-Hua Hu,et al.  Effect of surface roughness on determination of bulk tissue optical parameters. , 2003, Optics letters.

[9]  Xin-Hua Hu,et al.  Numerical study of reflectance imaging using a parallel Monte Carlo method. , 2007, Medical physics.

[10]  H. Rief.  Generalized Monte Carlo perturbation algorithms for correlated sampling and a second-order Taylor series approach , 1984 .

[11]  A. Welch,et al.  Determining the optical properties of turbid mediaby using the adding-doubling method. , 1993, Applied optics.

[12]  L Wang,et al.  MCML--Monte Carlo modeling of light transport in multi-layered tissues. , 1995, Computer methods and programs in biomedicine.

[13]  Jun Q. Lu,et al.  Determination of complex refractive index of polystyrene microspheres from 370 to 1610 nm , 2003, Physics in medicine and biology.

[14]  Huafeng Ding,et al.  Bulk optical parameters of porcine skin dermis at eight wavelengths from 325 to 1557 nm. , 2005, Optics letters.