International Journal of Biomedical Imaging
Volume 2018 (2018), Article ID 2815163, 7 pages
https://doi.org/10.1155/2018/2815163

Research Article

Localizing Optic Disc in Retinal Image Automatically with Entropy Based Algorithm

Lamia AbedNoor Muhammed

Computer Sciences and IT College, University of Al-Qadisiyah, Al-Qadisiyah, Iraq

Correspondence should be addressed to Lamia AbedNoor Muhammed; lamia.abed@qu.edu.iq

Received 13 September 2017; Revised 17 December 2017; Accepted 10 January 2018; Published 6 February 2018

Academic Editor: Guowei Wei

Copyright © 2018 Lamia AbedNoor Muhammed. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Examining retinal image continuously plays an important role in determining human eye health; with any variation present in this image, it may be resulting from some disease. Therefore, there is a need for computer-aided scanning for retinal image to perform this task automatically and accurately. The fundamental step in this task is identification of the retina elements; optical disk localization is the most important one in this identification. Different optical disc localization algorithms have been suggested, such as an algorithm that would be proposed in this paper. The assumption is based on the fact that optical disc area has rich information, so its entropy value is more significant in this area. The suggested algorithm has recursive steps for testing the entropy of different patches in image; sliding window technique is used to get these patches in a specific way. The results of practical work were obtained using different common data set, which achieved good accuracy in trivial computation time. Finally, this paper consists of four sections: a section for introduction containing the related works, a section for methodology and material, a section for practical work with results, and a section for conclusion.

1. Introduction

Glaucoma is a chronic eye disease, which can be controlled but cannot be cured. If left untreated, loss of vision occurs gradually, potentially leading to blindness. The detection and diagnosis of glaucoma are related to tracing the changes in the optic cup which is a portion of optic disc (OD). To perform this detection in an automated system, the optical disc region must be extracted from retinal image through segmentation process. However, localization of optic disc is an important step in simplifying this segmentation. Different methods have been proposed for localization of OD [1]. They have exploited some of OD region features, such as its yellow color, having more brightness, having high grey intensity, and containing a network of convergence vessels. They have been applied through different techniques.

Simple operations were adopted through the works such as Akram et al. [2] who used average filtering and threshold. Aquino et al. [3] exploited morphological operations and edge detection. Whardana and Suciati [4] combined two techniques: morphological operators and clustering with -means method. Li and Chutatape [5] presented an algorithm that uses the clustering of image according to the bright pixel and the candidate regions are passed through principal component analysis (PCA) in order to locate the center of optic disc region. Padmanaban and Kannan [6] suggested the use of Fuzzy C-Means (FCM) clustering. Foracchia et al. [7] worked on tracing the vessels, matching their path with directional pattern in OD in originate image. Nergiz et al. [8] introduced a study using vasculature geometry character in optical disc, convergence to its center. Mendonace et al. [9] presented new methodology: entropy based on the information resulting from the distribution of vessel through the optic disc region.

Learning techniques has good opportunity in optic disc localization fields; however, OD template would be learned based on different features. While Wu et al. [10] worked to build a specific model for network vessels shape in OD region which has form as parabolic shape, Dehghani et al. [11] suggested building a histogram-template for colors which OD region contains. Akyol et al. [12] proposed an algorithm which consists of a multiple-steps algorithm, with induction classifier being one of these steps. The study by Ichim and Popescu [13] used adaptive local texture analysis to generate several features that pass through classification algorithm. Muangnak et al. [14] used in their work decision model, which is induced from direction vectors that were derived from the vessel network, and points of convergence and then used hybrid method. Sinthanayothin et al. [15] applied neural network in their work with data input which is extracted from principle components analysis of the image in question.

Frequency transform techniques have been used in other studies. Pallawala et al. [16] proposed using Daubechies wavelet transform. Jafariani and Tabatabaee [17] employed Fourier-Mellin transform in their work. While another study produced by Esmaeili et al. [18] used curvelete transform technique.

2. Materials and Methods

2.1. Optic Disc

Eye fundus plays an important role in detecting various eye diseases; it consists of retina which is the transparent, light-sensitive structure at the back of the eye, optic nerve disc, and blood vessels (retinal arteries and veins). Anatomically, the retina contains structures including macula which is the central area, rods which are photoreceptor cells that surround the macula, the optic nerve which carries signals, and blood vessels as shown in Figure 1 [19].

Figure 1: Retinal fundus image [19].

One of the problems that eye suffers is the vision loss resulting from the diabetes complications that can be noticed through changes occurring in the retina components. In relation to this, fundus image is used continuously to examine and assess the eye health conditions [2].

