Ocular fundus reference images from optical coherence tomography
Introduction
Several pathologic conditions of the retina have been shown to be correlated to abnormal vascular patterns [1]. The retina is regarded as a window through which to reflect on the vascular system, and changes seen in the retinal microvasculature have been found to be related to several cardiac diseases [2], [3]. Hence, the detection and quantitative characterization of the human retina vascular network are important research areas. Properties such as arteriolar diameter at bifurcation and tortuosity have been associated with death from ischemic heart disease [4], and arteriolar narrowing has been associated with increased mortality from stroke [4]. Furthermore, several characteristics of the retinal vascular network have been shown to be connected to the state of the cerebral vascular system [5], [6]. Also, studies have reported consistent associations between vascular abnormalities and several disorders such as diabetes [7] and hypertension [8], [9].
Extensive work has been done in this field, mainly based on two widely used imaging modalities: color fundus photography (CFP) and fluorescein angiography (FA). To facilitate ocular fundus imaging and improve CFP image quality, eyedrops are used to dilate the pupil, making this imaging modality minimally invasive. On the other hand, sodium fluorescein (NaFl, a fluorescent dye) is administered intravenously before the imaging to highlight the vascular network, making FA invasive.
Optical coherence tomography (OCT) is an imaging modality capable of noninvasively imaging tissue microstructure in vivo and in situ. It is an important tool in the diagnosis of ocular pathologic conditions and has been used extensively in research and in clinical practice. This modality uses the interferometer principle and is extensively described in the literature [10], [11].
The introduction of the commercial high-definition spectral-domain OCT (SD-OCT) models has made it possible to obtain large, high-resolution 3D scans of the retina (Fig. 1) with short acquisition time. An accurate fundus reference image computed from OCT data is still necessary, even though several OCT systems use a fundus image from a scanning laser ophthalmoscope (SLO) as a surrogate reference. Although highly important, the SLO image is not necessarily aligned with the OCT volume, as SLO images are taken after the OCT scan and not simultaneously. In addition, the SLO does not convey the same information as the OCT data because of the different working principle. By adding up the 3D OCT data depth-wise, a 2D reference image can be computed. From a vascular network point of view, this image correlates to the fundus images obtained from CFP or FA. As opposed to SLO, these OCT fundus images are created from OCT data and, as such, are intrinsically aligned to the volume.
An accurate reference image is highly important, as shown by the several approaches that have been proposed to achieve the co-registration of OCT scans [12], [13], [14]. These approaches have been performed to co-register OCT scans taken at different points in time and to develop complementary imaging modalities (expected to lead to a deeper insight into the pathophysiology of retinal lesions) [12], [13], [14].
From the 2D fundus reference image of the OCT 3D data, one can also segment important retinal features that complement the 3D information, such as the fovea, the optic disk, and the retinal vascular network, similarly to CFP [14], [15], [16], [17]. On standard SD-OCTs, the retinal blood vessels are not directly visible. Instead, hemoglobin absorbs the light at the wavelength used by these systems. Consequently, the signal is attenuated at the structures below the vessels [10], [18], [19] (Fig. 2). This effect is well known and has been used to obtain the 2D vascular network from the 3D OCT data [20], [21].
In this work, we propose a method to generate a fundus reference image from 3D OCT data with improved vascular network extension and contrast. It combines data from three distinct 2D images computed from the 3D OCT data – the mean value fundus (MVF), the expected value fundus (EVF), and the error to local median (ELM) – to create the principal component fundus (PCF) reference image. To assess the extent of vascular network provided by the proposed OCT fundus reference image, we compared it with those provided by raw CFP and FA. This comparison is based on the visible vascular network as detected by a human grader who manually segmented the vascular network in the three imaging modalities. FA is the gold standard for ocular vascular network imaging and is used in this work to establish the total vascular network within the imaged region.
Section snippets
Fundus reference images
A few approaches are used to compute a fundus reference image from the 3D OCT data, and they all share the idea of using depth-wise averaging (total or partial) of each individual A-scan [20], [22], [21], [23]. In this work, we compute a partial depth-wise average, the MVF, and propose two new ones (the EVF and the ELM). Finally, we combine all three to compute the PCF reference image. The parameters used throughout this section are defined in Section 2.3.
As a preprocessing step, all the
Results
The four fundus reference images computed from the OCT 3D data (MVF, EVF, ELM, and PCF) were qualitatively and quantitatively evaluated. Then, the vascular network from the best image was compared to that of CFP to assess its similarity (in extension).
Comments and conclusions
In this paper, we have presented a method for computing an ocular fundus reference image from an OCT volumetric scan of the human macula. While of major importance, as it allows for the exact localization of the scan within the human macula, several OCT systems perform a SLO scan, immediately after the OCT scan, to be used as a surrogate fundus reference. The non-simultaneous acquisition is prone to misalignments between the fundus reference and the OCT data, in the one hand and, the
Acknowledgements
The authors would like to thank Dr. Maria da Luz Cachulo and Dr. Isabel Pires for their support in the evaluation of the fundus images, and Telmo Miranda for performing the manual segmentations.
This work was supported by FCT under the projects PTDC/SAU-ENB/111139/2009 and PEST-C/SAU/UI3282/2013, and by the COMPETE programs FCOMP-01-0124-FEDER-015712 and FCOMP-01-0124-FEDER-037299.
References (29)
- et al.
Retinal microvascular abnormalities and 10-year cardiovascular mortality: a population-based case-control study
Ophthalmology
(2003) - et al.
Retinal microvascular abnormalities and their relationship with hypertension, cardiovascular disease, and mortality
Survey of Ophthalmology
(2001) - et al.
A coarse-to-fine strategy for automatically detecting exudates in color eye fundus images.
Computerized Medical Imaging and Graphics
(2010) - et al.
The use of radial symmetry to localize retinal landmarks
Computerized Medical Imaging and Graphics
(2013) - et al.
Locating the fovea center position in digital fundus images using thresholding and feature extraction techniques
Computerized Medical Imaging and Graphics
(2013) - et al.
Digital ocular fundus imaging: a review
Ophthalmologica
(2011) - et al.
Abnormalities of retinal microvascular structure and risk of mortality from ischemic heart disease and stroke
Hypertension
(2006) - et al.
Retinal vascular image analysis as a potential screening tool for cerebrovascular disease: a rationale based on homology between cerebral and retinal microvasculatures
Journal of Anatomy
(2005) - et al.
Retinal arterial changes correlate with cerebral small-vessel disease
Neurology
(2002) - et al.
Retinal vessel diameters and the incidence of gross proteinuria and renal insufficiency in people with type 1 diabetes
Diabetes
(2004)