# A Hybrid Unsupervised Approach for Retinal Vessel Segmentation

Authors: Khan Bahadar Khan; Muhammad Shahbaz Siddique; Muhammad Ahmad; Manuel Mazzara
Journal: BioMed Research International
 (2020)

---

## Abstract

Retinal vessel segmentation (RVS) is a significant source of useful information for monitoring, identification, initial medication, and surgical development of ophthalmic disorders. Most common disorders, i.e., stroke, diabetic retinopathy (DR), and cardiac diseases, often change the normal structure of the retinal vascular network. A lot of research has been committed to building an automatic RVS system. But, it is still an open issue. In this article, a framework is recommended for RVS with fast execution and competing outcomes. An initial binary image is obtained by the application of the MISODATA on the preprocessed image. For vessel structure enhancement, B-COSFIRE filters are utilized along with thresholding to obtain another binary image. These two binary images are combined by logical AND-type operation. Then, it is fused with the enhanced image of B-COSFIRE filters followed by thresholding to obtain the vessel location map (VLM). The methodology is verified on four different datasets: DRIVE, STARE, HRF, and CHASE_DB1, which are publicly accessible for benchmarking and validation. The obtained results are compared with the existing competing methods.


---

## Body

## 1. Introduction

The most essential sensory system for gathering information, navigation, and learning is the human visual system [1]. The retina is the sensitive part of the eye that contains fovea, light receptors, Optic disk, and macula. The retina is a layered tissue, coating the interior of the eye, which is an initial sensor of the communication system and gives a sense of sight. Moreover, it allows understanding the colors, dimensions, and shape of objects by processing the amount of light it reflects or emits. Retina image of an eye is captured with a fundus camera [2]. RGB photographs of the fundus are the protrusion of the internal surface of an eye. Imaging of the retina has emerged swiftly and now one of the most common practices in healthcare and for screening the patients suffering from ophthalmologic or systemic diseases. For identify ing numerous ophthalmologic diseases, the ophthalmologist uses vessel condition as an indicator which is a vital component in retinal fundus images.Critical diagnostic to eye diseases in human retinal images can be indicated by its shape analysis, its appearance, blood vessels, morphological features, and tortuosity [3]. Structure of RVS is also used for screening of brain and heart stock diseases [4, 5]. Retinal vessel structures play a significant role among other structures in fundus images. RVS is the elementary phase utilized for the examination of retina images [6]. Vascular-related diseases are diagnosed with the help of vessel delineation which is an important component of medical image processing. Additionally, ongoing research in the area of deep learning suggested multiple approaches with emphasis on the separation and the delineation of the vasculature.The inadequate number of images and having low-contrast in publicly available retina datasets is challenging for deep learning-based research. A dataset having a large number of retina images captured with a different imaging system and under diverse environmental conditions is required to train the supervised network. Deep learning-based methods will aid to control blindness, timely and precise identification of diseases for successful remedy, and thus vividly increase the life quality of patients with eye ailments [7]. RVS is a very difficult task due to many reasons:
(1)
The structure and formation of retinal vessels are very complex and there is a prominent dissimilarity in various local parts regarding the shape, size, and intensity in vessels.(2)
Some structures have the same intensity and shape as vessels, e.g., hemorrhage. Moreover, there are also thin microvessels, whose width is normally between ranges from one to a few pixels and which can be easily mixed with the background. There are irregular illumination in the images and having low-varying contrast [7, 8]. Typically, noise in fundus images is added by the image-capturing procedure such as artifact on the lens or movement of the patient [9]. It is hard to differentiate vessels from other structures that are similar or noises in the retina image. In other words, thicker vessels are more prominent in comparison to the thinner ones as shown in Figure 1(3)
Different manual graders have different segmentation results. Manual RVS is also a very hard and tedious task. Over the recent two decades, automatic RVS has caught noteworthy attention and numerous such techniques are developed but they have performance degradation with the change of datasets. Some of the techniques are not fully automatic while others are incapable to handle pathological images. Some of these methods are evaluated on the datasets having a limited number of images while others have problems of oversegmentation or undersegmentation with abnormal images [10]. Hence, the dilemma of perfect RVS is still not answered.Figure 1
Graphical representation of thinner and thicker vessels.Automated RVS techniques provide incredible support to the ophthalmologist in terms of identification and medication of numerous ophthalmological abnormalities. In this article, an automatic unsupervised approach is developed for RVS that consists a combination of the preprocessing steps, segmentation, vessel structure-based enhancement, and postprocessing steps. The preprocessing steps aim at exterminating noise and improving the contrast of the fundus image. Segmentation is performed by using the Modified Iterative Self Organizing Data Analysis Technique (MISODATA) to acquire a binary image that is fused with the segmented image of the Combination Of Shifted Filter Responses (B-COSFIRE). Then, the fused image is multiplied with the enhanced image of the B-COSFIRE to obtain the initial vessel location map (VLM). Lastly, the VLM and the fused image are combined by logical OR-type operators to obtain final results. In a nutshell, the main contributions of this research are the following:(1)
A mask image is not provided with all retina datasets. Automatic masking creation is proposed for each image to extract ROI which suppresses the false positive rate (FPR).(2)
The proposed efficient denoising process (preprocessing steps) improves the selection of a suitable threshold.(3)
The basic ISODATA algorithm only one-time process the retina image locally and then globally, which sometimes makes it unable to find an optimal threshold. The modified ISODATA technique is introduced to find the global threshold of the entire image which is compared and equated with the individual local threshold of each segment in order to find the optimal threshold for more precise detection of vessels.(4)
The vessel location map (VLM) is a new scheme to achieve better performance. In this scheme, the background noise eradication and vessel enhancement are accomplished independently.(5)
A distinctive postprocessing steps (AND-type and OR-type operations) to reject misclassified foreground pixels.


## 2. Related Works

Numerous methodologies for RVS have been developed in literature [4, 10]. These methodologies are arranged into two sets: supervised and unsupervised procedures. Supervised techniques utilizing a trained classifier for pixel classification into the foreground or background. Supervised techniques utilized various classifiers, for instance, adaptive boosting (AdaBoost), support vector machines (SVM), neural networks (NN), Gaussian mixture models (GMM) and k-nearest neighbors (k-NN).A RVS method utilizing a supervisedk-NN classifier for isolation of foreground and background pixels was recommended by Niemeijer et al. [11], with a feature vector (FV) formation based on a multiscale (MS) Gaussian filter. Staal et al. [12] projected an equivalent RVS methodology using an FV generated based on a ridge detector. A feed-forward NN built classifier was applied by Marin et al. [13], using 7-D FV generated based on moment-invariant.An SVM-based approach was presented by Ricci et al. [14], utilizing FV constructed through a rotation-invariant linear operator and pixel intensity. An AdaBoost classifier was suggested by Lupascu et al. [15], utilizing a 41−D feature set. An ensemble-based RVS system applying a simple linear iterative clustering (SLIC) algorithm was presented by Wang et al. [16]. A GMM classifier-based scheme was recommended by Roychowdhury et al. [17], utilizing 8−D FV extracted from the pixel neighborhood on first and second-order gradient images.Zhu et al. [18] offered an extreme learning machine(ELM) based RVS scheme utilizing a 39−D FV generated by morphological and local attributes combined with attributes extracted from phase congruency, Hessian, and divergence of vector fields (DVF). Tang et al. [19] recommended an SVM-based RVS scheme utilizing an FV created based on MS vessel filtering and the Gabor wavelet features. A random forest classifier-based RVS system was proposed by Aslani et al. [20], utilizing a 17−D FV created based on MS and the multiorientation Gabor filter responses and intensity feature combined with feature extracted from vesselness measure and B-COSFIRE filter.A directionally sensitive vessel enhancement-based scheme combined with NN derived from the U-Net model was presented in [21]. Thangaraj et al. [22] constructed a 13−D FV from the Gabor filter responses, Frangi’s vesselness measure (1D), local binary pattern feature (1D), Hu moment invariants (7D), and grey-level cooccurrence matrix features (3D) for RVS utilizing NN-based approach. Memari et al. [23] recommended an arrangement of various enhancement techniques with the AdaBoost classifier to segregate foreground and background pixels.A three-stage (thick vessel extraction, thin vessel extraction, and vessel fusion-based) deep learning approach were proposed in [24]. Guo et al. [25] suggested an MS deeply supervised network with short connections (BTS-DSN) for RVS. Local intensities, local binary patterns, a histogram of gradients, DVF, higher-order local autocorrelations, and morphological transformation features were used for RVS in [26]. Random forests were used for the selection of feature sets which were utilized in combination with the hierarchical classification methodology to extract the vessels.Alternatively, unsupervised systems are categorized based on matched filtering (MF), mathematical morphology (MM), and multiscale-based approaches. In matched filtering approaches, thick and thin vessels are extracted by the selection of large and small filter kernels, respectively. However, the application of large kernels can accurately detect major vessels with the misclassification of thin vessels by increasing its width. Similarly, smaller kernels can accurately extract thin vessels along with the extraction of thick vessels in reduced widths. To obtain a complete vascular network, a conventional MF technique can be applied with a large number of diverse filter masks in various directions.Similar methods were employed using MF [27–32], combined filters [33], COSFIRE filters [3, 5, 34–36], Gaussian filters [37], wavelet filters [38], and Frangi’s filter [39]. The MM-based approaches are utilized for isolating retinal image segments such as optic disk, macula, fovea, and vasculature. Morphological operators utilized the application of structuring elements (SE) to images for extraction and representation of region contours. A morphological operation for detecting particular structures has the benefit of speed and noise elimination. But they are unable to achieve the known vessel cross-sectional shape. Moreover, there is an issue to extract extremely tortuous vessels in case of superimposing large SE. Morphological operations were utilized for both enhancement and RVS [2, 40–44]. On the other hand, retinal blood vessels of variable thickness at various scales were obtained by multiscale approaches [45–50].


