Measurement of retinal vessel widths from fundus images based on 2-D modeling

Changes in retinal vessel diameter are an important sign of diseases such as hypertension, arteriosclerosis and diabetes mellitus. Obtaining precise measurements of vascular widths is a critical and demanding process in automated retinal image analysis as the typical vessel is only a few pixels wide. This paper presents an algorithm to measure the vessel diameter to subpixel accuracy. The diameter measurement is based on a two-dimensional difference of Gaussian model, which is optimized to fit a two-dimensional intensity vessel segment. The performance of the method is evaluated against Brinchmann-Hansen's half height, Gregson's rectangular profile and Zhou's Gaussian model. Results from 100 sample profiles show that the presented algorithm is over 30% more precise than the compared techniques and is accurate to a third of a pixel

Assessment of total retinal blood flow using Doppler Fourier Domain Optical Coherence Tomography during systemic hypercapnia and hypocapnia.

The purpose of this study was to investigate changes in total retinal blood flow (RBF) using Doppler Fourier Domain Optical Coherence Tomography (Doppler FD-OCT) in response to the manipulation of systemic partial pressure of CO2 (PETCO2). Double circular Doppler blood flow scans were captured in nine healthy individuals (mean age ± standard deviation: 27.1 ± 4.1, six males) using the RTVue(™) FD-OCT (Optovue). PETCO2 was manipulated using a custom-designed computer-controlled gas blender (RespirAct(™)) connected to a sequential gas delivery rebreathing circuit. Doppler FD-OCT measurements were captured at baseline, during stages of hypercapnia (+5/+10/+15 mmHg PETCO2), return to baseline and during stages of hypocapnia (-5/-10/-15 mmHg PETCO2). Repeated measures analysis of variance (reANOVA) and Tukeys post hoc analysis were used to compare Doppler FD-OCT measurements between the various PETCO2 levels relative to baseline. The effect of PETCO2 on TRBF was also investigated using linear regression models. The average RBF significantly increased by 15% (P < 0.0001) with an increase in PETCO2 and decreased significantly by 10% with a decrease in PETCO2 (P = 0.001). Venous velocity significantly increased by 3.11% from baseline to extreme hypercapnia (P < 0.001) and reduced significantly by 2.01% at extreme hypocapnia (P = 0.012). No significant changes were found in the average venous area measurements under hypercapnia (P = 0.36) or hypocapnia (P = 0.40). Overall, increased and decreased PETCO2 values had a significant effect on RBF outcomes (P < 0.002). In healthy individuals, altered end-tidal CO2 levels significantly changed RBF as measured by Doppler FD-OCT

Joint segmentation and classification of retinal arteries/veins from fundus images

Objective Automatic artery/vein (A/V) segmentation from fundus images is required to track blood vessel changes occurring with many pathologies including retinopathy and cardiovascular pathologies. One of the clinical measures that quantifies vessel changes is the arterio-venous ratio (AVR) which represents the ratio between artery and vein diameters. This measure significantly depends on the accuracy of vessel segmentation and classification into arteries and veins. This paper proposes a fast, novel method for semantic A/V segmentation combining deep learning and graph propagation. Methods A convolutional neural network (CNN) is proposed to jointly segment and classify vessels into arteries and veins. The initial CNN labeling is propagated through a graph representation of the retinal vasculature, whose nodes are defined as the vessel branches and edges are weighted by the cost of linking pairs of branches. To efficiently propagate the labels, the graph is simplified into its minimum spanning tree. Results The method achieves an accuracy of 94.8% for vessels segmentation. The A/V classification achieves a specificity of 92.9% with a sensitivity of 93.7% on the CT-DRIVE database compared to the state-of-the-art-specificity and sensitivity, both of 91.7%. Conclusion The results show that our method outperforms the leading previous works on a public dataset for A/V classification and is by far the fastest. Significance The proposed global AVR calculated on the whole fundus image using our automatic A/V segmentation method can better track vessel changes associated to diabetic retinopathy than the standard local AVR calculated only around the optic disc.Comment: Preprint accepted in Artificial Intelligence in Medicin

Tram-Line filtering for retinal vessel segmentation

The segmentation of the vascular network from retinal fundal images is a fundamental step in the analysis of the retina, and may be used for a number of purposes, including diagnosis of diabetic retinopathy. However, due to the variability of retinal images segmentation is difficult, particularly with images of diseased retina which include significant distractors. This paper introduces a non-linear filter for vascular segmentation, which is particularly robust against such distractors. We demonstrate results on the publicly-available STARE dataset, superior to Stare’s performance, with 57.2% of the vascular network (by length) successfully located, with 97.2% positive predictive value measured by vessel length, compared with 57% and 92.2% for Stare. The filter is also simple and computationally efficient

Three-Dimensional Spectral-Domain Optical Coherence Tomography Data Analysis for Glaucoma Detection

