1 Introduction
Fractal theory has presented a growing interest in many applied areas in the last decades, for instance, in Medicine [1, 23, 19], Physics [16, 39, 35], Computer Science [6, 14, 27], Engineering [36, 22, 37], among many other fields.
Most of these applications employ the fractal dimension to describe objects that should
be classified or simply described in some manner. Fractal dimension provides important information about the object. While in a
mathematical fractal, the dimension measures the “fractality”, in a realworld object, it expresses the spatial occupation of the structure. More practically, this implies that fractal dimension is capable of capturing important physical and visual attributes, like roughness, luminance or the repetition of geometrical patterns.Despite its importance as a powerful descriptor, fractal dimension is still limited in the representation of more complex structures. This is true mainly in the analysis of real world objects, when the level of “fractality” varies along the same object. The literature shows some approaches to solve this problem, such as multifractals [18, 8, 3], multiscale fractal dimension [9, 26] and fractal descriptors [5, 12, 31, 13]. This work is focused on the fractal descriptors solution.
Here, we propose a novel fractal descriptor based on the triangular prism fractal dimension [7] and apply it to the discrimination and retrieval of texture images. In this approach, the fractal dimension is estimated at different scales of observation and a weight parameter is introduced as an exponent in the sum of the areas, changing the influence of each scale in the final result and ensuring a more complete and flexible description of the image.
The proposed technique is tested over two wellknown texture datasets used for benchmark purposes. The results of the classification and retrieval of such data sets are compared to other classical and recent texture analysis methods. The results confirm that the proposed descriptor is a valuable tool for image analysis tasks.
2 Related Works
Texture analysis is a paradigm where the image is described in terms of statistical patterns formed by spatial arrangements of pixel intensities. The first known systematic study on this topic was carried out by Haralick [17] and his cooccurrence matrices. Since then, a large number of methods on texture analysis have been proposed in the literature. Among the most successful approaches one can mention local binary patterns [30], bagoffeatures [34], scaleinvariant feature transform [20], spatial pyramid matching [21], invariants of scattering transforms [32], fast features invariant to rotation and scale of texture [33], and others.
During the last decades, another branch of methods that have presented interesting results in texture analysis, especially on natural images, are those based on fractal geometry, particularly multifractals [38], multiscale fractal dimension [26] and fractal descriptors [5]. In this context, this work proposes the study and application of fractal descriptors based on the estimative of the fractal dimension using a tessellation of triangular prisms [7].
Our proposal has some particular characteristics that distinguish it from other approaches in the literature. First, rather than preselecting preferable regions in the image as in [20, 21], here all pixels and scales are equally important a priori, which simplifies the modelling and interpretation of the texture descriptors. Another difference from methods such as those in [20, 32, 33] is that image invariances are not treated explicitly, although the underlying model and multiscale process ensures that such effect is attenuated in practice. This is confirmed here in the experimental analysis and avoids the use of cumbersome strategies when in many cases invariances are not a critical issue or even when, for example, a rotated texture should be interpreted as a different object. Finally, an important distinction should be done from approaches such as those in [17, 34, 30] where direct relations are established based on te pixel values. Here there is a complete and welldefined physical model behind the statistics extracted from the image, causing it to be more precise in most cases and more robust to deformations usually found in natural structures.
3 Fractal Theory
Fractal geometry has been applied to many diverse areas [6, 39, 19, 36]. This is motivated mainly by the flexibility of fractal theory in modelling natural objects, which usually cannot be precisely represented through conventional Euclidean geometry.
Fractal theory also presents a concise and powerful framework to describe and identify a natural object, based primarily on the fractal dimension concept. Roughly speaking, fractal dimension measures the complexity of a structure. In this case, complexity is related to the property of presenting details at different scales of observation. In this way, fractal dimension is of particular importance because it is strongly related to fundamental physical and visual features of the object, such as roughness, luminance, distribution of colors, and others.
The following section describes in a few words some important aspects of fractal geometry theory and its application to texture analysis.
3.1 Fractal Dimension
Fractal dimension is formally defined as being the HausdorffBesicovitch dimension of a geometrical set of points that composes the fractal object.
In order to define the HausdorffBesicovitch dimension, initially the Hausdorff measure should be defined:
(1) 
where states for the diameter of and is a cover of iff with , for all .
Analyzing the behavior of against , one observes that jumps from to at a particular real nonnegative value of . This value is the HausdorffBesicovitch or fractal dimension of .
Although the above definition is consistent and can be applied to any set of points immersed in the Euclidean space, it shows to be difficult or even impracticable in many situations where the fractal dimension of an object has to be estimated. With the aim of simplifying the computation in such situations, an approximate discrete version of the HausdorffBesicovitch dimension, known as similarity dimension , can be defined by:
(2) 
where is the number of “rulers” with length used to cover the fractal object. Actually, is a metric that can be generalized to a large sort of measures, both in spatial and frequency domain. This gives rise to a lot of methods for estimation of fractal dimension [11], like BouligandMinkowski, boxcounting, Fourier, etc. The triangular prism dimension employed in this work is an example of method derived from the similarity dimension.
3.2 Triangular Dimension
This method for the estimation of the fractal dimension of objects represented in a graylevel image, proposed in [7], is based on the relation between the surface area of a triangular tessellation of the graylevel map and the dimension of the base of each triangle.
The image is divided into a grid of squares with sidelength . For each square, a triangular prism is constructed using the pixel intensities in each corner as the heights and a central point whose height is given by the average of the corner heights. Thus in a graylevel image , let , , and be the pixel intensities delimiting the grid square, such that:
(3) 
The center of the prism has height given by the simple average:
(4) 
Figure 1 depicts a scheme of the prism construction.
The set of prisms composes a tessellation surface, such that the total area of this surface can be computed by using some geometrical procedures (Heron formulas). Thus for each prism, the semiperimeter of each face , , and is given by:
(5) 
where , , , , , , and are the segments as labelled in Figure 1 (c):
(6) 
The total area of each triangular prism is provided by the sum of the areas of the four faces:
(7) 
where the area of each face is given by:
(8) 
The total area of the surface is computed by summing the area of each prism in the grid with step :
(9) 
where is the set of points in the grid with step .
This procedure is repeated for a range of values of and in each step the area is estimated. The fractal dimension is extracted from the loglog relation between and , such that, , where is the slope of a straight line fit to the curve .
4 Fractal Descriptors
Fractal descriptors are an extension of fractal dimension concept. Actually, although fractal dimension is an important descriptor it is still insufficient to represent more complex systems. We can easily observe distinct fractals with the same fractal dimension despite their completely different appearance. Such situation is even more complicated when we deal with objects from the real world, which are not real fractals. In these structures, we find different levels of “fractality” according to the observed scale or even to the spatial region analyzed. In this context we need a tool capable of modelling the object in all its extension.
Fractal descriptors [5, 12, 31, 13] constitute a solution to fill this gap, making possible a richer analysis of fractal characteristics present in the object. Figure 2 illustrates two distinct texture images whose fractal dimensions (FD) are similar. Hence using only the FD estimation is not enough to distinguish the images. For the same images, Figure 2 shows the normalized fractal descriptor curves, which demonstrates to be capable of discriminating the textures in a straightforward manner.
Essentially, the purpose of fractal descriptors is to estimate the fractal dimension under different scales, providing information of different patterns and arrangements present in the structure. A natural candidate to allow this analysis is the powerlaw relation intrinsic to the fractal dimension . From the similarity dimension we can write:
(10) 
where is any measure, related to the spatial or frequency distribution of the object and is the scale parameter. The powerlaw relation may be stated in a quite simple fashion:
(11) 
Fractal descriptors consist in extracting features from the following function:
(12) 
The function may be used in different manners either directly, as in [2], or after a multiscale transform [5]
or by the application of Principal Component Analysis or still Functional Data Analysis as in
[13], among many other possibilities.5 Proposed Method
This work proposes a novel fractal descriptor based on the fractal dimension estimated by triangular prisms. In this case, the fractality function corresponds to the area function . Following the general idea of fractal descriptors, the proposed features are provided by .
Furthermore, to improve the ability of identifying multiscale patterns along the texture, the proposed method introduces an exponent weight to the area sum , which now has the following expression:
(13) 
where is defined as in Equation 7.
For image analysis purposes, a range of values is chosen empirically and the following values are taken into account to compose the descriptors :
(14) 
Finally, since the above expression generates a too much large set of features, a reduction of dimensionality is necessary. Such procedure is carried out by a KarhunenLoève (KL) transform. Thus let
be the feature matrix, that is, a matrix where each line corresponds to the vector
of a sample and each column is a descriptor. The covariance matrix is defined by:(15) 
where is the number of descriptors in , is the column of and is the average of the vector. In the following a second matrix
is defined as having in each column the eigenvectors of
, sorted according to decreasing values of eigenvalues. Finally, the KL transform
of is obtained by(16) 
Now, is the new matrix of features. It has the same dimensions of and each line in contains the final descriptors of the respective sample. Based on the KL transform theory, the first descriptors are the most meaningful for the analysis. Here, the descriptors are considered in this order and, as it will be described in Results section, the number of such descriptors is varied between 1 and a predefined maximum, to find out the best configuration empirically.
Figure 3 summarizes the described steps in a visual diagram.
The proposed method combines two types of information to provide a rich and reliable image descriptor. The first is the complexity, measured by the fractal dimension. Such property is strongly correlated with physical characteristics of the object. The second one is the multiscale analysis that is accomplished by the parameter . Using the curve of ensures that the complexity information is expressed along a range of scales. Finally, the exponent provides an empirical weight for each scale, making the proposed descriptors more flexible to address the description of so diverse objects. All this combination yields a method capable of identifying and discriminating objects even in severe situations as when there is high variability among elements of a same class and/or high similarity among samples of different classes. Figure 4 visually illustrates the discrimination of two classes of textures by the proposed method.
5.1 Motivation
5.1.1 Fractal Modelling
To better understand how and why tringular prisms work on the image, we should look at its mathematical interpretation. The classical analysis of the relation between triangular prisms and fractal theory is that adopted in [7]. It relies on the idea of extending the “walkingdividers” approach [24] to a twodimensional manifold. Here we propose a second interpretation based on the theory of fractional Brownian motion (fBm) [11].
Our first objective is to establish the relation between the prism areas and the pixel gray levels in the image. We illustrate the computation in face . As can be written as a function of we can rewrite and :
(17) 
We can also plug the definition of into the area :
(18) 
Replacing , and with their respective representation in terms of , , and we end up with the following expression:
(19) 
Rearranging terms we can write the above expression as
(20) 
The expression enclosed by the square root contains three squared terms and the leading one is . The points in the image where the intensities are and are separated by a distance . Therefore the distribution of this term is classically related to the fBm process [11].
An fBm process
is a Gaussian nonstationary stochastic process, with mean zero and variance
, whose covariance statistics satisfies(21) 
where is the expected value and is a parameter named Hurst exponent. Another way of defining the same process is to write
(22) 
Seminal works on fractal geometry in images, such as those of Mandelbrot [24] and Pentland [29], already demonstrated that fractal characteristics in images reflect fractality in the physical process originating the pictured object. Pentland even carried out a survey to discover that such relation can be perceptually confirmed by human beings. Since then, a number of studies have been based on this assumption and fBm has been considered the canonical representation for such physical processes giving rise to fractal images.
In this context, the expression (22) has been paramount. If the image obeys the statistics in (22) with an acceptable value the image can be analyzed as a statistical fractal. Furthermore, as demonstrated in [11], the fBm has fractal dimension (in the sense of HausdorffBesicovitch) equals to
with probability
.In practice it is usual to verify the statistics (22) in an image by checking the curve of where stands for the average and and are points separated by . This is in essence what is expressed by the term in (20). The second squared term () can be seen as a correction taking into account the influence of the neighbor pixels in the fBm covariance. Similarly, the parameter appears as a weight expressing the scale size.
We can conclude from this that the triangular prisms are capable of representing two classical viewpoints in fractal geometry, to know, a geometrical interpretation such as “walkingdividers” and a statistical approach derived from fBm. Back to the context of fractal descriptors, the statistics (22) does not to fit perfectly the image. Rather what we expect to find is how the image is more or less close to a fractal in each scale. And this analysis highly depends on the richness of the measure provided by the method employed to estimate the fractal dimension. In this way, a method with a more complete description of the fractal process, both in a geometrical and a statistical sense, tends to be more appropriate than other classical methods in the literature, hence explaining the success achieved by the proposed approach.
5.1.2 Role of exponent
Now we turn our attention to the exponent in (13). Although a generalized version of the multinomial theorem with real exponents exists, it is not straightforward to be interpreted. Therefore we opted for employing MacLaurin series up to the second order. The general expression for a real multivariable function is:
(23) 
We will replace the first term in the summation (13) by and all the others by resulting in the function . Plugging the power function derivatives:
(24) 
Whereas the first summation term is related to the mean of the following ones relate to the correlation between groups of variables. The exponent gives weight to these terms and the higher the exponent the larger the influence of the correlation within larger groups of variables. Correlation is a wellestablished statistics in texture images as it is related to Haralick secondorder statistics [17] and quantify such property in terms of the prism areas and indirectly on the local fractality of the image.
Altogether the combination of an accurate fractal analysis with an intrinsic correlation is what allows for the analysis proposed here to be powerful even in the most challenging scenarios such as when we have high variability among images of a same group or when images from different groups can be confused by more traditional approaches in the literature, where the modelling may be not sufficiently solid and flexible at the same time neither metrics employed can express the parameters of the underlying model with the necessary accuracy.
6 Algorithm
The following pseudocode contains the algorithm for the triangular prism descriptors. Essentially it is not too much than a mere translation of the mathematical expressions in Section 3.2. A simplified version written in Matlab/Octave can also be found at ^{1}^{1}1www….;
7 Experiments
The proposed method is tested over two data sets, that is, Brodatz [4] and Vistex [28]. Brodatz is a graylevel texture database composed by 111 texture images divided each one into 16 windows with 128128 pixels without overlapping. Vistex is composed by 54 color texture images, each one divided into 16 128128 windows. The color images are converted into graylevel ones before the descriptor analysis. In both experiments, each entire image is a class and each window corresponds to a sample.
The classification is performed by a Linear Discriminant Analysis (LDA) method [10]. The classification of each data set is achieved by inputting the proposed descriptors and other wellknown texture descriptors to the classifier, using a 10fold scheme, where one tenth part of the database is used for testing and the remaining samples for training. The classification results are compared with other stateoftheart and classical texture descriptors in the literature, that is, Local Binary Patterns (LBP) [30], Gaborwavelets [25], Multifractal spectrum [38], Gray Level Cooccurrence Matrix (GLCM) [17] and Fourier descriptors [15].
8 Results
8.1 Parameter Settings
The proposed method depends on an exponential parameter that can be tuned to find the most suitable descriptors for each type of image. In this way, to find the best parameters is a fundamental task to ensure a robust solution. In real situations, these variables can be empirically evaluated over the training set. Figure 5(a) shows the precision in the classification of Brodatz data set when using exponents ranging between 1 and 2. Outside this range there is no significant gain. Based on this plot, the triangular descriptors with exponents between 1 and 1.8 were concatenated and summarized through PCA to provide the texture descriptors used in the following experiments. Ten descriptors are used for each exponent.
(a)  (b) 
A similar test was accomplished over Vistex database giving rise to the curve in Figure 5 (b). In this case, exponents with values between 0.4 and 1.2 were selected to compose the proposed descriptors.
8.2 Image Retrieval
Now the proposed method is applied to practical problems in image analysis and its performance is compared to other texture descriptors proposed in the literature. The first task discussed is image retrieval, when the user inputs a sample and tries to recall similar images (from the same class). Each attempt to recall the expected sample is named “query” and the aim is to obtain the maximum ratio of expected samples with the minimum possible of queries. A complete description of this tradeoff is given by the precision/recall curve. Precision is the ratio between the number of expected samples returned and the number of queries. Recall is the ratio between the number of right guesses and the total number of expected samples in the database.
Originally, precision/recall is defined only for twoclasses problems. As there are more classes here, the curve is computed for each class, assuming the current class as the positive prediction and all the remaining as the negative prediction. Thus an average precision/recall curve is presented for discussion. Figure 6 shows the curves for Brodatz textures. A good method is supposed to have a curve the closest possible of a constant curve with precision always being 1. On average, this profile is satisfied by the proposed descriptors, once although it has not so high precision for low recall, it outperforms the other approaches for higher recall. Such behavior implies that the triangular descriptors do better in more complicated problems where more queries are required to obtain the correct result.
To quantify the performance on image retrieval a global measure may also be necessary and the most commonly used is the area under the curve of precision/recall. Table 1 compares this measure for different methods on Brodatz textures. The proposed method has the largest computed area confirming its efficiency for image retrieval. Actually, the poor result for lower recalls is superimposed by the great precision when the number of queries increases.
Method  Area Under Curve 

