1 November 2005 Heterodyne swept-source optical coherence tomography for complete complex conjugate ambiguity removal
Author Affiliations +
Abstract
Fourier domain (FD) techniques have increasingly gained attention in optical coherence tomography (OCT). This is primarily due to their demonstrated sensitivity of two to three orders of magnitude over conventional time-domain techniques. FDOCT images are subject to two primary sources of artifacts. First, a complex conjugate ambiguity arises because the Fourier transform of the real-valued spectral interferometric signal is Hermitian symmetric. This ambiguity leads to artifactual superposition of reflectors at positive and negative pathlength differences between the sample and reference reflectors. Second, noninterferometric and sample autocorrelation terms appear at dc, obscuring reflectors at zero pathlength difference. We show that heterodyne detection in swept-source OCT (SSOCT) enables the resolution of complex conjugate ambiguity and the removal of noninterferometric and autocorrelation artifacts. We also show that complex conjugate ambiguity resolution via frequency shifting circumvents falloff induced by finite source linewidth in SSOCT when samples are shifted to large pathlength differences. We describe an efficient heterodyne SSOCT design that enables compensation of power losses from frequency-shifting elements. Last, we demonstrate this technique, coupled with wavenumber triggering and electronic demodulation, for in vivo imaging of the human anterior eye segment.

1.

Introduction

Fourier domain optical coherence tomography (FDOCT) techniques have reshaped the OCT landscape due to demonstrated sensitivities that are two to three orders of magnitude greater than their time domain counterparts.1, 2, 3 The shift to FDOCT, however, is not without costs. FDOCT techniques are limited in imaging depth and suffer from two important sources of artifacts. The first, called “complex conjugate ambiguity,” arises because the Fourier transform of the real-valued spectral interferometric signal is Hermitian symmetric. This results in sample reflectors at a positive displacement, +Δx , with respect to the reference reflector, being superimposed on those at a negative displacement, Δx . The second source of artifact, termed “dc artifact,” originates from the noninterferometric light and autocorrelation from sample reflectors, which transform to Δx=0 , and thereby obscure reflectors positioned at zero pathlength difference. These artifacts can be removed by retrieval of the complex interferometric signal. Additionally, retrieval of the complex interferometric signal effectively doubles the imaging depth by removal of the complex conjugate ambiguity.

Previous techniques for acquiring the complex FD interferometric signal have relied on collecting the in-phase and quadrature ( π2 -shifted) components generated by phase shifting interferometry4 or 3×3 interferometry.5, 6 These techniques are constrained to homodyne insofar as both the reference and sample arm optical fields have the same phase velocity. Homodyne detection is required for spectral domain OCT (SDOCT) systems that employ spectrometers coupled to charge accumulation detectors such as charge-coupled devices (CCDs) and photodiode arrays. In swept-source OCT (SSOCT), however, a current-generating photodiode is employed, which enables the spectral interferometric signal to be encoded with a characteristic heterodyne beat frequency. Phase-shifting interferometry requires acquisition of multiple images and therefore is not instantaneous. 3×3 interferometry is an instantaneous technique for acquisition of the complex interferometric signal. However, N×N fiber couplers are wavelength dependent, making it difficult to accurately reconstruct the full complex interferometric signal.

Two techniques for complex conjugate artifact removal using heterodyne SSOCT have previously been reported. The first described a system utilizing electro-optic phase modulators to introduce a heterodyne beat frequency in the reference arm.7 A drawback to this technique is that electro-optic modulators are polarization dependent and highly dispersive, thus requiring methods for controlling the polarization and compensation of the group-velocity dispersion. The second technique utilized acousto-optic modulators for heterodyne SSOCT.8 This system entails an interferometric topology in which the frequency-shifting components are placed between the sample and the receiver. Frequency-shifting elements, such as electro- and acousto-optic modulators, are typically lossy. Placing these elements after the sample prevents the ability to compensate power losses due to light exposure limitations on in vivo samples. Additionally, both published techniques required oversampling of the instrinsic signal bandwidth due to the heterodyne modulation, and the latter technique required a computationally intensive algorithm for resampling the acquired data into frequency space.

