Phase-noise analysis of swept-source optical coherence tomography systems

MCOCT: an experimentally and numerically validated, open-source Monte Carlo simulator for optical coherence tomography

Here, we present MCOCT, a Monte Carlo simulator for optical coherence tomography (OCT), incorporating a Gaussian illumination scheme and bias to increase backscattered event collection. MCOCT optical fluence was numerically compared and validated to an established simulator (MCX) and showed concordance at the focus while diverging slightly with distance to it. MCOCT OCT signals were experimentally compared and validated to OCT signals acquired in tissue-mimicking phantoms with known optical properties and showed a similar attenuation pattern with increasing depth while diverging beyond 1.5 mm and proximal to layer interfaces. MCOCT may help in the design of OCT systems for a wide range of applications.

Importance sampling-accelerated simulation of full-spectrum backscattered diffuse reflectance

The Monte Carlo (MC) method is one of the most widely used numerical tools to model the light interaction with tissue. However, due to the low photon collection efficiency and the need to simulate the entire emission spectrum, it is computationally expensive to simulate the full-spectrum backscattered diffuse reflectance (F-BDR). Here, we propose an acceleration scheme based on importance sampling (IS). We derive the biasing sampling function tailored for simulating BDR based on the two-term scattering phase function (TT). The parameters of the TT function at different wavelengths are directly obtained by fitting the Mie scattering phase function. Subsequently, we incorporate the TT function and its corresponding biased function into the redefined IS process and realize the accelerated simulation of F-BDR. Phantom simulations based on the Fourier-domain optical coherence tomography (FD-OCT) are conducted to demonstrate the efficiency of the proposed method. Compared to the original simulator without IS, our proposed method achieves a 373× acceleration in simulating the F-BDR of the multi-layer phantom with a relative mean square error (rMSE) of less than 2%. Besides, by parallelly computing A-lines, our method enables the simulation of an entire B-scan in less than 0.4 hours. To our best knowledge, it is the first time that a volumetric OCT image of a complex phantom is simulated. We believe that the proposed acceleration method can be readily applied to fast simulations of various F-BDR-dependent applications. The source codes of this manuscript are also publicly available online.

Monte Carlo-based full-wavelength simulator of Fourier-domain optical coherence tomography

Monte Carlo (MC) simulation has been widely used to study imaging procedures, including Fourier-domain optical coherence tomography (FD-OCT). Despite the broadband nature of FD-OCT, the results obtained at a single wavelength are often used in previous studies. Some wavelength-relied imaging applications, such as spectroscopic OCT (S-OCT), are unlikely to be simulated in this way due to the lack of information from the entire spectrum. Here, we propose a novel simulator for full-wavelength MC simulation of FD-OCT. All wavelengths within the emission spectrum of the light source will be simulated, and the optical properties derived from Mie theory will be applied. We further combine the inverse discrete Fourier transform (IDFT) with a probability distribution-based signal pre-processing to combat the excessive noises in the OCT signal reconstruction, which is caused by the non-uniform distribution of the scattering events at different wavelengths. Proof-of-concept simulations are conducted to show the excellent performance of the proposed simulator on signal reconstruction and optical properties extraction. Compared with the conventional method, the proposed simulator is more accurate and could better preserve the wavelength-dependent features. For example, the mean square error (MSE) computed between the backscattering coefficient extracted by the proposed simulator and the ground truth is 0.11, which is far less than the value (7.67) of the conventional method. We believe this simulator could be an effective tool to study the wavelength dependency in FD-OCT imaging as well as a preferred solution for simulating spectroscopic OCT.

Suppression of the conjugate signal for broadband computed imaging via synthetic phase modulation

We present synthetic-phase-modulated interferometric synthetic aperture microscopy (SPM-ISAM), a method to perform 3D object reconstructions from data acquired with confocal broadband interferometric microscopy (BIM) that reconstructs images virtually free of coherent and depth-dependent defocus artifacts. This is achieved by implementing a sinusoidal SPM method in combination with an ISAM reconstruction algorithm that uses relatively low-modulation frequencies compared with acquisition frequencies. A theoretical framework and numerical results are provided here.

Generalized Image Reconstruction in Optical Coherence Tomography Using Redundant and Non-Uniformly-Spaced Samples

