Breaking the resolution limit in photoacoustic imaging using non-negativity and sparsity

The spatial resolution achievable in photoacoustic imaging decreases with the imaging depth, resulting in blurred images for deeper structures. Apart from technical limitations, the ultimate resolution limit results from the second law of thermodynamics. The attenuation of the optically generated acoustic waves on their way from the imaged structure to the sample surface by scattering and dissipation leads to an increase of entropy. The resulting loss of spatial resolution for structures embedded in attenuating media can be compensated by numerical methods that make use of additional available information. In this article, we demonstrate this using experimental data from plane one-dimensional (1D) acoustic waves propagating in fat tissue. The acoustic waves are optically induced by nanosecond laser pulses and measured with piezoelectric transducers. The experimental results of 1D compensation are also relevant for photoacoustic imaging in 2D or 3D in an acoustically attenuating medium by dividing the reconstruction problem into two steps: First, the ideal signal, which is the solution of the un-attenuated wave equation, is determined by the proposed 1D attenuation compensation for each detector signal. In a second step, any ultrasound reconstruction method for un-attenuated data can be used for image reconstruction. For the reconstruction of a small step milled into a silicon wafer surface, which allows the generation of two photoacoustic pulses with a small time offset, we take advantage of non-negativity and sparsity and inverted the measured, frequency dependent acoustic attenuation of the fat tissue. We were able to improve the spatial resolution for imaging through 20 mm of porcine fat tissue compared to the diffraction limit at the cut-off frequency by at least a factor of two.