 Methodology
 Open access
 Published:
Twodimensional multifractal detrended fluctuation analysis for plant identification
Plant Methods volume 11, Article number: 12 (2015)
Abstract
Background
In this paper, a novel method is proposed to identify plant species by using the two dimensional multifractal detrended fluctuation analysis (2D MFDFA). Our method involves calculating a set of multifractal parameters that characterize the texture features of each plant leaf image. An index, I _{0}, that characterizes the relation of the intraspecies variances and interspecies variances is introduced. This index is used to select three multifractal parameters for the identification process. The procedure is applied to the Swedish leaf data set containing leaves from fifteen different tree species.
Results
The chosen three parameters form a threedimensional space in which the samples from the same species can be clustered together and be separated from other species. Support vector machines and kernel methods are employed to assess the identification accuracy. The resulting averaged discriminant accuracy reaches 98.4% for every two species by the 10 − fold cross validation, while the accuracy reaches 93.96% for all fifteen species.
Conclusions
Our method, based on the 2D MFDFA, provides a feasible and efficient procedure to identify plant species.
Introduction
The increasing interest in biodiversity and biocomplexity, together with the growing availability of digital images and image analysis algorithms, makes plant species identification and classification a topic that has attracted many researchers’ attention. In general, many parts of a plant such as flowers, seeds, roots, and leaves can be used to identify plant species [13]. In this paper, we focus on the usage of image of leaves as they are widely available. Leaf’s shape, color, vein properties, texture and contours are important features for plant identification. For example, leaf shapes were used in [46]; complex veins and contours of leaves were used in [7] and leaf texture was used in [811] for plant species identification. For plant species identification using digital morphometrics, we refer the reader to [1214] and the references therein.
Note that in [7], a monofractal method was used to extract plant leaf’s features from leaf images. This method was then used in [15,16]. It’s been recognized that the monofractal method cannot fully extract detailed information from the leaf image and therefore cannot be efficiently applied to process the images of the objects that are locally irregular [17]. To overcome this difficulty, several multifractal analysis (MFA) methods were proposed [1822]. For example, Backes et al. [18,19] used multiscale fractal dimensions to describe the texture property of leaf’s surface to identify plants, which turned out to be very efficient. Note that the classical MFA is based on capacity measurement or probability measurement and thus describes only stationary measurements [17]. For a leaf image, the surface itself is hardly stationary. Therefore, the multifractal detrended fluctuation analysis (MFDFA) method that can deal with nonstationary is a desirable method for leaf image analysis [23]. Though the MFDFA method has been successfully applied in many fields for nonstationary series and surfaces [2430], to the best of our knowledge, no work yet has applied the MFDFA on leaf images for plant identification and classification. In this paper, we attempt to identify plant species via leaf images by using the MFDFA. More precisely, we first adopt the MFDFA to extract important texture features from leaf images and obtain several key multifractal parameters, and then we apply the support vector machines and kernel methods (SVMKM) to distinguish leaves from different plant species. The widely used Swedish leaf data set [31] containing leaves from fifteen different Swedish tree species are used for our experiments. Our results show that the average accuracy is 98.4% for every two species by the 10 − fold cross validation; for the overall species, the average accuracy reaches 93.96% by the same validation criterion.
We organize the rest of this paper as follows: in Methods and materials we adopt the twodimensional (2D) MFDFA to calculate the multifractal parameters. In Results and discussion, we present and discuss our results. Our method is then further tested in Model test. A summary is provided in Conclusions.
Methods and materials
Multifractal detrended fluctuation analysis
We first adopt the 2D MFDFA method proposed in [32] to our setting as follows:

Step 1: Regard a leaf image as a selfsimilar surface and represent it by an M × N matrix X = (X(i, j)), i = 1, 2,…, M and j = 1, 2,…, N. Partition the surface into M _{ s } × N _{ s } nonoverlapping square subsurface of equal length s, where M _{ s } ≡ [M / s] and N _{ s } ≡ [N / s] are positive integers (Here [u] stands for the largest integer that is less than or equal to u). Each subsurface is denoted by X _{ m,n } = X _{ m,n }(i, j) with X _{ m,n }(i, j) = X(r + i, t + j) for 1 ≤ i, j ≤ s, where r = (m1)s and t = (n1)s. Note that M and N are not necessarily multiples of the length s, therefore, the subsurfaces in the upperright and the bottom may not be taken into consideration. We can then repeat the partitioning procedure starting from the other three corners.