In this paper, we use Frame Theory to develop a generalized OCT image reconstruction method using redundant and non-uniformly spaced frequency domain samples that includes using non-redundant and uniformly spaced samples as special cases. We also correct an important theoretical error in the previously reported results related to OCT image reconstruction using the Non-uniform Discrete Fourier Transform (NDFT). Moreover, we describe an efficient method to compute our corrected reconstruction transform, i.e., a scaled NDFT, using the Fast Fourier Transform (FFT). Finally, we demonstrate different advantages of our generalized OCT image reconstruction method by achieving (1) theoretically corrected OCT image reconstruction directly from non-uniformly spaced frequency domain samples; (2) a novel OCT image reconstruction method with a higher signal-to-noise ratio (SNR) using redundant frequency domain samples. Our new image reconstruction method is an improvement of OCT technology, so it could benefit all OCT applications.

Stabilizing the phase of swept-source optical coherence tomography by a wrapped Gaussian mixture model

The phase of an optical coherence tomography (OCT) signal carries critical information about particle micro-displacements. However, swept-source OCT (SSOCT) suffers from phase instability problems due to trigger jitters from the swept source. In this Letter, a wrapped Gaussian mixture model (WGMM) is proposed to stabilize the phase of SSOCT systems. A closed-form iteration solution of the WGMM is derived using the expectation-maximization algorithm. Necessary approximations are made for real-time graphic processing unit implementation. The performance of the proposed method is demonstrated through ex vivo, in vivo, and flow phantom experiments. The results show the robustness of the method in different application scenarios.

Effect of spectral leakage on the image formation of Fourier-domain optical coherence tomography

We report on the investigation of spectral leakage's impact on the reconstruction of Fourier-domain optical coherence tomography (FD-OCT). We discuss the shift-variant nature introduced by the spectral leakage and develop a novel spatial-domain FD-OCT image formation model. A proof-of-concept phantom experiment is conducted to validate our model. Compared with previous models, the proposed framework could better describe the image formation process, especially when the fineness of the axial structure approaches the theoretical resolution limit.

Simultaneous Morphological and Flow Imaging Enabled by Megahertz Intravascular Doppler Optical Coherence Tomography

We demonstrate three-dimensional intravascular flow imaging compatible with routine clinical image acquisition workflow by means of megahertz (MHz) intravascular Doppler Optical Coherence Tomography (OCT). The OCT system relies on a 1.1 mm diameter motorized imaging catheter and a 1.5 MHz Fourier Domain Mode Locked (FDML) laser. Using a post processing method to compensate the drift of the FDML laser output, we can resolve the Doppler phase shift between two adjoining OCT A-line datasets. By interpretation of the velocity field as measured around the zero phase shift, the flow direction at specific angles can be qualitatively estimated. Imaging experiments were carried out in phantoms, micro channels, and swine coronary artery in vitro at a speed of 600 frames/s. The MHz wavelength sweep rate of the OCT system allows us to directly investigate flow velocity of up to 37.5 cm/s while computationally expensive phase-unwrapping has to be applied to measure such high speed using conventional OCT system. The MHz sweep rate also enables a volumetric Doppler imaging even with a fast pullback at 40 mm/s. We present the first simultaneously recorded 3D morphological images and Doppler flow profiles. Flow pattern estimation and three-dimensional structural reconstruction of entire coronary artery are achieved using a single OCT pullback dataset.

Optical coherence tomography angiography in preclinical neuroimaging

Preclinical neuroimaging allows for the assessment of brain anatomy, connectivity, and function in laboratory animals, such as mice and this imaging field has been a rapidly growing aimed at bridging the translation gap between animal and human research. The progress in the animal research could be accelerated by high-resolution in vivo optical imaging technologies. Optical coherence tomography-based angiography (OCTA) estimates the scattering from moving red blood cells, providing the visualization of functional micro-vessel networks within tissue beds in vivo without a need for exogenous contrast agents. Recent advancement of OCTA methods have expanded its application to neuroimaging of small animal models of brain disorders. In this paper, we overview the recent development of OCTA techniques for blood flow imaging and its preclinical applications in neuroimaging. In specific, a summary of preclinical OCTA studies for traumatic brain injury, cerebral stroke, and aging brain on mice is reviewed.

Compressed sensing-enabled phase-sensitive swept-source optical coherence tomography

