“…In the present contribution, the estimation of the Hölder exponent is only performed using wavelet leaders L (γ) f (j, k), as preliminary contributions [27][28][29]32] report that wavelet leader based estimation outperforms those based on other multiresolution quantities. To indicate that the estimation of h and the segmentation procedures proposed in Section 3 below could be applied using any other multiresolution quantity, e.g., the modulus of the 2D-DWT coefficients, a generic notation X(j, k) is used instead of the specific L In the discrete setting, Hölder exponents are estimated for the locations x = 2k associated with the finest scale j = 1, and we make use of the notation · when we are dealing with matrices, rather than with matrix elements, e.g.,…”