Covariance Estimation in High Dimensions via Kronecker Product Expansions