Neutrosophic and Fuzzy C-Means Clustering for Breast Ultrasound Image Segmentation

Neutrosophic and Fuzzy C-Means Clustering for Breast Ultrasound Image Segmentation

Abstract— Breast ultrasound image segmentation is one of the most difficult tasks due to its speckle noise, poor quality and location of the breast nodule. In this research, we proposed normalization algorithm to enhance image contrast in order to be segmented using neutrosophic and fuzzy c-means clustering. At first, the input image is filtered using speckle reduction anisotropic diffusion to reduce speckle noise and normalized to increase the contrast. Secondly, the normalized image is transferred to neutrosophic domain with three membership subset T, I and F to define the nodule area. Finally, the fuzzy c-mean method is used to segment the nodule area from the background. To evaluate and compare the performance of the proposed method, this research uses several measurements, namely Area Metric and Boundary Metric. The result shows that implementation of normalization improves the performance of segmentation results.

Keywords—Breast ultrasound; neutrosophic; fuzzy c-mean; segmentation;

I.     Introduction

According to the statistics, the most frequent diagnose cancer and the highest cause of cancer death among women worldwide is breast cancer. It is estimated that there are 1.7 million cases and 521,900 deaths in 2012. A quarter of all cancer cases and 15% of all cancer deaths among female are  accounted for breast cancer. However, mortality can be reduced by early detection and treatment appropriately [1].

One of the imaging modalities commonly used to detect and to classify the mass abnormalities is ultrasound image. Compared to other modalities, ultrasound has several benefits such as no radiation, safer, cheaper, faster and possible to increase the number of detected nodules. However, it is highly dependent on the operator capability and the experience of the radiologist [2][3]. Computer aided diagnosis (CAD) has a potential to detect abnormal masses and it can be used as a tool to obtain second opinion for the radiologist to make an accurate diagnosis [4][5].

Segmentation remains an active field in machine vision and image processing research. Moreover, it is one of the hardest tasks in computer vision systems design. Segmentation often becomes the vital first step which must be taken successfully before subsequent tasks, e.g. feature extraction, classification, and description [6].

In the last decades, the segmentation methods in breast ultrasound images have been proposes whether automatically, semi-automatically or manually. Huang et al [7] categorizes segmentation methods into seven groups: thresholding-based, clustering-based, watershed-based, graph-based methods, active contour model, Markov random field and neural network. Then, Cheng et al  [8] classified segmentation methods into four class, such as histogram thresholding method, active contour model, Markov random field, and neural network.

Neutrosophy is a statistic method to handle indeterminate condition, such condition that cannot be handled by fuzzy logic. The image is transformed into neutrosophic set based on neutrosophic domain. Several research using neutrosophic methods have been proposed. Akhtar et al [9] used neutrosophic and K-means algorithm, Cheng et al [10] used neutrosophic and fuzzy c-means (FCM) clustering algorithm and Zhang et al used neutrosophic and watershed method  [11]. All of them used non-medical images. Meanwhile Anter et al [12] used neutrosophic and FCM clustering , Shan et al used neutrosophic and l-means clustering [13] and Guo [14] et al used neutrosophic and level set. All of them used medical images.

Based on the reviewed literature, the segmentation in breast ultrasound is still big challanges for futher research. In this paper, we proposed a normalization method to enhance breast ultrasound image contrast in order to be segmented using neutrosophic and FCM algorithm.

The paper is organized as follows: Section II describes the methods. Section III discusses the evaluation, results and comparison. The conclusion of this study is presented in Section IV.

II.    Methods

A.    Pre-processing

The stages in pre-processing are cropping the image to the area of ​​interest (RoI), converting the image into a gray scale image, use normalization to contrast the image and then filter the image using the SRAD filter.

1)   Normalization image is scaling of linear pixel values ​​using the entire grayscale level range to get a sharper image. Contrast represents light and dark spread in the image. Images are grouped into 3 contrast categories, ie low contrast, normal contrast, and high contrast. This category is usually distinguished intuitively. The normalization process of the histogram is as follows:

  • Calculate the maximum grayscale (max_f) and the minimum grayscale (min_f) of the image.
  • The absolute difference between minimum and maximum grayscale is calculated based on (1)
