  In most photoacoustic (PA) measurements, variations in speed-of-sound (SOS)
of the subject are neglected under the assumption of acoustic homogeneity.
Biological tissue with spatially heterogeneous SOS cannot be accurately
reconstructed under this assumption. We present experimental and image
reconstruction methods with which 2-D SOS distributions can be accurately
acquired and reconstructed, and with which the SOS map can be used subsequently
to reconstruct highly accurate PA tomograms. We begin with a 2-D iterative
reconstruction approach in an ultrasound transmission tomography (UTT) setting,
which uses ray refracted paths instead of straight ray paths to recover
accurate SOS images of the subject. Subsequently, we use the SOS distribution
in a new 2-D iterative approach, where refraction of rays originating from PA
sources are accounted for in accurately retrieving the distribution of these
sources. Both the SOS reconstruction and SOS-compensated PA reconstruction
methods utilize the Eikonal equation to model acoustic wavefront propagation,
which is solved using a high accuracy fast marching method (HAFMM). We validate
the new reconstruction algorithms using numerical phantoms. For experiments we
use the PER-PACT method which can be used to simultaneously acquire SOS and PA
data from subjects. We test the reconstruction algorithms using experimental
data acquired with the PER-PACT setup from challenging physical phantoms. The
results show that it is important to take SOS inhomogeneities into account. The
iterative reconstruction algorithms, that model acoustic refractive effects,
yield artifact-free highly accurate images. Our approach of using the hybrid
measurement method and the new reconstruction algorithms, is successful in
substantially improving the quality of PA images with a minimization of
blurring and artefacts.
