Resolution analysis of inverting the generalized N-dimensional Radon transform in \mathbb R^n from discrete data

DOI10.1007/S00041-022-09975-XarXiv2102.09035MaRDI QIDQ6360904FDOQ6360904

Authors: A. Katsevich 

Publication date: 17 February 2021

Abstract: Let $mathcalR$ denote the generalized Radon transform (GRT), which integrates over a family of $N$-dimensional smooth submanifolds $mathcal{S}_{ildey}subsetmathcalU$, $1leNlen-1$, where an open set $mathcalUsubsetmathbb{R}^n$ is the image domain. The submanifolds are parametrized by points $ildeysubsetildemathcalV$, where an open set $ildemathcalVsubsetmathbb{R}^n$ is the data domain. The continuous data are $g=mathcalRf$, and the reconstruction is $checkf=mathcal{R}^{\ast }mathcalBg$. Here $mathcal{R}^{\ast }$ is a weighted adjoint of $mathcalR$, and $mathcalB$ is a pseudo-differential operator. We assume that $f$ is a conormal distribution, $extsupp(f)subsetmathcalU$, and its singular support is a smooth hypersurface $mathcalSsubsetmathcalU$. Discrete data consists of the values of $g$ on a lattice $ilde{y}^j$ with the step size $O(epsilon)$. Let $check{f}_{e}psilon=mathcal{R}^{\ast }mathcalB{g}_{e}psilon$ denote the reconstruction obtained by applying the inversion formula to an interpolated discrete data $g_{e}psilon(ildey)$. Pick a generic pair $(x_0,ilde{y}_0)$, where $x_{0}inmathcalS$, and $mathcal{S}_{ilde{y}_0}$ is tangent to $mathcalS$ at $x_0$. The main result of the paper is the computation of the limit ${\displaystyle }$ f_0(check x):=lim_{epsilon o0}epsilon^kappa check f_epsilon(x_0+epsiloncheck x). ${\displaystyle }$ Here $kappage0$ is selected based on the strength of the reconstructed singularity, and $checkx$ is confined to a bounded set. The limiting function $f_0(checkx)$, which we call the discrete transition behavior, allows computing the resolution of reconstruction.

Mathematics Subject Classification ID

Radon transform (44A12) Numerical methods for integral transforms (65R10)

This page was built for publication: Resolution analysis of inverting the generalized $N$-dimensional Radon transform in $\mathbb R^n$ from discrete data

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q6360904)