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.

Alternative Arnoldi process for ill-conditioned tensor equations with application to image restoration

Discrete cosine transform LSQR and GMRES methods for multidimensional ill-posed problems

Tensor Arnoldi-Tikhonov and GMRES-type methods for ill-posed problems with a t-product structure

Tensor Krylov subspace methods via the T-product for large Sylvester tensor equations

A preconditioned Krylov subspace method for linear inverse problems with general-form Tikhonov regularization

The Order-p Tensor Linear Complementarity Problem for Images Deblurring

Krylov subspace methods for large multidimensional eigenvalue computation

Regularized Kaczmarz Algorithms for Tensor Recovery

A New Alternating Minimization Algorithm for Total Variation Image Reconstruction

Tensor Conjugate Gradient Methods with Automatically Determination of Regularization Parameters for Ill-Posed Problems with t-Product

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

Restarting for the Tensor Infinite Arnoldi method

Handling The Non-Smooth Challenge in Tensor SVD: A Multi-Objective Tensor Recovery Framework

Alternating direction algorithms for total variation deconvolution in image reconstruction

A new efficient alternating method for image restoration and texture extraction

Matrix Forms of Iterative Algorithms to Solve Large-Scale Discrete Ill-Posed Problems with an Application to Image Restoration.

Medical Image Reconstruction Using Kernel Based Methods

Simultaneous Rectification and Alignment via Robust Recovery of Low-rank Tensors.

A fast iterative regularization method for ill-posed problems

Third-Order Tensors as Operators on Matrices: A Theoretical and Computational Framework with Applications in Imaging

Tensor Regularized Total Least Squares Methods with Applications to Image and Video Deblurring