Derivation of a Monte Carlo method for modeling heterodyne detection in optical coherence tomography systems

A Monte Carlo (MC) method for modeling optical coherence tomography (OCT) measurements of a diffusely reflecting discontinuity embedded in a scattering medium is presented. For the first time to the authors' knowledge it is shown analytically that the applicability of an MC approach to this optical geometry is firmly justified, because, as we show, in the conjugate image plane the field reflected from the sample is delta-correlated from which it follows that the heterodyne signal is calculated from the intensity distribution only. This is not a trivial result because, in general, the light from the sample will have a finite spatial coherence that cannot be accounted for by MC simulation. To estimate this intensity distribution adequately we have developed a novel method for modeling a focused Gaussian beam in MC simulation. This approach is valid for a softly as well as for a strongly focused beam, and it is shown that in free space the full three-dimensional intensity distribution of a Gaussian beam is obtained. The OCT signal and the intensity distribution in a scattering medium have been obtained for several geometries with the suggested MC method; when this model and a recently published analytical model based on the extended Huygens-Fresnel principle are compared, excellent agreement is found.     

High-sensitivity determination of birefringence in turbid media with enhanced polarization-sensitive optical coherence tomography

Polarization-sensitive optical coherence tomography provides high-resolution cross-sectional characterization of birefringence in turbid media. Weakly birefringent biological tissues such as the retinal nerve fiber layer (RNFL) require advanced speckle noise reduction for high-sensitivity measurement of form birefringence. We present a novel method for high-sensitivity birefringence quantification by using enhanced polarization-sensitive optical coherence tomography (EPS-OCT) and introduce the polarimetric signal-to-noise ratio, a mathematical tool for analyzing speckle noise in polarimetry. Multiple incident polarization states and non-linear fitting of normalized Stokes vectors allow determination of retardation with +/-1 degrees uncertainty with invariance to unknown unitary polarization transformations. Results from a weakly birefringent turbid film and in vivo primate RNFL are presented. In addition, we discuss the potential of EPS-OCT for noninvasive quantification of intracellular filamentous nanostructures, such as neurotubules in the RNFL that are lost during the progression of glaucoma.     

Variables affecting polarization-sensitive optical coherence tomography imaging examined through the modeling of birefringent phantoms

Characteristics of tissue birefringence provide valuable diagnostic information when optical coherence tomography (OCT) imaging is performed. We mathematically model single-detector polarization-sensitive (PS) OCT (PS-OCT) for imaging an alternating layered phantom with birefringence. The theoretical analysis and the simulated data all reveal complexities associated with PS-OCT imaging. Various tissue parameters affect the retrieved PS-OCT signal and hence the image features. A better understanding of the mechanisms that control PS-OCT imaging would be useful for interpretation of tissue images and ultimately useful to provide more diagnostic information.     

Monte Carlo modeling of spatial coherence: free-space diffraction

We present a Monte Carlo method for propagating partially coherent fields through complex deterministic optical systems. A Gaussian copula is used to synthesize a random source with an arbitrary spatial coherence function. Physical optics and Monte Carlo predictions of the first- and second-order statistics of the field are shown for coherent and partially coherent sources for free-space propagation, imaging using a binary Fresnel zone plate, and propagation through a limiting aperture. Excellent agreement between the physical optics and Monte Carlo predictions is demonstrated in all cases. Convergence criteria are presented for judging the quality of the Monte Carlo predictions.     

Semi-analytical Monte Carlo simulation for time-resolved light propagating in multilayered turbid media

In this study, a Monte Carlo (MC) method for time-resolved light scattering from multilayered turbid media (SMCML) has been developed. This method is particularly suitable for simulating light backscattering from layered media and receiving the time-resolved signal in a finite sensor area, such as ocean detection, photomedicine and photobiology. The classical semi-analytical MC method requires the scattering events to be located in a single-layer medium. To address the multilayer problem, the energy loss mechanism of photons propagating in tissue was analyzed in this study. According to the energy contribution to the detector, only photons that contribute significantly were considered. Simulations were conducted for stochastic turbid media with different optical parameters. Temporal profiles of the echo signal were obtained with a satisfactory convergence. Compared to the classical MC method, the SMCML method can dramatically reduce the computation time by more than two orders of magnitude.     

Ensemble-averaged Shack-Hartmann wave-front sensing for imaging through turbid media

A technique is described for ensemble-averaging the light wave emerging from a turbid medium, enabling the recovery of optical information that is otherwise lost in a speckle pattern. The technique recovers both an amplitude and a phase function for a wave that has been corrupted by severe scattering, without the use of holography. With the phase estimated, an ensemble-averaged field is constructed that can be backprojected to form an image of the object obscured by the scattering medium. Experimental results suggest that the technique can resolve two object points whose signals are unresolved on the exit surface of a diffuser.     

Multiscan time-domain optical coherence tomography for retina imaging

