A fuzzy vessel tracking algorithm for retinal images based on fuzzy clustering

Accurate image analysis of the retina using hessian matrix and binarisation of thresholded entropy with application of texture mapping

In this paper, we demonstrate a comprehensive method for segmenting the retinal vasculature in camera images of the fundus. This is of interest in the area of diagnostics for eye diseases that affect the blood vessels in the eye. In a departure from other state-of-the-art methods, vessels are first pre-grouped together with graph partitioning, using a spectral clustering technique based on morphological features. Local curvature is estimated over the whole image using eigenvalues of Hessian matrix in order to enhance the vessels, which appear as ridges in images of the retina. The result is combined with a binarized image, obtained using a threshold that maximizes entropy, to extract the retinal vessels from the background. Speckle type noise is reduced by applying a connectivity constraint on the extracted curvature based enhanced image. This constraint is varied over the image according to each region's predominant blood vessel size. The resultant image exhibits the central light reflex of retinal arteries and veins, which prevents the segmentation of whole vessels. To address this, the earlier entropy-based binarization technique is repeated on the original image, but crucially, with a different threshold to incorporate the central reflex vessels. The final segmentation is achieved by combining the segmented vessels with and without central light reflex. We carry out our approach on DRIVE and REVIEW, two publicly available collections of retinal images for research purposes. The obtained results are compared with state-of-the-art methods in the literature using metrics such as sensitivity (true positive rate), selectivity (false positive rate) and accuracy rates for the DRIVE images and measured vessel widths for the REVIEW images. Our approach out-performs the methods in the literature.

Tracing retinal vessel trees by transductive inference

BACKGROUND

Structural study of retinal blood vessels provides an early indication of diseases such as diabetic retinopathy, glaucoma, and hypertensive retinopathy. These studies require accurate tracing of retinal vessel tree structure from fundus images in an automated manner. However, the existing work encounters great difficulties when dealing with the crossover issue commonly-seen in vessel networks.

RESULTS

In this paper, we consider a novel graph-based approach to address this tracing with crossover problem: After initial steps of segmentation and skeleton extraction, its graph representation can be established, where each segment in the skeleton map becomes a node, and a direct contact between two adjacent segments is translated to an undirected edge of the two corresponding nodes. The segments in the skeleton map touching the optical disk area are considered as root nodes. This determines the number of trees to-be-found in the vessel network, which is always equal to the number of root nodes. Based on this undirected graph representation, the tracing problem is further connected to the well-studied transductive inference in machine learning, where the goal becomes that of properly propagating the tree labels from those known root nodes to the rest of the graph, such that the graph is partitioned into disjoint sub-graphs, or equivalently, each of the trees is traced and separated from the rest of the vessel network. This connection enables us to address the tracing problem by exploiting established development in transductive inference. Empirical experiments on public available fundus image datasets demonstrate the applicability of our approach.

CONCLUSIONS

We provide a novel and systematic approach to trace retinal vessel trees with the present of crossovers by solving a transductive learning problem on induced undirected graphs.

Automatic segmentation and measurement of vasculature in retinal fundus images using probabilistic formulation

The automatic analysis of retinal blood vessels plays an important role in the computer-aided diagnosis. In this paper, we introduce a probabilistic tracking-based method for automatic vessel segmentation in retinal images. We take into account vessel edge detection on the whole retinal image and handle different vessel structures. During the tracking process, a Bayesian method with maximum a posteriori (MAP) as criterion is used to detect vessel edge points. Experimental evaluations of the tracking algorithm are performed on real retinal images from three publicly available databases: STARE (Hoover et al., 2000), DRIVE (Staal et al., 2004), and REVIEW (Al-Diri et al., 2008 and 2009). We got high accuracy in vessel segmentation, width measurements, and vessel structure identification. The sensitivity and specificity on STARE are 0.7248 and 0.9666, respectively. On DRIVE, the sensitivity is 0.6522 and the specificity is up to 0.9710.

Computer-aided diagnosis of diabetic retinopathy: a review

