Sparse Covariance Function Estimation in High Dimensions