1.

Introduction

Optical coherence tomography angiography (OCTA) was recently introduced for the imaging of microvascular networks in the human eye.¹ 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) ² 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,⁷ and macular telangiectasia type 2 (MacTel2).⁸ 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)¹ 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,⁹ split-spectrum amplitude decorrelation angiography,¹⁰ and phase variance,¹¹ 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.¹² Consequently, OMAG is being established as a valuable tool for the investigation of ocular diseases both qualitatively ¹³ and quantitatively.¹⁴

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.¹⁶ 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.

2.

Method

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 $\mathrm{\Delta }\lambda$ and a 68-KHz A-scan rate. The axial resolution is $\sim 7\mu \mathrm{m}$ while lateral resolution is $\sim 15\mu \mathrm{m}$. The scanning protocol in this system provided a field of view (FOV) of $2.4\mathrm{mm}\times 2.4\mathrm{mm}$, 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™¹⁸ 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 $6.72\mathrm{mm}\times 6.72\mathrm{mm}$, 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 $\mathrm{\Delta }\lambda$ and 100-kHz A-scan rate. The axial resolution is $\sim 5\mu \mathrm{m}$ while lateral resolution is $\sim 14\mu \mathrm{m}$. SS-OCT provided a FOV of $3\mathrm{mm}\times 3\mathrm{mm}$ 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.¹⁹ This image is processed into a binary image in MATLAB (R2015b, MathWorks, Inc.) using a global threshold, hessian filter,²⁰ and adaptive threshold.¹⁵ 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).

Fig. 1

Overview of the quantitative OMAG algorithm.

Fig. 2

Representative OMAG images to illustrate the quantitative analysis algorithm. (a) Original en face OMAG image. (b) Vessel area map, a binarized vasculature image using hessian filter and adaptive threshold. This image is used for VAD, VDI, and VCI quantification. (c) Vessel skeleton map, which is obtained by iteratively deleting the pixels in the outer boundary of the vessel area map until one pixel remained along the width direction of the vessels. This image is used for VSD and VDI quantification. (d)Vessel perimeter map, which is obtain by detecting the edge of vessels in the vessel area map and deleting pixels that are not on the edge of vessels. This image is used for VPI and VCI quantification.

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 maps

VSD is calculated as the ratio of the length occupied by the blood vessels to the total area in the skeletonized vessel map

VDI is calculated using both the vessel area map and the skeletonized vessel map to yield the averaged vessel caliber as follows:

VPI is calculated using the vessel perimeter map as the ratio of the vessel perimeter to the total area of the OCTA image

VCI is calculated using both the vessel perimeter map and the vessel area map, defined as

Fig. 3

Demonstration of VCI quantification. (a)–(j) are area images (a), (c), (e), (g), and (i) and perimeter images (b), (d), (f), (h), and (j), with increasing visual complexity. Plot (k) shows their corresponding quantification of VCI.

3.

Results

Five illustrative cases are shown to demonstrate the usefulness of all the quantitative metrics that we propose.

3.1.

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,¹⁹ 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.²² 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.

Fig. 4

Quantitative OMAG on a normal case. (a) OMAG en face image with montage scanning protocol, nine smaller cubes are stitched together to achieve a $6.72\mathrm{mm}\times 6.72\mathrm{mm}$ FOV. (b) Illustration of segmentation and color code strategy. Retinal vasculature was segmented into superficial slab, deep slab, and the avascular outer retina slab. Red represents superficial plexus, green represents deep plexus, and blue for outer retinal slab. All three slabs were integrated together and represented in a color image in (a), excluding retinal nerve fiber layer. (c) Vessel density map. This map shows the VAD of image (a), calculated using a kernel moving across the entire OMAG image. (d) Vessel diameter map. This map shows the local averaged VDI using the same calculation method as (c). (e) Vessel complexity map. This map shows the local averaged VCI using the same calculation method as (c). (f)–(h) Integration of quantitative maps (c)–(e) with binary vessel area map [Fig. 2(b)].

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)].

Table 1

Quantitative analysis of a normal case and a NPDR case.

Indices	VDI	VAD	VSD	VPI	VCI	Flow impairment zone (mm2)
Normal	24.421	0.518	0.146	0.349	17838.540	0.264
DR	24.665	0.370	0.101	0.245	12954.770	3.450

3.2.

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 $3\text{pixel}\times 3\text{pixel}$) 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.

Fig. 5

Quantitative OMAG on a NPDR case. (a) FA image. (b) Whole retina OMAG image, correspond to the red rectangular region in the FA image. Scanning area is $6.72\mathrm{mm}\times 6.72\mathrm{mm}$. (c) Flow impairment zone quantification. Capillary drop out region was detected from vessel skeleton image and presented in green color, superimposed on (b). (d) Vessel density map. (e) Vessel diameter map. (f) Vessel complexity map. (g)–(i) Integration of quantitative maps (d)–(f) with binary vessel area map.

3.3.

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.

Fig. 6

Quantitative OMAG analysis of a BRVO case with repeated scans. (a) and (d) OMAG image of occluded region, covering a FOV of $4.3\mathrm{mm}\times 4.6\mathrm{mm}$. (b) and (e) Vessel area map. (c) and (f) Vessel skeleton map. (g) and (j) Vessel density map integrated with vessel area map. (h) and (k) Vessel diameter map integrated with vessel area map. (i) and (l) Vessel complexity map integrated with vessel area map.

Table 2

Quantitative analysis and repeatability test of a BRVO case for two repeated scans.

	VDI	VAD	VSD	VPI	VCI
