Open Access
Articles  |   October 2020
Choroidal Biomarkers: A Repeatability and Topographical Comparison of Choroidal Thickness and Choroidal Vascularity Index in Healthy Eyes
Author Affiliations & Notes
  • Katharina Breher
    Institute for Ophthalmic Research, University of Tübingen, Tübingen, Germany
  • Louise Terry
    School of Optometry and Vision Sciences, Cardiff University, Cardiff, UK
  • Thomas Bower
    School of Engineering, Cardiff University, Cardiff, UK
  • Siegfried Wahl
    Institute for Ophthalmic Research, University of Tübingen, Tübingen, Germany
    Carl Zeiss Vision International GmbH, Aalen, Germany
Translational Vision Science & Technology October 2020, Vol.9, 8. doi:
  • Views
  • PDF
  • Share
  • Tools
    • Alerts
      This feature is available to authenticated users only.
      Sign In or Create an Account ×
    • Get Citation

      Katharina Breher, Louise Terry, Thomas Bower, Siegfried Wahl; Choroidal Biomarkers: A Repeatability and Topographical Comparison of Choroidal Thickness and Choroidal Vascularity Index in Healthy Eyes. Trans. Vis. Sci. Tech. 2020;9(11):8.

      Download citation file:

      © ARVO (1962-2015); The Authors (2016-present)

  • Supplements

Purpose: Choroidal thickness (ChT) and choroidal vascularity index (CVI) represent two important metrics in health-, disease-, and myopia-related studies. Wide-field swept-source optical coherence tomography (OCT) provides improved and extended imaging and extraction of choroidal variables. This study characterizes the topography and repeatability of these parameters in healthy eyes.

Methods: Swept-source OCT volume scans were obtained on 14 young adult patients on three separate days. ChT and CVI were automatically corrected for image magnification and extracted for different enface regions within an extended ETDRS grid of 10 mm diameter. Topographical distribution, correlation to ocular length, and intersession repeatability of both choroidal parameters were assessed.

Results: CVI showed little fluctuation between subfields, unlike ChT, which demonstrated thinning toward the peripheral choroid (coefficients of variation 5.92 vs. 0.89). ChT showed a consistent negative correlation with axial length (ρ = −0.05 to −0.61), although this was only statistically significant in the inner superior subfield (P = 0.02). There was no consistent or significant relationship between CVI and axial length or between CVI and ChT. The repeatability of CVI measurements (3.90%–5.51%) was more consistent between scan regions than ChT measurements (10.37–20.33 µm).

Conclusions: CVI values were consistent across the central 10 mm of the retina, while ChT reduced with eccentricity. The repeatability of both parameters is similar to the effect size reported in many studies using the choroid as a biomarker, which should be considered in the interpretation of findings.

Translational Relevance: This study provided normative as well as metrological information for the clinical interpretation of ChT and CVI in health and disease.

