Abstract
The choroid provides nutritional support for the retinal pigment epithelium and photoreceptors. Choroidal dysfunction plays a major role in several of the most important causes of vision loss including age-related macular degeneration, myopic degeneration, and pachychoroid diseases such as central serous chorioretinopathy and polypoidal choroidal vasculopathy. We describe an imaging technique using depth-resolved swept-source optical coherence tomography (SS-OCT) that provides full-thickness three-dimensional (3D) visualization of choroidal anatomy including topographical features of individual vessels. Enrolled subjects with different clinical manifestations within the pachychoroid disease spectrum underwent 15 mm × 9 mm volume scans centered on the fovea. A fully automated method segmented the choroidal vessels using their hyporeflective lumens. Binarized choroidal vessels were rendered in a 3D viewer as a vascular network within a choroidal slab. The network of choroidal vessels was color depth-encoded with a reference to the Bruch’s membrane segmentation. Topographical features of the choroidal vasculature were characterized and compared with choroidal imaging obtained with indocyanine green angiography (ICGA) from the same subject. The en face SS-OCT projections of the larger choroid vessels closely resembled to that obtained with ICGA, with the automated SS-OCT approach proving additional depth-encoded 3D information. In 16 eyes with pachychoroid disease, the SS-OCT approach added clinically relevant structural details, including choroidal thickness and vessel depth, which the ICGA studies could not provide. Our technique appears to advance the in vivo visualization of the full-thickness choroid, successfully reveals the topographical features of choroidal vasculature, and shows potential for further quantitative analysis when compared with other choroidal imaging techniques. This improved visualization of choroidal vasculature and its 3D structure should provide an insight into choroid-related disease mechanisms as well as their responses to treatment.
Keywords
Impact statement
The choroid plays an essential role in the pathogenesis of many retinal disorders. Traditional imaging techniques used to visualize the choroid, including indocyanine green angiography, are often invasive and typically lack depth-resolved information. Herein, swept-source optical coherence tomography (SS-OCT) imaging methods are applied to visualize full-thickness choroidal structure and investigate its topographical features and choroidal thickness. Our technique encodes full-thickness choroidal slabs obtained with SS-OCT into three-dimensional (3D) networks of choroidal vessels, enabling in vivo visualization of individual choroidal vessels that provides potential for more precise quantitative analysis when compared to other choroidal imaging techniques. This improved visualization of choroidal vasculature structure may prove useful in providing greater insight into the choroid’s role in several of the most important causes of vision loss including age-related and myopic degeneration, diabetic retinopathy, and pachychoroid diseases.
Introduction
Improved in vivo imaging of choroidal structure has been of great clinical interests for decades due to the choroid’s critical role in common ocular diseases causing vision loss.1,2 Located below the retina and anterior to the sclera, the choroid is a tissue layer composed of blood vessels and stroma that is mainly responsible for delivering nutrients and supporting the metabolic exchange between the systemic circulation and the retinal pigment epithelium (RPE) and photoreceptors. The choroid contributes blood flow to the retina and optic nerve and maintains ocular thermoregulation.3–5 Choroidal blood flow accounts for 85% of all ocular blood flow, and choroidal dysfunction is strongly implicated in many ocular diseases including age-related macular degeneration, myopic degeneration, inherited retinal degenerations, and central serous chorioretinopathy (CSC).5,6 Because the choroid is located beneath the RPE and Bruch’s membrane (BrM), several anatomical and physiological compartments need to be crossed before choroid can be visualized with optical imaging devices. 7 Traditional clinical imaging technologies using short optical wavelengths including fundus autofluorescence and fluorescein angiography fail to image the full-thickness choroid due to light absorption by pigment in the photoreceptors, RPE, and within the choroid.8,9 Multispectral imaging has been explored to selectively image retinal and choroidal structures in traditional fundus photography. Studies have shown that the retina is sensitive to visible light, while the choroid can be imaged with near infrared light.10–12 Indocyanine green angiography (ICGA) utilizes deeper penetrating longer wavelengths capable of imaging the choroid’s full thickness; however, ICGA provides 2D images lacking depth information. 13 Moreover, ICGA is an invasive test with rare, but potentially, serious adverse reactions to the intravenous dye. 14 Optical coherence tomography (OCT) is a non-invasive imaging technique providing high-resolution depth-resolved images of ocular structures. Recent advances in OCT technology permit the rapid acquisition of dense volume scans which, when rendered in three dimensional (3D), enable functional ocular imaging. 15 Spectral domain OCT (SD-OCT) utilizes a technique called enhanced-depth imaging to mitigate the sensitivity roll off limitation, which enhances the measurement sensitivity for the structures below the RPE/BrM complex. 16 Swept-source optical coherence tomography (SS-OCT) uses a frequency swept laser that enables interference measurements at different optical frequencies or sequentially over time. This significantly improves penetration of the sampling beam and signal fall-off compared to SD-OCT.17,18 By using a longer wavelength of light, SS-OCT has been shown to provide superior choroidal imaging than SD-OCT. 15
Despite recent progress using SS-OCT for choroidal imaging, methods for segmentation and visualization of individual choroidal vessels from SS-OCT scans remain limited.6,15 To date, studies using SS-OCT to study the choroid have measured choroid thickness on individual B-scans or over the entire scanning region.19–24 However, since choroid thickness may be influenced by factors such as interstitial fluid and inflammatory infiltration, choroidal thickness measurements may not accurately reflect choroidal blood volume.
Choroidal vessels appear as dark lumens in en face OCT choroidal slabs. Choroidal vascularity index (CVI) is typically calculated as the area ratio in a small subfoveal region at a limited number of locations or along a relatively small number of B-scans.25,26 Although the CVI is potentially useful to quantify choroidal vascular density, it cannot reveal topographical features of individual choroidal vessels.27,28 While in some diseases, choroidal thickness and morphology of individual choroidal vessels may be atypical throughout the entire choroid, in others there may be focal or regional differences in these features.29,30 3D visualization of the choroid would be more likely to identify and quantify this local variability.
Challenges in complicated preprocessing efforts, dependence on modeling, heavy computational labor, and availability of a user-friendly 3D viewer have prevented researchers and clinicians from accessing the benefits of 3D choroidal visualization.27,31–36 In our preliminary study, we pioneered an attenuation correction-assisted method to automatically quantify the choroid using the entire SS-OCT 3D dataset. 37 This strategy was shown to increase the inter-layer contrast, enabling more accurate segmentation of choroid and choroid vessels, and minimizing the shadowing effects from large retinal vessels.37,38 To achieve an easy-to-access in vivo visualization of choroidal vasculature, herein we introduce an innovative approach to segment choroidal vessels and automatically present the binarized choroidal vessel network in a 3D viewer. In addition to visualization of vessel anatomy, we show that additional information may be mapped on the vessel’s surface. Since choroid thickness is a key metric in evaluating choroid-related eye diseases, we added color-coded depth information indicating the distance from the choroidal vessel surface to BrM. The choroidal vessels are presented in a 3D reconstruction of the choroidal vasculature in which vascular depth can be directly visualized.
Materials and methods
Subjects and choroid imaging
Patients were enrolled in a prospective SS-OCT imaging study approved by Western Institutional Review Board (Olympia, WA). Informed consent was obtained from each subject before imaging. The study was performed in accordance with the tenets of the Declaration of Helsinki and complied with the Health Insurance Portability and Accountability Act of 1996.
SS-OCT and ICGA images of the eye with most evident clinical manifestations in subjects with pachychoroid disease were selected for the present analysis. Patients were diagnosed as having pachychoroid disease if presenting with characteristics referable to conditions within the pachychoroid disease spectrum in at least one eye, as evaluated using a complex of clinical and multimodal imaging parameters. 30 Widefield ICGA images were acquired (Optos California, Dunfermline, Scotland, UK) on the same day of the SS-OCT study. Early-to-mid-phase ICGA of the detailed choroidal vascular were screened, and the best angiograms were used for image analysis. OCT data were acquired by a commercial SS-OCT system (PlexElite 9000, Carl Zeiss Meditec Inc.). The central wavelength of the light source is 1050 nm with a bandwidth of 100 nm, yielding a full width of half maximum axial resolution of ∼5 μm in tissue and a lateral resolution of ∼20 μm at the surface of retina. The A-scan scanning rate is 100 kHz and the A-scan depth is 3.0 mm in tissue sampled with 1536 pixels. The scanning protocol is 15 mm × 9 mm, which consisted of 834 B-scans (slow 15 mm scan), and each B-scan consisted of 500 A-scans (fast 9 mm scan). This SS-OCT prototype is equipped with FastTrac™ motion correction software to reduce motion artifacts. 39 Scans with signal intensity less than 7 were excluded according to the manufacturer’s instructions.
Choroid vessel segmentation
Choroid vessels were segmented from the 3D SS-OCT dataset using a previously published method.
37
Attenuation correction was first applied to eliminate the projection artifacts from retinal vessels and to enhance the contrast of choroid layer (Figure 1(a) and (b)). On each B-scan, the choroidal slab was segmented from BrM to the outer boundary of choroidal vessels using an automated graph search-based algorithm. The candidate nodes were firstly searched at the inflexion points in the vertical direction where OCT signals change from dark to bright at the vessel bottoms and then selected for the shortest path across the bottom of choroidal vessels.
37
Choroidal vessels were further segmented from the choroidal slab using Otsu’s binarization method (Figure 1(c)). Values of pixels in the OCT dataset were assigned as

