Approximate joint photocount statistics for Gaussian light with arbitrary spectral shape