Stable carrier generation and phase-resolved digital data processing in optical coherence tomography

Chromatic dispersion based axial length estimation using retinal spectral domain optical coherence tomography

We propose a method for extracting the axial length of the human eye from high-resolution spectral domain optical coherence tomography (SD-OCT) retinal scans. The method evaluates the chromatic dispersion introduced by the anterior segment and the vitreous of the eye. By analyzing sub-spectral scans, we quantify the axial shift caused by dispersion and relate it to the thicknesses of the media passed by the OCT beam. The method depends on accurate k-linearization and pixel-to-wavenumber calibration. First, we demonstrate the feasibility of our approach using a model eye with adjustable water chamber length. Subsequently, the method is explored for in vivo retinal OCT scans. Challenges are inter-subject variability and limited availability of exact chromatic dispersion data for ocular tissues in the relevant spectral range for OCT imaging. By interpolating the refractive indices of an established eye model from visible wavelengths to the infrared wavelengths of the OCT system using the dispersion of water and estimating refined dispersion properties of the lens, we improve the model's agreement with in vivo measurements.

Tunable image-mapping optical coherence tomography

We present tunable image-mapping optical coherence tomography (TIM-OCT), which can provide optimized imaging performance for a given application by using a programmable phase-only spatial light modulator in a low-coherence full-field spectral-domain interferometer. The resultant system can provide either a high lateral resolution or a high axial resolution in a snapshot without moving parts. Alternatively, the system can achieve a high resolution along all dimensions through a multiple-shot acquisition. We evaluated TIM-OCT in imaging both standard targets and biological samples. Additionally, we demonstrated the integration of TIM-OCT with computational adaptive optics in correcting sample-induced optical aberrations.

Investigation of methods to extract confocal function parameters for the depth resolved determination of attenuation coefficients using OCT in intralipid samples, titanium oxide phantoms, and in vivo human retinas

The attenuation coefficient provides a quantitative parameter for tissue characterization and can be calculated from optical coherence tomography (OCT) data, but accurate determination requires compensation for the confocal function. We present extensive measurement series for extraction of the focal plane and the apparent Rayleigh length from the ratios of OCT images acquired with different focus depths and compare these results with two alternative approaches. By acquiring OCT images for a range of different focus depths the optimal focus plane difference is determined for intralipid and titanium oxide (TiO₂) phantoms with different scatterer concentrations, which allows for calculation of the attenuation coefficient corrected for the confocal function. The attenuation coefficient is determined for homogeneous intralipid and TiO₂ samples over a wide range of concentrations. We further demonstrate very good reproducibility of the determined attenuation coefficient of layers with identical scatter concentrations in a multi-layered phantom. Finally, this method is applied to in vivo retinal data.

Inscription and analysis of non-uniform diffraction gratings in azobenzene molecular glass thin films

Non-uniform surface relief diffraction gratings were laser-inscribed on azobenzene molecular glass thin films using a modified Lloyd's mirror interferometer. The azobenzene films were exposed to an adjustable interference pattern produced by the recombination of collimated and spherically divergent laser wave fronts. The localized pitch, grating vector orientation and depth of the resulting non-uniform gratings were measured using an atomic force microscope and a theoretical model was analytically developed to explain the experimental results. The fabricated gratings exhibited a chirping or pitch variation along the imposed X-axis as well as an angular change in the grating vector orientation along the imposed Y-axis. Studies were conducted on various non-uniform grating configurations having central pitches of 500 nm, 1000 nm, 1500 nm and 2000 nm.

Interferometric mapping of material properties using thermal perturbation

Rapid, accurate, and nondestructive mapping of material properties is of great interest in many fields, with applications ranging from detection of defects or other subsurface features in semiconductors to estimating temperature rise in various tissue layers during laser therapy. We demonstrate the speed and precision of two interferometric techniques, quantitative phase imaging and phase-resolved optical coherence tomography, in recording optical phase changes induced by energy deposition in various materials. Such phase perturbations can be used to infer sample properties, ranging from absorption and temperature maps to distribution of electric field or resistivity. We derive the theoretical sensitivity limits of such techniques and demonstrate their applicability to the mapping of absorption coefficients, temperature, and electric fields in synthetic and biological samples.