Depth-encoded 3D visualization of choroid vessels. Attenuation correction was applied to the original SS-OCT scan (a and b). Choroidal vessels were segmented, binarized (c) and color-coded using the distance to BrM (d). The reassembled 3D choroidal vessel network without (e to g) and with depth-encoding (h to j) was visualized using the ImageJ 3D viewer with top view, side view and bottom view from left to right, respectively. (A color version of this figure is available in the online journal.)
En face choroidal vasculature
En face image of choroidal vasculature was generated by sum projection of the binarized choroidal vessels in 3D dataset. Histogram was redistributed to adjust the contrast for better presentation of small vessels by
Depth-encoded choroidal vessel
The distance of choroidal vessel to the BrM was color-coded via
3D visualization of choroid vessels
The choroidal vasculature can be visualized at any angle in the ImageJ 3D viewer. Figure 1(e) to (j) shows representative images of views from the anterior vitreous (top view, Figure 1(e) and (h)), from the peripheral retina (side view, Figure 1(f) and (i)) or from the sclera (bottom view, Figure 1(g) and (j)).
Results
We initially analyzed the agreement of choroidal imaging between SS-OCT and ICGA. Sixteen eyes underwent SS-OCT and ICGA imaging, and each was imaged on the same day. Statistical analysis on the repeatability of segmentation and its correlation with aging in normal eyes has been studied in our previous work. 38 En face images of the choroid vasculature were generated via sum projection of the vessels and compared with the same region cropped from a corresponding ICGA image (Figure 2). Due to the complexity of comparing resolution and correcting distortion in the two types of images, quantitative comparison was not attempted. Instead, skeletons of the major vessels visualized in the en face SS-OCT image were manually drawn by experts (first by HZ, and then confirmed by RKW and KBF) overlaid on the ICGA image. We observed similar patterns for the major branches of choroidal vessels in both SS-OCT and ICGA imaging in all 16 eyes.

