Multiple scattering in wide-field optical coherence tomography

Optical Coherence Tomography (OCT), a well-established imaging method based on low-coherence interferometry, provides cross-sectional images of the internal structure of biological samples with a resolution in the micrometer range. OCT was successfully applied on various tissues such as for instance the retina, the skin or a tooth. In highly scattering tissues like the skin, probing depth is limited to approximately 2mm, mainly due to insufficient rejection of multiply scattered light. Presently, the contribution of multiple scattering in OCT is not fully understood. Therefore, there is a strong and urgent need to develop models allowing a reliable evaluation of the system's limitations as well as the improvement of the imaging capabilities. It is generally believed that a relevant model should account for loss of correlation between the reference and the sample field due to multiple scattering. We developed a new comprehensive model of OCT. Our preliminary study revealed that the reference and sample fields are actually fully correlated. This important result allowed us to model the OCT signal as a sum of stationary random phasors and treated it as a statistical signal. The mean of this signal can be calculated thanks to classical results of statistical optics and to a Monte Carlo simulation. Unlike other existing models, our model accounts for the source autocorrelation function. The model proved to be in excellent agreement with a whole range of experimental data gathered in a comprehensive study of cross-talk in wide-field OCT. Moreover, our results put in question the applicability of widely used models of OCT based on the "extended Huygens-Fresnel principle", which assume a partial correlation between interfering fields due to multiple scattering. The construction of conventional OCT images is based on lateral scanning of a beam focused within the sample. To increase image acquisition speed and eliminate the need for lateral scanning, wide-field OCT was recently developed. Our experimental and theoretical investigations of the potential and limitations of wide-field OCT revealed the crucial role played by the spatial coherence of the light source. Spatially coherent illumination generates considerable coherent optical cross-talk, which prevents shot-noise-limited detection and diffraction-limited imaging in scattering samples. The dependence on several parameters of the optical system and of the sample properties was investigated in a comprehensive study. Cross-talk increases with the wide-field diameter, numerical aperture, source coherence length, and sample optical density; and strongly depends on sample anisotropy. We showed that spatially incoherent illumination realized with a thermal light source permits cross-talk suppression in wide-field OCT, i.e. rejection of multiply scattered light to a level comparable to that of point scanning OCT. We performed a theoretical study which revealed that the power per spatial mode radiated by thermal light sources is too low to permit a high signal-to-noise ratio while maintaining a fast acquisition speed. Therefore, wide-field OCT realized with either spatially coherent or spatially incoherent illumination suffers from inherent fundamental limitation. This led us to investigate the possibility of exploiting a spatially incoherent light source brighter than a thermal light source. We came to the conclusion that such a "pseudothermal" light source can potentially lead to wide-field OCT systems devoid of cross-talk and an image acquisition speed higher than that of a point scanning OCT system. However, the attractive properties of pseudothermal light sources could be gained at the expense of the simplicity and the economical advantages offered by thermal light sources. Furthermore, fast acquisition speed also relies on a performing "smart pixel detector array". Presently, such detectors do not have sufficient sensitivity and their frequency readout is too low as shown in our feasibility study. We performed a theoretical investigation of the potential of thermal light sources in terms of axial resolution and power per mode. The former revealed that the maximum power per mode is radiated at a wavelength higher than the spectral peak of a blackbody radiator. This led to the important practical conclusion that, at 6000 K, the maximum power is collected in the therapeutic spectral window in OCT (600 - 1300 nm), while at 3000 K this peak is shifted out of the therapeutic window leading to significant power losses. More generally, our work provides a design tool for choosing the optimal thermal light source for a given therapeutic window in terms of signal-to-noise ratio. Currently available sources at 6000 K consist of high pressure gas arc lamps providing a spectrum endowed with spectral lines deleterious for OCT. By suppressing a portion of the spectrum devoid of spectral lines of a mercury arc lamp, we obtained amongst the highest axial resolution so far reported in OCT. Furthermore, the importance of the speckle statistics in OCT incited us to clarify the origin of a difference between two theoretical results reported in the literature. Indeed, two calculations of the amplitude distribution of speckles in OCT, each of them based on a different mathematical formulation, yield different results. We showed that a modification of an initial assumption in one of the formulation leads to equivalent results. In conclusion, this thesis provides a deeper understanding of the potential and limitations of widefield OCT, leading to important design rules. Moreover, it presents a new comprehensive model of OCT putting in question other widely used models.

Endoscopic low coherence interferometry applied to tri-dimensional contouring of upper airways

Yves Delacrétaz