## 3. Proposed Model

The complete structure of the proposed RVS framework is introduced in this section. The information and description of every stage are also presented in subsections.


### 3.1. Overview

The proposed framework consists of two major blocks to obtain a final binary image: (1) retina image denoising and segmentation and (2) vessel structure-based enhancement and segmentation. The key objective of this framework is to extract vasculature excellently along with the elimination of noise and supplementary disease falsifications. The complete structure of the proposed framework is labeled in Figure2. In which Block-I consists of the selection of suitable retina channel, contrast enhancement, noise filtering, region of interest (ROI) extraction, thresholding, and post processing steps. Block-II includes the application of B-COSFIRE filter, logical operations, and postprocessing steps. The initial binary vessel map of Block-I is fused with the B-COSFIRE filter segmented image in Block-II. Then, it is multiplied with the B-COSFIRE filter-enhanced image which is further thresholded. This output image is combined with the initial postprocessed image by the logical OR-type operation to obtain the final binary.Figure 2
Sequential sketch of the proposed framework.


### 3.2. Block-I: Retina Image Denoising and Segmentation

In the first block, the retina image is passed through selected techniques to extract the initial denoised vessel map. The green band of the RGB retina image is extracted and nominated for subsequent operation due to its noticeable contrast difference between the vessel and other retina structures. The RGB retina images generally have contrast variations, low resolution and noise. To avoid such variations and produce more appropriate image for further processing, the vessel light reflex elimination and background uniformity operations are performed. Retinal vessel structures have poor reflectance when equated to other retinal planes. Some vessels contain a bright stripe (light reflex) which runs down the central length of the vessel. To overcome this problem, a disc-shape opening operator with a 3-pixel width SE is used on the green plane. A minimal value of disc width is selected to avoid the absorption of close vessels. The background uniformity and smoothness of random salt-and-pepper noise are obtained by the application of a3×3 mean filter. Additional noise flattening is achieved with the application of a Gaussian kernel of size 9×9, mean=0, and variance 1.8.CLAHE [51, 52] is applied on the preprocessed green channel to make vessel structures prominent. The CLAHE operation divides the input image into blocks (size 8×8 in our case) with the constraint of contrast improvement which is set to 0.01. The clip limit suppresses the noise level and escalates the contrast. The effect of the CLAHE process (Iclahe) along with the green plane is displayed in Figure 3. Histogram-based graphical demonstration of the contrast improvement operations is displayed in Figure 4. An averaging filter of size 49×49 is applied for smoothness and elimination of anatomical regions (e.g., optic disk, macula, and fovea). Iavg symbolizes the output image of the averaging filter. The difference image (Id) is computed for all pixels as follows.
(1)Idm,n=Iavgm,n−Iclahem,n.Figure 3
(a) DRIVE dataset color image. (b) Green plane. (c) CLAHE output.
(a)(b)(c)Figure 4
Histogram pictorial effects. (a) Green plane. (b) CLAHE.
(a)(b)The extra regions of the retinal image are cropped by the utilization of the masking method to extract ROI which reduced the computational complexity. An automatic mask is created from the red band of the retinal image. The reason behind using the red channel for mask construction is that it has a good vessel-background dissimilarity. The automatic mask is created for all datasets because the mask image is not available in some datasets.Id is thresholded by the MISODATA algorithm. The subsequent procedure is used to compute the threshold level, and the application of MISODATA is shown in Algorithm 1.Algorithm 1: MISODATA algorithm.
1functionlevel=isodataId;2 Step 1: compute the mean intensity of image from histogram;3 setT=meanI;4counts,N=imhistId;5i⟵ iteration, let i=1;6mu=cumsumcounts;7Ti=sumN∗counts/muend;8Ti=roundTi;9 M⟵ mean intensity of Id utilizing histogram;10 Step 2:fori⟵1do11   Compute MAT⟵ mean above threshold using T from step 1;12   Compute MBT⟵ mean below threshold using T from step 1;13i=i+1;14mu2=cumsumcounts1:Ti;15MBT=sumN1:Ti∗counts1:Ti/mu2end;16mu3=cumsumcountsTi:end;17MAT=sumNTi:end∗countsTi:end/mu3end;18T=MAT+MBT/2;19ifTi≠Ti−1then20     go toStep 2;21else22level=T−1/Nend−1;23end24 Step 3: divide the image into square local regions as follows:25M,N=sizeId;26   0I_bnry=zerosM,N;27W=floorN/3,H=floorM/2; fori=1:H:M−H+1do28forj=1:W:N−W+1do29Temp=Idi:i+H−1,j:j+W−1;30lval=isodataTemp;31iflval>levelthen32set=level−0.015;33else iflval<=level−level/2then34set=level−0.04;35else36set=level−0.03;37end38I_bnryi:i+H−1,j:j+W−1=im2bwTemp,set−0.00939end40end41endThe isolated pixels with an area less than25 pixels in the image (Is1) are trimmed and fused with the B-COSFIRE filter segmented image of Block-II by AND-type operation. The physical stats (eccentricity and area) are utilized for the rejection of nonvessel structures. The vessel structures have a higher area and eccentricity as their pixels are linked and having an elongated structure. Figure 5 indicates the graphical results of the Iavg, Id, and Is1.Figure 5
(a) Average filter image. (b) Difference image. (c) Segmented initial vessel map.
(a)(b)(c)


### 3.3. Block-II: Vessel Structure-Based Enhancement and Segmentation

In Block-II, the masked image of the Block-I is used as an input for vessel structure-based enhancement and RVS. B-COSFIRE filter [5] is applied for contrast improvement of vessel structures that will enhance noise also along with the enhancement of vessel structures if the image is not preprocessed. Therefore, the masked image is used for further processing. B-COSFIRE filter produced two results: binary segmented image (IsC) and vessel structure-based enhanced image (IEC). The outputs of B-COSFIRE filter are displayed in Figure 6. The AND-type operation is used to combine Is1 with IsC that produced output image denoted by IAnd. The effect of AND-type operation is shown in Figure 7, which demonstrates that if an alternative operator like OR-type is utilized, it will introduce noise and misclassification. The advantage of using an AND-type operator is exposed in Figure 8 by displaying the visual results with and without using the AND-type operator. The IAND is postprocessed (Ip1) and multiplied with IEC which is further thresholded to obtain a segmented image (Is2). Pixel-by-pixel multiplication aims at ensuring the detection of vessels at their correct position. The logical OR-type operation is used to produce the final result by coupling of Ip1 and Is2. The visual effects of the OR-type operator are presented in Figure 9.Figure 6
(a) B-COSFIRE binary segmented image (b) Enhanced image based on symmetric-asymmetric filter responses. (c) AND-type operation output.
(a)(b)(c)Figure 7
Inspecting the results of AND-type operator. (a) Manual image. (b) AND-type. (c) OR-type.
(a)(b)(c)Figure 8
Inspecting the results of AND-type operator. (a) Manual image. (b) AND-type. (c) without AND-type.
(a)(b)(c)Figure 9
Inspecting the results of OR-type operator. (a) Manual image. (b) OR-type. (c) without OR-type.
(a)(b)(c)The B-COSFIRE filter application includes convolution with difference of Gaussian (DoG) filters, its blurring effects, shifting the blurred responses, and an approximate point-wise weighted geometric mean (GM). A DoG functionDoGσx,y is given by [5]
(2)DoGσx,y=12πσ2exp−x2+y22σ2−12π0.5σ2exp−x2+y220.5σ2,where σ is the standard deviation (SD) of the Gaussian function (GF) that decides the range of the boundary. 0.5σ is manually set SD value of the internal GF, and x,y symbolizes the pixel position of the image. Response of DoG filter Cσx,y with kernal function of DoGσx−x′,y−y′ has been estimated by convolution, where (x′,y′) denotes pixels intensity distribution.
(3)Cσx,y=I∗DoGσ+,where .+ represents the half-wave rectification process to reject negative values.In the B-COSFIRE filter, three factors (σi,ρi,∅i) are used to represent each point i, where σi=SD of the DoG filter, while ρi and ∅i denote the polar coordinates. This set of parameters is indicated by S=σi,ρi,∅i∣i=1,⋯,n, where n represents the figure of measured DoG responses. The blurring process indicates the calculation of the extreme limit of the weighted thresholded responses of a DoG filter. The blurring operation is shown as follows.
(4)σ′=σ0′+αρi,where σ0′ and α are constants. Each DoG-blurred outcome is moved in the reverse direction to ∅i by a gap ρi, and as a result, they can merge at the support center of the B-COSFIRE filter. Blurred and shifted responses of the DoG filter is indicated by Sσi,ρi,∅ix,y for every tuple σi,ρi,∅i in set S. The ith blurred and shifted response of the DoG filter is defined as
(5)Sσi,ρi,∅ix,y=maxx′,y′cσix−δxi−x′,y−δyi−y′Gσix′,y′,where −3σ′≤x′,y′≤3σ′. The output of the filter is shown as GM of all the blurred and shifted DoG responses.
(6)rSx,y=∏t=1SSσi,ρi,∅ix,yωi∑i=1Sωit,where ωi=exp−ρi2/2σ2 and .t symbolizes the thresholding response at t,0≤t≤1. Equation (6) represents the AND-type operation that is attained by the B-COSFIRE filter only when all DoG filter responses Sσi,ρi,∅i are larger than zero. The overall step-by-step visual results according to the block diagram (Figure 2) are portrayed in Figure 10.Figure 10
Step-by-step visual results of the proposed method using image from the DRIVE dataset. (a) Color input image. (b) Green plane. (c) CLAHE. (d) Low-pass filter. (e) Difference output. (f) MISODATA. (g) Binarized vessel map of B-COSFIRE. (h) AND operatorf,g. (i) B-COSFIRE enhanced image. (j) Final result.
(a)(b)(c)(d)(e)(f)(g)(h)(i)(j)


