An Improved FCM Medical Image Segmentation Algorithm Based on MMTD

Image segmentation plays an important role in medical image processing. Fuzzy c-means (FCM) is one of the popular clustering algorithms for medical image segmentation. But FCM is highly vulnerable to noise due to not considering the spatial information in image segmentation. This paper introduces medium mathematics system which is employed to process fuzzy information for image segmentation. It establishes the medium similarity measure based on the measure of medium truth degree (MMTD) and uses the correlation of the pixel and its neighbors to define the medium membership function. An improved FCM medical image segmentation algorithm based on MMTD which takes some spatial features into account is proposed in this paper. The experimental results show that the proposed algorithm is more antinoise than the standard FCM, with more certainty and less fuzziness. This will lead to its practicable and effective applications in medical image segmentation.


Introduction
Image segmentation is the procedure in which the original image is partitioned into homogeneous regions and plays an important role in medical image processing. As the imaging mechanism and the tissues of medical images are different, medical images are easily affected by noise, field migration effect, and tissue movement. Compared with the common images, medical images have more nonuniformity and fuzziness. Several popular clustering techniques for segmentation are available. Fuzzy c-means (FCM) is one such soft segmentation technique applicable for medical images. The performance of this method to obtain an optimal solution depends on the initial positions of the centers of the clusters, the measure of membership degree for each data point, and so on. In the standard FCM, the centers are initialized randomly and the measure of membership only uses the gray feature. This leads to be quite time-consuming and be sensitive to noise. For years, many research efforts have been made towards effective FCM image segmentation approaches. In order to accelerate the segmentation process, some approaches [1][2][3][4] focus on how to initialize the centers of required clusters. As noise is always emerged in the medical images, some FCM segmentation approaches [5,6] which have less sensitivity to noise are also presented. Moreover, in order to overcome the defect that the FCM is easy to fall into local optimal solution, some scholars combine the FCM and other mathematical approaches [7,8]. Nevertheless, due to their complexities and poor portabilities, effective medical image segmentation methods are yet to be seen. This paper introduces medium mathematics system which is employed to process fuzzy information for image segmentation. It establishes the medium similarity measure based on the measure of medium truth degree (MMTD) and uses the correlation of the pixel and its neighbors to define the medium membership function. An improved FCM medical image segmentation algorithm based on MMTD which takes some spatial features into account is proposed in this paper. The experimental results show that the proposed algorithm is more antinoise than the standard FCM, with more certainty and less fuzziness.

Computational and Mathematical Methods in Medicine
╕ ╕P + P ∼P Figure 1: Relation between numerical value areas and predication.
x(i, j) medium mathematics system [10], the corner stone of which is medium axiomatic sets.

Basic Symbols in Medium Mathematics
System. In medium mathematics system [10], predicate (conception or quality) is represented by and any variable is denoted as , with completely possessing quality being described as ( ). The "j" symbol stands for inverse opposite negative and it is termed as "opposite to. " The inverse opposite of predicate is denoted as j and then the concept of a pair of inverse opposite is represented by both and j . Symbol "∼" denotes fuzzy negative which reflects the medium state of "either or" or "both this and that" in opposite transition process. The fuzzy negative profoundly reflects fuzziness; "≺" is a truth-value degree connective which describes the difference between two propositions.

Measure of Medium Truth
Degree. According to the concept of super state [11,12], the numerical value area of generally applicable quantification is divided into five areas corresponding to the predicate truth scale, namely, j + , j , ∼ , , and + , as shown in Figure 1. In "True" numerical value area T, is standard scale of predicate ; in "False" numerical value area F, is standard scale of predicate j .
Individual truth scale in each numerical value area can be calculated by the distance ratio ℎ ( ( )) (or ℎ ( ( ))) [11] which relates to (or j ).

