On the joint estimation of the spectra, cospectrum and quadrature spectrum of a two-dimensional stationary Gaussian process