This thesis presents a novel approach for optical tri-dimensional contouring, requiring only one acquisition to extract a contour depth. It is based on an interferometric setup, using a reduced coherence light source, and is fully adaptable to endoscopic procedure commonly in used in ear-nose-throat (ENT) medicine. It has been demonstrated that the method can be successfully applied to imaging of porcine upper airways. Although this method uses interferometry to obtain depth range measurement, the phase signal is not directly evaluated, rather the location of the coherent superposition of two waves is extracted, and thus provides a robust technique even in harsh environment. The key feature is to extract the coherent area on the image by locating intentionally induced fringes created by off-axis superposition of a smooth reference wave and an object wave coming from a rough or diffusing sample through an optical system with limited aperture. A model for the propagation of the light through the optical system has been developed. It takes into account the broadened emission spectrum of the light source, as well as the statistical speckle effect induced by the illumination of a rough surface with partially coherent light. Simulated results have been carefully compared with measurements, both in the spatial domain and in the Fourier domain (spatial frequencies response), which has permitted to better understand the processes involved in the creation of the interference fringes, and improve the design of the device. Different filters have been investigated in order to extract the depth information from acquired interferograms. By using Gabor filterbanks, which do not need any a priori knowledge of the signal, selectivity of the extracted signal has been enhanced by a factor of 1.3 compared to a simple local Fourier spectrum evaluation, and time of the extraction procedure divided by 6, even when compared to fast iterative calculation in the Fourier spectrum. Moreover, by imposing the frequency and/or orientation of the signal to detect, which can be known by a simple calibration procedure, it has been proved that the selectivity can be further enhanced by a factor of 2.5, keeping the calculation time as short or even shorter. A compact prototype has been designed and built, with a specific multiple fibers arrangement for illumination. The system is separated in two parts: the first one is directly attached to the rigid endoscope, and contains imaging as well as recombining object/reference wave optics, and the detector. It fits into a box with overall dimensions of 115 × 83 × 94 mm. The second part contains the laser source, the scanning arm and the injection optics for the fibered system. It is designed to be left on a small cart. This device can potentially be brought to the bedside, and used during routine endoscopic check-up. Associated with the device, a complete framework has been elaborated in order to simplify the acquisition procedure as well as the signal processing. This results in a fully automated environment taking in charge hardware control, data storage and signal processing, up to the tri-dimensional scene rendering.

EPFL2009

Multiple scattering in optical coherence tomography. I. Investigation and modeling

Stéphane Bourquin, Boris Karamata, Theo Lasser, Markus Laubscher, Marcel Leutenegger

We present a new model of optical coherence tomography (OCT) taking into account multiple scattering. A theoretical analysis and experimental investigation reveals that in OCT, despite multiple scattering, the field backscattered from the sample is generally spatially coherent and that the resulting interference signal with the reference field is stationary relative to measurement time. On the basis of this result, we model an OCT signal as a sum of spatially coherent fields with random-phase arguments-constant during measurement time-caused by multiple scattering. We calculate the mean of such a random signal from classical results of statistical optics and a Monte Carlo simulation. OCT signals predicted by our model are in very good agreement with a depth scan measurement of a sample consisting of a mirror covered with an aqueous suspension of microspheres. We discuss other comprehensive OCT models based on the extended Huygens-Fresnel principle, which rest on the assumption of partially coherent interfering fields. (C) 2005 Optical Society of America.

2005

Low-coherence interferometry based on customized detector arrays

Stéphane Bourquin

Many biological tissues and a lot of technical objects of interest are to a certain degree transparent to light in the visible and near infrared wavelength range of the electromagnetic spectrum. Nevertheless, it has not been possible for a long time to create clear optical images of these objects due to the multiple scattering of a large part of the employed light. The recent development of low-coherence interferometry (LCI) forms the basis for selecting just the ballistic photons from the back-scattered or transmitted light, and hence enables optical low-coherence reflectometry (OLCR) and optical coherence tomography (OCT) for profiling or imaging with micrometer resolution. All published work until now on OLCR or OCT for cross-sectional imaging relies on point detectors, single-channel electronics and mechanical scanning for 2-dimensional or 3-dimensional imaging, limiting the image acquisition rate to 4-8 per second. This work is devoted to improving the state of the art by enhancing the acquisition speed of OLCR 2D and 3D imaging by systematically exploiting the parallelism offered by custom image sensors with smart pixels based on CMOS technology and by developing a new optical setup without moving parts. We replace the conventional single-channel approach by a multi-channel approach using an array of smart pixels, where each pixel consists of a high-sensitivity photodetector with an offset subtraction circuitry, an amplifier, and an envelope detector. Although dated 2 μm CMOS technology is used for the fabrication, our layout is so compact that a pixel period of only 110 μm is achieved. This makes it possible to integrate linear arrays as well as area arrays of smart pixels. We have successfully demonstrated a 64-pixel line sensor and a 58´58 pixel image sensor in an optical low-coherence reflectometry (OLCR) system, making use of versatile, programmable, high-speed readout electronics we developed. This allows the acquisition of OLCR pixel information in arbitrary sequences for efficient region-of-interest scans. When reading out at a speed of 1 MHz, we can acquire volumetric OLCR data at a repetition rate of 6 Hz, which can be easily improved to the video rate of 25 Hz. The simultaneous photosensing and heterodyne signal detection performed by our smart pixels is carried out close to the photon shot noise limit. In order to simplify OLCR and to reduce system cost, we developed a system based on a novel optical setup with a reflection echelon placed in the reference arm of the interferometer. In combination with the 2D customized high-sensitivity image sensor, this allows high-speed depth cross-sectional measurements without any moving parts, opening the way to cost-effective applications of OLCR in medicine, biology and industry. This thesis has successfully demonstrated that the developed parallel approach to OLCR imaging with CMOS-based smart pixel imagers and a novel optical setup without moving parts is highly promising, representing an important contribution to increasing the practical impact of OLCR imaging.