Optical phase changes induced by transient perturbations provide a sensitive measure of material properties. We demonstrate the high sensitivity and speed of such methods, using two interferometric techniques: quantitative phase imaging (QPI) in transmission and phase-resolved optical coherence tomography (OCT) in reflection. Shot-noise–limited QPI can resolve energy deposition of about 3.4 mJ/cm² in a single pulse, which corresponds to 0.8 °C temperature rise in a single cell. OCT can detect deposition of 24 mJ/cm² energy between two scattering interfaces producing signals with about 30-dB signal-to-noise ratio (SNR), and 4.7 mJ/cm² when SNR is 45 dB. Both techniques can image thermal changes within the thermal confinement time, which enables accurate single-shot mapping of absorption coefficients even in highly scattering samples, as well as electrical conductivity and many other material properties in biological samples at cellular scale. Integration of the phase changes along the beam path helps increase sensitivity, and the signal relaxation time reveals the size of hidden objects. These methods may enable multiple applications, ranging from temperature-controlled retinal laser therapy or gene expression to mapping electric current density and characterization of semiconductor devices with rapid pump–probe measurements.

Multi-shaping technique reduces sidelobe magnitude in optical coherence tomography

Shaping methods that are commonly used in Fourier-domain optical coherence tomography (FD-OCT) can suppress sidelobe artifacts in the axial direction, but they typically broaden the mainlobe of the point spread function (PSF) and reduce the axial resolution. To improve OCT image quality without this tradeoff, we have developed a multi-shaping technique that reduces the axial sidelobe magnitude dramatically and achieves better axial resolution than conventional shaping methods. This technique is robust and compatible in various FD-OCT imaging systems. Testing of multi-shaping in three experimental settings shows that it reduced the axial sidelobe contribution by more than 8 dB and improved the contrast to noise by at least 30% and up to three-fold. Multi-shaping enables accurate image analysis and is potentially useful in many OCT applications.

Twenty-five years of optical coherence tomography: the paradigm shift in sensitivity and speed provided by Fourier domain OCT [Invited]

Optical coherence tomography (OCT) has become one of the most successful optical technologies implemented in medicine and clinical practice mostly due to the possibility of non-invasive and non-contact imaging by detecting back-scattered light. OCT has gone through a tremendous development over the past 25 years. From its initial inception in 1991 [Science254, 1178 (1991)] it has become an indispensable medical imaging technology in ophthalmology. Also in fields like cardiology and gastro-enterology the technology is envisioned to become a standard of care. A key contributor to the success of OCT has been the sensitivity and speed advantage offered by Fourier domain OCT. In this review paper the development of FD-OCT will be revisited, providing a single comprehensive framework to derive the sensitivity advantage of both SD- and SS-OCT. We point out the key aspects of the physics and the technology that has enabled a more than 2 orders of magnitude increase in sensitivity, and as a consequence an increase in the imaging speed without loss of image quality. This speed increase provided a paradigm shift from point sampling to comprehensive 3D in vivo imaging, whose clinical impact is still actively explored by a large number of researchers worldwide.

Quantitative depth-resolved microcirculation imaging with optical coherence tomography angiography (Part Ι): Blood flow velocity imaging

The research goal of the microvascular network imaging with OCT angiography is to achieve depth-resolved blood flow and vessel imaging in vivo in the clinical management of patents. In this review, we review the main phenomena that have been explored in OCT to image the blood flow velocity vector and the vessels of the microcirculation within living tissues. Parameters that limit the accurate measurements of blood flow velocity are then considered. Finally, initial clinical diagnosis applications and future developments of OCT flow images are discussed.

Computational optical coherence tomography [Invited]

Optical coherence tomography (OCT) has become an important imaging modality with numerous biomedical applications. Challenges in high-speed, high-resolution, volumetric OCT imaging include managing dispersion, the trade-off between transverse resolution and depth-of-field, and correcting optical aberrations that are present in both the system and sample. Physics-based computational imaging techniques have proven to provide solutions to these limitations. This review aims to outline these computational imaging techniques within a general mathematical framework, summarize the historical progress, highlight the state-of-the-art achievements, and discuss the present challenges.

