Source author record

Herwig Wendt

Herwig Wendt appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

11works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

11 published item(s)

preprint2024arXiv

On the empirical spectral distribution of large wavelet random matrices based on mixed-Gaussian fractional measurements in moderately high dimensions

In this paper, we characterize the convergence of the (rescaled logarithmic) empirical spectral distribution of wavelet random matrices. We assume a moderately high-dimensional framework where the sample size $n$, the dimension $p(n)$ and, for a fixed integer $j$, the scale $a(n)2^j$ go to infinity in such a way that $\lim_{n \rightarrow \infty}p(n)\cdot a(n)/n = \lim_{n \rightarrow \infty} o(\sqrt{a(n)/n})= 0$. We suppose the underlying measurement process is a random scrambling of a sample of size $n$ of a growing number $p(n)$ of fractional processes. Each of the latter processes is a fractional Brownian motion conditionally on a randomly chosen Hurst exponent. We show that the (rescaled logarithmic) empirical spectral distribution of the wavelet random matrices converges weakly, in probability, to the distribution of Hurst exponents.

preprint2022arXiv

Wavelet eigenvalue regression in high dimensions

In this paper, we construct the wavelet eigenvalue regression methodology in high dimensions. We assume that possibly non-Gaussian, finite-variance $p$-variate measurements are made of a low-dimensional $r$-variate ($r \ll p$) fractional stochastic process with non-canonical scaling coordinates and in the presence of additive high-dimensional noise. The measurements are correlated both time-wise and between rows. Building upon the asymptotic and large scale properties of wavelet random matrices in high dimensions, the wavelet eigenvalue regression is shown to be consistent and, under additional assumptions, asymptotically Gaussian in the estimation of the fractal structure of the system. We further construct a consistent estimator of the effective dimension $r$ of the system that significantly increases the robustness of the methodology. The estimation performance over finite samples is studied by means of simulations.

preprint2016arXiv

Combining local regularity estimation and total variation optimization for scale-free texture segmentation

Texture segmentation constitutes a standard image processing task, crucial to many applications. The present contribution focuses on the particular subset of scale-free textures and its originality resides in the combination of three key ingredients: First, texture characterization relies on the concept of local regularity ; Second, estimation of local regularity is based on new multiscale quantities referred to as wavelet leaders ; Third, segmentation from local regularity faces a fundamental bias variance trade-off: In nature, local regularity estimation shows high variability that impairs the detection of changes, while a posteriori smoothing of regularity estimates precludes from locating correctly changes. Instead, the present contribution proposes several variational problem formulations based on total variation and proximal resolutions that effectively circumvent this trade-off. Estimation and segmentation performance for the proposed procedures are quantified and compared on synthetic as well as on real-world textures.

preprint2016arXiv

p-exponent and p-leaders, Part I: Negative pointwise regularity

Multifractal analysis aims to characterize signals, functions, images or fields, via the fluctuations of their local regularity along time or space, hence capturing crucial features of their temporal/spatial dynamics. Multifractal analysis is becoming a standard tool in signal and image processing, and is nowadays widely used in numerous applications of different natures. Its common formulation relies on the measure of local regularity via the Hölder exponent, by nature restricted to positive values, and thus to locally bounded functions or signals. It is here proposed to base the quantification of local regularity on $p$-exponents, a novel local regularity measure potentially taking negative values. First, the theoretical properties of $p$-exponents are studied in detail. Second, wavelet-based multiscale quantities, the $p$-leaders, are constructed and shown to permit accurate practical estimation of $p$-exponents. Exploiting the potential dependence with $p$, it is also shown how the collection of $p$-exponents enriches the classification of locally singular behaviors in functions, signals or images. The present contribution is complemented by a companion article developing the $p$-leader based multifractal formalism associated to $p$-exponents.

preprint2016arXiv

p-exponent and p-leaders, Part II: Multifractal Analysis. Relations to Detrended Fluctuation Analysis

Multifractal analysis studies signals, functions, images or fields via the fluctuations of their local regularity along time or space, which capture crucial features of their temporal/spatial dynamics. It has become a standard signal and image processing tool and is commonly used in numerous applications of different natures. In its common formulation, it relies on the Hölder exponent as a measure of local regularity, which is by nature restricted to positive values and can hence be used for locally bounded functions only. In this contribution, it is proposed to replace the Hölder exponent with a collection of novel exponents for measuring local regularity, the $p$-exponents. One of the major virtues of $p$-exponents is that they can potentially take negative values. The corresponding wavelet-based multiscale quantities, the $p$-leaders, are constructed and shown to permit the definition of a new multifractal formalism, yielding an accurate practical estimation of the multifractal properties of real-world data. Moreover, theoretical and practical connections to and comparisons against another multifractal formalism, referred to as multifractal detrended fluctuation analysis, are achieved. The performance of the proposed $p$-leader multifractal formalism is studied and compared to previous formalisms using synthetic multifractal signals and images, illustrating its theoretical and practical benefits. The present contribution is complemented by a companion article studying in depth the theoretical properties of $p$-exponents and the rich classification of local singularities it permits.

