In optoacoustic (photoacoustic) tomography, several parameters related to tissue and detector features are needed for image formation, but they may not be known a priori. An autofocus (AF) algorithm is generally used to estimate these parameters. However, the algorithm works iteratively and is therefore impractical for clinical imaging with planar geometry systems due to the long reconstruction times. We have developed a fast autofocus (FAF) algorithm for 3D optoacoustic systems with planar geometry. Such an algorithm exploits the symmetries of the planar geometry and a virtual source concept to reduce the dimensionality of the parameter estimation problem. The dimensionality reduction makes FAF much simpler computationally than the conventional AF algorithm. We show that the FAF algorithm required about 5 s to provide accurate estimates of the speed of sound in simulated data and experimental data obtained using an imaging system that is poised to enter the clinic. The applicability of FAF for estimating other image formation parameters is discussed. We expect the FAF algorithm to contribute decisively to the clinical use of optoacoustic tomography systems with planar geometry.
