Estimating Basis Functions in Massive Fields under the Spatial Mixed Effects Model