An improved exact inversion formula for solenoidal fields in cone beam vector tomography

Abstract

In this paper we present an improved inversion formula for the 3D cone beam transform of vector fields supported in the unit ball which is exact for solenoidal fields. It is well known that only the solenoidal part of a vector field can be determined from the longitudinal ray transform of a vector field in cone beam geometry. The inversion formula, as it was developed in Katsevich and Schuster (2013 An exact inversion formula for cone beam vector tomography Inverse Problems 29 065013), consists of two parts. The first part is of the filtered backprojection type, whereas the second part is a costly 4D integration and very inefficient. In this article we tackle this second term and obtain an improved formula, which is easy to implement and saves one order of integration. We also show that the first part contains all information about the curl of the field, whereas the second part has information about the boundary values. More precisely, the second part vanishes if the solenoidal part of the original field is tangential at the boundary. A number of numerical tests presented in the paper confirm the theoretical results and the exactness of the formula. Also, we obtain an inversion algorithm that works for general convex domains.