LBP  0.84 
GLCM  0.82 
Multifractal  0.81 
Gabor  0.81 
Fourier  0.74 
Proposed  0.90 
Figure 7 shows the precision/recall curves for Vistex database. Here, a similar phenomenon to that found in Figure 6 is repeated and the triangular descriptors have advantage for recalls close to 1. Now the precision for recalls close to 0 is acceptable yet it is not better than stateoftheart methods such as LBP and Multifractal Spectrum. In both cases, the proposed method behaves more regularly than the other approaches, making it a potential candidate to be applied in practical situations of image retrieval problems.
The areas under curves are showed in Table 2. Again, the triangular fractal descriptors have the largest area and have great potential to provide good outcomes in this kind of task.
Method  Area Under Curve 

LBP  0.82 
GLCM  0.79 
Multifractal  0.79 
Gabor  0.81 
Fourier  0.80 
Proposed  0.95 
Generally speaking, in both data sets, the proposed method presents better results than their counterparts after a number of retrieval attempts. This is a consequence of combining different weights given by the exponents to each scale in the triangular descriptors. Although the empirical combination may not be enough for an initial guess, it ensures more robustness when more queries are requested.
8.3 Image Classification
In this section, the textures from the benchmark sets are classified. Since the actual classes are known previously, the accuracy of each approach can be measured and compared by statistical metrics.
Figure 8 shows in a plot the relation between the number of descriptors employed in the classification and the rate of images correctly classified (success rate) obtained in Brodatz data set. Even though the stateoftheart LBP method is the best when using a few descriptors, the triangular fractal descriptors ensures a more accurate classification using more than 34 features, which is a reasonable number given that the error is kept within an acceptable range.
Table 3 shows the success rates, the respective errors and number of descriptors for each compared texture descriptors. The proposed descriptors achieved an advantage of more than 2% over the second best approach (LBP). Such gain is relevant given the size and number of classes in this database, which makes such rates close to the limit of best possible results. Moreover, if on the one hand more descriptors are necessary, on the other the error in the cross validation is lower, showing that no dimensionality curse is detected.
Method  Success rate  Number of descriptors 

