a study developed by american cancer society estimates that between one in eight and one in twelve will be diagnosed with breast cancer in their life time . mostly , 25% of all cases of breast cancer deaths occur if women were diagnosed between the age of 40 and 49 . although breast cancer incidence has increased over the last decade , breast cancer mortality has declined among women of all ages , thanks to the development of both breast cancer treatment and mammography screening . among the different imaging modalities used for the detection of breast cancer , demonstrate that the ratio of breast tumor detection in vermont ( usa ) , when applying screening mammography , increased from 2% to 36% between 1995 and 1999 . nowadays , the digital mammography gives the opportunity of increasing the use of the cad systems in order to help the expert radiologists in the interpretation and diagnosis of mammograms . however , the rapid improvement of full digital mammography has been accompanied by natural increase of such systems . the cad is a set of tools developed to help radiologists in the detection and interpretation of mammographic images . back in 2001 , freer and ulissey   have proposed an algorithm using substantial dataset containing 12,860 cases and have concluded that the application of cad in the analysis of screening mammograms may increase the malignancies detection at an early stage . the main disadvantage of existent cad systems is the lack of general algorithms producing good results for all cases and images . we will focus in this paper on mass - related lesions . in order to develop an improved computer - aided clinical decision classifying the tumor and identifying the stage of the cancer , we must ensure whether it is an area which contains a mass or not . our work focuses on classification of the tissue of the breast as mass or mass free . so , mass detection system is able to assist health professionals in finding out mass abnormalities in mammograms . however , some recent approaches have used multiple views , but this has three main drawbacks.image views must be properly registered to allow a correct comparison of regions.gray level values must also be correctly registered . there are some cases where comparison is not possible , because no correspondance between pixels can be done ( e.g. , the case of patients which have suffered from a previous breast surgery ) . there are some cases where comparison is not possible , because no correspondance between pixels can be done ( e.g. , the case of patients which have suffered from a previous breast surgery ) . it is important to note that algorithms typically working with one view can always be applied to multiple views . textural information has already been used to solve this problem and has been introduced in several works [ 6 , 10 ] . results have indicated that the use of lbp and its extensions has been promising in different comparative studies and has been applied in different texture analysis tasks [ 6 , 19 ] . however , lbp descriptors are not able to model mammogram texture because they are mapping only the differences of pixel gray level values . so , we will work here on a new approach taking into consideration the texture whole information , the local difference and local gray values as features , namely , gray level and local difference features ( glld ) . accordingly , we intend to investigate the efficiency of the glld based approach as a method of feature extraction . we perform our experiments on a set of 1000 rois obtained from the ddsm database . this paper , using a single view , proposes a new cad methodology in order to achieve better performances in terms of false negative and false positive using glld operator . section 2 shows some related works on mass detection in mammogram images . in section 3 , we present a brief review of lbp operator and an analysis of our glld based approach . in section 4 , a brief description of a set of classification methods is given , namely , support vector machine ( svm ) , k - nearest neighbors ( knns ) , and artificial neural network .   section 5 is reserved to the validation of the glld proposed technique simulation , and results and discussion are conducted . in the last section , we summarize the paper contribution and end our work by some concluding remarks and future work . several image processing techniques have been formulated as tools that can assist early automatic mass detection [ 20 , 21 ] . algorithms for mammographic mass detection using a single image view are based on a characteristic classifier scheme : for a given database consisting of known cases , the decision making system learns how to distinguish between the two kinds of rois ( mass and nonmass roi ) . thereafter , once the given system has been trained , a new roi can be rightly classified . among all these detection algorithms the first one includes the algorithms which extract features usually related to their texture from the roi and then trains a classifier . each new roi is compared to all the remaining roi images obtained from the database in order to be finally classified as mass or nonmass . table 2 shows different approaches applying this strategy . on the one hand , qian et al . have analyzed the implementation of an adaptive cad to develop a fully automatic procedure for mass segmentation and classification which consists in training a novel kalman - filtering neural network to classify features extracted from wavelet decomposition . on the other hand , christoyianni et al .   have extracted features based on independent component analysis ( ica ) , gray level and texture , in order to train the ann as a classifier . furthermore , they have applied the principal component analysis ( pca ) for the preprocessing step to overcome the problem of complexity and of increasing dimensionality . the latter is based on the translation of eigenfaces approach for face detection / classification problems to the mass detection one . as result of such transformation , they have obtained a vector which describes the contribution of each eigenrau for the representation of the corresponding image . they have used these vectors in the construction of the models for the step of training . in , oliver et al . have extended their proposed method based on pca by using the two - dimensional pca ( 2dpca ) technique . varela et al .   have proposed a methodology based on extracting gray level as well as morphological features and classifying , using ann , the new roi . leonardo et al .   have proposed an algorithm for the detection of masses in mammographic images . the technique is based on the use of textural and shape measures for k - means clustering algorithm and the svm , aiming at detecting masses in mammographic images .   as shown in table 2 , the proposed approach of chang et al .   and tourassi et al . they have proposed for classification purposes to undertake a comparison of the new roi with the remaining rois in the database composed of rois depicting masses . the difference between these works appears in the similarity measure function . as indicated in table 2 , there are only limited publications trying to detect masses using template - matching based methods . from the two tables , we can conclude that one of the main dissimilarity among these recent works is the ratio between the rois depicting abnormality and the total number of cropped images . it is important to note that when the number of normal rois increases , the number of rois wrongly classified is likely to increase . one should remember that the purpose of this work is the classification of mammographic masses and normal breast tissue . all the developed methods allow the tradeoff between the reduction of false positive fraction and the increase of false negative fraction . such trade off can be ensured when using the receiver operating characteristics ( roc )   in the performance evaluation step . the roc curve is a graphical curve representing the true positive rate ( sensitivity ) versus the false positive rate ( 100 specificity ) , extensively used in classifier performance evaluation . points representing roc curve correspond to sensitivity / specificity pair representing a particular decision threshold . the auc ( known as az ) is an information about the overall performance of the approach . furthermore , the latter is a metric which can be used to compare different features , and it allows the reduction of the roc curve to a single value summarizing expected performance . most of the approaches in the first strategy have the drawback that a large number of features need to be calculated but only the most discriminant will be selected [ 7 , 8 ] . besides , for the second strategy , the used similarity function measure for classifying needs to be recomputed for each element . in our paper , to overcome such limitations , the lbp operator has been investigated with the idea of performing gray scale invariant texture analysis . the latter has proved to be relevant in many applications knowing that the gray level information is of great importance in mammography , our approach will add to lbp absolute gray level information rather than gray level difference . we will focus in our approach on making use of small size feature vector as well as possible . in the following section , we will introduce our mass detection methodology based on the glld for the extraction of texture features obtained from the rois . texture classification is nowadays a challenging problem . it is an active topic in computer vision research . early methods of texture classification are based on statistical analysis of images with different textures . the most representative ones are the cooccurrence matrix method   and filtering for texture - classification methods . at an early stage , exploratory models were developed to investigate rotation invariance in texture classification , such as hidden markov model   and gaussian markov random field . varma and zisserman   have proposed to learn from a training set a rotation invariant texton and to classify the obtained texture according to its texton distribution . varma and zisserman   have later proposed to use the image patch in order to represent features directly . some recently proposed works have been developed for scale as well as affine invariant texture classification . later , ojala et al .   have proposed the lbp histogram application in order to achieve a rotation invariant texture classification . it is worth noting that the lbp is efficient in describing local image pattern and its performance in computer vision and pattern recognition is promising . however , it still needs to be improved for mammography texture modeling . in order to generate texton , ojala et al . have applied the absolute gray level difference ( agld ) between each pixel and its neighbors . . then , ojala et al .   have proposed the lbp using the sign of the difference for the representation of local patterns . in , ojala et al . suggested to use signed gray level difference ( sgld ) and its multidimensional distribution for the description of texture and considered lbp as a simplification of sgld . with such variants of lbp , there still remain questions that need answering , such as what information is lost in the considered code ? how to represent the missing information to obtain better texture modeling ? here , we propose a new feature extractor to improve the system performance , based on glld features . the lbp operator used eight neighboring pixels when considering the center gray value as threshold . this operator generates  1  if the considered neighbor value is greater or equal to that of the center . otherwise , it generates  0 .  accordingly and referring to figure 1 , lbp   code may be computed as follows : 
 ( 2)lbpp , r=p=0p1s(gpgc)2p , s(x)={1,x00,x<0 , 
 where gc corresponds to the gray value of the central pixel , gp corresponds to the value of its neighbors , and ( p = 0,1 , p  1 ) and p , r correspond to the number of neighbors and to the radius of the neighborhood , respectively . gp coordinates are ( ( r  cos(2p / p ) and r  sin(2p / p ) ) . if neighbors are not in the image grids , their gray values may be estimated by interpolation . after identifying lbp pattern of each pixel ( i , j ) , we associate lbp histogram to the whole image , with a given image size ( n1n2 ) as 
 ( 3)hlbp(k)=i=1  n1j=1n2f(lbpp , r(i , j),k ) ,  k[0,k],f(x , y)={1,x = y0,otherwise , 
 where k corresponds to the maximum gray level value . let u be a function corresponding to the value of an lbp pattern , it is defined as the number of transitions ( i.e. , change from 0 to 1 or 1 to 0 ) in the following pattern : 
 ( 4)u(lbpp , r)=|s(gp1gc)s(g0gc)|  + p=1p1|s(gpgc)s(gp1gc)| . patterns corresponding to limited transitions or discontinuities are with u  2 , in a binary presentation . the mapping from the original lbpp , r to lbpp , r , knowing that the superscript u2 refers to uniform patterns , may be implemented using a look - up table containing 2 elements . a local rotation invariant pattern is defined as follows : 
 ( 5)lbpp , rriu2={p=0p1|s(gpgc)|if  u(lbpp , r)2p+1otherwise . the mapping from lbpp , r to lbpp , r , knowing that the superscript riu2 corresponds to rotation invariant uniform patterns , may be implemented using a look up table . the main limitation using lbp code is that it may give the same results with two completely different gray levels when the differences with the neighbors are the same . knowing that for mammographic images , the gray level information is directly related to the breast tissue density , gray level and local difference are two important features of the texture which must be used together in order to have more accurate results . in our approach , we propose to calculate the average for each 3  3 neighborhood and to attribute it to the central pixel . given the new value of the central pixel gc  mean and its p circularly symmetric neighbor ( see figure 2 ) , the substraction of the value of gc  mean is presented as follows : 
 ( 6)t = t(gc  mean , g0gc  mean, thus , the difference between gc  mean and gp may be represented as diffp = gp  gc  mean and the local difference may be represented with a vector noted diffp knowing that diffp = [ diff0 ,  , diffp1 ] , diffp describes the local image structure around the gc  mean . because of its robustness and efficiency , the obtained vector diffp is decomposed of sign and modulus components in order to achieve much better performance in texture classification . in our proposal , sp corresponds to the sign of the differences , and it is obtained by thresholding with respect to the value of gc  mean as expressed in ( 8) . however , mp corresponds to the absolute value of diffp as expressed in ( 9 ) . we obtain , also , two vectors , the sign vector [ s0 ,  , sp1 ] and the modulus vector [ m0 ,  , mp1 ] , with 
 ( 7)diffp = mpsp , 
 knowing that 
 ( 8)sp={1,diffp01,diffp<0,(9)mp=|diffp| . aiming at recognizing robustly and efficiently the texture patterns , we should extract both absolute and relative features from pixel gray levels . the modulus component provides discriminant information to the sign component ; the intensity value of the central pixel corresponding to the mean value of its neighbors may also give us useful information [ 33 , 34 ] . it will also be seen that by coding the sign , the modulus , and the central gray level features into rotation invariant binary codes and fusing them , results may provide much better performance in mammogram texture classification than using each one by itself . this fusion provides useful information about local gray level which is so important in the stage of mass detection in mammographic images .    in this subsection , we present the gray level and local difference ( glld ) different processing steps to explore the proposed three features , which are illustrated in figure 4 . after that , in the selected roi , each central gray level corresponds to the mean of its neighbors and its local difference . the latter is decomposed into sign and modulus components as expressed in ( 7 ) . given a pixel in the image , the sign coding component is noted as ( sglld ) and is computed by comparing it with the values of its neighbors as follows : 
 ( 10)sglldp , r=p=0p1s(gpgc  mean)2p , 
 where s(x ) is defined by 
 ( 11)s(x)={1,x01,x<0 , 
 where gc  mean is the average value of the central pixel and its neighbors . inspired by the method of coding ( sglld ) , the coding of the magnitude component is noted as ( mglld ) and is defined as follows : 
 ( 12)mglldp , r=p=0p1t(mp , c)2p,(13)t(x , c)={1,xc0,x < c , 
 where c corresponds to a global gray level threshold which is determined adaptively . the new value of the central pixel , which expresses the gray level of the image , represents also a discriminant information . so , to make it consistent with the two previous operators sglld and mglld , we code it as 
 ( 14)cglldp , r = t(gc  mean , ci ) , 
 where t is already defined in ( 13 ) , ci corresponds to the threshold and is set as the mean gray level of the whole input image . figure 5 illustrates the image results after the application of the three operators and their fusion . for the three obtained codes , the rotation invariant version is defined to achieve rotation invariant classification . each code carries specific texture information , that is why we concatenate them to build the glld feature , which corresponds to a vector . the procedure consists in using the glld in order to build local descriptor of the obtained rois knowing that the concatenation leads to global description and the obtained global and local glld texture descriptor are , then , used as features for mass detection . the following images ( figure 7 ) illustrate the obtained histogram for differents rois . in this figure , we have considered three trivial examples as well as three challenging examples which were misclassified by a radiologist . all of the six roi examples where correctly classified by our glld texture features . in a further section , the last step of our proposal is mass classification . for the sake of generality and for doing a best choice of the classifier , an investigation of three classifiers will be undertaken , namely , support vector machine ( svm ) , k - nearest neighbors ( knns ) and artificial neural network . it learns how to discriminate between positive and negativ(in our case mass and non mass ) , by finding a hyperplane as a decision surface separating the classes . the svm uses an optimization method identifying the support vectors si , the weights ai , and the bias b which are used for the classification of the vectors x according to the following equation : 
 ( 15)c(x)=iai(si , x)+b , 
 where  corresponds to a kernel function . then , if c  0 , x is classified as a member belonging to the first class . otherwise , it is classified as a member belonging to the second class . since knn is memory based , no models need to be trained . for a given instance x , the knn first finds the k closest training points with respect to a particular distance metric . then , it uses its labels in order to classify the instance x by majority vote . in this study , we use the euclidian distance to determine the nearest neighbors of the query element , and k is used as a training parameter . for each element x , the output score corresponds to the ratio of the winning class elements among the total number of neighborhood in the corresponding dataset . ann has been widely used in many applications where the expert knowledge is not clearly defined . the idea of the ann has been inspired from the biological nervous system and has been successfully applied in medical imaging . this technique is based on the adjustment of weights between the neurons for any input - output function approximation . therefore ann , has been widely used in digital mammography to mimic this computational power and the perception capabilities of human brain . two basic types of ann , the multilayer perceptron ( mlp ) as well as the radial basis function network ( rbf ) , are frequently used in recent works . on the one hand , multilayer perceptrons ( mlps ) their principal advantage is the ease of use and the approximation capability of any input / output map . on the other hand , radial basis function ( rbf ) networks are nonlinear hybrid networks containing a single hidden layer of processing elements . this type of ann is in generally used when the number of samples is so small ( < 100 ) . so , the limitation of the rbf neural network is that it is very sensitive to the dimensionality and has more and more difficulties if the number of units is large . based on this assumption , and knowing that the glld feature size is of 1352 , we intend to investigate the mlp to exploit the results using the ann . let us consider x = ( x1 , x2 ,  , xd ) the input vector ,  = ( 1 , 2 ,  , d ) , the weight vector , and g(x ) = ( 1 + e ) the activation function which corresponds to a sigmoid function , and the network output is thus defined as follows : 
 ( 16)y = g(txb)=g(i=1dixib ) . 
 for each roi sample , glld features are computed and used in the classification step as inputs of the neural network . the evaluation of the effectiveness of the training is based on the measure of the network relative error as follows : 
 ( 17)e=i=1n(yt)number  of  applied  samples , 
 where y is the roi corresponding to masses or nonmasses resulting from ann and t corresponds to the target . after the training step , generalization errors may be evaluated for various features and network conditions . figure 9 maps the different steps of the proposed method .   in the following section , the obtained results for ann , svm , and knn classifiers will be illustrated for comparison purposes . we then do an investigation of the feature relevance , by using each of the proposed features ( sglld , mglld , cglld ) separately as input vector to the classifier . in further step , we made the classifier input a concatenated vector made up with different feature vectors . finally , a comparative study of our proposal to those in the state - of - the - art will be done for a fair evaluation . our approach has been evaluated based on publicly available database taken from the ddsm database . a case is the collection of all information to the mammography exam of one patient . each case in the ddsm database contains two images , of each breast , that is , in each case the mammograms include a craniocaudal and mediolateral oblique view ( cc and mlo , resp . ) . the ddsm database provides the metadata ( date of study , breast density , assessment categories , etc . ) of each abnormality using the breast imaging reporting and data system ( bi - rads ) lexicon , it provides , also , the corresponding chain codes of the suspicious regions . with these chain codes , the outlines of the abnormalities may be identified . the ddsm provide delineations of mass regions . however , precision of such delineations is not adequate for validation in our approach , since it was done on downscaled images of ddsm database ( by factor of 8) , see section 2 . therefore , we wore based in the extraction of roi 's on manual segmentation entertained by two expert radiologists of more than 30 years of clinical experience from the farabi imaging . we should also notice that all the considered masses in ddsm are biopsy proven ones . these rois were randomly selected and separated into two sets : 500 samples for training and 500 samples for tests . in the training set as well as in the setting set the evaluation of our mass detection algorithm is performed by applying a leave - one - out methodology , where the input roi is classified by using the appropriate classification method and the procedure is reapplied for all the remaining rois used as input .   from the results presented in table 4 , we can conclude that with ( p , r ) = ( 24,3 ) , the area under curve for the glld is increased from 0.93 to 0.95 . , glld24,3 will be used .    from the comparative study , as shown in table 5 , we note that the ann provides the best results . as illustrated in figure 10 and table 6 the cad system achieves better performance ( az = 0.93 ) when using the sign component than the modulus component . however , their fusion may provide much better results in texture classification than using either sign or modulus ( az = 0.95 ) . the auc of the glld24,3 after the fusion of the three operators sglld24,3 , mglld24,3 , and cglld24,3 using the ann as classifier is about ( az = 0.95 ) for the experimental set . as can be noted from table 6 , the glld24,3 feature provides useful information about local gray level which is the most significant one for mass detection in mammographic images . based on the size of the lesion , we use six different group of roi images , which is an important aspect for the correct classification of the masses . these classes correspond to the following specified mass sizes intervals : size 1 : < 10  mm , size 2 : ( 1060 )  mm , size 3 : ( 60120 )  mm , size 4  : ( 120190 )  mm , size 5 : ( 190,270 )  mm , and size 6 : > 270  mm . however , the used numbers of masses in each class size were , respectively , 28 , 32 , 37 , 57 , 69 , and 33 masses . table 7 illustrates the az values for each class of roi image sizes and the obtained mean az values . we include in this table a quantitative comparison with the work of oliver et al . [ 8 , 9 ] have used our database of roi ratio ( 1/3 ) , the same specified mass size intervals , and the same number of masses in each class size . the results presented in table 7 have shown that the glld24,3 features are effective for mass detection at different roi image sizes , and the latter is an important aspect for correct classification of the masses . our method proves its performance in the most difficult case , which correspond to the smaller masses . note that for this proposed ratio , better results are obtained for all the size intervals . table 8 shows the different az performance values for different approaches in the state - of - the - art presented in section 2 . such which were represented comparison demonstrates the effectiveness of the proposed glld operator in mass detection . for instance , [ 10 , 12 , 13 ] , which used ratio ( 1/1 ) the same as in our cases , obtained az values 0.83 , 0.90 , and 0.81 , respectively , as our obtained az = 0.95 . cad systems have been used and gained greater utility in recent years , as a second virtual reader for the medical images , contributing to increase an early detecting of breast cancer . the combined descriptors are , respectively , sglld24,3 , mglld24,3 , and cglld24,3 providing a final texture feature descriptor named glld , which will be used to classify the rois to masses and mass free . the ann classifier gives better performances in term of classification owing to its higher function approximation . such comparison illustrates that our proposed method leads to the best performance az = 0.95 . the specialists who i have collaborated with found that the proposed cad improved the sensitivity of mammography screening . in fact , cad system is useful in situations where there is a high interobserver variability , lack of trained observers , or impossibility to perform the double reading with two or more radiologists as stated in the birads categories . future work will be focused in the classification of masses into the four breast imaging - reporting and data system ( bi - rads ) categories .
during the last decade , several works have dealt with computer 
 automatic diagnosis ( cad ) of masses in digital mammograms . 
 generally , the main difficulty remains the detection of masses . 
 this work proposes an efficient methodology for mass detection 
 based on a new local feature extraction . 
 local binary pattern 
 ( lbp ) operator and its variants proposed by ojala are 
 a powerful tool for textures classification . 
 however , it has been 
 proved that such operators are not able to model at their own 
 texture masses . 
 we propose in this paper a new local pattern 
 model named gray level and local difference ( glld ) where 
 we take into consideration absolute gray level values as well as 
 local difference as local binary features . 
 artificial neural networks 
 ( anns ) , support vector machine ( svm ) , and k - nearest 
 neighbors ( knns ) are , then , used for classifying masses from 
 nonmasses , illustrating better performance of ann classifier . 
 we have used 1000 regions of interest ( rois ) 
 obtained from the digital database for screening mammography 
 ( ddsm ) . 
 the area under 
 the curve of the corresponding approach has been found to 
 be az = 0.95 for the mass detection step . 
 a comparative study with previous approaches proves that our approach offers the 
 best performances .