A method based on the maximum a posteriori probability (MAP) criterion is proposed to estimate the channel frequency response (CFR) matrix and interference- plus-noise spatial covariance matrix (SCM) for multiple input and multiple output orthogonal frequency division multiplexing (MIMO-OFDM) systems. An iterative solution is proposed to solve the MAP-based problem and an interference rejection combining (IRC) receiver is derived to suppress co-channel interference (CCI) based on the estimated CFR and SCM. Furthermore, considering the property of SCM, i. e., Hermitian and semi-definite, two schemes are proposed to improve the accuracy of SCM estimation. The first scheme is proposed to parameterize the SCM via a sum of a series of matrices in the time domain. The second scheme measures the SCM on each subcarrier as a low-rank model while the model order can be chosen through the penalized-likelihood approach. Simulation results are provided to demonstrate the effectiveness of the proposed method.