Citation: Cao K, Xu J, Zhao WQ. Artificial intelligence on diabetic
retinopathy diagnosis: an automatic classification method based on grey level
co-occurrence matrix and naive Bayesian model. Int J Ophthalmol 2019;12(7):1158-1162
DOI:10.18240/ijo.2019.07.17
·Clinical
Research·
Artificial intelligence on diabetic retinopathy diagnosis: an automatic classification method based on grey level co-occurrence matrix and naive Bayesian model
Kai Cao, Jie Xu, Wei-Qi Zhao
Beijing Institute of Ophthalmology, Beijing Tongren Hospital of Capital Medical University, Beijing 100005, China
Correspondence to: Kai Cao. Beijing Institute of Ophthalmology, Beijing Tongren Hospital of Capital Medical University, Beijing 100005, China. anzhen602@163.com
Received:
Abstract
AIM: To develop an automatic tool on screening diabetic retinopathy (DR) from diabetic patients.
METHODS: We extracted textures from eye fundus images of each diabetes subject using grey level co-occurrence matrix method and trained a Bayesian model based on these textures. The receiver operating characteristic (ROC) curve was used to estimate the sensitivity and specificity of the Bayesian model.
RESULTS: A total of 1000 eyes fundus images from diabetic
patients in which 298 eyes were diagnosed as DR by two ophthalmologists. The
Bayesian model was trained using four extracted textures including contrast,
entropy, angular second moment and correlation using a training dataset. The
Bayesian model achieved a sensitivity of 0.949 and a specificity of
CONCLUSION: Textures extracted by grey level co-occurrence can be useful information for DR diagnosis, and a trained Bayesian model based on these textures can be an effective tool for DR screening among diabetic patients.
KEYWORDS: grey level co-occurrence matrix; Bayesian; textures; artificial intelligence; receiver operating characteristic curve; diabetic retinopathy
DOI:10.18240/ijo.2019.07.17
Citation: Cao K, Xu J, Zhao WQ. Artificial intelligence on diabetic retinopathy diagnosis: an automatic classification method based on grey level co-occurrence matrix and naive Bayesian model. Int J Ophthalmol 2019;12(7):1158-1162
INTRODUCTION
Diabetic retinopathy (DR) is a common microvascular complication of diabetes, and is one of the main causes of preventable blindness in population within working age[1, 2]. According to an analysis on global causes of blindness, DR ranks the 7th leading reason and takes responsibility for 0.85% of all blindness worldwide[3]. As the life expectancy getting extended, more and more diabetic patients are diagnosed. China alone has 92.4 million diabetic patients[4]. Meanwhile, DR population is up to 15.2 million[2].
Such large crowd of DR patients poses a challenge to ophthalmologists, because reading eye fundus takes lots of time. Doctors are in sore need of an automatic tool to help them do preliminary screening, and they can save time for the intractable cases.
Recently, many researchers cast their eyes over recognition filed, in this field many machine learning or deep learning methods are widely used to make classification of diseases based on images, and neural network algorithm is one the most commonly used method[5, 6, 7]. Other algorithms like support vector machine[8], random forest method etc. are also have been applied and showed good performance[9].
Naive bayes is also one kind of machine learning algorithm. Unlike other classifiers, naive bayes is a classification algorithm based on probability theory. Naive bayes is simple but efficient, and has been widely used in many fields like cardiovascular disease risk’s level detection[10], biomarker selection and classification from genome-wide SNP data[11], diagnosis of Parkinson’s disease[12], glaucoma management[13] and so on. What is more, at the early days of DR, symptoms like retinal hard exudates, microaneurysm and cotton version are very few and very hard to identify macroscopically. Some patients may be mis-diagnosed and thus miss the chance to stop the progress. However, we all know that DR could be stopped if treated by laser photocoagulation at the early stage[14]. Under this circumstance, early detection and early diagnosis is quite important. Some researchers even tried to find out some biomarkers of DR to help with early detection[15, 16, 17], but compared with auto-detection through eye fundus using artificial intelligence, biomarkers cost much more money and is of low efficiency[18].
In our study, we used grey-level co-occurrence matrix firstly to help extract textures from eye fundus images. Grey-level co-occurrence matrix method is a useful method on analyzing textural features and was originally proposed by Haralick[19], grey-level cooccurrence matrix method has been widely used in many fields[20, 21, 22]. Specifically, for DR, textures of symptoms like microaneurysm, exudation etc. are quite different from textures of normal images, so we can use the texture information to help make diagnosis. In our study, we applied naive Bayesian method to train a prediction model based on these textures. At last we used an external group of diabetic subjects’ eye fundus images to make validation on the trained model.
SUBJECTS AND METHODS
Ethical Approval The study followed the tenets of the Declaration of Helsinki. Oral consent was obtained from each patient, without any financial compensation.
All eye fundus images are from diabetic patients at Outpatient in Beijing Tongren Hospital during 2017 January to 2017 August. Diabetic patients were diagnosed by doctors or they reported a diabetic history and were using oral hypoglycemic drugs. All eye fundus pictures were taken with a Topcon TRC-NW6S/7S (Topcon; www.global.topcon.com) camera, and dilated 45° digital stereoscopic color retinal photographs of Early Treatment DR standard field 1 (centred on the optic disc, stereoscopic) and standard field 2 (centred on the macula, non-stereoscopic) were taken for each eye by trained and certified photographers. The eye fundus images were in a size of 1536 pixel × 1248 pixel. The Bayesian model was trained based on readable fundus images, blurry images due to cataract or other diseases were excluded.
Gray Level Co-occurrence Matrix The gray level co-occurrence matrix (GLCM) is defined as the probability that the gray value at a point away from a certain fixed position. We use Pδ to represent gray level co-occurrence, Pδ = L×L (L represents the number of different grey levels in a picture). Thus actually, Pδ is the frequency of occurrences of two pixels (i and j) with spatial positional relationship and gray levels.
For intuitive interpretation, we draw Figure 1. There are 3 grey levels in Figure 1.
Figure 1 Grey level co-occurrence matrix The number shows the grey level.
Thus we get a 3×3 matrix
Pδ= 
After normalization, the matrix will be
Pδ= 
Characteristic Quantities of Matrix The GLCM describes textural features with the spatial relationships among the gray levels through examining all pixels within the image space and counting the gray-level dependencies between two pixels separated by a certain interval distance along a particular direction. Haralick et al[23] described 14 statistics that can be calculated from the co-occurrence matrix. However, Ulaby et al[24] pointed out that only 4 statistics (contrast, angular second moment, entropy and correlation) are uncorrelated. In our study, we used four features to make diagnosis of DR.
Contrast Contrast reflects[23] how the matrix is distributed and to what degree the local change is. It shows the sharpness of the image and the depth of the texture. The deeper the groove from the gray texture is, the greater contrast and the clearer visual effect will be. Conversely, a small contrast comes with a shallow groove and a blurred vision effect.
Con=
(1)
Angular Second Moment Angular Second Moment (ASM) is also called energy. ASM reflects the uniformity of image gray distribution and the thickness of the texture. If the element values of the GLCM were similar, the energy would be small; if some of the values were large, and other values were small, then the energy value would be large.
ASM=
(2)
Entropy The image contains a measure of the randomness of the amount of information. When all the values in the co-occurrence matrix are equal or the pixel values exhibit the greatest randomness, the entropy is the largest. The entropy value indicates the complexity of the gray distribution of the image. The larger the entropy value is, the more complex the image will be[25].
Ent=
(3)
Correlation Correlation is used to measure the similarity of the gray level of the image in the row or column direction. So it is worth to reflect the local gray correlation. The larger the value, the greater the correlation.
COR= 
(4)
μi and μj are respectively the means for the rows and columns of the normalized matrix, and σi and σj are the corresponding standard deviations.
=
(5)
=
(6)
=
(7)
=
(8)
Image Preparation Before extracting the texture feature, we need to convert the multi-band image (RGB image) into a grayscale image. Because the texture feature is a structural feature, thus the texture features obtained from each band is the same. We arbitrarily chose the G band in our study. The image preparation process is shown in Figure 2.
Figure 2 Figure preparation To save time on computation, we compressed the grey levels of images from 256 levels to 16 levels which helped make the matrix smaller, with no influence on the texture extracting. Then we cut each eye fundus image into small pieces by a unit of 128×128 pixels. For each small piece, we extracted four textures from four directions: 0°, 45°, 90°, 135°, and then calculated the mean value of each texture. A: Original eye fundus; B: G band of eye fundus.
Quality Control All DR patients are diagnosed by two ophthalmologists after training. Although GLCM method doesn’t require the ophthalmologists to mark DR-related symptoms on the eye fundus images, the ophthalmologists were still asked to offer a brief description of diagnostic basis. When there is divergence, the eye fundus images will be sent to a third senior ophthalmologist to make final decision. We calculated the kappa coefficient to evaluate the consistency between two ophthalmologists.
Statistical
Analysis We use mean value and standard
deviation to make statistical description on continuous variables, and we use
frequency and percentage to make statistical description on categorical
variables. A naive Bayesian model was applied to train prediction model and a
10-fold cross validation method was applied to evaluate the accuracy of the
model. We used receiver operating characteristic (ROC) curve to make evaluation
on the classification ability of our built model. All the analysis was done
with open source R-program (version
RESULTS
Among all the included 1000 diabetic subjects, there are 298 DR patients and 702 control subjects. Specifically, among DR patients, there were 253 non-proliferative DR cases, 35 proliferative DR cases, and 10 diabetic macular edema cases. Two ophthalmologists misdiagnosed 49 and 25 cases respectively, and the kappa value (Table 1) between two doctors is 0.817 (95%CI: 0.776, 0.857). We think the diagnosis is trustable.
Table 1 Agreement on diabetic retinopathy diagnose between two doctors
Doctor 2 |
Doctor 1 |
Kappa (95%CI) |
|
No |
Yes |
||
No |
683 |
42 |
0.817 (0.776, 0.857) |
Yes |
32 |
243 |
|
We firstly
split all cases into two datasets by a probability of 7:3. The dataset with
larger sample size was used to train Bayesian model and the smaller dataset was
used to test the Bayesian model. As shown in Table 2, the Bayesian model
performed quite good. We got a sensitivity of 0.949 and a specificity of
Table 2 Cross table of predicted diagnosis and real diagnosis
Predicted |
Real |
Sensitivity |
Specificity |
|
No |
Yes |
|||
No |
192 |
15 |
0.949 |
0.928 |
Yes |
4 |
74 |
||
Figure 3 Receiver operating characteristic curve.
At last, we further used a 10-fold cross validation method to evaluate the accuracy of the Bayesian model, and we got an average accuracy rate of 93.5%.
DISCUSSION
DR has been rising globally. It is estimated that in 1990s DR contributed to 0.85% of all blindness worldwide while in 2015 this percentage reached 1.05%[3].
DR in the early days is reversible while at the late stage, it is irreversible[26, 27]. Thus the screening of diabetic population is of great public health meaning.
Nowadays, researchers indeed have developed some auxiliary diagnosis system using artificial neural network[28], support vector machine and Gaussian mixture model[29]. But all these methods underutilize the grey level and texture information contained in eye fundus images. In 2012, Acharya et al[30] did apply grey level co-occurrence matrix to extract texture information from fundus photograph and they indeed achieve a good sensitivity of 0.989. However, the specificity is only 0.895. What is more, Acharya et al[30] tried to detect diabetic features namely vessels, exudates, and hemorrhages in fundus images first, and then use these features to diagnose DR. During this process, there must be some information loss, and this process ask ophthalmologists mark each DR feature on the images before training model, which adds too much work to them. It is not good for clinical practice, and I think that explains why only 450 fundus images were used in this study.
In our study, because the textures of DR patients are different from diabetes patients, we use textures of each fundus photograph directly to train our model. There is no information loss and could save lots of ophthalmologists’ time. Usually, when we use artificial intelligence to diagnosis DR, a model training process is necessary, and the training of the diagnose model requires the ophthalmologists to mark the symptoms of DR on the eye fundus images firstly. Then the model treat these marks as a reference and learn by itself. For a new eye fundus image, artificial intelligence will try to judge whether there is a similar symptom what it has learned, and then make diagnosis of DR. The image marking procedure takes a lot of ophthalmologists’ time. However, GLCM takes use of disease-related information in another way. GLCM doesn’t require any so-call gold-standard mark because GLCM extracts and transforms images’ textures into quantitative values. None diagnosis on DR-related symptoms needed to be done. The ophthalmologists only need to spend time on diagnosis of DR. Using GLCM is a much more cost-effective way to diagnose DR. The most important thing is that we got both high sensitivity and specificity.
The limitation of our study is that eye fundus images of low quality were excluded. A small part of patients, like cataract patients, were unable to be screening by GLCM and Bayesian model. Actually, low quality of image is a common barrier for automatic diagnosis. Image enhancement techniques might be a promising method to figure out this problem[31].
Textures extracted by grey level co-occurrence can be useful information for DR diagnosis, and a trained Bayesian model based on these textures could be an effective tool for DR screening among diabetic patients.
ACKNOWLEDGEMENTS
Authors’ contributions: Cao K did the analysis; Xu J collected the data; Zhao WQ wrote the manuscript.
Foundation: Supported by the Priming Scientific Research Foundation for the Junior Researcher in Beijing Tongren Hospital, Capital Medical University.
Conflicts of Interest: Cao K, None; Xu J, None; Zhao WQ, None.
REFERENCES