Overcomplete representation in a hierarchical Bayesian framework

TL;DR

This paper extends Bayesian hierarchical models with an iterative algorithm to efficiently find sparse solutions in overcomplete systems, benefiting large-scale inverse problems and machine learning applications.

Contribution

It introduces a hybrid IAS algorithm for overcomplete sparse coding, improving efficiency and effectiveness in identifying the most sparse solutions.

Findings

Effective in large-scale inverse problems

Suitable for dictionary learning in machine learning

Outperforms traditional methods in sparsity recovery

Abstract

A common task in inverse problems and imaging is finding a solution that is sparse, in the sense that most of its components vanish. In the framework of compressed sensing, general results guaranteeing exact recovery have been proven. In practice, sparse solutions are often computed combining $\ell_1$-penalized least squares optimization with an appropriate numerical scheme to accomplish the task. A computationally efficient alternative for finding sparse solutions to linear inverse problems is provided by Bayesian hierarchical models, in which the sparsity is encoded by defining a conditionally Gaussian prior model with the prior parameter obeying a generalized gamma distribution. An iterative alternating sequential (IAS) algorithm has been demonstrated to lead to a computationally efficient scheme, and combined with Krylov subspace iterations with an early termination condition, the approach is particularly well suited for large scale problems. Here the Bayesian approach to sparsity is extended to problems whose solution allows a sparse coding in an overcomplete system such as composite frames. It is shown that among the multiple possible representations of the unknown, the IAS algorithm, and in particular, a hybrid version of it, is effectively identifying the most sparse solution. Computed examples show that the method is particularly well suited not only for traditional imaging applications but also for dictionary learning problems in the framework of machine learning.