We report a Doppler optical cardiogram gating technique for increasing the effective frame rate of Doppler optical coherence tomography (DOCT) when imaging periodic motion as found in the cardiovascular system of embryos. This was accomplished with a Thorlabs swept-source DOCT system that simultaneously acquired and displayed structural and Doppler images at 12 frames per second (fps). The gating technique allowed for ultra-high speed visualization of the blood flow pattern in the developing hearts of African clawed frog embryos (Xenopus laevis) at up to 1000 fps. In addition, four-dimensional (three spatial dimensions + temporal) Doppler imaging at 45 fps was demonstrated using this gating technique, producing detailed visualization of the complex cardiac motion and hemodynamics in a beating heart.
© 2007 Optical Society of America
Optical coherence tomography (OCT) is a high resolution (1–10 μm) imaging technique capable of visualizing structures to a depth of 2-3mm in soft tissues. The high spatial resolution provided by OCT has been used to study the developing cardiovascular system of various models in developmental biology including Drosophila (fruit fly) , Xenopus laevis (African clawed frog tadpole) [2–5], chick [6, 7] and mouse embryos [7, 8]. Doppler OCT (DOCT) is a functional extension of OCT, analogous to Doppler ultrasound, that can provide a flow sensitivity down to 10-100 μm/s [9–12]. DOCT imaging of Xenopus embryos has been demonstrated at up to 16 frames per second (fps) . In rodent models the heart rate can range from 200 to 500 beats per minute, which makes high temporal resolution a necessity for accurate characterization of cardiovascular function and development. To date, the highest in vivo structural OCT frame rate has been 360 fps . High frame rate color Doppler flow imaging in OCT is very challenging due to the increased Doppler noise floor that accompanies the higher scanning optics velocities . In Doppler ultrasound, this effect is negated through the use of arrayed transducers and electronic scanning such that no physical scanning motion is required. The concept of fiber array OCT systems has previously been described , but its use in DOCT has yet to be demonstrated.
The first demonstration of artificially increasing the imaging frame rate for OCT was by Yandzafar et al, who used a retrospective gating technique on the structural signal  to generate two-dimensional structural and color Doppler images from a time domain OCT system with a slow delay line frequency of 8 Hz . Jenkins et al. have used a prospective gating technique based on electrical pacing of excised hearts of chick and mouse embryo models to allow ex vivo 4D structural imaging at 16 fps . Most recently, in the related field of high frequency ultrasound (~50 μm resolution), Cherin et al. have performed in vivo structural and colour Doppler flow imaging of the mouse heart at 1000 fps using a retrospective gating technique based on a simultaneously acquired electrocardiogram (ECG) .
The retrospective gating techniques demonstrated by both Yandzafar and Cherin have the potential to increase the frame rate for color Doppler flow imaging while maintaining an acceptable noise floor. In this paper, we demonstrate a retrospective gating technique based on a Doppler optical cardiogram as it proved difficult to obtain an ECG signal from the Xenopus. In addition, for in utero mouse embryos, obtaining ECG signal would be problematic due to the difficulty in electrode placement and interference from the maternal heart or surrounding sibling’s embryonic hearts. We demonstrate 1 kHz frame rates for both structural and Doppler images of the Xenopus cardiac dynamics as well as 4D (3D spatial + temporal) imaging at 45 fps. This was accomplished on a swept source (SS) OCT system, while using a separate Time Domain (TD) DOCT to provide the Doppler optical cardiogram.
2. Methods and materials
2.1 SS-DOCT system hardware
The Swept Source Doppler OCT (SS-DOCT) system is based on a commercially available frequency swept source (Thorlabs, SL1325-P16) that is centered at 1325 nm, with a FWHM of the optical spectrum of >100 nm, and an average output power of 12 mW. Further details about the swept source can be found in Ref. . The OCT interferometer is based on a 90/10 Mach-Zehnder interferometer (MZI) where 90% of the power goes to the sample arm and 10% to the reference arm. The interference fringe signal is measured by a 15 MHz balanced photodetector (Thorlabs, PDB140C). Both the MZI clock signal and the interference fringe signal are digitized at 50 MS/s with 14-bit A/D resolution. The interference fringe data is recalibrated by using the MZI clock signal as described in Ref.  yielding ~1024 points per laser sweep. The SS-OCT system also outputs two waveforms, through a 16 bit digital to analog card (National Instruments, PCI-6731), to drive the XY galvanometer and scanning mirror assembly (Cambridge technology, 6210). The TD-OCT system is triggered by the start of the scanning driver waveform from the SS-DOCT system and streams the optical cardiogram data continuously to the hard drive. The SS-DOCT system acquires time-stamped B-mode structural images at 25 fps and streams the raw data to disk at the same time. Real time Doppler OCT imaging is performed at 12 fps using unidirectional wavelength sweeps from the laser. This imaging rate allowed for real-time feedback, which was required for tadpole positioning and visualization of breathing motion artifacts. Each B-mode image consisted of 512 A-scans per image. Prior to in vivo Xenopus cardiac imaging, the phase stability of the SS-OCT system was measured.
2.2 Doppler optical cardiogram hardware
A schematic of the imaging setup is shown in Fig. 1. The TD-OCT system has been previously described with a Michelson interferometer and a rapid scanning optical delay (RSOD) line at a 12.95 kHz resonant frequency . The light source (Inphenix, IPSDM-1322) is centered at 1310 nm with ~55nm full width half maximum (FWHM) bandwidth and an output power of 22 mW. The interference signal was measured with a 10 MHz dual balanced photodetector (New Focus, 2117-FC). A 4.3 MHz electro-optic phase modulator is placed in the delay line providing stability for phase measurements while also providing the carrier frequency for in-phase (I) and quadrature (Q) demodulation. The I and Q signals were digitized with a 10 MHz data acquisition card (DAQ) (National Instruments, PCI-6115). The optical cardiogram was obtained via a previously described GRIN fiber imaging probe [18, 19].
2.3 Structural and Doppler signal processing
The results presented in this paper used only the forward wavelength (short to long) sweeps of the laser. Similar results can be produced from analyzing the backward wavelength sweeps of the laser. The recalibrated interference fringe pattern F(k) can be processed to obtain both the structural and Doppler information encoded in the raw data. The Fourier transform of F(k) gives the complex signal I(y)+iQ(y). The structural image can be obtained by calculating the magnitude of the complex signal. A sliding window spatial averaging mask for the structural images defined in Eq. (1), is used to improve signal-to-noise (SNR) as has been described previously . Briefly, M defines the size of the mask in the axial direction while N defines the size of the spatial averaging mask in the transverse direction. The B-mode structural images shown in this paper were generated using N=4 and M=2.
The frequency shift induced by moving scatterers in the direction of the incident laser beam can be estimated using the Kasai autocorrelation function . The Kasai autocorrelation function measures phase shifts between two adjacent A-scans and is shown in Eq. (2).
The spatial averaging mask values used for all B-mode Doppler images acquired by the SS-DOCT system were N=8 and M=2.
The phase shift induced by moving scatterers is displayed using a standard Doppler colormap with red and yellow indicating flow in the direction of the optical beam while blue and turquoise coloring indicate flow in the opposite direction in unaliased images. The maximum detectable phase shift without aliasing is ±π, due to the complicated blood flow profiles and the presence of noise no attempt is made to phase unwrap images. Furthermore, due to the difficulty in determining the Doppler angle in the complicated structures of the heart we do not relate the phase shift to the velocity of moving scatterers.
2.4 Doppler optical cardiogram and retrospective gating
The RSOD frequency was measured at the beginning of each imaging session, allowing for accurate calculation of the elapsed time during data acquisition by counting the number of Ascans acquired. An example of the M-mode structural image through one of the two great vessels leaving the tadpole heart is shown in Fig. 2(a). The corresponding M-mode Doppler image is shown in Fig 2(b), the size of the averaging mask used to generate this image was N=13 and M=1. The pulsatile nature of the blood flow is clearly discernible in this image. Averaging the phase shift across the vessel shown in Fig 2(b) generates the cardiogram signal shown in Fig. 2(c). A low pass finite impulse response filter was used to remove the high frequency noise with the result demonstrated in Fig. 2(d). One cardiac cycle was defined by the region between the vertical dotted lines. This allowed for the calculation of the duration in each cardiac cycle. The mean and standard deviation of the Xenopus cardiac cycle duration during the entire acquisition were calculated, and cardiac cycles with variations greater than 1.75 standard deviations were rejected.
The overall system timing accuracy was 1 ms, and thus the highest frame rate achievable was 1000 fps for 2D movies. To achieve this maximum frame rate each cardiac cycle within the 1.75 standard deviation was divided into a number of equally spaced temporal bins determined by the mean cardiac cycle length over the image acquisition time of 2-3 minutes. For 1000 fps movies shown in section 3.2 the mean cardiac cycle length was approximately 800 ms and thus each cardiac cycle was divided into 800 temporal bins. B-Mode frames acquired by the SS-DOCT system were placed into the correct temporal bin using the time stamp attached to each frame. If multiple frames were found to originate from the same bin only the first frame was kept. The frames were then re-ordered according to their temporal bin location to achieve the higher effective frame rate.
To reconstruct the 4D structural and Doppler models, the tadpole heart was imaged at approximately 150 different spatial locations. Each position was imaged for 15 seconds after which the 2D scanner stepped 10 μm and the process repeated. During each 15 second acquisition the cardiogram signal was also simultaneously collected with the TD-OCT system. An identical process to that described above was also followed for increasing the frame rate from 12 to 45 fps, using 36 temporal bins. The use of 36 temporal bin here was chosen as a trade off between good temporal resolution and performing manageable 4D reconstructions. The heart was then segmented in each of 36 different bins using the structural images. Masks were created to generate the surface reconstruction at each of the 36 bins as well as to segment the corresponding Doppler images.
2.5 Animal protocols
Xenopus laevis embryos were provided by Dr. Winklbauer in the Department of Zoology at the University of Toronto and were housed in standard conditions under an animal protocol approved by the institutional animal care committee at Princess Margaret Hospital, Toronto, Canada. Stage 51 tadpoles were anesthetized with 21°C (room temperature) Lidocaine solution diluted to 0.001% with tap water treated with Aqua Plus water conditioner. Xenopus laevis were allowed to sit in this solution for approximately 10 minutes before being placed in a shallow V-groove with a small amount of the solution. Images where then taken for approximately 1-2 hours and finally the animals were euthanized by anesthetic overdose.
In previous studies we have used 0.01% Tricaine as our anesthetic of choice. We chose Lidocaine in this study as a lower cost alternative. We compared its performance to Tricaine and found no significant degree of arrythmia under our experimental conditions since over an imaging time of 180 secs, the cardiac cycle variation was only 2%. This may be attributed to the species difference between mammalian and amphibian cardiac systems.
The stage 51 Xenopus heart is a three-chambered heart consisting of two atria and one ventricle. Unlike the four-chambered mammalian heart the amphibian heart mixes oxygyenated (left atrium) and deoxygenated (right atrium) blood inside the single ventricle. This mixture is pumped out of the ventricle to the systemic capillaries and the lungs through the common outflow tract known as the truncus arteriosis. The oxygenated blood from the lungs then returns to the left atrium while the deoxygenated blood from the systemic capillaries returns to the right atrium.
3.1 Phase stability results
Prior to Xenopus imaging it was necessary to characterize the phase noise properties of the SS-DOCT system. The results are summarized in Fig. 3 and indicate a factor of five improvement (for N=2) in the intrinsic phase noise as compared to our TD-OCT system , as shown in Fig. 3(a). However, the overall Doppler noise in the system was still dominated by speckle modulation and sample arm scanning noise, which increased as a function of frame rate as shown in Figs. 3(b) and 3(c). The TD-OCT system used in this study, exhibits phase stabilities as a function of frame rate, that are similar to those found in Ref. .
3.2 Color doppler visualization at 1000 fps
An example of the original real-time 12 fps movie acquired by the SS-DOCT system is shown in Fig. 4. Structural intensity based thresholding of the Doppler images were performed in real-time to minimize Doppler noise, resulting in the movie in Fig. 4(b). This provided necessary feed back for proper positioning of the sample, identification of blood flow, and monitoring the tadpole for breathing motion during data acquisition to generate higher effective frame rate as described in Section 2.2.
The resulting 1000 fps movies are shown in Fig. 5. Gated imaging at effectively 1000 fps gave us the option to implement a sliding window temporal averaging scheme over several frames to improve the SNR in both the structural and Doppler images. Such a high effective frame rate allowed detailed capture of the complex cardiac motion and hemodynamics during the heart cycle. Even more importantly, this technique will allow for proper imaging of the mouse cardiovascular system, where heartbeat can exceed 6 beats per second.
Furthermore, the left and right arches traverse through different angles (which changes during the cardiac cycle); these angles represent different incident Doppler angles to the optical axis. In addition, since the branching point moves into and out of the 2D imaging plane, with the downstream branches (especially of the right aortic arch) pointing in different directions, positive and negative Doppler frequencies are observed due to this diverging flow.
The difference in flow pattern between the right and left aortic arches further exemplifies the need of 4D imaging to fully visualize and comprehend the complex nature of cardiac hemodynamics, even in a simple amphibian model.
The 1000 fps B-mode structural and Doppler images can be reconstructed to provide details in different planes. An example is shown in Fig. 6, where the data is presented in the x-t plane for better visualization of the temporal flow pattern, and appreciation of the different phases of the cardiac cycle.
3.3 Three dimensional + temporal (4D) imaging of the tadpole heart
Starting from the head end of the tadpole heart, we acquired SS-DOCT and Doppler optical cardiogram data as described in section 2.2 at different slice positions towards the tail. Three examples of the 45 fps structural and Doppler movies generated are shown in Figs. 7–9.
The entire 4D data set can then be used for other visualization reconstructions to better illustrate the complex cardiac motion during a heart beat cycle. Figure 10 shows a surface reconstruction of the tadpole heart, rotating while beating, to demonstrate the dynamic anatomical relationships. For example, the truncus arteriosis (TA) oscillates over a significant distance relative to the right and left atria (RA and LA), which are situated more posterior. During each cardiac cycle, the TA not only changes position, but its shape and diameter also vary periodically to accommodate the output blood flow from the ventricle (V).
Internal to the heart, the blood flow pattern and the heart wall motion velocity can also be visualized at different positions of the same cardiac phase (Fig. 11) or at the same position during different phases (Fig. 12).
The demonstration above reports the highest effective frame rate in OCT and DOCT imaging of embryonic heart samples in which the motion is periodic over 2 to 3 minutes of time. It is also the first demonstration of combined 4D structural and blood flow OCT imaging, where the full complexity of the embryonic heart motion and blood flow can be visualized in detail. The Doppler optical cardiogram provides a novel method of tracking the highly periodic flow and motion in the cardiovascular system. In situations where conventional ECG signals are difficult to obtain, such as mouse embryos in utero, it provides a relatively easy way to accomplish gated 2D, 3D and 4D reconstructions. The techniques presented are suitable for in utero high speed 4D imaging and can be extended to visualize mouse embryo hearts at frame rates over 300 fps with low Doppler noise background.
Using a second OCT probe to provide the Doppler optical cardiogram also offers the possibility of imaging the peripheral vascular system, away from the heart where periodical structural motion may not be prominent. Hence, in organs such as the brain or kidney, single channel OCT with structural gating technique  will not be sufficient for 4D imaging. Acquisition of the gating signal on the same SS-DOCT system is an attractive option and thus we are currently in the process of modifying existing systems to perform this function. Since most commercially available DAQ cards have multiple analog input channels, systems can be modified to capture both the cardiogram and the B-mode images using a two-interferometer setup. The use of booster semi-conductor optical amplifiers may be necessary if the laser output power is not sufficient.
During data acquisition there are a number of effects that can lead to degradation in the temporal resolution of the reconstructions. This degradation effect, appears as “jittering”, tends to increase as the duration of data acquisition increases. We attribute this degradation to: (a) environmental effects such as evaporation of water in V-groove can cause a drift of the specimen; (b) small variations of heart beat period due to anesthetics; (c) breathing motion of the tadpole; and (d) de-correlation between the two DOCT systems. The de-correlation is limited by the accuracy of the two computer system clocks and fluctuations in the RSOD and SS sweeping frequencies. These effects limit the extent to which frame rates can be increased without introducing a significant amount of jitter. We found that the best compromise between jitter and high frame rates occurred at approximately 200-300 fps, corresponding to an acquisition time of approximately 20-30 seconds with the current setup. Using a two-channel setup, one could limit the de-correlation effect and allow higher effective frame rates in the future.
The Doppler noise floor in the 12 fps real-time acquisition and display mode of the SS-DOCT used in this study can be improved by using a swept wavelength source with higher sweeping frequency to increase the A-scan density (i.e., A-scans/mm) for a given frame rate. In addition, frequency domain mode locking  will also improve the phase stability and subsequently the accuracy of Doppler shift estimation.
The spatial-temporal artefacts observed in the high frame rate movies can be removed by re-sampling the data in the direction parallel to the equal-temporal lines. An alternate approach is to stop the lateral scanning and simply acquire M-mode data in a point-by-point fashion in the x-y plane. The advantage of such an approach is that the Doppler noise floor will be reduced to the intrinsic noise level of the SS-DOCT without contribution from the lateral scanning optics. The trade-off would be the loss of real-time monitoring capability of any breathing motion or other artefacts, and further post-processing detection of such effects would be required.
In conclusion, we have demonstrated a Doppler optical cardiogram gating technique to increase the effective frame rate for DOCT systems. This method allowed detailed Doppler imaging of the blood flow in the aortic arches of Xenopus laevis embryos at 1000 fps with 512 lines per frame on a SS-DOCT system. This image gating technique also yielded 4D imaging visualization of the complex 3D motion of the embryonic heart, and the hemodynamics therein, at 45 fps. These techniques may be used to image the cardiovascular system development in wild type and transgenic animal models, where phenotypic determination of cardiac morphology and function are important for studying the underlying genetics.
This work was supported by the Canadian Institutes of Health Research, Photonics Research Ontario, the Premier’s Research Excellence Award, the Gordon Lang and the Samuel B. McLaughlin Foundations, and the Canadian Cancer Society through a grant from the National Cancer Institute of Canada. The authors acknowledge assistance from Dr. Rudolf Winklbauer from the Department of Zoology, University of Toronto, ON, Canada, and Jurek Smolen and Xiao Jin from the Ontario Cancer Institute, Toronto, ON, Canada.
References and links
1. M. A. Choma, S. D. Izatt, R. J. Wessells, R. Bodmer, and J. A. Izatt, “Images in cardiovascular medicine: in vivo imaging of the adult Drosophila melanogaster heart with real-time optical coherence tomography,” Circulation. 114,e35–e36 (2006). [CrossRef] [PubMed]
2. S. A. Boppart, G. J. Tearney, B. E. Bouma, J. F. Southern, M. E. Brenzinski, and J. G. Fujimoto, “Noninvasive assessment of the developing Xenopus cardiovascular system using optical coherence tomography,” Proc. Natl. Acad. Sci. 94,4256–4261 (1997). [CrossRef] [PubMed]
4. V. XD. Yang, M. Gordon, E. Seng-Yue, S. Lo, B. Qi, J. Pekar, A. Mok, B. Wilson, and I. Vitkin, “High speed, wide velocity dynamic range Doppler optical coherence tomography (Part II): Imaging in vivo cardiac dynamics of Xenopus laevis,” Opt. Express 11,1650–1658 (2003). [CrossRef] [PubMed]
5. S. A. Boppart, G. J. Tearney, B. E. Bouma, J. F. Southern, M. E. Berzinski, and J. G. Fujimoto, “Noninvasive assessment of the developing Xenopus cardiovascular system using optical coherence tomography,” Proc. Natl. Acad. Sci. USA. 94,4256–4261 (1997). [CrossRef] [PubMed]
6. T. Mesud Yelbuz, M. A. Choma, L. Thrane, M. L. Kirby, and J. A. Izatt, “A new high-resolution imaging technology to study cardiac development in chick embryos,” Circulation. 106,2771 (2002). [CrossRef] [PubMed]
7. M. W. Jenkins, F. Rothenberg, D. Roy, V. P. Nikolski, Z. Hu, M. Watanabe, D. L. Wilson, I. R. Efimov, and A. M. Rollins, “4D embryonic cardiography using gated optical coherence tomography,” Opt. Express 14,736”748 (2006). [CrossRef] [PubMed]
8. D. L. Marks, T. S. Ralston, and S. A. Boppart, “Three-dimensional optical coherence tomography of the embryonic murine cardiovascular system,” J. Biomed. Opt. Special Issue on Cardiovascular Photonics ,11,021014 (2006).
9. S. Yazdanfar, A. M. Rollins, and J. A. Izatt, “Imaging and velocimetry of the human retinal circulation with color Doppler optical coherence tomography,” Opt. Lett. 251448–1450 (2000). [CrossRef]
11. A. W. Schaefer, J. J. Reynolds, D. L. Marks, and S. A. Boppart, “Real-time digital signal processing-based optical coherence tomography and Doppler optical coherence tomography,” IEEE Trans. Biomed. Eng. 51186–190, (2004). [CrossRef] [PubMed]
13. R. Huber, D. C. Adler, and J. G. Fujimoto, “Buffered Fourier domain mode locking: unidirectional swept laser sources for optical coherence tomography imaging at 370,000 lines/s,” Opt. Lett. 31,2975–2977 (2006). [CrossRef] [PubMed]
14. V. XD. Yang, M. Gordon, B. Qi, J. Pekar, S. Lo, E. Seng-Yue, A. Mok, B. Wilson, and I. Vitkin, “High speed, wide velocity dynamic range Doppler optical coherence tomography (Part I): System design, signal processing, and performance,” Opt. Express 11,794–809 (2003). [CrossRef] [PubMed]
15. V. XD. Yang, N. Munce, J. Pekar, M. L. Gordon, S. Lo, N. E. Marcon, B. C. Wilson, and I. A. Vitkin, “Micromachined array tip for multifocus fiber-based optical coherence tomography,” Opt. Lett. 29,1754–1756 (2004). [CrossRef] [PubMed]
16. E. Chérin, R. Williams, A. Needles, G. Liu, C. White, A. S. Brown, Y-Q. Zhou, and F. S. Foster, “Ultrahigh frame rate retrospective ultrasound microimaging and blood flow visualization in mice in vivo,” Ultrasound Med Biol. 32,683–691 (2006). [CrossRef] [PubMed]
17. R. Huber, M. Wojtkowski, J. G. Fujimoto, J. Y. Jiang, and A. E. Cable, “Three-dimensional and C-mode OCT imaging with a compact, frequency swept laser source at 1300 nm,” Opt. Express 13,10523–10538 (2005). [CrossRef] [PubMed]
18. W. A. Reed, M. F. Yan, and M. J. Schnitzer, “Gradient-index fiber-optic microprobes for minimally invasive in vivo low-coherence interferometry,” Opt. Lett. 27,1794–1796 (2002). [CrossRef]
19. H. Li, B. A. Standish, A. Mariampillai, N. R. Munce, Y. Mao, S. Chiu, N. E. Marcon, B. C. Wilson, I. A. Vitkin, and V. XD. Yang, “Feasibility of Interstitial Doppler Optical Coherence Tomography for in vivo detection of microvascular changes during photodynamic therapy,” Lasers Surg. Med. 38,754–761 (2006). [CrossRef] [PubMed]
20. C. Kasai, K. Namekawa, A. Koyano, and R. Omoto, “Real-time two-dimensional blood flow imaging using an autocorrelation technique,” IEEE Trans. Sonics. Ultrason. 32458–464 (1985). [CrossRef]
22. J. F. de Boer, B. Cense, B. H. Park, M. C. Pierce, G. J. Tearney, and B. E. Bouma, “Improved signal-to-noise ratio in spectral-domain compared with time-domain optical coherence tomography,” Opt. Lett. 282067–2069 (2003). [CrossRef] [PubMed]