A Novel Brain MRI Image Segmentation Method Using an Improved Multi-View Fuzzy c-Means Clustering Algorithm

Background: The brain magnetic resonance imaging (MRI) image segmentation method mainly refers to the division of brain tissue, which can be divided into tissue parts such as white matter (WM), gray matter (GM), and cerebrospinal fluid (CSF). The segmentation results can provide a basis for medical image registration, 3D reconstruction, and visualization. Generally, MRI images have defects such as partial volume effects, uneven grayscale, and noise. Therefore, in practical applications, the segmentation of brain MRI images has difficulty obtaining high accuracy.

Materials and Methods: The fuzzy clustering algorithm establishes the expression of the uncertainty of the sample category and can describe the ambiguity brought by the partial volume effect to the brain MRI image, so it is very suitable for brain MRI image segmentation (B-MRI-IS). The classic fuzzy c-means (FCM) algorithm is extremely sensitive to noise and offset fields. If the algorithm is used directly to segment the brain MRI image, the ideal segmentation result cannot be obtained. Accordingly, considering the defects of MRI medical images, this study uses an improved multiview FCM clustering algorithm (IMV-FCM) to improve the algorithm’s segmentation accuracy of brain images. IMV-FCM uses a view weight adaptive learning mechanism so that each view obtains the optimal weight according to its cluster contribution. The final division result is obtained through the view ensemble method. Under the view weight adaptive learning mechanism, the coordination between various views is more flexible, and each view can be adaptively learned to achieve better clustering effects.

Results: The segmentation results of a large number of brain MRI images show that IMV-FCM has better segmentation performance and can accurately segment brain tissue. Compared with several related clustering algorithms, the IMV-FCM algorithm has better adaptability and better clustering performance.

Introduction

The brain is the human body’s nerve center, controlling people’s thinking, memory, speech, and movement, and plays a role in regulating human organs. When the brain is healthy, it can work quickly and efficiently. However, when something goes wrong, the result can be devastating. In recent years, the problem of brain diseases has become increasingly prominent due to comprehensive factors such as high pressure in people’s lives, fast-paced activities, extreme tension in thoughts and emotions, frequent accidents, and serious aging in the population, which continue to threaten people’s health. In clinical medicine, doctors usually use MRI technology to diagnose brain diseases. MRI is a very advanced medical imaging technology. It visualizes the structure and function of the human body through radiology and is particularly suitable for brain tissue research. By using this technology, high soft tissue contrast can be obtained, and it has the advantages of noninvasiveness, nonradiation, and high-precision spatial resolution (Ji et al., 2011; Jiang et al., 2012; Thanh and Wu, 2013).

Neuroscience researchers need to segment brain MRI images to quantitatively study brain diseases. B-MRI-IS marks each pixel in the image as the corresponding brain tissue anatomical structure, such as the thalamus, hippocampus, and ventricles. Then, the segmented tissue size, shape, location, and other characteristics are used to evaluate and formulate medical plans. Experiments have shown that the abnormal shape or volume of certain anatomical regions of the brain is related to brain diseases, such as Alzheimer’s disease and Parkinson’s disease. The quality of B-MRI-IS determines the reliability of researchers’ assessment of brain diseases. The segmentation quality is determined by the segmentation method, so the segmentation method is very important in the entire diagnosis process. Usually, accurate segmentation results of brain tissue are obtained by manual segmentation by experienced brain experts. However, when faced with a large number of datasets, manual segmentation methods become quite expensive, time consuming, and impractical. Moreover, due to differences in experience and knowledge among experts, the segmentation results are not uniform. Therefore, the use of more accurate B-MRI-IS algorithms to quantitatively analyze the volume shape of each tissue in brain MRI images has become the focus of medical image research.

Traditional image segmentation methods include threshold-based methods, clustering-based methods, region-based methods, edge-based methods, and graph theory-based methods (Gu et al., 2019; Kumar et al., 2020). Considering the problems in the segmentation process, scholars have proposed their solutions. Among them, fuzzy clustering algorithms are widely used. The literature (Bezdek, 1981) proposed the FCM algorithm for the first time, introducing the concept of ambiguity into the clustering method (Li et al., 2014), and proposed an optimized bias field estimation and tissue segmentation method. Through the simultaneous iteration of bias field estimation and tissue segmentation, the final segmentation accuracy is improved, and the algorithm is very robust to initialization. The literature (Chuang et al., 2006) proposed the FCM algorithm based on spatial relationships. The algorithm integrates the spatial information of the sum of membership degrees of all points in the window into the traditional FCM. Experiments show that the improved algorithm has certain robustness to noise. Noise interference has always been one of the important interference factors in medical image segmentation. Due to the complexity of the human brain tissue structure and the imaging mechanism of MRI technology, MRI images will show phenomena such as strong noise, false images, and weak boundaries. For example, traditional active contour models often only use target edge information, which leads to premature stopping when processing strongly noisy images, and weak boundary leakage occurs when the boundary is blurred. Many studies use grayscale information and spatial information to suppress noise interference. The literature (Ahmed et al., 2002) introduces local spatial terms into the objective function of FCM to make the algorithm perform better. The literature (Chen and Zhang, 2004) replaces the traditional Euclidean distance with the core distance to make the algorithm more robust. The literature (Krinidis and Chatzis, 2010) proposed a new blur factor and introduced the kernel distance, which is more robust to noise. The literature (Gong et al., 2013) proposed a kernel-based adaptive regularization FCM algorithm. This method introduces adaptive parameters and uses mean filtering, median filtering, and custom filtering of images. There is no need to calculate parameters for each iteration, which greatly saves time and cost. Additionally, the robustness of the algorithm to noise is improved. The literature (Elazab et al., 2015) improves the quality and compression rate of the reconstructed image by compressing the image in blocks.

