  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.
