25 April 2013 In vivo lung microvasculature visualized in three dimensions using fiber-optic color Doppler optical coherence tomography
Author Affiliations +
J. of Biomedical Optics, 18(5), 050501 (2013). doi:10.1117/1.JBO.18.5.050501
For the first time, the use of fiber-optic color Doppler optical coherence tomography (CDOCT) to map in vivo the three-dimensional (3-D) vascular network of airway segments in human lungs is demonstrated. Visualizing the 3-D vascular network in the lungs may provide new opportunities for detecting and monitoring lung diseases such as asthma, chronic obstructive pulmonary disease, and lung cancer. Our CDOCT instrument employs a rotary fiberoptic probe that provides simultaneous two-dimensional (2- D) real-time structural optical coherence tomography (OCT) and CDOCT imaging at frame rates up to 12.5 frames per second. Controlled pullback of the probe allows 3-D vascular mapping in airway segments up to 50 mm in length in a single acquisition. We demonstrate the ability of CDOCT to map both small and large vessels. In one example, CDOCT imaging allows assignment of a feature in the structural OCT image as a large (∼1 mm diameter) blood vessel. In a second example, a smaller vessel (∼80 μm diameter) that is indistinguishable in the structural OCT image is fully visualized in 3-D using CDOCT.
Lee, Ohtani, MacAulay, McWilliams, Shaipanich, Yang, Lam, and Lane: In vivo lung microvasculature visualized in three dimensions using fiber-optic color Doppler optical coherence tomography



Vascular remodeling and angiogenesis are important processes in the pathogenesis of many lung diseases such as asthma, chronic obstructive pulmonary disease, and lung cancer.1 To date, structural information on vascular remodeling in the lung has been primarily investigated in ex vivo surgical and postmortem specimens which cannot be used to measure dynamic vasculature parameters such as blood flow speed and directionality. Also, methods that require tissue removal do not permit longitudinal studies to follow vascular changes due to disease progression or response to treatment.

Real-time minimally invasive techniques for mapping the three-dimensional (3-D) vascular network in the airway wall over time may provide new clinical tools for detecting and monitoring lung diseases. Current methods for the visualization of airway vasculature in vivo include computed tomography angiography and Doppler ultrasound. However, these methods lack the resolution and sensitivity to detect small vessels. Optical coherence tomography (OCT),2 has proven to have sufficient resolution for examining tissue morphology and microstructure. Fiber-optic probes have demonstrated in vivo endoscopic structural OCT imaging of the lungs.34.5

Several techniques based on OCT for visualizing vascular networks have been presented in the literature including several variants of Doppler OCT imaging,6,7 and structural methods such as speckle variance OCT,8 and speckle decorrelation.9 The vast majority of in vivo vascular detection techniques have been developed for geometries where access to the region of interest is external to the organism using large-scale galvanometer scanners. Fiber probe-based OCT vascular imaging is considerably more difficult than its galvo-scanned analog because artifacts caused by subject motion are compounded by the additional motion of the fiber probes and their scanning mechanisms. Nonetheless, two-dimensional (2-D) vascular detection inside internal organs using fiber optic based probes using color Doppler optical coherence tomography (CDOCT) has been demonstrated.10 Here we present the first 3-D CDOCT-generated vascular maps of in vivo human airways. In addition to structural information on the tissue morphology, it is hoped that inclusion of vascular maps from CDOCT imaging will provide additional diagnostic capabilities for lung diseases.


Methods and Materials