Diabetes mellitus may cause alterations in the retinal microvasculature leading to diabetic retinopathy. Unchecked, advanced diabetic retinopathy may lead to blindness. It can be tedious and time consuming to decipher subtle morphological changes in optic disk, microaneurysms, hemorrhage, blood vessels, macula, and exudates through manual inspection of fundus images. A computer aided diagnosis system can significantly reduce the burden on the ophthalmologists and may alleviate the inter and intra observer variability. This review discusses the available methods of various retinal feature extractions and automated analysis.

Application of morphological bit planes in retinal blood vessel extraction

The appearance of the retinal blood vessels is an important diagnostic indicator of various clinical disorders of the eye and the body. Retinal blood vessels have been shown to provide evidence in terms of change in diameter, branching angles, or tortuosity, as a result of ophthalmic disease. This paper reports the development for an automated method for segmentation of blood vessels in retinal images. A unique combination of methods for retinal blood vessel skeleton detection and multidirectional morphological bit plane slicing is presented to extract the blood vessels from the color retinal images. The skeleton of main vessels is extracted by the application of directional differential operators and then evaluation of combination of derivative signs and average derivative values. Mathematical morphology has been materialized as a proficient technique for quantifying the retinal vasculature in ocular fundus images. A multidirectional top-hat operator with rotating structuring elements is used to emphasize the vessels in a particular direction, and information is extracted using bit plane slicing. An iterative region growing method is applied to integrate the main skeleton and the images resulting from bit plane slicing of vessel direction-dependent morphological filters. The approach is tested on two publicly available databases DRIVE and STARE. Average accuracy achieved by the proposed method is 0.9423 for both the databases with significant values of sensitivity and specificity also; the algorithm outperforms the second human observer in terms of precision of segmented vessel tree.

Retina verification system based on biometric graph matching

This paper presents an automatic retina verification framework based on the biometric graph matching (BGM) algorithm. The retinal vasculature is extracted using a family of matched filters in the frequency domain and morphological operators. Then, retinal templates are defined as formal spatial graphs derived from the retinal vasculature. The BGM algorithm, a noisy graph matching algorithm, robust to translation, non-linear distortion, and small rotations, is used to compare retinal templates. The BGM algorithm uses graph topology to define three distance measures between a pair of graphs, two of which are new. A support vector machine (SVM) classifier is used to distinguish between genuine and imposter comparisons. Using single as well as multiple graph measures, the classifier achieves complete separation on a training set of images from the VARIA database (60% of the data), equaling the state-of-the-art for retina verification. Because the available data set is small, kernel density estimation (KDE) of the genuine and imposter score distributions of the training set are used to measure performance of the BGM algorithm. In the one dimensional case, the KDE model is validated with the testing set. A 0 EER on testing shows that the KDE model is a good fit for the empirical distribution. For the multiple graph measures, a novel combination of the SVM boundary and the KDE model is used to obtain a fair comparison with the KDE model for the single measure. A clear benefit in using multiple graph measures over a single measure to distinguish genuine and imposter comparisons is demonstrated by a drop in theoretical error of between 60% and more than two orders of magnitude.

An approach to localize the retinal blood vessels using bit planes and centerline detection

The change in morphology, diameter, branching pattern or tortuosity of retinal blood vessels is an important indicator of various clinical disorders of the eye and the body. This paper reports an automated method for segmentation of blood vessels in retinal images. A unique combination of techniques for vessel centerlines detection and morphological bit plane slicing is presented to extract the blood vessel tree from the retinal images. The centerlines are extracted by using the first order derivative of a Gaussian filter in four orientations and then evaluation of derivative signs and average derivative values is performed. Mathematical morphology has emerged as a proficient technique for quantifying the blood vessels in the retina. The shape and orientation map of blood vessels is obtained by applying a multidirectional morphological top-hat operator with a linear structuring element followed by bit plane slicing of the vessel enhanced grayscale image. The centerlines are combined with these maps to obtain the segmented vessel tree. The methodology is tested on three publicly available databases DRIVE, STARE and MESSIDOR. The results demonstrate that the performance of the proposed algorithm is comparable with state of the art techniques in terms of accuracy, sensitivity and specificity.

