About this Journal Submit a Manuscript Table of Contents
Journal of Biomedicine and Biotechnology
Volume 2012 (2012), Article ID 830252, 7 pages
http://dx.doi.org/10.1155/2012/830252
Research Article

A Hybrid Technique for Medical Image Segmentation

1School of Electrical Engineering, University of Ulsan, Building 7, Room No. 308, 93 Daehak-ro, Nam-gu, Ulsan 680-749, Republic of Korea
2School of Electronics and Computer Engineering, Chonnam National University, Building 7, Room No. 506, 77 Yongbong-ro, Buk-gu, Gwangju 500-757, Republic of Korea

Received 22 May 2012; Revised 11 July 2012; Accepted 12 July 2012

Academic Editor: Tai Hoon Kim

Copyright © 2012 Alamgir Nyma et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Medical image segmentation is an essential and challenging aspect in computer-aided diagnosis and also in pattern recognition research. This paper proposes a hybrid method for magnetic resonance (MR) image segmentation. We first remove impulsive noise inherent in MR images by utilizing a vector median filter. Subsequently, Otsu thresholding is used as an initial coarse segmentation method that finds the homogeneous regions of the input image. Finally, an enhanced suppressed fuzzy c-means is used to partition brain MR images into multiple segments, which employs an optimal suppression factor for the perfect clustering in the given data set. To evaluate the robustness of the proposed approach in noisy environment, we add different types of noise and different amount of noise to T1-weighted brain MR images. Experimental results show that the proposed algorithm outperforms other FCM based algorithms in terms of segmentation accuracy for both noise-free and noise-inserted MR images.

1. Introduction

Image segmentation plays a pivotal role in computer-aided diagnosis and therapy. The objective of image segmentation is to partition an image into nonoverlapping, constituent regions that are homogeneous with respect to some attributes such as intensity and texture [1]. The representative image segmentation methods are generally based on the following techniques: thresholding, clustering, and edge detection and region extraction [2]. In this study, a combination of both thresholding and fuzzy clustering techniques will be considered because thresholding finds out the uniform regions for fuzzy clustering that makes the clustering performance better.

For thresholding methods, it is necessary to convert an image to a binary image depending on threshold values that extract the objects from the background [3]. The most convenient and widely used technique is histogram thresholding that is based on the shape properties of the histogram. The image histogram has distinct peaks, with each peak corresponding to one distinct region, and the valleys as the threshold values for separating these regions [4, 5]. Thresholding-based segmentation algorithms are then generally efficient in terms of computational complexity when compared to other segmentation methods, and Otsu’s clustering-based thresholding [6] is being one of the most representative methods for image segmentation.

In addition, fuzzy clustering has been widely studied and successfully applied in image clustering and segmentation [7]. Fuzzy c-means (FCM) clustering [8, 9], which is one of the basic methods of fuzzy clustering, is a soft segmentation method that has been used extensively to improve the compactness of the regions with its cluster validity and simple implementation. FCM depends on the Euclidean distance between pixels based on the assumption that each feature is of equal importance. However, in most real-world applications, features are not considered equally important. Thus, this assumption may seriously affect clustering performance. To improve the performance of FCM, many techniques have been proposed, such as rival checked FCM and suppressed FCM (SFCM), which integrates the hard c-means (HCM) and FCM in order to improve the convergence speed and clustering performance [10, 11]. Based on the advantages of thresholding and fuzzy clustering algorithms for image segmentation, some authors have proposed hybrid techniques combining various FCM-based methods with thresholding. Histogram thresholding using fuzzy theory was proposed by Tobias and Seara [12] in which thresholding was done according to the similarity between grey levels, and this similarity was accessed through a fuzzy measure. However, the application of this algorithm is limited based on the assumption that objects and background must occupy nonoverlapping regions. Chaabane Ben et al. proposed a hybrid method that combines automatic thresholding with FCM [13]. However, this technique yielded good results such that significant peaks and valleys are identified properly. Another hybrid approach was introduced by Tan and Isa, and it provided a good solution to overcome the FCM’s sensitiveness to the initialization condition of cluster centroids and selection of the number of clusters by using the histogram thresholding [14]. This algorithm, however, has one drawback such that some flat portions of the histogram curves had been recognized as dominating peaks. To improve the performance of medical image segmentation, we propose a hybrid technique using Otsu thresholding and enhanced SFCM (EnSFCM). Furthermore, we reduce impulsive noise that is widely presented in magnetic resonance (MR) images by employing vector median filtering.

