In this paper we investigate systematic sampling in the image- synthesis context. Systematic sampling has been widely used in stereology to improve the efficiency of different probes in experimental design. These designs are theoretically based on estimators of 1-dimensional and 2-dimensional integrals. For the particular case of the characteristic function, the variance of these estimators has been shown to be asymptotically N −−3/2, which improves on the O(N −−1) behaviour of independent estimators using uniform sampling. Thus, when no a priori knowledge of the integrand function is available, like in several image synthesis techniques, systematic sampling efficiently reduces the computational cost.