The mathematics of computerized tomography
The mathematics of computerized tomography
Parallel Optimization: Theory, Algorithms and Applications
Parallel Optimization: Theory, Algorithms and Applications
On image reconstruction algorithms for binary electromagnetic geotomography
Theoretical Computer Science
Adaptive wavelet-Galerkin methods for limited angle tomography
Image and Vision Computing
Motion and positional error correction for cone beam 3D-reconstruction with mobile C-arms
MICCAI'07 Proceedings of the 10th international conference on Medical image computing and computer-assisted intervention - Volume Part I
A hybrid Kaczmarz-Conjugate Gradient algorithm for image reconstruction
Mathematics and Computers in Simulation
Image reconstruction by an alternating minimisation
Neurocomputing
Optimization for limited angle tomography in medical image processing
Pattern Recognition
| Hi-index | 0.00 |
The algebraic reconstruction technique (ART), based on the well known algorithm proposed by Kaczmarz in 1937, is one of the most important class of solution methods for image reconstruction problems. But unfortunately, almost all the methods from the ART class give satisfactory results only in the case of consistent problems. In the inconsistent case (and unfortunately this is what happens in real applications, because of measurement errors) they give only more or less "controllable" versions of the exact solutions. This is exactly the case that we analyze in the present paper. We start with a theoretical analysis of the classical Kaczmarz's projection method in the case of an inconsistent linear least-squares problem and we prove that the approximations so obtained are at a certain distance from the set of exact least-squares solutions. This distance is controlled by the component of the right hand side of the problem lying in the orthogonal complement of the range of problem's matrix, i.e. exactly the component that makes the problem inconsistent. For overcoming this difficulty we consider an extended version of Kaczmarz's algorithm, previously analyzed by one of the authors. In the numerical experiments described in the last part of the paper we compare the above mentioned extension with two well known (ART) type algorithms for image reconstruction in two electromagnetic geotomography problems. The results indicate that the extended Kaczmarz algorithm gives much better results than the other two.