Solving or resolving inadequate and noisy tomographic systems (Q1074329)
From MaRDI portal
scientific article
Language | Label | Description | Also known as |
---|---|---|---|
English | Solving or resolving inadequate and noisy tomographic systems |
scientific article |
Statements
Solving or resolving inadequate and noisy tomographic systems (English)
0 references
1985
0 references
The determination of seismic velocity leads to the inversion of the Radon transform; i.e., line integrals of the searched-for distribution are given. In contrast to medicals application the data cannot be collected uniformly, resulting in a highly ill-posed problem. Also the recovery is complicated by very noisy data. The projection of the solution onto a finite dimensional pixel space then leads to very ill-conditioned linear equations. The author compares different stabilization methods, namely singular value decomposition (SVD), conjugate gradients (CG) and the Dines-Lytle method, which is a version of the Landweber iteration. It turns out that the CG-method is cheaper than SVD and that it picks up the recoverable part of the solution after a few steps. The Dines-Lytle method is much slower. In a resolution analysis the trade-off between resolution and noise damping is studied.
0 references
tomography
0 references
regularisation by iterative methods
0 references
seismology
0 references
inversion
0 references
Radon transform
0 references
ill-posed problem
0 references
noisy data
0 references
stabilization methods
0 references
singular value decomposition
0 references
conjugate gradients
0 references
Dines-Lytle method
0 references
Landweber iteration
0 references