In this paper, we present a heterodyne SSOCT system that enables efficient detection of the complex interferometric signal, thus having the potential for increased signal-to-noise ratio (SNR). Additionally, we use wavenumber triggering,9, 10 thereby eliminating the necessity to resample the acquired data in software. We show that heterodyne detection in SSOCT enables complete resolution of complex conjugate ambiguity and the removal of noninterferometric signals. We also show that frequency shifting the reference arm optical field upshifts the cross-interferometric signal to a user-tunable frequency that corresponds to an electronic pathlength mismatch between the interferometer arms. This electronic pathlength mismatching recenters the A-scan at an offset that can be far from dc, which effectively resolves the complex conjugate ambiguity problem. While it is also possible, in principle, to upshift the cross-interferometric signal by placing the sample focus at a large pathlength mismatch, this physical pathlength mismatching leads to significant signal attenuation due to amplitude falloff secondary to finite source linewidths.4, 5 We show that frequency shifting provides a method for complex conjugate ambiguity resolution that circumvents signal falloff that occurs by placing samples at a large pathlength mismatch. Through electronic demodulation, we gain access to the in-phase and quadrature components of the interferometric signal and enables wavenumber triggering, which eliminates the necessity for oversampling.

2.

Theoretical Analysis

2.1.

Imaging Depth

The imaging depth achieved using FDOCT systems is limited by two mechanisms, the spectral sampling interval (which limits the maximum depth observable) and the system spectral resolution (which leads to a falloff of SNR with depth). The maximum imaging depth in FDOCT systems is described as2

1

Δzmax=14δsk,
where δsk is defined as the spectral sampling interval of the FDOCT system. In SDOCT systems, the spectral sampling interval is limited by the pixel spacing of the CCD. Reported maximum imaging depths for 1.3μm SDOCT systems are of the order of 2.0mm .11 For SSOCT systems, the spectral sampling interval is limited by the sampling rate of the temporally sweeping source frequency. Over 4.0mm maximum imaging depths have been achieved for 1.3μm SSOCT systems.12

The second parameter that limits the imaging depth of FDOCT systems is falloff. The sensitivity of FDOCT systems degrades as a function of imaging depth due to fringe washout. The 3dB falloff depth, derived from the analysis reported by Yun 11 is given by

2

Δz3dB=2ln2δrk,
where Δz3dB is the imaging depth at which the SNR is reduced by half, and δrk is the spectral resolution of the FDOCT system. In SDOCT, the spectral resolution is limited by spectrometer optics and/or the pixel width of the CCD. The preceding expression for the 3dB falloff point was derived assuming the spectral resolution of SDOCT systems is limited by the Gaussian beam profile in the spectrometer as opposed to the width of the CCD pixel. For SSOCT systems, the spectral resolution is defined by the instantaneous linewidth of the swept laser source. Since spectral sampling and spectral resolution are coupled in spectrometer-based SDOCT systems, they are more limited by falloff compared to SSOCT techniques. The 3dB imaging depth of a 1.3μm SDOCT system was reported to be 1.6mm .11 In comparison, the 3dB imaging depth for 1.3μm SSOCT is 3.7mm , more than double the distance achieved using SDOCT. Imaging depth capabilities over 4mm can be valuable for applications such as endoscopy, small animal imaging, and human anterior segment imaging.

2.2.

Heterodyne SSOCT

In FDOCT, the photocurrent signal generated by n reflectors is given by

3

i(k)S(k){RR+nRn+2RRnRncos[2k(zRzn)]+2nmn(RnRm)12cos[2k(znzm)]}.
Here, i(k) is the detector photocurrent as a function of optical wavenumber k ; S(k) is the source power spectral density; RR and Rn are the reflectivities of the reference and nth sample reflector, respectively; and zR and zn are the positions of the reference and the nth sample reflector, respectively. The first two terms represent noninterferometric spectral artifacts. The third term represents the cross-interferometric terms, and the fourth term represents the autocorrelation artifact.

In SSOCT systems, the wavenumber is parameterized in time t by the relationship k=k0+t(dkdt) , where k0 is the starting wavenumber, and dkdt is the nonlinear sweep velocity. This sweeping leads to the conversion of pathlength differences in the auto- and cross-terms to an electronic frequency in i(t) , the time-varying photocurrent. The cross-frequencies have instantaneous values of ωn=(dkdt)(zRzn) , while the autofrequencies have instantaneous values of ωnm=(dkdt)(znzm) . If the reference arm optical field is shifted by some beat frequency ωD , then the time-varying photocurrent becomes

4

