Covariance Estimation for Matrix-valued Data