A Variational Level Set Model Combined with FCMS for Image Clustering Segmentation

The fuzzy C means clustering algorithm with spatial constraint (FCMS) is effective for image segmentation. However, it lacks essential smoothing constraints to the cluster boundaries and enough robustness to the noise. Samson et al. proposed a variational level set model for image clustering segmentation, which can get the smooth cluster boundaries and closed cluster regions due to the use of level set scheme. However it is very sensitive to the noise since it is actually a hard Cmeans clusteringmodel. In this paper, based on Samson’s work, we propose a new variational level set model combined with FCMS for image clustering segmentation. Compared with FCMS clustering, the proposed model can get smooth cluster boundaries and closed cluster regions due to the use of level set scheme. In addition, a block-based energy is incorporated into the energy functional, which enables the proposed model to bemore robust to the noise than FCMS clustering and Samson’smodel. Some experiments on the synthetic and real images are performed to assess the performance of the proposed model. Compared with some classical image segmentation models, the proposed model has a better performance for the images contaminated by different noise levels.


Introduction
Image segmentation is separating the image domain into dissimilar homogeneous regions, which is the precondition and foundation of further image analysis and understanding.The quality of segmentation affects the result of the following analysis and processing directly.So, it is an important technique in image processing and has drawn much research attention at the theory and application.In recent years, variational level set method and clustering technology have been widely exploited in image segmentation because of their good experimental performance and sound theoretical foundation.
The variational level set model used in image segmentation is formulated as follows [1][2][3][4].The contours are first implicitly represented by the zero level set (ZLS) of a higher dimensional function, usually referred to as the level set function.And then one can obtain evolution partial differential equation (PDE) or partial differential equations (PDEs) for level set function in terms of minimizing an energy functional which typically includes the internal energy that smoothes the level set function and the external energy that aligns the ZLS with object boundaries.At last, the level set function evolves according to the evolution PDE or PDEs and thus achieves the goal for evolving the ZLS implied therein.Compared with the traditional level set method where the level set function is driven purely by PDE, in the variational level set method, the evolution PDE is obtained by minimizing energy functional.So, more prior information (e.g., texture and shape information) on the image can be conveniently taken into account into the energy functional, which makes the variational level set method have a good performance and extensive adaptability.Here, we present some classical variational level set models for image segmentation.The well-known Mumford and Shah (MS) [5] model for image segmentation has been successfully extended to a wide range of applications.But it cannot be solved directly in practice because of the nonconvexity of its functional.Recently, some improved models are proposed, such as piecewise constant model [6] proposed by Chan and Vese (CV) and region-based active contour model [7] (RBACM) proposed by Zhang et al..In order to enhance the quality of segmentation for image with inhomogeneous intensity, Li et al. [8] proposed an implicit active contour driven by local binary fitting energy (LBF), and then Zhang et al. [9] proposed an active contour 2 Mathematical Problems in Engineering driven by local image fitting energy (LIF) which has higher computing efficiency than LBF.The models mentioned above can only be utilized for two-phase partition.Some multiphase models are proposed, such as Vese and Chan [10] proposed multiphase CV model (MCV), Gao and Yan [11] proposed multiphase local CV model (MLCV) to improve the efficiency for noisy image segmentation.
Clustering is to partition a given input dataset or image pixels into  clusters with most similarities in the same cluster and most dissimilarities between different clusters.In the last decades, the fuzzy C means clustering (FCM) [12] has been widely used in image segmentation (e.g., [13,14]) due to its good performance and a well-grounded theory.Such a success chiefly attributes to the introduction of fuzzy membership relations between the image pixels and the cluster centers,   ( = 1, . . ., ).This allows the ability of FCM to be able to retain more image information than the hard C means clustering.The original FCM clustering has a good performance on segmenting the most noisefree image, but it fails to segment images contaminated by noise, outliers, and other imaging artifacts.Ahmed et al. [15] first considered the fuzzy C means with spatial constraints (FCMS), that is, incorporating the spatial information of image into the objective function, to overcome this difficulty.Afterward, some FCMS-based models were proposed to meet different research requirements.Such as Chen and Zhang [16] modified the FCMS objective function to reduce the computational complexity.They then replaced the Euclidean distance in FCMS by a kernel-induced distance and then proposed a Gaussian kernel version of FCMS, called GFCMS later.Yang and Tsai [17] proposed a generalized type of GFCMS in which the parameters can be automatically estimated under a learning scheme.Kannan et al. [18] proposed an effective FCMS for segmenting medical images.Liu et al. [19] proposed a fuzzy spectral clustering combined with spatial information.He et al. [20] proposed a new FCM clustering with total variation regularization for segmenting the images with noisy and incomplete data.
The above mentioned clustering algorithms are all based on the discrete data.They utilize the intensity, statistics properties, and spatial features of image pixels to perform pixels clustering.But they cannot obtain the smooth cluster boundaries and closed cluster regions due to the lack of the essential smoothing constraint to the cluster boundaries, while the variational level set method just can deal with the above problems.So, the image segmentation quality will be improved if the clustering algorithms are appropriately combined with the variational level set method.However, most of the current clustering algorithms are based on the discrete dataset, while variation method is based on continuous function.So, these two techniques cannot be combined together easily.Samson et al. [21] first solved this problem in 2000 by the use of level set method.As the first attempt of combining data clustering with variation method, Samson's model still has some drawbacks, such as (1)  The remainder of this paper is organized as follows: in Section 2, some backgrounds concerning the standard FCMS clustering (and its variants) and Samson's model are presented; some drawbacks of them are also mentioned.In Section 3, a variational level set model combined with FCMS (i.e., VFCMS) is proposed.In Section 4, we apply the proposed model to image clustering segmentation.The comparisons with some classical image segmentation models are also performed in this section.This paper is summarized in Section 5.