A 50.4 kHz swept source laser (SSOCT-1310, Axsun Technologies Inc., Billerica, Massachusetts) with a polarized output power of 20 mW centered at wavelength λ=1310nm and a bandwidth of 100nm (full width half maximum) was coupled 90% sample/10% reference into a single-mode Mach–Zehnder OCT interferometer. The sample arm was connected using a fiber-optic rotary joint (MJP-SAPB, Princetel Inc., Pennington, New Jersey) to a 0.9 mm diameter, rotationally driven, side-looking fiber optic probe (C7 Dragonfly Imaging Catheter, St. Jude Medical Inc., St. Paul, Minnesota). Light exits the probe at approximately 64 deg relative to the tip of the fiber optic probe (i.e., slightly forward looking). For in vivo lung imaging, the open-ended probe was inserted into a closed-ended 1.5 mm diameter sheath to prevent direct probe-patient contact and was guided to the imaging location using a bronchoscope. The probe was driven using a custom built motor capable of rotational speeds up to 100 Hz and pullbacks up to 50 mm in length. In this study, the probe was spun at speeds of 6.25 Hz (8064 A-lines/frame) and 12.5 Hz (4032 A-lines/frame). The interference signal was detected using a pair of 75 MHz balanced detectors (PDB420C, Thorlabs-Inc., Newton, New Jersey) using a polarization diversity detection scheme. A high speed digitizer (ATS9350, Alazar Technologies Inc., Pointe-Claire, Quebec) in “k-clock” acquisition mode collected the data. The data was processed to provide live-streaming 2-D structural OCT and CDOCT imaging using custom written data acquisition software. This study was approved by the Research Ethics Board of the University of British Columbia and the British Columbia Cancer Agency. All subjects gave written, informed consent.

The data processing steps to produce a single CDOCT image frame, based on previously reported algorithms7 are briefly described. The data is collected in the polar domain (r, θ) to generate the structural OCT image and raw CDOCT images using the Kasai velocity estimator:


where fD,(m,n) is the Doppler frequency shift of the pixel (m, n), fa is the A-scan frequency, and Im,n and Qm,n are the real and imaginary parts of the Fourier-transformed polar pixel (m, n). A high degree of oversampling in the θ-direction is required for reliable CDOCT processing for slower flow velocties deeper into the airway wall because fanning out of the A-lines in our rotational scanning geometry means that adjacent pixels get farther apart at large radii.11 To improve accuracy, the calculation in Eq. (1) is averaged over M pixels in the radial scan direction, and N pixels in the azimuthal scan direction. For our data sets, we chose N=32, and M=8 as a compromise between color Doppler SNR (which increases with higher values of N and M) and spatial resolution (which decreases with higher values of M and N). The mean velocity vm,n of the pixel (m, n) is v(m,n)=λ0fD,(m,n)/(2ntcosθ), where λ0 is the center wavelength, nt is the tissue refractive index, and θ is the Doppler angle. Assuming that blood vessels lie parallel to the axis of the airway (i.e., Dopplerangle=64deg) and nt=1.40, our nonaliased Doppler velocity measurement range is ±8.5mm/s. To reduce Doppler noise, a copy of the structural OCT image is first smoothed using an M by N rectangular kernel. The Doppler image is then masked to include only pixels for which the SNR’s of the corresponding smoothed structural OCT image pixels are greater than 2.0. There may be artifact-based Doppler signals following the preceding processing steps that are present due to a combination of three factors: (1) involuntary motion of the subject, (2) in-plane (r, θ) vibrational motion of the rotary probe, and/or (3) out-of-plane (z) pullback motion of the probe. These artifacts are largely removed using a bulk tissue motion correction algorithm that first creates a Doppler velocity histogram for the n’th A-line, Hν(n) from nw to n+w adjacent A-lines in the θ-dimension (in our case w=5), and specified ROI in the r-dimension. Then, the Doppler velocity of the peak in Hν(n) is subtracted from all Doppler velocities in the n’th A-line. This process is repeated for all A-lines in each frame. The CDOCT and structural OCT images are transformed into Cartesian coordinates to generate cross-sectional views of the airway. Finally, to display the location of the flow regions with respect to the structural image, velocities <1.7mm/s are masked out prior to overlay on top of the structural image. Raw data sets were reprocessed into movies using MATLAB (Mathworks, Natick, Massachusetts) and 3-D rendered using Amira (Visualization Sciences Group, Burlington, Massachusetts).


Results and Discussion

