Spatial regression model

Hello,

is there a package which allows to fit a spatial regression model  for 
any k-dimensional data?
The spatial regression model is composed of three components:

y=mu+Z+epsilon,   

where mu is a low order polynomial, Z is a realization of a second-order 
stochastic process with covariance matrix Sigma and epsilon is a 
realization of a white-noise process. The covariance function belongs to 
the class of  either isotropic or anisotropic functions. Mean and 
covariance parameters are estimated simultaneously by ML or REML.

I hope somebody can answer my question.
Thank you!

Nadine Henkenjohann