Maximum likelihood estimation of regularization parameters in high-dimensional inverse problems: An empirical bayesian approach part i: Methodology and experiments

Ana Fernandez Vidal, Valentin De Bortoli, Marcelo Pereyra, Alain Durmus

Research output: Contribution to journalArticlepeer-review

2 Citations (Scopus)

Abstract

Many imaging problems require solving an inverse problem that is ill-conditioned or ill-posed. Imaging methods typically address this difficulty by regularizing the estimation problem to make it well-posed. This often requires setting the value of the so-called regularization parameters that control the amount of regularization enforced. These parameters are notoriously difficult to set a priori and can have a dramatic impact on the recovered estimates. In this work, we propose a general empirical Bayesian method for setting regularization parameters in imaging problems that are convex w.r.t. the unknown image. Our method calibrates regularization parameters directly from the observed data by maximum marginal likelihood estimation and can simultaneously estimate multiple regular-ization parameters. Furthermore, the proposed algorithm uses the same basic operators as proximal optimization algorithms, namely gradient and proximal operators, and it is therefore straightfor-ward to apply to problems that are currently solved by using proximal optimization techniques. Our methodology is demonstrated with a range of experiments and comparisons with alternative approaches from the literature. The considered experiments include image denoising, nonblind image deconvolution, and hyperspectral unmixing, using synthesis and analysis priors involving the ℓ1, total-variation, total-variation and ℓ1, and total-generalized-variation pseudonorms. A detailed theoretical analysis of the proposed method is presented in our companion paper.

Original languageEnglish
Pages (from-to)1945-1989
Number of pages45
JournalSIAM Journal on Imaging Sciences
Volume13
Issue number4
DOIs
Publication statusPublished - 18 Nov 2020

Keywords

  • Empirical Bayes
  • Image processing
  • Inverse problems
  • Markov chain Monte Carlo methods
  • Proximal algorithms
  • Statistical inference
  • Stochastic optimization

ASJC Scopus subject areas

  • Mathematics(all)
  • Applied Mathematics

Fingerprint

Dive into the research topics of 'Maximum likelihood estimation of regularization parameters in high-dimensional inverse problems: An empirical bayesian approach part i: Methodology and experiments'. Together they form a unique fingerprint.

Cite this