The choroid is located between the retina and the sclera and consists of five layers—Bruch's membrane, choriocapillaris, Haller's layer, Sattler's layer, and the suprachoroid—with its main task to provide blood supply to the outer retina.1,2 The choroid is suggested to play an essential role in various fields of human eye research. Changes to the choroidal structure have been suggested as biomarkers for a variety of ocular pathologies; altered choroidal thickness (ChT) was previously observed in central serous chorioretinopathy,3 Vogt-Koyanagi-Harada disease,4 and polypoidal choroidal vasculopathy and exudative age-related macular degeneration.5 Moreover, bidirectional changes of ChT influence the regulation of eye growth due to its location in the signaling pathway from retina to sclera.2,6 
However, the exact physiologic mechanisms behind the observed ChT changes still require further clarification. It is assumed that the fast alterations are based on adjustments of the vascular rather than stromal components.2,710 Therefore, choroidal vascularity developed as a further subject of interest in human eye research. The term choroidal vascularity index (CVI) was introduced as the ratio of luminal area to total choroidal area. The higher the CVI, the more vascular tissue compared with stromal tissue is present in the choroid.11 CVI was previously investigated in a healthy study population,12 as well as in ocular pathologies, such as diabetic retinopathy,13 central serous chorioretinopathy,14 age-related macula degeneration,15 panuveitis,16 and retinitis pigmentosa.17 While there are only a few studies existent, they report only vague and somewhat inconsistent association between choroidal vascularity and myopia.1821 
ChT and CVI can be obtained from optical coherence tomography (OCT) images. After segmentation of the choroidal area, the hyporeflective blood vessel lumina are distinguished from the hyperreflective stromal tissue. The reliability of the results from graphical analysis strongly depends on the OCT technology (scan depth, scan resolution, and image contrast)11,22 but also on subject factors (absolute choroidal thickness and the pigmentation of the retinal pigment epithelium).23,24 Therefore, these parameters influence the repeatability of ChT and CVI measurements. The repeatability metric helps to differentiate between measurement noise and true measured effect, which is pertinent for interpreting study findings. This becomes essential when effect sizes are small, in combination with high intersubject variability, as is typical of studies of the choroid. 
Therefore, the purpose of this study is to evaluate the repeatability of ChT and CVI measurements metrologically and to put them into the context of commonly found effect sizes in choroidal research. Moreover, the topographic distribution across an extended scan area, as well as the relationship between choroidal parameters and axial length, is assessed. To our knowledge, this is the first study for this purpose using automated analysis of widefield volume scans from swept-source OCT. 
Materials and Methods
Study Participants
The prospective and cross-sectional study adhered to the Declaration of Helsinki and was approved by the ethics committee of the Faculty of Medicine of the University of Tübingen. Written informed consent was obtained from all participants prior to data collection. Fifteen patients without known ocular pathologies were enrolled in the study. One patient was excluded from analysis due to erroneous choroidal segmentation resulting from small pupil size and thus reduced image quality. The study participants included in the analysis had a mean age of 27 ± 3 years (range, 24–37 years) and an average axial eye length (AEL) of 24.83 ± 0.91 mm (range, 23.12–26.52 mm). 
OCT Scanning Protocol
Participants underwent nine OCT volume scans on their undilated right eye, using swept-source OCT (ZEISS PlexElite 9000; Carl Zeiss Meditec, Inc., Dublin, CA, USA). The OCT device uses a central wavelength between 1040 and 1060 nm and a sweep range from 980 to 1120 nm, resulting in an increased scan depth of 3 mm. The axial and lateral resolutions in tissue are 6.3 µm and 20 µm, respectively. Three consecutive scans were performed, at the same time of day on three separate days, to control for diurnal variations in choroidal structure.25,26 Moreover, participants were advised to avert caffeine27,28 and nicotine29 at least 1 hour prior to the measurement to avoid short-term changes of ChT. In order to wash out near work-induced accommodative choroidal thinning30,31 and to stabilize blood pressure,26 participants completed a 10-minute resting phase with distance viewing. For the current study, a 12 × 12-mm isometric volume scan centered on the fovea was chosen. This scan pattern consisted of 1024 B-scans with 1024 A-scans per B-scan. Eye tracking, enhanced-depth imaging, and follow-up mode were enabled to reduce motion artifacts, to enhance the visibility of the choroid, and to ensure a constant scanning area throughout all scans. The same experienced examiner obtained all scans. Furthermore, participants moved their head off of and back onto the chinrest and headrest between each scan. 
Extraction of Choroidal Thickness
Choroidal segmentation was performed automatically by image processing of each volume scan via the Advanced Retinal Imaging Network (ARI Network; Carl Zeiss Meditec, Inc.). The ARI Network represents a research portal provided by Carl Zeiss Meditec, Inc., which offers various algorithms for researchers using the PlexElite 9000 OCT device. The choroidal segmentation algorithm uses a graph-based approach with a combination of intensity, axial gradient, gradient magnitude, and positional information to segment the retinal pigment epithelium (RPE) and the choroidal-scleral interface. To save time and increase segmentation reliability in volume scans with multiple adjacent B-scans, a subset of B-scans is segmented first. The remaining B-scans are segmented considering the results from the first subset. Exemplary central and peripheral B-scans of the scan were checked for segmentation errors by the examiner.32 
The choroidal thickness map for the entire scan area was extracted as a two-dimensional matrix in pixels and was further processed in MATLAB (MATLAB 2020a; The MathWorks, Inc., Natick, MA, USA). An axial conversion factor of 1.9531 from pixel to microns was applied, which is based upon the ratio of true scan depth (3000 µm) to digital image scan depth (1536 pixels). Subsequently, a modified extended Early Treatment of Diabetic Retinopathy Study (ETDRS) grid33 was overlaid on the segmented choroidal thickness map (Fig. 1). It consisted of four rings with baseline radii of 0.5 mm (central), 1.5 mm (inner), 3.0 mm (outer), and 5.0 mm (extended). The three outer rings were further divided into nasal, temporal, superior, and inferior sectors. The ETDRS grid size was adjusted for the AEL (equation 2) of the participant (ZEISS IOLMaster 700; Carl Zeiss Meditec AG, Jena, Germany) to compensate for the individual ocular magnification. Data from the optic nerve head (ONH) was selected and subsequently excluded via a semiautomated approach due to inherent segmentation errors in this region. The ONH was excluded by presenting the examiner with an en-face stack image and prompting the user to select with their cursor two points that define the diameter of a circle. This circle was later used to define an exclusion zone for calculations of CVI. These steps are described in more detail in the following section. 
Figure 1.
Extended ETDRS grid, showing the optic nerve head exclusion zone for a right eye. The ring diameters are 1 mm, 3 mm, 6 mm, and 10 mm.
Figure 1.
Extended ETDRS grid, showing the optic nerve head exclusion zone for a right eye. The ring diameters are 1 mm, 3 mm, 6 mm, and 10 mm.
Extraction of Choroidal Vascularity Index
The novel CVI calculation algorithm is written in MATLAB (MATLAB 2019b; The MathWorks, Inc.) and makes use of the Parallel Computing Toolbox for parallelized “for” loops with slices of the volume sent to each worker node in turn for processing. The user interface for this code is in the form of Windows dialogues and figure plots, which request cursor click positions for various tasks. The code also relies on the preparation of a data file, listing the filenames of the volume scans, along with the AEL of each scan. The scan itself is stored as an 8-bit unsigned grayscale image. The procedure for calculating CVI for an image once it is loaded into MATLAB is adapted from the protocol of Sonoda et al.34,35 and is outlined below. 
The volume scan is extracted from the input file and stored in a three-dimensional matrix. The RPE segmentation and choroidal thickness map are stored in a matrix, and the segmentation for the choroidal-scleral boundary is calculated by summing the RPE segmentation matrix and the choroidal thickness matrix. The AEL is also extracted from the input files. To give the algorithm a baseline to quantify the grayscale value of a typical fluid, the average voxel value is taken in the two non-en-face planes of the scan, at 50% of the width of the slice, and 75% of the height of the slice (assumed to be within the vitreous). This value is referred to as the “liquid norm.” A mask is then created that specifies the region of interest (ROI) for analysis. From the RPE and choroidal segmentation data, a three-dimensional binary mask matrix is generated, holding a value of 1 within the ROI and 0 outside the ROI. For normalization of the image, the liquid norm is used as a lower threshold on the volume scan (i.e., all values below the threshold were set to zero), and voxel values were rescaled between 0 and 255. To reduce speckle and assist with thresholding, a median blur was applied to each slice in the scan using a 3 × 3 window. The volume scan was binarized using Niblack auto local thresholding36 over each slice of the scan, whereby a moving circular window of radius 20 pixels specified the binarization threshold. The threshold, T, was calculated from the local mean, m, and local standard deviation, σ, of the pixels in the window, along with a parameter k, which was set to 0.2.  
\begin{equation}T = m + k\; \times \;\sigma .\end{equation}
Any value above the threshold was set to 1, and values below were set to 0, thus producing a binarized image. This study made use of an extended ETDRS grid, centered on the fovea, with an extended diameter of 10 mm and exclusion of the optic nerve head region (Fig. 1). 
Careful lateral scaling was performed for each image, to account for the effects of ocular biometry on OCT magnification.37 In the lateral dimension, the length each pixel, represented in mm pmm, was given by  
\begin{equation}{p_{mm}} = {\rm{\;}}\frac{{\left({{l_a} - {\rm{\;}}{\Delta _{cp}}} \right){\rm{\;}}{\theta _{air}}}}{{{d_p}{n_B}}},\end{equation}
where la is the AEL, Δcp is the distance between the cornea and principal plane and is given as 1.82 mm,38 θair is the scan angle (in radians) of the OCT in air and is 40° for this study, dp is the dimension of the en-face scan in pixels, and nB is the bulk refractive index of the ocular humors, which is taken as 1.336. The 13 ETDRS grid regions were stored as two-dimensional Boolean mask matrices. Each matrix was also masked against the optic nerve head exclusion circle; typically, this intersected with the extended nasal grid cell. The CVI was calculated from each voxel column of an en-face view to create a map of CVI for the scan. Within the ROI of each column, the number of voxels representing the vessel lumen was counted, nv, and the total number of voxels in the ROI column was counted, nROI. The CVI for each column was calculated as  
\begin{equation}CV{I_{col}} = {\rm{\;}}\frac{{{n_v}}}{{{n_{RoI}}}}{\rm{\;}} \times 100.\end{equation}
This CVI map was stored and presented as an image with an ETDRS grid overlay. Finally, the CVI for each of the ETDRS regions was calculated by summing the total number of voxels representing the vessel lumen in the grid region within the ROI, Nv, and the total number of voxels in the ROI within each region, NROI. The CVI for each region was calculated as  
\begin{equation}CV{I_{reg}} = {\rm{\;}}\frac{{{N_v}}}{{{N_{RoI}}}}{\rm{\;}} \times 100.\end{equation}
Statistical Data Analysis
MATLAB (MATLAB 2020a; The MathWorks, Inc.) was used for statistical analysis. Normality of the data was tested with the Lilliefors test. For the analysis of topographical distribution, measurements from all nine scans were averaged per patient and are presented as a median and interquartile range. The quartile-based coefficient of variation39 (CV) was calculated (equation 5) where Q25 and Q75 denote the 25th and 75th quantiles of the distribution.  
\begin{equation}CV = \left({\frac{{{Q_{75}} - {\rm{\;}}{Q_{25}}}}{{{Q_{75}} + {\rm{\;}}{Q_{25}}}}} \right) \times 100.\end{equation}
   Correlations of the choroidal parameters with AEL were performed using Spearman rank correlation.40 The coefficient of repeatability (CR) from nonnormal data was derived for each ETDRS region as explained elsewhere.41,42 In short, the median of each of the three measurements per day and its difference to the overall mean across the 3 days were calculated. Subsequently, the resulting deltas were evaluated in a cumulative distribution function. The 95% limits of the cumulative distribution function were determined and multiplied by \(\sqrt {3/2} \) as a correction factor for centered data. The CR was considered half the length of the 95% interval. It is noteworthy that only intersession but not intrasession repeatability was considered in the analysis, to reflect most studies that investigate choroidal metrics in an intersession rather than intrasession (consecutive) manner. 
