Latent Gaussian random fields for large discrete areal data