## 3.1. Overview

The proposed framework consists of two major blocks to obtain a final binary image: (1) retina image denoising and segmentation and (2) vessel structure-based enhancement and segmentation. The key objective of this framework is to extract vasculature excellently along with the elimination of noise and supplementary disease falsifications. The complete structure of the proposed framework is labeled in Figure2. In which Block-I consists of the selection of suitable retina channel, contrast enhancement, noise filtering, region of interest (ROI) extraction, thresholding, and post processing steps. Block-II includes the application of B-COSFIRE filter, logical operations, and postprocessing steps. The initial binary vessel map of Block-I is fused with the B-COSFIRE filter segmented image in Block-II. Then, it is multiplied with the B-COSFIRE filter-enhanced image which is further thresholded. This output image is combined with the initial postprocessed image by the logical OR-type operation to obtain the final binary.Figure 2
Sequential sketch of the proposed framework.


## 3.2. Block-I: Retina Image Denoising and Segmentation

In the first block, the retina image is passed through selected techniques to extract the initial denoised vessel map. The green band of the RGB retina image is extracted and nominated for subsequent operation due to its noticeable contrast difference between the vessel and other retina structures. The RGB retina images generally have contrast variations, low resolution and noise. To avoid such variations and produce more appropriate image for further processing, the vessel light reflex elimination and background uniformity operations are performed. Retinal vessel structures have poor reflectance when equated to other retinal planes. Some vessels contain a bright stripe (light reflex) which runs down the central length of the vessel. To overcome this problem, a disc-shape opening operator with a 3-pixel width SE is used on the green plane. A minimal value of disc width is selected to avoid the absorption of close vessels. The background uniformity and smoothness of random salt-and-pepper noise are obtained by the application of a3×3 mean filter. Additional noise flattening is achieved with the application of a Gaussian kernel of size 9×9, mean=0, and variance 1.8.CLAHE [51, 52] is applied on the preprocessed green channel to make vessel structures prominent. The CLAHE operation divides the input image into blocks (size 8×8 in our case) with the constraint of contrast improvement which is set to 0.01. The clip limit suppresses the noise level and escalates the contrast. The effect of the CLAHE process (Iclahe) along with the green plane is displayed in Figure 3. Histogram-based graphical demonstration of the contrast improvement operations is displayed in Figure 4. An averaging filter of size 49×49 is applied for smoothness and elimination of anatomical regions (e.g., optic disk, macula, and fovea). Iavg symbolizes the output image of the averaging filter. The difference image (Id) is computed for all pixels as follows.
(1)Idm,n=Iavgm,n−Iclahem,n.Figure 3
(a) DRIVE dataset color image. (b) Green plane. (c) CLAHE output.
(a)(b)(c)Figure 4
Histogram pictorial effects. (a) Green plane. (b) CLAHE.
(a)(b)The extra regions of the retinal image are cropped by the utilization of the masking method to extract ROI which reduced the computational complexity. An automatic mask is created from the red band of the retinal image. The reason behind using the red channel for mask construction is that it has a good vessel-background dissimilarity. The automatic mask is created for all datasets because the mask image is not available in some datasets.Id is thresholded by the MISODATA algorithm. The subsequent procedure is used to compute the threshold level, and the application of MISODATA is shown in Algorithm 1.Algorithm 1: MISODATA algorithm.
1functionlevel=isodataId;2 Step 1: compute the mean intensity of image from histogram;3 setT=meanI;4counts,N=imhistId;5i⟵ iteration, let i=1;6mu=cumsumcounts;7Ti=sumN∗counts/muend;8Ti=roundTi;9 M⟵ mean intensity of Id utilizing histogram;10 Step 2:fori⟵1do11   Compute MAT⟵ mean above threshold using T from step 1;12   Compute MBT⟵ mean below threshold using T from step 1;13i=i+1;14mu2=cumsumcounts1:Ti;15MBT=sumN1:Ti∗counts1:Ti/mu2end;16mu3=cumsumcountsTi:end;17MAT=sumNTi:end∗countsTi:end/mu3end;18T=MAT+MBT/2;19ifTi≠Ti−1then20     go toStep 2;21else22level=T−1/Nend−1;23end24 Step 3: divide the image into square local regions as follows:25M,N=sizeId;26   0I_bnry=zerosM,N;27W=floorN/3,H=floorM/2; fori=1:H:M−H+1do28forj=1:W:N−W+1do29Temp=Idi:i+H−1,j:j+W−1;30lval=isodataTemp;31iflval>levelthen32set=level−0.015;33else iflval<=level−level/2then34set=level−0.04;35else36set=level−0.03;37end38I_bnryi:i+H−1,j:j+W−1=im2bwTemp,set−0.00939end40end41endThe isolated pixels with an area less than25 pixels in the image (Is1) are trimmed and fused with the B-COSFIRE filter segmented image of Block-II by AND-type operation. The physical stats (eccentricity and area) are utilized for the rejection of nonvessel structures. The vessel structures have a higher area and eccentricity as their pixels are linked and having an elongated structure. Figure 5 indicates the graphical results of the Iavg, Id, and Is1.Figure 5
(a) Average filter image. (b) Difference image. (c) Segmented initial vessel map.
(a)(b)(c)


## 3.3. Block-II: Vessel Structure-Based Enhancement and Segmentation

In Block-II, the masked image of the Block-I is used as an input for vessel structure-based enhancement and RVS. B-COSFIRE filter [5] is applied for contrast improvement of vessel structures that will enhance noise also along with the enhancement of vessel structures if the image is not preprocessed. Therefore, the masked image is used for further processing. B-COSFIRE filter produced two results: binary segmented image (IsC) and vessel structure-based enhanced image (IEC). The outputs of B-COSFIRE filter are displayed in Figure 6. The AND-type operation is used to combine Is1 with IsC that produced output image denoted by IAnd. The effect of AND-type operation is shown in Figure 7, which demonstrates that if an alternative operator like OR-type is utilized, it will introduce noise and misclassification. The advantage of using an AND-type operator is exposed in Figure 8 by displaying the visual results with and without using the AND-type operator. The IAND is postprocessed (Ip1) and multiplied with IEC which is further thresholded to obtain a segmented image (Is2). Pixel-by-pixel multiplication aims at ensuring the detection of vessels at their correct position. The logical OR-type operation is used to produce the final result by coupling of Ip1 and Is2. The visual effects of the OR-type operator are presented in Figure 9.Figure 6
(a) B-COSFIRE binary segmented image (b) Enhanced image based on symmetric-asymmetric filter responses. (c) AND-type operation output.
(a)(b)(c)Figure 7
Inspecting the results of AND-type operator. (a) Manual image. (b) AND-type. (c) OR-type.
(a)(b)(c)Figure 8
Inspecting the results of AND-type operator. (a) Manual image. (b) AND-type. (c) without AND-type.
(a)(b)(c)Figure 9
Inspecting the results of OR-type operator. (a) Manual image. (b) OR-type. (c) without OR-type.
(a)(b)(c)The B-COSFIRE filter application includes convolution with difference of Gaussian (DoG) filters, its blurring effects, shifting the blurred responses, and an approximate point-wise weighted geometric mean (GM). A DoG functionDoGσx,y is given by [5]
(2)DoGσx,y=12πσ2exp−x2+y22σ2−12π0.5σ2exp−x2+y220.5σ2,where σ is the standard deviation (SD) of the Gaussian function (GF) that decides the range of the boundary. 0.5σ is manually set SD value of the internal GF, and x,y symbolizes the pixel position of the image. Response of DoG filter Cσx,y with kernal function of DoGσx−x′,y−y′ has been estimated by convolution, where (x′,y′) denotes pixels intensity distribution.
(3)Cσx,y=I∗DoGσ+,where .+ represents the half-wave rectification process to reject negative values.In the B-COSFIRE filter, three factors (σi,ρi,∅i) are used to represent each point i, where σi=SD of the DoG filter, while ρi and ∅i denote the polar coordinates. This set of parameters is indicated by S=σi,ρi,∅i∣i=1,⋯,n, where n represents the figure of measured DoG responses. The blurring process indicates the calculation of the extreme limit of the weighted thresholded responses of a DoG filter. The blurring operation is shown as follows.
(4)σ′=σ0′+αρi,where σ0′ and α are constants. Each DoG-blurred outcome is moved in the reverse direction to ∅i by a gap ρi, and as a result, they can merge at the support center of the B-COSFIRE filter. Blurred and shifted responses of the DoG filter is indicated by Sσi,ρi,∅ix,y for every tuple σi,ρi,∅i in set S. The ith blurred and shifted response of the DoG filter is defined as
(5)Sσi,ρi,∅ix,y=maxx′,y′cσix−δxi−x′,y−δyi−y′Gσix′,y′,where −3σ′≤x′,y′≤3σ′. The output of the filter is shown as GM of all the blurred and shifted DoG responses.
(6)rSx,y=∏t=1SSσi,ρi,∅ix,yωi∑i=1Sωit,where ωi=exp−ρi2/2σ2 and .t symbolizes the thresholding response at t,0≤t≤1. Equation (6) represents the AND-type operation that is attained by the B-COSFIRE filter only when all DoG filter responses Sσi,ρi,∅i are larger than zero. The overall step-by-step visual results according to the block diagram (Figure 2) are portrayed in Figure 10.Figure 10
Step-by-step visual results of the proposed method using image from the DRIVE dataset. (a) Color input image. (b) Green plane. (c) CLAHE. (d) Low-pass filter. (e) Difference output. (f) MISODATA. (g) Binarized vessel map of B-COSFIRE. (h) AND operatorf,g. (i) B-COSFIRE enhanced image. (j) Final result.
(a)(b)(c)(d)(e)(f)(g)(h)(i)(j)