nx,y 

= normalization image

2)   Speckle Reduction Anisotropic Diffusion (SRAD) Filter

Speckle noise is noise that occurs during image acquisition that’s multiplicative and locally correlated noise. Filtering reduces the speckle noise  without damaging the important features. SRAD filter is proposed by Yu, et al [15] combine anisotropic diffusion (AD) method proposed by Frost and Adaptive Mean Filter proposed by Lee. this technique based on Partial Differential Equation (PDE) and Minimum Mean Square Error (MMSE) [16].

The gradient based edge detector in the original anisotropic diffusion with the instantaneous coefficient  of variation is replaced by SRAD filter method to be suitable for speckle filtering [17]. Gradient magnitude and Laplacian operators as edge detection are instantaneous coefficient functions.

If the pixel of the image is represented as

To control how much smoothing performed on a pixel using instantaneous coefficient in the diffusion function is given by (4) [17].

of d-dimensional measured data. 

In this neutrosophic approach, fuzzy c-means method is used to segment the image with the following steps [11]:

  • Map and decide {T,F}.
  • Enhancement.
  • Find the threshold in T and F.
  • Define homogeneity in intensity domain and decide {I}.
  • Convert the image to a binary image based on {T, I, F}.
  • Apply the fuzzy c-means clustering to the converted binary image.

C.    Evaluation performance

1)   Area metrics can assess how much the nodule area is covered by the generated nodule area correctly and how much is covered wrongly evaluated. The false positive (FP), the false negative (FN) and dice coefficient (Dice) are calculated based on (10), (11) and (12), respectively

HD evaluate the defective possible disagreement between two contours whilst MDSS evaluate the disagreement averaged over two boundaries.

III.   Evaluation and Result

This research uses 62 breast ultrasound images which are obtained from 3 hospital databases and completed with anatomical pathology results. The examples of  breast nodule is shown in Fig. 1

j12_ROI g2_ROI

Fig. 1.  Breast nodule ultrasound

A.    Pre-processing

RoI is the process of determining an area containing nodule based on mark from radiologist. Several ultrasound images contain label and marker, therefore important to be eliminated. Adaptive median filtering affective to solve it without causing any blur effect. The results of these processes is shown in Fig. 2

D:Backup S2paperDATA RENIData uji fix17.bmp
(a) (b)
(c) (d)

Fig. 2. (a) Original image (b) RoI image (c) Gray scale image (d) Unmark image

The purpose of image normalization is to increase image contrast. SRAD filter can reduce the noise speckle without damaging edge of the object. Fig. 3 shows the results of the normalization process and the SRAD filter.

(a) (b) (c)

Fig. 3. (a) Original image (b) SRAD image (c) SRAD image with normalization

B.    Segmentation

A hybrid technique based on neutrosophic set and FCM clustering method algorithm to breast ultrasound image segmentation is proposed. In this paper, objects are T and background is F. The blurry edges are gradually changed from objects to background, and there are no clear boundaries between the objects and edges or between the background and edges. The blurry boundaries are defined in I. Fig. 4 shows the result of neutrosophic algorithm based on FCM clustering for nodule segmentation.

Fig. 5 shows the example of segmented nodule regions by two methods. From the segmentation results, the proposed methods can separate the nodule area appropriately and preserved the important area of diagnosis, i.e. edges and the boundary area of nodule. The red lines on original image represent the nodule boundaries that are obtained with the different pre-processing methods. The results of neutrosophic based on FCM with normalization methods are more exactly and more appropriate to the gold standard than neutrosophic based on FCM without normalization..

There are five measurement parameters used to evaluate quantitatively. It compares the segmentation results of the proposed method with the segmentation results by the radiologist. The radiologist manually outline the nodule boundaries as a gold standard for evaluation of nodule area.

Total 62 images were selected from 2 cases. The nodule boundaries were manually outlined by the researcher based on the gold standard from the radiologist. it was used as a performance evaluation for segmentation results. All data represent breast cancer characteristics on our data set. The corresponding images calculated by the parameters Dice, FN, FP, HD, and MSSD.

