Tikhonov-type iterative regularization methods for ill-posed inverse problems: theoretical aspects and applications

Ill-posed inverse problems arise in many fields of science and engineering. The ill-conditioning and the big dimension make the task of numerically solving this kind of problems very challenging. In this thesis we construct several algorithms for solving ill-posed inverse problems. Starting from the classical Tikhonov regularization method we develop iterative methods that enhance the performances of the originating method. In order to ensure the accuracy of the constructed algorithms we insert a priori knowledge on the exact solution and empower the regularization term. By exploiting the structure of the problem we are also able to achieve fast computation even when the size of the problem becomes very big. We construct algorithms that enforce constraint on the reconstruction, like nonnegativity or flux conservation and exploit enhanced version of the Euclidian norm using a regularization operator and different semi-norms, like the Total Variaton, for the regularization term. For most of the proposed algorithms we provide efficient strategies for the choice of the regularization parameters, which, most of the times, rely on the knowledge of the norm of the noise that corrupts the data. For each method we analyze the theoretical properties in the finite dimensional case or in the more general case of Hilbert spaces. Numerical examples prove the good performances of the algorithms proposed in term of both accuracy and efficiency

Tikhonov-type iterative regularization methods for ill-posed inverse problems: theoretical aspects and applications

Ill-posed inverse problems arise in many fields of science and engineering. The ill-conditioning and the big dimension make the task of numerically solving this kind of problems very challenging. In this thesis we construct several algorithms for solving ill-posed inverse problems. Starting from the classical Tikhonov regularization method we develop iterative methods that enhance the performances of the originating method. In order to ensure the accuracy of the constructed algorithms we insert a priori knowledge on the exact solution and empower the regularization term. By exploiting the structure of the problem we are also able to achieve fast computation even when the size of the problem becomes very big. We construct algorithms that enforce constraint on the reconstruction, like nonnegativity or flux conservation and exploit enhanced version of the Euclidian norm using a regularization operator and different semi-norms, like the Total Variaton, for the regularization term. For most of the proposed algorithms we provide efficient strategies for the choice of the regularization parameters, which, most of the times, rely on the knowledge of the norm of the noise that corrupts the data. For each method we analyze the theoretical properties in the finite dimensional case or in the more general case of Hilbert spaces. Numerical examples prove the good performances of the algorithms proposed in term of both accuracy and efficiency

Diffusion equations and inverse problems regularization.

The present thesis can be split into two dfferent parts: The first part mainly deals with the porous and fast diffusion equations. Chapter 2 presents these equations in the Euclidean setting highlighting the technical issues that arise when trying to extend results in a Riemannian setting. Chapter 3 is devoted to the construction of exhaustion and cut-o_ functions with controlled gradient and Laplacian, on manifolds with Ricci curvature bounded from below by a (possibly unbounded) nonpositive function of the distance from a fixed reference point, and without any assumptions on the topology or the injectivity radius. The cut-offs are then applied to the study of the fast and porous media diffusion, of Lq-properties of the gradient and of the selfadjointness of Schrödinger-type operators. The second part is concerned with inverse problems regularization applied to image deblurring. In Chapter 5 new variants of the Tikhonov filter method, called fractional and weighted Tikhonov, are presented alongside their saturation properties and converse results on their convergence rates. New iterated fractional Tikhonov regularization methods are then introduced. In Chapter 6 the modified linearized Bregman algorithm is investigated. It is showed that the standard approach based on the block circulant circulant block preconditioner may provide low quality restored images and different preconditioning strategies are then proposed, which improve the quality of the restoration

Diffusion equations and inverse problems regularization.

The present thesis can be split into two dfferent parts: The first part mainly deals with the porous and fast diffusion equations. Chapter 2 presents these equations in the Euclidean setting highlighting the technical issues that arise when trying to extend results in a Riemannian setting. Chapter 3 is devoted to the construction of exhaustion and cut-o_ functions with controlled gradient and Laplacian, on manifolds with Ricci curvature bounded from below by a (possibly unbounded) nonpositive function of the distance from a fixed reference point, and without any assumptions on the topology or the injectivity radius. The cut-offs are then applied to the study of the fast and porous media diffusion, of Lq-properties of the gradient and of the selfadjointness of Schrödinger-type operators. The second part is concerned with inverse problems regularization applied to image deblurring. In Chapter 5 new variants of the Tikhonov filter method, called fractional and weighted Tikhonov, are presented alongside their saturation properties and converse results on their convergence rates. New iterated fractional Tikhonov regularization methods are then introduced. In Chapter 6 the modified linearized Bregman algorithm is investigated. It is showed that the standard approach based on the block circulant circulant block preconditioner may provide low quality restored images and different preconditioning strategies are then proposed, which improve the quality of the restoration

An Arnoldi-based preconditioner for iterated Tikhonov regularization

Many problems in science and engineering give rise to linear systems of equations that are commonly referred to as large-scale linear discrete ill-posed problems. These problems arise, for instance, from the discretization of Fredholm integral equations of the first kind. The matrices that define these problems are typically severely ill-conditioned and may be rank-deficient. Because of this, the solution of linear discrete ill-posed problems may not exist or be very sensitive to perturbations caused by errors in the available data. These difficulties can be reduced by applying Tikhonov regularization. We describe a novel "approximate Tikhonov regularization method" based on constructing a low-rank approximation of the matrix in the linear discrete ill-posed problem by carrying out a few steps of the Arnoldi process. The iterative method so defined is transpose-free. Our work is inspired by a scheme by Donatelli and Hanke, whose approximate Tikhonov regularization method seeks to approximate a severely ill-conditioned block-Toeplitz matrix with Toeplitz-blocks by a block-circulant matrix with circulant-blocks. Computed examples illustrate the performance of our proposed iterative regularization method

A conjugate-gradient-type rational Krylov subspace method for ill-posed problems

Conjugated gradients on the normal equation (CGNE) is a popular method to regularise linear inverse problems. The idea of the method can be summarized as minimising the residuum over a suitable Krylov subspace. It is shown that using the same idea for the shift-and-invert rational Krylov subspace yields an order-optimal regularisation scheme

A comparison of parameter choice rules for ℓp - ℓq minimization

Images that have been contaminated by various kinds of blur and noise can be restored by the minimization of an ℓp-ℓq functional. The quality of the reconstruction depends on the choice of a regularization parameter. Several approaches to determine this parameter have been described in the literature. This work presents a numerical comparison of known approaches as well as of a new one