Figure 1 shows structural OCT and CDOCT imaging of a 15 mm long airway segment from a 68-year-old male subject. These images were collected with 8064 A-lines/frame (framerate=6.25Hz) and pullback speed=1mm/s (framepitch=0.16mm). Figure 1(a) shows the 3-D rendering of the structural OCT data with a representative 2-D slice shown in Fig. 1(b). Figure 1(c) and 1(d) show the same structural information as Fig. 1(a) and 1(b), respectively, with the CDOCT imaging information overlaid. The color scale for the CDOCT overlay is taken from color Doppler ultrasound convention that dictates that flow towards the probe is shown in shades of red, while flow away from the probe is shown in shades of blue. As light exits the probe slightly forward looking, the color map shows distal-to-proximal blood flow in shades of red, and proximal-to-distal blood flow in shades of blue. A movie of the structural OCT image side by side with the CDOCT overlay of this airway from distal to proximal is presented in Video 1. Figure 1(b) shows a round structure approximately 1 mm in diameter at the 12 o’clock position. The CDOCT image [Fig. 1(d)] clearly shows this structure contains flow, and thus indicates the structure to be a blood vessel that runs alongside the airway. In this case, identifying the architecture of this large vessel using CDOCT imaging alone can be complicated by the fact that the flow within the vessel is pulsatile, and causes the vessel to “blink” on and off in the CDOCT overlay image during the 3-D pullback (Video 1). However, once identified as vasculature using CDOCT imaging, the feature in the structural OCT movie can be used to map the entire vessel. The complete vessel is readily visible as a dark shadow in the 3-D rendering of the airway and is outlined with a dashed red line in Fig. 1(a). This example demonstrates the necessity of using both structural OCT and CDOCT imaging modalities to identify and map vasculature of this nature. The Doppler signal from the vessel is predominantly blue indicating blood flow proximal to distal along the airway. Quantification of velocity in this case requires Doppler phase unwrapping approaches as Doppler aliasing is readily apparent in Video 1.

Fig. 1

