Mammogram Images Thresholding for Breast Cancer Detection Using Different Thresholding Methods

The purpose of this study is to apply different thresholding in mammogram images, and then we will determine which technique is the best in thresholding (extraction) malignant and benign tumors from the rest breast tissues. The used technique is Otsu method, because it is one of the most effective methods for most real world views with regard to uniformity and shape measures. Also, we present all the thresholding methods that used the concept of between class variance. We found from the experimental results that all the used thresholding techniques work well in detection normal breast tissues. But in abnormal tissues (breast tumors), we found that only neighborhood valley emphasis method gave best detection of malignant tumors. Also, the results demonstrate that variance and intensity contrast technique is the best in extraction the micro calcifications which represent the first signs of breast cancer.


Introduction
Normal cells of the human body grow and divide to generate new cells in order to meet the body needs.If these normal cells grow old or destroyed, they die.So the new cells take their place.Cancer occurs when these new cells are generated and the body does not need them; in addition, the old and destroyed cells do not die as they should, so these additional cells are constructed as a mass tissue known as a lump, growth or tumor [1].
One of the main lethal cancers is breast cancer [1,2].It occurs in both males and females, but the breast cancer in men is scarce [1,3].Around 25% of all cancers in woman are breast cancer, and approximately 20% of cancers causing death are breast cancers.It primarily sites in ducts (tubas that transfers milk to the nipple) and in the lobules (glands which make milk) in the breast by taking the form as micro-calcifications or masses [1].Commonly, breast cancer cannot be avoidable, but at least early detection gives a big opportunity of the treatment of this disease [2].One of the best tools for scanning the structure of the breast is mammogram [1,2].Mammogram images are classified depended on the categories as malignant, benign and normal.Where benign and malignant are abnormal; benign tumors can be treated (not cancerous).In contrast, to malignant tumors (cancerous), this can dam-age neighboring tissues and spread to the rest parts of body [1,2].Mammogram images are represented as very accurate and complex images to be interpreted [1,3].Moreover, many cases (about 25%) of cancers are not detected in the screen.Therefore, there is a necessary need to improve computer-aided diagnosis (CAD) systems to help a radiologist in his diagnosis for the cancer [3,4].For these reasons, our aim in this paper is to put the scope on applying different thresholding techniques on mammogram images as a try to determine which one presents best results.
This paper is organized as follows: Section 2 is a brief introduction of segmentation, and mention types of segmentation used in mammograms.It also defined the thresholding, and related works for thresholding used in mammogram images.Section 3 illustrated the formulation of the structure of an image.Section 4 describes Otsu method, and the techniques related to it.Section 5 is about the thresholding evaluation methods.Section 6 is a discussion for experimental results.Conclusions appear in Section 7.

Thresholding
Segmentation in image processing plays a central role in detection the region of interest from background.Its in-puts are images and the outputs are the properties were obtained from those images [1,2].Segmentation of intensity images (such as mammograms), contains four types of methods: thresholding techniques, boundary-based methods, region-based methods, and hybrid techniques that used both boundary and region criteria [1].
The thresholding in mammograms images is based on separated the histogram into background and breast tissues.Depending on the value of threshold all pixels less than the threshold are classified as background, and the reminder pixels are breast or vice versa [1][2][3][4][5].As a related work there are many researchers used (global or local) thresholds to split the breast region from background; like in [6] presents simple histogram thresholding.Others like in [7][8][9] used global thresholding to extract the regions of interest (ROI) from the background.In contrary, in [10] local thresholds are used; in addition to region growing and morphological filter.In [5] the researcher applied entropy threshold to segment microcalcifications in mammographic X-ray images.Finally, in [11] a comparison has been done of automatic thresholding methods including Otsu's method, Kapur 's method, Pun's method, Johannes's method Iterative Selection, Fuzzy Sets, and semi-supervised thresholding method-Black PCT.After the comparison among these methods is done, the experimental outcomes ensured that Otsu method is the best in determination tumor cells.
In general, the previous thresholding techniques have disadvantages including some need manual adjustment of thresholds, and others require local statistics to get the thresholds (local thresholds) [5].Therefore, we will use one of the best automatic thresholding selection techniques (Otsu method).Its concept based on selected the optimal threshold that maximizes the separability classes in gray levels [12].Furthermore, many thresholding techniques have been appeared to revise Otsu technique; like [13][14][15][16]; each one developed Otsu method in a specific way.This research focuses on implemented thresholding techniques (Otsu method and the techniques related to it) on mammograms.

