Singular value decomposition and its application to numerical inversion for ray transforms in 2D vector tomography

摘要

The operators of longitudinal and transverse ray transforms acting on vector fields on the unit disc are considered in the paper. The goal is to construct SVD-decompositions of the operators and invert them approximately by means of truncated decomposition for the parallel scheme of data acquisition. The orthogonal bases in the initial spaces and the image spaces are constructed using harmonic, Jacobi and Gegenbauer polynomials. Based on the obtained decompositions inversion formulas are derived and the polynomial approximations for the inverse operators are obtained. Numerical tests for data sets with different noise levels of smooth and discontinuous fields show the validity of the approach for the reconstruction of solenoidal or potential parts of vector fields from their ray transforms.