## 4. Experimental Outcomes and Deliberation

This section will provide the information about datasets, performance metrics, analysis of experimental results, and time complexity of the proposed method.


### 4.1. Datasets

The proposed system obtained remarkable results on the freely online available datasets: DRIVE [11, 12], STARE [53], HRF [54], and CHASE_DB1 [55]. The magnificence of the framework is justified in terms of assessment with state-of-the-art systems. The datasets used for endorsement of the suggested framework are encapsulated in Table 1. The manually labeled results in all datasets are utilized as a gold standard for performance assessment of the proposed framework.Table 1
Datasets comparison.
DatasetImage classificationImage sizeFormatDRIVETotal 40 images565×584JPEG20 test, 20 training7 abnormal, 33 normalSTARETotal 20 images700×605PPM10 normal10 abnormalHRFTotal 45 images3504×2336JPEG15 normal, 15 DR15 glaucomatousCHASE_DBITotal 28 images1280×960JPEG14 left eye14 right eye


### 4.2. Performance Judgment Parameters

The quantitative results are obtained by equating the proposed segmentation’s with the manual segmentation available on each dataset. There are numerous performance standards mentioned in the literature. The performance metrics used for evaluation of the proposed framework are visible in Table2. Six performance standards (Acc, Sn, Sp, AUC, MCC, and CAL) are selected for the justification of the proposed methodology. The Acc metric tells about the overall valuation of the proposed method. Sn is a measure of the quantity of correctly classified vessel pixels, while Sp is an assessment of the competency of differentiating nonvessel pixels. The AUC is the ratio of Sn and Sp. The MCC [5, 56] is a more appropriate indicator of the accuracy of binary categorization in the case of unbalanced structures. For a comprehensive judgment of the superiority of segmentation, the CAL metric [57, 58] is computed. This metric provides justification based on the properties (connectivity-area-length) of the segmented structures beyond the correctly classified image pixels.Table 2
Performance judgment criteria of the proposed model.
ParameterFormulationSensitivity (Sn)TPTP+FNSpecificity (Sp)1-FPR orTN/TN+FPAccuracy (Acc)TN+TPTN+FP+TP+FNArea under ROC curve (AUC)Sn+Sp/2Matthews correlation coefficient (MCC)TP/N−S×PP×S×1−S×1−PConnectivity-area-length (CAL)fC,A,L=C×A×LIn Table2, N=TN+TP+FN+FP,S=TP+FN/N and P=TP+FP/N [58]. The terms TP, TN, FP, and FN denote the true positive (exactly matched vessel pixels), true negative (exactly matched nonvessel pixels), false positive (invalidly predicated vessel pixels), and false negative (invalidly predicated nonvessel pixels), correspondingly.LetIS be the extracted final binary image and IG the corresponding manual segmented image. The considered metric evaluates the following [57, 58]:
(i)
Connectivity (C): it calculates the fragmentation grade ofIS with respect to the manual segmentation IG and penalizes fragmented segmentation. It is computed as(7)CIS,IG=1−min1,#CIG−#CIS#IG,where #C· sums the linked segments while #· measures the number of vessel pixels in the considered binary image
(ii)
Area (A): it estimates the intersecting area betweenIS and IG, based on the Jaccard coefficient. Let δε· be a morphological dilation that utilizes a disc structuring element (SE) with a radius of ε pixels. The magnitude A is calculated as follows:(8)AIS,IG=#δεIS∩IS∩IG∪IS∩δεIG#IG∪IS.The value of ε controls the tolerance to lines of various sizes. We set ε=2(iii)
Length (L): it determines the equivalent degree betweenIS and IG by computing the length of the two line networks:(9)LIS,IG=#φIS∩δβIG∪δβIS∩φIG#φIS∪φIG,where φ· is a skeletonization process and δβ· is a morphological dilation with a disc SE of β pixel radius. The value of β controls the tolerance to dissimilarity of the line tracing output. We set β=2. The final assessment parameter, named CAL, is demarcated as fC,A,L=C·A·AL.


### 4.3. Experimental Results and Inspection

