Effects of the surface boundary on the determination of the optical properties of a turbid medium with time-resolved reflectance

Fast method of retrieving the asymmetry factor and scattering albedo from the maximum time-resolved reflectance of participating media

This research presents a parametric study of the time-resolved hemispherical reflectance of a semi-infinite plane-parallel slab of homogeneous, nonemitting, absorbing, and anisotropic scattering medium exposed to a collimated Gaussian pulse. The one-dimensional transient radiative transfer equation was solved by using the finite volume method. The internal reflection at the medium-air interface caused by the mismatch of the refractive indices was considered. In particular, this work focused on the maximum diffuse hemispherical reflectance. Three different optical regions were identified according to the dimensionless pulsewidth βct_p. The correlation between the normalized maximum hemispherical reflectance and βct_p was conformed to the Boltzmann function. The coefficients in the correlating functions of the match and mismatch refractive index cases were fitted as polynomial fitting functions of the single scattering albedo ω and Henyey-Greenstein asymmetric factor g. Thus, ω and g can be simultaneously reconstructed by the semi-empirical correlations without solving the forward model. In conclusion, the proposed method can potentially retrieve the asymmetry factor and single scattering albedo of participating media accurately and efficiently.

Fluorescence time-resolved imaging system embedded in an ultrasound prostate probe

Ultrasound imaging (US) of the prostate has a low specificity to distinguish tumors from the surrounding tissues. This limitation leads to systematic biopsies. Fluorescent diffuse optical imaging may represent an innovative approach to guide biopsies to tumors marked with high specificity contrast agents and therefore enable an early detection of prostate cancer. This article describes a time-resolved optical system embedded in a transrectal US probe, as well as the fluorescence reconstruction method and its performance. Optical measurements were performed using a pulsed laser, optical fibers and a time-resolved detection system. A novel fast reconstruction method was derived and used to locate a 45 µL ICG fluorescent inclusion at a concentration of 10 µM, in a liquid prostate phantom. Very high location accuracy (0.15 cm) was achieved after reconstruction, for different positions of the inclusion, in the three directions of space. The repeatability, tested with ten sequential measurements, was of the same order of magnitude. Influence of the input parameters (optical properties and lifetime) is presented. These results confirm the feasibility of using optical imaging for prostate guided biopsies.

Real time optical method for localization of inclusions embedded in turbid media

A simple and fast time-domain method for localizing inclusions, fluorescent optical probes or absorbers, is presented. The method offers new possibilities for situations where complete tomographic measurements are not permitted by the examined object, for example in endoscopic examination of the human prostate or the oesophagus. Feasibility has been envisioned with a phantom study conducted on a point-like fluorochrome embedded in a diffusing medium mimicking the optical properties of biological tissues.

Chi2 analysis for estimating the accuracy of optical properties derived from time resolved diffuse-reflectance

Weighted residuals and the reduced chi(2) (chi(R) (2)) value are investigated with regard to their relevance for assessing optical property estimates using the diffusion equation for time-resolved measurements in turbid media. It is shown and explained, for all photon counting experiments including lifetime estimation, why chi(R) (2) increases linearly with the number of photons when there is a model bias. Only when a sufficient number of photons has been acquired, chi(R) (2) is a pertinent value for assessing the accuracy of mu(a) and mu(s)' estimates. It was concluded that chi(R) (2) is of particular interest for cases of small interfiber separation, low-level scattering, strong absorption and incorrect measurement of instrument response function. It was also found that chi(R) (2) is less pertinent for judging mu(a) in case of air boundary effects.

Double-layer estimation of intra- and extracerebral hemoglobin concentration with a time-resolved system

We present in vivo measurements of baseline physiology from five subjects with a four-wavelength (690, 750, 800, and 850 nm) time-resolved optical system. The measurements were taken at four distances: 10, 15, 25, and 30 mm. All distances were fit simultaneously with a two-layered analytical model for the absorption and reduced scattering coefficient of both layers. The thickness of the first layer, comprising the skin, scalp, and cerebrospinal fluid, was obtained from anatomical magnetic resonance images. The fitting procedure was first tested with simulations before being applied to in vivo measurements and verified that this procedure permits accurate characterization of the hemoglobin concentrations in the extra- and intracerebral tissues. Baseline oxyhemoglobin, deoxyhemoglobin, and total hemoglobin concentrations and oxygen saturation were recovered from in vivo measurements and compared to the literature. We observed a noticeable intersubject variability of the hemoglobin concentrations, but constant values for the cerebral hemoglobin oxygen saturation.

Approximations of the measurable quantity in diffuse optical problems: theoretical analysis of model deviations

Photon density and photon flux are widely used to model the measurable quantity in diffuse optical tomography problems. However, it is not these two quantities that are actually measured, but rather the radiance accepted by the detection system. We provide a theoretical analysis of the model deviations related to the choice of the measurable quantity-either photon density or flux. By using the diffusion approximation to the radiative transfer equation and its solution with extrapolated boundary conditions, an exact analytical expression of the measurable quantity has been obtained. This expression has been employed as a reference to assess model deviation when considering the photon density or the photon flux as the measurable quantity. For the case of semi-infinite geometry and for both continuous wave and time domains, we show that the photon density approximates the measurable quantity better than the photon flux. We also demonstrate that the validity of this approximation strongly depends on the optical parameters.

Estimation of optical properties of turbid media: experimental comparison of spatially and temporally resolved reflectance methods

We compare two methods for the optical characterization of turbid media. The estimates of the absorption and reduced scattering coefficients (mu(a) and mu(')(s)) by a spatially resolved method and a time-resolved method are performed on tissue-like phantoms. Aqueous suspension of microspheres and Intralipid are used as scattering media with the addition of ink as an absorber. mu(')(s) is first measured on weakly absorbing media. The robustness of these measurements is then tested with respect to a variation of mu(a). The spatially resolved method gave more accurate estimates for mu(')(s) whereas the time-resolved method gave better results for mu(a) estimates.

Analytical method for localizing a fluorescent inclusion in a turbid medium

We describe a novel method for localizing a fluorescent inclusion in a homogeneous turbid medium through the use of time-resolved techniques. Based on the calculation of the mean time of the fluorescence curves, the method does not require a priori knowledge of either the fluorescence lifetime or the mean time of the instrument response function since it adopts a differential processing approach. Theoretical expressions were validated and experiments for assessing the accuracy of localization were carried out on liquid optical phantoms with a small fluorescent inclusion. The illumination and detection optical fibers were immersed in the medium to achieve infinite medium geometry as required by the model used. The experimental setup consisted of a time-correlated single-photon counting system. Submillimeter accuracy was achieved for the localization of the inclusion.