Source author record

Jürgen Frikel

Jürgen Frikel appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

7works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

7 published item(s)

preprint2022arXiv

Regularization of Inverse Problems by Filtered Diagonal Frame Decomposition

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.

preprint2015arXiv

Artifacts in incomplete data tomography - with applications to photoacoustic tomography and sonar

We develop a paradigm using microlocal analysis that allows one to characterize the visible and added singularities in a broad range of incomplete data tomography problems. We give precise characterizations for photo- and thermoacoustic tomography and Sonar, and provide artifact reduction strategies. In particular, our theorems show that it is better to arrange Sonar detectors so that the boundary of the set of detectors does not have corners and is smooth. To illustrate our results, we provide reconstructions from synthetic spherical mean data as well as from experimental photoacoustic data.

preprint2015arXiv

Joint Image Reconstruction and Segmentation Using the Potts Model

We propose a new algorithmic approach to the non-smooth and non-convex Potts problem (also called piecewise-constant Mumford-Shah problem) for inverse imaging problems. We derive a suitable splitting into specific subproblems that can all be solved efficiently. Our method does not require a priori knowledge on the gray levels nor on the number of segments of the reconstruction. Further, it avoids anisotropic artifacts such as geometric staircasing. We demonstrate the suitability of our method for joint image reconstruction and segmentation. We focus on Radon data, where we in particular consider limited data situations. For instance, our method is able to recover all segments of the Shepp-Logan phantom from $7$ angular views only. We illustrate the practical applicability on a real PET dataset. As further applications, we consider spherical Radon data as well as blurred data.

preprint2015arXiv

Limited data problems for the generalized Radon transform in $\mathbb{R}^n$

We consider the generalized Radon transform (defined in terms of smooth weight functions) on hyperplanes in $\mathbb{R}^n$. We analyze general filtered backprojection type reconstruction methods for limited data with filters given by general pseudodifferential operators. We provide microlocal characterizations of visible and added singularities in $\mathbb{R}^n$ and define modified versions of reconstruction operators that do not generate added artifacts. We calculate the symbol of our general reconstruction operators as pseudodifferential operators, and provide conditions for the filters under which the reconstruction operators are elliptic for the visible singularities. If the filters are chosen according to those conditions, we show that almost all visible singularities can be recovered reliably. Our work generalizes the results for the classical line transforms in $\mathbb{R}^2$ and the classical reconstruction operators (that use specific filters). In our proofs, we employ a general paradigm that is based on the calculus of Fourier integral operators. Since this technique does not rely on explicit expressions of the reconstruction operators, it enables us to analyze more general imaging situations.

preprint2015arXiv

On Artifacts in Limited Data Spherical Radon Transform: Curved Observation Surface

In this article, we consider the limited data problem for spherical mean transform. We characterize the generation and strength of the artifacts in a reconstruction formula. In contrast to the third's author work [Ngu15b], the observation surface considered in this article is not flat. Our results are comparable to those obtained in [Ngu15b] for flat observation surface. For the two dimensional problem, we show that the artifacts are $k$ orders smoother than the original singularities, where $k$ is vanishing order of the smoothing function. Moreover, if the original singularity is conormal, then the artifacts are $k+\frac{1}{2}$ order smoother than the original singularity. We provide some numerical examples and discuss how the smoothing effects the artifacts visually. For three dimensional case, although the result is similar to that [Ngu15b], the proof is significantly different. We introduce a new idea of lifting the space.

preprint2014arXiv

A paradigm for the characterization of artifacts in tomography

We present a paradigm for characterization of artifacts in limited data tomography problems. In particular, we use this paradigm to characterize artifacts that are generated in reconstructions from limited angle data with generalized Radon transforms and general filtered backprojection type operators. In order to find when visible singularities are imaged, we calculate the symbol of our reconstruction operator as a pseudodifferential operator.

preprint2011arXiv

Sparse regularization in limited angle tomography

We investigate the reconstruction problem of limited angle tomography. Such problems arise naturally in applications like digital breast tomosynthesis, dental tomography, electron microscopy etc. Since the acquired tomographic data is highly incomplete, the reconstruction problem is severely ill-posed and the traditional reconstruction methods, such as filtered backprojection (FBP), do not perform well in such situations. To stabilize the reconstruction procedure additional prior knowledge about the unknown object has to be integrated into the reconstruction process. In this work, we propose the use of the sparse regularization technique in combination with curvelets. We argue that this technique gives rise to an edge-preserving reconstruction. Moreover, we show that the dimension of the problem can be significantly reduced in the curvelet domain. To this end, we give a characterization of the kernel of limited angle Radon transform in terms of curvelets and derive a characterization of solutions obtained through curvelet sparse regularization. In numerical experiments, we will present the practical relevance of these results.