Blood vessel segmentation methodologies in retinal images--a survey

Retinal vessel segmentation algorithms are a fundamental component of automatic retinal disease screening systems. This work examines the blood vessel segmentation methodologies in two dimensional retinal images acquired from a fundus camera and a survey of techniques is presented. The aim of this paper is to review, analyze and categorize the retinal vessel extraction algorithms, techniques and methodologies, giving a brief description, highlighting the key points and the performance measures. We intend to give the reader a framework for the existing research; to introduce the range of retinal vessel segmentation algorithms; to discuss the current trends and future directions and summarize the open problems. The performance of algorithms is compared and analyzed on two publicly available databases (DRIVE and STARE) of retinal images using a number of measures which include accuracy, true positive rate, false positive rate, sensitivity, specificity and area under receiver operating characteristic (ROC) curve.

EXTRACTION OF CORONARY ARTERIAL TREE USING CINE X-RAY ANGIOGRAMS

An efficient and robust method for identification of coronary arteries and evaluation of the severity of the stenosis on the routine X-ray angiograms is proposed. It is a challenging process to accurately identify coronary artery due to poor signal-to-noise ratio, vessel overlap, and superimposition with various anatomical structures such as ribs, spine, or heart chambers. The proposed method consists of two major stages: (a) signal-based image segmentation and (b) vessel feature extraction. The 3D Fourier and 3D Wavelet transforms are first employed to reduce the background and noisy structures in the images. Afterwards, a set of matched filters was applied to enhance the coronary arteries in the images. At the end, clustering analysis, histogram technique, and size filtering were utilized to obtain a binary image that consists of the final segmented coronary arterial tree. To extract vessel features in terms of vessel centerline and diameter, a gradient vector-flow based snake algorithm is applied to determine the medial axis of a vessel followed by the calculations of vessel boundaries and width associated with the detected medial axis.

A function for quality evaluation of retinal vessel segmentations

Retinal blood vessel assessment plays an important role in the diagnosis of ophthalmic pathologies. The use of digital images for this purpose enables the application of a computerized approach and has fostered the development of multiple methods for automated vascular tree segmentation. Metrics based on contingency tables for binary classification have been widely used for evaluating the performance of these algorithms. Metrics from this family are based on the measurement of a success or failure rate in the detected pixels, obtained by means of pixel-to-pixel comparison between the automated segmentation and a manually-labeled reference image. Therefore, vessel pixels are not considered as a part of a vascular structure with specific features. This paper contributes a function for the evaluation of global quality in retinal vessel segmentations. This function is based on the characterization of vascular structures as connected segments with measurable area and length. Thus, its design is meant to be sensitive to anatomical vascularity features. Comparison of results between the proposed function and other general quality evaluation functions shows that this proposal renders a high matching degree with human quality perception. Therefore, it can be used to enhance quality evaluation in retinal vessel segmentations, supplementing the existing functions. On the other hand, from a general point of view, the applied concept of measuring descriptive properties may be used to design specialized functions aimed at segmentation quality evaluation in other complex structures.

An automated blood vessel segmentation algorithm using histogram equalization and automatic threshold selection

This paper focuses on the detection of retinal blood vessels which play a vital role in reducing the proliferative diabetic retinopathy and for preventing the loss of visual capability. The proposed algorithm which takes advantage of the powerful preprocessing techniques such as the contrast enhancement and thresholding offers an automated segmentation procedure for retinal blood vessels. To evaluate the performance of the new algorithm, experiments are conducted on 40 images collected from DRIVE database. The results show that the proposed algorithm performs better than the other known algorithms in terms of accuracy. Furthermore, the proposed algorithm being simple and easy to implement, is best suited for fast processing applications.

Digital ocular fundus imaging: a review