LBP  87.330.02  15 
GLCM  86.480.02  70 
Multifractal  85.640.03  70 
Gabor  85.420.02  19 
Fourier  78.710.03  15 
Proposed  89.580.01  85 
Figure 9 shows the confusion matrices for the respective descriptors. Such matrices expressing the number of images assigned to each class is represented in a grayscale image. In these diagrams, an ideal method should present a black diagonal with a clear white background. Any light point in the diagonal and gray point outside corresponds to errors in the classification. The primary goal of this type of expression is to show the behavior of the descriptors in each class. For instance, both triangular and LBP methods come off badly in the classes 42/43. Those classes do not have enough texture information, but they are more suitable for a shapebased analysis. On the other hand, they completely differ in classes like 69 and 62. Images from class 69 exhibit different macropatterns that impair the fractalbased method, but not local approaches such as LBP. The opposite situation arises with class 62, where there are local illumination changes but the multiscale patterns are preserved.
Figure 10 shows the relation success rate/number of descriptors for Vistex. The general shape of each curve is similar to those in Figure 8, confirming the reproducibility of the experiment. The highest success rates are caused by the smaller number of classes, even though each Vistex image is more challenging than Brodatz, as it has a higher level of variances intraclasses.
Table 4 shows the success rates for Vistex. It is worth to notice that a number of triangular descriptors smaller than in Brodatz case are necessary to obtain the best result. This is also a consequence of the reduced number of classes, which makes the problem statistically simpler.
Method  Success rate (%)  Number of descriptors 