The rest of this paper is organized as follows. Section 2 introduces the proposed image segmentation approach and Section 3 presents experimental results of the proposed approach and well-known image segmentation algorithms. Finally, Section 4 concludes this paper.

2. Proposed Image Segmentation Framework

The proposed image segmentation approach consists of vector median filtering, Otsu thresholding, and EnSFCM, as shown in Figure 1.

830252.fig.001
Figure 1: Process of the proposed image segmentation approach.

In the first step, vector median filtering is utilized to reduce impulsive noise in medical images. Otsu thresholding is then employed for rough segmentation of brain MR images, and finally EnSFCM is applied to have well-segmented images. More details about the proposed approach are described in the following sections.

2.1. Vector Median Filter

It is easily to see impulsive noise in medical images which is independent and uncorrelated to the image pixels and is also randomly distributed over the image. In the case of magnetic resonance (MR) images, uncertainty is widely presented because of impulsive noise in the partial volume effects originating from the low resolution of sensors [11, 15]. With these reasons, impulsive noise inherent in brain MR images can result in low segmentation performance. Thus this paper utilizes a vector median filter (VMF) to remove impulsive noise in images. VMF is a vector processing operator that has been introduced as an extension of scalar median filter and preserves the image without getting blurred and no shifting of boundary [16]. It approaches the problem of noise reduction by searching the most robust vector in the processing window. The process to find the most robust vector is as follows.

Denote by 𝐼 an image to be processed, and let 𝑊 be the processing window centered on the pixel under processing of size 𝑁×𝑁, 𝑁=3,5,7,, and so forth. Let us consider that each input vector 𝑥𝑖 is associated with the distance measure, 𝐿𝑖=𝑁𝑗=1𝑥𝑖𝑥𝑗𝛾for𝑖=1,2,,𝑁,(1) where 𝛾 represents the selected norm. The distance between two samples can be defined by 𝑥𝑖𝑥𝑗𝛾=𝑚𝑘=1||𝑥𝑖𝑘𝑥𝑗𝑘||𝛾1/𝛾,(2) where 𝑚 is the distance of the vectors and 𝑥𝑖𝑘 is the 𝑘th element of 𝑥𝑖. Thus distance 𝐿𝑖 serves as an ordering criterion of 𝐿𝑜𝐿1𝐿𝑁1 which implies the same ordering of the input vectors. The VMF output of the set {𝑥0,𝑥1,,𝑥𝑁} is defined as the sample 𝑥VMF{𝑥0,𝑥1,,𝑥𝑁} that satisfies the following condition: 𝑁𝑖=1𝑥VMF𝑥𝑖𝛾𝑁𝑖=1𝑥𝑗𝑥𝑖𝛾,for𝑗=1,2,,𝑁.(3)

2.2. Otsu Thresholding

The Otsu’s algorithm is one of the most referenced thresholding methods to partition images by automatically selecting threshold values from the histogram of the image. To find threshold values, the Otsu’s algorithm utilizes the variance property of the image because variance is the measure of uniformity; the greater value of variance represents the greater difference between the background and the object [4, 6]. Initially, two regions are separated by the intensity threshold, and then the optimal threshold is determined by minimizing the within-class variance or maximizing the between-class variance.

Assume that {𝑃𝑖}𝐿1𝑖=0 are the probabilities of the gray-level image histogram of an image, where 𝐿 is the range of intensity levels. We can calculate the probabilities of background (𝑃𝐵) and the probabilities of object (𝑃𝑂) of the image with a threshold 𝑡 as follows: 𝑃𝐵(𝑡)=𝑡𝑖=0𝑃𝑖,𝑃𝑂(𝑡)=1𝑃𝐵=𝐿1𝑖=𝑡+1𝑃𝑖.(4)