A 15-mm-long segment of airway from a 68-year-old male subject oriented distal to proximal left to right. (a) Three-dimensional (3-D) rendered structural optical coherence tomography (OCT) volume and (c) color Doppler optical coherence tomography (CDOCT) overlay onto the structural OCT volume of the OCT 3-D-pullback data. (b) Two-dimensional (2-D) cross-sectional structural OCT image and (d) CDOCT overlay on the structural OCT image taken from the locations indicated in (a, c). The scale bars in (b, d) are 0.5 mm. The CDOCT scale bar shown in (d) spans ±8.5mm/s. A 2-D structural OCT and CDOCT overlay movie of this airway from distal to proximal is presented in Video 1 (MPG, 5.3 MB) [URL: http://dx.doi.org/10.1117/10.1117/1.JBO.18.5.050501.1].


Figure 2 shows structural OCT and CDOCT imaging of an 8 mm long airway segment from a 52-year-old male subject with subfigures presented as in Fig. 1 with the corresponding movie shown in Video 2. These images were collected with 4032 A-scans/frame (framerate=12.5Hz) and pullback speed 0.6mm/s (framepitch=0.048mm). In the 3-D rendering shown in Fig. 2(c), a larger, deeper blood vessel has been cropped to reveal the underlying smaller, approximately 80 μm diameter vessel shown. This vessel follows the airway and is joined by another vessel of similar size near the proximal end. The blue color of vessel in CDOCT indicates that flow in this vessel is also proximal to distal along the airway. In this case, the small vessel is continuously visualized by CDOCT in Video 2 as it does not exhibit the same magnitude of pulsatility as the larger vessel seen in F. 1. Unlike the previous airway, there is no obvious feature in the structural OCT images that corresponds to the 80-μm diameter blood vessel. However, CDOCT unequivocally identifies and maps this blood vessel along the entire length of the airway.

Fig. 2

An 8-mm-long segment of airway from a 52-year-old male subject oriented distal to proximal left to right. Subfigures are presented as in Fig. 1. A 2-D structural OCT and CDOCT overlay movie of this airway from distal to proximal is presented in Video 2 (MPG, 8.6 MB) [URL: http://dx.doi.org/10.1117/10.1117/1.JBO.18.5.050501.2].


Comparing the CDOCT imaging quality between the two airways presented, Video 1 (8064 A-lines/frame) shows fewer radial artifacts due to laser phase instability than Video 2 (4032 A-lines/frame). However, the faster frame rate in Video 2 (12.5 fps) compared to Video 1 (6.25 fps) should mitigate subject motion artifacts more effectively and allows for finer imaging resolution along the airway axis.

In comparison to other OCT vasculature detection modalities, CDOCT is ideally suited for endoscope-guided vasculature detection in the lung. As this implementation of CDOCT only relies on intraframe data processing to extract vasculature, it is less sensitive to ever-present subject and probe motion that can cripple techniques that require frame-to-frame image stability. Rotating fiber optic probes can also be plagued by nonuniform rotational distortion that can complicate inter-frame algorithms even in subject-motion-free environments.



This study demonstrates for the first time 3-D imaging of in vivo lung vasculature using fiber-optic CDOCT. For diseases where vascular remodeling and angiogenesis play central roles, 3-D mapping of the lung vasculature may provide important insight into disease detection and diagnosis, following disease progression/regression, and interventional monitoring. In combination with morphological information gleaned from structural OCT imaging, vascular imaging is expected to increase the diagnostic capabilities of lung OCT imaging.


We thank Rosa Lopez Lisbona, Hamid Pahlevaninezhad, Myles McKinnon, Chulho Hyun, and Barry Vuong for technical assistance and helpful discussions. This work was supported by the Canadian Institutes of Health Research, the Natural Sciences and Engineering Research Council of Canada, the Michael Smith Foundation for Health Research, the Canada Research Chairs and the Canada Foundation for Innovation.



N. F. Voelkelet al., “Angiogenesis in chronic lung disease,” Chest J. 131(3), 874–879 (2007).CHETBF0012-3692http://dx.doi.org/10.1378/chest.06-2453Google Scholar


D. Huanget al., “Optical coherence tomography,” Science 254(5035), 1178–1181 (1991).SCIEAS0036-8075http://dx.doi.org/10.1126/science.1957169Google Scholar


M. Tsuboiet al., “Optical coherence tomography in the diagnosis of bronchial lesions,” Lung Cancer 49(3), 387–394 (2005).0169-5002http://dx.doi.org/10.1016/j.lungcan.2005.04.007Google Scholar


S. Lamet al., “In vivo optical coherence tomography imaging of preinvasive bronchial lesions,” Clin. Cancer Res. 14(7), 2006–2011 (2008).CCREF41078-0432http://dx.doi.org/10.1158/1078-0432.CCR-07-4418Google Scholar


L. P. Haririet al., “Volumetric optical frequency domain imaging of pulmonary pathology with precise correlation to histopathology,” Chest J. 143(1), 64–74 (2013).CHETBF0012-3692Google Scholar


J. A. Izattet al., “In vivo bidirectional color Doppler flow imaging of picoliter blood volumes using optical coherence tomography,” Opt. Lett. 22(18), 1439–1441 (1997).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.22.001439Google Scholar


V. X. D. Yanget al., “High speed, wide velocity dynamic range Doppler optical coherence tomography (Part I): system design, signal processing, and performance,” Opt. Express 11(7), 794–809 (2003).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.11.000794Google Scholar


A. Mariampillaiet al., “Speckle variance detection of microvasculature using swept-source optical coherence tomography,” Opt. Lett. 33(13), 1530–1532 (2008).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.33.001530Google Scholar


X. Liuet al., “Quantitative transverse flow measurement using optical coherence tomography speckle decorrelation analysis,” Opt. Lett. 38(5), 805–807 (2013).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.38.000805Google Scholar


B. A. Standishet al., “Doppler optical coherence tomography monitoring of microvascular tissue response during photodynamic therapy in an animal model of Barrett’s esophagus,” Gastrointest. Endosc. 66(2), 326–333 (2007).GAENBQ0016-5107http://dx.doi.org/10.1016/j.gie.2007.02.040Google Scholar


C. Sunet al., “In vivo feasibility of endovascular Doppler optical coherence tomography,” Biomed. Opt. Express 3(10), 2600–2610 (2012).BOEICL2156-7085http://dx.doi.org/10.1364/BOE.3.002600Google Scholar

Anthony M. D> Lee, Keishi Ohtani, Calum MacAulay, Annette McWilliams, Tawimas Shaipanich, Victor X. D. Yang, Stephen Lam, Pierre Lane, "In vivo lung microvasculature visualized in three dimensions using fiber-optic color Doppler optical coherence tomography," Journal of Biomedical Optics 18(5), 050501 (25 April 2013). http://dx.doi.org/10.1117/1.JBO.18.5.050501
Submission: Received ; Accepted

Optical coherence tomography


Fiber optics

Doppler tomography


In vivo imaging

Doppler effect


Back to Top