Abstract:We consider the solution of the $\ell_1$ regularized image deblurring problem using isotropic and anisotropic regularization implemented with the split Bregman algorithm. For large scale problems, we replace the system matrix $A$ using a Kronecker product approximation obtained via an approximate truncated singular value decomposition for the reordered matrix $\mathcal{R}(A)$. To obtain the approximate decomposition for $\mathcal{R}(A)$ we propose the enlarged Golub Kahan Bidiagonalization algorithm that proceeds by enlarging the Krylov subspace beyond either a given rank for the desired approximation, or uses an automatic stopping test that provides a suitable rank for the approximation. The resultant expansion is contrasted with the use of the truncated and the randomized singular value decompositions with the same number of terms. To further extend the scale of problem that can be considered we implement the determination of the approximation using single precision, while performing all steps for the regularization in standard double precision. The reported numerical tests demonstrate the effectiveness of applying the approximate single precision Kronecker product expansion for $A$, combined with either isotropic or anisotropic regularization implemented using the split Bregman algorithm, for the solution of image deblurring problems. As the size of the problem increases, our results demonstrate that the major costs are associated with determining the Kronecker product approximation, rather than with the cost of the regularization algorithm. Moreover, the enlarged Golub Kahan Bidiagonalization algorithm competes favorably with the randomized singular value decomposition for estimating the approximate singular value decomposition.

Singular Value Decomposition Approximation via Kronecker Summations for Imaging Applications

Representing Images By Multiple Kronecker Product Sum

Limited Memory Block Krylov Subspace Optimization for Computing Dominant Singular Value Decompositions.

Split Bregman Isotropic and Anisotropic Image Deblurring with Kronecker Product Sum Approximations using Single Precision Enlarged-GKB or RSVD Algorithms to provide low rank truncated SVDs

A Hybrid Algorithm for Computing a Partial Singular Value Decomposition Satisfying a Given Threshold

High Order Numerical Methods To Approximate The Singular Value Decomposition

Accelerated singular value thresholding for matrix completion.

A Fast Spectral Relaxation Approach to Matrix Completion Via Kronecker Products.

Approximation of Images via Generalized Higher Order Singular Value Decomposition over Finite-dimensional Commutative Semisimple Algebra

The Singular Value Decomposition, Applications and Beyond.

Singular value decomposition: applications and computations

Better than best low-rank approximation with the singular value decomposition

A Power Method for Computing Singular Value Decomposition

Large-scale Nyström kernel matrix approximation using randomized SVD.

Hybrid Kronecker Product Decomposition and Approximation

Separable Compressive Imaging Method Based on Singular Value Decomposition

Fast Singular Value Shrinkage with Chebyshev Polynomial Approximation Based on Signal Sparsity

The Application of ACA and SVD Algorithms in the Matrix Low-Rank Approximation

Incremental algorithms for truncated higher-order singular value decompositions

Best subspace tensor approximations

Towards a more robust algorithm for computing the restricted singular value decomposition