Published in

IOP Publishing, Inverse Problems, 1(37), p. 014002, 2020

DOI: 10.1088/1361-6420/aba5ee

Links

Tools

Export citation

Search in Google Scholar

Error analysis for filtered back projection reconstructions in Besov spaces

Journal article published in 2020 by Matthias Beckmann, Peter Maass ORCID, Judith Nickel
This paper is made freely available by the publisher.
This paper is made freely available by the publisher.

Full text: Download

Green circle
Preprint: archiving allowed
Orange circle
Postprint: archiving restricted
Red circle
Published version: archiving forbidden
Data provided by SHERPA/RoMEO

Abstract

Abstract Filtered back projection (FBP) methods are the most widely used reconstruction algorithms in computerized tomography (CT). The ill-posedness of this inverse problem allows only an approximate reconstruction for given noisy data. Studying the resulting reconstruction error has been a most active field of research in the 1990s and has recently been revived in terms of optimal filter design and estimating the FBP approximation errors in general Sobolev spaces. However, the choice of Sobolev spaces is suboptimal for characterizing typical CT reconstructions. A widely used model are sums of characteristic functions, which are better modelled in terms of Besov spaces B q α , p ( R 2 ) . In particular B 1 α , 1 ( R 2 ) with α ≈ 1 is a preferred model in image analysis for describing natural images. In case of noisy Radon data the total FBP reconstruction error ‖ f − f L δ ‖ ⩽ ‖ f − f L ‖ + ‖ f L − f L δ ‖ splits into an approximation error and a data error, where L serves as regularization parameter. In this paper, we study the approximation error of FBP reconstructions for target functions f ∈ L 1 ( R 2 ) ∩ B q α , p ( R 2 ) with positive α ∉ N and 1 ⩽ p, q ⩽ ∞. We prove that the L p -norm of the inherent FBP approximation error f − f L can be bounded above by ‖ f − f L ‖ L p ( R 2 ) ⩽ c α , q , W L − α | f | B q α , p ( R 2 ) under suitable assumptions on the utilized low-pass filter’s window function W. This then extends by classical methods to estimates for the total reconstruction error.