The disruption of mechanisms that regulate perfusion of ocular tissues has been recognized as an important contributor to the pathogenesis and development of a variety of ocular diseases.1 Such diseases include: glaucoma,23.4.–5 age-related macular degeneration,6 and diabetic retinopathy.7 Over the past several decades, various imaging techniques, such as fluorescein angiography,8 laser Doppler flowmeter (LDF),2,9,10 scanning laser Doppler flowmetry (SLDF),1184.108.40.206.–16 laser speckle flowgraphy (LSFG),17,18 and Doppler optical coherence tomography (D-OCT),19 have been used to directly measure the ocular blood flow and vascular volume, or to indirectly quantify the parameters via perfusion time, flow velocity, flux of the blood cells, or vessel diameter changes.
Although several studies have shown that these techniques are capable of quantifying ocular tissue perfusion, large variations in their measurements limit their application and use in ophthalmology. With SLDF [Heidelberg retina flowmeter (HRF), Heidelberg Engineering, Heidelberg, Germany], the coefficients of variation (CVs) in the literature range from 4% to 28%.11,20,21 Pemp et al. used an LDF to measure the optic nerve head (ONH) blood flow by detecting a frequency shift of the scattered light according to the optical Doppler effect,9 and found that the CVs were 5.9% for normal subjects and 8.6% for glaucoma subjects.22 A later study from the same group measured the ocular hemodynamic parameters from 16 healthy subjects using SLDF, the blue field entoptic technique, LDF, laser Doppler velocimetry, retinal vessel analyzer, and color Doppler imaging. The variability from five consecutive measurements was analyzed with each technique; it was found that the CVs of the hemodynamic parameters were between 1.6% and 18.5%.23 Using LSFG (Kyushu Institute of Technology, Iizuka, Japan), Yaoeda et al. investigated the microcirculation in ONH from 60 normal volunteers by measuring the square blur rate, volume, flow, and velocity. They reported that the CVs for intervisit repeatability of the parameters ranged from 9.7% to 19.8%.17,18
We previously developed a noninvasive imaging technique capable of providing three-dimensional in vivo blood flow visualization within microcirculatory tissue beds in the eye. The imaging approach uses a novel optical microangiography (OMAG) technique based on Fourier-domain optical coherence tomography (FD-OCT) systems.24220.127.116.11.–29 By detecting the differences in the scattered light from moving particles, e.g., the red blood cells (RBCs) in vessels, OMAG enables the visualization of microcirculation within tissue beds in the eye.2829.–30 Results with the OMAG technique show promise in the visualization of vessel networks and neovascularization in disease states.24,26,27,29,30 Furthermore, we have shown that the OMAG signal is proportional to the concentration of blood flow, which can be referred to as blood flux, i.e., the number of blood cells passing through a vessel cross-sectional area per unit time,31,32 and therefore, the blood flux signal permits OMAG to quantify retinal blood. The purpose of this study was to evaluate the intravisit and intervisit repeatability as well as the interobserver reproducibility of ONH perfusion measurements using the OMAG technique.
This study was approved by the Institutional Review Board of the University of Washington (UW) and informed consent was obtained from all subjects. This study followed the tenets of the Declaration of Helsinki and was conducted in compliance with the Health Insurance Portability and Accountability Act.
Volunteers with normal optic discs were prospectively recruited at the UW Eye Institute. Inclusion criteria were best-corrected visual acuity of or better, refractive error between and spherical equivalent, and a normal retinal nerve fiber layer (RNFL) thickness on FD-OCT. Exclusion criteria were significant media opacity preventing high-quality imaging, previous diagnosis of migraine, or any ocular disease. The right eye from each subject was imaged and included in this study.
Image Acquisition and Scanning Protocol
All eyes were scanned using a 68-kHz Cirrus HD-OCT 5000-based OMAG prototype spectral-domain (SD-) OCT system with active motion-tracking capability (Carl Zeiss Meditec Inc., Dublin, California) to minimize motion artifacts.33 The center wavelength of the OCT system worked at 840 nm with a typical 3-dB bandwidth of 45 nm, axial resolution of , and lateral resolution of . Three OMAG scans centered on the ONH were acquired at each visit. Between each scan, patients were asked to sit back and rest for about 30 s to 1 min. Participants returned for a second visit within a 6-week time frame, at which time the same scanning protocol was repeated by the same operator.
The OMAG scan pattern generated a volumetric dataset over a area centered on the optic disc. Each single B-scan was composed of 245 A-scans. For each A-scan, 1024 sampling points were collected along a 2.0-mm axial scan depth. Four consecutive B-scans were acquired at each fixed transverse location (i.e., along the fast axis) before the scanning probe moved to the next transverse location on the optic disc. A total of 245 transverse locations were sampled, located apart. The time difference between two successive B-scans was , which corresponds to a B-scan acquisition rate of .
All acquired volumetric scans were processed with an OMAG algorithm to extract both the structural and blood flow signals. The details of the OMAG algorithm have been described elsewhere.28,29 Briefly, blood flow signals were calculated by subtracting complex OCT signals between each consecutive B-scan pair at each location.28,29 The intensity of the blood flow signal is proportional to the amount of RBCs passing through the vessels.31,32 The average of the differences at each location was calculated to generate the final blood flow signal, , as described in
Retinal Layer Segmentation and Vascular En Face Images
Using the structural images, boundaries of the inner limiting membrane (ILM), the anterior surface of lamina cribrosa (aLC), and the visible posterior surface of lamina cribrosa (pLC) were delineated using a laboratory-developed semiautomatic retinal layer segmentation program.34 The segmentation software automatically detected the boundaries of the retinal layers from the structural OCT cross-sectional images by measuring the gradient of OCT signals.34 When the changes of intensity were too indistinct for the algorithm to detect a boundary, the operator could interrupt the automatic algorithm and manually guide the system to identify the correct boundaries. Three retinal layers were defined based on the segmented retinal boundaries: the prelaminar tissue [preLC, between ILM and aLC, Fig. 1(a)], lamina cribrosa (LC, between aLC and pLC), and full ONH (between ILM and pLC). Since vascular images were derived from the structural images, the boundaries detected based on the structural image could then be applied to OMAG vascular images [Fig. 1(b)] to provide precise correspondence between structural and OMAG retinal layer locations.
Retinal layer segmentation was performed on all the images by a single operator [Chieh-Li Chen (CLC)]. To further test the interobserver reproducibility, two independent operators (CLC and CX) performed the layer segmentation in a masked manner with the same program on one randomly selected scan for each eye.
Maximum intensity projection analyses, which detect the blood flow signal with the highest intensity value along each A-scan, were performed within each segmented retinal layer to generate vascular en face images of preLC, LC, and the full ONH. The boundary of the optic disc margin to be used for later analysis was manually delineated by identifying the edge of Bruch’s membrane on the structural B-scans.
Flux, Vessel Area Density, and Normalized Flux Measurements
To quantify ONH perfusion, we created three quantitative measurements: flux, vessel area density, and normalized flux.35
The OMAG technique measured the variations in OCT signals between consecutive B-scans by calculating the differences in complex OCT signals. The differences in the complex OCT signals contained both differences from phase and amplitude information. The phase differences resulted from the motion of RBCs flowing through the vessels, whereas the amplitude differences came from the variation of the scattered OCT signal and was related to the number of particles passing through the vessels. Combining both information of flow velocity and the number of particles (concentration), the intensity of flow signal is an indicator of flux (the number of particles passing through the cross section per unit time), and we first measured the flux using
Vessel area density was the ratio of vessel-containing area within the ONH to the total ONH area [Eq. (3)]. Vascular en face images of the blood vessels in the ONH were detected by a multiscale Hessian filter developed based on the method of Frangi et al.36 The Hessian filter measures the local curvature by calculating the second-order partial derivatives of local structure from the vascular en face images. Based on the calculated two-dimensional eigenvalues and , where , the structural information was converted into probability-like estimates of vessel presence36
As less nerve head tissue requires less ONH perfusion, we further calculated the normalized flux to avoid bias from the reduced vessel area, e.g., in the situation of comparing ONH perfusion between glaucomatous and normal eyes. In Eq. (4), the flow signals accumulated from the vessels were divided by the vessel area, indicating the evaluation of flux in the vessels only:
Demographic and clinical characteristics were summarized by standard descriptive statistics. For each eye, two CVs were calculated from three scans (acquired at the same visit) and six scans (from two visits) to evaluate the intravisit and intervisit repeatability, respectively, for all three ONH perfusion metrics. Mean CVs and standard deviations of 10 subjects were further computed to represent the final intravisit and intervisit repeatability for each ONH perfusion metric. To evaluate agreement between the observers’ measurements for interobserver reproducibility, the intraclass correlation coefficient (ICC) two-way mixed model with their 95% confidence interval on absolute agreement was used to analyze measurement reliability.37 The values of the ICC can range from 0 to 1, with a higher value indicating better reliability. ICC less than 0.40 was considered as poor, values between 0.4 and 0.6 as fair, values between 0.61 and 0.8 as good, and values higher than 0.8 as excellent.38 IBM SPSS Statistics for Windows, Version 23.0 (IBM Corp., Armonk, New York) was used for statistical analyses.
Ten eyes from 10 healthy volunteers (five men and five women) were recruited at the UW Eye Institute. The mean age was -years old (range 28.8 to 50.0-years old). The mean RNFL thickness of the normal subjects was . All the scans were included in the analysis.
Figure 2 shows an example of the vascular en face images and vessel maps detected using a Hessian filter of preLC, LC, and full ONH from a 28-year-old healthy woman. Three scans acquired at the same visit are displayed in adjacent images. Excellent intravisit repeatability of vascular networks in the ONH can be observed. In addition to the central retinal arteries and veins, the capillaries inside and surrounding the ONH were also visualized.
Analysis of the measurements of three ONH perfusion metrics evaluating intravisit repeatability and the corresponding CVs for the three retinal layers are listed in Table 1. In all retinal layers, vessel area density showed the best intravisit repeatability, followed by normalized flux, and finally by flux.
Intravisit repeatability of ONH perfusion metrics (flux, vessel area density, and normalized flux) using OCTA in normal eyes.
|Flux||CV (%)||Vessel area density||CV (%)||Normalized flux||CV (%)|
Note: CV, coefficient of variation. ONH perfusion metrics and CVs are presented as mean±standard deviation.
Table 2 summarizes the intervisit repeatability for the three ONH perfusion metrics. In preLC and full ONH, vessel area density demonstrated the highest intervisit repeatability, followed by normalized flux, and then by flux, whereas in the LC, normalized flux showed the highest intervisit repeatability, followed by flux; vessel area density had the lowest intervisit repeatability.
Intervisit repeatability of ONH perfusion metrics (flux, vessel area density, and normalized flux) using OCTA in normal eyes.
|Flux||CV (%)||Vessel area density||CV (%)||Normalized flux||CV (%)|
Note: CV, coefficient of variation. ONH perfusion metrics and CVs are presented as mean±standard deviation.
Interobserver reproducibility of ONH perfusion measurements was analyzed as ICC and its 95th percentile confidence intervals, which are presented in Table 3. For the interobserver reproducibility, all three ONH perfusion metrics showed excellent agreement between observers in each of the three retinal layers: all measurements were above 0.950.38
Interobserver reproducibility of ONH perfusion metrics (flux, vessel area density, and normalized flux) using OMAG in normal eyes.
|Flux||ICC||Vessel area density||ICC||Normalized flux||ICC|
|(0.975, 0.998)||(0.951, 0.997)||(0.978, 0.999)|
|(0.965, 0.998)||(0.973, 0.998)||(0.871, 0.991)|
|(0.998, 1.000)||(0.972, 0.998)||(0.971, 1.000)|
Note: ICC, intraclass correlation coefficient. ONH perfusion metrics are presented as mean±standard deviation. For ICC, 95% confidence interval is shown in the parentheses.
We demonstrate that OMAG enables visualization of the microcirculation in the ONH and allows the quantification of ONH perfusion in normal eyes with high repeatability and reproducibility. By detecting the differences in the backscattered OCT signals caused by RBCs between adjacent cross-sectional frames at the same location, OMAG is able to generate blood flow signals with high resolution in a noninvasive fashion. Since the blood flow signals are proportional to the number of RBCs passing through the vessels,31,32 the quantification of ONH perfusion becomes possible.
Various imaging techniques have been developed to measure and quantify the ONH perfusion. However, both repeatability and reproducibility have been considered problematic with all techniques used to date thus limiting their value in clinical practice as well as in clinical trials.11,17,18,2021.22.–23 In our study, using the OCTA with OMAG technique for quantifying ONH perfusion, the intravisit and intervisit repeatability were and (CVs for the full ONH). Our results provide ONH perfusion measurements with higher repeatability and reproducibility values than those from other reports of noninvasive imaging techniques. Unlike LSFG, which detects the accumulated speckle changes and measures the blurriness of the images, or LDF, which detects the flow signal by measuring the Doppler frequency shift, OMAG-based OCTA directly measures the localized differences in consecutive OCT signals. In this way, OMAG-based OCTA is more sensitive to the microcirculation in the tissue, and, therefore, is able to generate blood flow images with higher resolution, better contrast, and, as a consequence, less variability.
Two studies have assessed the intravisit and intervisit repeatability of ONH perfusion using a split-spectrum amplitude-decorrelation angiography (SSADA) algorithm, another OCTA technique.39,40 High intravisit ( mean CV) and intervisit repeatability ( mean CV) of ONH perfusion measurements were detected using OCTA with the SSADA algorithm, further supporting and strengthening the proposal that OCTA-based methods provide reproducible and reliable measurement of ONH perfusion and may be suitable for clinical practice.
A number of factors could be responsible for the slightly higher CV values reported in our study compared to others who used SSADA-based OCTA.39,40 In the SSADA algorithm, an extra image registration algorithm was required to register 4 orthogonal raster scanned volumes and to remove the motion between B-scan positions along the slow scan axis.39 The angle and location of the acquisition can vary between both same day and different day scans. It is conceivable that the additional image registration algorithm used to register the images between scans may unintentionally introduce an apparent reduction in between-scan and between-visit variability. Our OMAG processing compensates for the variation between consecutive B-scans at the same location but not for the variation in the sampling location in the longitudinal direction. Nonetheless, our results demonstrate excellent intravisit and intervisit repeatability and interobserver reproducibility.
Another possible reason for the slightly higher CV values seen with OMAG-based OCTA is that three subjects were unable to be scanned at the same time of day at two consecutive visits. Since ONH perfusion fluctuates during the day,21 the time difference may contribute to a larger variation in our results. Additionally, differences in the methods of blood flow signal extraction and image segmentation could contribute to variability among studies.
With the use of a light source with a central wavelength at 840 nm in our prototype SD-OCT system, overlying prelaminar tissue may reduce the visibility of microcirculation in the LC. The intact prelaminar tissue in the ONH may block the light signal and thus limit our ability to visualize deep tissues. Such a blockage of the light signal may cause increased variation in ONH perfusion measurements in the LC, especially for vessel area density assessment. The limited penetration issue might be resolved by using a light source with a longer center wavelength, or using swept-source OCT at 1060-nm wavelength.41,42
One limitation of this study is that we involved only healthy subjects to test the repeatability and reproducibility of ONH perfusion measurements using OCTA with the OMAG technique. Since this was the first study testing the repeatability and reproducibility to quantify ONH perfusion using OMAG, we started with normal subjects to make sure that the flow measurements were not affected by disease progression. Other challenges of having diseased eyes in the repeatability and reproducibility study are that the repeatability and reproducibility may be affected by the larger variations in the flow signal in the defect areas, or receiving treatment. Nevertheless, further investigation with diseased eyes, such as glaucoma or diabetic retinopathy, to reflect the actual clinical situation and to provide better estimation is required.
In conclusion, OMAG provides a noninvasive method to visualize and quantify ONH perfusion in human eyes with excellent repeatability and reproducibility. Our results show that OMAG may add additional insights into ONH perfusion. Such insights also have the potential to provide a better understanding of the mechanisms and etiology of perfusion abnormalities in eye diseases, monitoring of their progression and also improve decisions in managing their care.
This study was supported in part by the National Institutes of Health contract NEI R01-EY024158, Carl Zeiss Meditec Inc., and Research to Prevent Blindness (New York, New York). C.L.C., Q.Z., and R.K.W. received research support from Carl Zeiss Meditec, Inc. M.A.J. is a consultant for Healonics, Sensimed, Ivantis, Cascade Ophthalmics, and receives royalties from intellectual property owned by Allergan. R.K.W. has significant financial interest in the intellectual property of OMAG technology, owned by Oregon Health & Science University. K.D.B., C.X., J.C.W., D.G., R.C.M., and P.P.C.: None.
Chieh-Li Chen received her BS degree in electrical engineering and her MS degree in electronics engineering from the National Taiwan University in 2008 and 2010, respectively, and her PhD in bioengineering from the University of Pittsburgh in 2014. She is a senior research fellow at the University of Washington, Department of Bioengineering. Her research interests include medical image processing and data analysis, especially for optical coherence tomography (OCT) and OCT-based angiography images.