Accurate and efficient solutions of the three dimensional radiative transport equation were derived in all domains for the case of layered scattering media. Index mismatched boundary conditions based on Fresnel’s equations were implemented. Arbitrary rotationally symmetric phase functions can be applied to characterize the scattering in the turbid media. Solutions were derived for an obliquely incident beam having arbitrary spatial profiles. The derived solutions were successfully validated with Monte Carlo simulations and partly compared with analytical solutions of the diffusion equation.

An algorithm for the simulation of image formation in Fourier domain optical coherence tomography (OCT) for an infinitely long cylinder is presented. The analytical solution of Maxwell’s equations for light scattering by a single cylinder is employed for the case of perpendicular incidence to calculate OCT images. The A-scans and the time-resolved scattered intensities are compared to geometrical optics results calculated with a ray tracing approach. The reflection peaks, including the whispering gallery modes, are identified. Additionally, the Debye series expansion is employed to identify single peaks in the OCT A-scans. Furthermore, a Gaussian beam is implemented in order to simulate lateral scanning over the cylinder for two-dimensional B-scans. The fields are integrated over a certain angular range to simulate a detection aperture. In addition, the solution for light scattering by layered cylinders is employed and the various layers are identified in the resulting OCT image. Overall, the simulations in this work show that OCT images do not always display the real surface of investigated samples.

A method to correct for surface scattering in spatial frequency domain imaging (SFDI) is presented. The use of a modified analytical solution of the radiative transfer equation allows calculation of the reflectance and the phase of a rough semi-infinite geometry so that both spatial frequency domain reflectance and phase can be applied for precise retrieval of the bulk optical properties and the surface scattering. For validation of the method, phantoms with different surface roughness were produced. Contrarily, with the modified theory, it was possible to dramatically reduce systematic errors due to surface scattering. The evaluation of these measurements with the state-of-the-art theory and measuring modality, i.e., using crossed linear polarizers, reveals large errors in the determined optical properties, depending on the surface roughness, of up to ≈100  %  . These results were confirmed with SFDI measurements on a phantom that has a structured rough surface.

No abstract

No abstract

The propagation of different focused beams (e.g., Gaussian or quasi-Bessel beams) through scattering media is studied. The finite-difference time-domain method, a numerical solution of Maxwell's equations, is applied to propagate the light beams in two dimensions. The focused beams are modeled by applying the angular spectrum of the plane waves method. The results show that weakly focused beams exhibit comparable performance to strongly focused beams in delivering focused light deep into scattering media.

No abstract

We report a theoretical study on the determination of three optical properties from spatially resolved reflectance calculations. In particular, the reduced scattering coefficient μs′, the absorption coefficient μa, and the recently defined phase function parameter σ are identified. The solution of the inverse problem is based on the principal component analysis of a large set of reflectance profiles that were calculated using an analytical solution of the radiative transfer equation. Different phase function types were studied to test the method in the range of 0.63  mm−1≤μs′≤4.2  mm−1 and 0.002  mm−1≤μa≤0.1  mm−1. For curves impaired with noise, we were able to reconstruct μs′ and μa with relative median errors of 2.5% and 12%, respectively, and σ with an absolute median error of 0.04.