Identification of optic disc in fundal landmark is important as reference coordinates to locate anatomical components in retinal images, for vessel tracking as a reference length for measuring distances in retinal images, and for observing any change in the optic disc which may result from a disease [21].

With the development of digital imaging and computing power, the potential to use these technologies in ophthalmology analysis and computer vision techniques also increases. Optic disk localization with this technique has been performed successfully according to its features that can be used in image analysis such as brightness, high contrast, and yellowish disk, where the blood vessels and optic nerves pass through it [2]. Therefore optic disc area contains more detailed information that can be exploited in its identification.

2.2. Entropy

An information theory is associated with information measure that has played an important role in different application fields, such as image analysis. The entropy of a probability distribution can be interpreted not only as a measure of uncertainty but also as a measure of information. As a matter of fact, the amount of information acquired from the observation of the result of an experiment (depending on chance) can be taken numerically equal to the amount of uncertainty concerning the outcome of the experiment before carrying it out [22].

Formally, let be a discrete random variable with alphabet and probability mass function , . The Shannon entropy of is defined as [23]

where ,  (−log ) is an information association for occurrence , and .

2.3. Entropy in Image Analysis

A digital image consists of small units that represent brightness of a particular position in the image, called picture elements (pixels). The variation of the pixels values carries information in an image, so it can be measured by entropy metric. However, the image has different distributed brightness values; entropy of image () can be computed as

where is the distribution of brightness value indexed by in image ; is total of brightness levels in an interested image ().

In image processing, entropy measure generates a value which represents new feature that can be exploited in image analysis such as texture analysis. Low values of entropy refer to smoothing texture, while texture with more details has higher entropy values [24]. Thus, it can be used to generate a new feature to measure the smoothing of the texture of images.

As related, the optic disc area in fundus image contains more details such as nerves and vessels passing through that means its texture is not smooth, so it is expected that its entropy value would be higher than other regions in the fundus image.

2.4. Sliding Window Approach

This is an approach that has been used in locating object in an image. However the image is partitioned into subimages (regions); they will be evaluated separately according to quality function. So the region with maximum score will be candidate to be Region of Interest [25].

Definition 1. Let be an image that would be partitioned into regions ():
.
However, is a region labeled with (top, bottom, left, and right) coordinates.
is quality function; is the candidate region with

The proposed algorithm in this paper used sliding window approach with two methods: nonoverlap and overlap. However, in nonoverlap method as shown in Figure 2, each image pixel must not belong to more than one region, so

Meanwhile, overlap method as shown in Figure 3 permits image pixel to belong to more than one region.

, are regions in the image.

Figure 2: Nonoverlap sliding window.

Figure 3: Overlap sliding window.

3. Practical Work and Results

3.1. Proposed Algorithm

As mentioned in the last section, there are some concepts that would be exploited in this work as shown in Algorithm 1.

(1)However, entropy value in the optic disc area is significant; it would be used to find this area. So the algorithm uses searching technique; greedy method was used in this work to find the maximum entropy value.(2)Sliding window technique was used in order to find specific area from the whole image; the image area would be partitioned into patches. This technique was executed in two methods: nonoverlap sliding window in order to find the approximate optic disc area and then using this area with another method; overlap sliding window which is of importance in order to find the location of optic disc accurately.(3)Finally, we used green channel image in this work according to results from previous related experimental work about this image. It seems from subjective information about image channels (red, green, and blue) that the best one is green channel image; however it gives maximum contrast [26]. Therefore its texture is more distinguished than other channel images as shown in Figure 4.

Input: Retina Image

Output: Retina Image with marking the Optic Disc location

Step: Extract Green channel Image ()

Step: Divide Image () into () patches ()

Stepfor each patch

        (i) Compute entropy value

StepPick the patchwith the maximum value

StepFormulate new regionthat is:

StepDivideinto () overlap patches ()

Stepfor each patch

        (i) Compute entropy value

StepPick the patch with maximum value

StepMarking the center location of this patch with small square into produce new image

Algorithm 1: Proposed algorithm: localization of optic disk.

(a) Original color image

(b) Red channel image

(c) Green channel image

(d) Blue channel image

Figure 4: Original image: (a) colors image and (b) red, (c) green, and (d) blue channel image.

3.2. Experimental Work

The proposed algorithm was executed computationally using the following tools: MATLAB R2014a programming language and computer with specification Intel®Core™i5-3320M CPU@2.60 GHz; tests were applied with the known data sets: DRIVE, CHASEDB, DRIONS-DB, and DIARETDBI. In addition, a special data set is used [20], whose results will be given in detail in the following paragraphs.