LBP  91.550.03  13 
GLCM  88.210.03  70 
Multifractal  88.310.03  76 
Gabor  90.390.01  17 
Fourier  84.490.02  15 
Proposed  92.940.01  50 
Figure 11 shows the confusion matrices for Vistex database. The best three descriptors (Gabor, LBP and Triangular) have matrices visually similar, but the behavior in each class differs among the approaches. For example, the class 34 has different global patterns that can be confused by methods like Gabor and Triangular, while the class 5 has large homogeneous regions that impair the performance of local methods like LBP.
8.4 Invariances
Another important test to assess the robustness of any image descriptor is to verify its invariance to changes like rotation and addition of noise.
The first evaluated invariance is to noise following Gaussian distribution, a situation commonly found in practice. The levels of noise (ratio of affected pixels) ranges between 0.01 and 0.05. Figure
12 shows examples of such images affected by the different levels of noise in our tests.Table 5 presents results comparing with other descriptors in Brodatz data set. A similar behavior is expected for other textures. Although Gabor obtained the best results, mainly for more severe noises, the proposed method does better than the stateoftheart LBP method and for a small amount of noise it has a performance similar to Gabor.
Method  Noise ratio  

0.01  0.02  0.03  0.04  0.05  
LBP  71.23  69.88  65.14  62.39  59.79 
GLCM  80.46  78.10  78.21  78.04  77.65 
Multifractal  76.58  71.00  66.10  65.48  63.79 
Gabor  83.67  83.73  83.22  82.77  82.38 
Fourier  78.99  78.71  78.43  78.15  77.59 
Proposed  83.33  81.64  79.05  77.36  76.13 
A second invariance experiment evaluates textures rotated by predefined angles. To ensure the same dimension for all textures, a procedure to extract only the central region of each image is accomplished, as illustrated in Figure 13. In this way, if the original texture image has dimensions , the extracted region after the rotation is the central part with dimension .
Figure 14 shows the example of a texture rotated by the angles considered in this experiment, that is, 0, 30, 60 90, 120, 150 and 180 degrees.
Table 6 presents the correctness rates for the classification of Brodatz database after the rotation of the textures. The lowest rates are caused by the use of a cropped region of the original texture, leading to an obvious loss of information. The proposed approach outperforms the other methods with higher success rates and a low variability in the results when the angle is changed. The maximum difference is only a little over 1%, while LBP, for example, shows a difference of almost 10% in the result when the textures are rotated by 150 and 180 degrees.
Method  Rotation angle (degrees)  