Ocular fundus imaging plays a key role in monitoring the health status of the human eye. Currently, a large number of imaging modalities allow the assessment and/or quantification of ocular changes from a healthy status. This review focuses on the main digital fundus imaging modality, color fundus photography, with a brief overview of complementary techniques, such as fluorescein angiography. While focusing on two-dimensional color fundus photography, the authors address the evolution from nondigital to digital imaging and its impact on diagnosis. They also compare several studies performed along the transitional path of this technology. Retinal image processing and analysis, automated disease detection and identification of the stage of diabetic retinopathy (DR) are addressed as well. The authors emphasize the problems of image segmentation, focusing on the major landmark structures of the ocular fundus: the vascular network, optic disk and the fovea. Several proposed approaches for the automatic detection of signs of disease onset and progression, such as microaneurysms, are surveyed. A thorough comparison is conducted among different studies with regard to the number of eyes/subjects, imaging modality, fundus camera used, field of view and image resolution to identify the large variation in characteristics from one study to another. Similarly, the main features of the proposed classifications and algorithms for the automatic detection of DR are compared, thereby addressing computer-aided diagnosis and computer-aided detection for use in screening programs.

A new supervised method for blood vessel segmentation in retinal images by using gray-level and moment invariants-based features

This paper presents a new supervised method for blood vessel detection in digital retinal images. This method uses a neural network (NN) scheme for pixel classification and computes a 7-D vector composed of gray-level and moment invariants-based features for pixel representation. The method was evaluated on the publicly available DRIVE and STARE databases, widely used for this purpose, since they contain retinal images where the vascular structure has been precisely marked by experts. Method performance on both sets of test images is better than other existing solutions in literature. The method proves especially accurate for vessel detection in STARE images. Its application to this database (even when the NN was trained on the DRIVE database) outperforms all analyzed segmentation approaches. Its effectiveness and robustness with different image conditions, together with its simplicity and fast implementation, make this blood vessel segmentation proposal suitable for retinal image computer analyses such as automated screening for early diabetic retinopathy detection.

Automatic segmentation of coronary angiograms based on fuzzy inferring and probabilistic tracking

BACKGROUND

Segmentation of the coronary angiogram is important in computer-assisted artery motion analysis or reconstruction of 3D vascular structures from a single-plan or biplane angiographic system. Developing fully automated and accurate vessel segmentation algorithms is highly challenging, especially when extracting vascular structures with large variations in image intensities and noise, as well as with variable cross-sections or vascular lesions.

METHODS

This paper presents a novel tracking method for automatic segmentation of the coronary artery tree in X-ray angiographic images, based on probabilistic vessel tracking and fuzzy structure pattern inferring. The method is composed of two main steps: preprocessing and tracking. In preprocessing, multiscale Gabor filtering and Hessian matrix analysis were used to enhance and extract vessel features from the original angiographic image, leading to a vessel feature map as well as a vessel direction map. In tracking, a seed point was first automatically detected by analyzing the vessel feature map. Subsequently, two operators [e.g., a probabilistic tracking operator (PTO) and a vessel structure pattern detector (SPD)] worked together based on the detected seed point to extract vessel segments or branches one at a time. The local structure pattern was inferred by a multi-feature based fuzzy inferring function employed in the SPD. The identified structure pattern, such as crossing or bifurcation, was used to control the tracking process, for example, to keep tracking the current segment or start tracking a new one, depending on the detected pattern.

RESULTS

By appropriate integration of these advanced preprocessing and tracking steps, our tracking algorithm is able to extract both vessel axis lines and edge points, as well as measure the arterial diameters in various complicated cases. For example, it can walk across gaps along the longitudinal vessel direction, manage varying vessel curvatures, and adapt to varying vessel widths in situations with arterial stenoses and aneurysms.

CONCLUSIONS

Our algorithm performs well in terms of robustness, automation, adaptability, and applicability. In particular, the successful development of two novel operators, namely, PTO and SPD, ensures the performance of our algorithm in vessel tracking.

General retinal vessel segmentation using regularization-based multiconcavity modeling

