A Direct Algorithm For Optimization Problems With The Huber Penalty

Jingyan Xu, Frederic Noo, Benjamin Tsui

Research output: Contribution to journalArticle


We present a direct (noniterative) algorithm for one dimensional (1-D) quadratic data fitting with neighboring intensity differences penalized by the Huber function. Applications of such an algorithm include 1-D processing of medical signals, such as smoothing of tissue time concentration curves in kinetic data analysis or sinogram preprocessing, and using it as a subproblem solver for 2-D or 3-D image restoration and reconstruction. Dynamic programming (DP) was used to develop the direct algorithm. The problem was reformulated as a sequence of univariate optimization problems, for k = 1,…,N, where N is the number of data points. The solution to the univariate problem at index k is parameterized by the solution at k+1, except at k = N. Solving the univariate optimization problem at k = N yields the solution to each problem in the sequence using backtracking. Computational issues and memory cost are discussed in detail. Two numerical studies, tissue concentration curve smoothing and sinogram preprocessing for image reconstruction, are used to validate the direct algorithm and illustrate its practical applications. In the example of 1-D curve smoothing, the efficiency of the direct algorithm is compared with four iterative methods: the iterative coordinate descent, Nesterov’s accelerated gradient descent algorithm, FISTA, and an off-the-shelf second order method. The first two methods were applied to the primal problem, the others to the dual problem. The comparisons show that the direct algorithm outperforms all other methods by a significant factor, which rapidly grows with the curvature of the Huber function. The second example, sinogram preprocessing, showed that robustness and speed of the direct algorithm are maintained over a wide range of signal variations, and that noise and streaking artifacts could be reduced with almost no increase in computation time. We also outline the application of the proposed 1-D solver for imaging applications.

Original languageEnglish (US)
JournalIEEE Transactions on Medical Imaging
StateAccepted/In press - Oct 4 2017


  • Algorithm design and analysis
  • denoising
  • dynamic programming
  • Heuristic algorithms
  • Huber penalty
  • Imaging
  • Minimization
  • robust estimation
  • Signal processing algorithms
  • smoothing
  • Smoothing methods
  • total variation
  • TV

ASJC Scopus subject areas

  • Software
  • Radiological and Ultrasound Technology
  • Computer Science Applications
  • Electrical and Electronic Engineering

Fingerprint Dive into the research topics of 'A Direct Algorithm For Optimization Problems With The Huber Penalty'. Together they form a unique fingerprint.

  • Cite this