FCMS Clustering and Its Variants.
In [15], the authors proposed FCMS clustering to partition the discrete dataset {x  }  =1 into  clusters.The main contribution of FCMS is that the spatial information of the discrete dataset was incorporated into the objective function, which can increase the robustness to noise.The objective function of FCMS is defined as where   is a fuzzy membership matrix and  is weighting exponent on each fuzzy membership; {k  }  =1 are the cluster centers.  is the set of neighbors falling into the window centered at x  and   is its cardinality; the parameter  is the weighting coefficient of the spatial constraints.In essence, the spatial constraint (the second term in (1)) aims at keeping the continuity on the neighboring data values around x  .By minimizing in a way similar to the standard FCM algorithm, the necessary conditions on   and k  for (1) to be at a local minimum are The procedure of FCMS clustering is as follows.First, based on the prior information of dataset, predefine the cluster number  and the initial cluster centers {k  }  =1 .And then update the fuzzy membership matrix   and cluster centers    {k  }  =1 successively by (2) until ‖k new − k old ‖ < .Then the output k new = (V 1 , . . ., V  ) is the final clustering centers.
In [16], Chen and Zhang studied the FCMS clustering and pointed out a shortcoming of its update equations (2), that is, computing the neighborhood terms will take much more time than the classical FCM.They noticed that where x  is the median or mean of neighboring data lying within a window around x  , and then proposed a modified FCMS objective function by replacing     Similarly, minimizing the J FCMS (   , k  ) in (4) with respect to   and k  , we can obtain where x  can be computed in advance.Obviously, updating (5) is simpler than (2).So the clustering time can be saved.For convenience of notation, the authors of [16] named the clustering algorithm using (5) with median and mean filtering FCMS 1 and FCMS 2 , respectively.Although FCMS clustering and its variants (e.g., FCMS 1 and FCMS 2 ) have the benefits that it is simple and easy to manipulate, it cannot obtain the smooth cluster boundaries and the closed cluster regions for the lack of the essential smoothing constraints for the cluster boundaries.In addition, although the spatial constraint is incorporated into the objective function, FCMS cannot achieve good clustering result when the dataset is contaminated by strong noise.