The mean associated with the background and the object can be further calculated using the following equations: 𝜇𝐵(𝑡)=𝑡𝑖=0𝑖𝑃𝑖𝑃𝐵,𝜇(𝑡)𝑂(𝑡)=𝐿1𝑖=𝑡+1𝑖𝑃𝑖𝑃𝑂(.𝑡)(5)

By using values of mean, we can compute variance as follows: 𝜎𝐵(𝑡)=𝑡𝑖=0𝑖𝜇𝐵2𝑃𝑖𝑃𝐵,𝜎(𝑡)𝑂(𝑡)=𝐿1𝑖=𝑡+1𝑖𝜇𝑂2𝑃𝑖𝑃𝑂(.𝑡)(6)

The between-class variance (𝜎between-class(𝑡)) which is the weighted variance of the cluster means around the overall mean is defined as follows: 𝜎between-class(𝑡)=𝑃𝐵𝜇𝐵(𝑡)𝜇2+𝑃𝑂𝜇𝑂(𝑡)𝜇2,(7) where 𝜇=𝐿1𝑖=0𝑖𝑃𝑖 is the global mean of the image. Furthermore, the within-class variance (𝜎within-class(𝑡)) can be expressed as follows: 𝜎within-class(𝑡)=𝑃𝐵(𝑡)×𝜎𝐵(𝑡)+𝑃𝑂(𝑡)×𝜎𝑂(𝑡).(8)

Finally, the optimal threshold value, 𝑡opt, can be determined by maximizing the between-class variance or equivalently minimizing the with-class variance as follows: 𝑡opt=argmax1𝑖𝐿1𝜎between-class,(𝑡)=argmin1𝑖𝐿1𝜎within-class(.𝑡)(9)

2.3. Enhanced Suppressed Fuzzy C-Means

Fuzzy c-means (FCM) is a data clustering algorithm in which data can belong to several groups based on the membership value, and it is one of the most prominent clustering techniques for image segmentation, feature extraction, and pattern recognition [9, 17]. However, the conventional FCM sometimes leads to accuracy degradation in segmentation. To solve this problem, suppressed FCM (SFCM) was proposed to improve the clustering performance as well as convergence speed, which introduces the suppression factor for establishing a relationship between hard c-means (HCM) and FCM clustering algorithms [10]. However, there is no general consensus that the suppression factor for SFCM is optimal in any sense or suitable for any given purpose. Thus, we propose an enhanced suppressed FCM (EnSFCM) to automatically select the suppression factor for SFCM. The selection procedure of the suppression factor, 𝛼, is based on the separation strength between clusters that is a time variant suppression rate. More details about EhSFCM are as follows.

Let 𝑋={𝑥1,𝑥2,𝑥3,,𝑥𝑛} where 𝑛 is the number of image pixels. The conventional FCM algorithm sorts the data set 𝑋 into 𝑐 clusters. The standard FCM objective function is defined as follows: 𝐽𝑚(𝑈,𝑉)=𝑐𝑛𝑖=1𝑘=1𝜇𝑚𝑖𝑘𝑑2𝑥𝑘,𝑣𝑖,(10) where 𝑑2(𝑥𝑘,𝑣𝑖) is the Euclidian distance between the data point 𝑥𝑘 and the centroid 𝑣𝑖 of the 𝑖th cluster and 𝜇𝑖𝑘 is the degree of membership of the data 𝑥𝑘 to the 𝑘th cluster. The parameter 𝑚 which is called the fuzzy factor controls the fuzziness of the resulting partition, (𝑚1), and 𝑐 is the total number of clusters. FCM clustering is an iteration-based clustering technique that produces an optimal number of 𝑐 classes by minimizing the objective function 𝐽𝑚(𝑈,𝑉) with updated values of 𝜇𝑖𝑘 and 𝑣𝑖 according to the following equations: 𝜇𝑖𝑘=𝑐𝑗=1𝑑2𝑥𝑘,𝑣𝑖𝑑2𝑥𝑘,𝑣𝑗1/(𝑚1)1,𝑣𝑖=𝑛𝑘=1𝜇𝑚𝑖𝑘𝑥𝑘𝑛𝑘=1𝜇𝑚𝑖𝑘.(11)