A versatile time-domain optical coherence tomography system is presented that can generate cross-sectional images by using either transverse priority or depth priority scanning. This is made possible by using a transmissive scanning delay line compatible with balance detection operating at a speed similar to that of the transverse scanner used to scan the beam across the target. In vivo images from the retina are generated and shown using the same system switched to either transverse or depth priority scanning regime, by using the scanning delay line either in slow or fast scanning modes, respectively. A comparative analysis of different scanning regimes depending on image size to fit different areas to be imaged is presented. Safety thresholds due to the different continuous irradiation time per transverse pixel in different scanning regimes are also considered. We present the maximum exposure level for a variety of scanning procedures, employing either A scanning (depth priority) or T scanning (transverse priority) when generating cross-sectional images, en face images, or collecting 3D volumes.     

Polarization-sensitive optical coherence tomography for imaging human atherosclerosis

Polarization-sensitive optical coherence tomography (PS-OCT) combines the advantages of OCT with image contrast enhancement, which is based on its ability to detect phase retardation and the fast-axis angle. Both PS-OCT images and histopathology have demonstrated similar features that allowed differentiation of atherosclerotic structures (i.e., plaques) from normal tissue. Moreover, the picrosirius polarization method was used to confirm PS-OCT assessment of collagen in the fibrous cap of atherosclerotic plaques, and high-frequency (40 MHz) ultrasound images were used to identify calcium in the vessel wall. Our preliminary ex vivo investigation of human aortic specimens indicated that PS-OCT might help to identify atherosclerotic lesions.     

Characterization of optical fiber imaging bundles for swept-source optical coherence tomography

Fiber imaging bundles have been investigated for use in endoscopic optical coherence tomography (OCT) systems, to obviate the requirement for scanning components within the endoscope probe section. Images have been acquired using several optical configurations, two of which are common path in design. Configurations have been selected as having potential for miniaturization and inclusion in endoscopic-type systems, since the advantages of employing imaging bundles are most clearly seen in this type of system. The various types of bundle available are described, and the properties of the leached bundles used here are discussed in detail, with reference to their effect upon the performance of OCT systems. Images are displayed from measurements made on a range of samples.     

Ultrasound-modulated optical tomography of absorbing objects buried in dense tissue-simulating turbid media

Continuous-wave ultrasonic modulation of scattered laser light was used to image objects buried in tissue-simulating turbid media. The buried object had an absorption coefficient greater than the background turbid medium. The ultrasonic wave that was focused into the turbid media modulated the laser light that passed through the ultrasonic field. The modulated laser light that was collected by a photomultiplier tube reflected the local mechanical and optical properties in the zone of ultrasonic modulation. Objects buried in the middle plane of 5-cm-thick dense turbid media were imaged with millimeter resolution through the scanning and detecting alterations of the ultrasound-modulated optical signal. The optical properties of the dense turbid media included an absorption coefficient of 0.1 cm(-1) and a reduced scattering coefficient of 10 cm(-1) and were comparable with those of biological tissues in the visible and near-IR ranges. The dependence of the ultrasound-modulated optical signal on the off-axis distance of the detector from the optic axis and the area of the detector was studied as well.     

Wide-field optical coherence tomography: imaging of biological tissues

We describe a two-dimensional optical coherence tomography technique with which we were able to obtain multiple longitudinal slices of a biological sample directly in a single Z scan. The system is based on a femtosecond Cr4+:forsterite laser and an infrared camera for wide-field imaging of the sample with a depth resolution of 5 microm. With this imaging apparatus we were able to investigate human skin and mouse ear samples and to observe the different constitutive tissues.     

Cancellation of coherent artifacts in optical coherence tomography imaging

Coherent artifacts in optical coherence tomography (OCT) images can severely degrade image quality by introducing false targets if no targets are present at the artifact locations. Coherent artifacts can also add constructively or destructively to the targets that are present at the artifact locations. This constructive or destructive interference will result in cancellation of the true targets or in display of incorrect echo amplitudes of the targets. We introduce the use of a nonlinear deconvolution algorithm, CLEAN, to cancel coherent artifacts in OCT images of extracted human teeth. The results show that CLEAN can reduce the coherent artifacts to the noise background, sharpen the air-enamel and enamel-dentin interfaces, and improve the image contrast.     

Scaling method for fast Monte Carlo simulation of diffuse reflectance spectra from multilayered turbid media

A scaling Monte Carlo method has been developed to calculate diffuse reflectance from multilayered media with a wide range of optical properties in the ultraviolet-visible wavelength range. This multilayered scaling method employs the photon trajectory information generated from a single baseline Monte Carlo simulation of a homogeneous medium to scale the exit distance and exit weight of photons for a new set of optical properties in the multilayered medium. The scaling method is particularly suited to simulating diffuse reflectance spectra or creating a Monte Carlo database to extract optical properties of layered media, both of which are demonstrated in this paper. Particularly, it was found that the root-mean-square error (RMSE) between scaled diffuse reflectance, for which the anisotropy factor and refractive index in the baseline simulation were, respectively, 0.9 and 1.338, and independently simulated diffuse reflectance was less than or equal to 5% for source-detector separations from 200 to 1500 microm when the anisotropy factor of the top layer in a two-layered epithelial tissue model was varied from 0.8 to 0.99; in contrast, the RMSE was always less than 5% for all separations (from 0 to 1500 microm) when the anisotropy factor of the bottom layer was varied from 0.7 to 0.99. When the refractive index of either layer in the two-layered tissue model was varied from 1.3 to 1.4, the RMSE was less than 10%. The scaling method can reduce computation time by more than 2 orders of magnitude compared with independent Monte Carlo simulations.     

