Slicing of Radial Functions: a Dimension Walk in the Fourier Space
Nicolaj Rux, Michael Quellmalz, Gabriele Steidl
TL;DR
The paper addresses recovering a univariate radial slice $f$ from a radial high-dimensional function $F$ linked by $F(|x|)=\mathbb{E}_{\xi}[f(|\langle x,\xi\rangle|)]$. It develops a Fourier-space dimension-walk using the rotation $\mathcal{R}_d$, averaging $\mathcal{A}_d$, and multiplication $\mathcal{M}_d$ operators to reduce from $d$ to 1D, and extends the theory from $L^1$ radial functions to tempered distributions and positive definite kernels via Bochner’s theorem. The main contributions include explicit inversion formulas (e.g., $f=\tfrac{\omega_{d-1}}{2}(\mathcal{F}_1\circ\mathcal{M}_d)[\rho]$ with $\rho$ tied to $F$), a distributional slicing framework for radial tempered distributions, and a positive definite correspondence between radial $F$ and $f$. These results enable rigorous, efficient handling of high-dimensional radial convolutions in applications while providing theoretical guarantees on smoothness and positivity properties. The work integrates Abel-type relations, Radon-adjoint connections, and fractional derivative perspectives to unify both function and distributional settings.
Abstract
Computations in high-dimensional spaces can often be realized only approximately, using a certain number of projections onto lower dimensional subspaces or sampling from distributions. In this paper, we are interested in pairs of real-valued functions $(F,f)$ on $[0,\infty)$ that are related by the projection/slicing formula $F (\| x \|) = \mathbb E_ξ \big[ f \big(|\langle x,ξ\rangle| \big) \big]$ for $x\in\mathbb R^d$, where the expectation value is taken over uniformly distributed directions in $\mathbb R^d$. While it is known that $F$ can be obtained from $f$ by an Abel-like integral formula, we construct conversely $f$ from given $F$ using their Fourier transforms. First, we consider the relation between $F$ and $f$ for radial functions $F(\| \cdot\| )$ that are Fourier transforms of $L^1$ functions. Besides $d$- and one-dimensional Fourier transforms, it relies on a rotation operator, an averaging operator and a multiplication operator to manage the walk from $d$ to one dimension in the Fourier space. Then, we generalize the results to tempered distributions, where we are mainly interested in radial regular tempered distributions. Based on Bochner's theorem, this includes positive definite functions $F(\| \cdot\| )$ and, by the theory of fractional derivatives, also functions $F$ whose derivative of order $\lfloor d/2\rfloor$ is slowly increasing and continuous.