Figure 1 shows the motivation for this research. To reduce the influence of noise in the segmentation process, the accuracy of B-MRI-IS is improved. This article uses an improved multiview FCM algorithm. The work in this paper is summarized as follows:

FIGURE 1

Figure 1. The motivation for this research.

1. The B-MRI-IS process is proposed. First, the original brain MRI image is preprocessed, such as denoising. Second, the histogram of oriented gradient (HOG) feature, entropy feature, gradient feature, and contrast feature of the original image are extracted. Each feature data point is used as view data to construct multiview data (MVD). Third, the constructed MVD is input into the clustering model to obtain the final clustering result. The clustering result is the image segmentation result.

2. Based on multiview FCM, the IMV-FCM algorithm is proposed. IMV-FCM uses a view weight adaptive learning mechanism so that each view obtains the optimal weight according to its cluster contribution. The final division result is obtained through the view ensemble method. Under the view weight adaptive learning mechanism, the coordination between various views is more flexible, and each view can be adaptively learned to achieve better clustering effects.

3. To demonstrate the effectiveness of the clustering model used in this paper, single-view clustering and multiview clustering algorithms are used for experimental comparison. To demonstrate the antinoise of the algorithm, after adding noise to the experimental image, each algorithm is used for clustering comparison. The results show that the algorithm used in this paper is the most robust to noise compared with other comparison algorithms and can accurately segment brain tissue.

Backgrounds

MRI Technology

In medicine, human brain tissue includes the cerebrum, cerebellum, and brainstem, each of which has relatively independent characteristics. At the same time, the skull structure can be divided into white matter (WM), gray matter (GM), and cerebrospinal fluid (CSF) (Wang et al., 2015). GM contains WM, and CSF fills the crooked groove formed by the GM folds. The brain tissue images show that the gray levels of various brain tissues change slowly and are not constant. Coupled with the influence of various noises, it is easy to cause the gray values to cross and overlap the internal brain tissues. The only correspondence is that one gray value may correspond to several brain internal tissues. In addition, internal brain images with different imaging morphologies also provide different information, and they appear as different gray levels on the images. According to the image gray level, it can be divided into four types, namely, white, off-white, gray, and black. Table 1 shows the characteristics of different gray levels.

TABLE 1

Table 1. Image grayscale characteristics of adult brain tissue.

The proton density ρ, longitudinal relaxation time T1, and transverse relaxation time T2 of brain tissue constitute the high-dimensional biophysical property space of the brain. The T1 and T2 parameters of normal and pathological tissues in different organs of the human body are fixed, and there are certain differences between them. Additionally, different tissues have different densities, and the differences in these properties between brain tissues are the biological basis of the solvability of brain tissue segmentation. Table 2 shows the relevant attributes of the main brain tissues in MRI images.

TABLE 2

Table 2. Related attributes of main brain tissues in MRI images.

Clustering Algorithm

These attributes are determined by factors such as the composition of the brain tissue and the local microstructure. Cluster analysis aims to make reasonable classifications according to the characteristics of the sample and divide the data points with similar characteristics into one category. The idea of the cluster analysis algorithm comes from taxonomy. As MRI images contain increasing information, the requirements for classification are increasing. People have introduced multivariate analysis into taxonomy, forming cluster analysis. The goal of cluster analysis is to collect data to classify when the data points are similar (Bezdek, 1980; Hall et al., 1992; Jiang et al., 2014; Qian et al., 2015, 2016, 2018; Zheng et al., 2015; Jiang et al., 2019; Jiang et al., 2020). Cluster analysis is a technique to study the relationship between sample data points. The clustering results reveal the internal connections and differences between sample data and provide an important basis for data processing. The cluster analysis process is shown in Figure 2.

FIGURE 2

Figure 2. Cluster analysis flowchart.

The clustering object refers to the grayscale collection of image pixels. The similarity measure is a key step used to measure the similarity between data points. Euclidean distance is often used as the similarity measure. The clustering algorithm uses similarity measures to select data points that are very close to each other and classify them into a cluster. Additionally, the result of clustering will affect the selection of subsequent objects and the similarity measurement, so there is a feedback process. The last step of clustering is to evaluate the clustering results. The importance of evaluating the clustering results will not be lower than the importance of the previous steps. For example, in a cluster analysis method based on a threshold value, it is necessary to judge the classification result to guide the selection of the threshold value. In the cluster analysis method with the function as the criterion, the judgment or evaluation of the classification result can be used to guide the selection of the appropriate criterion function.

