R: Compute the multifractal indicators hmin, c1 and c2
fit-method
R Documentation
Compute the multifractal indicators hmin, c1 and c2
Description
This function estimates the multifractal indicators of an image. It takes as input a class object "leader" and a range of scales and returns a vector containing hmin, c1 and c2 .
Concretely the function fit uses "df.mf" slot and performs a linear regression of hmin, c1 and c2 columns on the choosen scales.
hmin characterizes the uniform smoothness of the image, c1 corresponds to the value of h where L(h) is maximal, and c2 explains the strength
of the multifractality. hmin and c1 are positive and c2 negative. The multifractal spectrum is approximated by :
L(h) = 2 + c2/2 ((h-c1)/c2)^2
hmin is the minimum value of h such that L(h) is greater than 0.
To make comparable analyzes , we substract frac to c1.
If the analysis is limited, we get only hmin. If the estimate of hmin is negative, strictly speaking one should repeat the analysis with an index of fractional integration "frac" greater than -hmin. hmin, t1 and t2 are calculated on the original wavelet coefficients (without the fractional integration).
Usage
## S4 method for signature 'Leader'
fit(object,scales)
Arguments
object
A object of "Leader" class
scales
The scales of range of the estimate. This range must be contained in 1:J
Value
A vector containing 1 (hmin) or 3 (hmin,c1,c2) values.
Author(s)
Francois Semecurbe
References
Herwig Wendt, Stephane Roux, Stephane Jaffard, Patrice Abry. Wavelet leaders and bootstrap
for multifractal analysis of images. Signal Processing, Elsevier, 2009, 6 (89), pp.1100-1114.
Patrice Abry, Herwig Wendt, Stephane Jaffard. When Van Gogh meets Mandelbrot: Multifractal Classification of Painting's Texture. Signal Processing, Elsevier, 2013, 93 (3), pp.554-572.