The wave equation and associated spherical means are a widespread model in modern imaging modalities like photoacoustic tomography. We consider a discretization for the Cauchy problem of the two dimensional wave equation by plane waves. The considered frequencies lie on a Cartesian or on a polar grid which gives rise to efficient algorithms for the computation of the spherical means. The theoretical findings are illustrated by a some numerical experiments.