I Introduction
Fractal theory plays a fundamental role as an auxiliary tool in the solution of problems in areas as different as Medicine TWZ07 ; LC10 ; LSSMDB10 , Physics SMS10 ; HWZ08 ; CCWH10 , Engineering CDHLAB03 ; W08 ; DAGGR09 , among many others. Particularly, in tasks involving texture analysis, fractal geometry is a powerful modelling tool, achieving interesting results in the description and discrimination of such textures.
In the last two decades, some different fractal approaches to deal with texture analysis have arisen, for instance, multifractals H01 ; LRAJ08 ; LGS00 , the multiscale fractal dimension MCSM02 ; CC00 , the fractal descriptors BPFC08 ; BCB09 ; PPFVOB05 ; FCB10 , among others. Here, we are interested in the fractal descriptors approach.
The main idea of fractal descriptors is to extract a set of features from the estimation of fractal dimension under different scales. Generally, the fractal dimension is based on a powerlaw relation which expresses the fractality of a structure as a function of measure scale. Unlike the fractal dimension which is a single value, the fractal descriptors are computed over the whole powerlaw curve FlorindoBCB12; FlorindoB11.
An example that illustrates the power of fractal descriptors is showed in BCB09
. In that solution, the values in the powerlaw of BouligandMinwkowski fractal dimension are used to compose a feature vector to discriminate among plant leaf textures. Actually, this method demonstrates to be successful in the discrimination of natural textures. Such kind of texture present an intrinsic selfsimilarity property which is notedly well represented by fractal modeling.
Despite their good results, conventional BouligandMinkowski fractal descriptors present still a limitation in the representation of textures, mainly when these textures present a higher degree of complexity. This limitation is due mainly to the fact that the descriptors are obtained from the global image, without a more specific treatment of local characteristics present in any real image. Thus, we can obtain more information by estimating those descriptors in different scales over the image.
Considering this assumption, the present work proposes a solution to extract fractal descriptors from a texture based on BouligandMinkowski method. Here, we propose the estimation of BouligandMinkowski descriptors at different scales (decomposition levels) of the image. The idea is to decompose recursively the image into 4 equal parts and, in each recursion step, we calculate an average and a deviation of the BouligandMinkowski descriptors. Thus, from both average and deviation descriptors, we extract entropy measures and compose the feature vector for the texture image.
The method is tested over wellknown benchmark texture datasets in a classification task and the results are compared to classical and stateoftheart texture features methods in the literature, like Gabor wavelets MM96 , Laws energy L84 , Gray Level Difference Matrix WDR76 , etc. The results confirmed the better accuracy of the proposed technique and pointed to the possibility of using the proposed method in a large number of problems involving the description and/or discrimination of textures.
This work id divided into 7 sections, including this introduction. The following provides mathematical background of fractal theory. The third section shows the original BouligandMinkowski fractal descriptors. The fourth presents the proposed method. The following explains the experiments. The sixth shows the results of experiments and the final section does the conclusions.
Ii Fractal Geometry
The literature shows some works applying fractal geometry to solve problems related to texture analysis in a large number of applications TWZ07 ; SMS10 ; CDHLAB03 . The importance of fractals in such kind of tasks is explained by the most flexible representation model provided by fractal theory. In this way, the fractal representation allows the extraction of measures which may describe more faithfully the original structure depicted in the texture image.
Most of such measures is based on the concept of fractal dimension (FD). The importance of FD is due to the fact that it captures the complexity of a fractal object or still, its spatial occupation. Furthermore, these properties are also related to the visual aspect of a texture. Thus, fractal geometry enables a link between the mathematical relations inside a pixel structure and the subjective concept of visual distinction. This link turns fractals into a particularly interesting tool for texture representation and description.
The original definition of fractal dimension is also known as HausdorffBesicovitch dimension . It is calculated for a set by
(1) 
in which is the dimensional Hausdorff measure, defined through:
(2) 
where corresponds to the diameter in , that is, .
In most practical applications, we are interested in calculating the fractal dimension of objects which are not exactly fractals. This is the case here where we are considering texture surfaces, which may be only approximated to fractal objects, if we weaken the infinite selfsimilarity criterion. A serious drawback in the above fractal dimension definition is that it cannot be applied to these cases where we have not a real fractal. To solve this issue, the literature presents a lot of estimation methods, which compute a fractal dimension value for the real world object, approximating the original fractal dimension concept. Most of such methods is based on the general expression:
(3) 
where is a specific measure (depending on the estimation method) of the object and is a scale parameter under which the measure is taken F86.
We may find a large number of fractal dimension estimation methods F86, like boxcouting, BouligandMinkowski, Fourier, etc. The present work is focused on BouligandMinkowski method.
Iii BouligandMinkowski Fractal Descriptors
The BouligandMinkowski fractal dimension is obtained by replacing in 3 by a dilation volume . So, in this approach, initially, we map the grayscale image onto a 3D surface:
(4) 
where:
(5) 
in which is the maximum pixel intensity.
Thus, the mapped surface is summited to a dilation process, by a dilation radius . Essentially, this operation consists in draw spheres with radius and with center in each point of . The dilation volume corresponds to the total amount of points pertaining to the union of the spheres. The BouligandMinkowski dimension is calculated from the curve . The Figure 1 illustrates the process.
Another possible interpretation for the dilation volume is that it corresponds to the number of points with a distance at most from the original object, in our case, the surface. In this way, the exact Euclidean distance transform (EDT) FCTB08 is used to optimize the computation of .
The EDT is a transform which maps each point in the 3D space to the distance of this point to a subset of the space. This subspace corresponds to the object of interest, defined by the user. Here, this subset is the mapped surface and the EDT for each point outside is defined through:
(6) 
where represents the euclidean distance.
In digital images like here, we use the exact EDT and the distances present discrete values :
(7) 
where
(8) 
The dilation volume is obtained through:
(9) 
where
(10) 
where:
(11) 
where
(12) 
Iv Proposed Method
Here, we propose the decomposition of the original texture image into decreasing cell sizes, followed by the calculus of BouligandMinkowski descriptors in each cell. The idea is in some way similar to that found in some classical multiscale approaches, like discrete wavelet transform or Gauss pyramid.
The essential idea is to divide recursively the image into 4 equal parts. Each step in this process constitutes a decomposition level. At each decomposition level, we take the average and the standard deviation of descriptors in each cell. Thus, we construct a feature vector from the entropy measure of such descriptors. Finally, we apply a simple attribute selection approach to the feature vector to compose the final descriptors.
Thus, we start with a digital image . This image is decomposed into levels , where is the maximum possible level in the image, given by . In each decomposition level, the image is partitioned into equal regions :
In each region , we apply the procedure described in the above section and obtain the BouligandMinkowski descriptors . For each level , we obtain the average descriptors and deviation descriptors :
In the following, we extract entropy features from both average and deviation descriptors in each level. Iniatially, for each component (index) of BouligandMinkowski average descriptors at all levels, we construct another vector , that is:
In the same fashion, we construct the vectors , from deviation descriptors:
Then, we compute one Shannon entropy value for each vector. In order to simplify the notation, we call a generic vector. The entropy is estimated through:
where is the length of .
The entropy feature vector is given by:
where is the number of components in average and deviation descriptors.
The final step consists in applying a basic feature selection approach to increase the performance and reduce the number of descriptors. We developed a selection based on the classifier method. Therefore, we computed the classification success rate
for each component of vector and sorted the success rates in descending order. We are particularly interested in the vector of sorting indices :Finally, we compose the proposed Multilevel Descriptors by indexing through the first indices in , where is the minimum necessary number of components to provide the best possible result. Formally, we have:
The Figure 3 shows a diagram depicting the whole process.
The Figure 4 illustrates the power of the proposed method in discriminating among textures from 3 classes in Brodatz benchmark dataset.
V Experiments
In order to verify the efficience of the proposed method, we compared its performance with other texture descriptors in the classification of well known texture data sets.
In the first experiment we used the Brodatz data set B66 . This is composed by 111 images photographed from an architecture book. In the database, each image is divided into 10 windows 200
200 and the images correspond to the classes. Brodatz data is broadly used as a benchmark set in computer vision and pattern recognition techniques, given the variety of characteristics found in its images, like variations in luminance, geometrical configurations, fidelity to real world textures, among others.
In the following, we classified the Vistex database Vistex , a set of natural color textures, composed by 54 classes with 16 samples in each class. Each sample is represented in a 128128 image. Here, we used the graylevel version of the texture images.
For a fair comparison, we extracted some different classical and stateoftheart texture descriptors found in the literature and applied to the classifier. The compared methods are Gabor wavelets MM96 , Cooccurrence matrix H67 , Laws energy L84 , Gray Level Difference Matrix (GLDM) WDR76 , multifractal spectrum H01 and original BouligandMinkowski fractal descriptors BCB09 . We used the Linear Discriminant Analysis (LDA) classifier DH00 and adopted the holdout technique as a statistical training and validating scheme.
Vi Results
The following Tables 2 and 3 show the results in terms of correctness rate and some other statistical metrics in the classification of the benchmark datasets. We employed the following metrics: Correctness Rate (CR), Kappa index (), Average Correctness Reliability (ACR), Average Error Reliability (AER), Average Error type 1 (AE1) and Average Error type 2 (AE2). The Table 1 shows a brief description of each metric. More details about each one may be found in BCB09 . We also show the number of descriptors (ND) employed by each compared approach.
Metric  What measures 