The success of the proposed framework is established by utilizing four freely obtainable datasets: DRIVE, STARE, HRF, and CHASE_DB1 for testing and evaluation. The average performance parameters results in Table3 are computed by processing 20 test images of the DRIVE and STARE datasets. The performance scores of HRF dataset (15 normal images, 15 DR, and 15 glaucomatous) and CHASE_DB1 are presented in Tables 4 and 5 and Table 6, respectively. The best and worst results within Tables 3–6 are highlighted in italic font. The best and worst image results from each dataset are selected based on their accuracy’s scores. Their pictorial results are shown in Figures 11–14.Table 3
Statistical scores achieved on the DRIVE and STARE datasets.
ImgDRIVESTAREAccSnSepAUCMCCCALAccSnSepAUCMCCCAL10.9610.8100.9760.8930.7640.7640.9930.4940.9970.7450.4680.55020.9660.7710.9880.8830.8020.7830.9960.8320.9770.9150.5940.41830.9480.8120.9610.8870.7120.6600.9880.8380.9880.9130.4760.53940.9670.7570.9870.8720.7860.7010.9950.4190.9990.7090.3830.25450.9670.7920.9830.8880.7800.7670.9890.9090.9890.9490.6250.57560.9570.7550.9780.8660.7410.6530.9950.9290.9960.9620.7350.66370.9590.8210.9690.8950.7210.6570.9960.7850.9970.8910.7090.68480.9470.8430.9540.8990.6600.6210.9970.7920.9980.8950.7140.67890.9550.7710.9710.8710.7090.6320.9980.8280.9990.9140.8430.655100.9620.8190.9730.8960.7360.7310.9930.9410.9930.9670.6950.576110.9580.7950.9730.8840.7370.6810.9990.9710.9990.9850.8920.617120.9550.8560.9640.9100.7390.7500.9990.9540.9990.9770.9610.677130.9550.7030.9830.8430.7360.6720.9980.8420.9990.9210.8190.551140.9520.8690.9590.9140.7170.6890.9970.8670.9990.9330.8480.613150.9470.7830.9600.8720.6660.5470.9980.7620.9990.8810.7730.685160.9620.7980.9770.8880.7620.7430.9950.6590.9980.8290.7140.489170.9540.8360.9630.8990.7160.6980.9990.9750.9990.9870.9610.612180.9610.8130.9760.8950.7710.7300.9990.7300.9990.8650.8290.531190.9610.7540.9840.8690.7750.7090.9990.7640.9990.8820.5380.550200.9570.7730.9760.8750.7460.7220.9960.5470.9980.7730.5620.397Avg0.9580.7970.9730.8850.7390.6960.9960.7920.9970.8950.7070.566Table 4
Average efficiency scores on the HRF dataset (normal and diabetic images).
ImgNormalDiabeticAccSnSepAUCMCCCALAccSnSepAUCMCCCAL10.9580.7010.9880.8440.7000.5740.9630.7680.9770.8450.6460.74220.9650.8060.9830.8940.7870.6710.9580.7550.9740.8430.6710.67230.9540.7130.9830.8480.7410.6320.9490.7750.9610.8600.6310.56340.9570.8100.9750.8820.7610.7190.9510.7580.9660.8360.6330.53150.9680.7650.9880.8760.7950.7540.9640.7250.9830.8230.7090.60960.9600.8320.9780.8880.7740.7590.9500.5640.9880.7280.6190.39470.9680.8170.9880.8780.7930.7410.9570.7370.9800.8340.7020.68880.9640.8060.9850.8790.7880.7290.9520.7760.9710.8410.6800.68890.9560.8080.9680.8880.6810.6320.9450.7580.9580.8580.6340.551100.9590.7650.9790.8630.7300.5610.9510.6640.9810.8090.6840.598110.9700.8390.9900.8860.8060.7530.9530.7000.9830.8080.7130.611120.9670.8520.9860.8950.8160.7880.9550.7010.9760.8130.6660.579130.9660.8370.9820.8940.7780.7110.9600.7530.9810.8290.7040.645140.9660.8350.9810.8900.7610.6990.9440.7260.9640.8290.6450.622150.9700.8270.9820.9050.7080.7680.9530.7300.9710.8340.6480.547Avg0.9630.8010.9820.8810.7610.6990.9540.7260.9740.8260.6660.603Table 5
Average efficiency scores on the HRF dataset (Glaucomatous images).
ImgGlaucomatousAccSnSepAUCMCCCAL10.9520.8130.9620.8870.6880.74320.9540.8050.9660.8860.7040.67930.9620.8300.9700.9000.6880.61040.9620.7880.9740.8810.7110.69550.9630.8330.9720.9020.7210.66360.9620.8110.9730.8920.7080.71370.9690.7730.9770.8750.7210.60980.9570.8350.9640.8990.7030.72090.9620.7740.9750.8750.7100.646100.9630.7800.9760.8780.7180.665110.9620.7490.9760.8630.7170.641120.9580.7990.9670.8830.7170.688130.9580.7530.9740.8630.6920.664140.9510.7640.9650.8650.6630.646150.9560.7480.9730.8600.6910.620Avg0.9590.7900.9790.8810.7030.667Table 6
Average performance achieved on the CHASE_DB1 dataset.
ImagesAccSnSepAUCMCCCAL01L0.9950.7530.9960.8740.5820.55201R0.9970.7050.9990.8520.7070.48502L0.9950.6470.9970.8220.6230.51602R0.9950.7440.9960.8700.6200.48803L0.9960.7830.9970.8900.6460.57103R0.9980.8480.9990.9230.7930.58204L0.9960.5880.9980.7930.6480.61604R0.9970.6510.9980.8250.6920.53505R0.9970.7580.9980.8780.7460.63406L0.9980.6200.9990.8090.6940.54806R0.9980.7050.9990.8520.7260.50307L0.9970.7220.9980.8600.6780.57407R0.9970.6300.9990.8150.6750.55308L0.9970.8500.9980.9240.6500.50708R0.9970.8670.9970.9320.6410.56609L0.9910.8160.9910.9030.3150.50909R0.9940.9150.9940.9550.3680.49310L0.9940.6720.9950.8340.4080.48310R0.9960.9000.9960.9480.5900.45811L0.9970.9330.9970.9650.6170.62211R0.9980.8780.9980.9380.6560.57012L0.9970.6990.9980.8480.5660.60912R0.9980.7700.9990.8850.6930.51313L0.9980.8330.9980.9160.6670.50313R0.9980.6500.9980.8240.5480.44514L0.9970.8130.9980.9050.6520.63914R0.9970.6920.9990.8450.6280.554Average0.9970.7570.9970.8770.6290.547Figure 11
Visual effects of the best and worst cases from the DRIVE dataset.Figure 12
Visual effects of the best and worst cases from the STARE dataset.Figure 13
Visual results of the best and worst cases from the HRF dataset.Figure 14
Best and worst cases visual results from the CHASE_DB1 dataset.The framework performs well on both healthy and pathological images of all selected datasets. The statistical results in Tables3–6 validates that the suggested system is robust and has the capability to handle the bright lesions images of the STARE dataset, higher resolution images of the HRF dataset, low resolution images of the DRIVE dataset, and left/right eyes images of the CHASE_DB1 dataset. The anatomical structures are also efficiently omitted to avoid any misclassification.The average statistical results of the proposed framework on all selected datasets are displayed in Table7, which reflects that the highest mean score of Acc 0.997, Sn 0.814, Sp 0.997, and AUC 0.905 is achieved on the CHASE_DB1 dataset. The lowest FPR is also observed using the same dataset. The highest value of MCC 0.761 and CAL 0.699 is recorded on the HRF dataset. The highest value of each parameter is italicized in the respective column of the Table 7.Table 7
Mean results comparison of the proposed method on different datasets.
DatasetsImagesAccSnSpAUCMCCCALDRIVE (1st observer)200.9540.7660.9720.8690.7210.690DRIVE (2nd observer)0.9580.7970.9730.8850.7390.696HRF_Normal150.9630.8010.9820.8810.7610.699HRF_DR0.9540.7260.9740.8260.6660.603HRF_Glaucoma0.9590.7900.9790.8810.7030.667HRF_Average450.9590.7720.9780.8630.7100.656CHASE_DB1 (1st observer)280.9970.7570.970.8770.6290.547CHASE_DB1 (2nd observer)0.9960.8140.9960.9050.5690.547STARE200.9960.7920.9970.8950.7070.566The average performance parameter scores of the proposed framework on the DRIVE and STARE datasets are compared with the existing literature in Table8, while Table 9 shows the result comparison of the HRF and CHASE_DB1 datasets. The Acc, Sn, and Sp results of all techniques in Tables 8 and 9 are acquired from their respective published articles while the AUC result is calculated by using the formula in Table 2.Table 8
Performance assessments of existing techniques on the DRIVE and STARE datasets.
MethodYearDRIVESTAREAccSnSpAUCAccSnSpAUCHuman observer0.9470.7790.9720.8740.9350.8950.9380.917Unsupervised techniquesChauduri [31]19890.877——0.788————Zana and Klein [42]20010.9380.697——————Martinez-Perez [45]20070.9340.7250.9650.8450.9410.7510.9550.853Zhang [27]20100.938———0.948———Bankhead [38]20120.9370.7030.9710.8370.9320.7580.9500.854Fraz [44]20120.9430.7150.9760.8450.9440.7310.9680.850Azzopardi [5]20150.9440.7660.9700.8680.9500.7720.9700.871Oliveira [33]20160.9460.8640.9560.9100.9530.8250.9650.895Khan [40]20160.9610.7460.9800.8630.9460.7580.9630.861Biswal [29]20170.9500.7100.9700.8400.9500.7000.9700.835Khan [32]20170.9440.7540.9640.8590.9480.7520.9560.854Soomro [63]20170.9430.7520.9760.8640.9610.7840.9810.883Badawi [3]20180.9550.7910.9710.8810.9530.8650.9610.913Yue [50]20180.9450.7530.9730.863————Soomro [9]20180.9480.7450.9620.8540.9510.7840.9760.880Soomro [64]20180.9530.7520.9760.8640.9670.7860.9820.884Fan [60]20180.9600.7360.9810.8580.9570.7910.9700.880Khan [65]20190.9510.7700.9650.8680.9510.7520.9810.867Memari [59]20190.9610.7610.9810.8710.9510.7820.9650.873Proposed20200.9580.7970.9730.8850.9960.7920.9980.895Supervised techniquesNiemeijer [11]20040.9420.6900.9700.830————Staal [12]20040.9440.7190.9770.8480.9520.6970.9810.839Ricci [14]20070.959———0.964———Lupascu [15]20100.9590.6730.9870.830————MarÃ­n [13]20110.9450.7070.9800.8440.9530.6940.9820.838Wang [16]20150.9770.8170.9730.8950.9810.8100.9790.894Roychowdhury [17]20150.9520.7250.9830.8540.9510.7720.9730.873Aslani [20]20160.9510.7540.9800.8670.9610.7550.9830.869Zhu [18]20170.9610.7140.9870.851————Thangaraj [22]20170.9610.8010.9750.8880.9430.8340.9540.893Memari [23]20170.9720.8720.9880.9300.9510.8090.9790.894Dharmawan [21]2018—0.8310.9720.902—0.7920.9830.887Yan [24]20180.9540.7630.9820.8730.9640.7740.9860.880Guo [25]20190.9550.7800.9810.8810.9660.8200.9830.902Khowaja [26]20190.9750.8180.9710.8950.9750.8240.9750.899Soomro [8]20190.9590.8020.9740.9480.9610.8010.9690.945Soomro [62]20190.9560.8700.9850.9860.9680.8480.9860.988Fan [61]20190.9660.7960.9820.8890.9740.8160.9870.901Table 9
Performance assessments of existing methods with the proposed model on the HRF and CHASE_DB1 datasets.
TechniqueYearHRFCHASE_DB1AccSnSpAUCAccSnSpAUCUnsupervised techniquesOdstrcilik [54]20130.9490.7740.9670.871————Azzopardi [5]2015————0.9390.7590.9590.859Zhang [66]20160.9570.7980.9740.8860.9460.7630.9680.866Biswal [29]2017————0.9400.7600.9700.865Rodrigues [67]20170.9480.7220.9640.843————Badawi [3]2018————0.9530.8000.9640.882Proposed20200.9600.7320.9790.8550.9960.7270.9970.862Supervised techniquesRoychowdhury [17]2015————0.9530.7200.9820.851Thangaraj [22]2017————0.9470.6290.9730.797Memari [23]2017————0.9480.8190.9590.889Dharmawan [21]2018—0.8130.9770.895————Yan [24]2018————0.9610.7640.9810.873Fan [60]2018————0.9510.6570.9730.815Guo [25]2019————0.9630.7890.9800.885Khowaja [26]2019————0.9520.7560.9760.866Soomro [62]20190.9620.8290.9620.9780.9760.8860.9820.985Fan [61]20190.9760.8240.9870.9050.9710.8020.9850.893In Table8, the obtained results of the framework are compared with 19 unsupervised and 18 supervised existing techniques. The proposed framework achieved the highest Acc result than all unsupervised methods on the DRIVE dataset except Khan et al. [40], Memari et al. [59] which is 0.003%, and Fan et al. [60] which is 0.002% better than ours. The supervised methods Ricci and Perfetti [14], Lupascu et al. [15], Wang et al. [16], Zhu et al. [18], Thangaraj et al. [22], Memari et al. [23], Khowaja et al. [26], and Fan et al. [61] show 0.001%, 0.001%, 0.019%, 0.003%, 0.003%, 0.014%, 0.017%, and 0.008% better results than the proposed method, respectively. But some of these methods are only validated on one dataset, which reflects that they are tuned for a single dataset. Some of these methods produce a very low AUC score, which is a trade-off between Sn and Sp. Moreover, supervised methods are computationally very expensive. In the case of the STARE dataset, the framework produced highest Acc scores than all other methods. Table 9 reflects that there are very few techniques that used both HRF and CHASE_DB1 datasets for validation. The Acc score of the framework is higher than both supervised and unsupervised approaches on the HRF and CHASE_DB1 datasets except Soomro et al. [62] and Fan et al. [61] which is slightly higher than ours on HRF dataset only. Fan et al. [61] showed higher Sp value than all other methods on the HRF dataset. The highest Sp value on CHASE_DB1 dataset is obtained by the proposed method. All the other supervised and unsupervised methods acquired a bit greater or equivalent values of Sn and AUC metric on the HRF and CHASE_DB1 datasets as compared to ours.In Table10, the MCC and CAL values are recorded by the proposed method and other existing supervised and unsupervised methods. The MCC and CAL values of Chauduri et al. [31], Niemeijer et al. [11], Hoover et al. [53], and B-COSFIRE [5] are calculated by utilizing their publicly accessible segmented images. The results of Fraz et al. [68, 69], RUSTICO [58], Yang et al. [70, 71], Vega et al. [72], FC-CRF [73], and UP-CRF [73] are extracted from their published articles.Table 10
Performance assessments of existing techniques on the four datasets.
MethodYearDRIVESTAREHRFCHASE_DB1MCCCALMCCCALMCCCALMCCCALUnsupervised techniquesChauduri [31]19890.4200.208——————Hoover [53]2000——0.6150.534————Fraz [68]20110.733—0.700—————Fraz [69]20130.736—0.691—————B-COSFIRE [5]20150.7190.7210.6980.7090.6860.5770.6560.608RUSTICO [58]20190.7290.7280.6980.7090.6910.5870.6630.620Proposed20200.7390.6960.7070.5660.7100.6560.6290.547Supervised techniquesYang [70]20190.736—0.704—0.712———Yang [71]20180.725—0.662—0.682———FC-CRF [73]20160.7560.7310.7270.6580.6900.5410.7040.622UP-CRF [73]20160.7400.6750.7260.6650.6770.4750.6890.571Vega [72]20150.662—0.640—————Niemeijer [11]20040.7220.659——————The average value of MCC attained by the proposed method is higher than all compared unsupervised approaches on the DRIVE, STARE, and HRF datasets, while it is statistically lower than the supervised methods (i.e., FC-CRF [73] and UP-CRF [73]) on the DRIVE, STARE, and CHASE_DB1 datasets. The CAL value of the proposed method is observed higher than all supervised and unsupervised methods on the HRF dataset, while it is statistically lower than or equivalent to CAL values of other methods on the DRIVE, STARE, and CHASE_DB1 datasets.


