The phase-space representation of stationary random electromagnetic fields is developed by using electromagnetic spatial coherence wavelets. The propagation of the field's power and states of spatial coherence and polarization results from correlations between the components of the field vectors at pairs of points in space. Polarization domains are theoretically predicted as the structure of the field polarization at the observation plane. In addition, the phase-space representation provides a generalization of the Poynting theorem. Theoretical predictions are examined by numerically simulating the Young experiment with electromagnetic waves. The experimental implementation of these results is a current subject of research.