A new recovery operator P : Q(n)(disc)(T) -> Q(n+1)(disc)(M) for discontinuous Galerkin is derived. It is based on the idea of projecting a discontinuous, piecewise polynomial solution on a given mesh T into a higher order polynomial space on a macro mesh M. In order to do so, we define local degrees of freedom using polynomial moments and provide global degrees of freedom on the macro mesh. We prove consistency with respect to the local L-2-projection, stability results in several norms and optimal anisotropic error estimates. As an example, we apply this new recovery technique to a stabilized solution of a singularly perturbed convection-diffusion problem using bilinear elements.