As most georeferenced data sets are multivariate and concern variables of different kinds, spatial mapping methods must be able to deal with such data. The main difficulties are the prediction of non-Gaussian variables and the dependence modelling between processes. The aim of this paper is to present a new approach that permits simultaneous modelling of Gaussian, count and ordinal spatial processes. We consider a hierarchical model implemented within a Bayesian framework. The method used for Gaussian and count variables is based on the generalized linear model. Ordinal variable is taken into account through a generalization of the ordinal probit model. We use the moving average approach of Ver Hoef and Barry to model the dependencies between the processes.