Detection of optic disc in retinal images by means of a geometrical model of vessel structure

We present here a new method to identify the position of the optic disc (OD) in retinal fundus images. The method is based on the preliminary detection of the main retinal vessels. All retinal vessels originate from the OD and their path follows a similar directional pattern (parabolic course) in all images. To describe the general direction of retinal vessels at any given position in the image, a geometrical parametric model was proposed, where two of the model parameters are the coordinates of the OD center. Using as experimental data samples of vessel centerline points and corresponding vessel directions, provided by any vessel identification procedure, model parameters were identified by means of a simulated annealing optimization technique. These estimated values provide the coordinates of the center of OD. A Matlab(R) prototype implementing this method was developed. An evaluation of the proposed procedure was performed using the set of 81 images from the STARE project, containing images from both normal and pathological subjects. The OD position was correctly identified in 79 out of 81 images (98%), even in rather difficult pathological situations.