Influence of size parameter and refractive index of the scatterer on polarization-gated optical imaging through turbid media

The influence of incident polarized light, refractive index, and size parameter of the scatterer on achievable resolution and contrast (image quality) of polarization-gated transillumination imaging in turbid media is reported here. Differential polarization detection led to significant improvement of image quality of an object embedded in a medium of small-sized scatterers (diameter Dor=lambda,g>or=0.7), the improvement in image quality was less pronounced using either linear or circular polarization gating when the refractive index of the scatterer was high (ns=1.59), but for a lower value of refractive index (ns=1.37), image quality improved with the differential circular polarization gating. We offer a plausible explanation for these observations.     

Primate retina imaging with polarization-sensitive optical coherence tomography.

Polarization-sensitive optical coherence tomography (PSOCT) is applied to determine the depth-resolved polarization state of light backreflected from the eye. The birefringence of the retinal nerve fiber layer (RNFL) was observed and measured from PSOCT images recorded postmortem in a Rhesus monkey. An image-processing algorithm was developed to identify birefringent regions in acquired PSOCT retinal images and automatically determine the thickness of the RNFL. Values of the RNFL thickness determined from histology and PSOCT were compared. PSOCT may provide a new method to determine RNFL thickness and birefringence for glaucoma diagnostics.     

Autocorrelation of scattered laser light for ultrasound-modulated optical tomography in dense turbid media

Based on measurement of the intensity autocorrelation function, a new method to determine the modulation depth of scattered laser light modulated by an ultrasonic wave in turbid media was applied to ultrasound-modulated optical tomography. Good signal-to-noise ratios and high sensitivities were demonstrated. Images of double optically absorbing objects buried in a highly optically scattering gel sample were obtained. The contrast was more than 10%, and the spatial resolution was approximately 2 mm.     

Evaluation of optical coherence quantitation of analytes in turbid media by use of two wavelengths

We introduce a novel method for determining analyte concentration as a function of depth in a highly scattering media by use of a dual-wavelength optical coherence tomography system. We account for the effect of scattering on the measured attenuation by using a second wavelength that is not absorbed by the sample. We assess the applicability of this technique by measuring the concentration of water in an Intralipid phantom, using a probe wavelength of 1.53 mum and a reference wavelength of 1.31 mum. The results of our study show a strong correlation between the measured absorption and the water content of the sample. The accuracy of the technique, however, was limited by the dominance of scattering over absorption in the turbid media. Thus, although the effects of scattering were minimized, significant errors remained in the calculated absorption values. More-accurate results could be obtained with the use of more powerful superluminescent diodes and a choice of wavelengths at which absorption effects are more significant relative to scattering.     

Transillumination optical tomography of tissue-engineered blood vessels: a Monte Carlo simulation

A Monte Carlo technique has been developed to simulate the transillumination laser computed tomography of tissue-engineered blood vessels. The blood vessel was modeled as a single cylinder layer mounted on a tubular mandrel. Sequences of images were acquired while rotating the mandrel. The tomographic image was reconstructed by applying a standard Radon transform. Angular discrimination was applied to simulate a spatial filter, which was used to reject multiply scattered photons. The simulation results indicated that the scattering effect can be overcome with angular discrimination because of the thin tissue thickness. However, any refractive-index mismatch among the tissue, the surrounding media, and the mandrel could produce significant distortions in the reconstructed image.     

Measurement of fluid-flow-velocity profile in turbid media by the use of optical Doppler tomography

Optical Doppler tomography is demonstrated to be a simple, accurate, and noncontact method for measuring the fluid velocity of laminar flow in small-diameter (~0.5-mm) ducts. Studies are described that utilize circular (square) plastic (glass) ducts infused with a moving suspension of polymer microspheres in air and buried in an optically turbid medium. The measurement of Doppler-shifted frequencies of backscattered light from moving microspheres is used to construct a high-resolution spatial profile of fluid-flow velocity in the ducts.     

Monte Carlo model and single-scattering approximation of the propagation of polarized light in turbid media containing glucose

We present a single-scattering model as well as a Monte Carlo model of the effect of glucose on polarized light in turbid media. Glucose alters the Mueller-matrix patterns of diffusely backscattered and forward-scattered light because glucose molecules rotate the polarization plane of linearly polarized light. For example, the angles of rotation in Mueller-matrix elements S21 and S12 are linearly related to the concentration of glucose and increase with the source-detector distance. In the nondiffusion regime, the two models agree well with each other. In the diffusion regime, the single-scattering model is invalid, but there still exists a linear relationship between the angles of rotation in the Mueller-matrix elements and the concentration of glucose, which is predicted by the Monte Carlo model.