Comparison of choroidal vasculature obtained from SS-OCT and ICGA. (a and b) 3D view of the choroid from a 75-year-old female patient with central serous chorioretinopathy. (c) En face images of choroidal vasculature generated via sum projection. (d) Cropped ICGA image at the same region. (e) Skeleton of major choroidal vessels (red) manually drawn on the en face image (c). (f) ICGA image with the vessel skeleton overlaid. (A color version of this figure is available in the online journal.)
The depth-encoded 3D visualization of choroidal vasculature provided: (1) the proportion of the choroid volume occupied by vessels; (2) a curved BrM surface; (3) curved choroid–scleral interface; (4) an appreciation of connectivity between individual vessels within the vascular network; (5) visualization of choroidal thickness throughout the volume. Figure 3 shows examples of 3D views of an eye with a thin choroid and an eye with a thick choroid. Vessels that were more than 150 µm away from BrM were presented in yellow to red color at the bottom. Choroidal vessels of medium size found mostly in Sattler’s layer were seen best in the top view, while larger choroidal vessels typically present in Haller’s layer were best seen from the bottom view. Side views provided the curvature of both BrM and choroid–scleral junction (Figure 3(b) and (e)). Branching vessels from the large vessels appeared more frequently in the top view than in the bottom view.

Examples of depth-encoded 3D visualization of a thin (Top row) and thick (Bottom row) choroid. Top views from the vitreous (a and d), side views from the peripheral (b and e), and bottom views from the sclera (c and f) were obtained from the 3D view of an eye with thin choroid (Top row), and an eye with thick choroid (Bottom row), respectively. (A color version of this figure is available in the online journal.)
Figure 4 shows a depth-encoded 3D visualization of choroid vessels in the left eye of a 71–year-old man with CSC. OCT B-scans show retinal thinning and outer retinal disruption subsequent to resorption of long-lasting subretinal fluid in the context of CSC (Figure 4). Most of the choroidal vessels in this eye are confined to a superficial thin layer, which is shown as a green network. A large vein located in the outer choroid of the macular region appears in red (bottom view, Figure 4(c)). A similar vascular morphology is observed in ICGA images; however, the depth of the large vein passing through macula and draining into the inferotemporal vortex vein cannot be determined (Figure 4(d)). From the corresponding SS-B-scans, we also observed the large dark lumen beneath the macula near the choroid–scleral junction that corresponds to the large vein (Figure 4(e)).

An example of depth-encoded 3D visualization of choroid vessels in the left eye of a 71-year-old male with resolved central serous chorioretinopathy. Top view from the vitreous (a), side view from the peripheral retina (b), and bottom view from the sclera (c) show the choroidal vessels segmented from SS-OCT dataset. (d) The corresponding ICGA image shows a large choroidal vessel passing through macula (arrow). (e) En face choroidal vasculature map with orthogonal SS-OCT images shown with and without attenuation correction. (A color version of this figure is available in the online journal.)
Figure 5 shows another depth-encoded 3D visualization of choroid vessels in an eye with pachychoroid disease (Figure 5). This patient is a 19-year-old man with fibrotic changes after chronic CSC with resolved subretinal fluid. An anomalous drainage pattern is observed at the macula in both Sattler’s layer from top view and Haller’s layer from the bottom view (Figure 5(a) and (c)). These vessels and their connectivity were poorly visualized on the ICGA (Figure 5(d)) and in the SS-OCT B-scans (Figure 5(e)).