The MMTDFCM Medical Image Segmentation Algorithm
3.1. Medium Similarity Measure of Two Pixels. Similarity degree between two pixels can be scaled by the difference or the ratio of the gray level between pixels. This paper firstly discusses a new similarity measure based on MMTD [11,12]. Given there are two pixels ( , ) and ( , ) in grey images whose range of gray level is 0 ∼ 255. The two pixels may be in the same image but have different locations or they belong to different images but have the same coordinate. First, we relate the gray level of pixels to a number axis. Second, predicate ( ( , ), ( , )) represents that ( , ) is similar to ( , ), j ( ( , ), ( , )) represents that ( , ) is different from ( , ) and ∼ ( ( , ), ( , )) transition, as shown in Figure 2. According to the measure of medium truth degree, we give the definition of the similarity degree between ( , ) and ( , ). where as the medium similarity measure of two pixels, where ( , ) is the Euclidean distance between and . Under the onedimensional circumstance, ( , ) = | − |. The larger the value is, the higher the similarity degree between ( , ) and ( , ) becomes. When it equals 1, it shows that the pixels ( , ) and ( , ) are similar. The smaller the value is, the lower the similarity degree between ( , ) and ( , ) becomes. When it equals 0, it shows that the pixels ( , ) and ( , ) are different.

Medium Membership Function.
There are abundant information reflected by gray, edge, texture, and space in an image. Images are always polluted during image acquisition and transmission, so the noise should be removed before the further image segmentation. Standard FCM image segmentation approach which only takes the gray feature into account and ignores the other features is very sensitive to noise. This leads to some wrong classifications. For example, some pixels which should belong to the homogeneous region are separated. In order to improve the antinoise and the effect of the segmentation, we introduce the spatial feature and the correlation between the pixel and its neighbors to image segmentation. Let × = [ ( , )] × and denote an image of × pixels. The image has gray levels; that is, = {0, 1, 2, . . . , − 1}. The grey level of pixel at coordinate ( , ) is expressed as ( , ). Take a 2 × 2 slip window whose center pixel is ( , ). Define the 2 × 2 neighbor region which includes the four regions ( −1, ), ( , −1), ( , +1), and ( +1, ) as shown in Figure 3 In standard FCM image segmentation approach, only related to gray feature is the fuzzy membership between the clustering sample pixel and the clustering centre V . In fact, no matter wherever the sample pixel is in the even region or region polluted by noise or the edge region, the influence of other pixels on the sample pixel should be considered. Because of the strong relationship between the sample pixel and its neighbors, we use the relationship to estimate the fuzzy membership in order to reduce the influence of noise and other uncorrelated pixels.
In grey image × = [ ( , )] × , take a 3 × 3 slip window and regard the center pixel ( , ) as the sample pixel. According to (1), we can get the medium similarity measure of the sample pixel ( , ) and its neighbors. The larger the value is, the higher the similarity degree between the sample pixel ( , ) and its neighbors becomes. The smaller the value is, the lower the similarity degree between the sample pixel ( , ) and its neighbors becomes. Let ℎ( , ) denote the mean similarity degree between the sample pixel ( , ) and its 8 neighbors, as follows:  ℎ( , ) reflects the similarity degree between the sample pixel ( , ) and its neighbors. When the sample pixel ( , ) is a normal pixel, the relationship between the sample pixel ( , ) and its neighbors is strong and the value of ℎ( , ) is high. When the sample pixel ( , ) is a noise pixel, the relationship between the sample pixel ( , ) and its neighbors is weak and the value of ℎ( , ) is low. The strong related pixels play the key role in the classification while the weak related pixels have relatively small impact on the classification. As a result, we use ℎ( , ) to modify membership function to reduce the influence of noise and uncorrelated pixels and improve the accuracy of classification.
Firstly, according to (3), compute and row the ℎ( , ) of every sample pixel ( , ), and then we get one-dimensional matrix [ℎ( )] × , where ℎ( ) represents the mean similarity degree between the th pixel and its 8 neighbors. Secondly row the gray image × = [ ( , )] × , and then use the classical FCM to calculate the fuzzy membership and the mean fuzzy membership between the th pixel and its neighbors: where ( ) represents the 3 × 3 neighbor region of center pixel .