i(t)S(t){RR+nRn+2RRnRncos[(ωn+ωD)t+ϕn]+2nmn(RnRm)12cos(ωnmt+ϕnm)},
where φn=k0(znzm) and φnm=k0(znzm) . After frequency shifting, the autocorrelation and source spectral terms remain centered at baseband, while the cross-interference terms are recentered around ωD . While the Fourier transform of i(t) remains Hermitian symmetric, the transform of fringes generated by pathlength differences of equal magnitude but opposite sign no longer overlap. This resolves complex conjugate ambiguity because positive displacements are above ωD , while negative displacements are below ωD as long as ωD is larger than the maximum ωn . If the wavenumber sweep is linear over a bandwidth sweep Δk that takes Δt seconds to complete, then ωD corresponds to a pathlength shift of zD=ωDΔt(2Δk) . This shift does not lead to signal falloff as falloff in SSOCT occurs because the interferometric signal is integrated over the source linewidth at the photodiode. If the source linewidth is of the order of 2π(zRzn) , then the linewidth spans an appreciable portion of the interferometric fringe. This decreases the fringe visibility, which decreases the peak height in the Fourier transform of i(t) . Frequency shifting creates a time-varying beat frequency that is independent of sweep speed or source linewidth and, as such, it is not susceptible to falloff.

The cross-interferometric signal can be recovered by bandpass filtering around ωD with a noise equivalent bandwidth of NEB=2zmaxΔkΔt . If demodulation is performed the bandpassed signal is electronically mixed with orthogonal local oscillators with frequency ωD . In this case, the in-phase (real) and quadrature (imaginary) parts of the complex interferometric signal can be recovered:

5.

iRe(t)=2S(t)RRRncos(ωnt+ϕn),
iIm(t)=2S(t)RRRnsin(ωnt+ϕn).
Additionally, after demodulation, the interferometric signals are dependent only on the time-varying frequency ωn , thereby enabling wavenumber triggering.

3.

Experimental Setup

We constructed the heterodyne SSOCT setup shown in Fig. 1 using a fiber-based swept laser source (Micron Optics, Inc., λ0=1310nm , Δλ=100nm , 250 Hz sweep rate), acousto-optic modulators (ACM-1002AA6 IntraAction, Corp.), and a New Focus balanced photoreceiver. The acousto-optic modulators (AOs) had a common center frequency of 100MHz , one having a user-adjustable offset (ωD) from that frequency. The diffraction efficiency of the AOs was measured to be 60%, and using the optical setup depicted in the figure, the maximum diffracted optical bandwidth recoupled into the circulator was 64nm . To trigger each line acquisition, a 250Hz clock was provided by the swept-laser source. Spectral interferogram samples evenly spaced in wavenumber were clocked into the data acquisition system by using a fiber Fabry-Perot interferometer (FFPI, Micron Optics, Inc.). Ten percent of the swept-laser source output was delivered to the FFPI and the output was detected using a New Focus 125 MHz photoreceiver. The detected signal was passed through a 2.1MHz low-pass filter and electronic gain to produce an electronic comb, where each 5V peak was evenly spaced in wavenumber. The wavenumber spacing was 0.1nm , given by the free spectral range of the FFPI.

Fig. 1

Heterodyne SSOCT system utilizing acousto-optic modulators. The swept-laser source has a center wavelength λ0=1310nm and bandwidth Δλ=100nm . A balanced photoreceiver is used for dc artifact suppression. FFPI, fiber Fabry-Pérot interferometer; PC, polarization controllers; AO, acousto-optic modulator; C, circulator; CASS, anterior segment scanner.

064005_1_007506jbo1.jpg

The heterodyne interferometric signal collected using the balanced photoreceiver was high-pass filtered (500kHz) to remove the spectral and autocorrelation artifacts, then demodulated (RF Micro Devices RF2713) by mixing with a local oscillator (LO) of frequency ωD . The in-phase and quadrature components were low-pass filtered ( 500kHz cutoff frequency) and digitized in dual analog-to-digital (A/D) channels (National Instruments PCI 6115) using the clocking signal from the wavenumber trigger. The output power of the swept-laser source was 500μW . There was approximately 6dB source power attenuation in the system prior to the sample ( 3dB through the AO and 2.9dB insertion loss into the fiber of the optical circulator), resulting in 60μW illumination on the sample (cornea).

The SNR of the system, using a 60dB calibrated reflector near zero pathlength difference, was measured to be 99dB . The predicted SNR, assuming same signal power, was 112dB . The heterodyne SSOCT setup described here, places the AO in front of the sample. This design enables compensation of optical power loss from the AO and, therefore, the SNR of the system could be increased by use of a higher power laser source.

4.

Results