Dispersion analysis and measurement of potassium tantalate niobate crystals by broadband optical interferometers

Electro-optic crystals, such as potassium tantalate niobate [KTa_1-xNb_xO₃(KTN)], are enabling materials for many optical devices. Their utility in broadband applications heavily depends on their dispersion property. To this end, an analysis of dispersion mismatch in broadband optical interferometers is first presented. Then a method utilizing polynomial phase fitting to measure the dispersion property of materials composing the arms of an interferometer is introduced. As a demonstration, an interferometry system based on optical coherence tomography (OCT) was built, where, for the first time, the group velocity dispersion of a KTN crystal around 1310 nm was measured and numerically compensated for OCT imaging. Several advantages over a widely used method in OCT, which is based on metric functions, are discussed. The results show the fitting method can provide a more reliable measurement with reduced computation complexity.

GPU-accelerated framework for intracoronary optical coherence tomography imaging at the push of a button

Frequency domain optical coherence tomography (FD-OCT) has become one of the important clinical tools for intracoronary imaging to diagnose and monitor coronary artery disease, which has been one of the leading causes of death. To help more accurate diagnosis and monitoring of the disease, many researchers have recently worked on visualization of various coronary microscopic features including stent struts by constructing three-dimensional (3D) volumetric rendering from series of cross-sectional intracoronary FD-OCT images. In this paper, we present the first, to our knowledge, "push-of-a-button" graphics processing unit (GPU)-accelerated framework for intracoronary OCT imaging. Our framework visualizes 3D microstructures of the vessel wall with stent struts from raw binary OCT data acquired by the system digitizer as one seamless process. The framework reports the state-of-the-art performance; from raw OCT data, it takes 4.7 seconds to provide 3D visualization of a 5-cm-long coronary artery (of size 1600 samples x 1024 A-lines x 260 frames) with stent struts and detection of malapposition automatically at the single push of a button.

Dispersion compensation in Fourier domain optical coherence tomography

In this work, we propose a numerical technique to compensate for errors due to dispersion effects in Fourier domain optical coherence tomography. The proposed technique corrects for errors in depth measurements and resolution loss due to dispersion. The results show that, by using this technique, errors in thickness measurement are reduced from about 5% to less than 0.1% depending on the sample length and the amount of dispersion. Also, an improvement in the resolution from about 50 μm to less than 10 μm is demonstrated.

Feasibility of capillary velocity assessment by statistical means using dual-beam spectral-domain Optical Coherence Tomography: a preliminary study

The assessment of vascular dynamics has been shown to yield both qualitative and quantitative metrics and thus play a pivotal role in the diagnosis and prognosis of various diseases, which may manifest as microcirculatory irregularities. Optical Coherence Tomography (OCT) is an established imaging modality which utilises the principle of optical interferometry to distinguish between spatial changes in refractive index and thus formulate a multi-dimensional representation of a specimen in vivo. Nonetheless, difficulties remain in obtaining accurate data (morphological and/or transient) in an environment which is subject to such large biological variability. In an effort to address the issue of angular dependence as with Doppler based analysis, a dual-beam Spectral-domain OCT system for quasi-simultaneous specimen scanning is described. A statistical based method of phase correlation is outlined which is capable of quantifying velocity values in addition to the ability to discern bidirectionality, without the necessity of angular computation.

Dispersion compensation in Fourier domain optical coherence tomography using the fractional Fourier transform

We address numerical dispersion compensation based on the use of the fractional Fourier transform (FrFT). The FrFT provides a new fundamental perspective on the nature and role of group-velocity dispersion in Fourier domain OCT. The dispersion induced by a 26 mm long water cell was compensated for a spectral bandwidth of 110 nm, allowing the theoretical axial resolution in air of 3.6 μm to be recovered from the dispersion degraded point spread function. Additionally, we present a new approach for depth dependent dispersion compensation based on numerical simulations. Finally, we show how the optimized fractional Fourier transform order parameter can be used to extract the group velocity dispersion coefficient of a material.

