Measurement of the coherent transfer function

It is important to characterize and quantify the imaging performance of optical coherence tomography (OCT) and especially optical coherence microscopy (OCM) systems for validation and calibration of the setup. We suggest using a simple rubber surface, presented to the setup at different focal positions, to measure the system’s lateral resolution and depth of field. Rubber is a nearly perfect scatterer and produces a fully developed speckle pattern. The speckle size and shape are directly defined by the optical system and its illumination and detection apertures. Analysis of the speckle pattern reveals the system’s coherent transfer function (CTF) and hence the point spread function (PSF) in a straightforward computation.

Coherent transfer functions and extended depth of field

Theo Lasser, Rainer Leitgeb, Christophe Pache, Martin Villiger

To preserve the speed advantage of Fourier Domain detection in Optical Coherence Microscopy (OCM), extended depth of field is needed. With a narrow probing volume that extends over a long axial range, tissue could be measured in vivo and at cellular resolution. To assess and improve the DOF and the lateral resolution, we analyzed the coherent transfer function (CTF) of OCM. Both the illumination and detection optics contribute equally to the overall imaging performance. In the Fourier domain detection, each pixel of the spectrometer has its specific CTF, sampling a different region of the object’s spatial frequency spectrum. For classical optics and increasing numerical apertures these regions start to overlap and bend, which limits the depth of field. Annular apertures, created with Bessel-like beams produced by axicon lenses or phase filters, circumvent these detrimental effects, but introduce strong side lobes. Decoupling the detection and the illumination apertures is needed to provide the flexibility in engineering a CTF that optimizes the lateral resolution and the DOF at the same time all while reducing these side lobes. We evaluated different combinations of Gaussian and Bessel-like illumination and detection optics, both theoretically and experimentally. Using Bessel-like beams as well in the illumination as in the detection paths, but with annular apertures of different lobe radii, we obtained a lateral resolution of 1.3 μm and an extended depth of field of more than 300 μm, which was completely decoupled from the numerical aperture and scalable to high lateral resolution.

The International Society for Optical Engineering Spie-Int Soc Optical Engineering2010

Multiple scattering in wide-field optical coherence tomography

Boris Karamata

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.

EPFL2004

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.