Variational Multiscale Nonparametric Regression: Algorithms and Implementation

Abstract

Many modern statistically efficient methods come with tremendous computational challenges, often leading to large-scale optimisation problems. In this work, we examine such computational issues for recently developed estimation methods in nonparametric regression with a specific view on image denoising. We consider in particular certain variational multiscale estimators which are statistically optimal in minimax sense, yet computationally intensive. Such an estimator is computed as the minimiser of a smoothness functional (e.g., TV norm) over the class of all estimators such that none of its coefficients with respect to a given multiscale dictionary is statistically significant. The so obtained multiscale Nemirowski-Dantzig estimator (MIND) can incorporate any convex smoothness functional and combine it with a proper dictionary including wavelets, curvelets and shearlets. The computation of MIND in general requires to solve a high-dimensional constrained convex optimisation problem with a specific structure of the constraints induced by the statistical multiscale testing criterion. To solve this explicitly, we discuss three different algorithmic approaches: the Chambolle-Pock, ADMM and semismooth Newton algorithms. Algorithmic details and an explicit implementation is presented and the solutions are then compared numerically in a simulation study and on various test images. We thereby recommend the Chambolle-Pock algorithm in most cases for its fast convergence. We stress that our analysis can also be transferred to signal recovery and other denoising problems to recover more general objects whenever it is possible to borrow statistical strength from data patches of similar object structure.

arXiv commentCodes are available at https://github.com/housenli/MIND

References in corpus (11)

• Beyond a Gaussian Denoiser: Residual Learning of Deep CNN for Image Denoising
• The Dantzig selector: Statistical estimation when $p$ is much larger than $n$
• Variance estimation in nonparametric regression via the difference sequence method
• A data-driven block thresholding approach to wavelet estimation
• Statistical Multiresolution Dantzig Estimation in Imaging: Fundamental Concepts and Algorithmic Framework
• General empirical Bayes wavelet methods and exactly adaptive minimax estimation
• Total variation regularization of multi-material topology optimization
• Extreme Value Analysis of Empirical Frame Coefficients and Implications for Denoising by Soft-Thresholding
• Variational Multiscale Nonparametric Regression: Smooth Functions
• A Globally Linearly Convergent Method for Pointwise Quadratically Supportable Convex-Concave Saddle Point Problems
• Parallel Statistical Multi-resolution Estimation

Cited by in corpus (1)

• A Variational View on Statistical Multiscale Estimation