To illustrate the behavior of the cross-interferometric term in Eq. 3 with the reference arm frequency shifted with respect to the sample arm, we show in Fig. 2 the fringe patterns of typical interferograms at two pathlength differences, centered around zero, for homodyne [Figs. 2(a) and 2(b)] and heterodyne [Figs. 2(c) and 2(d)] SSOCT setups. The homodyne setup was measured by setting ωD , the frequency difference between the two AOs, to zero. The fringe frequency for the cross-correlation term, when ωD=0 , was identical for equivalent positive [Fig. 2(a)] and negative [Fig. 2(b)] displacements, and it was thus ambiguous whether the pathlength difference was positive or negative. However, when ωD=20kHz , the fringe frequency for the positive displacement [Fig. 2(c)] was higher than for a negative displacement [Fig. 2(d)], as expected from Eq. 4, and therefore the positive and negative locations were resolved.

Fig. 2

Interferograms taken using heterodyne SSOCT: (a) and (b) taken from a reflector at +50μm and 50μm pathlength difference, respectively, with a 0kHz frequency shift, and (c) and (d) collected at the same pathlength difference, but with a 20kHz frequency shift. By using heterodyne detection, reflectors at positive positions produce higher fringe frequencies than those at negative positions and therefore are unambiguously resolved. The chirped signal is a result of wavenumber triggering without demodulation, reinforcing the importance of demodulation when acquiring heterodyne data using a wavenumber trigger.

064005_1_007506jbo2.jpg

Figure 3 demonstrates that the signal falloff remained centered at zero pathlength difference, even as the electronic frequency was shifted. Figure 3(a) shows falloff centered at zero pathlength for homodyne SSOCT. By shifting the zero pathlength difference frequency to ωD=1MHz , the image contents were moved away from dc and bandpass filtered to remove the spectral artifacts. After filtering and demodulation, the falloff was recentered around zero frequency and the 3.5dB signal attenuation out to a 4.0mm scan depth was preserved.

Fig. 3

Falloff measurements for (a) homodyne and (b) heterodyne SSOCTs acquired using setup in Fig. 1; the 3.5dB signal falloff was preserved using heterodyne SSOCT; 60dB calibrated reflector used in sample arm; 1MHz heterodyne frequency shift.

064005_1_007506jbo3.jpg

In Fig. 4 , complete removal of complex conjugate ambiguity and dc artifacts is demonstrated using heterodyne SSOCT. In the A-scan acquired using homodyne SSOCT [Fig. 4(a)], the complex conjugate ambiguity artifact is clearly shown. At zero pathlength difference a strong dc peak from the noninterferometric spectral artifact is also present. By shifting the frequency of the reference arm by ωD=1MHz , the complex conjugate ambiguity was completely resolved [Fig. 4(b)]. Also, the spectral artifact located at dc was removed. There is, however, still present a small dc artifact that originated from the demodulation board. The output signals from the demodulation board contained a small (<0.1V) dc offset. As a result, this dc offset appeared as a small peak located at zero pathlength difference in the A-scans. The dynamic range of the system was limited to 25dB by the input requirements of the demodulation board (<80mV) . As a result, we were able to demonstrate 25dB extinction of the complex conjugate ambiguity artifact. In theory, however, the capability of the system to remove this artifact is not limited.

Fig. 4

A-scans of a 60dB calibrated reflector at a 1.0mm pathlength difference. (a) A-scan showing unresolved complex conjugate ambiguity acquired using homodyne SSOCT. The complex conjugate ambiguity prevents distinguishing reflectors at a +1.0mm from those at 1.0mm pathlength difference. Using heterodyne SSOCT (b), the complex conjugate ambiguity is resolved. The spectral artifact located at dc is also removed, however, a small artifact resulting from a dc offset in the electronic demodulator is present.

064005_1_007506jbo4.jpg

In vivo images acquired using homodyne and heterodyne SSOCT techniques are shown in Figs. 5(a) and 5(b) , respectively. These images are of the anterior segment of a human eye and were constructed using 500 lines/image acquired at a rate of 250Hz . In the homodyne SSOCT image, the iris, cornea, and sclera were obscured by the mirror image arising from the complex conjugate ambiguity artifact. Using heterodyne SSOCT, however, the artifact was completely removed and the cornea, iris, and sclera were no longer obscured. It is evident that at least 6mm of imaging depth is necessary to visualize the entire anterior segment.

Fig. 5

In vivo images of human anterior segment using (a) homodyne and (b) heterodyne SSOCT techniques. Complete complex conjugate ambiguity removal and doubling of the imaging depth is shown to be necessary in order to visualize the entire anterior segment of the eye. C, cornea; S, sclera; I, iris.

064005_1_007506jbo5.jpg

5.

Conclusions