Clustering algorithms can be divided into partition-based (Jing et al., 2007; Zhu et al., 2009; Yu et al., 2012; Qian et al., 2017), hierarchy-based (Hall and Goldgof, 2011; Devi and Setty, 2018), density-based (Sunjana and Azizah, 2020; Yin et al., 2020), grid-based (Wang et al., 2019; Zhang et al., 2020), and model-based (Lee et al., 2018; Chrobak et al., 2020) algorithms. Due to the complexity of brain medical images, the boundaries and boundaries of different tissues in the brain are uncertain and unclear. This requires a method for reducing the unclearness. Fuzzy mathematics based on fuzzy sets can solve these uncertain and imprecise problems. The fuzzy clustering algorithm in the clustering algorithm is most suitable for B-MRI-IS. The typical fuzzy clustering algorithm is the FCM algorithm. The principle of image segmentation based on FCM is to use the uncertainty of the classification of image pixels. The degree of membership is used to describe this uncertainty, the distance relationship between pixels is described according to the objective function, and the best cluster center is selected. To describe the clustering algorithm more clearly, Table 3 shows the symbols used in the FCM algorithm.

TABLE 3

Table 3. Description of related symbols in FCM.

Assume the image pixel dataset is X = {x₁,x₂,…,x_N}, where x_i represents the gray value of the image pixel. FCM’s idea is to transform the image segmentation process into the optimization process of the feature function when the objective function converges to realize the fuzzy classification of pixel data. That is, the image segmentation problem is transformed into a clustering problem of dividing N pixels into C classes, and the cluster center of each class is expressed as Z = {z₁,z₂,…,z_N}. The functional expression of FCM is

$J_{FCM}=\sum _{i=1}^C\sum _{j=1}^N{u}_{ij}^m{d}^2(x_j,z_i)(1)$

$d^2(x_j,z_i)={||x_j-z_i||}^2(2)$

The size of m determines the clustering ambiguity. If the value of m is too large, the fuzziness of clustering will increase accordingly, which is not conducive to the defuzzification of the model. In contrast, the smaller the value of m is, the smaller the clustering ambiguity, and the segmentation result is close to hard segmentation. In practical applications, the value of m is usually taken as 2. Equations (3) and (4) give the iterative formula of cluster center z_i and membership degree u_ij.

$z_i=\frac{\sum _{j=1}^N{u}_{ij}^m{x}_j}{\sum _{j=1}^N{u}_{ij}^m}(3)$

$u_{ij}=\frac{1}{\sum _{r=1}^C\left[{\left(\frac{d(x_j,z_i)}{d(x_j,z_r)}\right)}^{\frac{2}{m-1}}\right]}(4)$

The steps of the FCM-based image segmentation algorithm are as follows: (1) Input the brain MRI image; (2) set the fuzzy factor m = 2, the maximum number of iterations ε (ε > 0), the number of clusters C, and initialize the cluster center and membership matrix randomly; (3) according to Eqs 3 and 4, update the clustering center z and membership degree u; (4) make the objective function converge until the cluster center no longer changes; and (5) output the segmented image. The flow of B-MRI-IS based on FCM is shown in Figure 3.

FIGURE 3

Figure 3. B-MRI-IS process based on FCM.

Multiview Clustering Algorithm

Traditional fuzzy clustering cannot achieve ideal results when processing multiview data. Based on this requirement, many multiview-related algorithms have been proposed. The existing multiview-related clustering algorithms include coclustering (Gu and Zhou, 2009), collaborative clustering algorithms (Pedrycz, 2002), and multitask clustering algorithms (Chen et al., 2013). There are many multiview clustering algorithms, the most classic of which is CoFKM. Next, we will briefly introduce the principle and implementation steps of the algorithm. Table 4 gives the definition of each symbol of the algorithm.

TABLE 4

Table 4. Description of related symbols in the multiview clustering algorithm.

Equation (5) gives the function expression of CoFKM:

$J_{C\text{o}FKM}=\sum _{v=1}^V\sum _{i=1}^C\sum _{j=1}^N{\tilde{u}}_{ij,v,\mathrm{\eta }}{||x_{j,v}-z_{i,v}||}^2(5)$

${\tilde{u}}_{ij,v,\mathrm{\eta }}=\left(1-\mathrm{\eta }\right)u_{ij,v}^m+\frac{\mathrm{\eta }}{V-1}\sum _{v^{\prime }=1,v^{\prime }\ne v}^V{u}_{ij,v^{\prime }}^m(6)$