As shown in TABLE I Dice, FP, HD and MSSD were improved from 82.0%, 43.7%, 47.5 pixel and 334.0 using neutrosophic based on FCM without normalization to 87.0%, 25.4%, 37.9 pixel, and 160.5 pixel using neutrosophic based on FCM with normalization, respectively. But On the contrary, FN was degraded performance from 1.4% using neutrosophic based on FCM without normalization to 4.0% using neutrosophic based on FCM with normalization. However, performance degradation in FN is offset by improved FP performance.

The comparison demonstrates that the neutrosophic based on FCM with normalization method can achieve better performance than neutrosophic based on FCM without normalization method. We can draw a conclusion that neutrosophic based on FCM clustering with normalization method is decent to segment nodule in the breast ultrasound images.

(a) (b) (c)
(d) (e) (f)
(g) (h) (i)

Fig. 4. (a) preprocessing image (b) T-domain (c) F-domain (d) indeterminate image (e) FCM for T-domain (f) FCM for F-domain (g) binary image based on T,I,F (h) binary segmentation (i) nodule segmentation

D:Backup S2Semester 4Neutrohasilsegmentdatabaru58 dataoundary 58 fcm SRADj25_ROI_filt_watershed.bmp D:Backup S2Semester 4Neutrohasilsegmentdatabaru58 dataoundary 58 fcm NSRADj25_ROI_normfilt_fcm.bmp
D:Backup S2Semester 4Neutrohasilsegmentdatabaru102 dataoundary 102 fcm SRADunMark_T_ROI_N Sumiyati  jinak  oval tegas hipo  posterior enhancement paralel_filt_fcm.bmp D:Backup S2Semester 4Neutrohasilsegmentdatabaru102 dataoundary 102 fcm NSRADSRAD_unMark_T_ROI_N Sumiyati  jinak  oval tegas hipo  posterior enhancement paralel_fcmnlm.bmp
(a) (b) (c)

Fig. 5. Example of nodule segmentation with different methods: (a) radiologist’s (b) result by neutrosophic based on FCM without normalization (c) ) result by neutrosophic based on FCM with normalization

  1. Performance Evaluation of Segmentation Methods
Methods Average area  metrics Average boundary metrics
Dice (%) FN (%) FP (%) HD (pixel) MSSD (pixel)
Neutrosophic FCM without normalization 82.0 1.4 43.7 47.5 334.0
Neutrosophic FCM with normalization 87.0 4.0 25.4 37.9 160.5

IV.   Conlusion

This research demonstrated that neutrosophic based on FCM clustering method with normalization improved the nodule segmentation. Automatic and exactly segmentation is a crucial step for many image processing and CAD application in breast ultrasound, and the proposed methods will find more implementations in this area.

Acknowledgment

We would like to express our gratitude to the department of Radiology RSUP Sardjito for the data set, the radiologists for the very helpful discussion and colleagues in Intelligent System research group at the Department of Electrical Engineering and Information Technology, Universitas Gadjah Mada.

References

[1] L. A. Torre, F. Bray, R. L. Siegel, J. Ferlay, J. Lortet-tieulent, and A. Jemal, “Global Cancer Statistics, 2012,” CA a cancer J. Clin., vol. 65, no. 2, pp. 87–108, 2015.

[2] J. M. Baran and J. G. Webster, “Design of Low-Cost Portable Ultrasound Systems : Review,” in 31st Annual International Conference of the IEEE EMBS, 2009, vol. 53706, pp. 792–795.

[3] H. A. Nugroho, Y. Triyani, and M. Rahmawaty, “Performance Analysis of Filtering Techniques for Speckle Reduction on Breast Ultrasound Images,” in International Electronics Symposium (IES), 2016, pp. 454–458.

[4] D. Chen, Y. Huang, and S. Lin, “Computer-aided diagnosis with textural features for breast lesions in sonograms,” Comput. Med. Imaging Graph., vol. 35, no. 3, pp. 220–226, 2011.