We demonstrated a heterodyne SSOCT technique using AOs for complex conjugate ambiguity and dc artifact removal. We also presented an efficient heterodyne SSOCT design that minimizes power losses after the sample, enabling an increased SNR. This technique enables the zero pathlength difference to be shifted away from electronic dc with no penalty due to finite laser linewidth falloff. Combined with coherent demodulation, this technique provides access to the entire complex interferometric signal without oversampling, as the method also enables data acquisition using wavenumber triggering, thus eliminating computationally intensive resampling algorithms. Although we had limited available laser source power and dynamic range, we presented a technique for heterodyne SSOCT with wavenumber triggering using electronic demodulation. We also demonstrated in vivo human anterior segment images where more than 6mm imaging depth is required.

Acknowledgments

This work was supported by National Institute of Health, R21 EB000243. We greatly appreciate the collaboration with K. Hsu at Micron Optics, Inc.

References

1. 

Leitgeb R., Hitzenberger C. K., Fercher A. F., “Performance of Fourier domain vs. time domain optical coherence tomography,” Opt. Express, 11 (8), 889 –894 (2003). 1094-4087 Google Scholar

2. 

Choma M. A., Sarunic M. V., Yang C., Izatt J. A., “Sensitivity advantage of swept source and Fourier domain optical coherence tomography,” Opt. Express, 11 (18), 2183 –2189 (2003). 1094-4087 Google Scholar

3. 

White B. R., Pierce M. C., Nassif N., Cense B., Park B. H., Tearney G. J., Bouma B. E., Chen T. C., de Boer J. F., “In vivo dynamic human retinal blood flow imaging using ultra-high-speed spectral domain optical Doppler tomography,” Opt. Express, 11 (25), 3490 –3497 (2003). 1094-4087 Google Scholar

4. 

Wojtkowski M., Kowalczyk A., Leitgeb R., Fercher A. F., “Full range complex spectral optical coherence tomography technique in eye imaging,” Opt. Express, 27 (16), 1415 –1417 (2002). 1094-4087 Google Scholar

5. 

Choma M. A., Yang C., Izatt J. A., “Instantaneous quadrature low-coherence interferometry with 3×3 fiber-optic couplers,” Opt. Lett., 28 (22), 2162 –2164 (2003). 0146-9592 Google Scholar

6. 

Sarunic M. V., Choma M. A., Yang C., Izatt J. A., “Instantaneous complex conjugate resolved spectral domain and swept-source OCT using 3×3 fiber couplers,” Opt. Express, 13 (3), 957 –967 (2005). 10.1364/OPEX.13.000957 1094-4087 Google Scholar

7. 

Zhang J., Nelson J. S., Chen Z., “Removal of a mirror image and enhancement of the signal-to-noise ratio in Fourier-domain optical coherence tomography using electro-optic phase modulator,” Opt. Lett., 30 (2), 147 –149 (2005). 10.1364/OL.30.000147 0146-9592 Google Scholar

8. 

Yun S. H., Tearney G. J., Boer J. F. d., Bouma B. E., “Removing the depth-degeneracy in optical frequency domain imaging with frequency shifting,” Opt. Express, 12 (20), 4822 –4828 (2004). 10.1364/OPEX.12.004822 1094-4087 Google Scholar

9. 

Choma M. A., Sarunic M. V., Yang C., Hsu K., Izatt J. A., “Sensitivity advantage of swept source and Fourier-domain optical coherence tomography,” (2004) Google Scholar

10. 

Choma M. A., Hsu K., Izatt J. A., “Swept source optical coherence tomography using an all-fiber 1300-nm ring laser source,” J. Biomed. Opt., 10 044 –009 (2005). 1083-3668 Google Scholar

11. 

Yun S. H., Tearney G. J., Bouma B. E., Park B. H., de Boer J. F., “High-speed spectral domain optical coherence tomography at 1.3μm wavelength,” Opt. Express, 11 (26), 3598 –3604 (2003). 1094-4087 Google Scholar

12. 

Maheshwari A., Choma M. A., Izatt J. A., “Heterodyne swept-source optical coherence tomography for complete complex conjugate ambiguity removal,” Proc. SPIE, 5690 91 –95 (2005). 0277-786X Google Scholar
© (2005) Society of Photo-Optical Instrumentation Engineers (SPIE)
Anjul Maheshwari, Anjul Maheshwari, Michael A. Choma, Michael A. Choma, Joseph A. Izatt, Joseph A. Izatt, "Heterodyne swept-source optical coherence tomography for complete complex conjugate ambiguity removal," Journal of Biomedical Optics 10(6), 064005 (1 November 2005). https://doi.org/10.1117/1.2136147 . Submission:
JOURNAL ARTICLE
6 PAGES


SHARE
Back to Top