Spherical means are a widespread model in modern imaging modalities like photoacoustic tomography. We develop Fourier based algorithms for an efficient computation of mean values. Furthermore we consider iterative reconstruction schemes, where we employ different regularization techniques.