SFCM modifies the membership function in FCM by utilizing the suppression factor as follows:𝜇𝑝𝑘=1𝛼𝑖𝑝𝜇𝑖𝑘=1𝛼+𝛼𝜇𝑝𝑘,𝜇𝑝𝑘=𝛼𝜇𝑝𝑘,𝑖𝑝,(12) where 𝜇𝑝𝑘 refers to data point, 𝑥𝑘 belongs to the largest cluster 𝑝, and 𝛼 is the suppression factor which ranges in the interval [0,1]. SFCM then updates the value of 𝑣𝑖 with the new membership function. The pixel clustering iterations are terminated when the termination measurement max1𝑖𝑐𝑣𝑡𝑖𝑣𝑖𝑡1<𝜀 is satisfied, where 𝑣𝑡𝑖 are the current centroids for 1𝑖𝑐, 𝑣𝑖𝑡1 are the previous centroids for 1𝑖𝑐, and 𝜀 is a predefined termination threshold. To match the perfect clustering of the given data set, it is necessary to select an optimal suppression factor for SFCM. We define a new exponential function to select the optimal suppression factor, and it is automatically updated at each iteration: 𝛼=expmin𝑖𝑗𝑣𝑖𝑣𝑗2𝑚,(13) where 𝑣𝑖 is the centroid of the 𝑖th cluster, 𝑣𝑗 is the centroid of the 𝑗th cluster, and 𝑚 is the degree of fuzzification. As mentioned before, selection of the suppression factor, 𝛼, is very essential for SFCM since higher values of 𝛼 indicate that FCM is superior to HCM whereas lower values represent vice versa. It is evident that good clustering performance exhibits compact clusters by giving more distance between centroids. Likewise, both the fuzzification parameter and the suppression rate 𝛼 influence the learning rate of the algorithms. Thus we intuitively selected min𝑖𝑗𝑣𝑖𝑣𝑗2/𝑚, illustrating the separation strength between the clusters and the fuzzy factor that signifies the fuzziness of the membership values for the clustered data points.

3. Experimental Results

This section evaluates the performance of the proposed approach. To evaluate the performance of the proposed approach, we set some empirical parameters. Bezdek et al. experimentally determined the optimal intervals for the degree of fuzzification (𝑚) and the termination threshold (𝜀), which ranged from 1.1 to 5 and 0.01 to 0.0001, respectively [17]. In this paper, we selected 𝑚=2 and the termination threshold to be 0.0001.

3.1. Segmentation Results for Gray Matter and White Matter

We evaluate the correctness of the segmentation using real brain scans with ground truth by expert segmentations obtained from the Internet Brain Segmentation Repository (IBSR) [18]. Figures 2(a)-2(b) and Figures 3(a)-3(b) show examples of two slices and the manual labeling provided by the IBSR, whereas the rest of the images illustrates the results of clustering by using FCM [9], MSFCM [11], FCMT [13], SFCM [19], and the proposed approach.

830252.fig.002
Figure 2: Comparison of the segmentation results on a simulated brain MR image. (a) original T1-weighted image, (b) manual class labeling of gray matter (GM) and white matter (WM) slice regions; results obtained with (c) FCM, (d) MSFCM, (e) FCMT, (f) SFCM, and (g) the proposed approach.
830252.fig.003
Figure 3: Comparison of the segmentation results on a simulated brain MR image. (a) original T1-weighted image, (b) manual class labeling of gray matter (GM) and white matter (WM) slice regions; results obtained with (c) FCM, (d) MSFCM, (e) FCMT, (f) SFCM, and (g) the proposed approach.