Samson's Model.
The mentioned above clustering algorithms are all based on the discrete data, so they cannot be solved by the variational method directly.It was the first time for Samson et al. to employ the variational method for data clustering by using level set method.In [21], they proposed a clustering model based on variational level set and then applied it to image clustering segmentation.Let image domain be Ω ∈ R 2 and let image function be (x) : Ω → R.
Image clustering segmentation is equivalent to solving the following minimization problem: inf The first term of the energy functional ( 6) aims to partition the image domain Ω into  subregions where the image intensity has a Gaussian distribution of mean   and of standard deviation   .The second term is the regularization energy; minimizing it is equivalent to minimizing the interface between clusters.Minimizing the third term leads to a solution where the formation of a vacuum (pixel with no labels) and overlapping (pixel with more than one labels) regions are penalized.
As the first attempt of data clustering manipulated by the use of variational method, Samson's model has an advantage over the traditional clustering algorithms in obtaining the smooth cluster boundaries and the closed cluster regions.However, there are still some drawbacks as follows.
In this paper, following Samson's work, we propose a variational level set model combined with FCMS clustering (VFCMS) for image clustering segmentation.Four schemes are introduced to resolve the above mentioned drawbacks of Samson's model (6).
(1) A block-based clustering energy and a spatial constraint are introduced into the energy functional.In addition, the fuzziness of belongingness of each pixel to the cluster centers is introduced.These improvements enable the new model to be more robust to noises than Samson's model and FCMS clustering.(2) A variational formulation is proposed for updating membership functions and cluster centers, which makes the new model more robust to the initial cluster centers and achieves a semisupervised clustering.(3) A regularization term based on a new edge stopping function is proposed, which enables the active contours to move quickly through the noise regions and reach the right boundaries of image.(4) A regularization term is introduced to eliminate the need of the costly reinitialization procedure.[21] is defined as

The Proposed VFCMS Model
Unfortunately, minimizing directly this energy functional to achieve the image clustering has two drawbacks: (1) it is sensitive to noises since it is point-based and (2) the ratio   /  of weighted parameter   to the standard deviation   can be seen as the membership between the image pixel (x) and the cluster center   .In [21], these two parameters are both chosen as constants.So, the clustering algorithm proposed by Based on the points discussed above, we introduce the following fuzzy and block-based clustering energy: where   (x) is a point spread function, (e.g., Gaussian function with standard deviation ).The constraint ∑  =1 (  (x)) = 1 is to penalize the overlapping and vacuum formation of the clustering regions and In what follows, we analyze this energy in the theory.Denote The point spread function   (x) satisfies that ∫   (x)x = 1.
Thus we have   in which From the last equation, we have If we ignore this term, the energy functional ( 9) can be rewritten as Since f(x) can be seen as a denoised image with a smoothing kernel   (x), minimizing energy ( 14) is equivalent to cluster the denoised image by FCM clustering.Thus, the proposed model is more robust to noise than the standard FCM clustering and Samson's model.In order to further increase the robustness to noises, we introduce the following spatial constraint into the energy functional.Note that here we adopt a modified form proposed by Chen and Zhang in [16] to save the computing time.
where  > 0 is a tuning parameter and (x) can be can be considered to be the mean or median of (x) supported on the disk centered at x. Similar to (4), the value of (x) can be computed in advance, thus the clustering time can be saved.Combining ( 9) and ( 15), the total external clustering energy is Similar to [16], for convenience of notation later, the proposed VFCMS model is renamed as VFCMS 1 and VFCMS 2 corresponding to (x) being median filtering and mean filtering, respectively.In addition, in what follows, we always write 3.1.2.The Optimal Membership Functions.Fixing level set functions {  }  =1 and clustering centers {  }  =1 , we seek the optimal membership functions   (x) which make the energy functional E  ( 1 ⋅ ⋅ ⋅   ,   ,   ) to converge to a local minimum.In order to improve the ability of active contour to capture the pixels belonging to its cluster, we compute the optimal membership functions that are supported on the whole image domain Ω, that is, computing the membership relations between each pixel x and each cluster center   .Firstly, we extend the fuzzy clustering energy E  ( 1 ⋅ ⋅ ⋅   ,   ,   ) to the whole image domain Ω, denoted as  respect to each   (x), we can obtain the optimal membership functions supported on total image domain Ω.The detail is stated as follows: min Using calculus of variation and Lagrange multiplier method, the necessary condition on   (x) for ( 18) to be at a local minimum is Solving each   (x) in the last equations, we can obtain the optimal membership functions In the experiments, (20) gives us the updating formula of the optimal membership functions.

The Optimal Cluster Centers. Fixing level set functions
{  }  =1 and membership functions {  }  =1 , we seek the optimal cluster centers {  }  =1 which make the energy functional E  ( 1 ⋅ ⋅ ⋅   ,   ,   ) to be a local minimum.That is, minimizing the energy Solving   in the last equation, we have The optimal cluster center   is actually the weighted mean of (x) supported on the Ω  .Equation (23) gives us the updating formula of the optimal cluster centers in the iterative process.

