We study the limited data problem of the spherical Radon transform in two and three-dimensional spaces with general acquisition surfaces. In such situations, it is known that the application of filtered-backprojection reconstruction formulas might generate added artifacts and degrade the quality of reconstructions. In this article, we explicitly analyze a family of such inversion formulas, depending on a smoothing function that vanishes to order k on the boundary of the acquisition surfaces. We show that the artifacts are k orders smoother than their generating singularity. Moreover, in two-dimensional space, if the generating singularity is conormal satisfying a generic condition then the artifacts are even orders smoother than the generating singularity. Our analysis for three-dimensional space contains an important idea of lifting up space. We also explore the theoretical findings in a series of numerical experiments. Our experiments show that a good choice of the smoothing function leads to a significant improvement of reconstruction quality.

To accelerate photoacoustic data acquisition, in [R. Nuster, G. Zangerl, M. Haltmeier, G. Paltauf (2010). Full field detection in photoacoustic tomography. Optics express, 18(6), 6288–6299] a novel measurement and reconstruction approach has been proposed, where the measured data consist of projections of the full 3D acoustic pressure distribution at a certain time instant T. Existing reconstruction algorithms for this kind of setup assume a constant speed of sound. This assumption is not always met in practice and thus can lead to erroneous reconstructions. In this paper, we present a two-step reconstruction method for full field detection photoacoustic tomography that takes variable speed of sound into account. In the first step, by applying the inverse Radon transform, the pressure distribution at the measurement time is reconstructed point-wise from the projection data. In the second step, a final time wave inversion problem is solved where the initial pressure distribution is recovered from the known pressure distribution at time T. We derive an iterative solution approach for the final time wave inversion problem and compute the required adjoint operator. Moreover, as the main result of this paper, we derive its uniqueness and stability. Our numerical results demonstrate that the proposed reconstruction scheme is fast and stable, and that ignoring sound speed variations significantly degrades the reconstruction.

Photoacoustic tomography (PAT) is a noninvasive imaging modality that requires recovering the initial data of the wave equation from certain measurements of the solution outside the object. In the standard PAT measurement setup, the used data consist of time-dependent signals measured on an observation surface. In contrast, the measured data from the recently invented full-field detection technique provide the solution of the wave equation on a spatial domain at a single instant in time. While reconstruction using classical PAT data has been extensively studied, not much is known for the full-field PAT problem. In this paper, we build mathematical foundations of the latter problem for variable sound speed and settle its uniqueness and stability. Moreover, we introduce an exact inversion method using time-reversal and study its convergence. Our results demonstrate the suitability of both the full-field approach and the proposed time-reversal technique for high-resolution photoacoustic imaging.