Optical coherence tomography-current technology and applications in clinical and biomedical research

Optical coherence tomography (OCT) is a noninvasive imaging technique that provides real-time two- and three-dimensional images of scattering samples with micrometer resolution. By mapping the local reflectivity, OCT visualizes the morphology of the sample. In addition, functional properties such as birefringence, motion, or the distributions of certain substances can be detected with high spatial resolution. Its main field of application is biomedical imaging and diagnostics. In ophthalmology, OCT is accepted as a clinical standard for diagnosing and monitoring the treatment of a number of retinal diseases, and OCT is becoming an important instrument for clinical cardiology. New applications are emerging in various medical fields, such as early-stage cancer detection, surgical guidance, and the early diagnosis of musculoskeletal diseases. OCT has also proven its value as a tool for developmental biology. The number of companies involved in manufacturing OCT systems has increased substantially during the last few years (especially due to its success in opthalmology), and this technology can be expected to continue to spread into various fields of application.

Artefact reduction for cell migration visualization using spectral domain optical coherence tomography

Visualization of cell migration during chemotaxis using spectral domain optical coherence tomography (OCT) requires non-standard processing techniques. Stripe artefacts and camera noise floor present in OCT data prevent detailed computer-assisted reconstruction and quantification of cell locomotion. Furthermore, imaging artefacts lead to unreliable results in automated texture based cell analysis. Here we characterize three pronounced artefacts that become visible when imaging sample structures with high dynamic range, e.g. cultured cells: (i) time-varying fixed-pattern noise; (ii) stripe artefacts generated by background estimation using tomogram averaging; (iii) image modulations due to spectral shaping. We evaluate techniques to minimize the above mentioned artefacts using an 800 nm optical coherence microscope. Effect of artefact reduction is shown exemplarily on two cell cultures, i.e. Dictyostelium on nitrocellulose substrate, and retinal ganglion cells (RGC-5) cultured on a glass coverslip. Retinal imaging also profits from the proposed processing techniques.

Signal processing for sidelobe suppression in optical coherence tomography images

In an optical coherence tomography system, the sidelobes of the point-spread function (PSF) introduced from the optical source reduce the A-scan imaging resolution and contrast of the images. A gradual iterative signal subtraction method based on the study of a point signal influenced by other points with different distances through the PSF is proposed in this paper. Comparing with the CLEAN algorithm and two typical deconvolution methods, the processed results demonstrate this algorithm can reduce sidelobes effectively with the least runtime. It is also found that it is insensitive to noise while slightly improving the longitudinal resolution, which shows this algorithm is good for improving image quality.

High speed optical coherence microscopy with autofocus adjustment and a miniaturized endoscopic imaging probe

Optical coherence microscopy (OCM) is a promising technique for high resolution cellular imaging in human tissues. An OCM system for high-speed en face cellular resolution imaging was developed at 1060 nm wavelength at frame rates up to 5 Hz with resolutions of < 4 microm axial and < 2 microm transverse. The system utilized a novel polarization compensation method to combat wavelength dependent source polarization and achieve broadband electro-optic phase modulation compatible with ultrahigh axial resolution. In addition, the system incorporated an auto-focusing feature that enables precise, near real-time alignment of the confocal and coherence gates in tissue, allowing user-friendly optimization of image quality during the imaging procedure. Ex vivo cellular images of human esophagus, colon, and cervix as well as in vivo results from human skin are presented. Finally, the system design is demonstrated with a miniaturized piezoelectric fiber-scanning probe which can be adapted for laparoscopic and endoscopic imaging applications.

Optical coherence tomography - development, principles, applications

This paper presents a review of the development of optical coherence tomography (OCT), its principles and important applications. Basic OCT systems are described and the physical foundations of OCT signal properties and signal recording systems are reviewed. Recent examples of OCT applications in ophthalmology, cardiology, gastroenterology and dermatology outline the relevance of this advanced imaging modality in the medical field.

Dispersion encoded full range frequency domain optical coherence tomography

