Abstract Modern approaches for the spatial simulation of categorical variables are largely based on multi-point statistical methods, where a training image is used to derive complex relationships using relevant patterns. In these approaches, simulated realizations driven by utilized, while statistics actual sample data ignored. This paper presents data-driven, high-order approach approximation ...