Automatic Identification of Retinal Arteries and Veins From Dual-Wavelength Images Using Structural and Functional Features

This paper presents an automated method to identify arteries and veins in dual-wavelength retinal fundus images recorded at 570 and 600 nm. Dual-wavelength imaging provides both structural and functional features that can be exploited for identification. The processing begins with automated tracing of the vessels from the 570-nm image. The 600-nm image is registered to this image, and structural and functional features are computed for each vessel segment. We use the relative strength of the vessel central reflex as the structural feature. The central reflex phenomenon, caused by light reflection from vessel surfaces that are parallel to the incident light, is especially pronounced at longer wavelengths for arteries compared to veins. We use a dual-Gaussian to model the cross-sectional intensity profile of vessels. The model parameters are estimated using a robust -estimator, and the relative strength of the central reflex is computed from these parameters. The functional feature exploits the fact that arterial blood is more oxygenated relative to that in veins. This motivates use of the ratio of the vessel optical densities (ODs) from images at oxygen-sensitive and oxygen-insensitive wavelengths () as a functional indicator. Finally, the structural and functional features are combined in a classifier to identify the type of the vessel. We experimented with four different classifiers and the best result was given by a support vector machine (SVM) classifier. With the SVM classifier, the proposed algorithm achieved true positive rates of 97% for the arteries and 90% for the veins, when applied to a set of 251 vessel segments obtained from 25 dual wavelength images. The ability to identify the vessel type is useful in applications such as automated retinal vessel oximetry and automated analysis of vascular changes without manual intervention.

[1]  Chia-Ling Tsai,et al.  The dual-bootstrap iterative closest point algorithm with application to retinal image registration , 2003, IEEE Transactions on Medical Imaging.

[2]  A.D. Hoover,et al.  Locating blood vessels in retinal images by piecewise threshold probing of a matched filter response , 2000, IEEE Transactions on Medical Imaging.

[3]  J. Beach,et al.  Hyperspectral imaging for measurement of oxygen saturation in the optic nerve head. , 2004, Investigative ophthalmology & visual science.

[4]  Theoretical Analysis of a Spectrophotometric Technique for Measuring Oxygen Saturation in Retinal Vessels , 2000 .

[5]  Langis Gagnon,et al.  Procedure to detect anatomical structures in optical fundus images , 2001, SPIE Medical Imaging.

[6]  Alfredo Ruggeri,et al.  A divide et impera strategy for automatic classification of retinal vessels into arteries and veins , 2003, Proceedings of the 25th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (IEEE Cat. No.03CH37439).

[7]  M. Goldbaum,et al.  Detection of blood vessels in retinal images using two-dimensional matched filters. , 1989, IEEE transactions on medical imaging.

[8]  Xiaohong W. Gao,et al.  A method of vessel tracking for vessel diameter measurement on retinal images , 2001, Proceedings 2001 International Conference on Image Processing (Cat. No.01CH37205).

[9]  Pascale Massin,et al.  A contribution of image processing to the diagnosis of diabetic retinopathy-detection of exudates in color fundus images of the human retina , 2002, IEEE Transactions on Medical Imaging.

[10]  J. Hickam,et al.  Studies of the Retinal Circulation in Man: Observations on Vessel Diameter, Arteriovenous Oxygen Difference, and Mean Circulation Time , 1966, Circulation.

[11]  William M. Campbell,et al.  Support vector machines for speaker verification and identification , 2000, Neural Networks for Signal Processing X. Proceedings of the 2000 IEEE Signal Processing Society Workshop (Cat. No.00TH8501).

[12]  Arthur Lompado,et al.  Multispectral confocal scanning laser ophthalmoscope for retinal vessel oximetry , 2000, BiOS.

[13]  Robert P. W. Duin,et al.  A Matlab Toolbox for Pattern Recognition , 2004 .

[14]  D. Roberts Analysis of vessel absorption profiles in retinal oximetry. , 1987, Medical physics.

[15]  Frans Vos,et al.  A model based method for retinal blood vessel detection , 2004, Comput. Biol. Medicine.

[16]  Massimiliano Pontil,et al.  Support Vector Machines for 3D Object Recognition , 1998, IEEE Trans. Pattern Anal. Mach. Intell..

[17]  K. Parker,et al.  Semiautomated computer analysis of vessel growth in preterm infants without and with ROP , 2003, The British journal of ophthalmology.

[18]  M. H. Smith,et al.  Oxygen Saturation Measurements of Blood in Retinal Vessels during Blood Loss. , 1998, Journal of biomedical optics.

[19]  L. Macaire,et al.  Fundus Images Segmentation by Unsupervised Classification , 1999 .

[20]  J. Beach,et al.  Oximetry of retinal vessels by dual-wavelength imaging: calibration and influence of pigmentation. , 1999, Journal of applied physiology.

[21]  Charles V. Stewart,et al.  Predictive scheduling algorithms for real-time feature extraction and spatial referencing: application to retinal image sequences , 2004, IEEE Transactions on Biomedical Engineering.

[22]  Shankar M. Krishnan,et al.  Detection and measurement of retinal vessels in fundus images using amplitude modified second-order Gaussian filter , 2002, IEEE Transactions on Biomedical Engineering.

[23]  Hong Shen,et al.  Optimal scheduling of tracing computations for real-time vascular landmark extraction from retinal fundus images , 2001, IEEE Transactions on Information Technology in Biomedicine.

[24]  Thorsten Joachims,et al.  Transductive Inference for Text Classification using Support Vector Machines , 1999, ICML.

[25]  Vladimir Vapnik,et al.  Statistical learning theory , 1998 .

[26]  M. Crawley Non‐linear Regression , 2007 .