The algorithm was initiated with extraction of green channel image () which was extracted from original image () as shown in Figure 4. So image () was divided into () patches () using nonoverlap sliding window as shown in Figure 5(a). The size of each patch is the same; it is approximate to size of optic disc area. From the experimental work here, this size is suitable for implying more of optic disc area in one of the patches, so significant patch that may contain optic disc comes out. Then, for each patch , entropy value would be generated as shown in Figure 5(b).

Figure 5: Image number (29) in data set with its nonoverlap patches entropies.

The patches are ranked according to their entropy values. The one with the highest entropy value would be selected as the area that may contain optic disc, as shown in Figure 5(b). The darker box referring to the highest value which is corresponding to the patch is addressed by as shown in Figure 6(a). New region () would be generated from patch , and half of each surrounded patch of . However, the existence of these halves in a new patch guarantees the possibility that some optic disc area may be in one of the surrounded patches. Then would be divided into overlap patches as (), (), or (). From the experimental work, this number of divisions is suitable; however, size of new patches is approximately equal to the size of optic disc area as shown in Figure 6(b). The new patches entropies would be measured. The patch with maximum value would be candidate to be the optic disc area. From the results shown in Figure 6(b), the patch that is indexed () has the highest entropy value. Finally, small black square is placed in the center location of this patch as shown in Figure 7.

(a) Candidate patch with maximum entropy

(b) Entropy values for overlap patches generated from region in (a)

Figure 6: New region () with its overlap patches entropies.

Figure 7: Original image with marking optic disc location.

3.3. Results and Findings

As related to data set [20] results, the computing entropies for candidate patch that is closer to optic disk area are more significant than other patches in image; this is true for all the images (35 images) as shown in Figure 8. However there is a comparison between candidate patch entropy value and the mean entropy values of the rest of the patches in an image; each one is represented by different line. It is obvious that candidate patch entropy value exceeds the mean patches entropies, for all tested images. So these images would be marked correctly. Sample of these marked images is shown in Figure 9. In addition, the experimental work would be applied with other data set and two measures were computed for comparison the results: accuracy of localization of OD in each data set image and execution time as shown in Table 1. There is a disparity in the accuracies values of OD localization, but the time execution is stable for different tested data because the executed steps are constant.

Table 1: Comparison of the results for different data sets.

Data set name
Measure

Number of images
Accuracy
Average time of execution (millisecond)

CHASEDB ()
28
75%
37.641

DRIONS-DB
110
77%
44.305

DRIVE
40
95%
36.9

DIARETDBI
89
89%
38.12

Special data set [20]
35
100%
37.955

Figure 8: Comparison of the entropy value of candidate patch versus the entropy mean of all patches in each of the 35 tested images from source [20].

Figure 9: Sample output marked images.

4. Conclusion

Researchers have focused on localization of optic disk using a computer, so more methods have been suggested with good results. In spite of this, the proposed algorithm in this paper gives significant results through simple computational steps that are executed in a short space of time; however there is no need for any preprocessing enhancement steps. Moreover, the proposed algorithm is attractive; it entails a simple technique, so it can be combined with other algorithms in order to be more effective as a supporting step.

Finally, there are more features that are significant in the OD area, so in future work combining them with the entropy feature in order to generate robust algorithm is suggested. Also the evolutionary algorithm can enhance this algorithm to give accurate localization.

Conflicts of Interest

The author declares that there are no conflicts of interest.

References

R. P. Apeksha and L. Bhailal, “A Review Paper on Detection of Optic Disc Damage using Retinal Images,” International Journal of Computer Applications (0975 – 8887), vol. 111, no. 13, pp. 1–4, 2015.

