Track new results across pure and applied mathematics.
Search papers, labs, and topics across Lattice.
This paper investigates the inversion of the geodesic ray transform using a finite number of local averages, proving Lipschitz stability for reconstructions within fixed finite-dimensional spaces under sufficient measurement density and geometric conditions. Convergent reconstruction algorithms are developed using Steepest Gradient Descent and Conjugate Gradients, implemented with piecewise constant finite element spaces. Numerical results in 2D and 3D highlight the impact of geometric properties and measurement discretization on reconstruction accuracy.
Finite local measurements can stably reconstruct functions from their geodesic ray transform, even in complex geometries.
We study the inversion of the geodesic ray transform from finitely many local averages of its data. Under suitable geometric assumptions, we prove Lipschitz stability on any fixed finite-dimensional reconstruction space when the measurement partition is sufficiently fine. We develop convergent reconstruction algorithms based on Steepest Gradient Descent and Conjugate Gradients and implement them using piecewise constant finite element spaces. Numerical experiments in 2D and 3D illustrate the influence of geometry and measurement discretization on reconstruction quality.