Scan 1	35.300	0.449	0.175	0.395	2862.719
Scan 2	34.196	0.451	0.182	0.393	2827.102
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.

Fig. 7

Repeatability of Quantitative OMAG on a normal case. (a) Whole retina OMAG image of a normal case, covering an area of $6.72\mathrm{mm}\times 6.72\mathrm{mm}$. (e) Repeated scan of the same subject at same location, whole retina OMAG image. (b) and (f) Vessel density map integrated with the binary vessel area map. (c) and (g) Vessel diameter map integrated with the binary vessel area map. (d) and (h) Vessel complexity map integrated with the binary vessel area map.

Table 3

Repeatability test of a normal case.

	VDI	VAD	VSD	VPI	VCI
Scan 1	23.206	0.511	0.151	0.353	19282.771
Scan 2	23.619	0.512	0.149	0.351	18979.230
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.

Fig. 8

(a) Illustration of region of interest selection. An original OCTA image with a FOV of $6.72\mathrm{mm}\times 6.72\mathrm{mm}$. (b) Illustration of quadratic analysis. Red represents superior quadrant, green represents nasal quadrant, blue represents in inferior quadrant, and yellow represents temporal quadrant. (c) ROI selection example of a circle. (d) ROI selection example of a rectangle. (e) ROI selection example of an ellipse. (f) ROI selection example of a circle.

Table 4

Quantitative analysis of a BRVO case with different choice of region of interest.

	VAD	VSD	VPI	VDI	VCI	ROI size (mm×mm)
Whole image	0.493	0.200	0.419	32.554	7234.772	44.890
Circle	0.491	0.195	0.418	33.272	396.946	10.521
Rectangle	0.491	0.193	0.415	33.633	457.660	11.389
Ellipse	0.519	0.220	0.454	31.28575	62.690	3.964
Ring	0.494	0.193	0.416	33.812	285.845	8.987

Table 5

Quadratic quantitative analysis of a BRVO case.

	VAD	VSD	VPI	VDI	VCI	ROI size (mm×mm)
Superior	0.493	0.197	0.411	33.942	239.2756	6.0934
Inferior	0.473	0.202	0.424	32.034	257.6289	6.0934
Nasal	0.478	0.200	0.424	32.838	251.9077	6.0934
Temporal	0.489	0.194	0.413	34.292	237.8684	6.0934
Whole image	0.493	0.200	0.419	32.554	7234.772	44.890

3.4.

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 plexus²³ 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.

Fig. 9

Quantitative OMAG analysis of a MacTel2 case. (a) FA image of the MacTel2 subject. (b) SS-OCT OMAG image of the whole retina, covering $3\mathrm{mm}\times 3\mathrm{mm}$ area, corresponding to the white rectangular region in FA image. (c) Vessel density map integrated with vessel area map on the whole retina. (d) Vessel diameter map integrated with vessel area map on the whole retina. (e) Vessel complexity map integrated with vessel area map on the whole retina. (f) OMAG image of the deep capillary plexus. (g) Vessel density map integrated with vessel area map of the deep capillary plexus. (d) Vessel diameter map integrated with vessel area map of the deep capillary plexus. (e) Vessel complexity map integrated with vessel area map of the deep capillary plexus.

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.

Fig. 10

Quantitative OMAG comparison between the MacTel2 and normal cases. (a) $3\mathrm{mm}\times 3\mathrm{mm}$ OMAG image of the whole retina of a MacTel2 subject. (b) Vessel density map. (c) Vessel diameter map. (d) Vessel complexity map. (e) $3\mathrm{mm}\times 3\mathrm{mm}$ OMAG image of the whole retina of a normal subject. (f) Vessel density map. (g) Vessel diameter map. (h) Vessel complexity map.

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.

Fig. 11

Three OCTA scans of a patient with MacTel2 monthly apart. a1–a3 shows the superficial retina for the first, second, and third scans, respectively. b1–b3 shows the deep retina for three scans, respectively. c1–c3 shows the outer retina for three scans, respectively. d1–d3 shows the vessel density map for three scans, respectively and e1–e3 shows the vessel diameter maps for three visits while f1–f3 shows the vessel complexity maps for all three visits.

Fig. 12

Illustration of the quadratic analysis. 2-mm circle is selected excluding the FAZ. Superior, inferior, temporal, and nasal quadrants are labeled in different colors as showed in text.

Table 6

Quantitative analysis of a MacTel2 case receiving anti-VEGF treatment.

Indices	Quadrant	1/9/2015	2/5/2015	3/31/2015
VDI	Superiora	21.397	19.704	18.995
	Temporala	20.228	18.857	18.063
	Inferiora	19.901	19.082	18.268
	Nasal	19.441	18.916	19.228
VAD	Superior	0.391	0.368	0.341
	Temporal	0.364	0.322	0.340
	Inferior	0.401	0.378	0.395
	Nasal	0.388	0.374	0.361
VCI	Superior	4.615	4.952	4.852
	Temporala	4.840	4.726	5.307
	Inferiora	5.313	5.378	5.830
	Nasal	5.112	5.371	5.329
VSD	Superior	0.097	0.095	0.092
	Temporal	0.092	0.086	0.096
	Inferior	0.105	0.103	0.107
	Nasal	0.097	0.098	0.094
VPI	Superior	0.225	0.221	0.213
	Temporal	0.236	0.221	0.241
	Inferior	0.264	0.257	0.271
	Nasal	0.250	0.254	0.248

^ap<0.05 in ANOVA test.

4.

Discussion

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.²⁴ 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 $6\mathrm{mm}\times 6\mathrm{mm}$ 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 ($\sim 10\mu \mathrm{m}$) 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.