Here we present a novel phase-sensitive swept-source optical coherence tomography (PhS-SS-OCT) system. The simultaneously recorded calibration signal, which is commonly used in SS-OCT to stabilize the phase, is randomly sub-sampled during the acquisition, and it is later reconstructed based on the Compressed Sensing (CS) theory. We first mathematically investigated the method, and verified it through computer simulations. We then conducted a vibrational frequency test and a flow velocity measurement in phantoms to demonstrate the system's capability of handling phase-sensitive tasks. The proposed scheme shows excellent phase stability with greatly discounted data bandwidth compared with conventional procedures. We further showcased the usefulness of the system in biological samples by detecting the blood flow in ex vivo swine left marginal artery. The proposed system is compatible with most of the existing SS-OCT systems and could be a preferred solution for future high-speed phase-sensitive applications.

Phase-stability optimization of swept-source optical coherence tomography

Phase-resolved imaging of swept-source optical coherence tomography (SS-OCT) is subject to phase measurement instabilities involved with the sweep variation of a frequency-swept source. In general, optically generated timing references are utilized to track the variations imposed on OCT signals. But they might not be accurately synchronized due to relative time delays. In this research, we investigated the impact of the signal delays on the timing instabilities and the consequent deviations of the measured phases. We considered two types of timing signals utilized in a popular digitizer operation mode: a sweep trigger from a fiber Bragg grating (FBG) that initiates a series of signal sampling actions clocked by an auxiliary Mach-Zehnder interferometer (MZI) signal. We found that significant instabilities were brought by the relative delays through incoherent timing corrections and timing collisions between the timing references. The best-to-worst ratio of the measured phase errors was higher than 200 while only the signal delays varied. Noise-limited phase stability was achieved with a wide dynamic range of OCT signals above 50 dB in optimized delays. This demonstrated that delay optimization is very effective in phase stabilization of SS-OCT.

Cerebral capillary flow imaging by wavelength-division-multiplexing swept-source optical Doppler tomography

Swept-source-based optical coherence tomography (SS-OCT) has demonstrated the unique advantages for fast imaging rate and long imaging distance; however, limited axial resolution and complex phase noises restrict swept-source optical Doppler tomography (SS-ODT) for quantitative capillary blood flow imaging in the deep cortices. Here, the wavelength-dividing-multiplexing optical Doppler tomography (WDM-ODT) method that divides a single interferogram into multiple phase-correlated interferograms is proposed to effectively enhance the sensitivity for cerebral capillary flow imaging. Both flow phantom and in vivo mouse brain imaging studies show that WDM-ODT is able to significantly suppress background phase noise and image cerebral capillary flow down to the vessel size of 5.6 μm. Comparison between the wavelength-division-multiplexing SS-ODT and the spectral-domain ultrahigh-resolution ODT (uODT) reveals that SS-ODT outpaces uODT by extending the capillary flow imaging depth to 1.6 mm in mouse cortex. Thus, for the first time, quantitative capillary flow imaging is demonstrated using SS-ODT in the deep cortex.

Multicontrast endomyocardial imaging by single-channel high-resolution cross-polarization optical coherence tomography

A single-channel high-resolution cross-polarization (CP) optical coherence tomography (OCT) system is presented for multicontrast imaging of human myocardium in one-shot measurement. The intensity and functional contrasts, including the ratio between the cross- and co-polarization channels as well as the cumulative retardation, are reconstructed from the CP-OCT readout. By comparing the CP-OCT results with histological analysis, it is shown that the system can successfully delineate microstructures in the myocardium and differentiate the fibrotic myocardium from normal or ablated myocardium based on the functional contrasts provided by the CP-OCT system. The feasibility of using A-line profiles from the 2 orthogonal polarization channels to identify fibrotic myocardium, normal myocardium and ablated lesion is also discussed.

Highly phase-stable 200 kHz swept-source optical coherence tomography based on KTN electro-optic deflector

The rapid advance in swept-source optical coherence tomography (SS-OCT) technology has enabled exciting new applications in elastography, angiography, and vibrometry, where both high temporal resolution and phase stability are highly sought-after. In this paper, we present a 200 kHz SS-OCT system centered at 1321 nm by using an electro-optically tuned swept source. The proposed system's performance was fully characterized, and it possesses superior phase stability (0.0012% scanning variability and <1 ns timing jitter) that is promising for many phase-sensitive imaging applications. Biological experiments were demonstrated within ex vivo human tracheobronchial ciliated epithelium where both the ciliary motion and ciliary beat frequency were successfully extracted.