#### 4.3.1. Processing Time

The proposed framework processes a single image in a very short time as equated to other approaches in Table11. The time values are computed on the single image taken from the DRIVE and STARE datasets.Table 11
Processing time evaluation of the systems.
MethodTimeHardware particularsRoychowdhury [17]3.11 secIntel Core i3 CPU 2.6 GHz, 2 GB RAMZhu [18]12.160 sec4.0 GHz Intel i7-4790K CPU and 32 GB RAMMemari [23]8.2 minsIntel i5-M480 CPU, 2.67 GHz, 4 GB RAMBiswal [29]3.3 secIntel i3 (4010U CPU) 1.7 GHz, 4 GB RAMBadawi [3]8 secCPU 2.7 GHz, 16 GB RAMYue [50]4.6 secIntel i5-6200U CPU 2.3 GHz, 8 GB RAMKhan [39]6.1 sec5∗Intel Core i3 CPU, 2.53 GHz, 4 GB RAMKhan [40]1.56 secAzzopardi [5]11.83 secVlachos [47]9.3 secBankhead [38]22.45 secProposed5.5 sec


## 4.1. Datasets

The proposed system obtained remarkable results on the freely online available datasets: DRIVE [11, 12], STARE [53], HRF [54], and CHASE_DB1 [55]. The magnificence of the framework is justified in terms of assessment with state-of-the-art systems. The datasets used for endorsement of the suggested framework are encapsulated in Table 1. The manually labeled results in all datasets are utilized as a gold standard for performance assessment of the proposed framework.Table 1
Datasets comparison.
DatasetImage classificationImage sizeFormatDRIVETotal 40 images565×584JPEG20 test, 20 training7 abnormal, 33 normalSTARETotal 20 images700×605PPM10 normal10 abnormalHRFTotal 45 images3504×2336JPEG15 normal, 15 DR15 glaucomatousCHASE_DBITotal 28 images1280×960JPEG14 left eye14 right eye


## 4.2. Performance Judgment Parameters

The quantitative results are obtained by equating the proposed segmentation’s with the manual segmentation available on each dataset. There are numerous performance standards mentioned in the literature. The performance metrics used for evaluation of the proposed framework are visible in Table2. Six performance standards (Acc, Sn, Sp, AUC, MCC, and CAL) are selected for the justification of the proposed methodology. The Acc metric tells about the overall valuation of the proposed method. Sn is a measure of the quantity of correctly classified vessel pixels, while Sp is an assessment of the competency of differentiating nonvessel pixels. The AUC is the ratio of Sn and Sp. The MCC [5, 56] is a more appropriate indicator of the accuracy of binary categorization in the case of unbalanced structures. For a comprehensive judgment of the superiority of segmentation, the CAL metric [57, 58] is computed. This metric provides justification based on the properties (connectivity-area-length) of the segmented structures beyond the correctly classified image pixels.Table 2
Performance judgment criteria of the proposed model.
ParameterFormulationSensitivity (Sn)TPTP+FNSpecificity (Sp)1-FPR orTN/TN+FPAccuracy (Acc)TN+TPTN+FP+TP+FNArea under ROC curve (AUC)Sn+Sp/2Matthews correlation coefficient (MCC)TP/N−S×PP×S×1−S×1−PConnectivity-area-length (CAL)fC,A,L=C×A×LIn Table2, N=TN+TP+FN+FP,S=TP+FN/N and P=TP+FP/N [58]. The terms TP, TN, FP, and FN denote the true positive (exactly matched vessel pixels), true negative (exactly matched nonvessel pixels), false positive (invalidly predicated vessel pixels), and false negative (invalidly predicated nonvessel pixels), correspondingly.LetIS be the extracted final binary image and IG the corresponding manual segmented image. The considered metric evaluates the following [57, 58]:
(i)
Connectivity (C): it calculates the fragmentation grade ofIS with respect to the manual segmentation IG and penalizes fragmented segmentation. It is computed as(7)CIS,IG=1−min1,#CIG−#CIS#IG,where #C· sums the linked segments while #· measures the number of vessel pixels in the considered binary image
(ii)
Area (A): it estimates the intersecting area betweenIS and IG, based on the Jaccard coefficient. Let δε· be a morphological dilation that utilizes a disc structuring element (SE) with a radius of ε pixels. The magnitude A is calculated as follows:(8)AIS,IG=#δεIS∩IS∩IG∪IS∩δεIG#IG∪IS.The value of ε controls the tolerance to lines of various sizes. We set ε=2(iii)
Length (L): it determines the equivalent degree betweenIS and IG by computing the length of the two line networks:(9)LIS,IG=#φIS∩δβIG∪δβIS∩φIG#φIS∪φIG,where φ· is a skeletonization process and δβ· is a morphological dilation with a disc SE of β pixel radius. The value of β controls the tolerance to dissimilarity of the line tracing output. We set β=2. The final assessment parameter, named CAL, is demarcated as fC,A,L=C·A·AL.


## 4.3. Experimental Results and Inspection

