Abstract
We propose a new model to characterize the phase noise in swept-source optical coherence tomography (SS-OCT). The new model explicitly incorporates scanning variability, timing jitter, and sample location in addition to intensity noise (shot noise). The model was analyzed and validated by using both Monte Carlo methods and experiments. We suggest that the proposed model can be used as a guideline for future SS-OCT experimental designs.
© 2017 Optical Society of America
There are several mathematical models to theoretically describe and analyze the phase noise of Fourier-domain optical coherence tomography (FD-OCT) systems. Park et al., proposed a phasor picture to understand the composition of the phase noise in spectral-domain OCT (SD-OCT): the noise was modeled as a vector orthogonal to the signal vector with a random magnitude [1]. Choma et al., proposed a similar model [2] and further reframed it as a constant phasor plus additive white Gaussian noise (AWGN) problem: the noise was considered as a phasor with a constant magnitude but with a uniformly distributed random angle [3]. Szkulmowski et al., on the other hand, modeled the same problem as a constant phasor plus a sum of random phasors [4]. Despite the differences among these models, similar conclusions were drawn: the phase stability (sensitivity) of the FD-OCT measurement was found inversely proportional to the square root of the system’s signal-to-noise ratio (SNR) if the system was shot noise limited [1–4]. In addition to the aforementioned studies based on SD-OCT, Vakoc et al., discussed the impact of the timing jitter on the phase performance of swept-source optical coherence tomography (SS-OCT) [5], and they reached a similar conclusion that the phase stability of SS-OCT was also SNR limited.
Currently, phase-resolved OCT or related phase-sensitive techniques relies on the previously mentioned frameworks to predict the system’s phase stability. Although the models generally fit the measurements in SD-OCT systems, their effectiveness was contested in SS-OCT platforms: the measured phase stability in SS-OCT could be orders of magnitude worse than the predictions [6,7]. In this Letter, we propose a modified phase noise model to incorporate two important factors that are unique to SS-OCT’s phase stability: the “scanning variability” and the “timing jitter.” The proposed model suggests that the phase stability of the SS-OCT, which is collectively determined by the newly introduced factors and the optical path length (OPL) between the reference mirror and the sample, is only linear with up to a certain level. Beyond the threshold, the system’s phase stability becomes saturated as the SNR further increases. We verified our model on a custom-built SS-OCT system, and it succeeded in fitting the experimental data with high accuracy.
The detected interference signal of single sample reflector from a typical FD-OCT system can be written as [8]
where is the power spectrum of the light source, is the wavenumber, is the responsivity of the photodetector, is the displacement from the sample reflector to the reference mirror, is the noise presented, and and are the reflectivity of the reference mirror and the sample reflector, respectively. Here we assume the noise is an AWGN, which is valid in a shot-noise-limited scenario [1,2].After being uniformly discretized times in the domain, the detected signal could be written as
where is the index of the spectral sampling point, is the sampling interval in the domain, , and . Without losing generality, we assumed to be a constant number . We then performed an -point inverse discrete Fourier transform on Eq. (2) against as where is a constant since is AWGN. The noise phase angle could be modeled by a uniform distribution over . In OCT, we are only interested in the positive “frequency” portion of Eq. (3). Notice that could be further decomposed into two parts: as the peak amplitude location of , and as the subaxial-sampling-interval displacement. Without considering the noise corruption at this moment, we could rewrite the phase angle of the positive “frequency” components at the peak amplitude location asIf the sample only presents subaxial-sampling-interval movement, the first term of Eq. (4) could be assumed constant during the imaging process. Then, a faithful demodulation of in the second term will rely on the condition that both coefficients, and , be constant. In SD-OCT, where the interferogram is measured by a spectrometer, this condition generally holds.
However, it is not guaranteed in SS-OCT. First, the output spectrum of the wavelength swept laser source could vary from scan to scan, as illustrated in Fig. 1(a), due to the mechanical or electrical scanning. The deviation could be characterized by a random variable , and we call it “scanning variability.” Without losing generality, we assume that all follow the same normal distribution, , with a zero mean and a standard deviation of . This assumption makes both terms, and , in Eq. (4) random variables rather than constants in SS-OCT. We acknowledge that various techniques have been proposed to reduce this effect [9,10]. For example, Gora et al., suggested to simultaneously record a reference signal to recalibrate the sample signal [10]. However, we do not believe this “scanning variability” could be completely eliminated and we keep this parameter for the completeness of the study.
Second, a timing jitter exists in SS-OCT, as shown in Fig. 1(b), due to the imperfect synchronization between the light source and the data acquisition (DAQ) board, the limited sampling rate of the DAQ, and the noise in the trigger signal. Since the timing jitter from the imperfect synchronization could be greatly suppressed by various methods [11,12], we focus on the jitter caused by the limited sampling rate and the trigger noise in this Letter, which inherently stems from the temporal measurement nature of SS-OCT. We then model the timing jitter as another random variable that is linearly proportional to , the sampling interval of the DAQ. Based on our assumption, a slower DAQ will lead to a larger timing jitter and a faster DAQ will result in a smaller timing jitter. Assume also follows a zero-mean normal distribution, and its standard deviation is written as . The standard deviation is determined by a coefficient and the sampling interval of the DAQ. The impact of the introduction of on the measured wavenumber is translated via the scanning speed of the light source . Here, is the A-line scanning period of the light source. Equation (4) could thus be written as
The first term in Eq. (5) corresponds to the structural information that is readily available in OCT; the second term represents the subaxial-sampling-interval phase term that techniques such as phase-resolved OCT intend to decouple; the third term is the phase modulation introduced by both the timing jitter and the scanning variability. The rest of the terms are considered higher order infinitesimals that could be ignored.A schematic of phasor interpretation of Eq. (5) is given in Fig. 2. The phase angle of the original signal is given by the second term of Eq. (5). The inclusion of the third term effectively rotates the phase angle of by . After the rotation, the phasor is further added by the random noise phasor , and the final resultant phasor has a phase angle of .
Since both and are random variables, the measured phase angle is also a random variable. Let us assume that the signal phasor has a constant amplitude 1, the random phasor has a constant amplitude , and a uniformly distributed phase angle . The probability density function of the phase angle of the measured phasor can be written as
whereThe standard deviation of the phase angle essentially represents how precisely an OCT could measure the phase of the sample. It can be calculated as the root mean square error (RMSE) of the multiple phase measurements. We first defined the following -sweep variation coefficient to conveniently characterize the system phase stability:
where is the wavenumber of the central wavelength of the light source. The larger the is, the higher phase stability the system possesses. The SNR of the system isWe then conducted a computer simulation by using the Monte Carlo method. The calculated at the sample location () was plotted in logarithmic scale as a color-coded image against two variables, and SNR, in Fig. 3(a). ranged from 1 to 200 dB and SNR varied from 1 to 100 dB. Both were in a 1 dB step. For each pair, 10 million signal photons were launched. The program was implemented in MATLAB and run on the high-performance computing cluster (Columbia University Yeti HPC cluster, USA). We plotted the RMSE of the measured phase at different SNRs against in Fig. 3(b). When the -sweep variation coefficient, , is very low, the system suffers large noises from the scanning variability and the timing jittering. The measured phase angle is random in the entire range of . In this regime, is bounded and is not changed with either SNR or . When rises above a certain threshold (30 dB), starts to reduce linearly against in logarithmic scale. If is further increased, will saturate at a certain level that is determined by both SNR and . We also plotted at different values of against SNR in Fig. 3(c). When SNR is large, the noise in Eq. (6) approaches 0. The probability distribution function shown in Eq. (6) is equivalent to a normal distribution with a standard deviation of . It is noted that tends to saturate at a higher level with higher . For large (), is mostly linear in the selected dynamic range. In fact, the curve will converge to the previously reported SD-OCT case [3] if goes to infinity.
Experiments were conducted to verify the proposed model. A custom-built SS-OCT system was developed by using an electrically tuned wavelength swept light source (NTT-AT, Japan). The source features a phase-locked loop-connected A-line sweeping and trigger output [13]. The central wavelength of the light source is 1329 nm, and the wavelength span is 101 nm. The A-line rate of the system is 200 kHz with a duty cycle of 50%. A dual-balance detection scheme was realized by using a broadband photoreceiver (PDB480C-AC, Thorlabs, U.S.), and a 12 bit high-speed data acquisition board (ATS9373-A3, AlazarTech, Canada) was used to digitize the analog interference signals. Samples were placed in the focal plane of the sample arm with the reference arm blocked to achieve a common-path interferometer configuration. A reference clock was used to resample the sample signal to minimize the scanning variability [10].
Since the scanning variability is inherent to a specific system, we first verified the effect of timing jitter with controlled . A piece of No. 1 glass coverslip was used, and we altered the by changing the sampling rate of the DAQ from 200 MS/s to 1.8 GS/s. For each sampling rate, we tuned the SNR of the OCT signal by a continuously variable neutral density filter (NDC-25C-4M, Thorlabs, U.S.). The results are plotted in Fig. 4(a). As a comparison, we also added the experimental data obtained from a 32 kHz SD-OCT system using a superluminescent diode. The SD-OCT measurements were very close to the theoretical limit as predicted in Ref. [3]. However, the SS-OCT measurements saturated at higher SNR; this behavior is predicted in our model, as shown in Fig. 3(c).
A family of curves was obtained based on Eq. (6) via Monte Carlo simulation and then fitted to the experimental data by using least squares approximation. The fitted are listed in Table 1.
After that, we used the fitted to obtain the coefficient and in Eq. (7). From Fig. 4(b), we can see that the fitted coefficient is equal to 0.5612, and the fitted is . We then fixed the sampling rate of the DAQ to be 800 MS/s and validated our model with different values of . We used No. 1 (110 μm), No. 1.5 (130 μm), and No. 2 (170 μm) glass coverslips as samples in a common-path configuration. Theoretical predictions were made by using the fitted parameters and . The experimental data and the theoretical predictions presented great agreements, as shown in Fig. 4(c). In summary, the fitted timing jitter of the system at 800 MS/s is 701.5 ps, while the fitted scanning variability is 0.0019% compared with .
Finally, based on the obtained system parameters ( and ), we plotted the achievable phase stability of our SS-OCT as a function of the sampling rate and sample location in Fig. 5. This plot could be used as a guideline for future SS-OCT experiment design; with required phase stability, available DAQ board, imaging depth, and data bandwidth in mind, the researchers could find a niche to achieve best performance.
In conclusion, we presented a new theoretical model for phase noise analysis of SS-OCT. Compared with the models previously conceived upon SD-OCT, the proposed model incorporates two new factors and fits the SS-OCT experimental data better. The new model also suggests that SS-OCT is fundamentally at a disadvantage in absolute phase measurements against SD-OCT in terms of phase stability/sensitivity due to its temporal measurement nature; the collection of the three factors affects as a multiplicative noise on the original signal in addition to the additive noise such as the shot noise and the excess noise. The sole increment on the SNR will no longer boost the phase stability in SS-OCT. However, there are merits for using SS-OCT in phase-related tasks. The faster A-line acquisition of SS-OCT could improve the temporal resolution of sectional or volumetric phase measurements. The variance on the absolute phase measurement could also be washed out if the absolute measurement was later Fourier transformed to obtain the vibrational information.
Funding
National Institutes of Health (NIH) (1DP2HL127776-01).
REFERENCES
1. B. H. Park, M. C. Pierce, B. Cense, S.-H. Yun, M. Mujat, G. J. Tearney, B. E. Bouma, and J. F. de Boer, Opt. Express 13, 3931 (2005). [CrossRef]
2. M. A. Choma, A. K. Ellerbee, C. Yang, T. L. Creazzo, and J. A. Izatt, Opt. Lett. 30, 1162 (2005). [CrossRef]
3. M. A. Choma, A. K. Ellerbee, S. Yazdanfar, and J. A. Izatt, J. Biomed. Opt. 11, 024014 (2006). [CrossRef]
4. M. Szkulmowski, A. Szkulmowska, T. Bajraszewski, A. Kowalczyk, and M. Wojtkowski, Opt. Express 16, 6008 (2008). [CrossRef]
5. B. J. Vakoc, S. H. Yun, J. F. de Boer, G. J. Tearney, and B. E. Bouma, Opt. Express 13, 5483 (2005). [CrossRef]
6. D. C. Adler, R. Huber, and J. G. Fujimoto, Opt. Lett. 32, 626 (2007). [CrossRef]
7. M. Bonesi, M. P. Minneman, J. Ensher, B. Zabihian, H. Sattmann, P. Boschert, E. Hoover, R. A. Leitgeb, M. Crawford, and W. Drexler, Opt. Express 22, 2632 (2014). [CrossRef]
8. W. Drexler and J. G. Fujimoto, eds., Theory of Optical Coherence Tomography (Springer, 2015), pp. 65–94.
9. D. C. Adler, Y. Chen, R. Huber, J. Schmitt, J. Connolly, and J. G. Fujimoto, Nat. Photonics 1, 709 (2007). [CrossRef]
10. M. Gora, K. Karnowski, M. Szkulmowski, B. J. Kaluzny, R. Huber, A. Kowalczyk, and M. Wojtkowski, Opt. Express 17, 14880 (2009). [CrossRef]
11. W. Choi, B. Potsaid, V. Jayaraman, B. Baumann, I. Grulkowski, J. J. Liu, C. D. Lu, A. E. Cable, D. Huang, J. S. Duker, and J. G. Fujimoto, Opt. Lett. 38, 338 (2013). [CrossRef]
12. G. Liu, O. Tan, S. S. Gao, A. D. Pechauer, B. Lee, C. D. Lu, J. G. Fujimoto, and D. Huang, Opt. Express 23, 9824 (2015). [CrossRef]
13. Y. Sasaki, M. Fujimoto, S. Yagi, S. Yamagishi, S. Toyoda, and J. Kobayashi, Proc. SPIE 8934, 89342Y (2014). [CrossRef]