Fast and customizable image formation model for optical coherence tomography.

Optical coherence tomography (OCT) is a technique that performs high-resolution, three-dimensional, imaging of semi-transparent scattering biological tissues. Models of OCT image formation are needed for applications such as aiding image interpretation and validating OCT signal processing techniques. Existing image formation models generally trade off between model realism and computation time.

On the inverse problem in optical coherence tomography.

We examine the inverse problem of retrieving sample refractive index information in the context of optical coherence tomography. Using two separate approaches, we discuss the limitations of the inverse problem which lead to it being ill-posed, primarily as a consequence of the limited viewing angles available in the reflection geometry. This is first considered from the theoretical point of view of diffraction tomography under a weak scattering approximation.

Application of a Taylor series approximation to the Debye-Wolf integral in time-domain numerical electromagnetic simulations.

Finite-difference time-domain (FDTD) and pseudospectral time-domain (PSTD) methods are numerical electromagnetic simulation techniques that have been employed to perform rigorous simulations of broadband illuminations in several contexts. However, the computational cost of calculating the incident source fields introduced into the FDTD/PSTD grid can be considerable. In some cases, this can exceed the computational cost of what might be considered the principal part of the FDTD/PSTD algorithm, which calculates the spatial derivative of fields throughout the computational grid.

Controlling the optical pathlength in continuous-wave reflectance spectroscopy using polarization.

We investigate potential improvements of continuous-wave diffuse reflectance spectroscopy within highly scattering media by employing polarization gating. Simulations are used to show the extent at which the effective optical pathlength varies in a typical scattering medium as a function of the optical wavelength, the total level of absorption, and the selected polarization channels, including elliptical and circular polarization channels. Experiments then demonstrate that a wavelength dependent polarization gating scheme may reduce the prior knowledge required to solve the problem of chromophore quantification.

Approximate image synthesis in optical coherence tomography.

Full-wave models of OCT image formation, which are based on Maxwell's equations, are highly realistic. However, such models incur a high computational cost, particularly when modelling sample volumes consistent with those encountered in practice. Here, we present an approximate means of synthesizing volumetric image formation to reduce this computational burden.

Efficient inversion strategies for estimating optical properties with Monte Carlo radiative transport models.

Significance: Indirect imaging problems in biomedical optics generally require repeated evaluation of forward models of radiative transport, for which Monte Carlo is accurate yet computationally costly. We develop an approach to reduce this bottleneck, which has significant implications for quantitative tomographic imaging in a variety of medical and industrial applications.

Aim: Our aim is to enable computationally efficient image reconstruction in (hybrid) diffuse optical modalities using stochastic forward models.

Designing phantoms to accurately replicate circular depolarization in biological scattering media.

We introduce an iterative method for designing optical phantoms that are able to replicate the depolarization profiles of various target media, including colloidal suspensions of Intralipid, bovine milk, and ex vivo samples of ovine kidney cortex tissue. The designed phantoms comprise spherical scattering particles with fine-tuned size distributions and are capable of simultaneously reproducing spatially resolved intensity measurements and depolarization measurements of target media when illuminated with circularly polarized light.

Characterizing the depolarization of circularly polarized light in turbid scattering media.

We investigate the effectiveness of various bulk optical parameters in characterizing the degree of circular polarization (DOCP) of light diffusely reflected from scattering media. It is demonstrated that the traditional set of bulk optical parameters (namely, the scattering and absorption coefficients and the scattering asymmetry parameter) fail to characterize the observed depolarization. However, we find that there exists an additional parameter connected to the circular polarization memory phenomenon that consistently relates to observations, even in media with widely varying refractive indices and particle size distributions.

Diffuse correlation tomography in the transport regime: A theoretical study of the sensitivity to Brownian motion.

Diffuse correlation tomography (DCT) uses the electric-field temporal autocorrelation function to measure the mean-square displacement of light-scattering particles in a turbid medium over a given exposure time. The movement of blood particles is here estimated through a Brownian-motion-like model in contrast to ordered motion as in blood flow. The sensitivity kernel relating the measurable field correlation function to the mean-square displacement of the particles can be derived by applying a perturbative analysis to the correlation transport equation (CTE).

Optimized diffusion approximation.

We show that the diffusion approximation (DA) to the radiative transport equation, which is commonly used in biomedical optics to describe propagation of light in tissues, contains a previously unexplored adjustable parameter. This parameter is related to the rate of exponential decay of the reduced intensity. In conventional theories, there are two distinct choices for this parameter.

Numerical investigation of polarization filtering for direct optical imaging within scattering media.

We investigate the ability of polarization filtering to improve direct imaging of absorbing objects which are buried within scattering environments. We extend on previous empirical investigations by exploiting an efficient perturbation-based formalism, which is applicable to arbitrarily arranged sources and detectors with arbitrary polarizations. From this approach, we are able in some cases to find certain non-trivial linear combinations of polarization measurement channels that maximize the object resolution and visibility.

Reciprocity relation for the vector radiative transport equation and its application to diffuse optical tomography with polarized light.

We derive a reciprocity relation for the 3D vector radiative transport equation that describes propagation of polarized light in multiple-scattering media. We then show how this result, together with translational invariance of a plane-parallel sample, can be used to efficiently compute the sensitivity kernel of diffuse optical tomography by Monte Carlo simulations. Numerical examples of polarization-selective sensitivity kernels are given.

Application of circularly polarized light for non-invasive diagnosis of cancerous tissues and turbid tissue-like scattering media.

Polarization-based optical techniques have become increasingly popular in the field of biomedical diagnosis. In the current report we exploit the directional awareness of circularly and/or elliptically polarized light backscattered from turbid tissue-like scattering media. We apply circularly and elliptically polarized laser light which illuminates the samples of interest, and a standard optical polarimeter is used to observe the polarization state of light backscattered a few millimeters away from the point of incidence.

Propagation of coherent polarized light in turbid highly scattering medium.

Within the framework of further development of unified Monte Carlo code for the needs of biomedical optics and biophotonics, we present an approach for modeling of coherent polarized light propagation in highly scattering turbid media, such as biological tissues. The temporal coherence of light, linear and circular polarization, interference, and the helicity flip of circularly polarized light due to reflection at the medium boundary and/or backscattering events are taken into account. To achieve higher accuracy in the results and to speed up the modeling, the implementation of the code utilizes parallel computing on NVIDIA graphics processing units using Compute Unified Device Architecture.