[5] Y. Triyani, H. A. Nugroho, M. Rahmawaty, I. Ardiyanto, and L. Choridah, “Performance Analysis of Image Segmentation for Breast Ultrasound Images,” in International Conference on Information Technology and Electrical Engineering (ICITEE ), 2016, pp. 415–420.

[6] C. Solomon and T. Breckon, Fundamentals of Digital Image Processing, 1st ed. John Wiley & Sons, Ltd, 2011.

[7] Q. Huang, Y. Luo, and Q. Zhang, “Breast ultrasound image segmentation: a survey.,” Int. J. Comput. Assist. Radiol. Surg., vol. 12, no. 3, pp. 493–507, 2017.

[8] H. D. Cheng, J. Shan, W. Ju, Y. Guo, and L. Zhang, “Automated breast cancer detection and classification using ultrasound images: A survey,” Pattern Recognit., vol. 43, no. 1, pp. 299–317, Jan. 2010.

[9] N. Akhtar, N. Agarwal, and A. Burjwal, “K-mean algorithm for Image Segmentation using Neutrosophy,” 2014 International Conference on Advances in Computing, Communications and Informatics (ICACCI). pp. 2417–2421, 2014.

[10] H. D. Cheng, Y. Guo, and Y. Zhang, “A Novel Image Segmentation Approach based on Neutrosophic Set and Improved Fuzzy C-Means Algorithm.,” New Math. Nat. Comput., vol. 7, no. 1, pp. 155–171, Mar. 2011.

[11] M. Zhang, L. Zhang, and H. D. Cheng, “A neutrosophic approach to image segmentation based on watershed method,” Signal Processing, vol. 90, no. 5, pp. 1510–1517, May 2010.

[12] A. M. Anter, A. E. Hassanien, M. A. A. ElSoud, and M. F. Tolba, “Neutrosophic Sets and Fuzzy C-Means Clustering for Improving CT Liver Image Segmentation BT  – Proceedings of the Fifth International Conference on Innovations in Bio-Inspired Computing and Applications IBICA 2014,” P. Kömer, A. Abraham, and V. Snášel, Eds. Cham: Springer International Publishing, 2014, pp. 193–203.

[13] J. Shan, H. D. Cheng, and Y. Wang, “A novel segmentation method for breast ultrasound images based on neutrosophic l-means clustering,” Med. Phys., vol. 39, no. 9, pp. 5669–5682, 2012.

[14] Y. Guo, A. Şengür, and J.-W. Tian, “A novel breast ultrasound image segmentation algorithm based on neutrosophic similarity score and level set,” Comput. Methods Programs Biomed., vol. 123, pp. 43–53, Jan. 2016.

[15] Y. Yu and S. T. Acton, “Speckle reducing anisotropic diffusion,” IEEE Transactions on Image Processing, vol. 11, no. 11. pp. 1260–1270, 2002.

[16] A. Balodi, M. L. Dewal, and A. Rawat, “Comparison of despeckle filters for ultrasound images,” in (INDIACom), 2015 2nd International Conference on “Computing for Sustainable Global Development,” 2015, pp. 1919–1924.

[17] C. P. Loizou and CS Pattichis, Despeckle Filtering for Ultrasound Imaging and Video Volume I : Algorithms and Software, Second. Cyprus: Morgan & Claypool Publishers, 2015.

[18] W. Jatmiko, P. Mursanto, B. Hardian, A. Bowolaksono, B. Wiweko, M. A. Akbar, I. P. Satwika, Z. Immadudin, M. sakti Alvissalim, I. Habibie, M. A. Ma’sum, and M. N. Kurniawan, Teknik Biomedis : Teori dan Aplikasi. Jakarta: UIP, 2013.

[19] J. C. Dunn, “A fuzzy relative of the ISODATA process and its use in detecting compact well-separated clusters,” 1973.

[20] J. C. Bezdek, Pattern Recognition with Fuzzy Objective Function Algorithms. Springer Berlin Heidelberg, 1981.

Professor

You must be logged in to post a comment