Over the last decade, various techniques have been reported to image microvascular networks, whose aims range from distinct research questions to the diagnosis and monitoring of specific pathologies. Color fundus photography, fluorescein angiography (FA), and indocyanine green angiography served as standard methods for imaging the retinal vascular structure.1 These imaging techniques require the injection of exogenous contrast agents into the circulation, which may cause allergic reaction. Moreover, these techniques cannot provide blood vessel depth information and are limited by the two-dimensional (2-D) nature of the images.
An attractive optical technique for imaging three-dimensional (3-D) microvasculature structure without the need of exogenous agents is vascular-sensitive optical coherence tomography (OCT). The OCT technique is based on the principle of low-coherence interferometry, and it offers the advantages of being noninvasive, contactless, and yields depth-resolved localization at high spatial and temporal resolutions cross-sectional imaging in biological systems.cross-sectional imaging in biological system.2,3 Moreover, intrinsic contrast mechanisms and depth resolution of tissue are advantages over confocal fluorescence microscopy. Unlike the sound echoes of ultrasound, OCT generates images based upon back-reflected light from the sample. Various OCT blood flow detection techniques have been developed for imaging vasculature structure, each of which has its own merits and drawbacks. For example, Doppler OCT22.214.171.124.126.96.36.199.13.14.–15 (also known as optical Doppler tomography), images are reconstructed based on phase changes or Doppler frequency shifts due to moving particles in the sample. Phase changes are calculated from multiple axial scans (A-scans) collected at the same position, or adjacent A-scans with sufficient spatial overlap. The DOCT method is useful in mapping blood vessels in various tissue such as skin16,17 and retina18 and can extract blood-flow information. Several useful strategies such as dual-beam DOCT1920.–21 and narrow bandwidth phase-reference OCT,22 are reported to improve the sensitivity of DOCT. All of these methods are phase-sensitive, and therefore a phase-stable system is necessary for obtaining high-contrast images.23
The power Doppler OCT and/or optical angiography/microangiography (OMAG) are other imaging techniques based on the Doppler effect that is able to isolate blood flow signals from the static tissue background.24188.8.131.52.184.108.40.206.33.34.35.–36 This method does not allow visualization of flow velocity directly; it provides a measure of the total backscattered optical signal attributable to moving targets with higher SNR37 compared to DOCT. Unlike phase-based approaches, OMAG directly analyzes and processes an image, which can minimize the phase instability noise. Different phase compensation methods3738.–39 were introduced to compensate motion artifacts from OMAG measurement. Overall real-time application of power Doppler or OMAG is computationally intensive and may require graphics processing unit (GPU)-based signal processing techniques to clinical suitability.
Two recently developed high-contrast in vivo 2-D/3-D microcirculation imaging techniques, with various implementations by research groups, are speckle-variance (SV) and phase-variance (PV) OCT. SV-OCT images microvasculature by calculating the interline or interframe speckle variance of the intensity-based structural OCT images.4041.42.43.44.45.–46 The calculation of the variance of the signal intensity is determined either by a spatial window40 or temporal window.41,42 The spatial window method uses the speckle of the conventional structural OCT images.40 However, the performance of this system is hindered by the computational complexity and slow frame rate in the initial demonstration.41 Developments of high-speed swept-source laser16,17 may yield sufficient imaging speed of this technique. The temporal window method41,42 is based on calculating the interframe speckle variations between successive structural (B-mode) OCT images. This method provides high-contrast microvasculature images with enhanced sensitivity and angle independence, without the need of a high-phase stability system. However, the images suffer from blood vessel shadowing artifacts47,48 and bulk tissue motion.41 Szkulmowski et al.49 averaged the multiple A-scans, and Lee et al.50 coded the discrete Fourier transform algorithm to improve the speckle contrast from the OCT images with the capability of real-time signal processing.
The PV-OCT method47,5220.127.116.11.18.104.22.168.22.214.171.124.64.65.–66 has been implemented to visualize the vasculature of zebrafish,59 ocular circulation of the mouse,60 and human retinal vasculature.47,6126.96.36.199.–66 This method identifies the phase difference between consecutive B-scans of the same transverse position and allows for the mapping of microvasculature, independent of the direction of flow relative to the imaging system.47,60 High-speed phase stable imaging systems and postimage processing can reduce motion artifacts, which may further advance this technique for use in clinical applications.
In this article, we report on the principles, algorithms, current state of the art, and future applications of microvascular imaging based on SV and PV techniques. Studies are ongoing for improving complete microvascular metrics, as well as assessing blood flow information by optimizing system acquisition parameters or by improving the OCT system.
The term speckle is used to describe the OCT signal acquired as a result of partial coherently backscattered light from biological tissues.6768.69.70.–71 If an OCT image is acquired in a stationary object, the speckle pattern is temporally stationary as well.44 Due to the presence of moving particles in biological tissues (such as red blood cells), the speckle pattern varies with time and can be quantified by speckle variance calculations using either interframe or interline comparisons.
Fundamentally, SV imaging is a decorrelative measurement technique that has the capability to visualize blood flow, independent of the Doppler angle. This method is suitable for detecting small blood vessel networks because of its angle independence and fast image acquisition capability. Decorrelation methods based on speckle intensity variation in B-scan images were previously introduced in ultrasound systems.72,73 In OCT, SV imaging has the capability of detecting microvasculature through intrinsic contrast in the scattering properties of red blood cells as they move through the vascular network. These speckle variance images () are created by evaluating variance in the structural pixel intensity () across the desired number () of B-scan images at the same spatial location using the following equation:
The SV estimate from Eq. (1) differs in solid and fluid components of biological tissue, as they have different time-dependent scattering properties, which provide the intrinsic contrast in speckle variance. For example, in the regions of solid tissue, assuming no bulk tissue motion (BTM), the pixel intensities have a Gaussian distribution profile [Fig. 1(a)]. Contrary to this are the fluid components of biological tissue, such as red blood cells traveling through a vessel, where pixel intensities change and lead to a Rayleigh distribution.74
As a consequence, the SV contrast between solid and fluid components of biological tissue can affect the accuracy of the measurement. To account for this, the SV SNR between a “fluid” and “solid” pixel is defined as1); and are the variance of and , respectively; and is the equal time-averaged intensity of “fluid” and “solid” pixel. To achieve higher SNR, the gate length () was optimized and was found for lower BTM scenarios.42
A major challenge with using the SV-OCT technique is the effect of blood vessel shadowing artifacts.75 The streaking and/or shadowing artifacts encountered as one moves through the blood vessel to the deeper layers of the tissue sample can lead to an overestimation of the vascular area (Fig. 2). To remove these artifacts from the SV image, a step-down exponential filter could be used, while preserving the overall vascular structure.11,47,48,75 In this de-shadowing process, the SV signal is attenuated by a factor from the top surface of the tissue and is proportional to the sum of the de-shadowed SV pixel above it. Mathematically, the de-shadowing process is expressed asFig. 2) and qualitatively found no loss of vascular information due to the overall attenuation of speckle variance intensity (Fig. 3).
Subsequent to the attenuation of the shadowing artifact, a depth-encoding colormap (Fig. 4) is used to preserve vascular depth information in 2-D.76 Each image in the stack is first encoded with an RGB color depending on its depth in the stack. The transparency (or alpha channel) for each pixel in the stack is then set to its SV intensity. Both RGB color and alpha channels are then combined to get color depth-encoded SV images.
An alternative approach for vasculature visualization considers the evaluation of phase change statistics as a contrast mechanism. Various phase-contrast methods, including Doppler imaging and power Doppler tomography (or angiography) and the PV method,51,53,55,62 were covered in the introduction. Here, we focus on PV-OCT principles and algorithms as it produces high-contrast microcirculation imaging performance over DOCT and optical microangiography.
PV-OCT utilizes the statistical variance of phase changes over successive OCT scans where phase changes are caused by motion of fluid in vessels or capillaries. It produces high-contrast imaging of 3-D capillary networks. The variance of phase changes can be calculated over successive A-scans during transverse scanning.24,36,5152.–53 The average phase shift () between sequential A-scans at the same location () is expressed as5152.–53
If denotes the standard deviation of the Doppler spectrum, the Doppler variance () can be derived from autocorrelation technique.25,51,77 Usually, averaging is done to improve SNR. The averaging can be performed either in the lateral direction51 or in both lateral (temporal) and axial (depth) directions.24,25,52 Compared to lateral averaging, only Eq. (5), averaging along both the lateral and depth directions [Eq. (6)], can reduce the background noise and improve the image quality:25,36
Phase variance can also be calculated from phase changes between B-scans, where multiple B-scans are required over the same transverse scan.61,78 Thus, multiple phase measurements are recorded over the same transverse point separated in time. The phase differences between adjacent frames of the same transverse points can then be obtained, and the variance of these changes can be calculated. The phase difference () between two consecutive B-scans at the same location is mathematically expressed as
There are three basic steps in PV-OCT measurement: OCT data acquisition, phase change calculation, and phase variance calculation. Fingler et al.59 proposed two different scan mechanisms, MB-scan and BM-scan, to acquire 2-D OCT images with a variety of time separations between phase measurements. MB-scan takes multiple M-scans over different transverse positions; on the other hand, multiple B-scans over time at the same transverse scan position are referred to as BM-scans. Densely sampled A-scans may be required to improve the transverse resolution. It was demonstrated that the phase contrast calculation through an MB-scan could allow cancellation of phase noise-induced error. Unfortunately, MB-scans are time consuming and limit the speed of 2-D phase contrast imaging. On the other hand, BM-scans are more time efficient and commonly used. Additional steps of phase error reduction are needed ultimately to improve the accuracy of PV measurement.59
To improve the measurement accuracy, image processing techniques are applied during the PV-OCT procedure. First, to improve the accuracy of phase change calculation, the bulk motion along the imaging direction–induced artifact needs to be compensated. Bulk motion within each frame can be removed by using the Doppler shifts between adjacent A-lines, whereas bulk motion between adjacent frames can be removed by correlation-based algorithms. Fingler et al.59 introduced a weighted mean technique that utilizes phase changes from all depths within a chosen region of the A-scan to calculate the bulk motion. The phase change is then corrected by subtracting the tissue motion as47,75,79 are often used to remove phase shifts caused by bulk motion (e.g., eye motion). Liu et al.36,80 showed an intensity-based Doppler variance method that is insensitive to bulk motion and can be used without the removal of bulk motion.
We have outlined the principles and algorithms of SV and PV techniques in the above sections. In the following sections, we outline specific applications of SV- and PV-OCT microvascular imaging, putting particular emphasis on the broad application of high-resolution methods to address specific research questions.
Over the past several years, SV- and PV-OCT has been applied to a variety of preclinical and clinical scenarios requiring high-resolution imaging of the microvasculature. In this section, we summarize these applications and highlight both the diversity of use and limitations of these technologies. To begin, we highlight the study of tumor biology, as this was the first reported use of SV-OCT in the literature. Following this introduction, we highlight a variety of applications, including imaging the embryo, applications for the study of neuroscience and ophthalmology, and imaging nanoparticles, and conclude with methodological developments that aim to improve certain aspects of SV-OCT imaging across disciplines. In the second part of this section, we highlight applications of the PV-OCT technique.
Applications of SV-OCT
Imaging of tumor biology
Mariampillai et al.41,42 successfully demonstrated the SV technique to image microvasculatures in animal skin using a window chamber model (Fig. 5). The technique provides Doppler angle-independent microvessel imaging, which overcame the limitation of the DOCT technique.9 In the initial work, the smallest detectable microvessels by SV imaging was in the range of 25 µm.41 One major limitation highlighted by the authors is bulk tissue motion, which prompted further work to improve SV imaging using a different scanning pattern to improve imaging speed and SV contrast.42 In addition to 2-D imaging, the SV vasculature map can be color-coded11 according to the depth of vessels in the tissue (Fig. 6). Unlike the structural OCT image [Fig. 6(a)], satellite metastases are clearly observed in SV images, as indicated with blue arrows in Fig. 6(b). Significant improvement of microvessel contrast was found when the number of frames () increased from 2 to 8. The SV approach showed better performance for detecting microvasculature over color DOCT, due to its sensitivity and angle independence; however, it still suffers from blood vessel shadowing artifacts.47,48
To detect locations of retinal microvasculature and assess blood oxygen level effectively, Liu et al.44 combined speckle variance with spectroscopic OCT (SOCT) techniques. The SV-OCT images were calculated [from Eq. (1)] and compared the effect of using different numbers of frames (, 100, and 500) of OCT images processed. Large numbers of frames required can increase image acquisition and data processing time. Simple thresholding when applied in the SV images could allow determination of the location and dimension of the capillary blood vessels [Fig. 7(d)]. The SOCT image was then combined with the thresholded SV-OCT image to determine the spectroscopic properties of blood vessels, such as blood oxygen level and speckle noise, based on the attenuation properties of the blood in different wavelength ranges and the results are shown in Fig. 7(e) and 7(f). The spatial variations of blood vessels were represented by the full dynamic range of color code; red and purple represent the minimum and maximum speckle variance in the blood vessels. The combined results showed the localized spectroscopic properties of the blood microvascular structure with future potential of functional OCT imaging applications.
A better understanding of structural and functional properties of tumor microcirculation can potentially contribute to monitoring, preventing, and improving the treatment of tumor or cancer therapies. Formation of new vessels, remodeling, and dropout of exiting vessels inside a tumor can continuously create temporal variations in blood flow and oxygenation.66 For the study of tumor hemodynamic, Skala et al.81 combined the hyperspectral imaging with Doppler and SV-OCT to noninvasively microvasculature, hemoglobin oxygen saturation level (), flow information, and sheer rate on the vessel wall with tumor growth (Fig. 8). Eight repeated B-scans were used to collect SV-OCT images, whereas for DOCT images, 10 repeated A-scans were used. Combining such measurements reveal potential for clinical applications requiring detection of flow velocity, vessel diameter and morphology, as well as oxygen saturation level inside tumors, which can be determined accurately without the use of contrast agents.
The SV-OCT technique can also be used to monitor the change of tissue properties during photodynamic therapy (PDT) and radiation-based cancer treatments.41,8283.–84 Mariampillai et al.41 showed the detection of microvascular changes induced by PDT before, during, and immediately after posttreatment. The smaller blood vessels disappeared approximately 1 min after the PDT laser exposure [Fig. 9(b)]. A transient complete shutdown of one of the larger vessels occurred between 5 min, 40 s, and 6 min, 35 s, with a partial reperfusion event at 6 min, 27 s [Fig. 9(c)–9(e)]. At the end of the PDT laser irradiation, the smaller vessels were completely closed, whereas the larger vessels had an approximately 30% reduction of diameter as measured by the SV-OCT technique [Fig. 9(f)].
Imaging oral tissue
Davoudi et al.85 reported 3-D microstructural morphology and compared images of healthy [Fig. 10(a)] and scar regions [Fig. 10(b)] in the human oral cavity using the SV technique. From initial case studies, the average volumetric vessel densities for scar and healthy region were equal to 12% and 8%, respectively. High concentrations of microvessels were observed in the scar region. In both images, the minimum detectable diameters of microvessels were found to be under clinical conditions.
Imaging embryological tissue
Imaging the microvascular structure of the mammalian embryos has great potential for the study of cardiovascular birth defects and potential treatments. Recently, Sudheendran et al.86 illustrated the vascular development in the cultured mouse embryos with the SV-OCT technique. The results showed that the embryonic yolk sac structure of the mouse undergoes significant changes at 8.5 and 9.5 days post coitus (dpc) similar to other embryonic cardiovascular studies.87 The vasculature was clearly visible at both stages with apparent morphological difference found in the primitive vascular plexus at 8.5 dpc [Fig. 11(a)] and the remodeled hierarchically branching vascular network at 9.5 dpc [Fig. 11(b)].
For comparison, both DOCT and SV-OCT were used to reconstruct the same embryonic vascular structures (Fig. 12). The DOCT technique could not reconstruct several large vessels and capillaries in the central part of the image, where the embryonic yolk sac is perpendicular to the laser beam [Fig. 12(a)]. One the other hand, the corresponding SV-OCT technique provided more accurate representation of the vasculature of the yolk sac regardless of the flow direction [Fig. 12(b)]. In contrast, both methods outline vasculature equally well when blood flow has significant axial component [Fig. 12(c) and 12(d)]. Various colors in Doppler images represent the different directions of the axial flow components.
Applications in neuroscience
One of the most recent applications of the SV-OCT technique has been to visualize the intramedullar vasculature of the rodent spinal cord, opening a plethora of opportunities for the neuroscience community. Cadotte et al.88 successfully demonstrated the 3-D microvascular structure of the rodent spinal cord using the SV-OCT technique. In the case of imaging the spinal cord of a rat (Fig. 13), a spinal stabilizer was used to reduce the bulk cardio-respiratory motion. During a breath-hold cycle, structural OCT data were acquired [Fig. 13(b)]. The vascular imaging depth was found to be approximately 200 µm. The dorsal vein and its extensive tributaries can prevent microvascular imaging at deeper depth. While imaging the spinal cord of a mouse, less cardio-respiratory motion was observed, likely related to the different animal model used. A simple spinal stabilization jig was used to avoid bulk tissue motion. The structural OCT technique was capable of imaging approximately 500 µm into the mouse spinal cord, visualizing the entire dorsal gray matter and a portion of the ventral gray matter [Fig. 14(a)]. The SV-OCT and depth-dependent color mapping images are shown in Fig. 14(b) and 14(c), respectively. The diameter of the resolving microvessels in the mouse spinal cord was approximately 10 to 20 µm. By imaging a 3-D volume of tissue, this high temporal (μsec) and spatial resolution () SV-OCT method offers the possibility of studying aspects of neurovascular coupling, whereby the optically invisible electrical events of neuronal signaling drive changes in the microvascular network to meet spatiotemporally isolated increases in metabolic requirements. While a host of other optical imaging techniques have been used to study the phenomenon of neurovascular coupling in the cerebral cortex, SV-OCT holds promise to uncover afferent or efferent stimulus-specific hemodynamic response functions in the spinal cord with high temporal resolution.
Another aspect of the SV technique is that it may be used to trace the diffusion of nanoparticles (i.e., Au nanorings) in a sample and observe its systematical diffusion behaviors.89 With resonant Au nanorings (Au-NRs), the distributions of nanorings were deeper inside the sample with time. However, with a high concentration of the Au-NRs solution, the overall localized surface plasmon resonance (LSPR) behavior became weaker than that of individual Au-NRs. As a result, the overall intensity of the OCT signal could appear attenuated, which limited further analysis from SVOCT. The technique’s potential remains to be proven in future studies under physiologically tolerable Au-NR concentrations.
Application in ophthalmology
Another application of the SV technique is in the field of ophthalmology. For this purpose, Jia et al.46 developed the split-spectrum amplitude-decorrelation (SSAD) algorithm, which has significantly improved the image contrast and allowed for the visualization of the microvasculature of the posterior part of the human eye. This method may be used to diagnose and/or treat blinding diseases such as glaucoma, diabetic retinopathy, and age-related macular degeneration. The algorithm consists of dividing the OCT spectrum into several spectral (wavenumber) bands and then calculates the inter-B-scan decorrelation of each of these spectral bands, both separately and averaged. The degradation in the axial resolution caused by splitting up the spectrum is a major drawback of this implementation. Figure 15 showed the in vivo optic nerve head (ONH) angiograms obtained by the reflectance intensity images and decorrelation images computed with SSAD algorithm. The en face reflectance intensity image [Fig. 15(a)] showed all major blood vessels; but capillaries and microvasculature of the retina, choroid, and optic disc were not visible during imaging. On the other hand, the en face decorrelation image [Fig. 15(c)] showed all major retinal vessels, as well as capillaries and the microvascular network of the ONH. Layers of tissue morphology, such as retina, choroid, and sclera, can be delineated in the OCT cross-section (B-mode) image [Fig. 15(b)]. In the decorrelation cross-section image [Fig. 15(d)], the disc blood vessels (marked with green arrows) were clearly visualized approximately 1.0 mm deep inside the disc. The SSAD algorithm was also tested on in vivo images of the human macula (Fig. 16). The vascular pattern and capillary networks and the foveal avascular zone (FAZ) were clearly visualized in the retinal circulation en face image [Fig. 16(a)]. In an ideal case, the vessel should not be present in the FAZ region; however, few false branches were found in this region due to noise.46 Previous studies suggested that the choriocapillaris layer develops a confluent overlapping plexus; therefore, an en face decorrelation image of the choroid circulation [Fig. 16(b)] showed confluent flow. The cross-sectional OCT images [Fig. 16(c) and 16(d)] of the retinal vessels [marked with upper and lower dashed lines in Fig. 16(a)] showed different tissue layers from the nerve fiber layer (NFL) to the outer plexiform layer.
Applications of PV-OCT
PV-OCT is widely used to create vascular maps of human eyes.47,6188.8.131.52.–66 Kim et al.61,64,65 demonstrated the feasibility of a high-speed PV-OCT system for 2-D visualization of the capillary network over the foveal and perifoveal regions, as well as 3-D representation of the vasculature in the human retina. A spectral domain OCT system (at center wavelength , with 2048 pixels) was used for the PV measurement. The scanning protocols used a series of BM-scans (here “BM-scans” means that it takes multiple B-scans over time for the same transverse scan region) acquired across the fovea, with each BM-scan consisting of three successive B-scans at the same location. The imaging acquisition time of PV-OCT is 3.5 s. Phase differences between sequential B-scans were extracted for PV contrast calculation. The imaging results were shown in Fig. 17.
The size of the smallest resolvable capillary detected from this acquisition scheme is 10 to 12 µm. PV-OCT shows details of capillary networks comparable to the gold standard fundus FA imaging. In addition, PV-OCT provides depth information of capillary locations. They also imaged multiple volumes and vascular maps to visualize a larger region of the retina, demonstrating potential for realistic clinical applications. The stitched vascular map and fundus FA are shown in Fig. 18. Depth-resolved color maps are shown in Fig. 18(b). Such PV-OCT imaging may be useful to improve the early diagnosis of retinal vascular diseases, including diabetic retinopathy and vascular-related macular degeneration.
Furthermore, Kim et al.78 also demonstrated en face visualization of human retinal microcirculation using average intensity OCT, SV-OCT, PV-OCT, and FA. The fundus FA photo of the magnified foveal region is shown in Fig. 19(a), and 2-D projection images generated by the average intensity projection of the retina is presented in Fig. 19(b). The bright spot (marked as a red dotted circle) in the center of the foveal region in Fig. 19(c) is the back-reflection signal from the central fovea. The en face projection image of the entire retina achieved by PV-OCT [Fig. 19(d)] shows similar results as the image from SV-OCT method [Fig. 19(c)], excluding two horizontal stripes (marked as red arrows).
There is an important need to enhance the imaging depth for capturing both retinal and choroidal microvasculatures. PV-OCT systems at a 850-nm central wavelength offer good resolution and sensitivity for retinal imaging, but longer wavelengths may improve imaging of the deeper choroidal vessels. Swept-source OCT systems at a 1000-nm central wavelength with a 50-kHz line-rate (A-scan) has been developed.62,63 With an OCT system at a 1,300-nm central wavelength, the phase contrast method can image vasculature in the brain and tumors, but it is not the ideal choice for ophthalmic imaging due to the significant absorption by the vitreous. Therefore, based on the intended purpose of imaging, different OCT systems should be applied to specific applications.
Unlike the phase processing methods, including DOCT and PV-OCT, the SV method allows for a simple acquisition and post-processing algorithm and shows less sample motion artifacts effects and vessel shadows in the axial direction;66 however, this method produces low-contrast images of the vascular network due to the bright NFL and inner retinal structures. The SV-OCT method demonstrates similar microvascular maps with the FA and the PV-OCT when scanned over the central foveal region. However, SV values from the microcapillaries in the NFL and below the thick NFL lacked sufficient SNRs for reliable visualization of microcirculation in these regions. This is in contrast to PV-OCT methods, which visualized microcapillaries and larger vessels from the OCT-scanned volumes without such limitations. PV-OCT with rapid image acquisition at a high-speed A-scan rate can help to reduce acquisition time and related eye motion artifacts.
Recently, Liu et al.36,80 developed a modified Doppler variance algorithm that is based on the intensity (or amplitude) information rather than on phase information. Performances of this intensity-based Doppler variance (IBDV) method were compared to the traditional phase-resolved Doppler variance (PRDV) and phase-resolved color Doppler (PRCD) methods while mapping vascular network and extraction of blood vessel information.36 PRCD is most sensitive when flow direction is along the probing beam. Both PRDV and IBDV are better to measure blood flow when the flow direction is nearly perpendicular to the probing beam. However, PRDV and IBDV show Doppler angle-dependent behavior above a certain angle (). In applications where microvascular mapping (such as angiogram) is more important than the flow velocity quantification, both IBDV and PRDV images are high-contrast images and can detect smaller blood vessels (capillary) than that of PRCD images. However, the IBDV method is less sensitive to bulk motion with relatively simple post-processing steps, which is preferred for fast data processing.
Another alternative method known as correlation mapping OCT (cmOCT) has been developed to map vasculature networks and the blood flow information from the OCT images.90,91 The technique is capable of extracting blood flow information without the need of phase information. In addition, the system does not require any modification of the exiting OCT hardware. Parameters such as capillary density, vessel location, and diameter can be detected using this technique. The minimum detectable diameter of the microvessel volar forearm was achieved at about 33 µm.90 Although the system can image the microvascular structures, the exact locations of these vessels within the tissue were found previously. Further studies are necessary to improve this technique.
Discussion and Conclusions
In conclusion, we have reviewed the principles, algorithms, and applications of two microcirculation imaging techniques: SV- and PV-OCT, and compared their performances to DOCT and other gold standard techniques. The PV-OCT method can visualize vasculature structure independent of the vessel orientation and the flow velocity. However, this method is phase-sensitive; therefore, a highly phase-stable optical system is necessary. With a high-speed A-scan acquisition rate (in the kilohertz to megahertz range), B-scan image size and imaging speed can be increased to realize 3-D phase-variance contrast imaging more efficiently. Several potential applications of the PV technique are in vivo high-contrast vasculature imaging of the human retina for diagnosis of the ocular diseases, brain tissue vasculature, and study of vasculature in tumors/cancer cells. Further improvements of this technique require the minimization of shadowing artifacts, necessary to full 3-D, large-volume imaging of high-resolution microvasculature networks.
The SV-OCT is another attractive technique of detecting, quantifying, and monitoring microcirculation. Unlike DOCT and PV-OCT, SV methods are not phase-resolved, instead relying on the comparisons of consecutive structural (A-scan or B-scan) OCT images. This method is advantageous for detecting microvasculature because of its sensitivity, Doppler angle-independent microvascular information with little additional computational complexity and fast imaging technique. Lower data storage requirements relative to PV microvascular imaging techniques add further utility to this technique. Despite promising results, this technique also suffers from multiple scattering-induced blood vessel shadowing artifacts47,48 and susceptible to bulk tissue motion,41 which causes high interframe variance. High A-scan imaging systems, which ensure a sufficient number of acquired frames for calculation of speckle decorrelation, can diminish the bulk motion artifacts.50,92 Further study is necessary to overcome the limitations of accurately quantifying blood flow information.81
We summarized both SV- and PV-OCT technologies in terms of algorithms and their applications for microvasculature imaging. In comparison to other vascular-sensitive OCT imaging methods, both of these techniques have unique strengths and weaknesses that make them useful in different preclinical and clinical situations. Neither is categorically superior to the other, and both require ongoing methodological improvements specific to the research question at hand.
The authors thank the Canadian Institutes of Health Research (CIHR), Natural Sciences and Engineering Council, and Ryerson Bio-photonics and Bioengineering Laboratory for funding contributions. Dr. David Cadotte is supported by a CIHR Health Professional Fellowship award.
J. Gass, Stereoscopic Atlas of Macular Diseases: Diagnosis and Treatment, 4th ed., Vol. 2, Mosby, St. Louis, Missouri (1997).Google Scholar
A. F. Fercheret al., “In vivo optical coherence tomography,” Am. J. Ophthalmol. 116(1), 113–114 (1993).AJOPAA0002-9394Google Scholar
S. YazdanfarM. KulkarniJ. Izatt, “High-resolution imaging of in vivo cardiac dynamics using color Doppler optical coherence tomography,” Opt. Express 1(13), 424–431 (1997).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.1.000424Google Scholar
Z. Chenet al., “Noninvasive imaging of in vivo blood flow velocity using optical Doppler tomography,” Opt. Lett. 22(14), 1119–1121 (1997).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.22.001119Google Scholar
S. YazdanfarA. M. RollinsJ. A. Izatt, “Imaging and velocimetry of the human retinal circulation with color Doppler optical coherence tomography,” Opt. Lett. 25(19), 1448–1450 (2000).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.25.001448Google Scholar
R. Leitgebet al., “Real-time assessment of retinal blood flow with ultrafast acquisition by color Doppler Fourier domain optical coherence tomography,” Opt. Express 11(23), 3116–3121 (2003).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.11.003116Google Scholar
B. R. Whiteet al., “In vivo dynamic human retinal blood flow imaging using ultra-high-speed spectral domain optical Doppler tomography,” Opt. Express 11(25), 3490–3497 (2003).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.11.003490Google Scholar
L. Wanget al., “Frequency domain phase-resolved optical Doppler and Doppler variance tomography,” Opt. Commun. 242(4–6), 345–350 (2004).OPCOB80030-4018http://dx.doi.org/10.1016/j.optcom.2004.08.035Google Scholar
B. J. Vakocet al., “Three-dimensional microscopy of the tumor microenvironment in vivo using optical frequency domain imaging,” Nat. Med. 15(10), 1219–1223 (2009).1078-8956http://dx.doi.org/10.1038/nm.1971Google Scholar
T. SchmollC. KolbitschR. A. Leitgeb, “Ultra-high-speed volumetric tomography of human retinal blood flow,” Opt. Express 17(5), 4166–4176 (2009).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.17.004166Google Scholar
I. V. Larinaet al., “Hemodynamic measurements from individual blood cells in early mammalian embryos with Doppler swept source OCT,” Opt. Lett. 34(7), 986–988 (2009).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.34.000986Google Scholar
M. BonesiS. MatcherI. Meglinski, “Doppler optical coherence tomography in cardiovascular applications,” Laser Phys. 20(6), 1491–1499 (2010).LAPHEJ1054-660Xhttp://dx.doi.org/10.1134/S1054660X10110034Google Scholar
I. V. Larinaet al., “Live imaging of blood flow in mammalian embryos using Doppler swept-source optical coherence tomography,” J. Biomed. Opt. 13(6), 060506 (2008).JBOPFO1083-3668http://dx.doi.org/10.1117/1.3046716Google 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
Y. Zhaoet al., “Phase-resolved optical coherence tomography and optical Doppler tomography for imaging blood flow in human skin with fast scanning speed and high velocity sensitivity,” Opt. Lett. 25(2), 114 (2000).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.25.000114Google Scholar
S. YazdanfarA. M. RollinsJ. A. Izatt, “In vivo imaging of human retinal flow dynamics by color Doppler optical coherence tomography,” Arch. Opthalmol. 121(2), 235–239 (2003).AROPAW0003-9950http://dx.doi.org/10.1001/archopht.121.2.235Google Scholar
S. Zotteret al., “Visualization of microvasculature by dual-beam phase-resolved Doppler optical coherence tomography,” Opt. Express 19(2), 1217–1227 (2011).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.19.001217Google Scholar
S. Makitaet al., “Comprehensive in vivo micro-vascular imaging of the human eye by dual-beam-scan Doppler optical coherence angiography,” Opt. Express 19(2), 1271–1283 (2011).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.19.001271Google Scholar
C. J. Pedersenet al., “Phase-referenced Doppler optical coherence tomography in scattering media,” Opt. Lett. 30(16), 2125–2127 (2005).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.30.002125Google Scholar
M. Szkulmowskiet al., “Flow velocity estimation using joint spectral and time domain optical coherence tomography,” Opt. Express 16(9), 6008–6025 (2008).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.16.006008Google Scholar
H. Renet al., “Power optical Doppler tomography imaging of blood vessel in human skin and M-mode Doppler imaging of blood flow in chick chrioallantoic membrane,” Proc. SPIE 4956, 225–231 (2003).PSISDG0277-786Xhttp://dx.doi.org/10.1117/12.479028Google Scholar
V. 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
H. Renet al., “Real-time in vivo blood-flow imaging by moving scatterer-sensitive, spectral-domain optical Doppler tomography,” Opt. Lett. 31(7), 927–929 (2006).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.31.000927Google Scholar
R. K. Wanget al., “Optical microangiography provides depth-resolved images of directional ocular blood perfusion in posterior eye segment,” J. Biomed. Opt. 15(2), 020502 (2010).JBOPFO1083-3668http://dx.doi.org/10.1117/1.3353958Google Scholar
H. M. Subhashet al., “Volumetric in vivo imaging of microvascular perfusion within the intact cochlea in mice using ultra-high sensitive optical microangiography,” IEEE Trans. Med. Imag. 30(2), 224–230 (2011).ITMID40278-0062http://dx.doi.org/10.1109/TMI.2010.2072934Google Scholar
S. YousefiZ. ZhiR. Wang, “Eigendecomposition-based clutter filtering technique for optical micro-angiography,” IEEE Trans. Biomed. Eng. 58(8), 2316–2323 (2011).IEBEAX0018-9294http://dx.doi.org/10.1109/TBME.2011.2152839Google Scholar
F. Jaillonet al., “Enhanced imaging of choroidal vasculature by high-penetration and dual-velocity optical coherence angiography,” Biomed. Opt. Express 2(5), 1147–1158 (2011).BOEICL2156-7085http://dx.doi.org/10.1364/BOE.2.001147Google Scholar
L. AnJ. QinR. K. Wang, “Ultrahigh sensitive optical microangiography for in vivo imaging of microcirculations within human skin tissue beds,” Opt. Express 18(8), 8220–8228 (2010).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.18.008220Google Scholar
R. K. Wanget al., “Depth-resolved imaging of capillary networks in retina and choroid using ultrahigh sensitive optical microangiography,” Opt. Lett. 35(9), 1467–1469 (2010).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.35.001467Google Scholar
R. K. WangL. An, “Doppler optical micro-angiography for volumetric imaging of vascular perfusion in vivo,” Opt. Express 17(11), 8926–8940 (2009).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.17.008926Google Scholar
L. AnR. K. Wang, “In vivo volumetric imaging of vascular perfusion within human retina and choroids with optical micro-angiography,” Opt. Express 16(15), 11438–11452 (2008).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.16.011438Google Scholar
L. AnT. T. ShenR. K. Wang, “Using ultrahigh sensitive optical microangiography to achieve comprehensive depth resolved microvasculature mapping for human retina,” J. Biomed. Opt. 16(10), 106013 (2011).JBOPFO1083-3668http://dx.doi.org/10.1117/1.3642638Google Scholar
J. BartonS. Stromski, “Flow measurement without phase information in optical coherence tomography images,” Opt. Express 13(14), 5234–5239 (2005).OPEXFF1094-4087http://dx.doi.org/10.1364/OPEX.13.005234Google 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
L. ConroyR. S. DacostaI. A. Vitkin, “Quantifying tissue microvasculature with speckle variance optical coherence tomography,” Opt. Lett. 37(15), 3180–3182 (2012).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.37.003180Google Scholar
X. Liuet al., “Spectroscopic-speckle variance OCT for microvasculature detection and analysis,” Biomed. Opt. Express 2(11), 2995–3009 (2011).BOEICL2156-7085http://dx.doi.org/10.1364/BOE.2.002995Google 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
Y. Jiaet al., “Split-spectrum amplitude-decorrelation angiography with optical coherence tomography,” Opt. Express 20(4), 3116–3121 (2012).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.20.004710Google Scholar
J. Fingleret al., “Volumetric microvascular imaging of human retina using optical coherence tomography with a novel motion contrast technique,” Opt. Express 17(24), 22190–22200 (2009).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.17.022190Google Scholar
V. J. Srinivasanet al., “Rapid volumetric angiography of cortical microvasculature with optical coherence tomography,” Opt. Lett. 35(1), 43–45 (2010).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.35.000043Google Scholar
M. Szkulmowski, et al., “Efficient reduction of speckle noise in optical coherence tomography,” Opt. Express 20(2), 1337–1359 (2012).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.20.001337Google Scholar
K. K. C. Leeet al., “Real-time speckle variance swept-source optical coherence tomography using a graphics processing unit,” Biomed. Opt. Express 3(7), 1557–1564 (2012).BOEICL2156-7085http://dx.doi.org/10.1364/BOE.3.001557Google Scholar
Y. Zhaoet al., “Doppler standard deviation imaging for clinical monitoring of in vivo human skin blood flow,” Opt. Lett. 25(18), 1358–1360 (2000).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.25.001358Google Scholar
H. Renet al., “Phase-resolved functional optical coherence tomography: simultaneous imaging of in situ tissue structure, blood flow velocity, standard deviation, birefringence, and Stokes vectors in human skin,” Opt. Lett. 27(19), 1702–1704 (2002).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.27.001702Google Scholar
B. Parket al., “Real-time fiber-based multi-functional spectral-domain optical coherence tomography at 1.3 μm,” Opt. Express 13(11), 3931–3944 (2005).OPEXFF1094-4087http://dx.doi.org/10.1364/OPEX.13.003931Google Scholar
Y. Zhaoet al., “Three-dimensional reconstruction of in vivo blood vessels in human skin using phase-resolved optical Doppler tomography,” IEEE J. Sel. Top. Quant. Electron. 7(6), 931–935 (2001).IJSQEN1077-260Xhttp://dx.doi.org/10.1109/2944.983296Google Scholar
H. Renet al., “Imaging and quantifying transverse flow velocity with the Doppler bandwidth in a phase-resolved functional optical coherence tomography,” Opt. Lett. 27(6), 409–411 (2002).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.27.000409Google Scholar
Z. Dinget al., “Real-time phase-resolved optical coherence tomography and optical Doppler tomography,” Opt. Express 10(5), 236–245 (2002).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.10.000236Google Scholar
J. ZhangZ. Chen, “In vivo blood flow imaging by a swept laser source-based Fourier domain optical Doppler tomography,” Opt. Express 13(19), 7449–7457 (2005).OPEXFF1094-4087http://dx.doi.org/10.1364/OPEX.13.007449Google Scholar
J. Fingleret al., “Mobility and transverse flow visualization using phase variance contrast with spectral domain optical coherence tomography,” Opt. Express 15(20), 12636–12653 (2007).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.15.012636Google Scholar
J. Fingleret al., “Phase-contrast OCT imaging of transverse flows in the mouse retina and choroid,” Invest. Ophthalmol. Vis. Sci., 49(11), 5055–5059 (2008).IOVSDA0146-0404http://dx.doi.org/10.1167/iovs.07-1627Google Scholar
D. Y. Kimet al., “In vivo volumetric imaging of human retinal circulation with phase-variance optical coherence tomography,” Biomed. Opt. Express 2(6), 1504–1513 (2011).BOEICL2156-7085http://dx.doi.org/10.1364/BOE.2.001504Google Scholar
R. MotaghiannezamS. Fraser, “In vivo human retinal and choroidal vasculature visualization using differential phase contrast swept source optical coherence tomography at 1060 nm,” Proc. SPIE 8213, 821304 (2012).PSISDG0277-786Xhttp://dx.doi.org/10.1117/12.910991Google Scholar
S. M. R. MotaghiannezamD. KoosS. E. Fraser, “Differential phase-contrast, swept-source optical coherence tomography at 1060 nm for in vivo human retinal and choroidal vasculature visualization,” J. Biomed. Opt. 17(2), 026011 (2012).JBOPFO1083-3668http://dx.doi.org/10.1117/1.JBO.17.2.026011Google Scholar
D. Y. Kimet al., “Noninvasive imaging of the foveal avascular zone with high-speed, phase-variance optical coherence tomography,” Invest. Ophthalmol. Vis. Sci. 53(1), 85–92 (2012).IOVSDA0146-0404http://dx.doi.org/10.1167/iovs.11-8249Google Scholar
D. Y. Kimet al., “Visualization of human retinal micro-capillaries with phase contrast high-speed optical coherence tomography,” Proc. SPIE 7889, 78890H (2011).PSISDG0277-786Xhttp://dx.doi.org/10.1117/12.873896Google Scholar
T. Schmollet al., “Imaging of the parafoveal capillary network and its integrity analysis using fractal dimension,” Biomed. Opt. Express 2(5), 1159–1168 (2011).BOEICL2156-7085http://dx.doi.org/10.1364/BOE.2.001159Google Scholar
A. F. Fercher, “Inverse scattering, dispersion, and speckle in optical coherence tomography,” in Biological and Medical Physics, Biomedical Engineering, pp. 119–146, Springer-Verlag, Berlin, Heidelberg (2008).Google Scholar
V. Yanget al., “High frequency ultrasound speckle flow imaging comparison with Doppler optical coherence tomography (DOCT),” in Proc. IEEE Ultrasonics Symposium, Vol. 1, pp. 453–456, IEEE, Ontario, Canada (2004).Google Scholar
W. Aoudiet al., “Estimation methods for flow imaging with high-frequency ultrasound,” Ultrasonics 44(Suppl. 1), e135–e140 (2006).ULTRA30041-624Xhttp://dx.doi.org/10.1016/j.ultras.2006.06.025Google Scholar
J. W. Goodman, Statistical Optics, Vol. 1, pp. 361–465, Wiley, New York (1985).Google Scholar
V. X. D. Yanget al., “Improved phase-resolved optical Doppler tomography using the Kasai velocity estimator and histogram segmentation,” Opt. Commun. 208(4–6), 209–214 (2002).OPCOB80030-4018http://dx.doi.org/10.1016/S0030-4018(02)01501-8Google Scholar
A. Mariampillai, “Development of a high-resolution microvascular imaging toolkit for optical coherence tomography,” PhD Thesis, University of Toronto (2010).Google Scholar
C. Kasaiet al., “Real-time two-dimensional blood flow imaging using an autocorrelation technique,” IEEE Trans. Sonics Ultrasonics 32(3), 458–464 (1985).IESUAU0018-9537http://dx.doi.org/10.1109/T-SU.1985.31615Google Scholar
D. Y. Kimet al., “Visualization of human retinal capillary networks: a comparison of intensity, speckle-variance and phase-variance optical coherence tomography,” Proc. SPIE 8213, 821307 (2012).PSISDG0277-786Xhttp://dx.doi.org/10.1117/12.906940Google Scholar
G. Liuet al., “Intensity-based modified Doppler variance algorithm: application to phase instable and phase stable optical coherence tomography systems,” Opt. Express 19(12), 11429–11440 (2011).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.19.011429Google Scholar
M. C. Skalaet al., “Combined hyperspectral and spectral domain optical coherence tomography microscope for noninvasive hemodynamic imaging,” Opt. Lett. 34(3), 289–291 (2009).OPLEDP0146-9592http://dx.doi.org/10.1364/OL.34.000289Google Scholar
H. A. Collinset al., “Blood-vessel closure using photosensitizers engineered for two-photon excitation,” Nat. Photon. 2(7), 420–424 (2008).1749-4885http://dx.doi.org/10.1038/nphoton.2008.100Google Scholar
B. Davoudiet al., “Noninvasive in vivo structural and vascular imaging of human oral tissues with spectral domain optical coherence tomography,” Biomed. Opt. Express 3(5), 826–839 (2012).BOEICL2156-7085http://dx.doi.org/10.1364/BOE.3.000826Google Scholar
D. W. Cadotteet al., “Speckle variance optical coherence tomography of the rodent spinal cord: in vivo feasibility,” Biomed. Opt. Express 3(5), 911–919 (2012).BOEICL2156-7085http://dx.doi.org/10.1364/BOE.3.000911Google Scholar
C.-K. Leeet al., “Characterizing the localized surface plasmon resonance behaviors of Au nanorings and tracking their diffusion in bio-tissue with optical coherence tomography,” Opt. Express 1(4), 1060–1073 (2010).OPEXFF1094-4087http://dx.doi.org/10.1364/BOE.1.001060Google Scholar
E. JonathanJ. EnfieldM. J. Leahy, “Correlation mapping method for generating microcirculation morphology from optical coherence tomography (OCT) intensity images,” J. Biophoton. 4(9), 583–587 (2011).JBOIBX1864-063Xhttp://dx.doi.org/10.1002/jbio.201000103Google Scholar
J. EnfieldE. JonathanM. Leahy, “In vivo imaging of the microcirculation of the volar forearm using correlation mapping optical coherence tomography (cmOCT),” Opt. Express 2(5), 1184–1193 (2011).OPEXFF1094-4087http://dx.doi.org/10.1364/BOE.2.001184Google Scholar