The success of the proposed framework is established by utilizing four freely obtainable datasets: DRIVE, STARE, HRF, and CHASE_DB1 for testing and evaluation. The average performance parameters results in Table3 are computed by processing 20 test images of the DRIVE and STARE datasets. The performance scores of HRF dataset (15 normal images, 15 DR, and 15 glaucomatous) and CHASE_DB1 are presented in Tables 4 and 5 and Table 6, respectively. The best and worst results within Tables 3–6 are highlighted in italic font. The best and worst image results from each dataset are selected based on their accuracy’s scores. Their pictorial results are shown in Figures 11–14.Table 3
Statistical scores achieved on the DRIVE and STARE datasets.
ImgDRIVESTAREAccSnSepAUCMCCCALAccSnSepAUCMCCCAL10.9610.8100.9760.8930.7640.7640.9930.4940.9970.7450.4680.55020.9660.7710.9880.8830.8020.7830.9960.8320.9770.9150.5940.41830.9480.8120.9610.8870.7120.6600.9880.8380.9880.9130.4760.53940.9670.7570.9870.8720.7860.7010.9950.4190.9990.7090.3830.25450.9670.7920.9830.8880.7800.7670.9890.9090.9890.9490.6250.57560.9570.7550.9780.8660.7410.6530.9950.9290.9960.9620.7350.66370.9590.8210.9690.8950.7210.6570.9960.7850.9970.8910.7090.68480.9470.8430.9540.8990.6600.6210.9970.7920.9980.8950.7140.67890.9550.7710.9710.8710.7090.6320.9980.8280.9990.9140.8430.655100.9620.8190.9730.8960.7360.7310.9930.9410.9930.9670.6950.576110.9580.7950.9730.8840.7370.6810.9990.9710.9990.9850.8920.617120.9550.8560.9640.9100.7390.7500.9990.9540.9990.9770.9610.677130.9550.7030.9830.8430.7360.6720.9980.8420.9990.9210.8190.551140.9520.8690.9590.9140.7170.6890.9970.8670.9990.9330.8480.613150.9470.7830.9600.8720.6660.5470.9980.7620.9990.8810.7730.685160.9620.7980.9770.8880.7620.7430.9950.6590.9980.8290.7140.489170.9540.8360.9630.8990.7160.6980.9990.9750.9990.9870.9610.612180.9610.8130.9760.8950.7710.7300.9990.7300.9990.8650.8290.531190.9610.7540.9840.8690.7750.7090.9990.7640.9990.8820.5380.550200.9570.7730.9760.8750.7460.7220.9960.5470.9980.7730.5620.397Avg0.9580.7970.9730.8850.7390.6960.9960.7920.9970.8950.7070.566Table 4
Average efficiency scores on the HRF dataset (normal and diabetic images).
ImgNormalDiabeticAccSnSepAUCMCCCALAccSnSepAUCMCCCAL10.9580.7010.9880.8440.7000.5740.9630.7680.9770.8450.6460.74220.9650.8060.9830.8940.7870.6710.9580.7550.9740.8430.6710.67230.9540.7130.9830.8480.7410.6320.9490.7750.9610.8600.6310.56340.9570.8100.9750.8820.7610.7190.9510.7580.9660.8360.6330.53150.9680.7650.9880.8760.7950.7540.9640.7250.9830.8230.7090.60960.9600.8320.9780.8880.7740.7590.9500.5640.9880.7280.6190.39470.9680.8170.9880.8780.7930.7410.9570.7370.9800.8340.7020.68880.9640.8060.9850.8790.7880.7290.9520.7760.9710.8410.6800.68890.9560.8080.9680.8880.6810.6320.9450.7580.9580.8580.6340.551100.9590.7650.9790.8630.7300.5610.9510.6640.9810.8090.6840.598110.9700.8390.9900.8860.8060.7530.9530.7000.9830.8080.7130.611120.9670.8520.9860.8950.8160.7880.9550.7010.9760.8130.6660.579130.9660.8370.9820.8940.7780.7110.9600.7530.9810.8290.7040.645140.9660.8350.9810.8900.7610.6990.9440.7260.9640.8290.6450.622150.9700.8270.9820.9050.7080.7680.9530.7300.9710.8340.6480.547Avg0.9630.8010.9820.8810.7610.6990.9540.7260.9740.8260.6660.603Table 5
Average efficiency scores on the HRF dataset (Glaucomatous images).
ImgGlaucomatousAccSnSepAUCMCCCAL10.9520.8130.9620.8870.6880.74320.9540.8050.9660.8860.7040.67930.9620.8300.9700.9000.6880.61040.9620.7880.9740.8810.7110.69550.9630.8330.9720.9020.7210.66360.9620.8110.9730.8920.7080.71370.9690.7730.9770.8750.7210.60980.9570.8350.9640.8990.7030.72090.9620.7740.9750.8750.7100.646100.9630.7800.9760.8780.7180.665110.9620.7490.9760.8630.7170.641120.9580.7990.9670.8830.7170.688130.9580.7530.9740.8630.6920.664140.9510.7640.9650.8650.6630.646150.9560.7480.9730.8600.6910.620Avg0.9590.7900.9790.8810.7030.667Table 6
Average performance achieved on the CHASE_DB1 dataset.
ImagesAccSnSepAUCMCCCAL01L0.9950.7530.9960.8740.5820.55201R0.9970.7050.9990.8520.7070.48502L0.9950.6470.9970.8220.6230.51602R0.9950.7440.9960.8700.6200.48803L0.9960.7830.9970.8900.6460.57103R0.9980.8480.9990.9230.7930.58204L0.9960.5880.9980.7930.6480.61604R0.9970.6510.9980.8250.6920.53505R0.9970.7580.9980.8780.7460.63406L0.9980.6200.9990.8090.6940.54806R0.9980.7050.9990.8520.7260.50307L0.9970.7220.9980.8600.6780.57407R0.9970.6300.9990.8150.6750.55308L0.9970.8500.9980.9240.6500.50708R0.9970.8670.9970.9320.6410.56609L0.9910.8160.9910.9030.3150.50909R0.9940.9150.9940.9550.3680.49310L0.9940.6720.9950.8340.4080.48310R0.9960.9000.9960.9480.5900.45811L0.9970.9330.9970.9650.6170.62211R0.9980.8780.9980.9380.6560.57012L0.9970.6990.9980.8480.5660.60912R0.9980.7700.9990.8850.6930.51313L0.9980.8330.9980.9160.6670.50313R0.9980.6500.9980.8240.5480.44514L0.9970.8130.9980.9050.6520.63914R0.9970.6920.9990.8450.6280.554Average0.9970.7570.9970.8770.6290.547Figure 11
Visual effects of the best and worst cases from the DRIVE dataset.Figure 12
Visual effects of the best and worst cases from the STARE dataset.Figure 13
Visual results of the best and worst cases from the HRF dataset.Figure 14
Best and worst cases visual results from the CHASE_DB1 dataset.The framework performs well on both healthy and pathological images of all selected datasets. The statistical results in Tables3–6 validates that the suggested system is robust and has the capability to handle the bright lesions images of the STARE dataset, higher resolution images of the HRF dataset, low resolution images of the DRIVE dataset, and left/right eyes images of the CHASE_DB1 dataset. The anatomical structures are also efficiently omitted to avoid any misclassification.The average statistical results of the proposed framework on all selected datasets are displayed in Table7, which reflects that the highest mean score of Acc 0.997, Sn 0.814, Sp 0.997, and AUC 0.905 is achieved on the CHASE_DB1 dataset. The lowest FPR is also observed using the same dataset. The highest value of MCC 0.761 and CAL 0.699 is recorded on the HRF dataset. The highest value of each parameter is italicized in the respective column of the Table 7.Table 7
Mean results comparison of the proposed method on different datasets.
DatasetsImagesAccSnSpAUCMCCCALDRIVE (1st observer)200.9540.7660.9720.8690.7210.690DRIVE (2nd observer)0.9580.7970.9730.8850.7390.696HRF_Normal150.9630.8010.9820.8810.7610.699HRF_DR0.9540.7260.9740.8260.6660.603HRF_Glaucoma0.9590.7900.9790.8810.7030.667HRF_Average450.9590.7720.9780.8630.7100.656CHASE_DB1 (1st observer)280.9970.7570.970.8770.6290.547CHASE_DB1 (2nd observer)0.9960.8140.9960.9050.5690.547STARE200.9960.7920.9970.8950.7070.566The average performance parameter scores of the proposed framework on the DRIVE and STARE datasets are compared with the existing literature in Table8, while Table 9 shows the result comparison of the HRF and CHASE_DB1 datasets. The Acc, Sn, and Sp results of all techniques in Tables 8 and 9 are acquired from their respective published articles while the AUC result is calculated by using the formula in Table 2.Table 8
Performance assessments of existing techniques on the DRIVE and STARE datasets.
MethodYearDRIVESTAREAccSnSpAUCAccSnSpAUCHuman observer0.9470.7790.9720.8740.9350.8950.9380.917Unsupervised techniquesChauduri [31]19890.877——0.788————Zana and Klein [42]20010.9380.697——————Martinez-Perez [45]20070.9340.7250.9650.8450.9410.7510.9550.853Zhang [27]20100.938———0.948———Bankhead [38]20120.9370.7030.9710.8370.9320.7580.9500.854Fraz [44]20120.9430.7150.9760.8450.9440.7310.9680.850Azzopardi [5]20150.9440.7660.9700.8680.9500.7720.9700.871Oliveira [33]20160.9460.8640.9560.9100.9530.8250.9650.895Khan [40]20160.9610.7460.9800.8630.9460.7580.9630.861Biswal [29]20170.9500.7100.9700.8400.9500.7000.9700.835Khan [32]20170.9440.7540.9640.8590.9480.7520.9560.854Soomro [63]20170.9430.7520.9760.8640.9610.7840.9810.883Badawi [3]20180.9550.7910.9710.8810.9530.8650.9610.913Yue [50]20180.9450.7530.9730.863————Soomro [9]20180.9480.7450.9620.8540.9510.7840.9760.880Soomro [64]20180.9530.7520.9760.8640.9670.7860.9820.884Fan [60]20180.9600.7360.9810.8580.9570.7910.9700.880Khan [65]20190.9510.7700.9650.8680.9510.7520.9810.867Memari [59]20190.9610.7610.9810.8710.9510.7820.9650.873Proposed20200.9580.7970.9730.8850.9960.7920.9980.895Supervised techniquesNiemeijer [11]20040.9420.6900.9700.830————Staal [12]20040.9440.7190.9770.8480.9520.6970.9810.839Ricci [14]20070.959———0.964———Lupascu [15]20100.9590.6730.9870.830————MarÃ­n [13]20110.9450.7070.9800.8440.9530.6940.9820.838Wang [16]20150.9770.8170.9730.8950.9810.8100.9790.894Roychowdhury [17]20150.9520.7250.9830.8540.9510.7720.9730.873Aslani [20]20160.9510.7540.9800.8670.9610.7550.9830.869Zhu [18]20170.9610.7140.9870.851————Thangaraj [22]20170.9610.8010.9750.8880.9430.8340.9540.893Memari [23]20170.9720.8720.9880.9300.9510.8090.9790.894Dharmawan [21]2018—0.8310.9720.902—0.7920.9830.887Yan [24]20180.9540.7630.9820.8730.9640.7740.9860.880Guo [25]20190.9550.7800.9810.8810.9660.8200.9830.902Khowaja [26]20190.9750.8180.9710.8950.9750.8240.9750.899Soomro [8]20190.9590.8020.9740.9480.9610.8010.9690.945Soomro [62]20190.9560.8700.9850.9860.9680.8480.9860.988Fan [61]20190.9660.7960.9820.8890.9740.8160.9870.901Table 9
Performance assessments of existing methods with the proposed model on the HRF and CHASE_DB1 datasets.
TechniqueYearHRFCHASE_DB1AccSnSpAUCAccSnSpAUCUnsupervised techniquesOdstrcilik [54]20130.9490.7740.9670.871————Azzopardi [5]2015————0.9390.7590.9590.859Zhang [66]20160.9570.7980.9740.8860.9460.7630.9680.866Biswal [29]2017————0.9400.7600.9700.865Rodrigues [67]20170.9480.7220.9640.843————Badawi [3]2018————0.9530.8000.9640.882Proposed20200.9600.7320.9790.8550.9960.7270.9970.862Supervised techniquesRoychowdhury [17]2015————0.9530.7200.9820.851Thangaraj [22]2017————0.9470.6290.9730.797Memari [23]2017————0.9480.8190.9590.889Dharmawan [21]2018—0.8130.9770.895————Yan [24]2018————0.9610.7640.9810.873Fan [60]2018————0.9510.6570.9730.815Guo [25]2019————0.9630.7890.9800.885Khowaja [26]2019————0.9520.7560.9760.866Soomro [62]20190.9620.8290.9620.9780.9760.8860.9820.985Fan [61]20190.9760.8240.9870.9050.9710.8020.9850.893In Table8, the obtained results of the framework are compared with 19 unsupervised and 18 supervised existing techniques. The proposed framework achieved the highest Acc result than all unsupervised methods on the DRIVE dataset except Khan et al. [40], Memari et al. [59] which is 0.003%, and Fan et al. [60] which is 0.002% better than ours. The supervised methods Ricci and Perfetti [14], Lupascu et al. [15], Wang et al. [16], Zhu et al. [18], Thangaraj et al. [22], Memari et al. [23], Khowaja et al. [26], and Fan et al. [61] show 0.001%, 0.001%, 0.019%, 0.003%, 0.003%, 0.014%, 0.017%, and 0.008% better results than the proposed method, respectively. But some of these methods are only validated on one dataset, which reflects that they are tuned for a single dataset. Some of these methods produce a very low AUC score, which is a trade-off between Sn and Sp. Moreover, supervised methods are computationally very expensive. In the case of the STARE dataset, the framework produced highest Acc scores than all other methods. Table 9 reflects that there are very few techniques that used both HRF and CHASE_DB1 datasets for validation. The Acc score of the framework is higher than both supervised and unsupervised approaches on the HRF and CHASE_DB1 datasets except Soomro et al. [62] and Fan et al. [61] which is slightly higher than ours on HRF dataset only. Fan et al. [61] showed higher Sp value than all other methods on the HRF dataset. The highest Sp value on CHASE_DB1 dataset is obtained by the proposed method. All the other supervised and unsupervised methods acquired a bit greater or equivalent values of Sn and AUC metric on the HRF and CHASE_DB1 datasets as compared to ours.In Table10, the MCC and CAL values are recorded by the proposed method and other existing supervised and unsupervised methods. The MCC and CAL values of Chauduri et al. [31], Niemeijer et al. [11], Hoover et al. [53], and B-COSFIRE [5] are calculated by utilizing their publicly accessible segmented images. The results of Fraz et al. [68, 69], RUSTICO [58], Yang et al. [70, 71], Vega et al. [72], FC-CRF [73], and UP-CRF [73] are extracted from their published articles.Table 10
Performance assessments of existing techniques on the four datasets.
MethodYearDRIVESTAREHRFCHASE_DB1MCCCALMCCCALMCCCALMCCCALUnsupervised techniquesChauduri [31]19890.4200.208——————Hoover [53]2000——0.6150.534————Fraz [68]20110.733—0.700—————Fraz [69]20130.736—0.691—————B-COSFIRE [5]20150.7190.7210.6980.7090.6860.5770.6560.608RUSTICO [58]20190.7290.7280.6980.7090.6910.5870.6630.620Proposed20200.7390.6960.7070.5660.7100.6560.6290.547Supervised techniquesYang [70]20190.736—0.704—0.712———Yang [71]20180.725—0.662—0.682———FC-CRF [73]20160.7560.7310.7270.6580.6900.5410.7040.622UP-CRF [73]20160.7400.6750.7260.6650.6770.4750.6890.571Vega [72]20150.662—0.640—————Niemeijer [11]20040.7220.659——————The average value of MCC attained by the proposed method is higher than all compared unsupervised approaches on the DRIVE, STARE, and HRF datasets, while it is statistically lower than the supervised methods (i.e., FC-CRF [73] and UP-CRF [73]) on the DRIVE, STARE, and CHASE_DB1 datasets. The CAL value of the proposed method is observed higher than all supervised and unsupervised methods on the HRF dataset, while it is statistically lower than or equivalent to CAL values of other methods on the DRIVE, STARE, and CHASE_DB1 datasets.


