key: cord-1031136-ei8crm1q authors: Chakraborty, Shouvik; Mali, Kalyani title: A radiological image analysis framework for early screening of the COVID-19 infection: A computer vision-based approach date: 2022-02-03 journal: Appl Soft Comput DOI: 10.1016/j.asoc.2022.108528 sha: 8af890345c3349dcd62d648bf1b153e4e7118ab8 doc_id: 1031136 cord_uid: ei8crm1q Due to the absence of any specialized drugs, the novel coronavirus disease 2019 or COVID-19 is one of the biggest threats to mankind Although the RT-PCR test is the gold standard to confirm the presence of this virus, some radiological investigations find some important features from the CT scans of the chest region, which are helpful to identify the suspected COVID-19 patients. This article proposes a novel fuzzy superpixel-based unsupervised clustering approach that can be useful to automatically process the CT scan images without any manual annotation and helpful in the easy interpretation. The proposed approach is based on artificial cell swarm optimization and will be known as the SUFACSO (SUperpixel based Fuzzy Artificial Cell Swarm Optimization) and implemented in the Matlab environment. The proposed approach uses a novel superpixel computation method which is helpful to effectively represent the pixel intensity information which is beneficial for the optimization process. Superpixels are further clustered using the proposed fuzzy artificial cell swarm optimization approach. So, a twofold contribution can be observed in this work which is helpful to quickly diagnose the patients in an unsupervised manner so that, the suspected persons can be isolated at an early phase to combat the spread of the COVID-19 virus and it is the major clinical impact of this work. Both qualitative and quantitative experimental results show the effectiveness of the proposed approach and also establish it as an effective computer-aided tool to fight against the COVID-19 virus. Four well-known cluster validity measures Davies–Bouldin, Dunn, Xie-Beni, and [Formula: see text] index are used to quantify the segmented results and it is observed that the proposed approach not only performs well but also outperforms some of the standard approaches. On average, the proposed approach achieves 1.709792, 1.473037, 1.752433, 1.709912 values of the Xie-Beni index for 3, 5,7, and 9 clusters respectively and these values are significantly lesser compared to the other state-of-the-art approaches. The general direction of this research is worthwhile pursuing leading, eventually, to a contribution to the community. Automated computer-aided systems prove their effectiveness and real-life applicability in various scenarios. Automated systems have a diverse domain of applications and sometimes, these systems are inevitable to perform certain jobs efficiently and in a cost-effective and highly time-bound manner. This domain is evolving day-byday and continuous effort can be observed from various researchers to enhance this domain. Computer-assisted systems can be categorized in two ways. The first one is the supervised approach in which some properly annotated data are required to perform the classification and interpretation job [1, 2] . Therefore, these automated systems are dependent on the ground truth data typically produced by some domain experts. But, it may not be always possible to acquire the properly annotated ground truth data due to the involvement of human experts [3] . Sometimes, some cases are not well-defined or not seen earlier, and therefore, it is very difficult to get some ground truth data for those cases. Unsupervised systems can be helpful in this context because these systems are not dependent on the ground truth data and can automatically explore some patterns from the underlying dataset by utilizing the surrounding knowledge [4] [5] [6] [7] . So, the unsupervised approaches are helpful in those situations where a sufficient amount of properly annotated ground truth data are not available. The unsupervised computer-aided systems are widely applied in different domains of research [8, 9] . Biomedical image analysis is no exception and exploits the advantages of unsupervised automated systems in various phases. Radiology is one of the important and frequently used parts of the biomedical imaging domain which is serving as an important tool for noninvasive diagnostic systems. X-Ray, CT Scan, etc. are widely used nowadays, to study the internal functionalities and the present state of the different organs [10] [11] [12] . Automated systems are helpful to analyze and diagnose different patients automatically and automated radiological image analysis systems are also helpful in preparing precise and timely reports by reducing the human intervention and also reducing some unintentional human-made errors. Physicians, radiological technicians, and all other concerned domain experts can be significantly benefitted from the advancement in the field of computer-aided radiological image analysis systems. Apart from the automated analysis of the radiological images, computer-assisted systems can be helpful in parameter tuning of the image acquisition hardware, image preprocessing, quality control, selecting the appropriate level of radiation, and many more. Therefore, automated systems can act as a helping hand in the decision-making process. In table 1 some of the related biomedical image segmentation works of literature are discussed which is helpful in a better understanding of the current trend and status of the same. Apart from these works, some comprehensive studies can be found in [13] [14] [15] [16] [17] . This approach is used to segment the PET volumes based on varying source-to-background (S/B) ratios which are collected from the phantom of a body. The calibrated source-to-background curves are used to determine the volume using the iterative thresholding procedure. One major drawback of the system is that it cannot effectively measure small volumes. Wiemker et. al. [19] Thresholding CT scan images This approach segments the CT scan images to easily interpret and study the lung nodules. This work proposes a divergence theorem and histogram-based Ct image segmentation approach. This approach is can effectively and optimally isolate the lung nodules from the CT scan images. In this context, the optimality is defined in terms of the mean gradient of the iso-surface and the sphericity. Asari et. al. [20] Thresholding and differential region growing This work is targeted to extract gastrointestinal lumen from the endoscopic images. This algorithm is consisting of two stages where the first stage employs a global thresholding approach and in the second phase, the differential region growing is used to extract the gastrointestinal lumen from the endoscopic images. The dynamic hill-clustering approach is used to ascertain the effectiveness of the termination criteria and to look after the growth process. Yu-qian et. al. [21] Edge detection This work is addressing the problem of edge detection in the presence of noise. Traditional gradient-based edge detection approaches are susceptible to noise and therefore, this approach proposes a novel approach to detect edges of the lung CT scan images using mathematical morphology. This approach is tested on the CT images which are corrupted with the salt-and-pepper noise and its efficiency is proved by comparing this approach with some of the other standard approaches. It is observed that this approach can efficiently reduce the effect of noise and also can generate precise edges. Falcao et. al. [22] Shortest-path based method This work is based on the computation of the shortest path using Dijkstra's algorithm. This approach is highly dependent on the user intervention to efficiently determine the segmented regions and to define the objects. This approach is found to be 3 to 15 times faster compared to manual tracing. This approach can be applied almost independently to the applications. One main problem associated with this method is the difficulties associated with the choice of slabs and orthogonal slices which has a serious impact on the efficiency of this approach. Pan et. al. [23] highlighted as the paths of the bacteria. Although this approach performs well and comparative study shows the effectiveness of the proposed approach still, one problem of this approach is not very robust to noise. Noise can lead to crumpled edges. This approach is not also suitable to handle overlapped cells. Ji et. al. [24] Fuzzy Cmeans clustering Syntheic, MR images, natural images This work proposes a modification to the traditional fuzzy Cmeans clustering by addressing some of the problems. Traditional fuzzy C-means clustering approach does not consider the spatial information and less robust to noise. This work proposes a modification which is known as the weighted image patch-based FCM. In this work, pixels are replaced with the weighted patches which is helpful to incorporate spatial information in the segmentation process. It is helpful to increase the reliability of the overall segmentation process but it also increases the computational overhead drastically. Agrawal et. al. [25] Optimum boundary point detection This approach is devoted to segment the intracranial region from the magnetic resonance images. This work proposes a novel hybrid approach which is based on the genetic algorithm and the bacterial foraging algorithm. The combination of these two approaches is used to optimize the objective function of the fuzzy c-means clustering. The final cluster centers are obtained using a method called optimum boundary point detection. This approach cannot determine the optimal number of clusters automatically and produces inaccurate results if the predefined clusters and the actual number of clusters differ. Chaira [26] Fuzzy Cmeans clustering This work proposes a new approach to segment CT images. This approach is based on intuitionistic fuzzy set theory and it is known as the intuitionistic fuzzy C means clustering. In this work, a novel objective function which is known as intuitionistic fuzzy entropy is incorporated with the traditional fuzzy C-means clustering. This approach is applied to different CT scan images to prove its efficiency. Miao et. al. [27] dictionary learning and improved fuzzy Cmeans clustering Synthetic, MRI, CT Scan In this work, a noise-resistant version of the fuzzy c-means clustering algorithm is proposed and applied to segment the images. This approach can be divided into two phases where the first phase incorporates a dictionary learning method to handle the noise. In the second phase, this dictionary learning approach is hybridized with the improved fuzzy c-means clustering approach. The proposed approach is not efficient for medical images with inhomogeneous intensity distribution. Apart from these works, some of the most recent and relevant works can be found in [28] [29] [30] [31] [32] that can be referred to, to understand the further advancements of this domain. In this context, it is worth mentioning here that the active contour model is an effective way of image segmentation. There are several variations available of this approach. The traditional active contour approach was proposed in 1988 [33] . A modified version of the traditional active contour approach is proposed in [34] and it is known as geometric active contours. This approach uses gradient information of an image to construct the edge stop function. A region information-based approach is proposed in [35] . This approach is developed by Chan and Vese and this is a parametric representation. Some deep learning approaches are developed that use the loss function of the active contour model as their loss function [36] . Some recent developments like the MAC model [37] , SBGFRLS model [38] , LSAC model [39] , etc. can also be observed in this domain. The highly infectious coronavirus disease 2019 or COVID-19 creates a worldwide pandemic scenario. Although the mortality rate is not very high, the highly infectious nature of this virus is the main threat to society. Due to the absence of any specialized drug, it is very difficult to restrict the drastic spread of this virus. Apart from using various protective equipment, early detection and isolation can be very effective to combat the spread of this highly infectious virus. In the middle of this pandemic scenario, some vaccines are invented and are being applied to the people and it is a ray of hope to fight against this virus. As per the report of the world health organization, 239,437,517 numbers of confirmed cases can be observed in 216 countries and 4,879,235 people are already expired due to this disease as of 15 th October 2021, 4:32pm CEST [40] . From these statistics, it is clear that the worldwide mortality rate is approximately 2.0378% which is not a very large value. The major risk factor lies in the highly infectious nature of this virus. Hopefully, 6,495,672,032 vaccine doses have already been administered worldwide which may be helpful in reducing the mortality rate. Many countries are not prepared with the appropriate infrastructures to support COVID-19 infected patients. Moreover, many people from remote areas are not even able to arrange protective gear like masks, sanitizers, etc. The reverse transcription-polymerase chain reaction test i.e., RT-PCR test is the only test available to date to confirm the presence of the COVID-19 virus. Some researches show that CT scan images of the chest region are showing some signs of the early COVID-19 infection [41] . It is a quite inspiring finding because CT scan images can be used to isolate some suspected patients at an early phase and therefore, the drastic spread of this virus can be stopped to some extent. The CT scan images cannot replace the RT-PCR test because some false negatives are reported in [42, 43] . The screening of the COVID-19 positive patients using the CT scan images are recommended in [44] . The presence of some prominent features like ground-glass opacities, crazy paving, etc. (which are given in table 2) helps us to trace the initial presence of this infection, and image segmentation is an essential task to automate the screening process. Typically, the absence of properly annotated data makes the automated biomedical image analysis job difficult. These are the basic motivation behind proposing a novel radiological image segmentation approach SUFACSO (SUperpixel based Fuzzy Artificial Cell Swarm Optimization). As the name suggests, the proposed approach is based on the superpixels and type 2 fuzzy systems where the type 2 fuzzy objective function is modified to incorporate the advantages of superpixels to efficiently process a large amount of spatial information. The fuzzy objective function is optimized with the recently developed metaheuristic procedure i.e., artificial cell swarm optimization. The proposed method allows automated and efficient analysis of the CT scan images which is beneficial to enhance the computer-aided diagnostic systems to act as a tool against the COVID-19 virus. To summarize, the major contributions are as follows:1) A novel superpixel-based image segmentation technique is proposed that reduces the incurred computational cost for processing a high amount of spatial information, 2) Type-II fuzzy system is incorporated with the superpixel-based approach, 3) A recently developed metaheuristic procedure ACSO is further enhanced, 4) The conventional fitness function of the FCM clustering approach is enhanced to exploit the advantages of superpixel 5) The cluster centers are updated with the help of the proposed fuzzy ACSO approach. The remaining article is prepared in the following way: section 2 and 3 describes the artificial cell swarm optimization method and the type 2 fuzzy clustering framework respectively. Section 4 and 5 describe the proposed SUFACSO approach and the obtained results respectively. Section 6 discusses some of the relevant points and a brief conclusion is presented in section 7. Table 2 : Significant properties which are found in the CT scan images of the chest region of COVID-19 positive patients [45] . Percentage of the observed samples ground-glass opacities (GGO) 100% multilobe and posterior involvement 93% bilateral pneumonia 91% subsegmental vessel enlargement (> 3 mm) 89% This is a recently developed metaheuristic procedure that is inspired by the artificial cell division procedure. The artificial cell swarm optimization procedure mimics the artificial cells as the search agents. The actual artificial cell division approach [46] is slightly modified to design the optimization procedure. The incorporated modifications are listed below [47] : i. The artificial cells are not depending on the current state to participate in the cell division process. ii. The artificial cells only take part in the cell division process The hierarchical tree structure is formed throughout the generations due to the artificial cell division process. Swarms of artificial cells are considered in the optimization process to take part in the artificial cell division A huge number of swarms can significantly increase the fitness evaluations and a small number of swarms can increase time to converge and therefore, is essential to decide the swarm count moderately. In this work, the swarm count is considered is a fixed parameter. One artificial cell can produce some new cells and the production of new cells occur at a certain distance which is inversely dependent on the fitness of the producer cell as expressed in equation 2. 1 1 . The distance between the th k cell and any of the th l cell, which are produced from the same parent cell, must be same. Therefore, if a cell is near to global optima, then it can generate some other cells at a smaller distance and vice-versa. Smaller steps help to search the nearest portions of the global optima cautiously so that the global optima may not be missed accidentally. A cell does not have any effect on the population once its lifespan is over. This property helps to maintain the size of the population and prevents getting overpopulated. The successor cells of a cell can produce some other cells by the cell division process to maintain the population. The life span of a cell can belong it belongs to the near-optimal area. The quality of a population is evaluated using the lambda function which is given in equation 3. Algorithm 1 illustrates the artificial cell swarm optimization approach in brief [47] . The proposed approach adopts the type 2 fuzzy logic-based clustering approach to effectively model and handle the random uncertainties. In most real-life applications, the uncertainty cannot be predicted in advance. A wide range of input types can produce random uncertainties. Hence, it is essential to cope up with the random uncertainties in real-life scenarios. The fuzzy C-means clustering approach is one of the widely used clustering approaches which is suitable to various problems of different domains [48] [49] [50] [51] . The main reason behind the increasing popularity of fuzzy systems is the suitability of this approach in different scenarios where the crisp clustering approaches do not perform well. A single point can be a member of more than one cluster at the same time with some membership values. The total sum of all membership values for a certain point must be one. So, the value of the membership can be anything between 0 and 1. The dissimilarity function which is optimized by the fuzzy C-means clustering approach is given in equation 5. The type 2 fuzzy logic systems use separate sets of membership values that are also fuzzy in nature. This approach allows efficient modeling of dynamic input uncertainties by providing additional degrees of freedom. In this work, the type 2 fuzzy logic-based clustering approach is adopted to overcome some of the common problems of type 1 fuzzy systems like noise sensitivity, relative membership values, etc., and also to handle uncertainties well [52] . It is essential to improve the outcome of the segmentation process. The uncertainty of a point must be decided depending on the membership value i.e., if a point has a membership value of 1 then its uncertainty will be certainly nil. So, a lower membership value indicates higher uncertainty and vice-versa. Some of the basic reason behind the adoption of type 2 fuzzy system in this work is listed below [53] : a. A point with higher uncertainty has a lesser impact on the overall clustering process and vice-versa. It helps to achieve more realistic results. b. Better noise handling capability can be achieved The membership value in type 2 fuzzy systems can be calculated using equation 8 and the cluster centers can be updated using equation 9. The proposed approach does not require equation 9 and can update the cluster centers. The artificial cell swarm optimization process will guide the proposed approach to determine the optimal cluster centers. The accuracy clustering process can be determined by a small threshold value  . The type 2 fuzzy clustering system can be easily understood from algorithm 2 and the schematic diagram of the type 2 fuzzy system can be visualized from figure 1. The ever-growing technology allows us to increase the quality of the image acquisition hardware. High-quality biomedical images can be acquired from various biomedical image acquisition devices and it is helpful in a precise analysis of the biomedical images. Automated biomedical image analysis devices are facing some challenges due to the increasing quality of biomedical images. A high amount of spatial information creates severe problems for automated and computer-aided diagnostic systems because medical diagnostic systems demand quick and accurate results. Image segmentation plays a vital role in many automated computer-aided image analysis systems. It is essential to generate precise reports within the stipulated amount of time to provide accurate treatment to the patients. To handle this situation effectively and to accelerate the screening process of the COVID-19 infection, a superpixel-based novel approach is proposed in this work to segment the CT scan images. Superpixels are useful to represent a set of pixels in a computation-friendly manner. Different approaches can be found in the literature to find the superpixel image from an input image [54] [55] [56] . Some superpixel computation methods like mean shift [54] and watershed [56] produce irregular superpixels and some methods like SLIC [55] generate regular superpixels. Meanshift and watershed approaches are more useful due to the capability to generate irregular superpixels. The watershed approach is simpler to implement compared to the mean-shift approach but it is sensitive to the noise which is not at all desirable for the image segmentation approaches. In this work, the noise sensitivity of the watershed approach is removed with the help of the gradient image, which is generated using the approach, proposed in [57] . The obtained gradient image is processed using the morphological erosion and dilation-based reconstruction operations, which are given in equations 10 and 11 respectively. In these equations, In the above equations, img and img denotes the original image and the marker image and the img can be expressed by equation 15 and 16.  and  are the two operators to compute the point wise maximum and the minimum values. Here, se is the structuring element and it is an important parameter that controls the segmented outcome. The size of the structuring element is subjective and depends on the image under consideration. Practically, it is not possible to determine different structuring elements which are of various sizes, depending on the image. Therefore, the pointwise maximum value is computed (using equation 17) from the gradient images, which are generated by applying more than one structuring elements where the number of structuring elements is decided as per the range of the size controlling parameter A very small lower bound is not desirable because it will produce very small regions and some essential edge information can be lost. A small threshold value  is used to control the error rate and the upper threshold value as given in equation 18 . A higher value of  indicates a higher error rate but, a smaller upper bound that helps to achieve lesser computational overhead. So, the threshold value can be adjusted as per the requirement and depending on the available resources. The conventional fuzzy C-means clustering approach often overlooks some important spatial information that can be costly in terms of the segmentation performance. Some approaches try to solve this problem by considering and blending some local spatial information in the objective function but it increases the computational cost and therefore not suitable on many occasions. Superpixels can help in this context by over-segmenting an image in many small, perceptually uniform, and homogeneous regions. In this work, the CT images are first processed to determine the superpixels using the proposed approach and then the fuzzy artificial cell swarm optimization approach is used to determine the segmented image by finding the optimal clusters. As discussed earlier, the type 2 fuzzy system is used to perform the segmentation. The fuzzy objective function which is given in equation 5, is optimized by the artificial cell swarm optimization algorithm. To incorporate the advantages of the superpixel, it is necessary to modify the fuzzy objective function. The fuzzy objective function which is given in equation 5 deals with the pixel-wise spatial information. The concept of superpixel represents a group of a pixel using a single value   , as given in equation 19 where, cntPix  is the count of the pixels in the  region and cnt  number of such regions exists. The representative value is used in the objective function, and the modified objective function is given in equation 20 to completely exploit the advantage of the superpixels. In this modified objective function, one superpixel is considered as a single unit and the value   is used to represent the superpixel. proposed procedure is given in algorithm 3 and the schematic flow diagram is given in figure 4 . Input: The input image. Output: Optimal segmented output image 1: Apply the method described in [57] Compute the fitness of each cell. 9: Generate the successor cells using the ACD process by following the successor distance 10: Assign the value of LS for each cell. Modify the LS value of the parent and also update the successor distance. The performance evaluation and comparison of the proposed SUFACSO approach are presented in this section. As discussed earlier, the properly annotated ground truth segmented images may not be available always, and therefore, some standard intrinsic cluster evaluation methods are used here to evaluate the proposed approach quantitatively. Davies Table 3 gives a brief overview of the test images and the test images along with their histograms are given in figure 5 . Figure 5 : The CT scan images and their histograms. The experiments are performed in the MatLab R2014a on a computer that is equipped with an Intel i3 processor and 4GB main memory. The proposed method is compared with some metaheuristic optimization-based image segmentation approaches like modified genetic algorithm [67] , modified PSO [68] , improved bat algorithm [69] and modified cuckoo search method [70] in both qualitative and quantitative manner. The visual comparison is presented in figure 6 where the 01 Test image is considered. The segmented output images which are obtained by applying the proposed SUFACSO approach, are reported in figure 7 . From the qualitative and quantitative results, it can be observed that the proposed SUFACSO approach outperforms some state-of-the-art works and can produce realistic outputs that are certainly helpful for the interpretation of the real-life CT scan images and therefore, this approach can be helpful for the early screening purposes. At the end of each table, the average performance of the five approaches is reported which is beneficial to understand the overall performance of these methods for the different number of clusters and different cluster validity indices. In the case of average, the column-wise optimal values are highlighted instead of highlighting the row-wise optimal values. The row-wise highlighted values talk about the performance of the individual algorithm for the different number of clusters whereas the column-wise highlighted values help to understand the performance of the individual algorithms. It can be observed that the proposed approach outperforms other approaches for most of the number of clusters as well as for most of the validity indices. For example, on a total of 16 occasions (i.e., 4 validity indices x 4 different number of clusters), the proposed approach is found to perform better 11 times. These comparative results are graphically presented in figure 8 . The experiments are carried out for the different numbers of clusters. A particular approach may perform well for a particular cluster count. That is why the average values of all experiments are reported at the end of each table for better interpretation. It can be observed that the proposed approach can optimize different objective functions effectively. Actually, the experiments are carried out on 200 CT images (in the first phase) and 100 CT images (in the second phase). It is already mentioned in section 5.1. Results that are obtained from all images are not possible to report in this stipulated amount of space. Therefore, only some results that are obtained from some selected images are reported. Apart from these tests, the proposed approach is also compared with some of the active contour models based on some standard parameters like accuracy, precision, and recall. This comparison is performed by using the database that is available at [71] . This dataset contains 100 CT scan images with dimensions 512x512. This dataset is created by collecting sample images from 49 patients with age range 32-86 years. The obtained average results are reported in table 8. The rate of convergence is an important parameter to be studied. The performance evaluation remains incomplete without studying and comparing the convergence of different algorithms. The convergence analysis gives a clear view of the comparative performance of different algorithms for the different numbers of clusters. The graphical analysis of the convergence is presented in this subsection using the image 01 Test for the Dunn index. In figure 9 , five different plots are given for five different methods. In a single plot, four separate curves are indicating four different clusters. These curves show that the proposed approach can efficiently segment the images for a higher number of clusters. Moreover, the proposed approach also outperforms some other methods in terms of convergence besides quantitative and qualitative performance. The time complexity is an important aspect that is to be analyzed. From the detailed discussion of the proposed approach, it can be noticed that the proposed approach can be viewed as a two-phase procedure where the watershed-based computation approach is used to determine the superpixel image from the underlying image in the first phase and the optimal segmented outcome is computed in the second phase. The task of optimization is performed using the proposed fuzzy ACSO approach. The gradient information of an image is used to avoid the noise sensitivity of the water-shed based superpixel computation process. The watershed-based technique is a simple method to compute the superpixel and the implementation follows linear complexity [56] . It is quite inspiring and lucrative to adopt this approach on different occasions. In the optimization part, the fuzzy objective function is optimized by using the proposed fuzzy ACSO method. The ACSO approach is an effective and efficient approach that can be executed in linear time [47] . So, the proposed approach is efficient enough and can be effectively used in various real-life problem-solving scenarios. The proposed SUFACSO approach is basically an unsupervised clustering approach that is used for image segmentation purposes. This approach can effectively process high-quality images with the help of the proposed superpixel-based approach that is an essential quality for the real-life application of an image segmentation approach. This approach removes the dependency of choice of the initial cluster centers as well as the ACSO approach determines the optimal cluster centers by optimizing some validity indices. These advantages motivate us to apply the proposed approach to automatically segment the radiological images that will be certainly helpful in diagnosing some symptoms of COVID-19. The experimental outcomes show the efficiency of the proposed approach. Under this pandemic environment, this work is designed hoping that it can help physicians and other domain experts to some extent in the early diagnosis of the disease. Early diagnosis can prevent the drastic spread of this highly infectious virus. The quantitative outcomes of the proposed SUFACSO approach are useful to assess the comparative performance. Quantitative results do not have any direct implications in real-life diagnosis. The segmented outcomes are useful in the diagnosis process. Physicians can investigate the segmented outcomes to find some prominent and common features as mentioned in table 2. The segmented images will be helpful in the easy interpretation of the radiological images. The proposed SUFACSO approach is an efficient image segmentation approach that can effectively segment the radiological images that highly useful in the easy interpretation of these images. The presence of the COVID-19 virus can only be confirmed with the help of some standard tests and one of the most popular and gold-standard tests is the RT-PCR test. Typically, the test reports of the RT-PCR tests are generated within 2-4 days. There is a high possibility that a suspected patient can spread the disease in the community completely unwillingly. The proposed approach can reduce this chance because an initial screening can be performed by the physicians comfortable with the help of the proposed SUFACSO approach. It is worth mentioning here that the proposed approach is neither a replacement of the RT-PCR test nor it can confirm the presence of the virus accurately. However, this approach can be helpful in an initial screening at an early stage that will restrict the spread of this highly infectious virus by separating suspected patients from the rest of the community. 6. Discussion The obtained results indicate that the proposed approach is suitable for real-life scenarios and also performs efficiently. This approach can be easily adapted for the automated screening purposes of the COVID-19 infected patients. It is assumed the quality of the CT scan images is considerably high and the performance of the proposed approach is not verified against the presence of noise. It will be interesting to study the proposed approach in the presence of noise. The scalability of the proposed approach to different types of biomedical images can be explored in future studies. Missing manual annotations can jeopardize the generalizability of the proposed work. On the other hand, the obtained results are quite promising and encouraging. From the best of the knowledge of the authors, there is no publicly available manually annotated dataset for the chest CT scan images of the COVID-19 positive cases. Although the proposed approach is efficient enough to segment the CT scan images automatically and produces realistic segmented outcomes still, some important drawbacks can be observed in this proposed approach that can be addressed in the subsequent works. One important drawback of the proposed approach is that it cannot automatically determine the number of clusters and it can be overcome in future works. Automated estimation of the clusters can make this approach more realistic, robust, and application friendly. The proposed method can handle only a single objective at a time. Therefore, the proposed approach is not suitable for multi-objective optimization issues unless enhanced further. The number of images in the dataset is not very large. So, the proposed approach can also be tested on some additional CT images of COVID-19 infection as well as on some standard dataset of the biomedical images. The proposed SUFACSO approach is an unsupervised segmentation approach. It neither use any training dataset nor uses any pre-trained model. The proposed approach can effectively segment the radiological images that are collected from different patients i.e., not only COVID-19 infected samples but samples collected from patients with other infections as well as normal patients. It is to be clarified that this approach cannot take any decision about the type of disease automatically. For example, the proposed approach cannot automatically differentiate between COVID-19 related lung images and other lung diseases. This approach aims to help physicians in early and quick interpretation of the radiological images and diagnosis of the diseases without any manual delineations. This article proposes a novel, simple and elegant solution that uses some of the important features of the chest CT scan images to screen the COVID-19 suspected patients easily and at an early phase which can be considered as an effective tool to reduce the drastic spread of this virus. From figure 8 , it can be observed that the proposed approach works well in most situations and outperforms most of the other standard approaches. Both qualitative and quantitative study produces some satisfactory results which help to make the proposed approach trustworthy so that it can be reliably adapted in the real-world scenarios. From figure 9 , it can be observed that the proposed approach performs well in terms of convergence. The proposed approach initially performs a superpixel-based clustering using the proposed superpixel computation method which significantly reduces the computational overhead for the further clustering process by reducing a large amount of spatial information. Therefore, radiological images can be conveniently explicated with the application of the proposed method and the proposed approach is also helpful in the easy interpretation of the radiological images. The proposed work neither claims that the suggested approach is cent percent accurate in determining the COVID-19 infection nor claims that it can be a replacement of the RT-PCR test but, the proposed method can help detect some common characteristics from the CT scan images, that can help to isolate some suspected patients from the rest of the community. The proposed approach is helpful for the early screening of the COVID-19 besides being a significant contribution to the image segmentation literature.  The proposed SUFACSO approach can be used as a tool for early screening the suspected COVID-19 infected patients using CT scan images.  The proposed SUFACSO approach exploits the advantages of superpixel to effectively process spatial information.  A newly developed metaheuristic optimization procedure namely artificial cell swarm optimization is investigated.  There is no dependency on the choice of initial cluster centers.  Some standard metaheuristic methods are outperformed by the proposed SUFACSO approach in terms of both qualitative and quantitative manner.  The convergence of the SUFACSO approach is better than some other standard metaheuristic optimization methods.  Experiments are carried out using various cluster validity indices to test the efficiency and the obtained Learning full pairwise affinities for spectral segmentation An Overview of Biomedical Image Analysis From the Deep Learning Perspective A Study on Different Edge Detection Techniques in Digital Image Processing Machine learning applications in genetics and genomics Interrelated two-way clustering: An unsupervised approach for gene expression data analysis An Advanced Approach to Detect Edges of Digital Images for Image Segmentation SUFMACS: a machine learning-based robust image segmentation framework for covid-19 radiological image interpretation Application of Multiobjective Optimization Techniques in Biomedical Image Multi-Objective Optim An optimized intelligent dermatologic disease classification framework based on IoT An integrated interactive technique for image segmentation using stack based seeded region growing and thresholding Contrast Optimization using Elitist Metaheuristic Optimization and Gradient Approximation for Biomedical Image Enhancement An integrated method for automated biomedical image segmentation A survey on image segmentation using metaheuristicbased deformable models: State of the art and critical analysis A Survey of Graph Cuts/Graph Search Based Medical Image Segmentation Intelligent Computing in Medical Imaging: A Study Liver segmentation from CT scans: A survey A survey on deep learning in medical image analysis Segmentation of PET volumes by iterative image thresholding Optimal thresholding for 3D segmentation of pulmonary nodules in high resolution CT A fast and accurate segmentation technique for the extraction of gastrointestinal lumen from endoscopic images Medical images edge detection based on mathematical morphology A 3D generalization of user-steered live-wire segmentation Cell image segmentation using bacterial foraging optimization Fuzzy c-means clustering with weighted image patch for image segmentation A study on fuzzy clustering for magnetic resonance brain image segmentation using soft computing approaches A novel intuitionistic fuzzy C means clustering algorithm and its application to medical images Local segmentation of images using an improved fuzzy C-means clustering algorithm based on self-adaptive dictionary learning A morphology-based radiological image segmentation approach for efficient screening of COVID-19 An Unsupervised Fuzzy Clustering Approach for Early Screening of COVID-19 from Radiological Images Unsupervised domain adaptation based COVID-19 CT infection segmentation network Unsupervised segmentation of COVID-19 infected lung clinical CT volumes using image inpainting and representation learning Lung Ultrasound Segmentation and Adaptation Between COVID-19 and Community-Acquired Pneumonia Snakes: Active contour models Geodesic Active Contours Active contours without edges Learning active contour models for medical image segmentation Magnetostatic active contour model Active contours with selective local or global segmentation: A new formulation and level set method A fast and robust level set method for image segmentation using fuzzy clustering and lattice boltzmann method COVID-19) Dashboard | WHO Coronavirus (COVID-19) Dashboard With Vaccination Data Essentials for Radiologists on COVID-19: An Update-Radiology Scientific Expert Panel Chest CT Findings in Coronavirus Disease-19 (COVID-19): Relationship to Duration of Infection Correlation of Chest CT and RT-PCR Testing in Coronavirus Disease 2019 (COVID-19) in China: A Report of 1014 Cases Sensitivity of Chest CT for COVID-19: Comparison to RT-PCR Chest CT Features of COVID-19 in Artificial cell division Artificial Cell Swarm Optimization Image segmentation using fuzzy clustering: A survey Fuzzy C-means (FCM) clustering algorithm: A decade review from Fuzzy Electromagnetism Optimization (FEMO) and its application in biomedical image segmentation SuFMoFPA: A superpixel and meta-heuristic based fuzzy image segmentation approach to explicate COVID-19 radiological images Typical Characteristics-based Type-2 Fuzzy C-Means Algorithm A type-2 fuzzy C-means clustering algorithm Mean shift: A robust approach toward feature space analysis SLIC superpixels compared to state-ofthe-art superpixel methods Watershed superpixel Finding Contours of Hippocampus Brain Cell Using Microscopic Image Analysis COVID-19 pneumonia | Radiology Case | Radiopaedia.org A Cluster Separation Measure A validity measure for fuzzy clustering Well-separated clusters and optimal fuzzy partitions Segmentation of remotely sensed images with fuzzy thresholding, and quantitative evaluation COVID-19 pneumonia | Radiology Case | Radiopaedia.org COVID-19 pneumonia | Radiology Case | Radiopaedia.org COVID-19 pneumonia | Radiology Case | Radiopaedia.org COVID-19 pneumonia | Radiology Case | Radiopaedia.org A modified genetic algorithm for distributed scheduling problems Vehicle routing problem with uncertain demands: An advanced particle swarm algorithm Improved bat algorithm with optimal forage strategy and random disturbance strategy Modified cuckoo search algorithm in microscopic image segmentation of hippocampus COVID-19 -Medical segmentation The authors would like to express their gratitude and thank the editors, anonymous reviewers, and referees for their valuable comments and suggestions which are helpful in further improvement of this research work. The experiments are performed in a laptop PC which is equipped with 4 gigabytes of RAM, 1 gigabyte of dedicated graphics memory, Intel Core i3 -3217U processor with 1.8 gigahertz clock speed. The dedicated graphics memory is not utilized for any kind of processing purposes. The system in which the experiments are carried out is equipped with the Microsoft Windows 7 (64 bit) operating system. The proposed SUFACSO approach is coded in the Matlab R2014a environment. It is not at all essential to use the Matlab environment to implement the proposed approach. We have chosen Matlab due to the availability of some inbuilt functions which are helpful to reduce the coding complexity. Still, any other languages or platforms can be used to implement the same. The graphs are also generated in the Matlab environment where the' best' position is selected to accommodate the legend without interfering the actual plot. The proposed SUFACSO approach is applied to the CT scan images which are collected from the chest region of the COVID-19 infected patients. It is assumed that there are no manual annotations available and the proposed approach is capable to process the images without having any prior knowledge. The proposed SUFACSO approach produces the optimal segmented image by computing the optimal cluster centers. The final segmented images are constructed by assigning the superpixel to their corresponding cluster centers. These segmented images are helpful to interpret different features from these radiological images.