Learning sparsity-promoting regularizers for linear inverse problems
Giovanni S. Alberti, Ernesto De Vito, Tapio Helin, Matti Lassas, Luca Ratti, Matteo Santacesaria
TL;DR
The work addresses learning sparsity-promoting regularizers for linear inverse problems by formulating a bilevel optimization to select a synthesis operator $B$ that induces sparse representations of the solution. The inner problem enforces data fidelity with an $\ell^1$ sparsity penalty in the $B$-domain, producing $\hat{x}_B=B\hat{u}_B$, while the outer objective minimizes an expected reconstruction error via $L(B)=\mathbb{E}[\|R_B(y)-x\|_X^2]$. The authors establish well-posedness and stability for fixed $B$, provide finite-sample generalization bounds for learning $B$ from data, and demonstrate two concrete parameter classes: compact perturbations of a known operator and learning the mother wavelet, each with explicit covering-number bounds. The results enable data-driven incorporation of prior structure into infinite-dimensional sparsity-promoting regularization, with proven guarantees on existence, uniqueness, stability, and sample complexity. Overall, the framework blends variational regularization, operator learning, and statistical learning theory to promote sparse reconstructions in linear inverse problems.
Abstract
This paper introduces a novel approach to learning sparsity-promoting regularizers for solving linear inverse problems. We develop a bilevel optimization framework to select an optimal synthesis operator, denoted as $B$, which regularizes the inverse problem while promoting sparsity in the solution. The method leverages statistical properties of the underlying data and incorporates prior knowledge through the choice of $B$. We establish the well-posedness of the optimization problem, provide theoretical guarantees for the learning process, and present sample complexity bounds. The approach is demonstrated through examples, including compact perturbations of a known operator and the problem of learning the mother wavelet, showcasing its flexibility in incorporating prior knowledge into the regularization framework. This work extends previous efforts in Tikhonov regularization by addressing non-differentiable norms and proposing a data-driven approach for sparse regularization in infinite dimensions.
