MAP Image Recovery with Guarantees using Locally Convex Multi-Scale Energy (LC-MUSE) Model
Jyothi Rikhab Chand, Mathews Jacob
TL;DR
The paper addresses robust MAP reconstruction from undersampled MR data by introducing LC-MuSE, a locally convex multi-scale energy model parameterized by a CNN. By enforcing $m$-strong local convexity via a local monotonicity constraint on the energy gradient, the approach guarantees uniqueness and convergence within a neighborhood while remaining compatible with practical architectures. The energy is learned through a DSM-based objective with a Lipschitz regularization, and the MAP solution is obtained via a Majorization-Minimization update, yielding robust and reliable reconstructions. Empirical results on parallel MR imaging show LC-MuSE outperforms globally convex regularizers and matches or closely approaches state-of-the-art non-convex and end-to-end methods, demonstrating both theoretical guarantees and competitive practical performance across acquisition settings. The framework is broadly applicable to computational imaging problems beyond MR, offering a principled way to combine energy-based priors with deep networks without imposing overly restrictive architectural constraints.
Abstract
We propose a multi-scale deep energy model that is strongly convex in the local neighbourhood around the data manifold to represent its probability density, with application in inverse problems. In particular, we represent the negative log-prior as a multi-scale energy model parameterized by a Convolutional Neural Network (CNN). We restrict the gradient of the CNN to be locally monotone, which constrains the model as a Locally Convex Multi-Scale Energy (LC-MuSE). We use the learned energy model in image-based inverse problems, where the formulation offers several desirable properties: i) uniqueness of the solution, ii) convergence guarantees to a minimum of the inverse problem, and iii) robustness to input perturbations. In the context of parallel Magnetic Resonance (MR) image reconstruction, we show that the proposed method performs better than the state-of-the-art convex regularizers, while the performance is comparable to plug-and-play regularizers and end-to-end trained methods.