Detecting blood vessels in retinal images with the presence of bright and dark lesions is a challenging unsolved problem. In this paper, a novel multiconcavity modeling approach is proposed to handle both healthy and unhealthy retinas simultaneously. The differentiable concavity measure is proposed to handle bright lesions in a perceptive space. The line-shape concavity measure is proposed to remove dark lesions which have an intensity structure different from the line-shaped vessels in a retina. The locally normalized concavity measure is designed to deal with unevenly distributed noise due to the spherical intensity variation in a retinal image. These concavity measures are combined together according to their statistical distributions to detect vessels in general retinal images. Very encouraging experimental results demonstrate that the proposed method consistently yields the best performance over existing state-of-the-art methods on the abnormal retinas and its accuracy outperforms the human observer, which has not been achieved by any of the state-of-the-art benchmark methods. Most importantly, unlike existing methods, the proposed method shows very attractive performances not only on healthy retinas but also on a mixture of healthy and pathological retinas.

Detection of blood vessels in ophthalmoscope images using MF/ant (matched filter/ant colony) algorithm

Blood vessels in ophthalmoscope images play an important role in diagnosis of some serious pathologies on retinal images. Hence, accurate extraction of vessels is becoming a main topic of this research area. Matched filter (MF) implementation for blood vessel detection is one of the methods giving more accurate results. Using this filter alone might not recover all the vessels (especially the capillaries). In this paper, a novel approach (MF/ant algorithm) is proposed to overcome the deficiency of the MF. The proposed method is a hybrid model of matched filter and ant colony algorithm. In this work, the accuracy and parameters of the hybrid algorithm are also discussed. The proposed method shows its success using the well known reference ophthalmoscope images of DRIVE database.

Unsupervised fuzzy based vessel segmentation in pathological digital fundus images

Performing the segmentation of vasculature in the retinal images having pathology is a challenging problem. This paper presents a novel approach for automated segmentation of the vasculature in retinal images. The approach uses the intensity information from red and green channels of the same retinal image to correct non-uniform illumination in color fundus images. Matched filtering is utilized to enhance the contrast of blood vessels against the background. The enhanced blood vessels are then segmented by employing spatially weighted fuzzy c-means clustering based thresholding which can well maintain the spatial structure of the vascular tree segments. The proposed method's performance is evaluated on publicly available DRIVE and STARE databases of manually labeled images. On the DRIVE and STARE databases, it achieves an area under the receiver operating characteristic curve of 0.9518 and 0.9602 respectively, being superior to those presented by state-of-the-art unsupervised approaches and comparable to those obtained with the supervised methods.

A model-based consecutive scanline tracking method for extracting vascular networks from 2-D digital subtraction angiograms

We propose a new model-based algorithm for the automated tracking of vascular networks in 2-D digital subtraction angiograms. Consecutive scanline profiles are fitted by a parametric imaging model to estimate local vessel center point, radius, edge locations and direction. An adaptive tracking strategy is applied with appropriate termination criteria to track each vessel segment. When tracking stops, to prevent premature termination and to detect bifurcations, a look ahead detection scheme is used to search for possible continuation points of the same vessel segment or those of its bifurcated segments. The proposed algorithm can automatically extract the majority of the vascular network without human interaction other than initializing the start point and direction. Compared to other tracking methods, the proposed method highlights accurate estimation of local vessel geometry. Accurate geometric information and a hierarchical vessel network are obtained which can be used for further quantitative analysis of arterial networks to obtain flow conductance estimates.

A novel method for the automatic grading of retinal vessel tortuosity

