Saved in:
Bibliographic Details
Main Authors: Chen, Zixuan, Nagy, James, Xi, Yuanzhe, Yu, Bo
Format: Preprint
Published: 2019
Subjects:
Online Access:https://arxiv.org/abs/1901.00913
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866910753002356736
author Chen, Zixuan
Nagy, James
Xi, Yuanzhe
Yu, Bo
author_facet Chen, Zixuan
Nagy, James
Xi, Yuanzhe
Yu, Bo
contents In this paper, we propose an efficient numerical scheme for solving some large scale ill-posed linear inverse problems arising from image restoration. In order to accelerate the computation, two different hidden structures are exploited. First, the coefficient matrix is approximated as the sum of a small number of Kronecker products. This procedure not only introduces one more level of parallelism into the computation but also enables the usage of computationally intensive matrix-matrix multiplications in the subsequent optimization procedure. We then derive the corresponding Tikhonov regularized minimization model and extend the fast iterative shrinkage-thresholding algorithm (FISTA) to solve the resulting optimization problem. Since the matrices appearing in the Kronecker product approximation are all structured matrices (Toeplitz, Hankel, etc.), we can further exploit their fast matrix-vector multiplication algorithms at each iteration. The proposed algorithm is thus called structured fast iterative shrinkage-thresholding algorithm (sFISTA). In particular, we show that the approximation error introduced by sFISTA is well under control and sFISTA can reach the same image restoration accuracy level as FISTA. Finally, both the theoretical complexity analysis and some numerical results are provided to demonstrate the efficiency of sFISTA.
format Preprint
id arxiv_https___arxiv_org_abs_1901_00913
institution arXiv
publishDate 2019
record_format arxiv
spellingShingle Structured FISTA for Image Restoration
Chen, Zixuan
Nagy, James
Xi, Yuanzhe
Yu, Bo
Numerical Analysis
In this paper, we propose an efficient numerical scheme for solving some large scale ill-posed linear inverse problems arising from image restoration. In order to accelerate the computation, two different hidden structures are exploited. First, the coefficient matrix is approximated as the sum of a small number of Kronecker products. This procedure not only introduces one more level of parallelism into the computation but also enables the usage of computationally intensive matrix-matrix multiplications in the subsequent optimization procedure. We then derive the corresponding Tikhonov regularized minimization model and extend the fast iterative shrinkage-thresholding algorithm (FISTA) to solve the resulting optimization problem. Since the matrices appearing in the Kronecker product approximation are all structured matrices (Toeplitz, Hankel, etc.), we can further exploit their fast matrix-vector multiplication algorithms at each iteration. The proposed algorithm is thus called structured fast iterative shrinkage-thresholding algorithm (sFISTA). In particular, we show that the approximation error introduced by sFISTA is well under control and sFISTA can reach the same image restoration accuracy level as FISTA. Finally, both the theoretical complexity analysis and some numerical results are provided to demonstrate the efficiency of sFISTA.
title Structured FISTA for Image Restoration
topic Numerical Analysis
url https://arxiv.org/abs/1901.00913