Fractal descriptors based on the probability dimension: a texture analysis and classification approach

In this work, we propose a novel technique for obtaining descriptors of gray-level texture images. The descriptors are provided by applying a multiscale transform to the fractal dimension of the image estimated through the probability (Voss) method. The effectiveness of the descriptors is verified in a classification task using benchmark over texture datasets. The results obtained demonstrate the efficiency of the proposed method as a tool for the description and discrimination of texture images.


page 3

page 6


Texture Analysis And Characterization Using Probability Fractal Descriptors

A gray-level image texture descriptors based on fractal dimension estima...

Fractal Descriptors of Texture Images Based on the Triangular Prism Dimension

This work presents a novel descriptor for texture images based on fracta...

Multiscale Fractal Descriptors Applied to Nanoscale Images

This work proposes the application of fractal descriptors to the analysi...

Gabor wavelets combined with volumetric fractal dimension applied to texture analysis

Texture analysis and classification remain as one of the biggest challen...

Enhancing fractal descriptors on images by combining boundary and interior of Minkowski dilation

This work proposes to obtain novel fractal descriptors from gray-level t...

Fusing Multiscale Texture and Residual Descriptors for Multilevel 2D Barcode Rebroadcasting Detection

Nowadays, 2D barcodes have been widely used for advertisement, mobile pa...

Graph entropies in texture segmentation of images

We study the applicability of a set of texture descriptors introduced in...

I Introduction

Fractals have played an important role in many areas with applications related to computer vision and pattern recognition

SMS10 ; HWZ08 ; CDHLAB03 ; W08 ; TWZ07 ; LC10 , owing to their flexibility in representing structures usually found in nature. In such objects, we observe different levels of detail at different scales, which are described in a straightforward manner by fractals, rather than through classical Euclidean geometry.

Most fractal-based techniques are based on the concept of fractal dimension. Altough this concept was originally defined only for mathematical fractal objects, it contains some properties that make it a very interesting descriptor for any object in the real world. Indeed, fractal dimension measures how the complexity (level of detail) of an object varies with scale, an effective and flexible means of quantifying how much space an object ocupies, as well as important physical and visual properties of the object, such as luminance and roughness.

Fractal techniques include the use of Multifractals H01 ; LRAJ08 ; LGS00 , Multiscale Fractal Dimension MCSM02 ; CC00 and Fractal Descriptors BPFC08 ; BCB09 ; PPFVOB05 ; FBCB12 . Here we are focus on the last approach, which has demonstrated the best results in texture classification FB13 . The main idea of fractal descriptors theory is to provide descriptors of an object represented in a digital image from the relation among fractal dimensions taken at different observation scales, thus these values provide a valuable information on the complexity of the object, in the sense that they capture the degree of detail at each scale. In this way, fractal descriptors are capable of quantifying important physical characteristics of the structure, as the fractal dimension, but presenting a richer information than can be provided by a single number (fractal dimension).

Although fractal descriptors have demonstrated to be a promising technique, we observe that they are defined mostly on well-known methods to estimate the fractal dimension. Here, we propose fractal descriptors based on a less known definition of fractal dimension: the probability dimension. This is a statistical approach, which measures the distribution of pixel intensities along the image. In this way, such descriptors can express how the statistical arrangement of pixels in the image changes with the scale and how much such correlation approximates a fractal behavior. In this sense, our descriptor also measure the self-similarity and complexity of the image but upon a statistical viewpoint. This is a rich and not explored perspective, which is studied in depth in this work.

We use the whole power-law curve of the dimension and apply a time-scale transform to emphasize the multiscale aspect of the features. Finally, we test the proposed method over two well-known datasets, that is, Brodatz and Outex, comparing the results with another fractal descriptor approach showed in BCB09 and other conventional texture analysis methods. The results demonstrat that probability descriptors achieve a more precise classification than other classical techniques.

Ii Fractal Theory

In recent years, fractal geometry concepts have been applied to the solution of a wide range of problems SMS10 ; HWZ08 ; CDHLAB03 ; W08 ; TWZ07 ; LC10 , mainly because conventional Euclidean geometry has severe limitations in providing accurate measures of real-world objects.

ii.1 Fractal Dimension

The first definition of fractal dimension provided in M68 , is the Hausdorff dimension. In this definition, a fractal object is a set of points immersed in a topological space. Thus one can use results from Measure Theory to define a measure over this object. This is the Hausdorff measure expressed by


where denotes the diameter of , that is, the maximum possible distance among any elements of :


Here, a countable collection of sets , with , is a -cover of if .

Notice that also depends on a parameter , which expresses the scale at which the measure is taken. We can eliminate such dependence by applying a limit over , defining in this way the -dimensional Hausdorff measure:


