Michael Quellmalz

Clemens Kirisits, Michael Quellmalz, Eric Setterqvist

This paper concerns diffraction-tomographic reconstruction of an object characterized by its scattering potential. We establish a rigorous generalization of the Fourier diffraction theorem in arbitrary dimension, giving a precise relation in the Fourier domain between measurements of the scattered wave and reconstructions of the scattering potential. With this theorem at hand, Fourier coverages for different experimental setups are investigated taking into account parameters such as object orientation, direction of incidence and frequency of illumination. Allowing for simultaneous and discontinuous variation of these parameters, a general filtered backpropagation formula is derived resulting in an explicit approximation of the scattering potential for a large class of experimental setups.

Robert Beinert, Jonas Bresch, Michael Quellmalz

The Radon transform is a fundamental tool for analyzing data in tomographic imaging, optimal transport, crystallography, and geometric analysis. Numerical computations require an accurate discretization. To deal with voxelized images and objects, we derive a closed-form, piecewise polynomial expression for the Radon transform of an axis-aligned cube in arbitrary dimension $d$. Building on this formula, we propose a discrete Radon transform in $\mathbb{R}^d$ that is both analytically exact for voxelized data and computationally efficient. For improved numerical stability, we introduce a regularized variant replacing the Radon transform of a cube, i.e.\ the $(d-1)$-dimensional area of the intersection between that cube and a hyperplane, by the $d$-dimensional volume of the intersection between the cube and a thin slab around the hyperplane. Numerical experiments demonstrate the effectiveness of the proposed approach in several applications including 3D shape matching, classification, and sliced Wasserstein barycenters. The computational efficiency in higher dimensions is verified by a comparison with Monte Carlo integration.

Michael Quellmalz, Mia Kvåle Løvmo, Simon Moser et al.

We investigate the application of the common circle method for estimating sample motion in optical diffraction tomography (ODT) of sub-millimeter sized biological tissue. When samples are confined via contact-free acoustical force fields, their motion must be estimated from the captured images. The common circle method identifies intersections of Ewald spheres in Fourier space to determine rotational motion. This paper presents a practical implementation, incorporating temporal consistency constraints to achieve stable reconstructions. Our results on both simulated and real-world data demonstrate that the common circle method provides a computationally efficient alternative to full optimization methods for motion detection.

Nicolaj Rux, Michael Quellmalz, Gabriele Steidl

Negative distance kernels $K(x,y) := - \|x-y\|$ were used in the definition of maximum mean discrepancies (MMDs) in statistics and lead to favorable numerical results in various applications. In particular, so-called slicing techniques for handling high-dimensional kernel summations profit from the simple parameter-free structure of the distance kernel. However, due to its non-smoothness in $x=y$, most of the classical theoretical results, e.g. on Wasserstein gradient flows of the corresponding MMD functional do not longer hold true. In this paper, we propose a new kernel which keeps the favorable properties of the negative distance kernel as being conditionally positive definite of order one with a nearly linear increase towards infinity and a simple slicing structure, but is Lipschitz differentiable now. Our construction is based on a simple 1D smoothing procedure of the absolute value function followed by a Riemann-Liouville fractional integral transform. Numerical results demonstrate that the new kernel performs similarly well as the negative distance kernel in gradient descent methods, but now with theoretical guarantees.

Michael Quellmalz

The Funk-Radon transform, also known as the spherical Radon transform, assigns to a function on the sphere its mean values along all great circles. Since its invention by Paul Funk in 1911, the Funk-Radon transform has been generalized to other families of circles as well as to higher dimensions. We are particularly interested in the following generalization: we consider the intersections of the sphere with hyperplanes containing a common point inside the sphere. If this point is the origin, this is the same as the aforementioned Funk--Radon transform. We give an injectivity result and a range characterization of this generalized Radon transform by finding a relation with the classical Funk--Radon transform.

Ralf Hielscher, Daniel Potts, Michael Quellmalz

In spherical surface wave tomography, one measures the integrals of a function defined on the sphere along great circle arcs. This forms a generalization of the Funk--Radon transform, which assigns to a function its integrals along full great circles. We show a singular value decomposition (SVD) for the surface wave tomography provided we have full data. Since the inversion problem is overdetermined, we consider some special cases in which we only know the integrals along certain arcs. For the case of great circle arcs with fixed opening angle, we also obtain an SVD that implies the injectivity, generalizing a previous result for half circles in [Groemer, On a spherical integral transform and sections of star bodies, Monatsh. Math., 126(2):117--124, 1998]. Furthermore, we derive a numerical algorithm based on the SVD and illustrate its merchantability by numerical tests.