Widefield Topography of ChT and CVI
Presented ChT and CVI were sorted by retinal location and eccentricity, as presented in Table 1. In general, the ChT decreased with increasing eccentricity and was typically highest superiorly, followed by the temporal, inferior, and nasal areas. In contrast, CVI showed a relatively uniform pattern with almost no eccentricity- or region-dependent fluctuations. Quartile-based CV across the scanned retina were 5.92 for ChT and 0.89 for CVI. There were no significant correlations found between ChT and CVI in any ETDRS sector (all P > 0.05). However, there appears to be a general transition from negative to increasingly positive correlations between the parameters with increasing eccentricity. 
Table 1.
Topographical Distribution of ChT and CVI across the Scan Area and Their Correlation with Each Other
Table 1.
Topographical Distribution of ChT and CVI across the Scan Area and Their Correlation with Each Other
Figure 2 shows an example of the CVI algorithm output with typical topographic variations for ChT and CVI. As depicted in Figure 2b and Figure 2c, the structures of larger choroidal vessels are highly visible with a subsequently higher CVI. This relationship is not surprising, since these pixel columns will comprise predominantly vessel lumina, yielding a high CVI in these localized areas. 
Figure 2.
Output of the CVI algorithm for an example participant: (a) ChT segmentation as ROI. (b) Calculated number of vessel pixels. (c) CVI map without ETDRS grid. (d) CVI map with overlaying ETDRS grid, adjusted for image magnification. The white area of the ONH was excluded from analysis.
Figure 2.
Output of the CVI algorithm for an example participant: (a) ChT segmentation as ROI. (b) Calculated number of vessel pixels. (c) CVI map without ETDRS grid. (d) CVI map with overlaying ETDRS grid, adjusted for image magnification. The white area of the ONH was excluded from analysis.
Relationship between Axial Length and Choroidal Metrics
The relationship of both choroidal parameters with AEL was investigated to explore the potential effects of eye size on ChT and CVI. While ChT decreased with increasing AEL in all scan fields, CVI showed a more variable relationship with AEL (Table 2). The central and inner fields correlated positively with AEL, whereas negative correlations were found in the most peripheral regions. Despite the consistency in the pattern of ChT, the relationship with AEL was statistically significant only in the inner superior field (P = 0.02), and there were no significant correlations between CVI and AEL. 
Table 2.
Spearman Correlation Coefficient (ρ) and Associated P Value of Choroidal Metrics with Axial Eye Length
Table 2.
Spearman Correlation Coefficient (ρ) and Associated P Value of Choroidal Metrics with Axial Eye Length
Intersession Repeatability
The CR between the three measurement sessions for choroidal metrics is reported in Table 3. There were small regional fluctuations for ChT, ranging from 10.37 µm (inner temporal and extended inferior fields) to 20.33 µm (extended nasal field). CVI repeatability ranged between 3.90% in the outer inferior field and 5.51% in the temporal extended field. The quartile-based CVs for the repeatability across the total scan area were 11.11 and 9.66 for ChT and CVI, respectively. The last two columns of Table 3 represent the relative proportion between the repeatability coefficient and the absolute median values from Table 1. With exception of the nasal extended subfield, ChT exhibits a relatively better repeatability than CVI. 
Table 3.
CR for ChT and CVI across the ETDRS Sectors
Table 3.
CR for ChT and CVI across the ETDRS Sectors
The study undertaken investigated the widefield distribution of ChT and CVI, their intersession repeatability, and their relationship to axial eye length. Moreover, novel and advantageous methodology was applied as follows: 
  • (1) Swept-source OCT with an increased widefield scan area of 12 × 12 mm was used. A smaller region of interest with a diameter of 10 mm was evaluated due to potential information loss resulting from magnification correction using AEL. To our knowledge, this is the largest evaluated retinal area for CVI distribution, ChT and CVI repeatability, and axial length scaling correlation to date.
  • (2) The analysis of three-dimensional volume scans was possible thanks to enhanced scan depth and thus improved image quality compared to spectral-domain OCT. Previously, most CVI and ChT metrics were obtained from line scans or as single measurement points only.
  • (3) All algorithms used were optimized to run entirely automated (with the exception of the semiautomated exclusion of the optic nerve head) on all 1024 B-scans per volume. This includes image-processing steps for choroidal thickness and choroidal vascularity analysis, as well as the subsequent scaling of the ETDRS grid to account for image magnification. This approach reduces analysis time, examiner bias, and interobserver variation.
  • (4) The size of the ETDRS sectors was compensated for image magnification for each participant, from their AEL measurement. Moreover, the study measurements were controlled for choroidal diurnal rhythm by performing the measurements always at the same time of day for each participant. Eye movement artifacts during and between scans were minimized by using the inbuilt eye-tracking and follow-up mode of the device.