Formulation
To analyze and process any image we should know that an image is generated from a set of pixels denoted as n; for each image level there are a set of pixels denoted as n i .Therefore, the total number of pixels is defined as: Gray level histogram is normalized and regarded as a probability distribution: The gray level of an image is [0… L−1].Where the gr obability of occurrence of the two classes can be de The mean and variance of the foreground and gr i (4) hi (5) It worth to mention that in each image there is a speci

Otsu Method
d on established the optimal thres-ay level 0 is the darkest and the gray level L−1 is the lightest.
The pr noted as the following: backound are denoted respectively as the following:  fic thresholding algorithm used to get an optimal threshold, which separated the object from the background.
Otsu method is base hold that (minimizing intra-class variance) or maximizing inter-class variance.The following equations will represent the within-class variance, and the between class variance respectively.
T (7) The final form of between-class variance can also de Otsu method is known as: It's simple, effective, easier to optimal threshold is displayed in Fi

Valley Emphasis Method
be noted as the following: apply and choosing the optimal threshold is automatically and stably [12].
The position of the gure 1 where first class represents the background, while the second class represents the object.
The reason for formulating valley that it can establish an optimal threshold for both bi-    t (9) In this technique the smallest probability of occurrence h

Neighborhood Valley Emphasis Method
ourne al and unimodal histogram.In other word, it extract large and small objects from the background.This technique can be considered as a modified of Otsu method, which failed in detection small objects.Establishing the optimal threshold will be based on two factors:  It (valley in the grey level histogram). Making between class variance is lar possible like in Otsu method.The valley-emphasis equation i ).The largest weight confirms that the chosen threshold is in the valley between two peaks for (bimodal histogram) as seen in Figure 2(a), or at the bottom rim of a single peak for (unimodal histogram) as seen in Figure 2(b) [13].
To give more reliable segmentation result; a neighb hood valley emphasis technique is presented.It developed the prior technique by computing the grey values for both the valley point and its neighborhood.So that, selection the optimal threshold will be based on maximum weight  of the valley and its neighborhood.

If the imag one dimensional histogram h(i)
; the e has ighborhood grey value h(i) of the grey level i is denoted as : Equation ( 10) refers to the sum of probability of occu ique is de-rrence for the neighborhood valley point; within the range n = 2m + 1.Where n is denoted as neighborhood length; usually it is odd number.
The neighborhood valley emphasis techn noted as the following:  The optimal threshold which increasing between variance 2 2 ξ t ( ) 1 t h   class will be denoted as the following: Equation (11)  This technique is suitable for processing images with big diversity between the object variance variance [14].

Threshol Intensity Contrast
A presented technique used both within-class va ntensity contrast at the sam This technique succeeded in extraction an object from complex background distribution in an ima sed on attempt to make the within class variance is smallest, and the intensity contrast is largest simultaneously.In addition, there is a parameter λ that represents the important factor in the formula.It's a weight that balanced the contribution for the (within class variance and intensity contrast).
The proposed formula is defined as the following: where 1  (t), 2  (t) are the mean intensities of the object and ba ound ckgr . is formulated from the following Equation (14).

 
In Equation ( 14), the first term refers to the p of class ariance Discrepancy Method ented images robability occurrence and the standard deviation (variance) of the first class, while the second term refers to the probability of class occurrence and the standard deviation (variance) of the second class.
The optimal threshold will dichotomize the image into object and background to get best segmentation result [15].

V
The presented technique is used for segm that have large variance discrepancy between the object and background.It computes both the variance sum and variances discrepancy to get effective segmentation outcomes.
The introduced technique is formulated as the following: and w is a main factor that determines the contributions of variance sum and variance discrepancy.
The optimal threshold will tend to make the effect of class variance sum is smallest, and in riance discrepancy at the same time [16].

