Inference with Multivariate GMRFs on Regular Lattices with General Boundary Conditions