The Internal Energy.
In this section, we introduce two internal energies E  1 and E  2 , where E  1 is to penalize the singularities of level set functions and E  2 is to penalize the formation of vacuum and overlapping regions.
In the traditional variational level set method for image processing, in order to maintain the stability of the level set function during the evolution, the evolving level set function needs periodical reinitialization to keep it close to a signed distance function [22].Samson et al. achieved the reinitialization by periodically solving PDEs (7).Many serious problems remain such as when to apply the reinitialization and the computational complexity increases.In this paper, we adopt method proposed by Zhang et al. in [23] and introduce the following internal energy: to eliminate the need of the expensive reinitialization procedure.Here,  > 0 is a tuning parameter.The energy E  1 can be identified as a metric to measure the smoothness of the level set function.
The constraint ∑  =1 (  (x)) = 1 is to penalize the vacuum and overlapping regions.In this paper, similar to (6), we introduce the following internal energy to meet this constraint: where  > 0 is a tuning parameter.The value of ∑  =1 (  (x)) will trend to 1 in the process of minimizing the energy E  2 .So, the overlapping and vacuum cluster regions will decrease in  the process of clustering.In what follows, we write the total of internal energy as 3.3.The Regularization Energy.Samson et al. [21] introduced the following regularization energy to smooth the boundaries of the clusters: where the stopping function is defined as which is a decreasing function of the gradient module of image.The evolving velocity of ZLS is about one (((x)) → 1) at the smooth position of the image, since the gradient module is about zero (|∇(x)| → 0) at this position, which makes the ZLS to move quickly through the smooth position.At the position of the edge, the gradient module |∇(x)| → ∞ and the evolving velocity ((x)) → 0, which makes the ZLS to stay the edges.But if the data is very noisy, the evolving velocity of ZLS is about zero at the position of isolated noise, since the gradient module is about infinite at this position, which makes ZLS stay at the position of isolated noise and results in failed image segmentation.So, the traditional edge indicators based on the image gradient module cannot effectively distinguish between edges and isolated noises.In [24], the authors presented a new edge indicator based on the second derivatives, which is defined as In experiments, we choose   (x) Gaussian function with standard deviation , that is, The functional E still has a drawback from the practical point of view; that is, E is not Gateaux differentiable.So we have to regularize E. To do this, we use the following regularization Heaviside and Dirac function defined as to approximate the standard Heaviside and Dirac functions, respectively.In this paper, we choose that  = 5.

Algorithm (VFCMS Clustering)
Step 1.Given the number of the classes  and the initial level set functions {  }  =1 , choosing the initial cluster centers {  }  =1 .
Step 5.If max  |  −   | ≤ , end the algorithm; else let   =   and go back to Step 2.