preprint2015arXiv

Bayesian estimation of the multifractality parameter for image texture using a Whittle approximation

Texture characterization is a central element in many image processing applications. Multifractal analysis is a useful signal and image processing tool, yet, the accurate estimation of multifractal parameters for image texture remains a challenge. This is due in the main to the fact that current estimation procedures consist of performing linear regressions across frequency scales of the two-dimensional (2D) dyadic wavelet transform, for which only a few such scales are computable for images. The strongly non-Gaussian nature of multifractal processes, combined with their complicated dependence structure, makes it difficult to develop suitable models for parameter estimation. Here, we propose a Bayesian procedure that addresses the difficulties in the estimation of the multifractality parameter. The originality of the procedure is threefold: The construction of a generic semi-parametric statistical model for the logarithm of wavelet leaders; the formulation of Bayesian estimators that are associated with this model and the set of parameter values admitted by multifractal theory; the exploitation of a suitable Whittle approximation within the Bayesian model which enables the otherwise infeasible evaluation of the posterior distribution associated with the model. Performance is assessed numerically for several 2D multifractal processes, for several image sizes and a large range of process parameters. The procedure yields significant benefits over current benchmark estimators in terms of estimation performance and ability to discriminate between the two most commonly used classes of multifractal process models. The gains in performance are particularly pronounced for small image sizes, notably enabling for the first time the analysis of image patches as small as 64x64 pixels.

preprint2015arXiv

Multifractal analysis based on p-exponents and lacunarity exponents

Many examples of signals and images cannot be modeled by locally bounded functions, so that the standard multifractal analysis, based on the Hölder exponent, is not feasible. We present a multifractal analysis based on another quantity, the p-exponent, which can take arbitrarily large negative values. We investigate some mathematical properties of this exponent, and show how it allows us to model the idea of "lacunarity" of a singularity at a point. We finally adapt the wavelet based multifractal analysis in this setting, and we give applications to a simple mathematical model of multifractal processes: Lacunary wavelet series.

preprint2015arXiv

Multiscale reverse-time-migration-type imaging using the dyadic parabolic decomposition of phase space

We develop a representation of reverse-time migration in terms of Fourier integral operators the canonical relations of which are graphs. Through the dyadic parabolic decomposition of phase space, we obtain the solution of the wave equation with a boundary source and homogeneous initial conditions using wave packets. On this basis, we develop a numerical procedure for the reverse time continuation from the boundary of scattering data and for RTM migration. The algorithms are derived from those we recently developed for the discrete approximate evaluation of the action of Fourier integral operators and inherit from their conceptual and numerical properties.

preprint2013arXiv

Frequency estimation based on Hankel matrices and the alternating direction method of multipliers

We develop a parametric high-resolution method for the estimation of the frequency nodes of linear combinations of complex exponentials with exponential damping. We use Kronecker's theorem to formulate the associated nonlinear least squares problem as an optimization problem in the space of vectors generating Hankel matrices of fixed rank. Approximate solutions to this problem are obtained by using the alternating direction method of multipliers. Finally, we extract the frequency estimates from the con-eigenvectors of the solution Hankel matrix. The resulting algorithm is simple, easy to implement and can be applied to data with equally spaced samples with approximation weights, which for instance allows cases of missing data samples. By means of numerical simulations, we analyze and illustrate the excellent performance of the method, attaining the Cramér-Rao bound.

preprint2013arXiv

Multi-scale discrete approximations of Fourier integral operators associated with canonical transformations and caustics

We develop an algorithm for the computation of general Fourier integral operators associated with canonical graphs. The algorithm is based on dyadic parabolic decomposition using wave packets and enables the discrete approximate evaluation of the action of such operators on data in the presence of caustics. The procedure consists in the construction of a universal operator representation through the introduction of locally singularity-resolving diffeomorphisms, enabling the application of wave packet driven computation, and in the construction of the associated pseudo-differential joint-partition of unity on the canonical graphs. We apply the method to a parametrix of the wave equation in the vicinity of a cusp singularity.

preprint2012arXiv

Irregularities and Scaling in Signal and Image Processing: Multifractal Analysis

B. Mandelbrot gave a new birth to the notions of scale invariance, selfsimilarity and non-integer dimensions, gathering them as the founding corner-stones used to build up fractal geometry. The first purpose of the present contribution is to review and relate together these key notions, explore their interplay and show that they are different facets of a same intuition. Second, it will explain how these notions lead to the derivation of the mathematical tools underlying multifractal analysis. Third, it will reformulate these theoretical tools into a wavelet framework, hence enabling their better theoretical understanding as well as their efficient practical implementation. B. Mandelbrot used his concept of fractal geometry to analyze real-world applications of very different natures. As a tribute to his work, applications of various origins, and where multifractal analysis proved fruitful, are revisited to illustrate the theoretical developments proposed here.