We propose an iterative algorithm that exploits the dispersion mismatch between reference and sample arm in frequency-domain optical coherence tomography (FD-OCT) to effectively cancel complex conjugate mirror terms in individual A-scans and thereby generate full range tomograms. The resulting scheme, termed dispersion encoded full range (DEFR) OCT, allows distinguishing real structures from complex conjugate mirror artifacts. Even though DEFR-OCT has higher post-processing complexity than conventional FD-OCT, acquisition speed is not compromised since no additional A-scans need to be measured, thereby rendering this technique robust against phase fluctuations. The algorithm uses numerical dispersion compensation and exhibits similar resolution as standard processing. The residual leakage of mirror terms is further reduced by incorporating additional knowledge such as the power spectrum of the light source. The suppression ratio of mirror signals is more than 50 dB and thus comparable to complex FD-OCT techniques which use multiple A-scans.

[Optical coherence tomography by evaluation of the spectral modulation of white light interferometry]

Amplitude and phase factor of a wave scattered from an object relative to the amplitude and phase factor of the incident wave reveals the spatial distribution of the object's scattering power by means of Fourier transform. In this paper this relation is discussed and the dependence of the spatial resolution and the maximal range of depth on the physical parameters as well as the influence of dispersion are investigated. The signal-to-noise ratio yields a measure on the sensitivity of the method which is put into relation of the sensitivity of Time-Domain Optical Coherence Tomography (TD-OCT). A white light interferometer is discussed. The spectrum of the interference between the wave scattered at an object and the reference wave is observed and yields the spectral distribution of the relative amplitude and phase factor.

Adaptive spectral apodization for sidelobe reduction in optical coherence tomography images

Often in imaging systems, the bandpass of the system is not uniform. In temporal coherence imaging methods such as optical coherence tomography, one would like to achieve the most spatially confined impulse response possible with a given source spectrum, minimizing sidelobes that blur adjacent features together. Typically the spectrum of the source is controlled in order to remove sidelobes from the measured interferogram. However, the measured interferogram is not necessarily the best estimate of the scattering density of the object. In this work, a sidelobe supression method is proposed and demonstrated to achieve low sidelobes even with highly nonuniform, non-Gaussian spectra.

Ultrahigh-resolution, high-speed, Fourier domain optical coherence tomography and methods for dispersion compensation

Ultrahigh-resolution optical coherence tomography uses broadband light sources to achieve axial image resolutions on the few micron scale. Fourier domain detection methods enable more than an order of magnitude increase in imaging speed and sensitivity, thus overcoming the sensitivity limitations inherent in ultrahigh-resolution OCT using standard time domain detection. Fourier domain methods also provide direct access to the spectrum of the optical signal. This enables automatic numerical dispersion compensation, a key factor in achieving ultrahigh image resolutions. We present ultrahigh-resolution, high-speed Fourier domain OCT imaging with an axial resolution of 2.1 ìm in tissue and 16,000 axial scans per second at 1024 pixels per axial scan. Ultrahigh-resolution spectral domain OCT is shown to provide a ~100x increase in imaging speed when compared to ultrahigh-resolution time domain OCT. In vivo imaging of the human retina is demonstrated. We also present a general technique for automatic numerical dispersion compensation, which is applicable to spectral domain as well as swept source embodiments of Fourier domain OCT.

Phase noise and phase modulation in optical coherence tomography

The signal in optical coherence tomography is often modulated either in phase or by use of the Doppler modulation generated by a depth-scanning mechanism. The effect of each type of modulation on the signal's amplitude is evaluated. The advantages of each type of modulation in terms of immunity to phase noise and penetration depth are discussed in relation to two envelope detection schemes, i.e., lock-in detection and rms-to-dc conversion. Phase noise due to drifts and demodulation instabilities causes distortion of the signal envelope and can be responsible in part for the speckle appearance of the image.

In vivo dynamic human retinal blood flow imaging using ultra-high-speed spectral domain optical coherence tomography

