Abstract:Recovering a function from its integrals over circular cones recently gained significance because of its relevance to novel medical imaging technologies such emission tomography using Compton cameras. In this paper we investigate the case where the vertices of the cones of integration are restricted to a sphere in $n$-dimensional space and symmetry axes are orthogonal to the sphere. We show invertibility of the considered transform and develop an inversion method based on series expansion and reduction to a system of one-dimensional integral equations of generalized Abel type. Because the arising kernels do not satisfy standard assumptions, we also develop a uniqueness result for generalized Abel equations where the kernel has zeros on the diagonal. Finally, we demonstrate how to numerically implement our inversion method and present numerical results.
What problem does this paper attempt to address?
### What problem does this paper attempt to solve?
This paper aims to solve the problem of how to recover functions from conical integrals, especially when the vertex is located on the sphere and the axis of symmetry is orthogonal to the sphere. Specifically, the paper studies the invertibility of the **conical Radon transform** and develops an inversion method based on series expansion and reduction to a one - dimensional integral equation of the generalized Abel type.
#### Background and Motivation
1. **Medical imaging technology**: In particular, single - photon emission computed tomography (SPECT) using a Compton camera, which can significantly increase the number of photons recorded, thereby reducing the noise level.
2. **Data reconstruction challenges**: Traditional SPECT uses a collimator, which causes most photons to be blocked, so a new method is needed to reconstruct images from the data obtained by the Compton camera.
#### Main Contributions
1. **Theoretical analysis**:
- Studies the invertibility of the conical Radon transform when the vertex is located on the sphere and the axis of symmetry is orthogonal to the sphere.
- Develops a method based on spherical harmonic function decomposition, which simplifies the conical Radon transform into a set of one - dimensional integral equations given explicitly.
2. **Mathematical derivation**:
- Derives the uniqueness results of the generalized Abel equation, even when the kernel function has zeros on the diagonal.
- Proves the invertibility of the considered transform and gives an effective inversion method.
3. **Numerical implementation**:
- Proposes a numerical algorithm based on the discrete Fourier transform (FFT) and regularization methods (such as Tikhonov regularization).
- Verifies the effectiveness of the algorithm through simulation experiments and shows the results of image reconstruction from simulated data.
#### Formula Representation
- The conical Radon transform is defined as:
\[
R_m f(z, \psi) = \int_{C(z, \psi)} f(x) \|x - z\|^m dS(x)
\]
where \( C(z, \psi) = \{ z + r\omega \mid r \geq 0, \omega \in S^{n - 1}, -\omega \cdot z = \cos(\psi) \} \) is the right conical surface with vertex \( z \), central axis \( -z \), and half - open angle \( \psi \).
- The generalized Abel equation has the form:
\[
(R_m f)_{\ell,k}(\psi) = |S^{n - 2}| \sin(\psi)^{-m} \int_{\sin(\psi)}^1 f_{\ell,k}(\rho) \rho K_\ell(\psi, \rho) \sqrt{\rho^2 - \sin(\psi)^2} d\rho
\]
where \( K_\ell(\psi, \rho) \) is the kernel function.
#### Conclusion
The paper successfully solves the problem of recovering functions from conical integrals under specific geometric conditions, provides a theoretical basis and numerical methods, and lays the foundation for further research and practical applications.