Tortuosity is among the first alterations in the retinal vessel network to appear in many retinopathies, such as those due to hypertension. An automatic evaluation of retinal vessel tortuosity would help the early detection of such retinopathies. Quite a few techniques for tortuosity measurement and classification have been proposed, but they do not always match the clinical concept of tortuosity. This justifies the need for a new definition, able to express in mathematical terms the tortuosity as perceived by ophthalmologists. We propose here a new algorithm for the evaluation of tortuosity in vessels recognized in digital fundus images. It is based on partitioning each vessel in segments of constant-sign curvature and then combining together each evaluation of such segments and their number. The algorithm has been compared with other available tortuosity measures on a set of 30 arteries and one of 30 veins from 60 different images. These vessels had been preliminarily ordered by a retina specialist by increasing perceived tortuosity. The proposed algorithm proved to be the best one in matching the clinically perceived vessel tortuosity.

A novel vessel segmentation algorithm for pathological retina images based on the divergence of vector fields

In this paper, a method is proposed for detecting blood vessels in pathological retina images. In the proposed method, blood vessel-like objects are extracted using the Laplacian operator and noisy objects are pruned according to the centerlines, which are detected using the normalized gradient vector field. The method has been tested with all the pathological retina images in the publicly available STARE database. Experiment results show that the method can avoid detecting false vessels in pathological regions and can produce reliable results for healthy regions.

Detecting false vessel recognitions in retinal fundus analysis

Automatic tracking of blood vessels in images of retinal fundus is an important and non-invasive procedure for the diagnosis of many diseases. Tracking techniques often present a high rate of false positives. This paper presents six methods to discriminate false detections from true positives, each based on a different model of the vessel. They describe a candidate vessel in terms of its average geometric and grayscale properties considered along the full trajectory of the vessel itself. The rationale is that false vessels are caused by the small scale of the tracking algorithm necessary during the tracking phase. Once tracking has been completed, we can gather information from the full vessel trajectory and solve ambiguities that cannot be fixed during tracking. We apply Fisher linear discriminant analysis to these features to get the desired discrimination. Results on 28 images show satisfactory rejection of false positives and better results when using more complex models.

Retinal blood vessel segmentation using line operators and support vector classification

In the framework of computer-aided diagnosis of eye diseases, retinal vessel segmentation based on line operators is proposed. A line detector, previously used in mammography, is applied to the green channel of the retinal image. It is based on the evaluation of the average grey level along lines of fixed length passing through the target pixel at different orientations. Two segmentation methods are considered. The first uses the basic line detector whose response is thresholded to obtain unsupervised pixel classification. As a further development, we employ two orthogonal line detectors along with the grey level of the target pixel to construct a feature vector for supervised classification using a support vector machine. The effectiveness of both methods is demonstrated through receiver operating characteristic analysis on two publicly available databases of color fundus images.

Texton-based segmentation of retinal vessels

With improvements in fundus imaging technology and the increasing use of digital images in screening and diagnosis, the issue of automated analysis of retinal images is gaining more serious attention. We consider the problem of retinal vessel segmentation, a key issue in automated analysis of digital fundus images. We propose a texture-based vessel segmentation algorithm based on the notion of textons. Using a weak statistical learning approach, we construct textons for retinal vasculature by designing filters that are specifically tuned to the structural and photometric properties of retinal vessels. We evaluate the performance of the proposed approach using a standard database of retinal images. On the DRIVE data set, the proposed method produced an average performance of 0.9568 specificity at 0.7346 sensitivity. This compares well with the best-published results on the data set 0.9773 specificity at 0.7194 sensitivity [Proc. SPIE5370, 648 (2004)].

Segmentation of retinal blood vessels using a novel clustering algorithm (RACAL) with a partial supervision strategy

In this paper, segmentation of blood vessels from colour retinal images using a novel clustering algorithm with a partial supervision strategy is proposed. The proposed clustering algorithm, which is a RAdius based Clustering ALgorithm (RACAL), uses a distance based principle to map the distributions of the data by utilising the premise that clusters are determined by a distance parameter, without having to specify the number of clusters. Additionally, the proposed clustering algorithm is enhanced with a partial supervision strategy and it is demonstrated that it is able to segment blood vessels of small diameters and low contrasts. Results are compared with those from the KNN classifier and show that the proposed RACAL performs better than the KNN in case of abnormal images as it succeeds in segmenting small and low contrast blood vessels, while it achieves comparable results for normal images. For automation process, RACAL can be used as a classifier and results show that it performs better than the KNN classifier in both normal and abnormal images.