As indicated in the figures by a red mark, other FCM-based segmentation approaches miss considerable pixels in the red-marked region. On the other hand, the proposed approach successfully avoids this classification error. In addition to qualitative results, we quantitatively evaluate and calculate the segmentation accuracy of the proposed approach, which is described below.

3.2. Segmentation Accuracy

The performances of the proposed and other FCM-based approaches are compared with respect to the optimal segmentation accuracy [20], which is defined as the sum of the correctly classified pixels divided by the sum of the total number of pixels: SegmentationAccuracy=𝑐𝑖=1𝐴𝑖𝐶𝑖𝑐𝑗=1𝐶𝑗,(14) where 𝑐 is the number of clusters, 𝐴𝑖 represents the set of pixels belonging to the 𝑖th cluster by the segmentation algorithm, and 𝐶𝑖 represents the set of pixels belonging to the 𝑖th cluster in the reference segmented image. To evaluate the robustness of the proposed segmentation approach in noisy environment, we add different types of noise and different amount of noise to a T1-weighted brain MR image, shown in Figure 2(a). All test images are corrupted by 8% and 10% Gaussian noise, speckle noise, and 10% and 12% salt and pepper noise as shown in Figure 4. Figure 4 also shows segmentation results of the proposed and other FCM-based approaches.

fig4
Figure 4: Segmentation results in noisy environment. (a) 8% Gaussian noise, (b) 10% Gaussian noise, (c) speckle noise, (d) 10% salt and pepper noise, and (e) 12% salt and pepper.

Figure 5 presents the quantitative results of these segmentation approaches for noise-free and noise-inserted brain MR images. The proposed algorithm increases 3.92% (GM) and 3.87% (WM) of the segmentation accuracy for a noise-free brain MR image compared to the baseline performance, where we used the average value of the segmentation accuracy as a baseline. For noise-inserted images, the proposed algorithm achieves 6.79%–14.73% (GM) and 6.42%–17.03% (WM) improvements in the segmentation accuracy, which is significant in the field of image segmentation for diagnosis purpose.

fig5
Figure 5: Segmentation accuracy results for gray matter (GM) and white matter (WM) in noisy environment by using FCM, MSFCM, FCMT, SFCM, and the proposed segmentation approach.

4. Conclusions

This paper presented a hybrid medical image segmentation technique. In the preprocessing step, we utilized vector median filtering to reduce the influence of noise. We then produced rough-segmented images using the Otsu thresholding method. To obtain well-segmented images, we finally used the enhanced suppressed FCM algorithm of the proposed approach. To match the perfect clustering of the given data set, we defined a new exponential function to automatically select the suppression factor for the suppressed FCM. To evaluate the robustness of the proposed segmentation approach in noisy environment, we added different types of noise and different amount of noise to the T1-weighted brain MR image and calculated the segmentation accuracy. Experimental results showed that the proposed method outperforms other FCM-based algorithms in terms of segmentation accuracy for both noise-free and noise-inserted images.

Acknowledgments

This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korean government (MEST) (no. 2012-000496, 2012-0001727).