$\{\begin{array}{c}\hfill u_{ij,v}\in [0,1],1\le i\le C,1\le j\le N\hfill \\ \hfill \sum _{i=1}^C{u}_{ij,v}=1,1\le v\le V\hfill \end{array}(7)$

where ${\tilde{u}}_{ij,v,\mathrm{\eta }}$ represents the weighted fusion of the membership degree $u_{ij,v}^m$ under the current view and the membership degree $u_{ij,v^{\prime }}^m$ of the other views. According to Eq. 5, through the Lagrangian extremum solution method, the calculation formulas of the clustering center and membership degree are obtained. Equations (8) and (9) give the calculation formulas of the cluster center and membership matrix of the vth view.

$z_{i,v}=\frac{\sum _{j=1}^N{\tilde{u}}_{ij,v,\mathrm{\eta }}{x}_{i,v}}{\sum _{j=1}^N{\tilde{u}}_{ij,v,\mathrm{\eta }}},i=1,2,\mathrm{\dots },C(8)$

$u_{ij,v}=\frac{1}{\sum _{h=1}^C\left[\frac{\left(1-\mathrm{\eta }\right)d_{ij,v}^2+\frac{\mathrm{\eta }}{V-1}\sum _{v^{\prime }=1,v^{\prime }\ne v}^V{d}_{ij,v^{\prime }}^2}{\left(1-\mathrm{\eta }\right)d_{hj,v}^2+\frac{\mathrm{\eta }}{V-1}\sum _{v^{\prime }=1,v^{\prime }\ne v}^V{d}_{hj,v^{\prime }}^2}\right]}(9)$

Through Eqs 8 and 9, the division matrix for any view is calculated. The final division result is obtained using the ensemble method of the geometric mean. The expression of the integrated method is

$\stackrel{⌢}{U}=\sqrt[V]{\prod _{v=1}^V{U}_v}(10)$

Algorithm Evaluation Index

The goal of segmentation is to divide pixels with certain characteristic properties (such as grayscale and spatial position) in an image into a region. All the pixels in a region have similar properties, and the characteristic properties of pixels in different regions are quite different. The evaluation index of the effectiveness of image segmentation is used to describe the pros and cons of the segmentation results. The commonly used evaluation indicators include segmentation accuracy JS, Dice coefficient (DSC), similarity (KI), segmentation coefficient (V_pc), segmentation entropy (V_pe), and pixel error rate (ME). The specific descriptions of commonly used evaluation indicators are shown in Table 5.

TABLE 5

Table 5. Evaluation index descriptions.

Experimental Dataset

Common brain image databases are mainly as follows: Brain Perfusion Database, Allen Brain Atlas, BRATS database, simulated brain MR image data from the Brain Web website, and the Internet Brain Segmentation Repository (IBSR). The description of each database is shown in Table 6.

TABLE 6

Table 6. Brain image database.

Brain MRI Image Segmentation Based on Multiview Fuzzy Clustering

Brain MRI Image Segmentation Process Based on Multiview Clustering

MVD carries more information than single-view data. This MVD feature is beneficial to the improvement of cluster analysis. For clustering analysis of MVD, the traditional single-view clustering algorithm is to first separate each view data from the MVD. Second, the single-view clustering algorithm is used to process the data of each view, and the clustering results of the data of each view are obtained. Finally, a suitable ensemble strategy is used to fuse the clustering results of all views to obtain the final division result. The single-view clustering algorithm for processing MVD has the following shortcomings. When the clustering results of a certain view have obvious deviations or the clustering results between different views are very different, if each view is artificially separated for analysis in this way, the final division result obtained by integrated learning is likely to be poor or the performance of the algorithm is unstable. To make full use of the information carried by MVD and improve the accuracy of clustering, this research introduces multiview technology into the single-view clustering method. The multiview clustering algorithm enables the collaborative learning of various views in the clustering process, makes full use of the data information of each view, and improves the clustering performance and stability of the algorithm. The B-MRI-IS process based on multiview clustering is shown in Figure 4.

FIGURE 4

Figure 4. B-MRI-IS process based on multiview clustering.

As shown in Figure 4, first, the original brain image is preprocessed, such as denoising. Second, the HOG feature, entropy feature, gradient feature, and contrast feature of the original image are extracted. The feature data serve as view data. Third, all view data are input into the multiview clustering model to obtain the final clustering result. Due to the different amounts of information contained in different features, the contribution of the formed view data to the clustering result is also different. Therefore, how to effectively adjust the weight of each view is very important during the use of the multiview clustering algorithm.

IMV-FCM

The CoFKM introduced in the section “Multiview Clustering Algorithm” assigns the same weight to each view, which shows that the CoFKM believes that the contribution of each view data to the clustering is the same. This is obviously inconsistent with the actual situation. For MVD in real production and living environments, the contribution of each view data to the clustering is necessarily different. If the angle of view data with a high degree of contribution is given a high weight and the angle of view data with a low degree of contribution is given a low weight, the overall clustering performance can be improved. Therefore, a multiview fuzzy clustering method with adaptive viewing angle adjustment capability is used in this study, and its objective function is

$\begin{array}{c}J=\sum _{v=1}^V\sum _{i=1}^C\sum _{j=1}^N\sum _{t=1}^V{w}_{v,t}{u}_{ij,t}^m{d}_{ij,v}^2+\mathrm{\gamma }\sum _{v=1}^V\sum _{t=1}^V{w}_{v,t}\log (w_{v,t})\hfill \\ \sum _{i=1}^C{u}_{ij,v}=1,u_{ij,v}\in [0,1]\hfill \\ 1\le i\le C,1\le j\le N,1\le v\le V\hfill \\ \sum _{t=1}^V{w}_{v,t}=1,w_{v,t}\in [0,1],1\le v\le V\hfill \end{array}(11)$

where $d_{ij,v}^2={||x_{j,v}-z_{i,v}||}^2$ and U_v = [u_ij,v] is the division matrix corresponding to the ith view. ${\tilde{u}}_{ij,v}=\sum _{t=1}^V{w}_{v,t}{u}_{ij,t}^m$ is the division fusion item, which realizes the view fusion of different views in the vth view clustering task. At this time, the importance of each view is reflected by the weight w_v,t. w_v,t represents the importance of the tth view in the vth view cluster. W is the weight matrix of all views. The Lagrange multiplier method is used to solve the extreme value of Eq. 11, and the expression of each variable is obtained as follows:

$z_{i,v}=\frac{\sum _{j=1}^N\sum _{t=1}^V{w}_{v,t}{u}_{ij,t}^m{x}_{j,v}}{\sum _{j=1}^N\sum _{t=1}^V{w}_{v,t}{u}_{ij,t}^m}(12)$

$u_{ij,t}=\frac{1}{\sum _{h=1}^C{\left[\frac{\sum _{v=1}^V{w}_{v,t}{d}_{ij,v}^2}{\sum _{v=1}^V{w}_{v,t}{d}_{hj,v}^2}\right]}^{1/(m-1)}}(13)$

$w_{v,t}=\frac{\exp \left[\frac{-\sum _{i=1}^C\sum _{j=1}^N{u}_{ij,t}^m{d}_{ij,v}^2}{\mathrm{\gamma }}\right]}{\sum _{g=1}^V\exp \left[\frac{-\sum _{i=1}^C\sum _{j=1}^N{u}_{ij,g}^m{d}_{ij,v}^2}{\mathrm{\gamma }}\right]}(14)$

The following ensemble method is used to obtain the final partition matrix:

$\overline{U}=\max \left(\sum _{t=1}^V\left(\frac{{\sum }_{v=1}^V{w}_{v,t}}{{\sum }_{t=1}^V{\sum }_{v=1}^V{w}_{v,t}}{U}_t\right)\right)(15)$

Figure 5 shows the schematic diagram of the IMV-FCM algorithm. As shown in Figure 5, after the multiview data with v views are executed by the IMV-FCM algorithm, v membership matrices are obtained. Each membership matrix is assigned a corresponding weight w, and the final global membership matrix is obtained according to the fusion strategy of Eq. 15.

FIGURE 5

Figure 5. Schematic diagram of IMV-FCM algorithm.

The specific steps of the IMV-FCM algorithm are in Algorithm 1.

ALGORITHM 1

Algorithm 1. IMV-FCM algorithm.

Simulation Experiment Analysis

Experimental Background

The dataset used in this study was the data downloaded from BrainWeb. The 2D images of 89 cross sections (89 CS), 92 cross sections (95 CS), and 95 cross sections (95 CS) in the T1-weighted brain MRI image were selected for segmentation. The comparison algorithms used were FCM, CoFKM (Cai et al., 2019), two-layer automatic weighted clustering algorithm (TW-k-means) (Singh et al., 2020), multitask-based K-means (CombKM) (Chen et al., 2013), and collaborative clustering based on sample and feature space (coclustering) (Gu and Zhou, 2009). In the experiment, the iteration stop threshold ε of each algorithm is set to 0.001, and the maximum number of iterations l was set to 100. The parameter settings of each algorithm are shown in Table 7. The parameter selection methods of each comparison algorithm include the grid search method and clustering index. This article mainly uses the grid search method to confirm the optimal parameters. The reason for choosing the grid search method for parameter selection is that the performance of the model trained based on the parameters selected by this method is the best. The evaluation index adopts JS and DSC. All experiments were based on the MATLAB 2018a programming environment, which was simulated on a Lenovo PC configured with a Windows 10 operating system, 2.60 GHz CPU, and 8G memory Intel(R) Core(TM)i7 processor.

TABLE 7

Table 7. Parameter setting of various algorithms based on grid search.

Experimental Results and Analysis

To verify the superiority of the IMV-FCM algorithm, the experiment section gives a comparison of the effect of each algorithm on image segmentation. In analyzing the antinoise performance of the algorithm, this study adds different proportions of Gaussian noise to the image. The segmentation results of each algorithm are shown in Tables 8, 9. The symbol A represents WM, B represents GM, and C represents CSF.

TABLE 8

Table 8. JS indicators of each algorithm (100%).

TABLE 9

Table 9. DSC indicators of each algorithm (100%).

From the experimental data in Tables 8, 9, it can be determined that when processing original brain images, the segmentation performance of the multiview algorithm used in this study is significantly better than FCM. This shows that the introduction of the multiview mechanism can mine information from multiple views, thereby improving the clustering accuracy. Compared with other multiview clustering algorithms, the segmentation effect of the IMV-FCM algorithm in this study is better. This shows that the IMV-FCM algorithm with partition adaptive fusion capability requires the Co-FKM algorithm and CombKM algorithm to manually set the degree of partition fusion. The IMV-FCM algorithm can obtain better view weights during multiview learning and finally obtain better multiview clustering effects.

As the image noise increases, the clustering performance of all clustering algorithms begins to decline. However, even if the performance of each algorithm is declining, the performance of the IMV-FCM algorithm is still better than that of the other algorithms. This is manifested in two aspects. First, the clustering performance of IMV-FCM is better than that of the comparison algorithm. The second is that with the increase in image noise, the clustering performance of the IMV-FCM algorithm decreases the most slowly, as shown in Figures 6, 7. This fully shows that IMV-FCM is more robust to noise.

FIGURE 6

Figure 6. JS indicator drop rate of each algorithm.

FIGURE 7

Figure 7. DSC index drop rate of each algorithm.

Conclusion

To further improve the performance of B-MRI-IS, the IMV-FCM algorithm is applied in this research. The IMV-FCM algorithm is a multiview fuzzy clustering algorithm with adaptive view weight. Under the adaptive learning of the view fusion weight matrix, the coordination between various views is more flexible. Additionally, each view can be self-adapted to learn and then achieve a better clustering effect. Since the original brain MRI image is not MVD, this study extracts multiple feature data of the original image and uses one feature data as one view data to construct MVD. The experimental results prove that the segmentation method used in this paper optimizes the segmentation effect. With the increase in noise, the clustering effect of IMV-FCM decreases more slowly than the comparison method, which shows that the method used has stronger noise immunity. Although the effectiveness of IMV-FCM was verified, it still faces certain limitations. For example, the proposed method uses the classic FCM algorithm as a framework and uses Euclidean distance, which causes it to suffer from the dimensionality disaster problem when facing high-dimensional multiview clustering tasks. How to solve this problem will be the focus of future research. In addition, similar to most unsupervised learning algorithms, the selection of optimal parameters is an important issue. Different practical applications require different ranges of parameter values. Since the optimal parameters are usually determined by the application, for an exact application, the appropriate parameter range of the algorithm used can be determined through prior knowledge or available valid labeled datasets. The ensemble learning strategy can avoid the selection of optimal parameters to a certain extent, and it is planned to be further studied in follow-up work.

Data Availability Statement

Publicly available datasets were analyzed in this study. This data can be found here: https://brainweb.bic.mni.mcgill.ca.

Author Contributions

LH and YG conceived and developed the theoretical framework of the manuscript. YG and XG designed and performed the surgical procedure of the case presentation. LH, JX, and TN performed the data evaluation, analysis, and interpretation, designed the figures, designed and wrote the manuscript in consultation with YG and XG, who took the lead in writing the manuscript. All authors participated in the editing process.

Funding

This work was supported in part by the National Natural Science Foundation of China under Grants 81701793, 61702225, 61772241, 61873321, and61806026; in part by the Natural Science Foundation of Jiangsu Province under Grant BK20180956; in part by the 2018 Six Talent Peaks Project of Jiangsu Province under Grant XYDXX-127; in part by the Science and Technology Demonstration Project of Social Development of Wuxi under Grant WX18IVJN002; and in part by the Youth Foundation of the Commission of Health and Family Planning of Wuxi under Grant Q201654.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

Ahmed, M. N., Yamany, S. M., Mohamed, N., Farag, A. A., and Moriarty, T. (2002). Amodified fuzzy c-means algorithm for bias field estimation and segmentation of MRI data. IEEE Trans. Med. Imaging 21, 193–199. doi: 10.1109/42.996338

CrossRef Full Text | Google Scholar

Bezdek, J. C. (1980). A convergence theorem for the fuzzy ISODATA clustering algorithm. IEEE Trans. Pattern Anal. Mach. Intell. 2, 1–8. doi: 10.1109/TPAMI.1980.4766964

PubMed Abstract | CrossRef Full Text | Google Scholar

Bezdek, J. C. (1981). Pattern Recognition With Fuzzy Objective Function Algorithms. New York, NY: Plenum Press.

Google Scholar

Cai, T. T., Ma, J., and Zhang, L. (2019). Chime: clustering of high-dimensional Gaussian mixtures with EM algorithm and its optimality. Ann. Stat. 47, 1234–1267. doi: 10.1214/18-AOS1711

CrossRef Full Text | Google Scholar

Chen, S., and Zhang, D. (2004). Robust image segmentation using FCM with spatialconstraints based on new kernel-induced distance measure. IEEE Trans. Syst. Man Cybern. B Cyber. 34, 1907–1916. doi: 10.1109/TSMCB.2004.831165

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, X., Xu, X., Huang, J. Z., and Ye, Y. (2013). TW-k-means: automated two-level variable weighting clustering algorithm for multiview data. IEEE Trans. Knowl. Data Eng. 25, 932–944. doi: 10.1109/TKDE.2011.262

CrossRef Full Text | Google Scholar

Chrobak, M., Dürr, C., Aleksander, F., and Nilsson, B. J. (2020). Online clique clustering. Algorithmica 82, 938–965. doi: 10.1007/s00453-019-00625-1

CrossRef Full Text | Google Scholar

Chuang, K. S., Tzeng, H. L., Chen, S., Wu, J., and Chen, T. J. (2006). Fuzzy c-means clustering with spatial information for image segmentation. Comput. Med. Imaging Graph. Official J. Comput. Med. Imaging Soc. 30, 9–15. doi: 10.1016/j.compmedimag.2005.10.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Cleuziou, G., Exbrayat, M., Martin, L., and Sublemontier, J. H. (2009). “Co FKM: a centralized method for multiple-view clustering,” in Proceedings of the 9th IEEE Internaional Conference on Data Mining, (Miami FL: IEEE), 752–757. doi: 10.1109/ICDM.2009.138

CrossRef Full Text | Google Scholar

Devi, B. R., and Setty, S. P. (2018). Hybrid clustering algorithm ‘KCu’ for combining the features of K-means and CURE Algorithm for efficient outliers handling. Adv. Model. Anal. B 61, 76–79. doi: 10.18280/ama_b.610204

CrossRef Full Text | Google Scholar

Elazab, A., Wang, C., Jia, F., Wu, J., Li, G., and Hu, Q. (2015). Segmentation of brain tissues from magnetic resonance images using adaptively regularized kernel-based fuzzy c-means clustering. Comput. Math. Methods Med. 2015:485495. doi: 10.1155/2015/485495

PubMed Abstract | CrossRef Full Text | Google Scholar

Gong, M., Liang, Y., Shi, J., Ma, W., and Ma, J. (2013). Fuzzy c-means clustering with local information and kernel metric for image segmentation. IEEE Trans. Image Process. 22, 573–584. doi: 10.1109/TIP.2012.2219547

PubMed Abstract | CrossRef Full Text | Google Scholar

Gu, J., Cheng, T., Hua, L., Wang, J., Zhao, J., and Cao, Y. (2019). Overview of image segmentation and registration for spine biological modeling. J. Syst. Simul. 31, 167–173. doi: 10.16182/j.issn1004731x.joss.18-0806

CrossRef Full Text | Google Scholar

Gu, Q., and Zhou, J. (2009). Learning the shared subspace for multi-task clustering and transductive transfer classification. Proceedings of the 2009 9th IEEE International Conference on Data Mining, (Miami Beach, FL: IEEE), 159–168. doi: 10.1109/ICDM.2009.32

CrossRef Full Text | Google Scholar

Hall, L. O., Bensaid, A. M., and Clarke, L. P. (1992). A comparison of neural network and fuzzy clustering techniques in segmenting magnetic resonance images of the brain. IEEE Trans. Neural Netw. 3, 672–681. doi: 10.1109/72.159057

CrossRef Full Text | Google Scholar

Hall, L. O., and Goldgof, D. B. (2011). Convergence of the single-pass and online fuzzy C-means algorithms. IEEE Trans Fuzzy Syst. 19, 792–794. doi: 10.1109/TFUZZ.2011.2143418

CrossRef Full Text | Google Scholar

Ji, Z. X., Sun, Q. S., and Xia, D. S. (2011). A modified possibilisticfuzzy c-means clustering algorithm for bias field estimation and segmentation of brain MR image. Comput. Med. Imaging Graph. 35, 383–397. doi: 10.1016/j.compmedimag.2010.12.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, C. F., Chang, C. C., and Huang, S. H. (2012). Regions of interest extraction fromspect images for neural degeneration assessment using multimodality image fusion. Multidimen. Syst. Signal Process. 23, 437–449. doi: 10.1007/s11045-011-0162-3

CrossRef Full Text | Google Scholar

Jiang, Y., Bi, A., Xia, K., Xue, J., and Qian, P. (2020). Exemplar-based data stream clustering toward internet of things. J. Supercomput. 76, 2929–2957. doi: 10.1007/s11227-019-03080-5

CrossRef Full Text | Google Scholar

Jiang, Y., Chung, F. L., Wang, S., Deng, Z., Wang, J., and Qian, P. (2014). Collaborative fuzzy clustering from multiple weighted views. IEEE Trans. Cybern. 45, 688–701. doi: 10.1109/TCYB.2014.2334595

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, Y., Zhao, K., Xia, K., Xue, J., Zhou, L., Ding, Y., et al. (2019). A novel distributed multitask fuzzy clustering algorithm for automatic MR brain image segmentation. J. Med. Syst. 43:118. doi: 10.1007/s10916-019-1245-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Jing, L. P., Ng, M. K., and Huang, J. Z. (2007). An entropy weighting k-means algorithm for subspace clustering of high-dimensional sparse data. IEEE Trans. Knowl. Data Eng. 19, 1026–1041. doi: 10.1109/TKDE.2007.1048

CrossRef Full Text | Google Scholar

Krinidis, S., and Chatzis, V. (2010). A robust fuzzy local information C-means clustering algorithm. IEEE Trans. Image Process. 19, 1328–1337. doi: 10.1109/TIP.2010.2040763

PubMed Abstract | CrossRef Full Text | Google Scholar

Kumar, S. N., Lenin, F. A., and Sebastin, V. P. (2020). An overview of segmentation algorithms for the analysis of anomalies on medical images. J. Intell. Syst. 29, 612–625. doi: 10.1515/jisys-2017-0629

CrossRef Full Text | Google Scholar

Lee, K., Moon, C., and Nam, Y. (2018). Diagnosing vocal disorders using cobweb clustering of the jitter, shimmer, and harmonics-to-noise ratio. KSII Trans. Internet Inform. Syst. 12, 5541–5554. doi: 10.3837/tiis.2018.11.020

CrossRef Full Text | Google Scholar

Li, C., Gore, J. C., and Davatzikos, C. (2014). Multiplicative intrinsic component optimization(MICO) for MRI bias field estimation and tissue segmentation. Magn. Reson. Imaging 32, 413–439. doi: 10.1016/j.mri.2014.03.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, L., Kuang, L., and Ji, Y. (2020). Multimodal MRI brain tumor image segmentation using sparse subspace clustering algorithm. Comput. Math. Methods Med. 2020:8620403. doi: 10.1155/2020/8620403

PubMed Abstract | CrossRef Full Text | Google Scholar

Okita, Y., Miyata, A. H., Motomura, B. N., and Takamoto, S. (2015). A study of brain protection during total arch replacement comparing antegrade cerebral perfusion versus hypothermic circulatory arrest, with or without retrograde cerebral perfusion: analysis based on the Japan adult cardiovascular surgery database. J. Thorac. Cardiovasc. Surgery 149, 65–73. doi: 10.1016/j.jtcvs.2014.08.070

PubMed Abstract | CrossRef Full Text | Google Scholar

Pedrycz, W. (2002). Collaborative fuzzy clustering. Pattern Recogn. Lett. 23, 1675–1686.

Google Scholar

Qian, P., Jiang, Y., Deng, Z., Hu, L., Sun, S., Wang, S., et al. (2015). Cluster prototypes and fuzzy memberships jointly leveraged cross-domain maximum entropy clustering. IEEE Trans. Cybern. 46, 181–193. doi: 10.1109/TCYB.2015.2399351

PubMed Abstract | CrossRef Full Text | Google Scholar

Qian, P., Jiang, Y., Wang, S., Su, K. H., Wang, J., Hu, L., et al. (2016). Affinity and penalty jointly constrained spectral clustering with all-compatibility, flexibility, and robustness. IEEE Trans. Neural Netw. Learn. Syst. 28, 1123–1138. doi: 10.1109/TNNLS.2015.2511179

PubMed Abstract | CrossRef Full Text | Google Scholar

Qian, P., Zhao, K., Jiang, Y., Su, K. H., Deng, Z., Wang, S., et al. (2017). Knowledge-leveraged transfer fuzzy C-Means for texture image segmentation with self-adaptive cluster prototype matching. Knowl. Based Syst. 130, 33–50. doi: 10.1016/j.knosys.2017.05.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Qian, P., Zhou, J., Jiang, Y., Liang, F., and Muzic, R. F. Jr. (2018). Multi-view maximum entropy clustering by jointly leveraging inter-view collaborations and intra-view-weighted attributes. IEEE Access 6, 28594–28610. doi: 10.1109/ACCESS.2018.2825352

PubMed Abstract | CrossRef Full Text | Google Scholar

Rosati, P., Lupascu, C. A., and Tegolo, D. (2018). Analysis of low-correlated spatial gene expression patterns: a clustering approach in the mouse brain data hosted in the Allen Brain Atlas. IET Comput. Vision 12, 996–1006. doi: 10.1049/iet-cvi.2018.5217

CrossRef Full Text | Google Scholar

Singh, M., Vishal, V., Verma, A., and Sharma, N. (2020). Segmentation of MRI data using multi-objective antlion based improved fuzzy c-means. Biocybern. Biomed. Eng. 40, 1250–1266. doi: 10.1016/j.bbe.2020.07.001

CrossRef Full Text | Google Scholar

Sunjana Azizah, Z. (2020). Outlier detection of transaction data using DBSCAN algorithm. Int. J. Psychosoc. Rehabil. 24, 3232–3240. doi: 10.37200/IJPR/V24I2/PR200632

CrossRef Full Text | Google Scholar

Thanh, M. N., and Wu, Q. M. J. (2013). A fuzzy logic model based Markov random field for medical image segmentation. Evolv. Syst. 4, 171–181. doi: 10.1007/s12530-012-9066-1

CrossRef Full Text | Google Scholar

Wang, J., Schreiber, D. K., Bailey, N., Hosemann, P., and Toloczko, M. B. (2019). The application of the OPTICS algorithm to cluster analysis in atom probe tomography data. Microsc. Microanal. 25, 338–348. doi: 10.1017/S1431927618015386

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X. Y., Zhang, D. D., and Wei, N. (2015). Fractal image coding algorithm using particle swarm optimisation and hybrid quadtree partition scheme. Iet Image Process. 9, 153–161. doi: 10.1049/iet-ipr.2014.0001

CrossRef Full Text | Google Scholar

Weijer, J., and Gevers, T. (2006). Boosting color saliency in image feature. IEEE Trans. Pattern Anal. Mach. Intell. 28, 150–156. doi: 10.1109/TPAMI.2006.3

PubMed Abstract | CrossRef Full Text | Google Scholar

Yin, S., Li, H., Liu, D., and Karim, S. (2020). Active contour modal based on density-oriented BIRCH clustering method for medical image segmentation. Multimedia Tools Appl. 79, 31049–31068. doi: 10.1016/j.ijleo.2018.01.004

CrossRef Full Text | Google Scholar

Yu, S., Tranchevent, L., Liu, X., Glanzel, W., Suykens, J. A. K., De Moor, B., et al. (2012). Optimized data fusion for kernel k-means clustering. IEEE Trans. Pattern Anal. Mach. Intell. 34, 1031–1039. doi: 10.1109/TPAMI.2011.255

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, C., Churazov, E., and Zhuravleva, I. (2020). Pairs of giant shock waves (N-waves) in merging galaxy clusters. Monthly Notices R. Astronomic. Soc. 501, 1038–1045.

Google Scholar

Zheng, Y., Byeungwoo, J., Xu, D., Wu, Q. M. J., and Hui, Z. (2015). Image segmentation by generalized hierarchical fuzzy C- means algorithm. J. Intell. Fuzzy Syst. 28, 4024–4028. doi: 10.3233/IFS-141378

CrossRef Full Text | Google Scholar

Zhu, L., Chung, F. L., and Wang, S. T. (2009). Generalized fuzzy k-means clustering algorithm with improved fuzzy partitions. IEEE Trans Syst. Man Cybern. 39, 578–591. doi: 10.1109/TSMCB.2008.2004818

PubMed Abstract | CrossRef Full Text | Google Scholar