The plot of as a function of shows a similar behavior in any fractal object analyzed. The value of is for any and it is for any , where always is a non-negative real value. is the Hausdorff fractal dimension of . More formally,


In most practical situations, the Hausdorff dimension is difficult or even impossible to calculate. Thus assuming that any fractal object is intrinsically self-similar, the literature shows a simplified version, also known as the similarity dimension or capacity dimension:


where is the number of rules with linear length used to cover the object.

In practice, the above expression may be generalized by considering to be any kind of self-similarity measure and to be any scale parameter. This generalization has given rise many methods for estimating fractal dimension, with widespread applications to the analysis of objects that are not real fractals (mathematically defined) but that present some degree of self-similarity in specific intervals. An example of such a method is the probability dimension, used in this work and described in the following section.

ii.2 Probability Dimension

The probability dimension, also known as the information dimension, is derived from the information function. This function is defined for any situation in which we have an object occupying a physical space. We can divide this space into a grid of squares with side-length and compute the probability of points of the object pertaining to some square of the grid. The probability function is given by


where is the maximum possible number of points of the object inside a unique square. Here we use a generalization of teh above expression defined in the multifractal theory PV02 :


where is any real number.

The dimension itself is given as


When this dimension is estimated over a gray-level digital image , a common approach is to map it onto a three-dimensional surface as


In this case, we construct a three-dimensional grid of 3D cubes also with side-length . The probability is therefore given by the number of grid cubes containing points on the surface divided by the maximum number of points inside a grid cube.

Figure 1: Two estimates of the probability dimension. Above, the 2D version used for shape analysis. Below, the 3D version used for gray-level images. In the 3D case, the original image is represented by a surface.

Iii Fractal Descriptors

Fractal descriptors are values extracted from the relationship common to most methods of estimating fractal dimension. Actually, any fractal dimension method derived from the concept of the Hausdorff dimension obeys a power-law relation, which may be expressed as


where is a measure depending on the fractal dimension method and is the scale at which this measure is taken.

Therefore Fractal descriptors are provided from the function :


We call the independent variable to simplify the notation. Thus and our fractal descriptor function is denoted . For the probability dimension used in this work, we have


The values of may be directly used as decriptors of the analyzed image or may be post-processed by some kind of operation aimed at emphasizing some specifical aspects of that function. Here, we apply a multiscale transform to and obtain a bi-dimensional function , in which the variable is related to and is related to the scale at which the function is observed. A common means of obtaining is through a wavelet transform:


where is a wavelet basis function and is the scale parameterGM84 . Figure 2 shows an example where two textures with the same dimension, but visually distinct, provide different descriptors.

Figure 2: Two textures with the same dimension (FD), but with different descriptors.

Iv Proposed Method

This work proposes to obtain fractal descriptors from textures by using the probability fractal dimension, computing them from the curve in Eq. 8. Empirically, we obtained as the best value of in the Equation 7. Therefore we apply a multiscale transform to .

The multiscale process employs a wavelet transform of , as described in the previous section:


As the multiscale transform maps a one-dimensional signal onto a bi-dimensional function, it is a process that generates intrinsic redundancies. There are different approaches to elliminating such redundancies and keeping only the relevant information CC00 . Here, we adopt a simple method, fine-tuning smoothing, in which is projected onto a specific value of the Gaussian parameter. We tested values of ranging between and and used the values that provided the best performance in the training experiments.

Finally, we selected a specific region from to compose the descriptors. Empirically, we observed that the initial points in this curve provided better performance in our application. Then, we established a threshold after which all points in the convolution curve are disregarded and the values in the curve are taken as the proposed descriptors.

Figure 3: Discrimination of texture image by the proposed descriptors. We have six images from two classes and the respective descriptors. Notice that the classes are substantially separated by the curves.

V Experiments

In order to verify the efficiency of the proposed technique, we applied our probability descriptors to the classification of two benchmark datasets and compared our results to the performance of other well-known and state-of-the-art methods for texture analysis.

The first classification task used the Brodatz dataset, a classic set of natural gray-level textures photographed and assembled in an architecture book B66 . This dataset is composed by 111 classes with 10 textures in each class. Each image has a pixel dimension of 200200.

The second data set was Outex, a set of color textures extracted from natural scenes OMPVKH02 . Here, we used the first 20 classes, each one having 20 images with a 128128 pixel dimension, and converted them to gray-level images.

We compared our probability descriptors to six other techniques, namely, Local Binary Patterns (LBP) PHZA11 , Gabor-wavelets MM96 , Gray-Level Difference Method (GLDM) WDR76 , a multifractal approach described in PV02 and Bouligand-Minkowski fractal descriptors BCB09 ; FB13 .

Therefore we applied a Principal Component Analysis (PCA)


