Modeling and prediction for multivariate spatial factor analysis