Multiscale multifractal detrended-fluctuation analysis of two-dimensional surfaces

Two-dimensional (2D) multifractal detrended fluctuation analysis (MF-DFA) has been used to study monofractality and multifractality on 2D surfaces, but when it is used to calculate the generalized Hurst exponent in a fixed time scale, the presence of crossovers can bias the outcome. To solve this problem, multiscale multifractal analysis (MMA) was recent employed in a one-dimensional case. MMA produces a Hurst surface h(q,s) that provides a spectrum of local scaling exponents at different scale ranges such that the positions of the crossovers can be located. We apply this MMA method to a 2D surface and identify factors that influence the results. We generate several synthesized surfaces and find that crossovers are consistently present, which means that their fractal properties differ at different scales. We apply MMA to the surfaces, and the results allow us to observe these differences and accurately estimate the generalized Hurst exponents. We then study eight natural texture images and two real-world images and find (i) that the moving window length (WL) and the slide length (SL) are the key parameters in the MMA method, that the WL more strongly influences the Hurst surface than the SL, and that the combination of WL = 4 and SL = 4 is optimal for a 2D image; (ii) that the robustness of h(2, s) to four common noises is high at large scales but variable at small scales; and (iii) that the long-term correlations in the images weaken as the intensity of Gaussian noise and salt and pepper noise is increased. Our findings greatly improve the performance of the MMA method on 2D surfaces.