### 4.3.1. Processing Time

The proposed framework processes a single image in a very short time as equated to other approaches in Table11. The time values are computed on the single image taken from the DRIVE and STARE datasets.Table 11
Processing time evaluation of the systems.
MethodTimeHardware particularsRoychowdhury [17]3.11 secIntel Core i3 CPU 2.6 GHz, 2 GB RAMZhu [18]12.160 sec4.0 GHz Intel i7-4790K CPU and 32 GB RAMMemari [23]8.2 minsIntel i5-M480 CPU, 2.67 GHz, 4 GB RAMBiswal [29]3.3 secIntel i3 (4010U CPU) 1.7 GHz, 4 GB RAMBadawi [3]8 secCPU 2.7 GHz, 16 GB RAMYue [50]4.6 secIntel i5-6200U CPU 2.3 GHz, 8 GB RAMKhan [39]6.1 sec5∗Intel Core i3 CPU, 2.53 GHz, 4 GB RAMKhan [40]1.56 secAzzopardi [5]11.83 secVlachos [47]9.3 secBankhead [38]22.45 secProposed5.5 sec


## 4.3.1. Processing Time

The proposed framework processes a single image in a very short time as equated to other approaches in Table11. The time values are computed on the single image taken from the DRIVE and STARE datasets.Table 11
Processing time evaluation of the systems.
MethodTimeHardware particularsRoychowdhury [17]3.11 secIntel Core i3 CPU 2.6 GHz, 2 GB RAMZhu [18]12.160 sec4.0 GHz Intel i7-4790K CPU and 32 GB RAMMemari [23]8.2 minsIntel i5-M480 CPU, 2.67 GHz, 4 GB RAMBiswal [29]3.3 secIntel i3 (4010U CPU) 1.7 GHz, 4 GB RAMBadawi [3]8 secCPU 2.7 GHz, 16 GB RAMYue [50]4.6 secIntel i5-6200U CPU 2.3 GHz, 8 GB RAMKhan [39]6.1 sec5∗Intel Core i3 CPU, 2.53 GHz, 4 GB RAMKhan [40]1.56 secAzzopardi [5]11.83 secVlachos [47]9.3 secBankhead [38]22.45 secProposed5.5 sec


## 5. Conclusion

Vessel extraction is momentous for inspecting abnormalities inside and around the retinal periphery. The retinal vessel segmentation is a challenging task due to the existence of pathologies, unpredictable dimensions and contour of the vessels, nonuniform clarification, and structural inconsistency between subjects. The proposed methodology is consistent, faster, and completely automated for isolation of retinal vascular network. The success of the proposed framework is evidently revealed by the RVS statistics on the DRIVE, STARE, HRF, and CHASE_DB1 datasets. The eradication of anomalous structures prior to enhancement boosted the efficiency of the proposed method. The application of logical operators avoids misclassification of foreground pixels which enhances the accuracy and makes the method robust. Pictorial representation validates that the framework is able to segment both healthy and unhealthy images. Furthermore, the method does not include any hand-marked data by experts for training, which makes it computationally fast.



---

*Source: 8365783-2020-12-12.xml*