Efficient Computation of Blue Noise Point Sets through Importance Sampling

Dart‐throwing can generate ideal Poisson‐disk distributions with excellent blue noise properties, but is very computationally expensive if a maximal point set is desired. In this paper, we observe that the Poisson‐disk sampling problem can be posed in terms of importance sampling by representing the available space to be sampled as a probability density function (pdf). This allows us to develop an efficient algorithm for the generation of maximal Poisson‐disk distributions with quality similar to naïve dart‐throwing but without rejection of samples. In our algorithm, we first position samples in one dimension based on its marginal cumulative distribution function (cdf). We then throw samples in the other dimension only in the regions which are available for sampling. After each 2D sample is placed, we update the cdf and data structures to keep track of the available regions. In addition to uniform sampling, our method is able to perform variable‐density sampling with small modifications. Finally, we also propose a new min‐conflict metric for variable‐density sampling which results in better adaptation of samples to the underlying importance field.

[1]  Li-Yi Wei,et al.  Parallel Poisson disk sampling , 2008, ACM Trans. Graph..

[2]  Eugene Fiume,et al.  Hierarchical Poisson disk sampling distributions , 1992 .

[3]  Thouis R. Jones Efficient Generation of Poisson-Disk Sampling Patterns , 2006, J. Graph. Tools.

[4]  Greg Humphreys,et al.  A spatial data structure for fast Poisson-disk sample generation , 2006, ACM Trans. Graph..

[5]  Ares Lagae,et al.  A Comparison of Methods for Generating Poisson Disk Distributions , 2008, Comput. Graph. Forum.

[6]  S. P. Lloyd,et al.  Least squares quantization in PCM , 1982, IEEE Trans. Inf. Theory.

[7]  Robert L. Cook,et al.  Stochastic sampling in computer graphics , 1988, TOGS.

[8]  Peter Shirley Nonuniform Random Point Sets via Warping , 1992, Graphics Gems III.

[9]  Oliver Deussen,et al.  Wang Tiles for image and texture generation , 2003, ACM Trans. Graph..

[10]  Michael Balzer,et al.  Capacity-constrained point distributions: a variant of Lloyd's method , 2009, ACM Trans. Graph..

[11]  Li-Yi Wei Parallel Poisson disk sampling , 2008, SIGGRAPH 2008.

[12]  Adrian Secord,et al.  Weighted Voronoi stippling , 2002, NPAR '02.

[13]  Robert Ulichney,et al.  Digital Halftoning , 1987 .

[14]  Dani Lischinski,et al.  Recursive Wang tiles for real-time blue noise , 2006, ACM Trans. Graph..

[15]  Ares Lagae,et al.  A procedural object distribution function , 2005, TOGS.

[16]  O. Deussen,et al.  Capacity-constrained point distributions: a variant of Lloyd's method , 2009, SIGGRAPH 2009.

[17]  V. Ostromoukhov,et al.  Fast hierarchical importance sampling with blue noise properties , 2004, SIGGRAPH 2004.

[18]  Alexander Keller,et al.  Tiled Blue Noise Samples , 2001, VMV.

[19]  Chi-Wing Fu,et al.  Anisotropic blue noise sampling , 2010, ACM Trans. Graph..

[20]  Steve C. Maddock,et al.  Accurate multidimensional Poisson-disk sampling , 2009, TOGS.

[21]  Mark A. Z. Dippé,et al.  Antialiasing through stochastic sampling , 1985, SIGGRAPH.

[22]  K.B. White,et al.  Poisson Disk Point Sets by Hierarchical Dart Throwing , 2007, 2007 IEEE Symposium on Interactive Ray Tracing.

[23]  Anton Alstes Wang Tiles for Image and Texture Generation , 2004 .