Purpose: To develop a new three-dimensional (3D) spectral-domain optical coherence tomography (SD-OCT) data analysis method using a machine learning technique based on variable-size super pixel segmentation that efficiently utilizes full 3D dataset to improve the discrimination between early glaucomatous and healthy eyes. Methods: 192 eyes of 96 subjects (44 healthy, 59 glaucoma suspect and 89 glaucomatous eyes) were scanned with SD-OCT. Each SD-OCT cube dataset was first converted into 2D feature map based on retinal nerve fiber layer (RNFL) segmentation and then divided into various number of super pixels. Unlike the conventional super pixel having a fixed number of points, this newly developed variable-size super pixel is defined as a cluster of homogeneous adjacent pixels with variable size, shape and number. Features of super pixel map were extracted and used as inputs to machine classifier (LogitBoost adaptive boosting) to automatically identify diseased eyes. For discriminating performance assessment, area under the curve (AUC) of the receiver operating characteristics of the machine classifier outputs were compared with the conventional circumpapillary RNFL (cpRNFL) thickness measurements. Results: The super pixel analysis showed statistically significantly higher AUC than the cpRNFL (0.855 vs. 0.707, respectively, p = 0.031, Jackknife test) when glaucoma suspects were discriminated from healthy, while no significant difference was found when confirmed glaucoma eyes were discriminated from healthy eyes. Conclusions: A novel 3D OCT analysis technique performed at least as well as the cpRNFL in glaucoma discrimination and even better at glaucoma suspect discrimination. This new method has the potential to improve early detection of glaucomatous damage. © 2013 Xu et al

Optical Coherence Tomography Angiography Vessel Density in Healthy, Glaucoma Suspect, and Glaucoma Eyes.

PurposeThe purpose of this study was to compare retinal nerve fiber layer (RNFL) thickness and optical coherence tomography angiography (OCT-A) retinal vasculature measurements in healthy, glaucoma suspect, and glaucoma patients.MethodsTwo hundred sixty-one eyes of 164 healthy, glaucoma suspect, and open-angle glaucoma (OAG) participants from the Diagnostic Innovations in Glaucoma Study with good quality OCT-A images were included. Retinal vasculature information was summarized as a vessel density map and as vessel density (%), which is the proportion of flowing vessel area over the total area evaluated. Two vessel density measurements extracted from the RNFL were analyzed: (1) circumpapillary vessel density (cpVD) measured in a 750-μm-wide elliptical annulus around the disc and (2) whole image vessel density (wiVD) measured over the entire image. Areas under the receiver operating characteristic curves (AUROC) were used to evaluate diagnostic accuracy.ResultsAge-adjusted mean vessel density was significantly lower in OAG eyes compared with glaucoma suspects and healthy eyes. (cpVD: 55.1 ± 7%, 60.3 ± 5%, and 64.2 ± 3%, respectively; P < 0.001; and wiVD: 46.2 ± 6%, 51.3 ± 5%, and 56.6 ± 3%, respectively; P < 0.001). For differentiating between glaucoma and healthy eyes, the age-adjusted AUROC was highest for wiVD (0.94), followed by RNFL thickness (0.92) and cpVD (0.83). The AUROCs for differentiating between healthy and glaucoma suspect eyes were highest for wiVD (0.70), followed by cpVD (0.65) and RNFL thickness (0.65).ConclusionsOptical coherence tomography angiography vessel density had similar diagnostic accuracy to RNFL thickness measurements for differentiating between healthy and glaucoma eyes. These results suggest that OCT-A measurements reflect damage to tissues relevant to the pathophysiology of OAG

Retinal Vessel Segmentation Using the 2-D Morlet Wavelet and Supervised Classification

We present a method for automated segmentation of the vasculature in retinal images. The method produces segmentations by classifying each image pixel as vessel or non-vessel, based on the pixel's feature vector. Feature vectors are composed of the pixel's intensity and continuous two-dimensional Morlet wavelet transform responses taken at multiple scales. The Morlet wavelet is capable of tuning to specific frequencies, thus allowing noise filtering and vessel enhancement in a single step. We use a Bayesian classifier with class-conditional probability density functions (likelihoods) described as Gaussian mixtures, yielding a fast classification, while being able to model complex decision surfaces and compare its performance with the linear minimum squared error classifier. The probability distributions are estimated based on a training set of labeled pixels obtained from manual segmentations. The method's performance is evaluated on publicly available DRIVE and STARE databases of manually labeled non-mydriatic images. On the DRIVE database, it achieves an area under the receiver operating characteristic (ROC) curve of 0.9598, being slightly superior than that presented by the method of Staal et al.Comment: 9 pages, 7 figures and 1 table. Accepted for publication in IEEE Trans Med Imag; added copyright notic

Portable dynamic fundus instrument

A portable diagnostic image analysis instrument is disclosed for retinal funduscopy in which an eye fundus image is optically processed by a lens system to a charge coupled device (CCD) which produces recordable and viewable output data and is simultaneously viewable on an electronic view finder. The fundus image is processed to develop a representation of the vessel or vessels from the output data