Depth-encoded 3D visualization of the right eye of a 19-year-old male patient with resolved central serous chorioretinopathy. Top view from the vitreous (a), side view from the peripheral retina (b), and bottom view from the sclera (c) of the choroidal vessels segmented from an SS-OCT dataset. (d) The corresponding ICGA image shows shadows at macula (arrow). (e) En face choroidal vasculature map with orthogonal SS-OCT images with and without attenuation correction. (A color version of this figure is available in the online journal.)
Discussion
The 3D visualizations of the choroidal vasculature described in this work were reconstructed from binarized choroidal vessels segmented from SS-OCT volume scans. This segmentation method has been shown to have good repeatability 37 and has been used to quantify choroid thickness and CVI in normal and diseased eyes from other 3D SS-OCT datasets.38,41 By representing the choroidal vessel network in 3D, we were able to observe morphological features not appreciated with regional thickness maps or standard measures of CVI. Furthermore, the depth-encoded visualization provided information of choroid thickness and the depth of each vessel below BrM. This information will provide opportunity to better assess the choroid’s role in disease mechanisms such as in pachychoroid disease in which compressive effects overlying anteriorly displaced congested veins are believed to induce localized ischemia of the RPE and outer retina. 30
By comparing the reconstructed choroidal vessels from SS-OCT with that obtained from ICGA, we found that SS-OCT was able to provide comparable en face representations of the choroidal vasculature. Furthermore, the depth-encoded information provided additional features that can distinguish certain vessels of interest. We also found that SS-OCT was able to present choroidal vessels in greater detail compared to ICGA, especially in eyes presenting with profound anatomic alterations affecting dye distribution and visualization of sub-RPE structures, and to identify choroidal vascular structures that were not identified by ICGA, which could represent non-perfused vessels.
Although 3D visualization of choroidal vessels has been explored using SD-OCT and SS-OCT previously, most studies analyzed views from vitreous or sclera.27,31,32,34–36,42,43 By encoding the distance to BrM for each vessel, we were able to distinguish the large choroidal vessels in deeper choroid layers especially in eyes with thick choroids. This feature should provide better assessment of choroid in eyes where choroid is significantly thickened by diseases such as CSC. 44 In addition, we found that the side view of the entire volume revealed features such as regional curvature of the surfaces of BrM and sclera which should be useful to assess pathologies such as choroidal tumors and myopic degeneration.
There are some limitations of this work. Firstly, the en face choroidal vasculature was compared with that from ICGA qualitatively. Quantitative analysis is challenging because the vessel segmentations from either OCT images or ICGA images would depend on image resolution and binarization thresholding. We manually overlaid the skeleton of major vessels and observed consistency of both imaging methodologies. Secondly, we reported features of the depth-encoded 3D visualization of choroidal vessels based on limited number of cases. Advantages and disadvantages of this 3D visualization method remain to be explored in more eyes with various pathologies.
Conclusions
We successfully segmented the choroidal vessels from the 3D OCT dataset and color-coded the distance to BrM on the surface of choroidal vessels. The choroidal vasculatures extracted from SS-OCT were compared with those from ICGA and showed comparable patterns of vessel network. Moreover, the 3D network generated from SS-OCT was able to simultaneously provide choroid thickness and vessel depth below BrM information which should prove useful for the study of choroid-related disease.
Footnotes
AUTHORS’ CONTRIBUTIONS
All authors participated in the design, interpretation of the studies and analysis of the data and review of the manuscript. TB conducted the patient enrollment and image acquiring; HZ conducted the image processing and data analysis; HZ and RKW wrote the manuscript.
DECLARATION OF CONFLICTING INTERESTS
The author(s) declared the following potential conflicts of interest with respect to the research, authorship, and/or publication of this article: KBF is a consultant for Optovue, Zeiss, Allergan, Genentech, Heidelberg Engineering, Novartis, and Roche. RKW discloses intellectual property owned by the Oregon Health and Science University and the University of Washington. RKW also receives research support from Tasso Inc., Moptim Inc., and Colgate Palmolive Company. He is a consultant to Insight Photonic Solutions and Carl Zeiss Meditec. HZ and TB have no disclosures.
FUNDING
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This research is supported by grants from the National Eye Institute (R01EY028753), Carl Zeiss Meditec. The funding organization had no role in the design or conduct of this research.
