An approach to the spherical mean Radon transform with detectors on a line

Rafik Aramyan
DOI: https://doi.org/10.48550/arXiv.1705.10690
2017-05-28
Abstract:The article suggests a new approach what is called a consistency method for the inversion of the spherical Radon transform in 2D with detectors on a line. It is known that there is not an exact inversion formula in 2D. By means of the method was proved that the reconstruction has a local description and found a new iteration formula which give an practical algorithm to recover an unknown function supported completely on one side of a line $L$ from its spherical means over circles centered on the line $L$. Such an inversion is required in problems of thermo- and photo-acoustic tomography.
Numerical Analysis
What problem does this paper attempt to address?
The problem that this paper attempts to solve is how to reconstruct an unknown function from the spherical mean Radon transform in two - dimensional space when the detector is located on a straight line. Specifically, given a function \( f \) with support on only one side, and its average value \( M_f(p, r) \) on circles centered at points on the straight line \( L \), how to recover the original function \( f \) through these average values. ### Main problems and background 1. **Problem description**: - The author studies how to recover \( f \) from its average value \( M_f(p, r) \) on circles centered at points on the straight line \( L \) in the case of a continuous real - valued function \( f \). - Here, \( M_f(p, r) \) represents the average value of the function \( f \) on a circle with center \( p\in L \) and radius \( r > 0 \). 2. **Background and motivation**: - This problem is very important in the field of medical imaging (such as thermoacoustic tomography, TAT). In thermoacoustic tomography, electromagnetic pulses pass through biological tissues, triggering thermoacoustic responses, and the generated sound waves can be measured by externally placed transducers. These measurement data are actually the integrals of the function \( f \) on different spherical surfaces. - The exact inversion formula for the spherical Radon transform is currently known only for the boundaries of some special domains (such as spheres, cylinders, and planes). For the general case, especially when the detector is located on a straight line, new methods are required to achieve inversion. ### Main contributions of the paper 1. **Proposed a new method**: - The author introduced a new method called the "consistency method" for the inversion of the spherical Radon transform in two - dimensional space. - Through this method, the author found a new iterative formula, which provides a practical algorithm for recovering the unknown function \( f \) from the spherical average values. 2. **Proved that the reconstruction has a local description**: - The author proved that the reconstruction process has local descriptiveness, that is, for any point \((x, y)\), the value of the function \( f(x, y)\) depends only on the values of \( M_f \) near \( p=(x, 0)\) and in the range of \( 0\leq r\leq y \). 3. **Provided a specific inversion formula**: - The author gave an inversion formula, which approximates the original function \( f \) through a series of standard polynomials \( Z_{n,i} \). The specific formula is as follows: \[ f(x, y)=\lim_{n\rightarrow\infty}\left(2(n + 1)M_f((x,0),y)+\sum_{i = 0}^{n}\int_{0}^{y}y^{2i-1}Z_{n,i}\left(\frac{u}{y}\right)\left(M_f((x,0),u)\right)_{x}^{(2i)}\,du\right) \] where \(\left(M_f((x,0),u)\right)_{x}^{(2i)}\) represents the \( 2i\)-th derivative with respect to the variable \( x \). ### Summary This paper solves the inversion problem of the spherical Radon transform in two - dimensional space by introducing the consistency method and provides a practical inversion formula. This result is of great significance for application fields such as thermoacoustic tomography.