Refine
Year of publication
Document Type
- Article (18)
- conference proceeding (article) (6)
- Preprint (6)
- conference proceeding (presentation, abstract) (2)
- Part of a Book (1)
- Doctoral Thesis (1)
- Report (1)
Has Fulltext
- no (35)
Is part of the Bibliography
- no (35)
Keywords
- limited angle tomography (5)
- Fourier integral operators (4)
- computed tomography (4)
- microlocal analysis (4)
- Lambda tomography (3)
- Radon transform (3)
- Curvelets (2)
- Dimensionality reduction (2)
- IMAGE-RECONSTRUCTION (2)
- Image reconstruction (2)
Institute
Begutachtungsstatus
- peer-reviewed (22)
This article provides a mathematical analysis of singular (nonsmooth) artifacts added to reconstructions by filtered backprojection (FBP) type algorithms for X-ray computed tomography (CT) with arbitrary incomplete data. We prove that these singular artifacts arise from points at the boundary of the data set. Our results show that, depending on the geometry of this boundary, two types of artifacts can arise: object-dependent and object-independent artifacts. Object-dependent artifacts are generated by singularities of the object being scanned, and these artifacts can extend along lines. They generalize the streak artifacts observed in limited-angle tomography. Object-independent artifacts, on the other hand, are essentially independent of the object and take one of two forms: streaks on lines if the boundary of the data set is not smooth at a point and curved artifacts if the boundary is smooth locally. We prove that these streak and curve artifacts are the only singular artifacts that can occur for FBP in the continuous case. In addition to the geometric description of artifacts, the article provides characterizations of their strength in Sobolev scale in certain cases. The results of this article apply to the well-known incomplete data problems, including limited-angle and regionof-interest tomography, as well as to unconventional X-ray CT imaging setups that arise in new practical applications. Reconstructions from simulated and real data are analyzed to illustrate our theorems, including the reconstruction that motivated this work a synchrotron data set in which artifacts appear on lines that have no relation to the object.
This article provides a mathematical classification of artifacts from arbitrary incom-plete X-ray tomography data when using the classical filtered backprojection algorithm. Usingmicrolocal analysis, we prove that all artifacts arise from points at the boundary of the data set.Our results show that, depending on the geometry of the data set boundary, two types of artifactscan arise: object-dependent and object-independent artifacts. The object-dependent artifacts aregenerated by singularities of the object being scanned and these artifacts can extend all along lines.This is a generalization of the streak artifacts observed in limited angle CT. The article also char-acterizes two new phenomena: the object-independent artifacts are caused only by the geometryof the data set boundary; they occur along lines if the boundary of the data set is not smooth andalong curves if the boundary of the data set is smooth. In addition to the geometric descriptionof artifacts, the article also provides characterizations of their strength in Sobolev scale in certaincases. Moreover, numerical reconstructions from simulated and real data are presented illustratingour theorems.This work is motivated by a reconstruction we present from a synchrotron data set in whichartifacts along lines appeared that were independent of the object.The results of this article apply to a wide range of well-known incomplete data problems, in-cluding limited angle CT and region of interest tomography, as well as to unconventional x-ray CTimaging setups. Some of those problems are explicitly addressed in this article, theoretically and numerically.
The characteristic feature of inverse problems is their instability with respect to data perturbations. In order to stabilize the inversion process, regularization methods have to be developed and applied. In this work we introduce and analyze the concept of filtered diagonal frame decomposition which extends the standard filtered singular value decomposition to the frame case. Frames as generalized singular system allows to better adapt to a given class of potential solutions. In this paper, we show that filtered diagonal frame decomposition yield a convergent regularization method. Moreover, we derive convergence rates under source type conditions and prove order optimality under the assumption that the considered frame is a Riesz-basis.
Translation invariant diagonal frame decomposition of inverse problems and their regularization
(2022)
Solving inverse problems is central to a variety of important applications, such as biomedical image reconstruction and non-destructive testing. These problems are characterized by the sensitivity of direct solution methods with respect to data perturbations. To stabilize the reconstruction process, regularization methods have to be employed. Well-known regularization methods are based on frame expansions, such as the wavelet-vaguelette (WVD) decomposition, which are well adapted to the underlying signal class and the forward model and allow efficient implementation. However, it is well known that the lack of translational invariance of wavelets and related systems leads to specific artifacts in the reconstruction. To overcome this problem, in this paper we introduce and analyze the concept of translation invariant diagonal frame decomposition (TI-DFD) of linear operators. We prove that a TI-DFD combined with a regularizing filter leads to a convergent regularization method with optimal convergence rates. As illustrative example, we construct a wavelet-based TI-DFD for one-dimensional integration, where we also investigate our approach numerically. The results indicate that filtered TI-DFDs eliminate the typical wavelet artifacts when using standard wavelets and provide a fast, accurate, and stable solution scheme for inverse problems.
In a variety of tomographic applications, data cannot be fully acquired, leading to severely underdetermined image reconstruction. Conventional methods result in reconstructions with significant artifacts. In order to remove these artifacts, regularization methods have to be applied that incorporate additional information. An important example is TV reconstruction which is well known to efficiently compensate for missing data and well reduces reconstruction artifacts. At the same time, however, tomographic data is also contaminated by noise, which poses an additional challenge. The use of a single regularizer within a variational regularization framework must therefore account for both the missing data and the noise. However, a single regularizer may not be ideal for both tasks. For example, the TV regularizer is a poor choice for noise reduction over different scales, in which case ℓ1 curvelet regularization methods work well. To address this issue, in this paper we introduce a novel variational regularization framework that combines the advantages of two different regularizers. The basic idea of our framework is to perform reconstruction in two stages, where the first stage mainly aims at accurate reconstruction in the presence of noise, and the second stage aims at artifact reduction. Both reconstruction stages are connected by a data proximity condition. The proposed method is implemented and tested for limited-view CT using a combined curvelet-TV approach. We define and implement a curvelet transform adapted to the limited view problem and demonstrate the advantages of our approach in a series of numerical experiments in this context.
Most of the common model-based reconstruction schemes in magnetic particle imaging (MPI) use idealized assumptions, e.g., of an ideal field-free-line (FFL) topology. However, the magnetic fields that are generated in real MPI scanners have distortions and, therefore, model-based approaches often lead to inaccurate reconstructions and may contain artifacts. In order to improve the reconstruction quality in MPI, more realistic MPI models need to be derived. In the present work, we address this problem and present a hybrid model for MPI that allows us to incorporates real measurements of the applied magnetic fields. We will explain that the measurements, that are needed to setup a model for the magnetic fields, can be obtained in a novel calibration procedure that is independent of the resolution and which is much less time-consuming than the one employed in measurement-based MPI reconstructions.We will also present a discretization strategy for this model, that can be used in context of algebraic reconstructions. The presented approach was validated on simulated data in [1], however, its evaluation on real data is a topic for future research.
In various fields of applications, inverse problems are characterized by their sensitivity to data perturbations which can cause severe reconstruction errors. Hence, regularization procedures are employed in order to ensure stability and reconstruction quality. To overcome limitations of classical approaches such as the filtered singular value decomposition (SVD), frame based
diagonalization methods have been studied in the recent years, e.g., wavelet-vagulette (WVD) decomposition. While these methods can be well adapted to the problem at hand, it is well-known, that the lack of translation invariance in multiscale systems can cause specific artifacts in the recovered object. Thus, to overcome these drawbacks we use the translation invariant diagonal frame decomposition (TI-DFD) of linear operators. For illustration, we construct a TI-WVD for one-dimensional integration operator, and confirm our theoretical findings by numerical simulations.
We study samples with full and partial occlusion causing streak artifacts, and propose two mod-ifications of filtered backprojection for artifact removal. Data is obtained by the SPring-8 synchrotron using a monochromatic parallel-beam scan [1]. Thresholding in the sinogram segments the metal, resulting in edges on which we apply 1) a smooth transition, or 2) a Dirichlet boundary condition.
We present two methods that combine image reconstruction and edge detection in computed tomography (CT) scans. Our first method is as an extension of the prominent filtered backprojection algorithm. In our second method we employ ℓ1-regularization for stable calculation of the gradient. As opposed to the first method, we show that this approach is able to compensate for undersampled CT data.