Estimation of multivariate self-similarity exponents

Multivariate self-similarity exponents across multivariate time series, modeled by an operator fractional Brownian motion (ofBm), are classically estimated by a linear regression of wavelet eigenvalues across scales. A bias corrected estimation is performed by the computation of wavelet eigenvalues from the same number of wavelet coefficients. A clustering strategy is built from successive pairwise tests between close estimates.

OFBM Toolbox (GitHub)
Clustering strategy: C.-G. Lucas, P. Abry, H. Wendt and G. Didier, 2022 (pdf)
Estimation: C.-G. Lucas, G. Didier, H. Wendt and P. Abry, 2024 (pdf)

Automated image reconstruction and contour detection

Joint piecewise smooth image reconstruction and contour detection from an observed noisy image is performed by the minimization of the Discrete Mumford-Shah (D-MS) functional, which enforces constraints related to smoothing over image and sparsity over contours to obtain a piecewise smooth reconstructed image `u` and sparse estimated contours `e` from observed noisy images `z`: $$\underset{u,e}{\rm minimize} \quad \Vert u - z \Vert_2 + \beta \Vert (1-e) \odot Du \Vert_2 + \lambda \Vert e \Vert_1,$$ where `\odot` denotes the component-wise product, `D` is a discrete difference operator, and `\beta > 0` and `\lambda > 0` are regularization parameters. A Stein-like strategy providing optimal hyperparameters `\beta` and `\lambda` is designed, based on the minimization of an unbiased quadratic risk estimate. Efficient and automated minimization of the risk estimate relies on an unbiased estimate of the risk's gradient with respect to hyperparameters. The resulting parameter-free procedure, not requiring prior image processing expertise, is thus amenable to real-world applications.

Here is a comparison with the state-of-the-art SUGAR T-ROF (see Cai and Steidl, 2013 and Deledalle et al., 2014) run using gsugar on a BSD69 dataset image:

SUGAR D-MS Toolbox (GitHub)
C.-G. Lucas, B. Pascal, N. Pustelnik and P. Abry, 2021 (pdf)