over the data to elliminate or at least attenuate the correlation among the features. Finally we classified each descriptor by a K-fold process, with

, using the Support Vector Machine (SVM) method

V99 and compared the results.

Vi Results

Table 1 shows the correctness rate in the classification of the Brodatz dataset using the compared descriptors. The proposed method obtained the best result, outperforming the powerful Bouligand-Minkowski fractal descriptors and taking substantial advantage over other state-of-the art techniques such as Gabor and LBP. For this result we used and a threshold

. A particularly important aspect of our method with this data set is the reduced number of descriptors needed to provide a precise classification. This point is especially important in large databases, for which computational performance is more relevant. Furthermore, the small number of features avoids the curse of dimensionality, which impairs the reliability of the global result.

Results for the Brodatz dataset Method Correctness Rate (%) Number of descriptors LBP 82.5 0.2 10 Gabor 86.8 0.1 8 GLDM 79.6 0.2 8 Multifractal 54.9 0.5 9 Bouligand-Minkowski 90.8 0.1 10 Proposed method 93.0 0.1 8

Results for the Outex dataset Method Correctness Rate (%) Number of descriptors LBP 98.2 0.0 7 Gabor 98.5 0.0 8 GLDM 90.5 0.1 9 Multifractal 82.0 0.2 10 Bouligand-Minkowski 97.2 0.0 10 Proposed method 99.2 0.0 7

Table 1: Correctness rate for Brodatz and Outex datasets.

Table 1 shows the results for the Outex textures. In this case, we obtained the best result by using and . Again, the proposed approach provided the greatest success rate, despite the challenge of applying a gray-scale-based method to color analysis. In fact, Outex textures exhibit nuances which are better expressed in the color information, such as the changes in the lighting perspective and the images from different classes presenting similarities in the intensity distribution, though distiguished by color. Based on this result, our method demonstrates that although it does not use any color information, it is powerful also for color image analysis.

Figure 4

shows how the success rate varies according to the number of descriptors used in both datasets. The graphs show a well-known property of Karhunen-Loève transform. The most expressive information is concentrated in the initial descriptors, so that the success curves show a quick growing and then tend to stabilize at a constant rate. The larger size and the native gray-scale format of Bradatz data leads to a clearer advantage of probability descriptors in that database. In Outex, the first descriptors, corresponding to the PCA components with higher variance, do not have as much significance for the classification purpose. However, the sum of all of them provide the best result. This is a specific property of fractal descriptors, as can be observed in the Bouligand-Minkowski descriptors for the Brodatz data as well. Fractal descriptors are tightly correlated among themselves, thus we do not have a large significance carried only in a few descriptors.

Figure 4: Success rate against the number of descriptors in each dataset. a) Brodatz. b) Outex.

Finally, Figures 5 and 6 show the confusion matrices of the methods with the best performances. In this kind of representation, a good descriptor must produce a matrix with a diagonal as lighter and continuous as possible and the minimum of dark points outside the diagonal.

As can be seen, in Brodatz data, the probability descriptors clearly presented these characteristics, with almost no “gap” in the diagonal and with a few dark points outside. Both gaps and gray points indicate the confusion of the classifier, that is, elements classified incorrectly in some way. This confusion is caused mostly by the high similarity inter-class and low similarity intra-class. A precise descriptor, like the proposed, avoids such confusion by providing measures capable of faithfully representing the most complex structures.

Figure 5: Confusion matrices in Brodatz dataset. In such figures, we have the predicted classes in the rows and the expected ones in the columns. The number of images expected/predicted in each class is given by the gray-level at each point (lighter points correspond to large number of images). a) Gabor. b) Co-occurrence. c) LBP. d) Proposed method.

In the case of Vistex, the diagonal gaps are not so clear, given the small number of classes. Thus the advantage of the proposed method can be seen in the reduced number of gray squares outside the diagonal. Such squares correspond to the confused classes. We observe that, particularly, the last classes have some discrimination difficulties. The elements of those classes are often assigned to other classes as depicted in the matrices. However, even in these cases, the proposed descriptors showed the expected robustness, assigning the elements correctly.

Figure 6: Confusion matrices in Vistex dataset. a) Gabor. b) Co-occurrence. c) LBP. d) Proposed method.

An overall analysis of the results demonstrates that the proposed method outperformed the compared ones in both datasets, using a small number of descriptors. Such results were expected from fractal theory given its wide applicability to the analysis of natural textures. Actually, fractal geometry presents a remarkable flexibility in the modeling of objects that cannot be well represented by Euclidean rules. The fractal dimension is a powerful metric for the complex patterns and spatial arrangements usually found in nature. Fractal descriptors provide a way of capturing multiscale variations and nuances that could not be measured by conventional methods. More specifically, the probability descriptors proposed here combine a statistical approach with fractal analysis, comprising a framework that supports a precise and reliable discrimination technique, as confirmed in the above results.