0  30  60  90  120  150  180  
LBP  81.98  73.70  72.86  81.98  73.98  72.75  82.15 
GLCM  77.70  75.28  75.50  77.93  75.50  76.23  77.53 
Multifractal  76.63  76.24  73.99  77.25  76.12  76.18  76.63 
Gabor  80.13  80.41  79.06  80.01  80.13  79.11  80.29 
Fourier  70.32  72.07  71.56  69.98  72.07  71.78  70.21 
Proposed  83.33  84.69  84.18  83.67  84.58  83.84  83.56 
8.4.1 Discussion
The proposed method differs from other approaches to texture analysis in the literature in that the physical process underlying the generation of the image is taken into account to a high extent, as we demonstrated its relation with fractional Brownian processes. These processes are known to be widely present in nature and its association with the way that our brain perceives materials ans scenarios around us is well established. Parallel to this, the triangular prisms also encompass a geometrical analysis similar to “walkingdividers”, intimately related to the primitive idea of fractals and also commonly observed in natural structures. In itself, such complementary and dual viewpoint of fractality is an important contribution and novelty with regards to other fractal approaches such as multifractals [38], BouligandMinkowski [2], and others.
The consideration of the physical process also showed to be more advantageous than simply quantifying relations between neighbor pixels without accounting for the semantic involved in the represented object, as in LBP and GLCM descriptors. Even though the interpretation of these approaches can be considered more straightforward, the lack of a more realistic model makes them insufficient in more complex cases with larger and more heterogeneous databases as those presented here.
Finally, we also compare to other methods where the image is observed beyond the simple pixel values, such as in Fourier and Gabor descriptors. These descriptors, however, are ultimately based on linear filtering and are not the ideal solution to explain the nonlinearities present in many natural images, mostly caused by the chaotic behavior associated to fBm processes.
Generally speaking, the great performance of triangular fractal descriptors in the presented results confirm what is expected from the fractal analysis of natural texture images. Indeed, fractal geometry constitutes a reliable tool to model those structures, which cannot be welldescribed in an Euclidean framework. This analysis is especially complete here as two fractal viewpoints (statistical and geometrical) are merged to provide an even more accurate description of the image. Furthermore, fractal descriptors approach enhances the conventional fractal geometry analysis by extracting relevant information of complexity under different scales, by the multiscale analysis implicit when the dimension is estimated for different values of . The combination of a solid modeling with a multiscale context also explains the robustness to noise and rotation, as noise is more active in local scales and deformations like the image rotation do not alter model parameters at a substantial level. In this way, these descriptors are capable of providing features representing precise and rich measures of the spatial structure of complex textures like those in benchmark data, as used here, as well as in practical application involving the identification and discrimination of intricate patterns in objects and scenarios represented in a digital image.
9 Conclusion
This work proposed a novel graylevel texture descriptor based on a particular method of fractal dimension estimation, named triangular prism method. The descriptors were obtained by combining values of area of a triangular tessellation of the texture, using different steps for the tessellation grid as well as different values for an exponential parameter used as a weight for the area within each cell in the grid.
Our study demonstrated that the triangular prisms are capable of extracting fractal characteristics of the image under two complementary perspectives: geometrical (walkingdividers) and statistical (fractional Brownian motion).
These descriptors were compared to classical and stateoftheart texture descriptors in a task of classification of two wellknown benchmark texture datasets, e.g., Brodatz and Vistex.
The classification results illustrated that the novel approach is a valuable descriptor, achieving remarkable results in image classification and retrieval tasks and suggesting its application to a wide range of problems comprising tasks of pattern recognition in digital images.
Acknowledgements
O. M. Bruno gratefully acknowledges the financial support of CNPq (National Council for Scientific and Technological Development, Brazil) (Grant #307797/20147 and Grant #484312/20138) and FAPESP (The State of São Paulo Research Foundation) (Grant # 14/080261). J. B. Florindo gratefully acknowledges the financial support of FAPESP Proc. 2013/222053 and 2012/191433.
References
 [1] M. Z. C. Azemin, D. K. Kumar, T. Y. Wong, R. Kawasaki, P. Mitchell, and J. J. Wang. Robust Methodology for Fractal Analysis of the Retinal Vasculature. IEEE Transactions on Medical Imaging, 30(2):243–250, FEB 2011.

[2]
A. R. Backes, D. Casanova, and O. M. Bruno.
Plant leaf identification based on volumetric fractal dimension.
International Journal of Pattern Recognition and Artificial Intelligence (IJPRAI)
, 23(6):1145–1160, 2009.  [3] R. Benzi, G. Paladin, G. Parisi, and A. Vulpiani. On the Multifractal Nature of Fullydeveloped Turbulence and Chaotic Systems. Journal of Physics AMathematical and General, 17(18):3521–3531, 1984.
 [4] P. Brodatz. Textures: A photographic album for artists and designers. Dover Publications, New York, 1966.
 [5] O. M. Bruno, R. de Oliveira Plotze, M. Falvo, and M. de Castro. Fractal dimension applied to plant identification. Information Sciences, 178(12):2722–2733, 2008.

[6]
H.C. Chen, F.C. Gu, and M.H. Wang.
A novel extension neural network based partial discharge pattern recognition method for highvoltage power apparatus.
Expert Systems with Applications, 39(3):3423–3431, FEB 15 2012.  [7] K. C. Clarke. Computation of the fractal dimension of topographic surfaces using the triangular prism surface area method. Comput. Geosci., 12:713–722, September 1985.
 [8] L. V. Constantin and D. A. Iordache. Study of the Fractal and Multifractal Scaling Intervening in the Description of Fracture Experimental Data Reported by the Classical Work: Nature 308, 721722(1984). Mathematical Problems in Engineering, 2012.
 [9] L. da F. Costa and R. M. Cesar, Jr. Shape Analysis and Classification: Theory and Practice. CRC Press, December 2000.
 [10] R. O. Duda and P. E. Hart. Pattern Classification and Scene Analysis. Wiley, New York, 1973.
 [11] K. J. Falconer. The Geometry of Fractal Sets. Cambridge University Press, New York, NY, USA, 1986.
 [12] J. B. Florindo and O. M. Bruno. Fractal descriptors in the fourier domain applied to color texture analysis. Chaos, 21(4):043112, 2011.
 [13] J. B. Florindo, M. De Castro, and O. M. Bruno. Enhancing Multiscale Fractal Descriptors Using Functional Data Analysis. International Journal of Bifurcation and Chaos, 20(11):3443–3460, NOV 2010.
 [14] K. Gdawiec and D. Domanska. Partitioned Iterated Function Systems with Division and a Fractal Dependence Graph in Recognition of 2D Shapes. International Journal of Applied Mathematics and Computer Science, 21(4):757–767, DEC 2011.
 [15] R. C. Gonzalez and R. E. Woods. Digital Image Processing (2nd Edition). Prentice Hall, Upper Saddle River, N.J., January 2002.

[16]
Q. Guo, J. Guo, Z. Liu, and S. Liu.
An adaptive watermarking using fractal dimension based on random fractional Fourier transform.
Optics and Laser Technology, 44(1):124–129, FEB 2012.  [17] R. M. Haralick. Statistical and structural approaches to texture. Proceedings of the IEEE, 67(5):786–804, 1979.
 [18] D. Harte. Multifractals: theory and applications. Chapman and Hall/CRC, 2001.
 [19] P.W. Huang and C.H. Lee. Automatic Classification for Pathological Prostate Images Based on Fractal Analysis. IEEE Transactions on Medical Imaging, 28(7):1037–1050, JUL 2009.
 [20] S. Lazebnik, C. Schmid, and J. Ponce. A sparse texture representation using local affine regions. IEEE Transactions on Pattern Analysis and Machine Intelligence, 27:1265–1278, 2005.

[21]
S. Lazebnik, C. Schmid, and J. Ponce.
Beyond bags of features: Spatial pyramid matching for recognizing
natural scene categories.
In
Proceedings of the 2006 IEEE Computer Society Conference on Computer Vision and Pattern Recognition  Volume 2
, CVPR ’06, pages 2169–2178, Washington, DC, USA, 2006. IEEE Computer Society.  [22] C. Liu, R. L. Panetta, and P. Yang. The Influence of Water Coating on the Optical Scattering Properties of Fractal Soot Aggregates. Aerosol Science and Technology, 46(1):31–43, 2012.
 [23] R. Lopes and N. Betrouni. Fractal and multifractal analysis: A review. Medical Image Analysis, 13(4):634–649, AUG 2009.
 [24] B. B. Mandelbrot. The Fractal Geometry of Nature. Freeman, 1982.
 [25] B. Manjunath and W. Ma. Texture features for browsing and retrieval of image data. IEEE Transactions on Pattern Analysis and Machine Intelligence, 18:837–842, 1996.
 [26] E. T. M. Manoel, L. da Fontoura Costa, J. Streicher, and G. B. Müller. Multiscale fractal characterization of threedimensional gene expression data. In SIBGRAPI, pages 269–274. IEEE Computer Society, 2002.
 [27] G. Min, J. Hu, and M. E. Woodward. Performance Modelling and Analysis of the TXOP Scheme in Wireless Multimedia Networks with Heterogeneous Stations. IEEE Transactions on Wireless Communications, 10(12):4130–4139, DEC 2011.
 [28] MIT. Mit vistex texture database, 2011.
 [29] A. P. Pentland. Fractalbased description of natural scenes. PAMI6(6):661–674, 1984.
 [30] M. Pietikäinen, A. Hadid, G. Zhao, and T. Ahonen. Computer vision using local binary patterns, 2011.
 [31] R. O. Plotze, J. G. Padua, M. Falvo, M. L. C. Vieira, G. C. X. Oliveira, and 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):287–301, 2005.
 [32] L. Sifre and S. Mallat. Rotation, scaling and deformation invariant scattering for texture discrimination. In CVPR, pages 1233–1240. IEEE Computer Society, 2013.
 [33] M. Sulc and J. Matas. Fast Features Invariant to Rotation and Scale of Texture, pages 47–62. Springer International Publishing, Cham, 2015.
 [34] M. Varma and A. Zisserman. A statistical approach to texture classification from single images. International Journal of Computer Vision, 62(12):61–81, 2005.
 [35] B. B. Wang, G. B. Dong, and X. Z. Xu. Carbon fractals grown from carbon nanotips by plasmaenhanced hot filament chemical vapor deposition. Applied Surface Science, 258(5):1677–1681, DEC 15 2011.
 [36] Y. Wu, Q. Lin, Z. Chen, W. Wu, and H. Xiao. Fractal analysis of the retrogradation of rice starch by digital image processing. Journal of Food Engineering, 109(1):182–187, MAR 2012.
 [37] H. P. Xie, J. F. Liu, Y. Ju, J. Li, and L. Z. Xie. Fractal property of spatial distribution of acoustic emissions during the failure process of bedded rock salt. International Journal of Rock Mechanics and Mining Sciences, 48(8):1344–1351, DEC 2011.
 [38] Y. Xu, H. Ji, and C. Fermüller. Viewpoint invariant texture description using fractal analysis. Int. J. Comput. Vision, 83(1):85–100, June 2009.
 [39] L. Yaomin, L. Zhongliang, and H. Lingyan. Experimental and theoretical investigations of the fractal characteristics of frost crystals during frost formation process. Experimental Thermal and Fluid Science, 36:217–223, JAN 2012.