English

Hessian Schatten-Norm Regularization for Linear Inverse Problems

Optimization and Control 2014-02-20 v3 Computer Vision and Pattern Recognition Numerical Analysis

Abstract

We introduce a novel family of invariant, convex, and non-quadratic functionals that we employ to derive regularized solutions of ill-posed linear inverse imaging problems. The proposed regularizers involve the Schatten norms of the Hessian matrix, computed at every pixel of the image. They can be viewed as second-order extensions of the popular total-variation (TV) semi-norm since they satisfy the same invariance properties. Meanwhile, by taking advantage of second-order derivatives, they avoid the staircase effect, a common artifact of TV-based reconstructions, and perform well for a wide range of applications. To solve the corresponding optimization problems, we propose an algorithm that is based on a primal-dual formulation. A fundamental ingredient of this algorithm is the projection of matrices onto Schatten norm balls of arbitrary radius. This operation is performed efficiently based on a direct link we provide between vector projections onto q\ell_q norm balls and matrix projections onto Schatten norm balls. Finally, we demonstrate the effectiveness of the proposed methods through experimental results on several inverse imaging problems with real and simulated data.

Keywords

Cite

@article{arxiv.1209.3318,
  title  = {Hessian Schatten-Norm Regularization for Linear Inverse Problems},
  author = {Stamatios Lefkimmiatis and John Paul Ward and Michael Unser},
  journal= {arXiv preprint arXiv:1209.3318},
  year   = {2014}
}

Comments

15 pages double-column format. This manuscript will appear in IEEE Transactions on Image Processing

R2 v1 2026-06-21T22:05:21.872Z