CR  Percentage of elements correctly classified 
Precision gain relative to a hypothetical random classification  
ACR  Average a posterioriprobability of correctly classified elements 
AER  Average a posteriori probability of misclassified elements 
AE1  Probability of elements classified as being from any class when pertain to the class 
AE2  Probability of elements classified as being from the class when pertain to any other class 
Method  ND  CR (%)  ACR  AER  AE1  AE2  

Gabor  20  90.09  0.90  0.98  0.84  0.08  0.10 
Cooccurrence  84  92.07  0.92  0.99  0.91  0.07  0.08 
GLDM  20  84.14  0.84  0.94  0.79  0.15  0.16 
Laws  15  87.03  0.87  0.93  0.68  0.11  0.13 
Multifractal  101  37.48  0.37  0.92  0.86  0.63  0.62 
BouligandMinkowski  85  98.92  0.99  1.00  0.88  0.01  0.01 
Proposed method  62  99.28  0.99  1.00  1.00  0.01  0.01 
The first point to be observed is that the proposed technique has overcome all the classical and stateoftheart descriptors. Although in Brodatz data the margin to improve the classification rate is small, we may notice that the proposed technique provided a more robust result and with a significant advantage over the other methods. Relative to the other statistical measures, we notice that they confirm the correctness efficiency. Particularly, the proposed method also presents a minimum error (1 and 2) and a perfect reliability (until the significance level adopted), both in correct and wrong classifications. We see that, although BouligandMinkowski shows similar values, it presents a significantly lower AER value, implying that the classifier confuses descriptors from a relevant number of classes.
Vistex dataset is a more complicated case once it was developed for color analysis approaches. Thus, the use of gray level descriptors is waited to present a defficient result. Nevertheless, even not using color properties, the multilevel descriptors achieved a good classification result, mainly relatively to other gray level texture approaches. Again, as in Brodatz set, the other statistical metrics support the correctness rate. Also, again, the AER in Gabor descriptors is meaningly lower than that of MLD descriptors.
Method  ND  CR (%)  ACR  AER  AE1  AE2  