Thresholding Evaluation Methods
The quality of thresholding technique is a critica various depending on the type of the thresholdin nique and the kind of image.In order to analyze the performance of the thresholding techniques, there are different evaluation methods used to measure their robustness and efficiency.In our study we used two evaluation methods Region Non-Uniformity (NU) and Inter-Region Contrast (GC) [17].Then, we compare the results of the five thresholding techniques to determine which technique is the best in determination the region of interest (object) from the background.

Region Non-Uniformity (NU)
This method measures the ability to dis the background and object in the thres good thresholded image should contain higher intra region uniformity, which is related to the similarity attribute about region element In the following NU Equation ( 17): σ 2 (t) denotes to the variance of the whole image, while σ 0 2 (t) denotes to the variance of the object (foreground).w 0 (t) denotes to the probability of occurrence of the object.NU equal to zero denotes to well thresholded image, but NU = 1 denotes to incorrect thresholded image [17].

Inter-Region Contrast (GC)
This method is very important in measure the contrast good thresholded degree in the thresholded image.A image should have higher contrast across adjacent regions.In the following GC Equation (18) the object average gray-level is known as μ 0 (t), and the background average gray-level is known as μ b (t) [17].

Experimental Results
In the following examples we applied the thresholding images in order to techniques on negative mammogram bring the part of interest from bright region (second mode) to dark region (first mode) and then determine the threshold values.It is worth to mention that the radiologist diagnosis the first Figure 3(a1) normal breast masses, Figure 4(a2) is an image for micro calcifications, while Figure 5(a3) represents breast tissues with small malignnant tumor.In the following examples we will show the performance of each thresholding techniques.
In the first example Figure 3(a1) is the original image of normal breast tissues.In general, and as we see all the five thresholding techniques detect the normal breast tissues.In Figure 3(b), the result of Otsu method that separated the object from the background successfully with T = 157.Figure 3(c), Valley emphasis method also extracted the object from the background with T = 150.Figure 3(d) The third method neighborhood valley emphasis method gave the same result of segmentation of the valley emphasis method with T = 149.Figure 3(e), Variance and intensity contrast method gave T = 141 to separate the object from the background.Finally, Figure 3(f) is the variance discrepancy method with T = 141 that also extracted the normal tissues from the background.The best region non uniformity NU = 0.112841 is introduced from Otsu method, on the other side the best inter region contrast GC = 0.556721 is get from both variance and intensity contrast technique and variance discrepancy technique.
Most early breast cancer can be diagnosed by detecting micro calcifications in the mammogram image.Otsu method did not detect the micro calcifications, and its optimal threshold T = 53.In Figure 4(c) valley emphasis method failed in isolated the micro calcifications from the background; its optimal threshold T = 49.The third method is neighborhood valley emphasis method with threshold T = 78 gives bad results as shown in Figure 4(d).Variance and intensity method succeed in isolation the micro calcifications from the complex intensity background.Its optimal threshold is get from λ = 0.05.Moreover, this method presents the smallest value of region non uniformity NU = 0.00799638 and the smallest     5(d) the threshold that is obtained from neighborhood valley emphasis method succeeded in detection the tumor from the background with optimal threshold T = 114 NU = 0.614376, GC = 0.606957, and the average = 0.610667.Moreover, the radiologist diagnoses the tumor with all its protrusions as malignant tumor.Table 1 displays the results of the five thresholding methods including the optimal threshold, region non uniformity, inter region contrast and the average.

Conclusion
In co   Otsu method.The ve big re peri e that wh f brig detected t bject insi he bright tis e breast, Further fi more, one of our efficient results represents that the micro calcifications are detected by using variance and intensity contrast.This contribution is very useful in detection breast cancer early.Then, this gives a big opportunity for the women patients to handle the breast cancer early.Also, we found, for breast tumors, only neighborhood valley emphasis method can detect the small region of bright tissue inside the breast, and it can help in other cases of diagnosis breast tumors.

σ
Refers to the square root of within-class variance; the   t w σ

Fig- ure 4
(a) is an example for the micro calcifications in the mammogram image.As we see in Figure 4(b)

Figure 5 (
b) Otsu method with T = 190 did not detect the tumor; moreover, the objec inder methods like Figure 5(c) valley emphasis method failed to extract the tumor from the rest normal tissues.

Figure 5 (
e) variance and intensity contrast method also failed in detection the tumor, and it reported incorrect threshold with white image.
l Research,