Experimental Results
In this section, we show the experimental results of image segmentation on several synthetic and real images.There are a total of twelve models used in this section, that is, (1) clustering models (FCM, FCMS 1 , FCMS 2, and GFCMS), (2) variational level set models (CV, IVC, LIF, MCV, and MLCV), and (3) the integration of variational level set and clustering (Samson's model, VFCMS 1, and VFCMS 2 ).

The Choice of Parameters.
For the choice of the initial level set function, under the combined effects of the clustering energy and internal energy, the choice of the initial level set function is very flexible.In our experiments, the initial level set functions are all chosen as   (x) = 1 ( = 1, . . ., ).
We also tested our model by using other initial level set functions, such as   (x) = 0, signed distance function, and piecewise constant function.We found that our model using these different initial level set functions can all get correct image segmentation.We adopt a simple method to choose the initial cluster centers, which is stated as follows.Let  = min x∈Ω (x),  = max x∈Ω (x), and ℎ = ( − )/; then where Ω  = {x ∈ Ω :  + ( − 1)ℎ ≤ (x) <  + ℎ}.
The parameter  is a weighting parameter to measure the penalization for the formation of a vacuum and regions overlapping, which has a very important effect on the performance of the proposed model.So the value of  should not be very small.In all experiments, we set  = 50.We also find that for a wide range of  over 50 (e.g., from 50 to 100), there seem to be no apparent changes in the results.The parameter  adjusts the degree of smoothing.If it is too small, we cannot obtain good segmentation results due to the effect of the noises.Conversely, if  is too large, the clustering boundaries will deviate from the image boundaries due to the oversmoothing of the image.How to choose an optimal smoothing parameter  is still an "open question." In this paper, we adopted the "trial and error" technique to determine the value of smoothing parameter .In the experiments, if the image is corrupted by mixed 0.5% or 1% Salt and Pepper, Gaussian, and Speckle noise, we set  = 2 and if the image is corrupted by 1.5% or 2% mixed noise, we set  = 2.In order to obtain the optimal segmentation under the case of  being fixed, the choice of parameter  is very important in the proposed model.For the FCMS 1 , FCMS 2 , GFCMS, VFCMS 1, and VFCMS 2 clustering models, the parameter  controls the weight of the spatial constraints.In the experiments, we still adopt "trial and error" technique to determine the value of parameter .The detail of the choice of  is presented in Section 4.2.For the other parameters, we set  = 0.1 and  = 1.

Image Clustering Segmentation
Example 1 (the clustering segmentation of noisy synthetic images).Figure 1 shows the segmentation results of a synthetic image corrupted by 1% Salt and Pepper, Gaussian, and Speckle noise simultaneously.This image with 128 × 128 pixels contains two clusters (i.e., object and background).Experiment results show that FCM clustering, CV model, LIF model, and Samson's model cannot segment image well due to the effect of the noise.Since RBACM model utilizes the global image intensities inside and outside the contour, it can get better segmentation result.FCMS 1 , FCMS 2, and GFCMS model also can get good segmentation due to the introduction of spatial information.However, we find that FCMS 2 and GFCMS clustering have fewer pixels being misclassified than FCMS 1 .This is because FCMS 1 is implemented by median filtering, while FCMS 2 is implemented by mean filtering and GFCMS replaces the European distance with Gaussian kernel-induced distance.From Figures 1(n)-1(q), we can see that the proposed VFCMS 1 and VFCMS 2 models both can give us good segmentation.This is because we introduce the block-based term in the functional, except for spatial constraints.The last two plots show the surface plot for the convergent level set functions of VFCMS 2 model, respectively.We take a set of values for  to test its performance in FCMS 1 , FCMS 2 , GFCMS, VFCMS 1, and VFCMS 2 clustering models.Figure 2(a) shows the comparisons of classification errors of these models under different values of  on the synthetic two-phase image shown in Figure 1(a).From Figure 2(a), as  increases, the numbers of misclassified pixels of all five models firstly reduce.FCMS 1 , FCMS 2 , and GFCMS reach minima between  = 2.8 and  = 3.2, and VFCMS 1 and VFCMS 2 reach minima when  in the interval of [1.4,3.0].From Figure 2(a), we note that as  continues to increase, the numbers of misclassified pixels of all five models will increase.In the first experiment (i.e., Figure 1), we set  = 2.2 in VFCMS 1 and VFCMS 2 , and set  = 2.8 in FCMS 1 , FCMS 2, and GFCMS.
Figure 3 shows segmentation results of the synthetic two-phase image corrupted by 2% Salt and Pepper, Gaussian, and Speckle noise simultaneously.For the choice of , similarly, we take a set of values for  to test its performance.Figure 2(b) shows the comparisons of classification errors of FCMS 1 , FCMS 2 , GFCMS, VFCMS 1, and VFCMS 2 clustering models under different values of .We obtain similar conclusions to the first experiment (see Figure 2(a)).From Figure 2(b), we set  = 3 in VFCMS 1 and VFCMS 2 , and set  = 3.2 in FCMS 1 , FCMS 2, and GFCMS.From Figure 3, we can clearly see that compared with the other models, the proposed models can achieve obvious predominance when image is corrupted by strong noise.
Figure 4 shows the segmentation results of a synthetic three-phase image corrupted by 0.5% Salt and Pepper, Gaussian, and Speckle noise.This image with 128 × 128 pixels contains three clusters.In the experiments, the parameters  Figure 5 shows the segmentation results of the synthetic three-phase image corrupted by 1.5% Salt and Pepper, Gaussian, and Speckle noise.We set  = 2.6 in VFCMS 1 and VFCMS 2 , and  = 3 in FCMS 1 , FCMS 2, and GFCMS.From experiment results, we can see that FCM clustering, CV model, LIF model, and Samson's model still do not segment image well.It is observed that FCMS 1 also cannot achieve satisfactory result in the case of image being corrupted by strong noise.This is because it is implemented by median filtering.Similar to Figure 4, FCMS 2 , GFCMS, MCV, MLCV, VFCMS 1, and VFCMS 2 can achieve better segmentation.
Table 1 gives us the segmentation accuracy (SA) of the models on the images shown in Figures 1(a), 3(a), 4(a), and 5(a), where SA is defined as the total number of correctly classified pixels by the total number of all pixels [16].By quantitative comparison of SA, the proposed VFCMS 1 and VFCMS 2 can achieve more accurate segmentation than the other models under different noise levels.
Example 2 (the clustering segmentation of noisy plane image).Figure 6 presents the comparison results (clustering into 2 clusters) on a plane image corrupted by 2% Salt and Pepper, Gaussian, and Speckle noise simultaneously.It is obvious that FCM, CV, LIF, and Samson's model cannot segment image well.FCMS 1 , FCMS 2 , GFCMS, and GBACM can achieve better segmentation than FCM, CV, LIF, and Samson's model.Among these ten models, the proposed VFCMS 1 and VFCMS 2 can obtain the superior segmentation, and only a few pixels are misclassified.Furthermore, we find that VFCMS 2 can give us a better segmentation than VFCMS 1 , This is because VFCMS 1 is implemented by median filtering, and VFCMS 2 is implemented by mean filtering.We conclude that if the test data is corrupted by strong noise, VFCMS 2 can get better segmentation than VFCMS 1 in general.The last two plots show the surface plot of the convergent level set functions of VFCMS 2 .
Figure 7 presents the comparison results (clustering into 3 clusters) on the noisy plane image corrupted by 2% mixed Salt and Pepper, Gaussian, and Speckle noises.This image contains two low-contrast regions (i.e., background region and the shadow region of plane).Furthermore, these two regions are contaminated by strong noises.We would like to segment noisy plane image into 3 parts and discriminate

Conclusions
In this paper, based on the Samson's work and FCMS clustering, we proposed a new variational level set model combined with FCMS for image clustering segmentation.In addition, a block-based energy was incorporated into the energy functional, which enables the proposed models robust to the noise.Some synthetic and real noisy images with different noise levels were employed to compare the performance of 12 models.Experimental results show that the proposed model has a superior performance among these methods.And different kinds of real noisy image were also used to show the practicability and validity of the proposed model.The experimental results reported in this paper show that the proposed VFCMS model is very effective for noise image clustering segmentation.This model can also be improved by incorporating other FCMS-based clustering algorithms, for example, Kernel-induced FCMS proposed by Chen and Zhang [16] or Gaussian Kernel-induced FCMS proposed by Yang and Tsai [17], and so forth.We note that in passing our model still has some drawbacks, such as (1) it is a semisupervised image clustering segmentation and needs to predetermine the clustering numbers; (2) Key parameters such as  weighting spatial term and  controlling Gaussian smoothing are determined by human-machine interaction; (3) the proposed model is more complex and takes more computing time than some other classical variational level set model, such as CV, LBF, LIF, and MCV.So, the proposed model has major drawback of weak real-time performance.Our further works will include autoselecting parameters, adaptive determination of the clustering number, and improving real-time capability of VFCMS.

Figure 1 :
Figure 1: Comparison of segmentation results on a synthetic image with mixed 1% Salt and Pepper, Gaussian, and Speckle noise (cluster number  = 2).

Figure 2 :
Figure 2: Comparison of classification errors on two-phase synthetic image with different mixed noise levels under different values of alpha.

Figure 3 :
Figure 3: Comparison of segmentation results on a synthetic image with mixed 2% Salt and Pepper, Gaussian, and Speckle noise (cluster number  = 2).
The surface plot of  3

Figure 4 :
Figure 4: Comparison of segmentation results on a synthetic image with mixed 0.5% Salt and Pepper, Gaussian, and Speckle noise (cluster number  = 3).

( 1 )
Samson's model is actually a hard C means clustering which lacks the ability to retain abundant information from the original image and is also very sensitive to noise.(2) The external clustering energy (the first term of ( 1 ⋅ ⋅ ⋅   )) is point-based, which makes the clustering results sensitive to noise and outliers.(3) It is a supervised image classification model.The cluster number  and the cluster centers {  }  =1 must be given by the preclustering.In addition, the updating schemes for the cluster centers are not introduced by Samson et al.So, the clustering result is very sensitive to the choice of the initial cluster centers.(4) The evolving level set function   (x, ) needs periodical reinitialization to keep it close to a signed distance function during its evolution by solving the following PDE:   (x, )  = sign (  (x, )) (1 −     ∇  (x, )     ) .

Figure 6 :
Figure 6: Comparison of segmentation results on a plane image with mixed 2% Salt and Pepper, Gaussian, and Speckle noise (cluster number  = 2).

Figure 7 :
Figure 7: Comparison of segmentation results on a plane image with mixed 2% Salt and Pepper, Gaussian, and Speckle noise (cluster number  = 3).

Figure 8 :
Figure 8: Segmentation results on a sun image with mixed 2% Salt and Pepper, Gaussian, and Speckle noise (cluster number  = 2).

Figure 9 :
Figure 9: Segmentation results on a palm image with mixed 2% Salt and Pepper, Gaussian, and Speckle noise (cluster number  = 2).

Figure 10 :
Figure 10: Segmentation results on a light image with mixed 2% Salt and Pepper, Gaussian, and Speckle noise (cluster number  = 2).

Figure 11 :
Figure 11: Segmentation results on satellite image with mixed 2% Salt and Pepper, Gaussian, and Speckle noise (cluster number  = 3).

Figure 12 :
Figure 12: Segmentation results on a butterfly image with mixed 2% Salt and Pepper, Gaussian, and Speckle noise (cluster number  = 3).

3. 4 .
The Numerical Implementation.Combined with external energy E  , internal energy E  and regularization energy E  , image clustering segmentation is equivalent to minimize the following energy functional: (a) Original image and its noisy version (b) VFCMS 1 (c) VFCMS 2

Figure 13 :
Figure 13: Segmentation results on a panda image with mixed 2% Salt and Pepper, Gaussian, and Speckle noise (cluster number  = 3).

Figure 14 :
Figure 14: Segmentation results on a CT image with mixed 2% Salt and Pepper, Gaussian, and Speckle noise (cluster number  = 4).

Figure 15 :
Figure 15: Segmentation results on a CT image with mixed 2% Salt and Pepper, Gaussian, and Speckle noise (cluster number  = 4).

Figure 16 :
Figure 16: Segmentation results on a CT image with mixed 2% Salt and Pepper, Gaussian, and Speckle noise (cluster number  = 4).
it is a hard C means clustering; (2) it is very sensitive to the noise; and (3) it is a supervised clustering model.To solve these drawbacks, we propose a variational level set model combined with FCMS clustering (called VFCMS later) in this paper.Compared with Samson's model, the proposed model has the following advantages: (1) it is a fuzzy clustering model; (2) it is very robust to the noise; and (3) it is a semisupervised clustering model.
and   represent the second directional derivatives in the direction of the gradient ∇ and in the perpendicular direction of ∇, respectively.| ⋅ | denotes the absolute value.The performance of the new edge indicator is as follows: (1) for the edges, |  | is large and |  | is small, so  is large and (2) for the isolated noises, |  | and |  | are both large and almost equal, so  is small.According to the analysis mentioned above, edges and isolated noises can be well distinguished based on the value of .
which is based on the edge indicator (29).Mathematical Problems in Engineering

Table 1 :
SA % of twe12e models on noisy synthetic image.-contrastregions.From this figure, we can see that FCM, FCMS 1 , FCMS 2 , GFCMS, MCV, MLCV, and Samson's model cannot distinguish these two low-contrast regions well.While the proposed VFCMS 1 and VFCMS 2 model can obtain the good segmentation results on these two low-contrast regions due to the use of block-based energy and variational level set scheme.The last three plots show the surface plot of the convergent level set functions of VFCMS 2 .Example 3 (the clustering segmentation of real images corrupted by mixed noise).Finally, to show the practicability and validity of the proposed model, different kinds of real images corrupted by mixed Salt and Pepper, Gaussian, and Speckle noise are tested.Compared with the synthetic images, these real images contain much more complex boundaries, weak boundaries, inhomogeneous regions, and low-contrast regions.Figures8, 9, and 10 show the segmentation results of the images which are segmented into two clusters.Figures 11,12, and 13 show the segmentation results of the images which are segmented into three clusters.The segmentation results of the images that are segmented into four clusters are shown in Figures 14, 15, and 16.From Figures 8-16, we can clearly see that the proposed model can overcome the influence of noise and obtain the excellent segmentation results for different kinds of noisy real images.