Multifeature Detection of Microaneurysms Based on Improved SSA

: The early diagnosis of retinopathy is crucial to the prevention and treatment of diabetic retinopathy. The low proportion of positive cases in the asymmetric microaneurysm detection problem causes preprocessing to treat microaneurysms as noise to be eliminated. To obtain a binary image containing microaneurysms, the object was segmented by a symmetry algorithm, which is a combination of the connected components and SSA methods. Next, a candidate microaneurysm set was extracted by multifeature clustering of binary images. Finally, the candidate microaneurysms were mapped to the Radon frequency domain to achieve microaneurysm detection. In order to verify the feasibility of the algorithm, a comparative experiment was conducted on the combination of the connected components and SSA methods. In addition, PSNR, FSIM, SSIM, ﬁtness value, average CPU time and other indicators were used as evaluation standards. The results showed that the overall performance of the binary image obtained by the algorithm was the best. Last but not least, the accuracy of the detection method for microaneurysms in this paper reached up to 93.24%, which was better than that of several classic microaneurysm detection methods in the same period.


Introduction
Diabetic retinopathy (DR) is one of the four leading causes of blindness in the world; it causes blindness mainly in working-age adults around the world [1].The other three leading causes are cataract, glaucoma and maculopathy [2].A global assessment of the prevalence of DR in a 2012 study, which used data from population-based studies worldwide, found that DR patients accounted for 34.6% of diabetics [3].Previous studies have demonstrated significant differences between diagnosed and undiagnosed DR in prevalence.An estimated 100 million people in China and 80 million people in India have diabetes [4,5].This means that the early diagnosis of DR is important.The early lesions of DR mainly include microaneurysms (MAs) and rigid exudations, among others.DR is clinically diagnosed with several lesions (MA, hard exudate and soft exudate) at the initial stages [6].Some studies have shown that the presence of even one or two dynamic MAs increases the rate of accelerated disease progression [7,8].Moreover, the number of MAs is significantly different among different graders, but the identification of MAs is relatively consistent after image processing [9].Common fundus diagnostic methods include RGB fundus angiography and fundus fluorescein angiography (FFA).In this paper, the research object was FFA, because FFA images have a unique role in the diagnosis of fundus diseases, especially for the diagnosis of vascular diseases (MA is a disease that acts on fundus blood vessels) [10].MA is a kind of fundus oculi lesion that can be accurately diagnosed using FFA [11].FFA is an angiographic method in which sodium fluorescein is injected into the patient's blood vessels to make the angiography more visible and more easily differentiated from surrounding tissue.FFA was first proposed by Novotny and Alvis in 1961, and it has since contributed to the development of MA detection [12].Aiming to achieve high contrast of MAs in FFA images, Quan Qi segmented the MAs in FFA images by OTSU [13].Manual segmentation of fundus images is low in efficiency, costly and difficult to perform.Many scholars introduced metaheuristic algorithms for multiple optimal threshold values of fitness functions to solve the problem mentioned above.In order to improve the segmentation speed, fundus images were segmented by whale optimization algorithm (WOA), genetic algorithm (GA) and ant colony optimization (ACO) algorithm [14][15][16].The above metaheuristic algorithms (in the literature [14][15][16]) effectively explore from global thresholds to local thresholds.Metaheuristic algorithms are used widely as evolutionary algorithms and swarm intelligence optimization algorithms [17].Intelligent optimization algorithms have excellent performance in the field of image processing, but there is no perfect algorithm to solve all optimization problems [18].Considering the engineering perspective, a faster and more accurate multithreshold image segmentation method is urgently needed.In this paper, DR FFA images were segmented by the salp swarm algorithm (SSA) technique, an algorithm initially proposed by Mirjalili et al. [19].The SSA has advantages in real-time medical image processing because it iterates quickly and easily converges and obtains the global optimal solution.The biggest problem of common population optimization algorithms is that they can easily fall into local optimization.An advantage of the SSA in regard to this problem is that the exploration and exploitation ratio of the algorithm is adjustable, which makes it possible to change update by adjusting the exploration and exploitation ratio.Literally, threshold segmentation is the clustering of pixels according to the threshold value T. Threshold techniques fall into two categories: global threshold methods and local threshold methods.In the former, a common threshold is shared by all pixels of the image [20].In the latter, a uniform threshold is shared by all pixels in a region separated from the image according to certain rules.In many cases, the two methods are combined to segment an image [21].Local techniques are commonly used for complex inputs and irregular lighting and noise distribution.The local threshold method is faster than the global method, but it is not suitable for complex inputs as described above.Therefore, we extended the research field to multiple dimensions.In order to verify the performance of the connected components with SSA (CSSA) for the segmentation of fundus images, a comparative experiment was performed with extreme optimization (EO), simplified swarm optimization (SSO), WOA, particle swarm optimization (PSO) and GA.The experimental results showed that the binary image created by CSSA contained more information, which contributed to the accuracy of subsequent MA detection [22][23][24][25][26]. MA testing is divided into fully supervised learning and semisupervised learning (SSL).The field of fully supervised learning for MA positioning has grown rapidly in recent years.However, because of the limitations of FFA image datasets, fully supervised learning is applied mainly to MA localization of color fundus images.Norah Asir proposed a dual-stream red lesion detection system based on faster RCNN for MA detection [27].Cheng Wan proposed EAD-Net, incorporating the attention mechanism for MA detection [28].SSL is a learning method that combines supervised learning with unsupervised learning.SSL uses a large amount of unlabeled data.It is increasingly important because it requires fewer people to do the work while delivering higher accuracy compared to fully supervised learning.Because FFA has few public datasets, many scholars have used SSL.Gwénolé Quellec located MAs by matching the subbands of a wavelet transform image and a mathematical model of MAs [29].In the literature [30], generative adversarial networks and thermal regression have been used to implement multifeature detection for MA.Marwan D. Saleh segmented h-maxima transformed images with multiple features to achieve MA detection [31].
The main content of this paper is as follows: First, to obtain an asymmetric binary image containing MAs, an FFA image was multidimensionally segmented by combining the connected components with the CSSA.Second, the binary images were multifeature clustered by the MSA-K algorithm to obtain a set of candidate MAs.Third, a preprocessing operation was performed on the candidate MA windows extracted from the original image.Finally, the windows were mapped to the Radon domain with multiple angles θ i (i =1, 2, • • • , m ) to obtain Radon(MA − window j , and the detection of MA was achieved by clustering the standard deviation values between the central band of Radon(MA − window j .
The method of this paper can be implemented on a simple and basic operating platform, which means that the algorithm has low requirements on the operating environment and low economic costs.The multithreshold processing of diabetic retinal images by CSSA has greatly reduced the workload of doctors and reduced the cost of diabetes diagnosis.

Materials and Methods
This paper did not involve human-related experiments.All the pictures in the simulation section were obtained from the Internet.
As shown in Figure 1, the asymmetric MA detection in this paper included binary images, sets of candidate MAs and locations of MAs.MAs were highlighted as dots in FFA images, as shown in Figure 2 [32].The method of this paper can be implemented on a simple and basic operating platform, which means that the algorithm has low requirements on the operating environment and low economic costs.The multithreshold processing of diabetic retinal images by CSSA has greatly reduced the workload of doctors and reduced the cost of diabetes diagnosis.

Materials and Methods
This paper did not involve human-related experiments.All the pictures in the simulation section were obtained from the Internet.
As shown in Figure 1, the asymmetric MA detection in this paper included binary images, sets of candidate MAs and locations of MAs.MAs were highlighted as dots in FFA images, as shown in Figure 2 [32].

SSA
The SSA is a new swarm optimization algorithm that simulates salp foraging behavior in the deep ocean.The mathematical model of the salp chain includes leaders and followers.The former contributes to exploration and the latter contributes to exploitation.The fronts of salp chains are defined as leaders, and the salps for the rest of the chain are defined as followers [19].As this sentence implies, leaders lead the salp chain forward and  The method of this paper can be implemented on a simple and basic operating platform, which means that the algorithm has low requirements on the operating environment and low economic costs.The multithreshold processing of diabetic retinal images by CSSA has greatly reduced the workload of doctors and reduced the cost of diabetes diagnosis.

Materials and Methods
This paper did not involve human-related experiments.All the pictures in the simulation section were obtained from the Internet.
As shown in Figure 1, the asymmetric MA detection in this paper included binary images, sets of candidate MAs and locations of MAs.MAs were highlighted as dots in FFA images, as shown in Figure 2 [32].

SSA
The SSA is a new swarm optimization algorithm that simulates salp foraging behavior in the deep ocean.The mathematical model of the salp chain includes leaders and followers.The former contributes to exploration and the latter contributes to exploitation.The fronts of salp chains are defined as leaders, and the salps for the rest of the chain are defined as followers [19].As this sentence implies, leaders lead the salp chain forward and  The former contributes to exploration and the latter contributes to exploitation.The fronts of salp chains are defined as leaders, and the salps for the rest of the chain are defined as followers [19].As this sentence implies, leaders lead the salp chain forward and explore the search space, and followers lead each other.A salp chain is shown in Figure 3. explore the search space, and followers lead each other.A salp chain is shown in Figure 3.The SSA is a stochastic algorithm to avoid falling into local optima stagnation.The weights of SSA exploration and exploitation can be changed by changing the proportion of leaders and followers.As in other swarm algorithms, elements in the SSA explore the N-dimensional search space.It is assumed that two-dimensional matrixes store the position of the salp, called X, and that of the swarm's target in the search space, represented by the food source, called F.
The equation for the position update is as follows: where i n X is the position of the th i salp follower, N is the number of variables covered with the objective problem, n F is the position of the food for the nth variable and n lb and n ub indicate the lower and upper limitations of the th n variable, respectively.The parameter 1 c balances exploration and exploitation by changing the ratio of leaders to followers as follows: The SSA is a stochastic algorithm to avoid falling into local optima stagnation.The weights of SSA exploration and exploitation can be changed by changing the proportion of leaders and followers.As in other swarm algorithms, elements in the SSA explore the N-dimensional search space.It is assumed that two-dimensional matrixes store the position of the salp, called X, and that of the swarm's target in the search space, represented by the food source, called F.
The equation for the position update is as follows: where X i n is the position of the i th salp follower, N is the number of variables covered with the objective problem, F n is the position of the food for the n th variable and lb n and ub n indicate the lower and upper limitations of the n th variable, respectively.The parameter c 1 balances exploration and exploitation by changing the ratio of leaders to followers as follows: where l is the current iteration number of the algorithm, L is the total iteration number of the algorithm and c 3 determines the forward direction of the chain.When c 3 ≥ 0.5, the chain moves towards +∞, and when c 3 < 0.5, the chain moves towards −∞.c 2 determines the length of the chain's movement.A schematic diagram of a salp chain update is shown in Figure 4.
The initial salp position in the SSA is random, which means that the SSA can explore a problem in multiple directions.The positions of leaders are updated according to Equations ( 1) and ( 2), and the positions of followers are updated according to Equation (3).The SSA is better at exploration than most ACO-like swarm algorithms because the position of its followers is determined only by that of the previous salp.In addition, the location of leaders depends on the food F n .F n is the position of the salp with the best fitness value in the current iteration l.The higher the percentage of leaders, the more random the chain is.In this paper, 30 salps were selected, among which the first 15 were leaders and the last 15 were followers, to avoid falling into local optimal solutions.is the optimal solution; is the swarm's optimal n F for the current iteration l .
The initial salp position in the SSA is random, which means that the SSA can explore a problem in multiple directions.The positions of leaders are updated according to Equations (1) and ( 2), and the positions of followers are updated according to Equation (3).The SSA is better at exploration than most ACO-like swarm algorithms because the position of its followers is determined only by that of the previous salp.In addition, the location of leaders depends on the food n F .n F is the position of the salp with the best fitness value in the current iteration l .The higher the percentage of leaders, the more random the chain is.In this paper, 30 salps were selected, among which the first 15 were leaders and the last 15 were followers, to avoid falling into local optimal solutions.

CSSA
The segmentation threshold of globally connected components depends on the global probability of pixel distribution.In general, global threshold segmentation results in the loss of image details.Experiments have shown that the lesions in DR images are extracted by global threshold methods in a distorted manner.In particular, some blood vessel endings became broken blood vessels under global threshold methods.In order to solve this problem, symmetry CSSA is proposed herein.Each connected component identified by global threshold methods is subdivided into several subregions, different thresholds are calculated for different subregions and subregions are processed according to local thresholds.The SSA is used to segment subgraphs based on multiple thresholds.In this paper, we defined the number of individuals in the salp group K = 30, among which the first 15 were defined as leaders and the last 15 were defined as followers.In order to is the optimal solution; is the swarm's optimal n F for the current iteration l .
The initial salp position in the SSA is random, which means that the SSA can explore a problem in multiple directions.The positions of leaders are updated according to Equations (1) and ( 2), and the positions of followers are updated according to Equation (3).The SSA is better at exploration than most ACO-like swarm algorithms because the position of its followers is determined only by that of the previous salp.In addition, the location of leaders depends on the food n F .n F is the position of the salp with the best fitness value in the current iteration l .The higher the percentage of leaders, the more random the chain is.In this paper, 30 salps were selected, among which the first 15 were leaders and the last 15 were followers, to avoid falling into local optimal solutions.

CSSA
The segmentation threshold of globally connected components depends on the global probability of pixel distribution.In general, global threshold segmentation results in the loss of image details.Experiments have shown that the lesions in DR images are extracted by global threshold methods in a distorted manner.In particular, some blood vessel endings became broken blood vessels under global threshold methods.In order to solve this problem, symmetry CSSA is proposed herein.Each connected component identified by global threshold methods is subdivided into several subregions, different thresholds are calculated for different subregions and subregions are processed according to local thresholds.The SSA is used to segment subgraphs based on multiple thresholds.In this paper, we defined the number of individuals in the salp group K = 30, among which the first 15 were defined as leaders and the last 15 were defined as followers.In order to is the optimal solution; is the swarm's optimal n F for the current iteration l .
The initial salp position in the SSA is random, which means that the SSA can explore a problem in multiple directions.The positions of leaders are updated according to Equations (1) and ( 2), and the positions of followers are updated according to Equation (3).The SSA is better at exploration than most ACO-like swarm algorithms because the position of its followers is determined only by that of the previous salp.In addition, the location of leaders depends on the food n F .n F is the position of the salp with the best fitness value in the current iteration l .The higher the percentage of leaders, the more random the chain is.In this paper, 30 salps were selected, among which the first 15 were leaders and the last 15 were followers, to avoid falling into local optimal solutions.

CSSA
The segmentation threshold of globally connected components depends on the global probability of pixel distribution.In general, global threshold segmentation results in the loss of image details.Experiments have shown that the lesions in DR images are extracted by global threshold methods in a distorted manner.In particular, some blood vessel endings became broken blood vessels under global threshold methods.In order to solve this problem, symmetry CSSA is proposed herein.Each connected component identified by global threshold methods is subdivided into several subregions, different thresholds are calculated for different subregions and subregions are processed according to local thresholds.The SSA is used to segment subgraphs based on multiple thresholds.In this paper, we defined the number of individuals in the salp group K = 30, among which the first 15 were defined as leaders and the last 15 were defined as followers.In order to is the swarm's optimal F n for the current iteration l.

CSSA
The segmentation threshold of globally connected components depends on the global probability of pixel distribution.In general, global threshold segmentation results in the loss of image details.Experiments have shown that the lesions in DR images are extracted by global threshold methods in a distorted manner.In particular, some blood vessel endings became broken blood vessels under global threshold methods.In order to solve this problem, symmetry CSSA is proposed herein.Each connected component identified by global threshold methods is subdivided into several subregions, different thresholds are calculated for different subregions and subregions are processed according to local thresholds.The SSA is used to segment subgraphs based on multiple thresholds.In this paper, we defined the number of individuals in the salp group K = 30, among which the first 15 were defined as leaders and the last 15 were defined as followers.In order to transform the image threshold segmentation problem into a multidimensional solution problem, OTSU was selected as the fitness function of the SSA.
The OTSU algorithm is an efficient image segmentation method for binary images [33].It assumes that an image, with size M × N and gray levels belonging to [0, 1, . . ., L − 1], is divided into two parts (foreground and background) with level T [34].
The variance between classes is: For CSSA, the first step is the extraction of connected components.An FFA image f (m, n) includes the foreground (object region) and background, and the gray value of the pixels belonging to these two ranges is significantly different [35].Threshold processing can be treated as a function: where p(m, n) is the characteristic of the points distributed around the point (m, n) and f (m, n) is the gray level of point (m, n).Extracting connected components avoids this problem.This approach assumes that the representation of the connection components contained in set C is Y and that the point p(m, n) in Y is known.Then, all elements of Y are generated using the Equation expressed iteratively as follows: In this paper, connected components were extracted by the 8 connected neighboring pixels.Connected component sets with sizes of length = [min(C), max(C)] were segmented by different dimensions.To avoid unnecessary calculation, connected components were divided by the 3-dimensional algorithm at most.
The food F L n obtained by the SSA is a 1-3-dimensional solution, based on the F L n segmentation image shown in Figure 5.
Extract the connected components of a DR FFA image Q i (i=1,2, . . .,T); Determine the segmentation dimension based on the size of Q i For n = 1: T; Initialize the salp population While the termination condition is not satisfied; Calculate the fitness of each X i j based on (5); The food n L F obtained by the SSA is a 1-3-dimensional solution, based on the n L F segmentation image shown in Figure 5.

Extraction of the Candidate MA Sets
The main flow of this part is shown in Figure 6.

MSA-K Algorithm
The MSA-K algorithm is a quickly convergent version of the clustering algorithm proposed by Liwei at the 39th CCC International Conference.The steps of the MSA-K algorithm are as follows [36]: First, calculate the centers by salp algorithm according to the adjacent principle; Calculate the distance between object and cluster center; Assign each object to the nearest category; Update the cluster center; Stop the iteration.

Multifeature Clustering
This section mainly concerns locating lesions in binary images obtained by CSSA segmentation.A candidate lesion set was obtained by MSA-K clustering, including the minimum length-width ratio of connected components bounding box (T); the connected component, including the area of the connected component (P), the farthest distance within the connected component (L) and the radius-to-diameter ratio (RD).The window of Radon transforms was then clustered to get the final candidate lesion set.
The blood vessels of binary images obtained by CSSA were easier to eliminate.More accurate candidate lesion sets were extracted by MSA-K clustering of T, P, L and RD.
T, P, L and RD are defined as follows:

The Location of MAs
A 50 50 × template was used to capture MA windows from the original image.Make an inner circle of the window with the pixel outside the circle being 0. This operation is to reduce the interference of the edge pixels in the subsequent multiangle Radon transform.
Radon transform is a mapping that projects an image in the direction of angle θ .The transformation is as follows [37]: Figure 7 contains 4 typical lesions and the corresponding Radon transformation.In particular, (a) is a window that includes only an MA; the MA in (b) is near the blood vessel; and (c) and (d) are severed blood vessels wrongly identified as MAs.
As shown in Figure 7, it was found that the MA corresponded to a uniform and central light band in the Radon varying-frequency band.Therefore, the detection of MAs was achieved by the clustering of the standard deviations of the central Radon band.
In order to verify this point of view, the 3D graphs of a center circle and a diagonal line after Radon transformation are shown in Figure 8.As shown in Figure 8, the Radon transform of the circle was evenly distributed, and the height of the spine was close to unity.The window containing the straight line was mapped to the Radon domain through 18 angles (0°, 10°, …, 180°).As shown in Figure 8b, the spine region was unevenly distributed and had mutations.This phenomenon was consistent with the theorem that the line integral of the same angle as the line is the largest.

The Location of MAs
A 50 × 50 template was used to capture MA windows from the original image.Make an inner circle of the window with the pixel outside the circle being 0. This operation is to reduce the interference of the edge pixels in the subsequent multiangle Radon transform.
Radon transform is a mapping that projects an image in the direction of angle θ.The transformation is as follows [37]: As shown in Figure 7, it was found that the MA corresponded to a uniform and central light band in the Radon varying-frequency band.Therefore, the detection of MAs was achieved by the clustering of the standard deviations of the central Radon band.
In order to verify this point of view, the 3D graphs of a center circle and a diagonal line after Radon transformation are shown in Figure 8.As shown in Figure 8, the Radon transform of the circle was evenly distributed, and the height of the spine was close to unity.The window containing the straight line was mapped to the Radon domain through 18 angles (0 • , 10 • , . . ., 180 • ).As shown in Figure 8b, the spine region was unevenly distributed and had mutations.This phenomenon was consistent with the theorem that the line integral of the same angle as the line is the largest.In order to prove the feasibility of the experiment, the process was reproduced for 20 MA windows and 20 noise windows.The results are shown in Figure 9.In Figure 9, the boundary between the standard deviations of true and false MAs is clear.In order to prove the feasibility of the experiment, the process was reproduced for 20 MA windows and 20 noise windows.The results are shown in Figure 9.In Figure 9, the boundary between the standard deviations of true and false MAs is clear.In order to prove the feasibility of the experiment, the process was reproduced for 20 MA windows and 20 noise windows.The results are shown in Figure 9.In Figure 9, the boundary between the standard deviations of true and false MAs is clear.

Evaluation Parameters and Data Analysis
Peak signal-to-noise ratio (PSNR) was the basic criterion to evaluate the quality of the images [38]: In order to evaluate whether there was damage in the structure of the segmented image, the SSIM parameter was selected as an evaluation criterion in this paper [39].
For images X and Y, the parameter SSIM is defined as follows: FSIM is a series of improvements on the basis of SSIM.The evaluation criterion considers that different pixels in an image have different weights.For example, the pixels at the edge of an object are more important than other pixels in defining the structure of the image [40].FSIM is defined as follows:

Simulation Results
Table 1 is a comparison of SSIM, PSNR and FSIM in the different algorithms.In order to make the correlation among the data clearer, Table 1 was translated into Figure 10.
Table 2 shows the fitness values of the different algorithms.In order to make the data more intuitive, Table 2 can be interpreted with a radar chart.Figure 11 is the radar chart of SSIM, PSNR and FSIM.
The chord graph is mainly used to show the relationships between multiple objects; the connection between two points on the circle shows the relationship between the two corresponding objects, and the width of the line indicates the degree of the proportional relationship between the two objects.The same color in Figure 11a represents the fitness value of a particular algorithm for different dimensions.The same color in Figure 11b represents the fitness value of the same dimension for different algorithms.

Evaluation Parameters and Data Analysis
Peak signal-to-noise ratio (PSNR) was the basic criterion to evaluate the quality of the images [38]: In order to evaluate whether there was damage in the structure of the segmented image, the SSIM parameter was selected as an evaluation criterion in this paper [39].
For images X and Y, the parameter SSIM is defined as follows: FSIM is a series of improvements on the basis of SSIM.The evaluation criterion considers that different pixels in an image have different weights.For example, the pixels at the edge of an object are more important than other pixels in defining the structure of the image [40].FSIM is defined as follows:

Simulation Results
Table 1 is a comparison of SSIM, PSNR and FSIM in the different algorithms.In order to make the correlation among the data clearer, Table 1 was translated into Figure 10.
Table 2 shows the fitness values of the different algorithms.In order to make the data more intuitive, Table 2 can be interpreted with a radar chart.Figure 11 is the radar chart of SSIM, PSNR and FSIM.
The chord graph is mainly used to show the relationships between multiple objects; the connection between two points on the circle shows the relationship between the two corresponding objects, and the width of the line indicates the degree of the proportional relationship between the two objects.The same color in Figure 11a represents the fitness value of a particular algorithm for different dimensions.The same color in Figure 11b represents the fitness value of the same dimension for different algorithms.Since the CSSA needs to calculate the threshold value of each connected component, it is meaningless to compare the CPU time (mean run time) of the CSSA with those of other global threshold algorithms.The computing speed of the CSSA can be evaluated by comparing the CPU time of the SSA with those of other algorithms.In this paper, the run time of all algorithms was the average run time (1-3D) of the same standard (same nonrepresentative code, iterated 1000 times) and segmentation standard on the same equipment (CPU: AMD A8-5550 M; graphics card: AMD Radeon HD 8670 M + AMD Radeon HD 8550 G; RAM: RY: 4 GB).The results are shown in Table 3.
In conclusion, the standard deviation of the window by Radon transformation can be used as the standard for dividing MA.The final result is shown in Figure 12.    3.
In conclusion, the standard deviation of the window by Radon transformation can be used as the standard for dividing MA.The final result is shown in Figure 12.

Original Image
Image by CSSA MA Detection

Discussion
The above experimental results showed that compared with the other two classic MA detection methods, the average accuracy of asymmetric MA detection in this paper was  Table 4 shows the accuracy of the algorithm in this paper and several other algorithms with the highest accuracy in locating fundus lesions.

Discussion
The above experimental results showed that compared with the other two classic MA detection methods, the average accuracy of asymmetric MA detection in this paper was higher.The method of asymmetric MA detection in this article consisted of two parts.The first part was the improvement of the SSA to obtain a binary image containing MAs.This improvement solves the difficulty in medical image recognition that MAs are treated as noise to be canceled in preprocessing and contributes to subsequent MA detection.The second part was clustering the binary image based on T, P, L and RD to extract the candidate MA set.Then MA set was mapped to the Radon frequency domain based on the morphological characteristics of MAs.Finally, MSA-K clustering Radon (MA) was used to detect MAs.Therefore, the experimental verification also contained two parts.
The first was the evaluation of the binary image obtained by the improved CSSA.In order to objectively and quantitatively verify the performance of the binary image obtained by this method, the method was compared with the classical algorithms.The CSSA is an improved multidimensional algorithm based on the SSA for obtaining binary images that contain MAs.Therefore, the performance of the CSSA can be represented by the SSA. Figure 10 visualizes the excellence of the binary images obtained by CSSA in terms of SSIM, PSNR and FSIM.
In this part of the experiment, a radar chart and a chord chart were used to clearly interpret the experimental data.By observing the proportion of the target in Figure 11a, it was found that the fitness value of the three-dimensional algorithm was the best at 66.98%.Correspondingly, the SSA had the highest proportion of overall dimensions in Figure 11b.This means that the performance of the binary images obtained by CSSA was comprehensive optimization.
The other part of the experiment was a comparative experiment between the method developed in this paper and several classic MA detection methods for accuracy.The results showed that the accuracy of the asymmetric MA detection algorithm in this paper was the highest and that the average accuracy could reach 93.24%.
Limitations: compared with similar swarm optimization algorithms, the SSA has faster convergence speed and relatively excellent exploration and exploitation ability, but the limitation of the SSA is that c 1 = 2e −( 4L l ) 2 is related to the preset number of iterations L, which makes the minimum number of iterations for the algorithm to obtain the optimal solution proportional to L. This means that, in order to make the algorithm converge faster, it is necessary to determine the complexity of the research problem in advance and determine the maximum number of iterations L. This necessitates adding one more step: problem preprocessing before solving problems by SSA.This will be addressed in further work.In addition, the detection of MA in color fundus images will be studied in the future.

Conclusions
Since the number of MA pixels has an extremely low global scale, inherent image preprocessing regards MAs as noise to be canceled.Therefore, asymmetric MA detection has always been a difficulty in fundus image processing.In this paper, we solved this problem in three steps.First, we merged multiple connected components and the SSA to obtain a binary image containing MAs.Second, this binary image was multifeature clustered using the MSA-K algorithm to obtain the set of candidate MAs.Finally, asymmetric MA detection was achieved by clustering the standard deviation of the center band in the Radon(window j .
The semisupervised learning method based on the SSA proposed in this paper achieved high accuracy in MA detection; the average accuracy of this method could reach 93.24%.

Figure 1 .
Figure 1.Block diagram of asymmetric MA detection.

Figure 2 .
Figure 2. Enlarged view of an MA.

Figure 1 .
Figure 1.Block diagram of asymmetric MA detection.

Figure 1 .
Figure 1.Block diagram of asymmetric MA detection.

Figure 2 .
Figure 2. Enlarged view of an MA.

Figure 2 .
Figure 2. Enlarged view of an MA.

2. 1 .
Extraction of Binary Images 2.1.1.SSA The SSA is a new swarm optimization algorithm that simulates salp foraging behavior in the deep ocean.The mathematical model of the salp chain includes leaders and followers.

Figure 4 .
Figure 4. Schematic diagram of a salp chain searching for optimal solutions in two dimensions.is the optimal solution;is the swarm's optimal n

Figure 4 . 16 Figure 4 .
Figure 4. Schematic diagram of a salp chain searching for optimal solutions in two dimensions.

16 Figure 4 .
Figure 4. Schematic diagram of a salp chain searching for optimal solutions in two dimensions.is the optimal solution;is the swarm's optimal n

Figure 7
Figure 7 contains 4 typical lesions and the corresponding Radon transformation.In particular, (a) is a window that includes only an MA; the MA in (b) is near the blood vessel; and (c) and (d) are severed blood vessels wrongly identified as MAs.As shown in Figure7, it was found that the MA corresponded to a uniform and central light band in the Radon varying-frequency band.Therefore, the detection of MAs was achieved by the clustering of the standard deviations of the central Radon band.In order to verify this point of view, the 3D graphs of a center circle and a diagonal line after Radon transformation are shown in Figure8.As shown in Figure8, the Radon transform of the circle was evenly distributed, and the height of the spine was close to unity.The window containing the straight line was mapped to the Radon domain through 18 angles (0 • , 10 • , . . ., 180 • ).As shown in Figure8b, the spine region was unevenly distributed and had mutations.This phenomenon was consistent with the theorem that the line integral of the same angle as the line is the largest.

Figure 7 .Figure 8 .
Figure 7. Candidate lesion windows and Radon transforms.It contains 4 typical lesions and the corresponding Radon transformation.In particular, (a) is a window that includes only an MA; the MA in (b) is near the blood vessel; and (c,d) are severed blood vessels wrongly identified as MAs.

Figure 7 . 17 TrueFigure 7 .Figure 8 .
Figure 7. Candidate lesion windows and Radon transforms.It contains 4 typical lesions and the corresponding Radon transformation.In particular, (a) is a window that includes only an MA; the MA in (b) is near the blood vessel; and (c,d) are severed blood vessels wrongly identified as MAs.

Figure 8 .
Figure 8. Eighteen-angle Radon transform 3D diagram: (a) is the Radon transform of a circle; (b) isthe Radon transform of an image containing a diagonal line.

Symmetry 2021 , 16 Figure 9 .
Figure 9.The standard deviation distribution of 18-angle Radon transformation, where o represents true MA and * represents false MA.

Figure 9 .
Figure 9.The standard deviation distribution of 18-angle Radon transformation, where o represents true MA and * represents false MA.

Figure 11 .
Figure 11.The chord graph of Table 2. (a) is the chord diagram of each algorithm for three dimensions, and (b) is the chord diagram of each dimension for the algorithms.
Since the CSSA needs to calculate the threshold value of each connected component, it is meaningless to compare the CPU time (mean run time) of the CSSA with those of other global threshold algorithms.The computing speed of the CSSA can be evaluated by comparing the CPU time of the SSA with those of other algorithms.In this paper, the run time of all algorithms was the average run time (1-3D) of the same standard (same nonrepresentative code, iterated 1000 times) and segmentation standard on the same equipment (CPU: AMD A8-5550 M; graphics card: AMD Radeon HD 8670 M + AMD Radeon HD 8550 G; RAM: RY: 4 GB).The results are shown in Table

Figure 12 .
Figure 12.Locations of lesions in DR images.Left, the original image; middle, the binary image obtained by the method of this paper; right, the location of the final lesion.
Symmetry 2021, 13, x FOR PEER REVIEW 3 of 16 the connected components with the CSSA.Second, the binary images were multifeature clustered by the MSA-K algorithm to obtain a set of candidate MAs.Third, a preprocessing operation was performed on the candidate MA windows extracted from the original image.Finally, the windows were mapped to the Radon domain with multiple angles

Table 1 .
Comparison of SSIM, PSNR and FSIM in the different algorithms.

Table 1 .
Comparison of SSIM, PSNR and FSIM in the different algorithms.

Table 1 .
Comparison of SSIM, PSNR and FSIM in the different algorithms.

Table 2 .
Comparison of the fitness values of each algorithm (shows 4 figures after the decimal point).

Table 2 .
Comparison of the fitness values of each algorithm (shows 4 figures after the decimal point).

Table 3 .
The mean run time of algorithm (shows 4 figures after the decimal point; unit: minutes).

Table 4
shows the accuracy of the algorithm in this paper and several other algorithms with the highest accuracy in locating fundus lesions.

Table 4 .
Accuracy of lesion localization.

Table 3 .
The mean run time of algorithm (shows 4 figures after the decimal point; unit: minutes).

Table 4 .
Accuracy of lesion localization.