Automated detection algorithm for arteriolar narrowing on fundus images

We have developed a computer-aided diagnosis system (CAD) to detect abnormalities in fundus images. In Japan, ophthalmologists usually detect hypertensive changes by identifying arteriolar narrowing and focal arteriolar narrowing. The purpose of this study is to develop an automated method for detecting arteriolar narrowing and focal arteriolar narrowing on fundus images. The blood vessel candidates were detected by the density analysis method. In blood vessel tracking, a local detection function was used to determine the centerline of the blood vessel. A direction comparison function using three vectors was designed to optimally estimate the next possible location of a blood vessel. After the connectivity of vessel segments was adjusted based on the recognized intersections, the true tree-like structure of the blood vessels was established. The blood vessels were recognized as arteries or veins by hue of HSV color space and their diameters. The arteriolar narrowing was detected by the ratio of diameters (artery vs. vein; A/V ratio). Focal arteriolar narrowing was detected by measuring the diameter of an artery. By applying this method to 100 fundus images, the detection sensitivity for arteriolar narrowing was found to be 76% when the specificity was 91%. Furthermore, by applying this method to 70 other different fundus images, the detection sensitivity for the focal arteriolar narrowing was 75% with 2.9 false positives per image. The number of some false positives is planned to be reduced during the next stage of development. Such an automated detection of abnormal vessels could help ophthalmologists in diagnosing ocular diseases.

A new tracking system for the robust extraction of retinal vessel structure

Identification and measurement of blood vessels in retinal images could allow quantitative evaluation of clinical features, which may allow early diagnosis and effective monitoring of therapies in retinopathy. A new system is proposed for the automatic extraction of the vascular structure in retinal images, based on a sparse tracking technique. After processing pixels on a grid of rows and columns to determine a set of starting points (seeds), the tracking procedure starts. It moves along the vessel by analyzing subsequent vessel cross sections (lines perpendicular to the vessel direction), and extracting the vessel center, calibre and direction. Vessel points in a cross section are found by means of a fuzzy c-means classifier. When tracking stops because of a critical area, e.g. low contrast, bifurcation or crossing, a "bubble technique" module is run. It grows and analyzes circular scan lines around the critical points, allowing the exploration of the vessel structure beyond the critical areas. After tracking the vessels, identified segments are connected by a greedy connection algorithm. Finally bifurcations and crossings are identified analyzing vessel end points with respect to the vessel structure. Numerical evaluation of the performances of the system compared to human expert are reported.

Segmentation of blood vessels from red-free and fluorescein retinal images

The morphology of the retinal blood vessels can be an important indicator for diseases like diabetes, hypertension and retinopathy of prematurity (ROP). Thus, the measurement of changes in morphology of arterioles and venules can be of diagnostic value. Here we present a method to automatically segment retinal blood vessels based upon multiscale feature extraction. This method overcomes the problem of variations in contrast inherent in these images by using the first and second spatial derivatives of the intensity image that gives information about vessel topology. This approach also enables the detection of blood vessels of different widths, lengths and orientations. The local maxima over scales of the magnitude of the gradient and the maximum principal curvature of the Hessian tensor are used in a multiple pass region growing procedure. The growth progressively segments the blood vessels using feature information together with spatial information. The algorithm is tested on red-free and fluorescein retinal images, taken from two local and two public databases. Comparison with first public database yields values of 75.05% true positive rate (TPR) and 4.38% false positive rate (FPR). Second database values are of 72.46% TPR and 3.45% FPR. Our results on both public databases were comparable in performance with other authors. However, we conclude that these values are not sensitive enough so as to evaluate the performance of vessel geometry detection. Therefore we propose a new approach that uses measurements of vessel diameters and branching angles as a validation criterion to compare our segmented images with those hand segmented from public databases. Comparisons made between both hand segmented images from public databases showed a large inter-subject variability on geometric values. A last evaluation was made comparing vessel geometric values obtained from our segmented images between red-free and fluorescein paired images with the latter as the "ground truth". Our results demonstrated that borders found by our method are less biased and follow more consistently the border of the vessel and therefore they yield more confident geometric values.

