Search

Photoacoustic computed tomography (PACT), also known as optoacoustic tomography, is an emerging imaging modality that has great potential for a wide range of biomedical imaging applications. In this note, we derive a hybrid reconstruction formula that is mathematically exact and operates on a data function that is expressed in the temporal frequency and spatial domains. This formula explicitly reveals new insights into how the spatial frequency components of the sought-after object function are determined by the temporal frequency components of the data functionmeasured with a circular or spherical measurement geometry in two- and three-dimensional implementations of PACT, respectively. The structure of the reconstruction formula is surprisingly simple compared with existing Fourier-domain reconstruction formulae. It also yields a straightforward numerical implementation that is robust and two orders ofmagnitude more computationally efficient than filtered backprojection algorithms. [ABSTRACT FROM AUTHOR]

This paper deals with image restoration problems where the data are nonuniform samples of the Fourier transform of the unknown object. We study the inverse problem in both semidiscrete and fully discrete formulations, and our analysis leads to an optimization problem involving the minimization of the data discrepancy under nonnegativity constraints. In particular, we show that such a problem is equivalent to a deconvolution problem in the image space. We propose a practical algorithm, based on the gradient projection method, to compute a regularized solution in the discrete case. The key point in our deconvolution-based approach is that the fast Fourier transform can be employed in the algorithm implementation without the need of preprocessing the data. A numerical experimentation on simulated and real data from the NASA RHESSI mission is also performed. [ABSTRACT FROM AUTHOR]