M. U. Akram, A. Khan, K. Iqbal, and W. H. Butt, “Retinal Images: Optic Disk Localizationand Detection,” in Proceedings of the International Conference Image Analysis and Recognition (ICIAR '10), vol. LNCS 6112, Part II, pp. 40–49, Springer-Verlag Berlin Heidelberg, Berlin, Germany, 2010.

A. Aquino, M. E. Gegúndez-Arias, and D. Marín, “Detecting the optic disc boundary in digital fundus images using morphological, edge detection, and feature extraction techniques,” IEEE Transactions on Medical Imaging, vol. 29, no. 11, pp. 1860–1869, 2010.

A. K. Whardana and N. Suciati, “A Simple Method for Optic Disk Segmentation from Retinal Fundus Image,” International Journal of Image, Graphics and Signal Processing, vol. 11, pp. 36–42, 2014.

H. Li and O. Chutatape, “Automatic location of optic disk in retinal images,” in Proceedings 2001 International Conference on Image Processing (Cat. No.01CH37205), vol. 2, pp. 837–840, IEEE, Thessaloniki, Greece, 2001.

K. Padmanaban and R. J. kannan, “Localization of Optic Disc Using Fuzzy C Means Clustering,” in International Conference on Current Trends in Engineering and Technology, ICCTET'13, Technical Sponsor, pp. 571–573, IEEE, Coimbatore, India, 2013.

M. Foracchia, E. Grisan, and A. Ruggeri, “Detection of optic disc in retinal images by means of a geometrical model of vessel structure,” IEEE Transactions on Medical Imaging, vol. 23, no. 10, pp. 1189–1195, 2004.

M. Nergiz, M. Akin, and S. Ari, “A fast Optic Disc Localization Alogorithm based on Geometry of vessel distribution,” in Proceedings of the International Conference on Natural Science and Engineering (ICNSE'16), kilis, Turkey, 2016.

A. M. Mendonace, A. Sousaand, L. Mendonca, and A. Campilho, “Automatic localization of the optic disc by combining vascular and intensity information,” Computerized Medical Imaging and Graphics, vol. 37, no. 5-6, pp. 409–417, 2013.

X. Wu, B. Dai, and W. Bu, “Optic disc localization using directional models,” IEEE Transactions on Image Processing, vol. 25, no. 9, pp. 4433–4442, 2016.

A. Dehghani, H. A. Moghaddam, and MS. J Moin, “Optic disc localization in retinal images using histogram matching,” EURASIP Journal on Image and Video Processing, vol. 19, 2012.

K. Akyol, B. Şen, and Ş. Bayir, “Automatic Detection of Optic Disc in Retinal Image by Using Keypoint Detection, Texture Analysis, and Visual Dictionary Techniques,” Computational and Mathematical Methods in Medicine, vol. 2016, Article ID 6814791, 10 pages, 2016.

L. Ichim and D. Popescu, “Optic Disc Localization Based on Feature Sorting,” Memoirs of the Scientific Sections of the Romanian Academy Tome XXXIX, pp. 63–70, 2016.

N. Muangnak, P. Aimmanee, and S. Makhanov, “Automatic optic disk detection in retinal images using hybrid vessel phase portrait analysis,” Medical & Biological Engineering & Computing, pp. 1–16, 2017.

C. Sinthanayothin, J. F. Boyce, H. L. Cook, and T. H. Williamson, “Automated localisation of the optic disc, fovea, and retinal blood vessels from digital colour fundus images,” Journal of Medical Genetics, vol. 83, no. 8, http://bjo.bmj.com/.

P. M. D. S. Pallawala, W. Hsu, M. L. Lee, and K. A. Eong, “Automated optic disc localization and contour detection using ellipse fitting and wavelet transform,” in Computer Vision - ECCV 2004, T. Pajdla and J. Matas, Eds., vol. 3022 of Lecture Notes in Computer Science, pp. 139–151, Springer, Berlin, Germany, 2004.

H. Jafariani and H. Tabatabaee, “Retinal identification system using fourier-mellin transform and fuzzy clustering,” Indian Journal of Science and Technology, vol. 7, no. 9, pp. 1289–1296, 2014.

M. Esmaeili, H. Rabbani, and A. M. Dehnavi, “Automatic optic disk boundary extraction by the use of curvelet transform and deformable variational level set model,” Pattern Recognition, vol. 45, no. 7, pp. 2832–2842, 2012.

S. J. Garg, http://www.merckmanuals.com/home/eye-disorders/retinal-disorders/overview-of-retinal-disorders.

http://misp.mui.ac.ir/en/Color%20Fundus%20Images%20with%20Exudates.

N. Patton, T. M. Aslam, T. MacGillivray et al., “Retinal image analysis: concepts, applications and potential,” Progress in Retinal and Eye Research, vol. 25, no. 1, pp. 99–127, 2006.

A. Renyi, “On measures of entropy and information,” in Proceedings of the Fourth, Berkeley Symposium, pp. 547–561, University of California Press, Berkeley, calif, USA, 1961.

C. Wang and H. Shen, “Information theory in scientific visualization,” Entropy, vol. 13, pp. 254–273, 2011.

S. Theodoridis and K. Kourtroumbas, Pattern Recognition, Acadimemic Press, 4th edition.

C. H. Lampert, M. B. Blaschko, and T. Hofmann, “Efficient subwindow search: A branch and bound framework for object localization,” IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 31, no. 12, pp. 2129–2142, 2009.

M. Abdullah, M. M. Fraz, and S. A. Barman, “Localization and segmentation of optic disc in retinal images using circular Hough transform and grow-cut algorithm,” PeerJ, vol. 4, no. e2003, 2016, http://doi.org/10.7717/peerj.2003.