Vii Conclusion

We have proposed a novel method for extracting descriptors by applying a multiscale transform over the power-law relation of the fractal dimension estimated by the probability method.

We tested the efficiency of the proposed technique in the classification of two well-known benchmark texture datasets and compared its performance to that of other classical texture analysis methods. The results demonstrated that probability fractal descriptors are a powerful tool for modeling such textures. The proposed method achieved a high success rate in the classification of the benchmark data sets, using fewer than 10 descriptors in this task. These results demonstrate that the proposed method is capable of combining precision, low computational cost and robustness.

As a consequence, our method offers a reliable approach to solve a large class of problems involving the analysis of texture images.


J. B. Florindo gratefully acknowledges the financial support of FAPESP Proc. 2012/19143-3. O. M. Bruno gratefully acknowledges the financial support of CNPq (National Council for Scientific and Technological Development, Brazil) (Grant #308449/2010-0 and #473893/2010-0) and FAPESP (The State of São Paulo Research Foundation) (Grant # 2011/01523-1).


  • (1) E. I. Scarlat, M. Mihailescu, A. Sobetkii, Spatial frequency and fractal complexity in single-to-triple beam holograms, Journal of Optoelectronics and Advanced Materials 12 (1) (2010) 105–109.
  • (2) D. Han, M. Wang, J. Zhou, Fractal analysis of self-mixing speckle signal in velocity sensing, Optics Express 16 (5) (2008) 3204–3211.
  • (3) 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.
  • (4) 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 American-Physical-Society, New Orleans, LA, MAR 10, 2008.
  • (5) L. Tian-Gang, S. Wang, N. Zhao, Fractal Research of Pathological Tissue Images, Computerized Medical Imaging and Graphics 31 (8) (2007) 665–671.
  • (6) S. Lorthois, F. Cassot, Fractal analysis of vascular networks: Insights from morphogenesis, Journal of Theoretical Biology 262 (4) (2010) 614–633.
  • (7) D. Harte, Multifractals: theory and applications, Chapman and Hall/CRC, 2001.
  • (8) 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.
  • (9) 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 (1-2) (2000) 49–82.
  • (10) E. T. M. Manoel, L. da Fontoura Costa, J. Streicher, G. B. Müller, Multiscale fractal characterization of three-dimensional gene expression data, in: SIBGRAPI, IEEE Computer Society, 2002, pp. 269–274.
  • (11) L. da F. Costa, R. M. Cesar, Jr., Shape Analysis and Classification: Theory and Practice, CRC Press, 2000.
  • (12) 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.
  • (13)

    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.

  • (14) 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 Botany-Revue Canadienne de Botanique 83 (3) (2005) 287–301.
  • (15) 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.
  • (16) J. B. Florindo, O. M. Bruno, Texture analysis by multi-resolution fractal descriptors, Expert Systems with Applications 40 (10) (2013) 4022–4028.
  • (17) B. B. Mandelbrot, The Fractal Geometry of Nature, Freeman, 1968.
  • (18)

    T. Parrinello, R. A. Vaughan, Multifractal analysis and feature extraction in satellite imagery, International Journal of Remote Sensing 23 (9) (2002) 1799–1825.

  • (19) A. Grossmann, J. Morlet, Decomposition of hardy functions into square integrable wavelets of constant shape, SIAM Journal on Mathematical Analysis 15 (4) (1984) 723–736.
  • (20) P. Brodatz, Textures: A photographic album for artists and designers, Dover Publications, New York, 1966.
  • (21) T. Ojala, T. Maenpaa, M. Pietikainen, J. Viertola, J. Kyllonen, S. Huovinen, Outex - new framework for empirical evaluation of texture analysis algorithms, in: 16TH International Conference on Pattern Recognition, Vol I, Proceedings, International Conference on Pattern Recognition, IEEE Computer SOC, 10662 Los Vaqueros Circle, PO BOX 3014, Los Alamitos, CA 90720-1264 USA, 2002, pp. 701–706.
  • (22) M. Pietikäinen, A. Hadid, G. Zhao, T. Ahonen, Computer vision using local binary patterns (2011).
  • (23) 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.
  • (24) J. Weszka, C. Dyer, A. Rosenfeld, A comparative study of texture measures for terrain classification, SMC 6 (4) (1976) 269–286.
  • (25) R. O. Duda, P. E. Hart, Pattern Classification and Scene Analysis, Wiley, New York, 1973.
  • (26)

    V. Vapnik, The Nature of Statistical Learning Theory (Information Science and Statistics), 2nd Edition, Springer, New York, 1999.