References

  1. R. Gonzalez and R. Woods, Digital Image Processing, Prentice-Hall, Upper Saddle River, NJ, USA, 2nd edition, 2002.
  2. K. S. Fu and J. K. Mui, “A survey on image segmentation,” Pattern Recognition, vol. 13, no. 1, pp. 3–16, 1981. View at Scopus
  3. D. Jayadevappa, S. Srinivas Kumar, and D. S. Murthy, “A Hybrid segmentation model based on watershed and gradient vector flow for the detection of Brain tumor,” International Journal of Signal Processing, Image Processing and Pattern Recognition, vol. 2, no. 3, pp. 29–42, 2009.
  4. H. S. Prasantha, H. L. Shashidhara, K. N. B. Murthy, and L. G. Madhavi, “Medical image segmentation,” International Journal on Computer Science and Engineering, vol. 2, no. 4, pp. 1209–1218, 2010.
  5. Z. K. Huang, Y. M. Xie, D. H. Liu, and L. Y. Hou, “Using fuzzy c-means cluster for histogram-based color image segmentation,” in Proceedings of the International Conference on Information Technology and Computer Science (ITCS '09), pp. 597–600, July 2009. View at Publisher · View at Google Scholar · View at Scopus
  6. N. Otsu, “A Threshold selecting method from gray-level histograms,” IEEE Transactions on Systems, Man, and Cybernetics, vol. 9, no. 1, pp. 62–66, 1979. View at Scopus
  7. M. Tabakov, “A fuzzy clustering technique for medical image segmentation,” in Proceedings of the International Symposium on Evolving Fuzzy Systems (EFS '06), pp. 118–122, September 2006. View at Publisher · View at Google Scholar · View at Scopus
  8. J. C. Dunn, “A fuzzy relative of the ISODATA process and its use in detecting compact well-separated clusters,” Journal of Cybernetics, vol. 3, no. 3, pp. 32–57, 1973. View at Scopus
  9. J. C. Bezdek, Pattern Recognition With Fuzzy Objective Function Algorithms, Advanced Applications in Pattern Recognition, Springer, Berlin, Germany, 1981.
  10. J. L. Fan, W. Z. Zhen, and W. X. Xie, “Suppressed fuzzy c-means clustering algorithm,” Pattern Recognition Letters, vol. 24, no. 9-10, pp. 1607–1612, 2003. View at Publisher · View at Google Scholar · View at Scopus
  11. W. L. Hung, M. S. Yang, and D. H. Chen, “Parameter selection for suppressed fuzzy c-means with an application to MRI segmentation,” Pattern Recognition Letters, vol. 27, no. 5, pp. 424–438, 2006. View at Publisher · View at Google Scholar · View at Scopus
  12. O. J. Tobias and R. Seara, “Image segmentation by histogram thresholding using fuzzy sets,” IEEE Transactions on Image Processing, vol. 11, no. 12, pp. 1457–1465, 2002. View at Publisher · View at Google Scholar · View at Scopus
  13. S. Chaabane Ben, M. Sayadi, F. Fnaiech, and E. Brassart, “Color image segmentation using automatic thresholding and the fuzzy c-means techniques,” in Proceedings of the IEEE Mediterranean Electrotechnical Conference (MELECON '08), pp. 857–861, May 2008. View at Publisher · View at Google Scholar · View at Scopus
  14. K. S. Tan and N. A. M. Isa, “Color image segmentation using histogram thresholding Fuzzy C-means hybrid approach,” Pattern Recognition, vol. 44, no. 1, pp. 1–15, 2011. View at Publisher · View at Google Scholar · View at Scopus
  15. M. H. A. Hadi, K. Mortexza, and S. Y. Hadi, “Least mean square algorithm tuned by fuzzy C-means for impulsive noise suppression of gray-level images,” International Journal of Signal Processing, Image Processing and Pattern Recognition, vol. 3, no. 4, pp. 55–66, 2010.
  16. J. Astola, P. Haavisto, and Y. Neuvo, “Vector median filters,” Proceedings of the IEEE, vol. 78, no. 4, pp. 678–689, 1990. View at Publisher · View at Google Scholar · View at Scopus
  17. J. C. Bezdek, J. Keller, R. Krisnapuram, and N. Pal, Fuzzy Models and Algorithms For Pattern Recognition and Image Processing, Springer, Berlin, Germany, 2005.
  18. Internet Brain Segmentation Repository (IBSR).
  19. M. A. Ali, G. C. Karmakar, and L. S. Dooley, “Image segmentation using fuzzy clustering incorporating spatial information,” in Proceedings of the International Conference on Computer and Information Technology, pp. 878–881, February 2004. View at Scopus
  20. L. Chen, C. L. P. Chen, and M. Lu, “A multiple-Kernel Fuzzy C-means algorithm for image segmentation,” IEEE Transactions on Systems, Man, and Cybernetics, Part B, vol. 41, no. 5, pp. 1263–1274, 2011. View at Publisher · View at Google Scholar