Numerical Simulation of Non-Gaussian Random Fields with Prescribed   Marginal Distributions and Cross-Correlation Structure II: Multivariate   Random Fields
  We provide theoretical procedures and practical recipes to simulate non-Gaussian correlated, homogeneous random fields with prescribed marginal distributions and cross-correlation structure, either in a N-dimensional Cartesian space or on the celestial sphere. We illustrate our methods using far-infrared maps obtained with the Infrared Space Observatory. However, the methodology presented here can be used in other astrophysical applications that require modeling correlated features in sky maps, for example, the simulation of multifrequency sky maps where backgrounds, sources and noise are correlated and can be modeled by random fields. 
R. Vio, P. Andreani, L. Tenorio & W. Wamsteker