Step 2: For each subdomain X _{ m,n }, find its cumulative sum
$$ {G}_{m,n}\left(i,j\right)={\displaystyle \sum_{k_1=1}^i{\displaystyle {\sum_{k_2=1}^j}_{m,n}Xm\left({k}_1,{k}_2\right)}}, $$(1) 
where 1 ≤ i, j ≤ s, m = 1, 2, …, M _{ s } and n = 1, 2, …, N _{ s }. Then G _{ m,n } = G _{ m,n }(i, j) (i, j = 1, 2, · · ·, s) itself is a surface.

Step 3: For each surface G _{ m,n }, obtain a local trend G ^{~} _{ m,n } by fitting it with a prechosen bivariate polynomial function. In this paper, we choose the trending function as
$$ {\tilde{G}}_{m,n}\left(i,j\right)= ai+bj+c, $$(2) 
a.where 1 ≤ i, j ≤ s and a, b and c are free parameters to be determined by the leastsquares method. The residual matrix is then given by y _{ m,n } = y _{ m,n }(i, j) with
$$ {y}_{m,n}\left(i,j\right)={G}_{m,n}\left(i,j\right){\tilde{G}}_{m,n}\left(i,j\right). $$(3) 
Step 4: Define the detrended fluctuation function F(m, n, s) for the segment X _{ m,n } as follows:
$$ {F}^2\left(m,n,s\right)=\frac{1}{s^2}{\displaystyle \sum_{i=1}^s{\displaystyle \sum_{j=1}^s{y}_{m,n}{\left(i,j\right)}^2}} $$(4) 
and the qthorder fluctuation function
$$ {F}_q(s)={\left[\frac{1}{M_s{N}_s}{\displaystyle \sum_{m=1}^{M_s}{\displaystyle \sum_{n=1}^{N_s}{\left[F\left(m,n,s\right)\right]}^q}}\right]}^{1/q},q\ne 0. $$(5)$$ {F}_q(s)= \exp \left\{\frac{1}{M_s{N}_s}{\displaystyle \sum_{m=1}^{M_s}{\displaystyle \sum_{n=1}^{N_s} \ln \left[F\left(m,n,s\right)\right]}}\right\},\ q=0. $$(6) 
Step 5: Vary the value of s ranging from 6 to min(M, N)/4. If there is longrange powerlaw correlation for large values of s, then
$$ {F}_q(s)\propto {s}^{h(q)}. $$
This allows us to obtain the scaling exponent h(q) via linearly regressing lnF _{ q }(s) on lns. Note that h(2) is the so called Hurst index of the surface, we then call h(q) the generalized Hurst index of the surface. For each q, the corresponding classical multifractal scaling exponent τ(q) is given by:
where D _{ f } is the fractal dimension of the geometric support of the multifractal measure, and takes the value of D _{ f } = 2 in our work. The generalized multifractal dimension D _{ q } is then given by
In the case where q = 1, D _{1} can be obtained via a linear regression of \( {\displaystyle {\sum}_{m=1}^{M_s}{\displaystyle {\sum}_{n=1}^{N_s}{P}_{m,n} \ln {P}_{m;n}}} \) against lns, where
The other two indicators characterizing the singularity strength of the multifractal surface are the Hölder exponent α(q) and the singularity spectrum f (α), which are given by
Here α(q) characterizes the local singularity of an image texture, and f (α) measures the global singularity of an image texture. Varying the value of q in the range from −15 to 15 determines ∆α and ∆f as follows:
α _{max} = max{α(q), q∈[−15,15]} and α _{min} = min{α(q), q∈[−15,15]}. Note that the index ∆α is considered as an indicator to measure the absolute magnitude of the gray scale volatility. The larger value of ∆α, the smaller even distribution of probability measure and the more roughness image surface will be expected. The index ∆f is the Hausdorff dimension of the measure object, which measures the degree of confusion. Therefore both ∆α and ∆f are important multifractal parameters in describing the characteristics of an image in our study.
Experiment materials
To demonstrate our method of identifying plant species by using the leaf texture, we use the Swedish leaf data set [31] for our experiment, which is widely employed in computer vision and pattern recognition fields [4,33,34], plant taxon fields [1] and image processing fields [6,35]. This leaf data set has images of 15 species of leaves with 75 sample images per species. We label the fifteen species by MI, MII, · · ·, MXV (See Figure 1).
We first transform the color image to gray scale so that each image can be viewed as a three dimensional surface with the first two coordinates (i, j) denoting the 2D position and the third coordinate z denoting the gray level of the corresponding pixel.
Multifractal nature of image surfaces
Each image is stored as a 2D matrix in 256 grey levels. This allows us to follow the procedure introduced in Multifractal detrended fluctuation analysis to calculate the associated h(q) and τ(q). If τ(q) is nonlinear in q, that is h(q) is not independent of q, then the image possesses the multifractal nature.
For the Swedish leaf data set, we find that the leaf images all possess the multifractal nature. Figure 2 and Figure 3 demonstrate the multifractal nature of two randomly chosen leaf images, namely, image MIV004 and image MX017, the former has 1793 × 979 pixels and the latter has 2934 × 1771 pixels. In each the left panel illustrates the dependence of the detrended fluctuation function F _{ q }(s) as a function of the scale s for different q. The well fitted straight lines indicate the evident power law scaling of F _{ q }(s) versus s. The right panel shows that τ(q) is nonlinear in q, indicated by the fact that h(q) depends on q.
Results and discussion
For each image, we can calculate the generalized Hurst exponents h(q) and six other multifractal parameters including α _{max}, α _{min}, ∆α, ∆f, D _{1} and D _{2}. For each tree species, we take the averaged value over the 75 samples and report our calculated values in Figures 4 and 5. Their standard deviations are given in Figures 6 and 7, respectively.
As seen in Figure 4, comparing with h(2) and h(3), the estimations of h(−3), h(−2), h(−1) and h(1) vary in relatively wider dynamic ranges and thus demonstrate better abilities to distinguish textures among different species. Yet, one notes that there are relatively large variations in the standard deviations among the 75 samples for the h(q) exponents in Figure 5. This suggests that this indicator alone may not be adequate to identify the fifteen tree species. Also as seen in Figure 6 that the three parameters, α _{max}, ∆α, and ∆f admit wider dynamic ranges than the other three parameters do. The variations among the 75 samples in the same tree species are notably large as shown in Figure 7.
For species i (i = I, II, · · ·, XV), with respect to each calculated multifractal parameter, we denote the standard deviation of the 75 samples by σ _{ in }(i) and define σ _{ in } as
which represents the intraspecies variance. Note also that for each indicator, we can calculate its value corresponding to each species and there are 15 values in total for those 15 species. We define σ _{ bet.} as the standard deviation of these 15 calculated values. Then the term σ _{ bet.} represents the interspecies variance for each multifractal indicator. We now define an index, I _{0}, as
From the definition, we note that the multifractal parameter with larger I _{0} serves better as an indicator to distinguish species. We present the calculated values of I _{0} in Table 1.
We choose the combination of three multifractal parameters with larger I _{0} values, namely, {h(−3), α _{min}, Δα}, as the feature descriptors for our classification purpose and apply the support vector machines and kernel methods (SVMKM) with the heavytailed radial basis function’htrfb’ as the kernel [36]. It is worth mentioning that the combination of 4 or more parameters does not lead to significant higher accuracies, but at a cost with much longer computational time and with no visual advantages. In this sense, the combination of the above three parameters is optimal. For the total sample set containing 75 × 15 = 1125 samples, we use the K − fold cross validation to evaluate the learning performance. This means that 100 (K − 1)/K% samples are randomly chosen as a training set and the remaining 100/K% samples are considered as a test set. The calculation process is then repeated 10 times to eliminate the impact of randomness.
In our first identification experiment, we test the proposed method through examining the distinguishing effect for every two species. To this end, we form a threedimensional parameter space with components given by the above chosen feature descriptors {h(−3), α _{min}, Δα}. In this space, one point represents a leaf sample image. In Figure 8(a)(d), we plot the corresponding points for Ulmus carpinifolia versus Alnus incana, Salix aurita versus Salix alba Sericea, Salix sinerea versus Tilia and Sorbus aucuparia versus Fagus silvatica, respectively. As shown in these plots, the samples from the same tree species are clustered together reasonably well.
In addition, we calculate the discriminant accuracies of every two tree species by SVMKM using the K − fold cross validation with different K values. The average accuracies of 10 trials are shown in Figure 9(a). To display the applicability of identifying different tree species by our proposed method, as an example, we plot the accuracy of identifying species MI (Ulmus carpinifolia) versus other 14 species with K = 10 in Figure 9(b). As expected, the average accuracy of every two species is increasing with respect to K. The obtained best accuracy is 98.40%, higher than 96.82% reported in [35], which requires a very complex preprocessing process for leaf images. It is seen from Figure 9(b) that there are accuracy variations between species Ulmus carpinifolia and the other 14 species. Five species, namely, Salix aurita, Betula pubescens, Ulmus glabra, Salix sinerea and Fagus silvatica, have accuracies below the average accuracy. This suggests that species Ulmus carpinifolia has high similarity with the above mentioned five species, which agrees with the observation from Figure 1.
For each species, the averaged {h(−3), α _{min}, Δα} of the 75 samples is represented by a single point in the threedimensional parameter space (see Figure 10) in which different points representing different species may be clustered into several groups. We use the calculated discriminant accuracy of every two species as the distance between these two points (species). This allows us to conduct a cluster analysis for all samples of the 15 species by the method of hierarchical clustering [37]. The result is given in Figure 10(b), which suggests that the 15 tree species’ leaf samples can be clustered into five groups: (i) {Ulmus carpinifolia, Salix aurita, Ulmus glabra, Salix sinerea, Fagus silvatica}; (ii) {Betula pubescens, Populus, Sorbus intermedia}; (iii) {Quercus, Alnus incana, Salix alba Sericea, Populus tremula}; (iv) {Acer, Tilia} and (v) {Sorbus aucuparia}. This is consistent with visualizing the images directly from Figure 1 showing our proposed approach is applicable.
As another important aspect of identification experiment, we next test our method through calculating the identification accuracies for different numbers of species. The averaged accuracy result calculated when K = 10 is shown in Figure 11(a). Note that the average accuracy is decreasing as the number of tree species increases. This is due to the increasing probability of incorrect classification. However, under the worst situation, all 75 × 15 = 1125 sample leaf images are well mixed together, which gives the lowest average accuracy: 93.96%. This is still very convincing that our approach is feasible. We calculate the identification accuracy also when K = 10 for each species and report the result in Figure 11(b), while the identification result for each species is displayed in Table 2. The best three accuracies reach 98.67%, 97.33% and 96%, and the corresponding species are Sorbus aucuparia, Sorbus intermedia and Tilia. As is seen in Figure 1, these three species are clearly distinct from the other species in leaf shapes and textures. This again shows that our method is effective and feasible.
We remark that the sample size of each species has little effect on the average discriminant accuracy. To justify this, we randomly choose n (n ≤ 75) leaf samples for each species and run the procedure. Then repeat the process 10 times and take the average accuracy, which is reported in Figure 12. It can be seen from Figure 12 that as the number of samples changes from 40 to 75, the accuracy changes only 0.73%.
Model test
In this section, we test our proposed method to demonstrate its efficiency. More precisely, we test the validity of the optimal multifractal parameter combination {h(−3), α _{min}, Δα}. To this end, we choose other four combinations composed by three multifractal parameters to construct four threedimensional spaces from Table 1. These four choices are {h(−3), Δf, D _{1}}, {h(2), h(3), α _{min}}, {h(2), Δα, Δf} and {h(1), h(2), Δf}. One notes that each of the first three combinations contains one multifractal parameter from {h(−3), α _{min}, Δα} and the fourth combination consists of the three parameters that produce the three smallest I _{0} values. As in the procedure proposed in the previous subsection, we place the 1125 leaf samples into the four new threedimensional spaces and also use the SVMKM to distinguish them. Under the K − fold cross validation, the discriminant accuracies with increasing K are shown in Figure 13. Obviously, the highest accuracy still comes from the combination {h(−3), α _{min}, Δα} for each K and the lowest accuracy comes from the combination {h(1), h(2), Δf}. This again suggests that the index I _{0} successfully indicates the optimal multifractal parameter combination.
Conclusions
In this paper we have adopted the 2D MFDFA method proposed in [32] to extract important texture features from leaf images. This allow us to calculate the generalized Hurst exponents, h(q), and several other multifractal parameters including α _{max}, α _{min}, ∆α, ∆f, D _{1} and D _{2}. By defining an index, I _{0}, which examines the variation of the interspecies variances and the intraspecies variances, we are able to find an optimal combination of the multifractal parameters that best characterizes the key features of plant species allowing high accuracy in plant species identification. For the Swedish leaf data set which contains 15 species and 75 × 15 = 1125 samples in total [31], the combination of {h(−3), α _{min}, Δα} turns out to be optimal compared to other combinations of parameters. We have obtained 98.4% of averaged discriminant accuracy for every two species by SVMKM with the 10 − fold cross validation, while the accuracy reaches 93.96% for the overall 15 species. Software based on our work can be designed and coded, for that purpose, we provided the corresponding flow chart in the Figure 14.
We should point out that most of the existing work on texture image recognition focuses mainly on the standard multifractal analysis. Our work has shown that the MFDFA is of particular practice for plant leaf identification as the MFDFA multifractal parameters can be combined to distinguish similar but different leaf textures.
References
Agarwal G, Belhumeur P, Feiner S, Jacobs D, Kress WJ, Ramamoorthi R, et al. First steps toward an electronic field guide for plants. Taxon. 2006;55(3):597–610.
Judd W, Campbell CS, Kellog EA, Stevens PF. Plant Systematics: A Phylogenetic Approach. Massachusetts: Sinauer Associates; 1999.
Kumannn MH, Hemsley AR. The Evolution of Plant Architecture. Kew, London: Royal Botanic Gardents; 1999.
Ling H, Jacobs DW. Shape classification using the innerdistance. IEEE trans Pattern Analysis and Machine Intelligence. 2007;29(2):286–99.
Neto JC, Meyer GE, Jones DD, Samal AK. Plant species identification using elliptic Fourier leaf shape analysis. Comput Electron Agr. 2006;50:121–34.
Wang B, Gao Y. Hierarchical string cuts: a translation, rotation, scale, and mirror invariant descriptor for fast shape retrieval. IEEE Trans Image Process. 2014;23(9):4101–11.
Bruno OM, de Oliverira Plotze R, Falvo M, de Castro M. Fractal dimension applied to plant identification. Inform Sciences. 2008;178:2722–33.
Backes AR, Bruno OM. Plant Leaf Identification Using MultiScale Fractal Dimension, in International Conference On Image Analysis and Processing. Berlin/ Heidelberg: Springer; 2009. p. 143150.
Backes AR, Goncalves WN, Martinez AS, Bruno OM. Texture analysis and classification using deterministic tourist walk. Pattern Recogn. 2010;43:685–94.
Cope JS, Remagnino P, Barman S, Wilkin P. Plant Texture Classification Using Gabor coOccurrences, in International Symposium on Visual Computing. Las Vegas Nevada, USA: SpringerVerlag; 2010.
Liu J, Zhang S, Deng S. A method of plant classification based on wavelet transforms and support vector machines. Emerging Intelligent Computing Technology and Applications. 2009;5754:253–60.
Cope JS, Corney D, Clark JY, Remagnino P, Wilkin P. Plant species identification using digital morphometrics: a review. Expert Syst Appl. 2012;39:7562–73.
Du JX, Wang XF, Zhang GJ. Leaf shape based plant species recognition. Appl Math Comput. 2007;185:883–93.
Zhang SW, Lei YK, Dong TB, Zhang XP. Label propagation based supervised locality projection analysis for plant leaf classification. Pattern Recogn. 2013;46:1891–7.
Kong Y, Wang SH, Ma CW, Li BM, Yao YC. Effect of image processing of a leaf photograph on the calculated fractal dimension of leaf veins. Computer and computing technologies in agriculture. 2008;259:755–60.
Rossatto DR, Casanova D, Kolb RM. Fractal analysis of leaftexture properties as a tool for taxonomic and identification purposes: a case study with species from Neotropical Melastomataceae (Miconieae tribe). Plant Syst Evol. 2011;29:103–16.
Chaudhuri BB, Sarkar N. Texture segmentation using fractal dimension. IEEE Trans Pattern Anal Mach Intell. 1995;17:72–7.
Backes AR, Bruno OM. Plant leaf identification using color and multiscale fractal. Dimension Image and Signal Processing. 2010;6134:463–70.
Backes AR, Bruno OM. Plant leaf identification using multiscale fractal dimension. Image Analysis and Processing. 2009;5716:143–50.
Lopes R, Dubois P, Bhouri I, Bedoui MH, Maouche S, Betrouni N. Local fractal and multifractal features for volumic texture characterization. Pattern Recogn. 2011;44:1690–7.
Wang F, Liao GP, Wang XQ, Li JH, Li JW, Shi W. Feature description for nutrient deficiency rape leaves based on multifractal theory. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE). 2013;29(24):181–9 (in Chinese with English abstract).
Xia Y, Feng DG, Zhao RC. Morpholgybased multifractal estimation for texture segmentation. IEEE trans image processing. 2006;15:614–22.
Kantelhardt JW, Zschiegner SA, Bunde EK, Havlin S, Bunde A, Stanley HE. Multifractal detrended fluctuation analysis of nonstationary time series. Phys A. 2002;316:87–114.
Du G, Ning X. Multifractal properties of Chinese stock market in Shanghai. Physica A. 2008;387:261–9.
Hedayatifar L, Vahabi M, Jafari GR. Coupling detrended fluctuation analysis for analyzing coupled nonstationary signals. Phys Rev E. 2011;84:021138.
Wang F, Liao GP, Li JH, Zhou TJ, Li XC. Multifractal detrended fluctuation analysis for clustering structures of electricity price periods. Physica A. 2013;392(22):5723–34.
Wang F, Liao GP, Zhou XY, Shi W. Multifractal detrended crosscorrelation analysis for power markets. Nonlinear Dynamics. 2013;72(1–2):353–63.
Wang F, Li JW, Shi W, Liao GP. Leaf image segmentation method based on multifractal detrended fluctuation analysis. J Appl Phys. 2013;114(21):214905.
Wang F, Li ZS, Liao GP. Multifractal detrended fluctuation analysis for image texture feature representation. Int J Pattern Recognit Artif Intell. 2014;28(3):1455005.
Wang F, Li ZS, Li JW. Local multifractal detrended fluctuation Analysis for Nonstationary Image’s Texture Segmentation. Appl Surf Sci. 2014;233:116–25.
Söderkvist OJO. Computer Vision Classification of Leaf from Swedish Trees, Master Thesis, Linköping University, 2001
Gu GF, Zhou WX. Detrended fluctuation analysis for fractals and multifractals in higher dimensions. Phys Rev E. 2006;74:061104.
Lei YK, Zou JW, Dong T, You ZH, Yuan Y. Orthogonal locally discriminant spline embedding for plant leaf recognition. Comput Vis Image Underst. 2014;119:116–26.
Wu JX, Rehg JM. A visual descriptor for scene categorization. IEEE trans Pattern Analysis and Machine Intelligence. 2010;33(8):1489–501.
Zhang SW, Zhang CL, Cheng L. Plant leaf image classification based on supervised orthogonal locality preserving projections. Transactions of the Chinese Society of Agricultural Engineering. 2013;29(5):125–31.
Canu S, Grandvalet Y, Guigue V, Rakotomamonjy A. SVM and kernel methods matlab toolbox, perception systemes et information. Rouen, France: INSA de Rouen; 2005.
Székely GJ, Rizzo ML. Hierarchical clustering via joint between within distances: extending ward’s winimum variance method. J Classif. 2005;22:151–83.
Acknowledgments
The authors wish to thank the two reviewers and the editorinChief Prof. Brian G. Forde for their comments and suggestions, which led to a great improvement to the presentation of this work.
FW was supported by the Young Scholar Project Funds of the Education Department of Hunan Province (14B087). JL and GL were supported by the Major project funds of Science and technology program of Hunan Province (2013FJ10064). DL was partially supported by a grant from the Forestry Department of Hunan Province, China.
Author information
Authors and Affiliations
Corresponding author
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
FW initiated the study, carried out modeling, programming and drafted the manuscript. DL participated in the design of the study and performed the statistical analysis. JL participated in programming. GL participated in revising the final version. All authors read and approved the final manuscript.
Rights and permissions
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Wang, F., Liao, Dw., Li, Jw. et al. Twodimensional multifractal detrended fluctuation analysis for plant identification. Plant Methods 11, 12 (2015). https://doi.org/10.1186/s1300701500497
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s1300701500497