An ultra-high-speed spectral domain optical Doppler tomography (SD-ODT) system is used to acquire images of blood flow in a human retina in vivo, at 29,000 depth profiles (A-lines) per second and with data acquisition over 99% of the measurement time. The phase stability of the system is examined and image processing algorithms are presented that allow accurate determination of bi-directional Doppler shifts. Movies are presented of human retinal flow acquired at 29 frames per second with 1000 A-lines per frame over a time period of 3.28 seconds, showing accurate determination of vessel boundaries and time-dependent bi-directional flow dynamics in artery-vein pairs. The ultra-high-speed SD-ODT system allows visualization of the pulsatile nature of retinal blood flow, detects blood flow within the choroid and retinal capillaries, and provides information on the cardiac cycle. In summary, accurate video rate imaging of retinal blood flow dynamics is demonstrated at ocular exposure levels below 600 microW.

Phase-slope and group-dispersion calculations in the frequency domain by simple optical low-coherence reflectometry

We report a new method by which phase slope and group dispersion can be calculated with a simple optical low-coherence reflectometer to quantify physiological conditions. A discrete-time signal processing algorithm based on the first and second derivatives of the phase with respect to wave number was developed from discrete-time Fourier properties. The algorithm avoids the 2pi ambiguity associated with most phase unwrapping. Experimental data collected by use of well-characterized optical materials validated the algorithm, which was minimally sensitive to phase noise. The group dispersion of bovine cornea was measured at various hydrations and was significantly dependent on hydration. The results suggest that group dispersion is an indicator of corneal alterations.

High speed, wide velocity dynamic range Doppler optical coherence tomography (Part I): System design, signal processing, and performance

Improvements in real-time Doppler optical coherence tomography (DOCT), acquiring up to 32 frames per second at 250 x 512 pixels per image, are reported using signal processing techniques commonly employed in Doppler ultrasound imaging. The ability to measure a wide range of flow velocities, ranging from less than 20 microm/s to more than 10 cm/s, is demonstrated using an 1.3 microm DOCT system with flow phantoms in steady and pulsatile flow conditions. Based on full implementation of a coherent demodulator, four different modes of flow visualization are demonstrated: color Doppler, velocity variance, Doppler spectrum, and power Doppler. The performance of the former two, which are computationally suitable for real-time imaging, are analyzed in detail under various signal-to-noise and frame-rate conditions. The results serve as a guideline for choosing appropriate imaging parameters for detecting in vivo blood flow.

Real-time multi-functional optical coherence tomography

We demonstrate real-time acquisition, processing, and display of tissue structure, birefringence, and blood flow in a multi-functional optical coherence tomography (MF-OCT) system. This is accomplished by efficient data processing of the phase-resolved inteference patterns without dedicated hardware or extensive modification to the high-speed fiber-based OCT system. The system acquires images of 2048 depth scans per second, covering an area of 5 mm in width x 1.2 mm in depth with real-time display updating images in a rolling manner 32 times each second. We present a video of the system display as images from the proximal nail fold of a human volunteer are taken.

Real-time dispersion compensation in scanning interferometry

We propose and demonstrate a method of real-time dispersion compensation suitable for scanning interferometry and optical coherence tomography. Static grating tilt in a scanning frequency-domain optical delay line is shown to produce dispersion that is linearly proportional to scan position, and we use this property to achieve depth-dependent dispersion compensation during an interferometric scan through a dispersive sample.

Simultaneous intensity, birefringence, and flow measurements with high-speed fiber-based optical coherence tomography

We demonstrate that tissue structure, birefringence, and blood flow can be imaged simultaneously by use of techniques of polarization-sensitive optical coherence tomography and phase-resolved optical Doppler tomography. An efficient data-acquisition procedure is implemented that optimizes the concurrent processing and display of all three image types. Images of in vivo human skin acquired with a high-speed fiber-based system are presented.

Spectral shaping for non-Gaussian source spectra in optical coherence tomography

We present a digital spectral shaping technique to reduce the sidelobes (ringing) of the axial point-spread function in optical coherence tomography for non-Gaussian-shaped source spectra. The spectra of two superluminescent diodes were combined to generate a spectrum with significant modulation. Images of onion cells demonstrate the improved image quality in a turbid biological sample. A quantitative analysis of the accompanying penalty in signal-to-noise ratio is given.