Optical coherence tomography angiography (OCTA) was recently introduced for the imaging of microvascular networks in the human eye.1 Recent investigations using OCTA have been mainly focused on showcasing its applications in visualizing the flow characteristics of microvascular diseases within the macula, such as age-related macular degeneration (AMD) 2 and diabetic retinopathy (DR),3,4 the two leading causes of blindness worldwide.5,6 In addition, OCTA has documented microvascular abnormalities in other macular diseases such as branch retinal vein occlusion (BRVO), central retinal vein occlusion,7 and macular telangiectasia type 2 (MacTel2).8 Compared with conventional dye-based angiography, OCTA offers a safer, faster, and more cost effective alternative for ocular imaging of the macular vasculature. Overall, the combination of OCT and OCTA can present integrated structural and flow information of the human eye in vivo, opening new opportunities for both qualitative and quantitative analysis of ocular diseases.
Optical microangiography (OMAG)1 is one of the many OCTA approaches that utilizes the intrinsic properties of particles’ (e.g., red blood cells) motion, to highlight the contrast between signals due to red blood cells (RBCs) and signals due to static tissues. Unlike other OCTA approaches such as speckle variance,9 split-spectrum amplitude decorrelation angiography,10 and phase variance,11 OMAG is able to harness motion information to the fullest extent by exploiting both amplitude and phase information contained within the OCT signals. Therefore, OMAG is able to produce OCT angiography with better vascular connectivity, higher signal-to-noise ratio, and higher sensitivity to capillary blood flows.12 Consequently, OMAG is being established as a valuable tool for the investigation of ocular diseases both qualitatively 13 and quantitatively.14
To objectively assess the OCT angiograms, a number of studies have introduced vessel area density (VAD),15,16 vessel skeleton density (VSD),15,17 and nonperfusion area.16 These parameters have been reported useful in the analyses of disease progression in DR and AMD. While promising, there are limitations when evaluating indices individually. VAD provides information about the function of vascular network by calculating the percentage of area occupied by vessels on OCTA images. Thus, a change in VAD is coupled with changes in both vessel length and vessel size. Since the decrease of microvascular perfusion is sometimes accompanied with vessel dilation, VAD alone is only capable of providing partial information about vascular abnormalities. Similarly, VSD falls short because it only evaluates the vessel length without taking the vessel size information into account. Moreover, compared with fluorescein angiography (FA), it is difficult for OCTA to convey information about vascular leakage since it detects the existence of moving particles while FA senses the diffusion of injected dyes into the extravascular space. To our best knowledge, there has not been an attempt to extract vascular information other than vessel density from OCT angiograms.
In this study, we introduce a comprehensive quantitative OMAG analysis for OCTA vascular images. We propose to use the following metrics to quantitatively describe the OCTA images: VAD, VSD, vessel diameter index (VDI), vessel perimeter index (VPI), and vessel complexity index (VCI). We illustrate the usefulness of these metrics by using five clinical cases. We also assess the repeatability of these metrics by using OCTA scans captured from the same patients at different times. With our comprehensive quantitative OMAG strategy, we aim to provide clinicians and physicians with a rapid strategy for interpreting OCT angiograms from multiple perspectives and to help simplify the process of identifying pathological vascular features so as to enable the objective assessment of disease progression and treatment monitoring.
An automated tool with a user interactive interface was developed to quantify the abnormalities in the retinal vasculature from multiple perspectives. Included within this tool are the quantitative parameters of VAD, VSD, VPI, VDI, and VCI. To demonstrate the usefulness of these parameters, we used OMAG en face images generated by using both spectral domain OCT (SD-OCT) and swept-source OCT (SS-OCT) angiography prototype instruments.
SD-OCT data was acquired using a Cirrus HD 5000 SD-OCT angiography system prototype (Carl Zeiss Meditec Inc., Dublin, California), with an 840-nm central wavelength, 45-nm and a 68-KHz A-scan rate. The axial resolution is while lateral resolution is . The scanning protocol in this system provided a field of view (FOV) of , with 245 lines in each B-scan and 245 locations in the slow axis (where four repeated B-scans were acquired at each location). Equipped with FastTrac™18 technology, a custom sized montage scan could be acquired to give a widefield image with high resolution and high definition. In this study, customized scans consisting of nine motion-tracked smaller cubes, with a total FOV of , were collected from each subject.
SS-OCT data were acquired using a modified CIRRUS prototype containing a swept source laser provided by Carl Zeiss Meditec Inc., which had a 1050-nm central wavelength, 100-nm and 100-kHz A-scan rate. The axial resolution is while lateral resolution is . SS-OCT provided a FOV of with 300 A-lines per B-scan and a total of 300 B-scans centered at the fovea. To achieve OMAG imaging of the vascular networks, a previously described repeated B-scan protocol was implemented to acquire the volumetric OCTA data.13,18
All participants in this study were consented before OCTA imaging, following the tenets of Declaration of Helsinki, in compliance with the Health Insurance Portability and Accountability Act. The study was approved by the Institutional Review Board of the University of Washington (UW), Seattle, as well as the University of Miami Miller School of Medicine, Miami. All SD-OCT scans were obtained at the University of Washington Eye Institute, while all SS-OCT scans were obtained at the Bascom Palmer Eye Institute, University of Miami Miller School of Medicine. Subjects with the diagnosis of DR and BRVO were prospectively enrolled at the UW Medicine Eye Institute, and subjects with diagnosis of MacTel2 were enrolled in the MacTel project and imaged at the Bascom Palmer Eye Institute.
Figure 1 shows a flowchart providing an overview of the processes involved in our proposed quantitative approach to OCTA metrics. First, a two-dimensional en face OCTA/OMAG image [Fig. 2(a)] is generated from the volumetric data using a semiautomatic segmentation software.19 This image is processed into a binary image in MATLAB (R2015b, MathWorks, Inc.) using a global threshold, hessian filter,20 and adaptive threshold.15 The global threshold is determined by the noise level in the foveal avascular zone (FAZ) that is selected with a user interactive interface by asking the operator to click on the center and edge of FAZ. This global threshold is then applied to the entire image. Then, hessian filter and adaptive threshold are combined to generate a binary vessel map [Fig. 2(b)]. In doing so, the information of vessel area can be obtained. Second, a skeletonized vessel map is created [Fig. 2(c)], where every vessel, regardless of its size or diameter, is represented by a single pixel line. This way, the information of vessel length can be obtained. Last, the vessel perimeter or the edges of vessels are detected and contoured as the vessel perimeter map [Fig. 2(d)] to give the information about vessel perimeter. All five quantitative indices are calculated based on these three binary vessel maps, i.e., vessel area map, skeletonized vessel map, and vessel perimeter map. Note that there are only two values, 1 and 0, in the binary maps where 1 represents white pixel (i.e., existence of blood vessel) and 0 the black pixel (i.e., absence of blood vessel).
VAD is calculated as a unit less ratio of the total image area occupied by the vasculature to the total image area in the binary vessel mapsFig. 2(b)] and all the pixels in the vessel area map [all pixels on Fig. 2(b)]. are the pixel coordinates in the OCTA image (assuming with a size of pixel array). With this definition, VAD provides the best estimate of real vessel density as it takes both vessel length and vessel diameter into consideration. However, a complication would be that VAD remains unchanged when there are decreased perfusion and vessel dilation happening at the same time, potentially resulting in false negatives on the vascular abnormality. Furthermore, other than computing a single VAD for the whole OCT angiogram, we could also produce a spatially resolved VAD map. A down sampled block map was generated with (pixels) block size, the map was resized to original size and Gaussian filtered with a (pixels) kernel. All maps were generated in the same fashion.
VSD is calculated as the ratio of the length occupied by the blood vessels to the total area in the skeletonized vessel mapFig. 2(c)] and represents all the pixels in the skeletonized vessel map [all pixels on Fig. 2(c)]. The unit for VSD is pixel/pixel.2 VSD quantifies the vessel density by only considering whether or not the vessel exists, which is simply an assessment of the vessel length density regardless of the vessel diameters. Since every vessel is presented as a single pixel line, large vessels and small capillaries contribute equally to the VSD quantification. Consequently, compared with VAD, VSD would be more sensitive to the perfusion changes at the capillary levels.
VDI is calculated using both the vessel area map and the skeletonized vessel map to yield the averaged vessel caliber as follows:Fig. 2(b)] and represents the pixels registered as vessel length [white pixels on Fig. 2(c)]. The unit of VDI is pixel and can be converted into micron by multiplying VDI with pixel size, all VDI data in this article has been converted into microns. VDI calculates the averaged vessel caliber within the image; it does not reflect the change in vessel density. VDI presents the vessel size information regardless of the vessel length; therefore, it is sensitive to vascular dilation in the OCTA images. With the VDI map, localized vascular dilation would be easily identified and serve as a marker of vascular abnormalities.
VPI is calculated using the vessel perimeter map as the ratio of the vessel perimeter to the total area of the OCTA imageFig. 2(d)] and represents all the pixels in the vessel perimeter map [Fig. 2(d)]. The unit for VPI is pixel/pixel.2 Contrary to VSD, VPI uses the perimeter to present every vessel, which includes both the length and the diameter, similar to VAD.
VCI is calculated using both the vessel perimeter map and the vessel area map, defined asFig. 2(d)] and represents all the pixels of the vessel area [white pixels on Fig. 2(b)]. VCI is also a unit less parameter and derived from the digital image processing discipline.21 With the vessel complexity map, specific locations where the vascular morphology is less complicated could be identified. In most cases, decreased perfusion would result in less complex vasculature. To the best of our knowledge, this is the first time VCI has been introduced to quantify morphological characteristics of the retinal vasculature, though fractal dimension has been previously used to describe vascular tortuosity.15 To better illustrate the concept of VCI, an example study was conducted to test the robustness of this quantitative parameter. As shown in Fig. 3, five images were selected as examples for the calculation of VCI. Both binary area images and perimeter images are shown in Figs. 3(a)–3(j), and the VCI results for each image are plotted in Fig. 3(k), with the exact values denoted above each data point. We find a positive relationship between VCI and its visual complexity; as VCI increases, the visual complexity of the image also increases. Theoretically the VCI value should be 1 for a circle; however, due to the discrete nature of computing, computed VCI value for a circle would always be larger than 1, depends on sampling rate. All computed VCIs have been normalized according to this 1.5711 to 1 ratio.
Five illustrative cases are shown to demonstrate the usefulness of all the quantitative metrics that we propose.
Retinal Vasculature and Quantitative Maps in Normal Case
The three-dimensional (3-D) nature of OCT allows for segmentation of vascular plexi located at different physiological depths within the retina. Using a validated semiautomatic segmentation software,19 we separated the retinal microvasculature into three distinct physiological layers: the superficial vascular network encompassing from ganglion cell layer to inner plexiform layer, the deep vascular network including the outer plexiform layer, and the avascular outer retinal layer extending from outer nuclear layer to retinal pigment epithelium as shown in Fig. 4(b). We excluded the retinal nerve fiber layer in the following examples since the vasculature in this layer is morphologically different from the vasculature in the deeper retinal layers, especially for large FOV OCTA scans.22 However, this layer can be included in the analysis if needed. Figure 4(a) shows the en face OCTA image with maximum projection, with the color scheme for segmented vascular plexus shown in Fig. 4(b).The vessel density map, vessel diameter map, and vessel complexity map are shown in Figs. 4(c)–4(e). Figures 4(f)–4(h) were obtained by combining Figs. 4(c)–4(e) with the vessel map [Fig. 2(b)], which would make it easier to localize the particular vessels of interests. Good correlations can be observed in these three maps with the anatomic features of the retinal vasculature. For both vessel density map and vessel diameter map, large vessels in the superficial layer appear red, which represents higher density/diameter; in the vessel complexity map, large vessels are blue since they are less complex in morphology compared to capillaries. Figure 4(h) shows the overlay of vessel complexity map and deep retinal vascular layer en face maximum projection image, where we have observed that the red regions in the vessel complexity map correspond to the radial fan-like capillaries that originate from vertical anastomosis connecting superficial and deep network.
With the algorithm, VAD, VSD, VPI, VDI, and VCI can all be quantified from the OMAG images. For illustrative purposes, an example of five quantitative indices has been given in Table 1 based on the entire retinal image [Fig. 4(a)].
Quantitative analysis of a normal case and a NPDR case.
|Indices||VDI||VAD||VSD||VPI||VCI||Flow impairment zone (mm2)|
Microvascular Abnormality and Flow Impairment in Diabetic Retinopathy
Quantitative OMAG analysis of an eye with nonproliferative diabetic retinopathy (NPDR) is shown in Fig. 5. A late phase fluorescein angiographic image shown in Fig. 5(a) corresponds well to the OMAG image shown in Fig. 5(b). The areas of flow impairment are automatically detected according to our vessel skeleton map and superimposed onto OMAG en face image in green [Fig. 5(c)]. For this example, the negation of the vessel skeleton map was first obtained, and then the morphological operations were conducted to close smaller unconnected regions (600 pixels) that come naturally between vessels. The size of this region was determined according to normal OCTA scans and might differ among different OCT systems. Last, a Gaussian filter (with a kernel size of ) was used to smooth the edges of the detected regions of flow impairment. Vessel density, diameter, and complexity maps have been quantified as shown in Figs. 5(d)–5(f). Figures 5(g)–5(i) show the results by overlying the quantitative maps onto binary vessel area map. It can be observed from the quantitative maps that the regions of capillary flow impairment are often associated with an increase of vessel diameter and a decrease in VAD and vessel complexity, which agrees well with clinical features observed from DR patients. Specific quantitative results are shown in Table 1.
Repeatability of Quantitative Optical Microangiography in Branch Retinal Vein Occlusion and Normal Cases
Figure 6 shows a side-by-side comparison of two scans from the same subject diagnosed with BRVO, but scanned at time points one month apart. This case was clinically diagnosed as stable over the time period as judged by fundus examination, FA and visual acuity. OMAG results demonstrated that the detailed characteristics of BRVO such as vascular dilation and capillary nonperfusion are repeatable between these two scans. Quantitative maps [Figs 6(g)–6(l)] also show great repeatability, where the increase in vessel diameter and decrease in vessel density and vessel complexity are present in regions of the vein occlusion in both scans. VAD, VSD, VDI, VPI, and VCI have also been quantified on two scans. Results show high repeatability as the coefficients of variance are all smaller than 0.031 (Table 2). In doing so, each map was resized to 100 superpixels and repeatability between different scans was quantified based on these 100 samples.
Quantitative analysis and repeatability test of a BRVO case for two repeated scans.
|Coefficient of variance||0.017||0.014||0.013||0.031||0.014|
Similarly, two repeated scans of a normal subject (Fig. 7) show great repeatability of not only the OMAG en face images, but also three quantitative maps. Coefficients of variance have also been quantified for five quantitative indices and all are no greater than 0.021 (Table 3). Therefore, our OCTA quantification algorithm should be useful for monitoring disease progression and treatment response.
Repeatability test of a normal case.
|Coefficient of variance||0.018||0.009||0.021||0.011||0.017|
Moreover, we have implemented this algorithm in a software package where the operator can flexibly select their regions of interest on the en face image with the following options of e.g., a rectangle, circle, ring, or ellipse. Figure 8 shows an example of our region of interest (ROI) selection, Figure 8(a) is the original OMAG en face image, Figure 8(b) shows the quadratic analysis with a ring ROI, and Figs. 8(c)–8(f) shows the circle, rectangle, ellipse, and circle ROI selection, respectively. Table 4 shows the specific quantitative values when different ROIs are selected. If the selection is a ring, it can be divided into four quadrants to enable quadrant specific quantification. Table 5 shows results of quadratic analysis of this BRVO case. Note that in the scans centered at fovea, the operator can choose to center their ring ROI at the fovea centralis.
Quantitative analysis of a BRVO case with different choice of region of interest.
|VAD||VSD||VPI||VDI||VCI||ROI size (mm×mm)|
Quadratic quantitative analysis of a BRVO case.
|VAD||VSD||VPI||VDI||VCI||ROI size (mm×mm)|
Vascular Dilation in Macular Telangiectasia Type 2
Figure 9 shows an example of an early stage MacTel2 case, where the FA image [Fig. 9(a)] shows the signature juxtafoveal leakage that is more prominent temporally than nasally. The en face OMAG image and quantitative maps are given for the entire retina and deep capillary plexus in Figs. 9(b)–9(e) and 9(f)–9(i), respectively. The vessel diameter map and vessel complexity map demonstrate good accordance with the dye leakage on FA, and all demonstrate the telangiectatic changes in the microvasculature. With the dilation of the juxtafoveal microvasculature, local vessel diameter increases and the local vessel complexity decreases significantly, when compared with the relatively normal regions. The comparison between the entire retinal vasculature and deep capillary plexus also validates that the early telangiectatic changes mainly happen in deep capillary plexus23 as the quantitative maps manifest more obvious changes. With the segmentation strategy and high resolution OMAG images, a detailed depiction of the deep retinal microvasculature can be produced, as shown in Figs. 9(g)–9(i), which cannot be achieved using traditional FA. Good correlations can be observed between the leaking regions on FA and the abnormal regions on the vessel diameter map and the vessel complexity map on the OMAG images. This observation may be clinically important because if true, then quantitative OMAG may present the ability not only to detect the same vascular abnormalities as FA, but also to provide information about the areas of potential leakage without the need for FA. Additional patients are needed to validate this observation.
Figure 10 shows the side-by-side comparison of a normal subject and the MacTel2 subject. Subtle differences can be noticed in the vessel density maps as the MacTel2 case shows a lower global VAD. Distinct differences can be found on both vessel diameter maps and vessel complexity maps. The MacTel2 case also shows a significantly increased value in the juxtafoveal temporal and nasal regions on the vessel diameter map when compared with the normal case and an extensively decreased vessel complexity value when compared with the normal eye in the same regions of interest.
Moreover, our five-index quantitative analysis may be useful in the treatment monitoring. Figure 11 shows results from a patient who received monthly injections of a vascular endothelial growth factor (VEGF) inhibitor from January 2015 to March 2015. OCTA scans were performed at each visit with the first scan before treatment and the subsequent two scans performed one month after each treatment. OCTA volume data were segmented into superficial retinal layer, deep retinal layer, and outer retinal layer. En face images were generated for each layer using maximum projection for all three scans. A1-C1 are the superficial retinal layer, deep retinal layer, and outer retinal layer OCT angiograms, respectively, for the first visit before treatment; A2-C2 are for the second visit one month after the first treatment; and A3-C3 are for the third visit one month after the second treatment. For each visit, the deep retinal layer angiograms were chosen for the quantification of vessel density maps, vessel diameter maps, and vessel complexity maps, which are shown in D1-D3, E1-E3, and F1-F3, respectively. By observing these quantitative maps after the first injection, we found that the VDI decreased between the first and second visits, and the VCI has increased. In order to test whether our quantitative analysis has the ability to monitor the changes in vascular abnormalities after treatment, we performed quadrant-specific statistical analysis on this case for all three visits. A 2-mm circle was selected for statistical analysis as showed in Figure 12, and superior, inferior, temporal, and nasal quadrants are labeled in different colors. In each quadrant, all data points were equally divided into ten subsets and the averaged values were used to represent each subset. Therefore, we have ten samples of each quantitative index for each quadrant. The averaged value of each quadrant for each index are shown in Table 6. ANOVA was performed and the VDI showed significant decrease in superior, temporal, and inferior quadrants after treatment while VCI were found significantly increased in temporal and inferior quadrants after treatment. No significant differences were found in any quadrants for VAD, VSD, or VPI. FA was only performed at the first visit therefore we do not have all three corresponding FA images to compare with OCTA images. While these results are intriguing, a larger study is clearly needed to fully validate the usefulness of our quantitative analysis for monitoring treatment response after antiVEGF therapy.
Quantitative analysis of a MacTel2 case receiving anti-VEGF treatment.
p<0.05 in ANOVA test.
Building Database for Normal Subjects
Looking forward, we are in the process of building a normal OCTA database for quantification. We have currently enrolled five normal subjects (age matched) and obtained montage OCTA scans from each of them. OMAG en face images have been produced and used for quantification. Descriptive statistics of normal subjects are shown in Table 7. With the establishment of a normal database, future studies that compare diseased cases with normal cases or that determine disease progression or classification would be feasible.
Statistics of quantitative indices for normal subjects.
|Mean||Standard deviation||Lower 95%||Upper 95%|
In this report, we proposed five metrics, VAD, VSD, VDI, VPI, and VCI, for the quantitative assessment of OCT angiograms. There are several important features in our analytical approach. First, the quantitative OCTA analysis using the five metrics is highly repeatable as the coefficients of variance are all smaller than 0.031, suggesting that the quantitative assessment at follow-up and when used for therapeutic monitoring should provide a high level of accuracy when using OCTA. We showcased repeatability using both a normal case and a BRVO case, with different FOV and scanning location. Second, the quantitative metrics are useful to provide clinicians and physicians with a rapid and comprehensive strategy for interpreting OCT angiograms. Unlike previous studies that only provided a single quantitative index such as VAD or skeleton density, our proposed quantitative analysis includes five indices, e.g., VAD, VSD, VDI, VPI, and VCI. Along with three quantitative maps, the vessel density map, the vessel diameter map, and the vessel complexity map, the information about vascular abnormalities could be visualized and extracted from OCTA images based on multiple perspectives. Vessel length, size, area, perimeter, and morphological complexity can be made all available in one step of image analysis, potentially helping clinicians understand OCTA images better. For example, in the case that we studied here, the five metrics showed its promise to distinguish BRVO from NPDR. Both of these two cases demonstrate VAD decrease. However, for BRVO, the decrease of VAD is dominated by skeleton decline; whereas for NPDR, both VDI and VSD are increased. Third, the quantitative metrics can provide additional information when presented alongside with the qualitative OCTA images.
One commonly acknowledged shortcoming of OCTA is that it is difficult to show hyperpermeability or leakage from abnormal vasculature as seen with FA. This difference between OCTA and FA originates from the fact that OCTA detects the movement of RBCs while FA detects the presence and leakage of a small molecular weight dye that remains mostly unbound to protein and capable of diffusing out from abnormal vessels. Though fluorescein leakage registers as an important marker of a vascular abnormality, it also prevents detailed visualization of microvasculature. In this study, we showed that with our quantitative maps, OMAG might also detect the regions of abnormal microvasculature that might correspond to areas of dye leakage on FA images. This abnormal vasculature possesses many abnormal characteristics that are probably associated with increased permeability, such as vascular dilation, a decrease of vascular complexity, and capillary drop out. While it is the current dogma that OCTA is not able to detect vascular hyper-permeability, this lack of leakage detection may be mitigated to some extent by our quantitative analysis because our results suggest that the regions of vascular dilation or decreasing complexity correlate well with leakages on FA. However, this finding is based on a limited number of clinical cases. Clearly, additional studies with larger sample sizes are needed to confirm this finding.
Moreover, our quantitative analysis is versatile in terms of selecting the ROIs. In this study, the analysis was reported on the entire image. Our algorithm also allows users to freely select only part of the OCT angiograms as a region of interest with a variety of useful choices, such as rectangle, circle, ring, or ellipse regions as shown in Figure 8. This feature may be useful when user wants to compare the same ROI of different scans as the scanned regions might vary. Furthermore, we have also demonstrated the software’s ability in conducting quadratic analysis in both a BRVO case and a MacTel2 case. Even though previous studies have shown the quantification of vessel density in normal and DR human subjects,16,17 to the best of our knowledge, this is the first study that vessel diameter and vessel complexity have been quantified on OCTA images with capillary level resolution. Our three quantitative maps together are useful for providing an integrated analysis on various ocular diseases. In addition, it should also be possible to combine all our quantitative maps together to yield an integrated “vascular abnormality indicator map”. In order to do so, we would need to complete the normal database, upon which to combine the vessel density map, vessel diameter map, and vessel complexity map together, with appropriate weighting factors if necessary, to serve as a quick indicator to locate a region of vascular abnormality from patients’ scans.
There are several limitations in our quantitative OMAG technology. First, the OCTA data shows projection artifacts in the deep retinal layers on both cross-sectional images and en face images. When photons interact with moving particles, the phase of the light wave changes as a result of the Doppler effect. Therefore, when the light travels through blood vessels in living tissue, only part of the light would reflect back during its first interaction with RBCs while the majority of the photons carry on that are influenced by their first interaction. As a result, OCTA signals in the deeper layers also carry the shadows of blood vessels from the superficial plexus. These projection artifacts from the retinal vasculature projected onto deeper layers of the retina are ubiquitous and difficult to notice unless the OCTA data is segmented into superficial and deep plexus and the superficial signals are removed from the deeper flow signals. The effects of the projection tails in retinal vasculature on quantitative analysis remain underexplored, but the inherent physical characteristics of these artifacts suggest that the effects should be consistent. As a consequence, we did not perform any correction on the projection tails in this study. Zhang et al.24 recently proposed a practical and automatic projection artifact removal algorithm, which has been demonstrated to successfully remove these projection artifacts efficiently for visualizing macular neovascularization. Work is currently in progress to incorporate this projection artifact removal algorithm into our quantification assessment in order to provide more accurate quantification of the retinal microvasculature in the deep retinal layers.
The second limitation of this study is the inherent trade-off between scanning time and sampling density that makes high resolution wide-field OCTA scans more time-consuming. Although there is also a single scanning protocol available with a lower sampling spacing between adjacent A-scans, we chose the montage scan mode over the single scan mode because the montage scan provides a much higher sampling spacing () and would not lose detailed information of the microvasculature. Last, we have not performed a reproducibility study on different OCT systems with the same subject. Since some of our quantitative indices are pixel based values, the comparison of different images will require the same pixel size.
In conclusion, we have demonstrated a comprehensive quantitative OCTA analysis package for normal and diseased eyes. Collectively, the VAD, VSD, VPI, VSD, and VCI can detect abnormal microvasculature and also provide repeatable quantitative results. Our vessel density map, vessel diameter map, and vessel complexity map jointly could be useful for detecting regions with vascular abnormalities that could also leak fluorescein on FA images. Overall, quantitative OCTA has promising future applications in research and clinical management of ocular diseases.
Research supported in part by grants from National Eye Institute (R01EY024158), Carl Zeiss Meditec, Inc. (Dublin, CA), an unrestricted grant from the Research to Prevent Blindness, Inc., New York, NY.
E. Moult et al., “Ultrahigh-speed swept-source OCT angiography in exudative AMD,” Ophthalmic Surg. Lasers Imaging Retina 45(6), 496–505 (2014).http://dx.doi.org/10.3928/23258160-20141118-03Google Scholar
R. F. Spaide, “Volume-rendered optical coherence tomography of diabetic retinopathy pilot study,” Am. J. Ophthalmol. 160(6), 1200–1210 (2015).AJOPAA0002-9394http://dx.doi.org/10.1016/j.ajo.2015.09.010Google Scholar
A. H. Kashani et al., “Optical coherence tomography angiography of retinal venous occlusion,” Retina 35(11), 2323–2331 (2015).RETIDX0275-004Xhttp://dx.doi.org/10.1097/IAE.0000000000000811Google Scholar
Q. Zhang et al., “Swept source optical coherence tomography angiography of neovascular macular telangiectasia type 2,” Retina 35(11), 2285–2299 (2015).RETIDX0275-004Xhttp://dx.doi.org/10.1097/IAE.0000000000000840Google Scholar
A. Mariampillai et 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
Y. Jia et al., “Split-spectrum amplitude-decorrelation angiography with optical coherence tomography,” Opt. Express 20(4), 4710–4725 (2012).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.20.004710Google Scholar
B. Braaf et al., “Angiography of the retina and the choroid with phase-resolved OCT using interval-optimized backstitched B-scans,” Opt. Express 20(18), 20516–20534 (2012).OPEXFF1094-4087http://dx.doi.org/10.1364/OE.20.020516Google Scholar
A. Zhang et al., “Methods and algorithms for optical coherence tomography-based angiography: a review and comparison,” J. Biomed. Opt. 20(10), 100901 (2015).JBOPFO1083-3668http://dx.doi.org/10.1117/1.JBO.20.10.100901Google Scholar
R. K. Wang et 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
C. L. Chen et al., “Optic nerve head perfusion in normal eyes and eyes with glaucoma using optical coherence tomography-based microangiography,” Quant. Imaging. Med. Surg. 6(2), 125–133 (2016).http://dx.doi.org/10.21037/qims.2016.03.05Google Scholar
R. Reif et al., “Quantifying optical microangiography images obtained from a spectral domain optical coherence tomography system,” J. Biomed. Imaging 2012, 9 (2012).http://dx.doi.org/10.1155/2012/509783Google Scholar
Y. Jia et al., “Quantitative optical coherence tomography angiography of vascular abnormalities in the living human eye,” Proc. Natl. Acad. Sci. U. S. A. 112(18), E2395–E2402 (2015).http://dx.doi.org/10.1073/pnas.1500185112Google Scholar
S. A. Agemy et al., “Retinal vascular perfusion density mapping using optical coherence tomography angiography in normals and diabetic retinopathy patients,” Retina 35(11), 2353–2363 (2015).RETIDX0275-004Xhttp://dx.doi.org/10.1097/IAE.0000000000000862Google Scholar
Q. Zhang et al., “Wide-field imaging of retinal vasculature using optical coherence tomography-based microangiography provided by motion tracking,” J. Biomed. Opt. 20(6), 066008 (2015).JBOPFO1083-3668http://dx.doi.org/10.1117/1.JBO.20.6.066008Google Scholar
X. Yin, J. R. Chao and R. K. Wang, “User-guided segmentation for volumetric retinal optical coherence tomography images,” J. Biomed. Opt. 19(8), 086020 (2014).JBOPFO1083-3668http://dx.doi.org/10.1117/1.JBO.19.8.086020Google Scholar
A. F. Frangi et al., “Multiscale vessel enhancement filtering,” in Medical Image Computing and Computer-Assisted Interventation—MICCAI’98, pp. 130–137, Springer Verlag, Berlin, Germany (1998).Google Scholar
M. R. Thorell et al., “Swept-source OCT angiography of macular telangiectasia type 2,” Ophthalmic Surg. Lasers Imaging Retina 45(5), 369–380 (2014).http://dx.doi.org/10.3928/23258160-20140909-06Google Scholar