Retinal vessel centerline extraction using multiscale matched filters, confidence and edge measures

Motivated by the goals of improving detection of low-contrast and narrow vessels and eliminating false detections at nonvascular structures, a new technique is presented for extracting vessels in retinal images. The core of the technique is a new likelihood ratio test that combines matched-filter responses, confidence measures and vessel boundary measures. Matched filter responses are derived in scale-space to extract vessels of widely varying widths. A vessel confidence measure is defined as a projection of a vector formed from a normalized pixel neighborhood onto a normalized ideal vessel profile. Vessel boundary measures and associated confidences are computed at potential vessel boundaries. Combined, these responses form a six-dimensional measurement vector at each pixel. A training technique is used to develop a mapping of this vector to a likelihood ratio that measures the "vesselness" at each pixel. Results comparing this vesselness measure to matched filters alone and to measures based on the Hessian of intensities show substantial improvements, both qualitatively and quantitatively. The Hessian can be used in place of the matched filter to obtain similar but less-substantial improvements or to steer the matched filter by preselecting kernel orientations. Finally, the new vesselness likelihood ratio is embedded into a vessel tracing framework, resulting in an efficient and effective vessel centerline extraction algorithm.

Retinal vessel segmentation using the 2-D Gabor 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 nonvessel, based on the pixel's feature vector. Feature vectors are composed of the pixel's intensity and two-dimensional Gabor wavelet transform responses taken at multiple scales. The Gabor 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. 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 (Staal et al., 2004) and STARE (Hoover et al., 2000) databases of manually labeled images. On the DRIVE database, it achieves an area under the receiver operating characteristic curve of 0.9614, being slightly superior than that presented by state-of-the-art approaches. We are making our implementation available as open source MATLAB scripts for researchers interested in implementation details, evaluation, or development of methods.

Segmentation of retinal blood vessels by combining the detection of centerlines and morphological reconstruction

This paper presents an automated method for the segmentation of the vascular network in retinal images. The algorithm starts with the extraction of vessel centerlines, which are used as guidelines for the subsequent vessel filling phase. For this purpose, the outputs of four directional differential operators are processed in order to select connected sets of candidate points to be further classified as centerline pixels using vessel derived features. The final segmentation is obtained using an iterative region growing method that integrates the contents of several binary images resulting from vessel width dependent morphological filters. Our approach was tested on two publicly available databases and its results are compared with recently published methods. The results demonstrate that our algorithm outperforms other solutions and approximates the average accuracy of a human observer without a significant degradation of sensitivity and specificity.

Automatic grading of retinal vessel caliber

New clinical studies suggest that narrowing of the retinal blood vessels may be an early indicator of cardiovascular diseases. One measure to quantify the severity of retinal arteriolar narrowing is the arteriolar-to-venular diameter ratio (AVR). The manual computation of AVR is a tedious process involving repeated measurements of the diameters of all arterioles and venules in the retinal images by human graders. Consistency and reproducibility are concerns. To facilitate large-scale clinical use in the general population, it is essential to have a precise, efficient and automatic system to compute this AVR. This paper describes a new approach to obtain AVR. The starting points of vessels are detected using a matched Gaussian filter. The detected vessels are traced with the help of a combined Kalman filter and Gaussian filter. A modified Gaussian model that takes into account the central light reflection of arterioles is proposed to describe the vessel profile. The width of a vessel is obtained by data fitting. Experimental results indicate a 97.1% success rate in the identification of vessel starting points, and a 99.2% success rate in the tracking of retinal vessels. The accuracy of the AVR computation is well within the acceptable range of deviation among the human graders, with a mean relative AVR error of 4.4%. The system has interested clinical research groups worldwide and will be tested in clinical studies.