Iterative image reconstruction in photoacoustic tomography using Kaiser-Bessel windows

Iterative image reconstruction algorithms for photoacoustic computed tomography (PACT) have the ability to improve image quality over analytic algorithms due to their ability to mitigate artifacts from incomplete data, incorporate the relevant imaging physics, and model the instrument response. In this work, Kaiser-Bessel functions are employed as the basis functions in an iterative reconstruction algorithm for PACT. Kaiser-Bessel functions, or blobs, can be made arbitrarily smooth (differentiable to arbitrary order), have finite spatial support, and can be made quasi-bandlimited in the spatial Fourier domain. Closed-form solutions exist in the time-domain and the temporal-frequency domain for the pressure signal generated by blobs.