CVI Reveals Less Topographical Variations Than ChT
ChT showed a higher variability of topographical distribution across the retina compared to CVI. ChT generally decreased toward the periphery (particularly in the nasal sector), as reported previously.4346 CVI distribution patterns were more uniform with no obvious influence of retinal location and thus lower variability metrics. Regional consistency is also reflected in the low quartile-based CV of 0.89 for CVI, compared to 5.92 for ChT. This is consistent with earlier studies, although these found generally lower absolute CVI values of around 45% to 65%.12,18,47,48 These differences likely arise from the use of different OCT devices, scan patterns, image-processing algorithms, and participant age and ethnicity distributions. However, it has been stated that central CVI (measured in two dimensions) is not significantly influenced by the choice of scan pattern49 or OCT device.50 Moreover, ChT and CVI showed a tendency of eccentricity-dependent association with each other, with a shift from negative to positive correlations toward the peripheral retina. This suggests that the thicker central and inner choroid results from a higher stromal rather than vascular component. This is in contrast to previously reported positive correlations of ChT and CVI in the central retina up to 1.5 mm diameter,12 which is comparable to the central and inner fields defined in the current study. 
Eccentricity-Dependent Correlations of AEL with CVI
When establishing and adjusting normative databases, the interactions of choroidal metrics with AEL are an important factor to be considered. Like this, it is possible to differentiate physiologic from truly pathologic deviations. ChT is known to thin with increasing AEL.12,5153 The lack of statistical significance in the current study, except in the inner superior field, most likely arises from the relatively small sample size. CVI again exhibits mixed associations with AEL; centrally, longer ocular length leads to an increased vascular component and/or reduced stromal component. In contrast, higher AEL comes along with lower CVI in the peripheral retina. Interestingly, this eccentricity-dependent correlation pattern is the inverse of the ChT and CVI correlation described above. Previous correlation findings appear somewhat inconclusive; low negative correlations of CVI with AEL have been reported,12,18,47 while slightly increased vascular components with myopia have also been observed.20,21 However, these findings are all limited to the central or maximally the outer scan fields and do not all reach statistical significance. 
Repeatability of Choroidal Metrics Is Similar to Reported Effect Sizes
The analysis of intersession repeatability was included for two reasons: first, to evaluate whether there are more metrologically preferred retinal areas for measurements in future studies. These are areas that exhibit lower CR values. Second, the repeatability of these choroidal metrics is put into context of commonly found effect sizes. To address the first question, the nasal extended area manifests the highest CR with 20.33 µm compared to the other ETDRS sectors (10.37–16.35 µm). One earlier publication reported an average CR of 30 µm for ChT with automated segmentation in swept-source volume scans but without differentiating between ETDRS areas.54 Moreover, another methodologically similar study found the opposite for spectral-domain OCT scans, with the best repeatability in the nasal sectors due to thinner choroids and thus advantageous imaging properties.42 However, comparisons to other past literature are limited, since the range of methodologies used and analyses reported are broad. In contrast to ChT, repeatability of CVI measurements did not show any regional variations, surprisingly not even in the extended nasal subfield. Therefore, CVI can be measured irrespective of the retinal region from a metrologic perspective. There are few published reports on the repeatability of CVI measurements and all limited to spectral-domain and line scan analysis, making comparisons to the present study difficult. Moreover, these commonly evaluate the repeatability of luminal and stromal areas separately, which causes difficulties when translating it to repeatability results of final and summarized CVI measurements. However, the somewhat most comparable study found 95% limits of agreements ranging from approximately −4% to +3%.49 From a relative perspective, however, ChT exhibits a better repeatability compared to the absolute thickness than CVI (Table 3). This finding is in line with another study in healthy participants.12 There has not been a definite theory for this relation yet. However, one hypothesis suggests that the repeatability of CVI is dependent on the repeatability of ChT as a first and inherent processing step. Therefore, this finding can be expected. The second step was to compare reported effect sizes with measurement variability. Only if the CR does not exceed the effect size can measurement noise be neglected as a source of error. Typical reported differences between healthy and disease conditions are around 2% to 6% for CVI (review: Agrawal et al.11). For ChT, the effect size is highly dependent on the condition studied. In retinal disease, ChT alterations of around 100 µm are commonly found.35 However, in myopia research, the measurement reliability is similar to or even exceeds published effect sizes of maximally 20 µm (review: Read et al.55), even when using swept-source OCT technology with improved choroidal imaging. In summary, measurement repeatability should be carefully considered when interpreting “true” change or difference in both choroidal parameters and distinguishing from measurement variation. 
A relatively small sample size of n = 14 was used for this study. Although axial length was normally distributed and covered the range of usual emmetropic and myopic ocular sizes in the normal population, the small sample size could affect the homogeneity of the sample. It is therefore prudent to view the correlation results as indicative rather than definitive, and these relationships warrant further investigation with a larger sample with a wider range of axial lengths and stronger homogeneity. All included participants were young adults without ocular pathology and from Caucasian descent (except for one Indian participant). Therefore, the findings should be applied cautiously to other age groups and ethnicities, as well as in pathology. Moreover, participants could not be fully controlled for all factors shown to affect the choroid, such as smoking and caffeine. By enrollment into the study, participants were advised to avoid caffeine and nicotine intake at least for 1 hour before the measurements, but this could not be verified by the examiner. Blood pressure as another confounding factor was not measured specifically, as it was assumed to be stable within one individual because of the 10-minute resting and distance viewing period prior to the measurements. Despite imaging advances with swept-source and long-wavelength OCT, there were still B-scans in the evaluated subset with limited visibility of either the choroidal-scleral interface or vessels. The limited contrast of choroidal structures occurred due to OCT signal roll-off with increasing scan depth and therefore measurement noise. Retinal vessel shadows (areas of hyporeflectivity in the choroid from overlying structures) represent another measurement artifact and could lead to a misjudgment of CVI. However, these artifacts are limited to the extended nasal, inferior, and superior analysis regions. While this factor might have a small influence on CVI in those regions, it should not affect the repeatability results per se, as the retinal vessels are not transient. Therefore, these artifacts were acknowledged in these three out of 13 analysis regions but were not removed from the images due to their expected limited effect on the outcomes. The results of the extended nasal sector might have been additionally influenced by the semiautomated exclusion of the optic nerve head, leading to a reduced number of analyzed pixels in that region, which varied between participants based on the size of that feature. Despite the image size compensation for magnification, caused by axial length differences of the participants, they were not corrected for ocular curvature. Subsequently, this might lead to small biases in measurements of absolute ChT, especially in the periphery of the scan and with increasing curvature. However, the magnitude of error usually tends to be relatively small with up to 6.0 µm at a scan eccentricity of almost 7 mm.56 Moreover, no highly myopic eyes were included into the study and only affects the absolute measures of ChT but not the repeatability measurements, as these were analyzed in a relative manner to each other. Therefore, the total error in the current study is existent but most probably has only very limited effect on the final outcomes and conclusion. 
This study provides evidence that choroidal vascularity does not follow a clear and distinct topographical pattern across the central retina, unlike ChT. Moreover, mixed correlations between ChT and CVI, as well as CVI and axial length, were found, whereas choroidal thinning with increasing axial length was confirmed. It was further observed that the repeatability of CVI measurements does not vary across scan areas, in contrast to ChT measurements. However, the repeatability of both metrics is similar to or even exceeds previously reported effect sizes in many studies using choroidal parameters. Measurement repeatability should be considered in the interpretation of results from studies involving the choroid as a biomarker. 
The authors thank Focke Ziemssen and Laura Kühlewein for their support and cooperation in order to allow successful study measurements. 
Supported by Eberhard-Karls-University Tuebingen (ZUK 63) as part of the German Excellence initiative from the Federal Ministry of Education and Research (BMBF), Deutsche Forschungsgemeinschaft and Open Access Publishing Fund of University of Tübingen, and Carl Zeiss Vision International GmbH (SW). The funders did not have any additional role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript. 
Disclosure: K. Breher, None; L. Terry, None; T. Bower, None; S. Wahl, Carl Zeiss Vision International GmbH (E, F) 
Hogan MJ, Alvarado JA, Wedell JE. Histology of the Human Eye. Philadelphia: Saunders; 1971.
Nickla DL, Wallman J. The multifunctional choroid. Prog Retin Eye Res. 2010; 29(2): 144–168. [CrossRef] [PubMed]
Maruko I, Iida T, Sugano Y, Ojima A, Ogasawara M, Spaide RF. Subfoveal choroidal thickness after treatment of central serous chorioretinopathy. Ophthalmology. 2010; 117(9): 1792–1799. [CrossRef] [PubMed]
Maruko I, Iida T, Sugano Y, et al. Subfoveal choroidal thickness after treatment of Vogt–Koyanagi–Harada disease. Retina. 2011; 31(3): 510–517. [CrossRef] [PubMed]
Chung SE, Kang SW, Lee JH, Kim YT. Choroidal thickness in polypoidal choroidal vasculopathy and exudative age-related macular degeneration. Ophthalmology. 2011; 118(5): 840–845. [CrossRef] [PubMed]
Wallman J, Wildsoet C, Xu A, et al. Moving the retina: choroidal modulation of refractive state. Vision Res. 1995; 35(1): 37–50. [CrossRef] [PubMed]
Wang D, Chun RKM, Liu M, et al. Optical defocus rapidly changes choroidal thickness in schoolchildren. PLoS One. 2016; 11(8): e0161535. [CrossRef] [PubMed]
Fitzgerald MEC, Wildsoet CF, Reiner A. Temporal relationship of choroidal blood flow and thickness changes during recovery from form deprivation myopia in chicks. Exp Eye Res. 2002; 74(5): 561–570. [CrossRef] [PubMed]
Zhang S, Zhang G, Zhou X, et al. Changes in choroidal thickness and choroidal blood perfusion in guinea pig myopia. Invest Ophthalmol Vis Sci. 2019; 60(8): 3074–3083. [CrossRef] [PubMed]
Li Z, Cui D, Hu Y, Ao S, Zeng J, Yang X. Choroidal thickness and axial length changes in myopic children treated with orthokeratology. Cont Lens Anterior Eye. 2017; 40(6): 417–423. [CrossRef] [PubMed]
Agrawal R, Ding J, Sen P, et al. Exploring choroidal angioarchitecture in health and disease using choroidal vascularity index. Prog Retin Eye Res. 2020;77: 100829.
Agrawal R, Gupta P, Tan K-A, Cheung CMG, Wong T-Y, Cheng C-Y. Choroidal vascularity index as a measure of vascular status of the choroid: measurements in healthy eyes from a population-based study. Sci Rep. 2016; 6: 21090. [CrossRef] [PubMed]
Tan K-A, Laude A, Yip V, Loo E, Wong EP, Agrawal R. Choroidal vascularity index—a novel optical coherence tomography parameter for disease monitoring in diabetes mellitus? Acta Ophthalmol. 2016; 94(7): e612–e616. [CrossRef] [PubMed]
Agrawal R, Chhablani J, Tan K-A, Shah S, Sarvaiya C, Banker A. Choroidal vascularity index in central serous chorioretinopathy. Retina. 2016; 36(9): 1646–1651. [CrossRef] [PubMed]
Wei X, Ting DSW, Ng WY, Khandelwal N, Agrawal R, Cheung CMG. Choroidal vascularity index: a novel optical coherence tomography based parameter in patients with exudative age-related macular degeneration. Retina. 2017; 37(6): 1120–1125. [CrossRef] [PubMed]
Agrawal R, Salman M, Tan K-A, et al. Choroidal vascularity index (CVI)—a novel optical coherence tomography parameter for monitoring patients with panuveitis? PLoS One. 2016; 11(1): e0146344. [CrossRef] [PubMed]
Tan R, Agrawal R, Taduru S, Gupta A, Vupparaboina K, Chhablani J. Choroidal vascularity index in retinitis pigmentosa: an OCT study. Ophthalmic Surg Lasers Imaging Retina. 2018; 49(3): 191–197. [CrossRef] [PubMed]
Li Z, Long W, Hu Y, Zhao W, Zhang W, Yang X. Features of the choroidal structures in myopic children based on image binarization of optical coherence tomography. Invest Ophthalmol Vis Sci. 2020; 61(4): 18. [CrossRef] [PubMed]
Alshareef RA, Khuthaila MK, Januwada M, Goud A, Ferrara D, Chhablani J. Choroidal vascular analysis in myopic eyes: evidence of foveal medium vessel layer thinning. Int J Retin Vitr. 2017; 3(1): 28. [CrossRef]
Alshareef RA, Khuthaila MK, Goud A, Vupparaboina KK, Jana S, Chhablani J. Subfoveal choroidal vascularity in myopia: evidence from spectral-domain optical coherence tomography. Ophthalmic Surg Lasers Imaging Retina. 2017; 48(3): 202–207. [CrossRef] [PubMed]
Gupta P, Thakku SG, Saw S-M, et al. Characterization of choroidal morphologic and vascular features in young men with high myopia using spectral-domain optical coherence tomography. Am J Ophthalmol. 2017; 177: 27–33. [CrossRef] [PubMed]
Iovino C, Pellegrini M, Bernabei F, et al. Choroidal vascularity index: an in-depth analysis of this novel optical coherence tomography parameter. J Clin Med. 2020; 9(2): 595. [CrossRef]
Kong M, Eo DR, Han G, Park SY, Ham D-I. Error rate of automated choroidal segmentation using swept-source optical coherence tomography. Acta Ophthalmol. 2016; 94(6): e427–e431. [PubMed]
Kong M, Choi DY, Han G, et al. Measurable range of subfoveal choroidal thickness with conventional spectral domain optical coherence tomography. Transl Vis Sci Technol. 2018; 7(5): 16. [PubMed]
Lee SW, Yu S-Y, Seo KH, Kim ES, Kwak HW. Diurnal variation in choroidal thickness in relation to sex, axial length, and baseline choroidal thickness in healthy Korean subjects. Retina. 2014; 34(2): 385–393. [PubMed]
Tan CS, Ouyang Y, Ruiz H, Sadda SR. Diurnal variation of choroidal thickness in normal, healthy subjects measured by spectral domain optical coherence tomography. Invest Ophthalmol Vis Sci. 2012; 53(1): 261–266. [PubMed]
Dervisougullari MS, Totan Y, Yüce A, Kulak AE. Acute effects of caffeine on choroidal thickness and ocular pulse amplitude. Cutan Ocul Toxicol. 2016; 35(4): 281–286. [PubMed]
Zengin MO, Cinar E, Karahan E, Tuncer I, Kucukerdonmez C. The effect of caffeine on choroidal thickness in young healthy subjects. Cutan Ocul Toxicol. 2015; 34(2): 112–116. [PubMed]
Zengin MO, Cinar E, Kucukerdonmez C. The effect of nicotine on choroidal thickness. Br J Ophthalmol. 2014; 98(2): 233–237. [PubMed]
Woodman EC, Read SA, Collins MJ. Axial length and choroidal thickness changes accompanying prolonged accommodation in myopes and emmetropes. Vision Res. 2012; 72: 34–41. [PubMed]
Woodman-Pieterse EC, Read SA, Collins MJ, Alonso-Caneiro D. Regional changes in choroidal thickness associated with accommodation. Invest Ophthalmol Vis Sci. 2015; 56(11): 6414–6422. [PubMed]
Bagherinia H, Gregori G, Rosenfeld PJ, et al. A method for automated choroidal-scleral interface segmentation in optical coherence tomography. Invest Ophthalmol Vis Sci. 2019; 60(9): 143.
Early Treatment Diabetic Retinopathy Study Research Group. Photocoagulation for diabetic macular edema—Early Treatment Diabetic Retinopathy Study Report No. 1. Arch Ophthalmol. 1985; 103: 1796–1806. [PubMed]
Sonoda S, Sakamoto T, Yamashita T, et al. Luminal and stromal areas of choroid determined by binarization method of optical coherence tomographic images. Am J Ophthalmol. 2015; 159(6): 1123–1131. [PubMed]
Sonoda S, Sakamoto T, Yamashita T, et al. Choroidal structure in normal eyes and after photodynamic therapy determined by binarization of optical coherence tomographic images. Invest Ophthalmol Vis Sci. 2014; 55(6): 3893–3899. [PubMed]
Niblack W . An Introduction to Digital Image Processing. Englewood Cliffs: Prentice Hall; 1986: 115–116.
Terry L, Cassels N, Lu K, et al. Automated retinal layer segmentation using spectral domain optical coherence tomography: evaluation of inter-session repeatability and agreement between devices. PLoS One. 2016; 11(9): e0162001. [PubMed]
Bennett AG, Rudnicka AR, Edgar DF. Improvements on Littmann's method of determining the size of retinal features by fundus photography. Graefes Arch Clin Exp Ophthalmol. 1994; 232(6): 361–367. [PubMed]
Zwillinger D, Kokoska S. CRC Standard Probability and Statistics Tables and Formulae. Boca Raton, FL: CRC Press; 1999.
Myers L, Sirois MJ. Spearman correlation coefficients, differences between. In Kotz S, Read CB, Balakrishnan N, Vidakovic B, Johnson NL, eds. Encycl Stat Sci. 2006, doi:10.1002/0471667196.ess5050.pub2.
ISO/IEC Guide 100:2008. Evaluation of Measurement Data—Guide to the Expression of Uncertainty in Measurement. Geneva, Switzerland: Joint Committee for Guides in Metrology; 2008.
Breher K, Ohlendorf A, Wahl S. A metrological approach to the analysis of choroidal thickness by optical coherence tomography 3D scans in myopia research. Sci Rep. 2019; 9(1): 20322. [PubMed]
Mohler KJ, Draxinger W, Klein T, et al. Combined 60 wide-field choroidal thickness maps and high-definition en face vasculature visualization using swept-source megahertz OCT at 1050 nm. Invest Ophthalmol Vis Sci. 2015; 56(11): 6284–6293. [PubMed]
Hoseini-Yazdi H, Vincent SJ, Collins MJ, Read SA, Alonso-Caneiro D. Wide-field choroidal thickness in myopes and emmetropes. Sci Rep. 2019; 9(1): 1–10. [PubMed]
Breher K, García García M, Ohlendorf A, Wahl S. The effect of the optical design of multifocal contact lenses on choroidal thickness. PLoS One. 2018; 13(11): e0207637. [PubMed]
Rasheed MA, Singh SR, Invernizzi A, et al. Wide-field choroidal thickness profile in healthy eyes. Sci Rep. 2018; 8(1): 1–7. [PubMed]
Goud A, Singh SR, Sahoo NK, et al. New insights on choroidal vascularity: a comprehensive topographic approach. Invest Ophthalmol Vis Sci. 2019; 60(10): 3563–3569. [PubMed]
Singh SR, Invernizzi A, Rasheed MA, et al. Wide-field choroidal vascularity in healthy eyes. Am J Ophthalmol. 2018; 193: 100–105. [PubMed]
Agrawal R, Wei X, Goud A, Vupparaboina KK, Jana S, Chhablani J. Influence of scanning area on choroidal vascularity index measurement using optical coherence tomography. Acta Ophthalmol. 2017; 95(8): e770–e775. [PubMed]
Agrawal R, Seen S, Vaishnavi S, et al. Choroidal vascularity index using swept-source and spectral-domain optical coherence tomography: a comparative study. Ophthalmic Surg Lasers Imaging Retina. 2019; 50(2): e26–e32. [PubMed]
Hirata M, Tsujikawa A, Matsumoto A, et al. Macular choroidal thickness and volume in normal subjects measured by swept-source optical coherence tomography. Invest Ophthalmol Vis Sci. 2011; 52(8): 4971–4978. [PubMed]
Sanchez-Cano A, Orduna E, Segura F, et al. Choroidal thickness and volume in healthy young white adults and the relationships between them and axial length, ammetropy and sex. Am J Ophthalmol. 2014; 158(3): 574–83.e1. [PubMed]
Jin P, Zou H, Zhu J, et al. Choroidal and retinal thickness in children with different refractive status measured by swept-source optical coherence tomography. Am J Ophthalmol. 2016; 168: 164–176. [PubMed]
Mansouri K, Medeiros FA, Tatham AJ, Marchase N, Weinreb RN. Evaluation of retinal and choroidal thickness by swept-source optical coherence tomography: repeatability and assessment of artifacts. Am J Ophthalmol. 2014; 157(5): 1022–1032. [PubMed]
Read SA, Fuss JA, Vincent SJ, Collins MJ, Alonso-Caneiro D. Choroidal changes in human myopia: insights from optical coherence tomography imaging. Clin Exp Optom. 2019; 102(3): 270–285. [PubMed]
Kuo AN, McNabb RP, Chiu SJ, et al. Correction of ocular shape in retinal optical coherence tomography and effect on current clinical measures. Am J Ophthalmol. 2013; 156(2): 304–311. [PubMed]
Figure 1.
Extended ETDRS grid, showing the optic nerve head exclusion zone for a right eye. The ring diameters are 1 mm, 3 mm, 6 mm, and 10 mm.
Figure 1.
Extended ETDRS grid, showing the optic nerve head exclusion zone for a right eye. The ring diameters are 1 mm, 3 mm, 6 mm, and 10 mm.
Figure 2.
Output of the CVI algorithm for an example participant: (a) ChT segmentation as ROI. (b) Calculated number of vessel pixels. (c) CVI map without ETDRS grid. (d) CVI map with overlaying ETDRS grid, adjusted for image magnification. The white area of the ONH was excluded from analysis.
Figure 2.
Output of the CVI algorithm for an example participant: (a) ChT segmentation as ROI. (b) Calculated number of vessel pixels. (c) CVI map without ETDRS grid. (d) CVI map with overlaying ETDRS grid, adjusted for image magnification. The white area of the ONH was excluded from analysis.
Table 1.
Topographical Distribution of ChT and CVI across the Scan Area and Their Correlation with Each Other
Table 1.
Topographical Distribution of ChT and CVI across the Scan Area and Their Correlation with Each Other
Table 2.
Spearman Correlation Coefficient (ρ) and Associated P Value of Choroidal Metrics with Axial Eye Length
Table 2.
Spearman Correlation Coefficient (ρ) and Associated P Value of Choroidal Metrics with Axial Eye Length
Table 3.
CR for ChT and CVI across the ETDRS Sectors
Table 3.
CR for ChT and CVI across the ETDRS Sectors

This PDF is available to Subscribers Only

Sign in or purchase a subscription to access this content. ×

You must be signed into an individual account to use this feature.