Implicit regularization of the incomplete oblique projections method
The aim of this paper is to improve the performance of the incomplete oblique projections method (IOP), previously introduced by the authors for solving inconsistent linear systems, when applied to image reconstruction problems. That method employs incomplete oblique projections onto the set of solu...
Guardado en:
Publicado: |
2009
|
---|---|
Materias: | |
Acceso en línea: | https://bibliotecadigital.exactas.uba.ar/collection/paper/document/paper_09696016_v16_n4_p525_Scolnik http://hdl.handle.net/20.500.12110/paper_09696016_v16_n4_p525_Scolnik |
Aporte de: |
id |
paper:paper_09696016_v16_n4_p525_Scolnik |
---|---|
record_format |
dspace |
spelling |
paper:paper_09696016_v16_n4_p525_Scolnik2023-06-08T15:59:01Z Implicit regularization of the incomplete oblique projections method Computerized tomographies Image reconstruction Incomplete projections Least squares problems Minimum norm solution Regularizing The aim of this paper is to improve the performance of the incomplete oblique projections method (IOP), previously introduced by the authors for solving inconsistent linear systems, when applied to image reconstruction problems. That method employs incomplete oblique projections onto the set of solutions of the augmented system Ax−r=b, and converges to a weighted least squares solution of the system Ax=b. Many tomographic image reconstruction problems are such that the limitation of the range of rays makes the model underdetermined, the discretized linear system is rank-deficient, the nullspace is non-trivial, and the minimal norm least squares solution may be far away from the true image. In a previous paper, we have added a quadratic term reflecting neighboring pixel information to the standard least squares model for improving the quality of the reconstructed images. In this paper we replace the quadratic function by a more general regularizing function avoiding the modification of the original system. The key idea is to perform a joint optimization of the norm of the residual and of the regularizing function in each iteration. The theoretical properties of this new algorithm are analyzed, and numerical experiments are presented comparing its performance with other well-known methods. They show that the new approach improves the quality of the reconstructed images. © 2009 International Federation of Operational Research Societies. 2009 https://bibliotecadigital.exactas.uba.ar/collection/paper/document/paper_09696016_v16_n4_p525_Scolnik http://hdl.handle.net/20.500.12110/paper_09696016_v16_n4_p525_Scolnik |
institution |
Universidad de Buenos Aires |
institution_str |
I-28 |
repository_str |
R-134 |
collection |
Biblioteca Digital - Facultad de Ciencias Exactas y Naturales (UBA) |
topic |
Computerized tomographies Image reconstruction Incomplete projections Least squares problems Minimum norm solution Regularizing |
spellingShingle |
Computerized tomographies Image reconstruction Incomplete projections Least squares problems Minimum norm solution Regularizing Implicit regularization of the incomplete oblique projections method |
topic_facet |
Computerized tomographies Image reconstruction Incomplete projections Least squares problems Minimum norm solution Regularizing |
description |
The aim of this paper is to improve the performance of the incomplete oblique projections method (IOP), previously introduced by the authors for solving inconsistent linear systems, when applied to image reconstruction problems. That method employs incomplete oblique projections onto the set of solutions of the augmented system Ax−r=b, and converges to a weighted least squares solution of the system Ax=b. Many tomographic image reconstruction problems are such that the limitation of the range of rays makes the model underdetermined, the discretized linear system is rank-deficient, the nullspace is non-trivial, and the minimal norm least squares solution may be far away from the true image. In a previous paper, we have added a quadratic term reflecting neighboring pixel information to the standard least squares model for improving the quality of the reconstructed images. In this paper we replace the quadratic function by a more general regularizing function avoiding the modification of the original system. The key idea is to perform a joint optimization of the norm of the residual and of the regularizing function in each iteration. The theoretical properties of this new algorithm are analyzed, and numerical experiments are presented comparing its performance with other well-known methods. They show that the new approach improves the quality of the reconstructed images. © 2009 International Federation of Operational Research Societies. |
title |
Implicit regularization of the incomplete oblique projections method |
title_short |
Implicit regularization of the incomplete oblique projections method |
title_full |
Implicit regularization of the incomplete oblique projections method |
title_fullStr |
Implicit regularization of the incomplete oblique projections method |
title_full_unstemmed |
Implicit regularization of the incomplete oblique projections method |
title_sort |
implicit regularization of the incomplete oblique projections method |
publishDate |
2009 |
url |
https://bibliotecadigital.exactas.uba.ar/collection/paper/document/paper_09696016_v16_n4_p525_Scolnik http://hdl.handle.net/20.500.12110/paper_09696016_v16_n4_p525_Scolnik |
_version_ |
1768544190155718656 |