Definition 2. For pixel , we define
as the medium membership function. The which considers the spatial information, the similarity degree, and the correlation between the sample pixel and its neighbors can help classify the sample pixel into the correct region.
When the sample pixel has a high correlation with pixels in its 3 × 3 neighbor region, medium membership grade is mainly judged by the mean fuzzy membership . It shows that the sample pixel has a high medium membership to class . Otherwise, when the sample pixel has a low correlation with pixels in its 3 × 3 neighbor region, medium membership grade is mainly judged by the fuzzy membership grade . This is more reasonable. The classification approach which uses the medium membership instead of fuzzy membership to get the cluster centers by iterative computation is more likely to get the correct classification. medical image segmentation algorithm is implemented in the following steps.
(4) According to (4), compute the mean fuzzy membership between the pixel and its neighbors. (5) According to (5), compute the medium membership of every pixel. (6) According to (7), update the cluster centre matrix ( +1) : where is the medium membership of sample pixel.
where represents the classification which the th pixel belongs to.

Experimental Results.
Peak-value signal-to-noise (PSNR), partition coefficient V , and partition entropy V are chosen as the criterions to evaluate the performance of the MMTDFCM method: where ( , ) is the grey level of pixel in original image and ( , ) is the gray level of pixel in processed image; max is the maximum gray level of the image; and are the sizes of the image.
For a given cluster number and membership matrix [ ] × , partition coefficient V and partition entropy V are defined as follows: When it gets the best cluster result, partition coefficient V will be the maximum and partition entropy V will be the minimum.
An artificial image, an MR image, and an ROI image are chosen as experimental samples to verify the MMTDFCM algorithm described in this paper. Standard FCM is used in comparison with the experiment in order to evaluate the performance of the new algorithm in medical image segmentation. As the most common noises that emerged in medical images are Gaussian noise and salt and pepper noise, images with salt and pepper noise and Gaussian noise are chosen to test the antinoise ability of MMTDFCM. Fuzzy weighting exponent controls the share degree of clusters, and, in practice, the best value range of is considered as [1.5,2.5]. This also accords with conclusion of Pal's experiments [13]. Set the fuzzy weighting exponent = 2 and iterative threshold = 0.00001. The experimental results are shown in Figure 4, Figure 5, and Tables 1, 2, and 3.
The artificial image (128 × 128) includes three gray levels: 0, 160, and 255, and we set the cluster number = 3. The MR image (256 × 256) and the ROI image (128 × 128) include the background, white matter, gray matter, and cerebrospinal fluid, so we set the cluster number = 4.
The performance of the new algorithm can be evaluated through both subjective visual and objective quality.   As it can be seen, due to the influence of the noise, standard FCM incorrectly classifies some pixels and the segmentation images by the standard  FCM include much noise. The MMTDFCM which considers the spatial feature and the correlation between the pixel and its neighbors can get the higher quality of segmentation to the polluted images. MMTDFCM correctly classifies the background, white matter, gray matter, and cerebrospinal fluid. Moreover, the segmentation images have little noise. The segmentation effect is obviously better than standard FCM. Moreover, standard FCM is more dependent on the initial positions of the centers of the clusters.
(2) Better Objective Quality. Table 1 to Table 3 show the segmentation accuracy of images with noise. From the tables, it can be seen that the PSNR of new algorithm is higher than the standard FCM. It shows that the MMTDFCM has higher antinoise performance. The results also show that, in many cases, the MMTDFCM has a larger partition coefficient V and a smaller partition entropy V than the standard FCM. It represents that the segmentation result by MMTDFCM has more certainty and less fuzziness.

Conclusion
The standard FCM image segmentation algorithm is sensitive to noise because of not taking into account the spatial information of the image. This paper proposes an improved FCM medical segmentation algorithm based on MMTD. It uses medium similarity measure and correlation of the pixel and its neighbors and defines the medium membership function. Both the gray information and the spatial information are applied in MMTDFCM. As a result, the MMTDFCM can classify the pixels more correctly. The experimental results show that the proposed algorithm is more antinoise than the standard FCM image segmentation algorithm, with more certainty and less fuzziness. This will lead to its practicable and effective applications in medical image segmentation.