Gabor  20  88.19  0.88  0.96  0.79  0.11  0.12 
Cooccurrence  24  79.63  0.79  0.95  0.79  0.17  0.20 
GLDM  20  67.36  0.67  0.84  0.67  0.31  0.33 
Laws  15  84.03  0.84  0.87  0.63  0.15  0.16 
Multifractal  101  32.41  0.32  0.89  0.82  0.63  0.68 
BouligandMinkowski  85  86.81  0.87  0.98  0.91  0.12  0.13 
Proposed method  101  92.82  0.93  0.99  0.94  0.06  0.07 
Generally speaking, we observe in both datasets that the proposed descriptors have presented a performance even better than BouligandMinkowski approach. BouligandMinkowski method has already demonstrated to be an efficient tool for the discrimination of natural textures. Such performance is explained by the dilation process in that method. As the surface points are dilated, some wavefronts start to emerge. The distribution of these wavefronts provide a rich description of the original arrangement of pixel intensities in the image. Moreover, it still gives information about physical characteristics, like luminance, roughness and even material composition. On the other hand, here we proposed an improvement to the original BouligandMinkowski approach. The present technique shows three main additions to the original method. The first is the multiscale decomposition. In this way, we are able to capture localized spatial details with a greater accuracy. The second is the use of deviation descriptors. Actually, BouligandMinkowski descriptors are in some sense similar in their general aspect. The deviation allows a more highlighted representation of patterns embedded in the original descriptors and, as a consequence, increases the power of the classifier technique. A third point is the use of entropy measure. Shannon entropy is a classical method to measure the information content in a data. In this way, here, the entropy plays the role of attenuating possible redundancies and describe the original texture with the minimum necessary number of descriptors.
Finally, the Figures 5 and 6 exhibit the confusion matrices for the 4 methods which presented the greater correctness rates in each dataset. In both figures, the matrix is represented through a surface in which the heights correspond to the number of samples classified to class and pertaining to class . In these surfaces, a good method is represented through a diagonal with a wall aspect (with a minimum amount of “holes”) and with a minimum number of peaks outside the diagonal. In this sense, for Brodatz dataset, we observe that the MLD matrix shows only one protuberant peak in the lower half of the matrix. In the case of Vistex, we also see a reduced number of peaks outside the diagonal.
Vii Conclusions
This work develops and study a new approach for texture descriptors based on fractal geometry, more specifically, on BouligandMinkowski fractal descriptors. The technique computes the BouligandMinkowski descriptors of an image under different decomposition levels and for each level we estimate the average and deviation descriptors. Thus, we extract statistical measures for each average and deviation. These measures compose the feature vector of the texture.
The method was tested on a classification task of benchmark texture datasets and compared to other classical texture descriptors approaches. The results demonstrated the higher accuracy of multilevel descriptors in this task. These results also illustrates the importance of a decomposition step in the application of BouligandMinkowski descriptors and the relevance of the deviation descriptors in the discrimination of more complex texture datasets.
Finally, this outcome suggests that the present technique is a powerful approach to describe and discriminate graylevel textures. This also points to the possibility of using the proposed method in a large number of problems involving texture analysis and related issues.
Viii Acknowledgements
Odemir M. Bruno gratefully acknowledges the financial support of CNPq (National Council for Scientific and Technological Development, Brazil) (Grant #308449/20100 and #473893/20100) and FAPESP (The State of São Paulo Research Foundation) (Grant # 2011/015231). João B. Florindo is grateful to CNPq(National Council for Scientific and Technological Development, Brazil) for his doctorate grant.
References
 (1) L. TianGang, S. Wang, N. Zhao, Fractal Research of Pathological Tissue Images, Computerized Medical Imaging and Graphics 31 (8) (2007) 665–671.
 (2) S. Lorthois, F. Cassot, Fractal analysis of vascular networks: Insights from morphogenesis, Journal of Theoretical Biology 262 (4) (2010) 614–633.
 (3) R. Lopes, M. Steinling, W. Szurhaj, S. Maouche, P. Dubois, N. Betrouni, Fractal features for localization of temporal lobe epileptic foci using SPECT imaging, Computers in Biology and Medicine 40 (5) (2010) 469–477.
 (4) E. I. Scarlat, M. Mihailescu, A. Sobetkii, Spatial frequency and fractal complexity in singletotriple beam holograms, Journal of Optoelectronics and Advanced Materials 12 (1) (2010) 105–109.
 (5) D. Han, M. Wang, J. Zhou, Fractal analysis of selfmixing speckle signal in velocity sensing, Optics Express 16 (5) (2008) 3204–3211.
 (6) W. Y. Chen, S. J. Chang, M. H. Weng, C. Y. Hung, Design of the fractalbased dualmode bandpass filter on ultra thin liquidcrystalpolymer substrate, Journal of electromagnetic waves and applications 24 (23) (2010) 391–399.
 (7) D. Chappard, I. Degasne, G. Hure, E. Legrand, M. Audran, M. Basle, Image analysis measurements of roughness by texture and fractal analysis correlate with contact profilometry, Biomaterials 24 (8) (2003) 1399–1407.
 (8) R. P. Wool, Twinkling Fractal Theory of the Glass Transition, Journal of Polymer Science Part B  Polymer Physics 46 (24) (2008) 2765–2778, Annual Meeting of the AmericanPhysicalSociety, New Orleans, LA, MAR 10, 2008.
 (9) I. Das, N. R. Agrawal, S. K. Gupta, S. K. Gupta, R. P. Rastogi, Fractal Growth Kinetics and Electric Potential Oscillations during Electropolymerization of Pyrrole, Journal of Physical Chemistry A 113 (18) (2009) 5296–5301.
 (10) D. Harte, Multifractals: theory and applications, Chapman and Hall/CRC, 2001.
 (11) B. Lashermes, S. G. Roux, P. Abry, S. Jaffard, Comprehensive multifractal analysis of turbulent velocity using the wavelet leaders, European Physical Journal B 61 (2) (2008) 201–215.
 (12) S. Lovejoy, P. Garrido, D. Schertzer, Multifractal absolute galactic luminosity distributions and the multifractal Hubble 3/2 law, Physica A  Statistical Mechanics and its Applications 287 (12) (2000) 49–82.
 (13) E. T. M. Manoel, L. da Fontoura Costa, J. Streicher, G. B. Müller, Multiscale fractal characterization of threedimensional gene expression data, in: SIBGRAPI, IEEE Computer Society, 2002, pp. 269–274.
 (14) L. da F. Costa, R. M. Cesar, Jr., Shape Analysis and Classification: Theory and Practice, CRC Press, 2000.
 (15) O. M. Bruno, R. de Oliveira Plotze, M. Falvo, M. de Castro, Fractal dimension applied to plant identification, Information Sciences 178 (12) (2008) 2722–2733.

(16)
A. R. Backes, D. Casanova, O. M. Bruno, Plant leaf identification based on volumetric fractal dimension, International Journal of Pattern Recognition and Artificial Intelligence (IJPRAI) 23 (6) (2009) 1145–1160.
 (17) R. O. Plotze, J. G. Padua, M. Falvo, M. L. C. Vieira, G. C. X. Oliveira, O. M. Bruno, Leaf shape analysis by the multiscale minkowski fractal dimension, a new morphometric method: a study in passiflora l. (passifloraceae), Canadian Journal of BotanyRevue Canadienne de Botanique 83 (3) (2005) 287–301.
 (18) J. B. Florindo, M. De Castro, O. M. Bruno, Enhancing Multiscale Fractal Descriptors Using Functional Data Analysis, International Journal of Bifurcation and Chaos 20 (11) (2010) 3443–3460.
 (19) B. Manjunath, W. Ma, Texture features for browsing and retrieval of image data, IEEE Transactions on Pattern Analysis and Machine Intelligence 18 (1996) 837–842.
 (20) M. Pietikäinen, A. Hadid, G. Zhao, T. Ahonen, Computer Vision Using Local Binary Patterns, Computational Imaging and Vision, Springer, 2011.
 (21) K. I. Laws, Goaldirected texturedimage segmentation, Tech. Rep. 334, AI Center, SRI International, 333 Ravenswood Ave., Menlo Park, CA 94025 (Sep 1984).
 (22) J. Weszka, C. Dyer, A. Rosenfeld, A comparative study of texture measures for terrain classification, SMC 6 (4) (1976) 269–286.
 (23) J. B. Florindo, A. R. Backes, M. de Castro, O. M. Bruno, A comparative study on multiscale fractal dimension descriptors, PATTERN RECOGNITION LETTERS 33 (6) (2012) 798–806.
 (24) J. B. Florindo, O. M. Bruno, Fractal descriptors in the fourier domain applied to color texture analysis, CHAOS 21 (4) (2011) 1–10.
 (25) P. Brodatz, Textures: A photographic album for artists and designers, Dover Publications, New York, 1966.

(26)
MIT, Mit
vistex texture database (2011).
URL http://vismod.media.mit.edu/vismod/imagery/VisionTexture/vistex.html  (27) R. M. Haralick, Statistical and structural approaches to texture, Proceedings of the IEEE 67 (5) (1979) 786–804.
 (28) R. O. Duda, P. E. Hart, Pattern Classification, John Wiley and Sons, New York, 2000.