[27]  Frédéric Zana,et al.  Segmentation of vessel-like patterns using mathematical morphology and curvature evaluation , 2001, IEEE Trans. Image Process..

[28]  Federico Girosi,et al.  Training support vector machines: an application to face detection , 1997, Proceedings of IEEE Computer Society Conference on Computer Vision and Pattern Recognition.

[29]  O. Chutatape,et al.  Retinal blood vessel detection and tracking by matched Gaussian and Kalman filters , 1998, Proceedings of the 20th Annual International Conference of the IEEE Engineering in Medicine and Biology Society. Vol.20 Biomedical Engineering Towards the Year 2000 and Beyond (Cat. No.98CH36286).

[30]  P J Derhaag,et al.  Quantitative analysis of retinal vascular changes in essential and renovascular hypertension , 1995, Journal of hypertension.

[31]  Joseph Jy-Haw Yu,et al.  Automatic Recognition Of Retinopathy From Retinal Images , 1990, [1990] Proceedings of the Twelfth Annual International Conference of the IEEE Engineering in Medicine and Biology Society.

[32]  Non - invasive retinal oximetry in normal human sub - jects , 2002 .

[33]  B. Roysam,et al.  Algorithms for automated oximetry along the retinal vascular tree from dual-wavelength fundus images. , 2005, Journal of biomedical optics.

[34]  R A Laing,et al.  Photographic measurements of retinal blood oxygen saturation: falling saturation rabbit experiments. , 1975, Investigative ophthalmology.

[35]  R. Klein,et al.  Retinal arteriolar narrowing and risk of coronary heart disease in men and women. The Atherosclerosis Risk in Communities Study. , 2002, JAMA.

[36]  Max A. Viergever,et al.  Ridge-based vessel segmentation in color images of the retina , 2004, IEEE Transactions on Medical Imaging.

[37]  E. Stefánsson,et al.  Laser treatment and the mechanism of edema reduction in branch retinal vein occlusion. , 2000, Investigative ophthalmology & visual science.

[38]  E. Stefánsson,et al.  Oxygenation and vasodilatation in relation to diabetic and other proliferative retinopathies. , 1983, Ophthalmic surgery.

[39]  Stephen R. Aylward,et al.  Initialization, noise, singularities, and scale in height ridge traversal for tubular object centerline extraction , 2002, IEEE Transactions on Medical Imaging.

[40]  Bjarne K. Ersbøll,et al.  Quantitative measurement of changes in retinal vessel diameter in ocular fundus images , 2000, Pattern Recognit. Lett..

[41]  Yunqian Ma,et al.  Practical selection of SVM parameters and noise estimation for SVM regression , 2004, Neural Networks.

[42]  Hong Shen,et al.  Rapid automated tracing and feature extraction from retinal fundus images using direct exploratory algorithms , 1999, IEEE Transactions on Information Technology in Biomedicine.

[43]  F. Delori Noninvasive technique for oximetry of blood in retinal vessels. , 1988, Applied optics.

[44]  J. Tukey,et al.  The Fitting of Power Series, Meaning Polynomials, Illustrated on Band-Spectroscopic Data , 1974 .

[45]  José A.F. Machado,et al.  Robust Model Selection and M-Estimation , 1993, Econometric Theory.

[46]  J. C. Ross,et al.  A Study of Retinal Venous Blood Oxygen Saturation in Human Subjects by Photographic Means , 1963, Circulation.

[47]  Chia-Ling Tsai,et al.  Automated Model-Based Segmentation, Tracing, and Analysis of Retinal Vasculature from Digital Fundus Images , 2003 .

[48]  Hideki Kuga,et al.  A computer method of understanding ocular fundus images , 1982, Pattern Recognit..

[49]  D. Schweitzer,et al.  In vivo measurement of the oxygen saturation of retinal vessels in healthy volunteers , 1999, IEEE Transactions on Biomedical Engineering.

[50]  Jae S. Lim,et al.  A new method for estimation of coronary artery dimensions in angiograms , 1988, IEEE Trans. Acoust. Speech Signal Process..

[51]  Tien Yin Wong,et al.  Retinal arteriolar narrowing and risk of diabetes mellitus in middle-aged persons. , 2002, JAMA.

[52]  Xiaoyi Jiang,et al.  Adaptive Local Thresholding by Verification-Based Multithreshold Probing with Application to Vessel Detection in Retinal Images , 2003, IEEE Trans. Pattern Anal. Mach. Intell..

[53]  Mong-Li Lee,et al.  A piecewise Gaussian model for profiling and differentiating retinal vessels , 2003, Proceedings 2003 International Conference on Image Processing (Cat. No.03CH37429).

[54]  Badrinath Roysam,et al.  Robust model-based vasculature detection in noisy biomedical images , 2004, IEEE Transactions on Information Technology in Biomedicine.

[55]  C. Riva,et al.  Hemoglobin oxygen saturation (So2) in the human ocular fundus measured by reflectance oximetry: preliminary data in retinal veins. , 2002, Klinische Monatsblatter fur Augenheilkunde.

[56]  Chris J. Wild,et al.  Computational Methods for Nonlinear Least Squares , 2005 .

[57]  Liang Zhou,et al.  The detection and quantification of retinopathy using digital angiograms , 1994, IEEE Trans. Medical Imaging.

[58]  P. Holland,et al.  Robust regression using iteratively reweighted least-squares , 1977 .

[59]  J S Tiedeman,et al.  Retinal oxygen consumption during hyperglycemia in patients with diabetes without retinopathy. , 1998, Ophthalmology.

[60]  G Valli,et al.  An algorithm for real-time vessel enhancement and detection. , 1997, Computer methods and programs in biomedicine.