1 Introduction
Nonperiodic tilings have emerged as an important mathematical tool in a variety of computer graphics applications [28]. They have proven especially useful in the design of sampling algorithms, where they serve to direct the spatial distribution of rendering primitives by enforcing spatial uniformity while precluding regular repetition. Recently, Wang tilings [4, 20, 24, 25, 26], Penrose tilings [38], Socolar tilings [39] and polyominoes [40] have been used to generate point sets for nonperiodic sampling. In one of the earliest applications of nonperiodic tilings in computer graphics, Penrose tilings [13, 15, 45] were employed by RangelMondragon and Abas [43] in the design of decorative patterns inspired by Islamic art. They had effectively reinvented the medieval trade secrets of the craftsmen of fifteenth century Islamic mosques [29] who created by hand highly intricate mosaics closely resembling quasicrystal tilings only discovered by modern science in the late twentieth century. Wang tilings [13, 15] were first introduced by Stam [49] in order to enable wave texture patches to cover water surfaces of arbitrary size without the appearance of regularly repeating artifacts. Further computer graphics applications of nonperiodic tilings include texture mapping and synthesis [4, 24, 25, 26, 49, 50], photorealistic rendering using environmental maps [25, 38], and nonphotorealistic rendering using stippling [24, 25]. For computer graphics applications, nonperiodic tilings have usually been generated by geometric constructs, such as matching rules and hierarchical substitution [12, 15]. In this work, we present the cutandproject method of generating quasicrystals as an alternative algebraic approach to the production of nonperiodic tilings and point sets (Fig. 1).
Point set  Voronoi diagram 
Delaunay graph  Delaunay triangulation 
This algebraic approach has the advantages of being straightforward to implement, easy to calculate, and readily amenable to rigorous mathematical analysis. Moreover, it may be directly extended to higher dimensions as well as adaptive sampling applications, although this is outside the scope of our present work. We choose to base our method on the algebra of the golden ratio , as its geometrical properties have been previously successfully exploited in computer graphics in the context of spatial sampling [38] and color quantization [34] techniques that rely on the Fibonacci number system. For an introduction to the theory of quasicrystals, consult Senechal’s comprehensive textbook [45], while a more advanced treatment of the cutandproject method may be found in surveys by Patera [41] as well as Chen, Moody, and Patera [3].
In the evaluation of the effectiveness of quasicrystals as a nonadaptive image sampling strategy, our work is motivated by the use of image sampling in multiresolution image representation and progressive image rendering. In particular, we base our experimental investigations on our experience with the development of a pointbased rendering approach to multiresolution image representation for digital photography [16, 17]
based on scattered data interpolation techniques
[2], which has been shown to support a secure and compact image encoding suitable for both photorealistic image reconstruction and nonphotorealistic image rendering. A thorough discussion of the standard image sampling strategies can be found in Glassner’s textbook [11]. Their effectiveness has been extensively investigated for use in photorealistic computer graphics applications [48], such as Monte Carlo integration in 3D ray tracing. For nonadaptive sampling, the key tradeoff is between aliasing and noise, as exemplified by the regular structure of periodic sampling using a square grid and the irregular clustering of random sampling using a uniform distribution. The classic compromise strategies are jittered sampling
[5, 7, 23], which disrupts the regularity of periodic sampling by randomly perturbing the sample sites, and quasirandom sampling [42], which avoids the irregularity of random sampling by ensuring a consistent density of sampling is maintained. In our experiments, we demonstrate that quasicrystal sampling can permit more accurate photorealistic image reconstruction than either standard jittered sampling or standard quasirandom sampling. For photorealistic image reconstruction [11], the ideal strategy is generally considered to be the Poisson disk distribution [5, 7], random point sets conditioned on a minimum distance between the points, while for nonphotorealistic image rendering [6, 18, 19, 21, 44], a popular strategy relies on centroidal Voronoi diagrams [8], optimized point sets with every point placed at the centroid of its Voronoi polygon. As both of these sampling strategies prove time consuming to compute exactly, a variety of approximation techniques have been proposed to produce sampling patterns that have similar properties in the frequency domain, in particular those that exhibit a blue noise Fourier power spectrum characteristic of a Poisson disk distribution. Historically, blue noise sampling strategies relied on slow, trial and error, stochastic procedures involving dart throwing algorithms that approximate the Poisson disk distribution by rejecting prospective locations for new sample sites whenever they are deemed to be too close to the preceding samples [5, 7, 31, 33]. Improved performance of blue noise sampling can be obtained through the use of efficient geometric data structures [9, 10, 22, 52] and parallel processing GPU hardware [51]. Alternatively, one can readily generate a blue noise sampling pattern using a nonperiodic tiling composed of a suitable set of tiles, where each tile contains a precomputed optimal arrangement of sample sites [4, 20, 24, 25, 26, 38, 40]. A detailed evaluation of the spectral properties of various blue noise sampling algorithms can be found in the survey by Lagae and Dutre [27]. In practical applications, it can often be quite difficult to visually distinguish between renditions produced by blue noise sampling patterns generated using different algorithms. Hence, for the purpose of our evaluation, we relied on farthest point sampling [10] since in previous work we have shown this blue noise sampling technique to be highly suitable for multiresolution image representation [17]. In general, blue noise sampling algorithms tend to have higher requirements for either computational processing, data storage, or implementation complexity than simpler sampling strategies, such as quasicrystal sampling. Naturally, simpler sampling strategies cannot replicate all of the desirable qualities of blue noise sampling. Nevertheless, as demonstrated by our evaluation, quasicrystal sampling is shown to be proficient at supporting a uniform sampling density, centroidal Voronoi regions, accurate image reconstruction, and progressive image rendering, despite having only a small number of local sampling configurations arranged in an anisotropic manner incompatible with a blue noise Fourier power spectrum. Hence, as a potential alternative to periodic sampling in image representation, sampling using cutandproject quasicrystals can deterministically guarantee minimum and maximum distances between nearest neighbors in a uniformly spacefilling sampling pattern without the overhead of geometric data structures or tiling lookup tables for tracking their placement.2 Method
While a periodic point set is characterized by its translational symmetries, a nonperiodic point set admits no translational symmetries. For use in image sampling, we focus on nonperiodic point sets that are determined by their inflation symmetries. In such a nonperiodic point set, a fixed configuration of sample sites can be repeated at different scales to generate a selfsimilar pattern. The simplest way of producing nonperiodic point sets is to use hierarchical substitution tilings [14, 15]. For instance, hierarchical substitution can be readily applied to the famous Penrose tiling [13, 15, 38, 43, 45]. The strategy starts with a small set of polygonal tiles. The tile set is carefully designed so that each tile can be decomposed by geometric subdivision into smaller instances of itself and the other tiles. A hierarchy is formed whereby an existing tile becomes the parent of new child tiles. Starting with an initial configuration of the tiles scaled to cover the image plane, the tiling is refined through an iterative process of deflation and substitution. The tile vertices or centroids are used to derive a point set from the tiling. The choice of the initial configuration appears mirrored in the global structure and symmetry properties of the tiling and the resulting point set.
Phase function contour map  Phase function level set 
Quasicrystal phase function: defined for regular decagon vertices and golden ratio with dot product and complex conjugate .
Evolution of a 1D quasicrystal point set
from to .
Evolution of a Delaunay graph of a 2D quasicrystal,
where 10 new sites are added at each generation.
Evolution of a Delaunay triangulation of a 2D quasicrystal, where the acceptance window is a decagon centered at the origin that is expanded by a factor of , from radius on the left and to radius on the right.
We focus on a more general class of nonperiodic point sets corresponding to cutandproject quasicrystals [3, 41, 45]. The cutandproject method was originally introduced by Meyer [32] in the context of harmonic analysis and it was later adapted for generating quasicrystals by Moody and Patera [35]. The Fibonacci chain and Penrose tilings can be regarded as special cases of such quasicrystals. In our work, we employ the standard root lattice of the noncrystallographic Coxeter group , a group of reflections taken from Lie algebra theory. This approach to quasicrystals can be used to relate discrete, nonperiodic point sets and tilings (Fig. 1) with the level sets of continuous, nonperiodic functions (Fig. 2). To produce a 2D cutandproject quasicrystal, a 4D periodic lattice is projected on a suitable 2D plane that is irrationally oriented with respect to the lattice. Using this method, we obtain a dense subspace consisting of integer coefficient linear combinations of the vertices of a regular decagon centered on the origin, which are the roots of the noncrystallographic Coxeter group
. This construction ensures that the coordinates of all its elements can be expressed using only integers and an irrational number, the golden ratio. To obtain a finite 2D quasicrystal, we select only those elements of the dense subspace contained in a specified bounded region, called the viewing window, that are mapped by an everywhere discontinuous algebraic transformation, called the star map, to another specified bounded region, called the acceptance window. Through the gradual expansion of a rotationally symmetric acceptance window, a quasicrystal sequence can be uniquely ordered by radial distance and angle from then center of the acceptance window in order to produce a uniformly spacefilling point set in the viewing window. As an important practical consequence, this property directly enables progressive sampling. It could also potentially enable adaptive sampling by varying the radius of the acceptance window according to an application dependent importance map defined for the viewing window. A 2D quasicrystal can also be expressed as a subset of a Cartesian product of two 1D quasicrystals, in accordance with the fact that the points that lie on any straight line through a 2D quasicrystal correspond to some linearly transformed 1D quasicrystal. Therefore, in practice, a 2D quasicrystal (Fig.
4) can be generated from a 2D lattice of 1D quasicrystals. Meanwhile, a 1D quasicrystal (Fig. 3) is produced by taking a strip of a 2D periodic lattice, having finite width and irrational slope, and orthogonally projecting its points onto a line of the same slope. The resulting 1D quasicrystal is composed of at most three distinct tiles. It is easy to generate 1D quasicrystal points using an iterative numerical algorithm. Alternatively, it is possible to exploit the selfsimilar structure of a 1D quasicrystal, viewing it as the fixed point of a set of substitution rules that act recursively on a finite alphabet of possible tile arrangements.Based on the geometry and algebra of the golden ratio, these quasicrystal point sets exhibit some remarkable properties. They display pentagonal and decagonal rotational symmetries, which cannot occur in any periodic point set. Originally, the theory of quasicrystals was motivated by solid state physics as a model of the nonperiodic geometric structures that describe the symmetries of a new kind of longrange atomic order discovered in certain metallic alloys [47]. While translational symmetries define periodic crystals, inflation symmetries can be used to describe quasicrystals based on algebraic irrational numbers, such as the golden ratio. When the acceptance window is a convex region, every point of a quasicrystal can be viewed as a center of inflation symmetry. A quasicrystal can have no translational symmetries and no periodic subsets. Moreover, it can be partitioned into subsets such that each subset forms a valid quasicrystal. Consider the local configurations of tiles in an infinite quasicrystal mosaic formed by a Voronoi diagram or a Delaunay triangulation. According to the repetitivity property implied by suitable regularity conditions, each fragment is repeated infinitely many times in the mosaic. Yet no single fragment is ever sufficient to determine the structure of the whole mosaic because every finite fragment, no matter its size, occurs in an uncountable infinity of nonequivalent mosaics. Furthermore, according to the Delaunay point set property, quasicrystals are both uniformly discrete and relatively dense, creating spacefilling tilings.
As Delaunay point sets, quasicrystals are particularly well suited to image sampling. They enforce both a minimal and a maximal distance between each sample site and its closest neighboring site. In quasicrystal sampling, we rely on the golden ratio to ensure symmetry and self similarity, which are generally absent when other than algebraic irrational numbers are used with the cutandproject method to produce nonperiodic point sets. By taking this approach, we can endow a rendition with a decorative symmetry, which viewers may find attractive in the context of nonphotorealistic rendering. Compared with the regular grids of periodic point sets, the selfsimilar, spacefilling structure of nonperiodic point sets (Fig. 5) appears less monotonous, especially during progressive image rendering. In effect, the geometric structure of quasicrystal sampling eliminates the possibility of aliasing artifacts regularly repeating in the rendition. However, in quasicrystal sampling, fixed local configurations of sample sites can be repeated at multiple places and orientations, albeit not at regular intervals, with the potential to yield some recurring, anisotropic aliasing artifacts. Although we did do so in this work, for photorealistic image reconstruction, it is preferable to avoid inducing global rotational symmetry in the sampling pattern, which is done by ensuring the viewing window does not contain the origin when the acceptance window is symmetric with respect to the origin.
3 Evaluation
We now present a qualitative evaluation of quasicrystal sampling in the context of nonadaptive sampling strategies for use in image representation. In this application, nonadaptive sampling strategies serve as building blocks for interactive sampling, adaptive sampling, and importance sampling techniques. In effect, they dictate the placement of sample sites in image regions sampled at a uniform resolution. For this purpose, a nonadaptive sampling strategy should satisfy several image representation objectives. The sample sites should be distributed in a manner that fairly and accurately represents the image. The sampling pattern should be evenly spacefilling in order to enable progressive image rendering. Without any preconceptions about the distribution of visually salient features in the image, the same amount of information should be devoted to capturing each part of the image. Hence, the number of sample sites placed in any region of the image should be proportional to its area, so the sampling density remains the same throughout the image. Both globally and locally, the placement of sample sites should be uniform and isotropic while still allowing for a variety of different sample site configurations. To prevent aliasing, the sample sites should not be arranged into fixed configurations that visibly repeat locally or globally. To prevent clustering, a minimum distance between sample sites should be maintained throughout the image. Assuming that the correlation between pixels decreases with distance, a sample site should be placed close to the centroid of its Voronoi polygon to enable its sampled color to be most representative of its region of influence when the image is reconstructed using a local interpolation technique. Naturally, the relative importance of these considerations depends on the requirements of a particular application. We have formulated these priorities based on our previous work on a multiresolution image representation [17] designed to simultaneously support both photorealistic image reconstruction and nonphotorealistic image rendering. While previous qualitative surveys [11] and quantitative comparisons [48] of image sampling have focused on applications in photorealistic image reconstruction, they did not cover quasicrystal sampling and farthest point sampling. Furthermore, their tests were not carried out on an image representation of digital photographs and they did not specifically address the needs of nonphotorealistic image rendering.
We compared quasicrystal sampling to a number of standard nonadaptive image sampling strategies [11], reflecting different approaches to the inherent tradeoff between aliasing and noise. For our evaluation, we chose approaches that exemplify divergent aims in sampling. For each approach, we selected a representative implementation. As noted below, alternative implementations are certainly possible but they are likely to produce similar qualitative results. From the deterministic to the stochastic, we tested a range of sampling strategies (Fig. 7):

Periodic sampling [37] aims for global regularity. Our implementation relies on a square lattice refined in scan line order. An alternative implementation could use a hexagonal lattice, the densest periodic lattice in the plane.

Quasicrystal sampling [45] aims for local regularity. Our implementation relies on the cutandproject method applied using the golden ratio. An alternative implementation could use a Penrose tiling produced using a hierarchical substitution algorithm.

Farthest point sampling [10] aims for spatial uniformity. Our implementation relies on the principle of progressively sampling at the point of least information, placing each new sample site at the point farthest from any preceding sample site, which is necessarily a vertex of the Voronoi diagram of the preceding sample sites. An alternative implementation could position sample sites to conform to a centroidal Voronoi diagram so that each sample site is placed at the centroid of its Voronoi polygon.

Jittered sampling [23] aims for local variability. Our implementation relies on a full random displacement of a square lattice refined in scan line order. An alternative implementation could use a partial random displacement of a hexagonal lattice.

Quasirandom sampling [42] aims for low discrepancy. Our implementation relies on the Halton sequence. An alternative implementation could use a Sobol sequence.

Random sampling [11] aims for global variability. Our implementation relies on a uniform distribution. An alternative implementation could use a random walk on a unit square with toroidal boundary conditions.
To perform a qualitative evaluation of the image sampling strategies, we applied a number of computer visualization techniques. For each nonadaptive sampling strategy, we visualize its sample sites (Fig. 7) in the spatial domain using a Voronoi diagram (Fig. 8) and in the frequency domain using a Fourier power spectrum (Fig. 9). We examined the visual effects of applying the image sampling strategies in the context of various image rendering techniques that are used in multiresolution image representations. For photorealistic image reconstruction [2], we tested the accuracy of Shepard’s interpolation (Fig. 13), an inverse distance weighted interpolation method that applies the Voronoi diagram to determine the color of each pixel based on its four nearest sample sites, as well as Gouraud shading (Fig. 14), a piecewise linear interpolation method that applies the Delaunay triangulation to determine the color of each pixel based on its three surrounding sample sites. To quantitatively assess the results of these widely used local interpolation techniques, we relied on the peak signaltonoise ratio (PSNR). This standard image fidelity metric [46]estimates the accuracy of the rendition according to the negative logarithm of the mean squared error between the rendered and actual RGB color values of the pixels. Hence, higher peak signaltonoise ratio scores are considered better. For nonphotorealistic image rendering [17], we experimented with a simple “paint strokes” rendering style (Fig. 12), which applies geometric subdivision to the Delaunay triangulation of the sample sites and then uses linear and nonlinear interpolation to emphasize the transitions between the sampled colors.
Our mosaic rendering style (Fig. 6) offers a new computer visualization tool for evaluating the spatial properties of point sets, such as the sample sites produced by the various image sampling strategies. To produce a mosaic rendering, we apply geometric subdivision to the Delaunay triangulation of the sample sites. The midpoints of the edges of each Delaunay triangle are joined to form three outer triangles and one inner triangle. Each outer triangle is rendered with the color sampled at its vertex of the original Delaunay triangle, while the central triangle is colored black. As each sample site is represented by a starshaped polygonal tile, the resulting mosaic appears packed as tightly as possible, with the black central triangles serving as grout between the tiles. As a sample site’s local neighborhood (Fig. 6, far left) comprises the surrounding sample sites connected to it by edges in the Delaunay triangulation (Fig. 6, center left), the sample site’s mosaic tile (Fig. 6, far right) is shaped to reflect the star of its surrounding Delaunay triangles (Fig. 6, center right). For instance, a sample site’s mosaic tile is a convex or concave polygon according to whether its neighboring sites are arranged in a convex or concave configuration. As a visualization tool, the advantage of mosaic rendering is that the layout of the starshaped mosaic tiles makes the spatial properties of a sampling strategy easier to see at a glance than the triangles of a Delaunay triangulation or the convex polygons of a Voronoi diagram. The sizes of the mosaic tiles are indicative of the uniformity of sampling, as coarsely sampled regions give rise to large tiles and finely sampled regions give rise to small tiles, making common defects such as clustering, undersampling, and oversampling easy to detect. The orientations of the mosaic tiles are indicative of the isotropy of sampling, as the preferred directions of the sampling are revealed in the preferred rotations of the tiles, making global or local grid structures easy to detect. The shapes of the mosaic tiles are indicative of the heterogeneity of sampling, as the local configurations of neighboring sites uniquely determine the tile polygons, making repetitive patterns easy to detect. For instance, farthest point sampling produces tiles of uniform size and similar shape to create the appearance of a pebble mosaic, while quasicrystal sampling yields a decorative tiling with just a small set of possible tile shapes.
Our qualitative evaluation of image sampling strategies uses seven criteria (Fig. 10) known to affect the visual quality of photorealistic image reconstruction and nonphotorealistic image rendering [17]:

Accurate reconstruction requires the rendition to faithfully represent the likeness of the original image. It is a necessary but not sufficient condition of success in both photorealistic and nonphotorealistic image rendering. This objective appears to be closely associated with uniform coverage and centroidal regions. It is assessed by measuring the peak signaltonoise ratio for the results of photorealistic image reconstruction (Figs. 13 and 14). Its most pronounced effects can also be observed in the results of nonphotorealistic image rendering (Fig. 12). When the resolution of sampling is uniform, periodic sampling yields the most accurate image interpolation (in Fig. 14, for the regular square grid, this takes place when there are , , and samples). However, when regions of varying resolution arise during progressive refinement, the accuracy of periodic sampling can substantially deteriorate. Similar behavior is observed in jittered sampling since it applies random perturbations to a periodic point set. Farthest point sampling produces image reconstructions that are nearly as accurate as periodic sampling, but its performance does not diminish during progressive refinement. Intermediate accuracy is offered by quasicrystal sampling, which appears to be slightly more accurate than quasirandom sampling. The least accurate reconstructions are produced by jittered and random sampling. Given its popularity in computer graphics implementations, the poor performance of jittered sampling is rather disappointing. Of course, the accuracy of jittered sampling can always be made closer to that of periodic sampling by reducing the amount of random displacement, which risks reintroducing the aliasing artifacts of periodic sampling. In general, accuracy is improved by uniformity and reduced by randomness, an effect that can be readily seen as producing tight or loose image stylization.
Figure 7: Nonadaptive sampling strategies: periodic sampling (top left), quasicrystal sampling (top center), farthest point sampling (top right), jittered sampling (bottom left), quasirandom sampling (bottom center), and random sampling (bottom right). Sampling starts with the dark blue sites and finishes with the light green sites. 
Progressive refinement requires the sample sites to smoothly fill the available space, avoiding abrupt changes in appearance as new sample sites are sequentially added to the rendition. This objective serves to enable a multiresolution image representation to support progressive rendering of compressed images based on an incremental sampling of the image data. It is assessed by examining the spatial layout of the sequence of sample sites (Fig. 7). Under ideal circumstances, progressive refinement should yield a smooth curve for the peak signaltonoise ratio (Fig. 14). The best progressive refinement results are produced by farthest point sampling and quasicrystal sampling, as these methods maintain a uniform sampling density by ensuring that new sample sites are placed in the largest empty spaces between the existing sample sites. Quasirandom sampling proves slightly less proficient, as it places some sample sites very close together while keeping others far apart. Random sampling is even less effective due to its tendency to locally cluster sample sites. The regular grids used in periodic and jittered sampling are not suitable for smooth progressive refinement, especially when they are refined in scan line order. While other refinement schemes can be applied to regular grids, such as refinement in random order, their intrinsic symmetry makes it difficult to smoothly increase the sampling density throughout the image.

Uniform coverage requires the sample sites to be evenly distributed regardless of position, avoiding configurations that place sample sites too close or too far from their nearest neighbors. This objective is assessed using Voronoi diagrams (Fig. 8) as well as mosaic renderings (Fig. 11). Its effect determines the sizes of brush strokes in nonphotorealistic image rendering (Fig. 12). Uniform coverage is associated with a Fourier power spectrum (Fig. 9) that displays an empty ring around the central spike, as low frequencies are attenuated in favor of a threshold frequency corresponding to the most commonly observed nearest neighbor distance between sample sites. Although a blue noise spectrum can ensure uniform coverage, it is not a necessary condition. When the resolution of sampling is uniform, periodic sampling generates uniform coverage, as its mosaic tiles are all exactly the same size. However, periodic sampling cannot sustain uniform coverage during progressive refinement. By design, farthest point sampling maintains uniform coverage at all times, as its mosaic tiles are all approximately the same size. Quasicrystal sampling maintains nearly as uniform coverage, as its mosaic tiles are limited to just a few comparable sizes. While quasirandom and jittered sampling strive to uphold a uniform density of sampling, they nevertheless are less effective at providing uniform coverage, as their mosaic tiles come in many sizes. In the case of jittered sampling, uniform coverage can be improved by reducing the amount of random displacement. Random sampling does not give uniform coverage, as its mosaic tiles exhibit the greatest range of different sizes.

Isotropic distribution requires the sample sites to be evenly distributed regardless of orientation, avoiding configurations that align sample sites along globally or locally preferred directions. This objective is assessed using Voronoi diagrams (Fig. 8) as well as mosaic renderings (Fig. 11). Its effect is to determine the orientations of brush strokes in nonphotorealistic image rendering (Fig. 12). An isotropic distribution produces a Fourier power spectrum (Fig. 9) that displays a rotational symmetry around the central spike, as the power at each frequency does not depend on its orientation. Although a blue noise spectrum can ensure isotropic distribution, it is not a necessary condition. Random sampling is the most isotropic, as its sample sites are both locally and globally uncorrelated. Farthest point and jittered sampling are nearly as isotropic, as their sample sites can exhibit slight local alignment. Farthest point samples can appear to be placed in roughly hexagonal local configurations. Jittered samples can appear to retain some of the structure of the underlying square grid, as the isotropy of jittered sampling reflects the amount of random displacement used to generate the sampling. Quasirandom sampling has intermediate isotropy, as its sample sites can exhibit slight global alignment, which can be verified in the lack of radial symmetry in its Fourier power spectrum. Periodic sampling and quasicrystal sampling do not have isotropic distribution since their sample sites are globally aligned along predetermined axes.
Figure 9: Fourier power spectrums of image sampling strategies: periodic sampling (top left), quasicrystal sampling (top center), farthest point sampling (top right), jittered sampling (bottom left), quasirandom sampling (bottom center), and random sampling (bottom right). 
Blue noise spectrum requires the sample sites to be distributed similarly to a Poisson disk distribution, a random point field conditioned on a minimum distance between the points. According to this objective, for an image sampling strategy to provide effective antialiasing for image rendering, it should attempt to mimic the idealized distribution of photoreceptors in the human eye. Usually implying both uniform coverage and isotropic distribution, a blue noise spectrum is highly desirable in many computer graphics applications, particularly photorealistic image reconstruction. It is assessed by examining the Fourier power spectrums of the sampling strategies (Fig. 9) for a radially symmetric profile that concentrates noise in the high frequencies while attenuating the power of the low frequencies, thereby eliminating the aliasing artifacts associated with low frequency patterns that can appear distracting to the eye. In effect, a blue noise spectrum exhibits a disk of low power around the origin, surrounded by roughly constant power at the higher frequencies. Its effects can be judged according to the amount of aliasing present in photorealistic image reconstruction (Fig. 13) and nonphotorealistic image rendering (Fig. 12
). Farthest point sampling has a Fourier power spectrum that is closest to a blue noise spectrum. Jittered sampling attempts to replicate the blue noise spectrum, but it does not clearly exhibit the threshold frequency ripple around the central spike. Quasirandom is even less successful because its Fourier power spectrum lacks radial symmetry. Periodic, quasicrystal, and random sampling have Fourier power spectrums that do not resemble the blue noise spectrum. The Fourier power spectrums of periodic and quasicrystal sampling reflect the spatial structures and directional symmetries of the lattices used to place the sample sites. On the other hand, the white noise spectrum of random sampling assigns roughly the same power to all frequencies.

Centroidal regions require sample sites to be well centered with respect to their Voronoi polygons, approximating a centroidal Voronoi diagram. Typically associated with uniform coverage and accurate reconstruction, this objective is popular in nonphotorealistic image rendering. Sampling strategies, such as periodic sampling, that produce centroidal regions can still be prone to aliasing artifacts since centroidal regions do not guarantee a blue noise spectrum. Centroidal regions can be readily assessed using the Voronoi diagrams (Fig. 8). The effects can also be observed in the shapes of tiles in mosaic rendering (Fig. 11) and brush strokes in nonphotorealistic image rendering (Fig. 12). Periodic sampling, placing each sample site at the same distance from all of its nearest neighbors, generates exact centroidal Voronoi regions. Quasicrystal and farthest point sampling produce approximately centroidal Voronoi regions. To place sample sites close to the center of their Voronoi polygons, quasicrystal sampling relies on local symmetries while farthest point sampling relies on nearest neighbor distance. Jittered sampling and quasirandom sampling have difficulty ensuring centroidal regions because they are less effective at maintaining a minimal nearest neighbor distance. Finally, random sampling can only produce centroidal regions by chance.

Heterogeneous configurations require sample sites to be placed in a variety of different local arrangements, avoiding regularly or randomly repeating the same sampling patterns. While this objective is not traditionally a concern in photorealistic image reconstruction, it helps to prevent nonphotorealistic image rendering from appearing too perfect, seemingly mechanical and artificial. For instance, it helps to give a vibrant appearance to brush stroke rendering. Typically, when sampling strategies yield centroidal regions, they also tend to produce homogeneous configurations and vice versa, illustrating an apparent tradeoff between these competing objectives. Heterogeneous configurations are assessed using the Voronoi diagrams (Fig. 8) and mosaic renderings (Fig. 11). Their effect is also visible in the arrangement of brush strokes in nonphotorealistic image rendering (Fig. 12). Random sampling produces the most heterogeneous local configurations, as its mosaic tiles exhibit the greatest variety of shapes. Jittered and quasirandom sampling are nearly as heterogeneous, since their mosaic tiles are almost as widely varied, though few of them are exceptionally large in size. Quasicrystal and farthest point sampling are far less heterogeneous. By upholding local symmetries, quasicrystal sampling causes sample sites to have only a few possible local configurations, resulting in mosaic tiles that have only a few possible shapes. By upholding nearest neighbor distance, farthest point sampling causes sample sites to have similar local configurations, resulting in mosaic tiles that look very much alike, mostly convex and rounded. Based on repetitions of a single local configuration, periodic sampling is entirely homogeneous.
Sampling  Accurate  Progressive  Uniform  Isotropic  Blue Noise  Centroidal  Heteroeneous 

Strategies  Reconstruction  Refinement  Coverage  Distribution  Spectrum  Regions  Configurations 
Periodic  
Quasicrystal  
Farthest Point  
Jittered  
Quasirandom  
Random 
Superior Good Fair Poor
Our qualitative analysis (Fig. 10) indicates that quasicrystal sampling offers a useful compromise between the ordered behavior of standard periodic sampling using a regular square lattice and the disordered behavior of standard Monte Carlo sampling using jittered, quasirandom, or random sampling. Compared with periodic sampling, quasicrystal sampling displays a greater variety of local sample site configurations resulting in smoother progressive refinement, although its sampling patterns are somewhat less uniform, leading to lower accuracy of image reconstruction. Compared with jittered, quasirandom, and random sampling, quasicrystal sampling displays more uniform coverage resulting in better accuracy of image reconstruction, although its sampling patterns are anisotropic, exhibiting significantly less variety of local sample site configurations. By virtue of its deterministic construction, quasicrystal sampling does not suffer from the variability that can affect the results of random sampling, jittered sampling, and, to a much lesser degree, farthest point sampling. Nevertheless, its lack of a blue noise power spectrum renders it rather susceptible to aliasing artifacts. Research on quasicrystal sampling based on the Penrose tiling [38] suggests that it may be possible to partially alleviate this problem by taking advantage of the symmetries and the repetitions of the local sample site configurations in order to systematically displace the sample sites in a manner that improves the spectral properties of the sampling pattern.
Based on our qualitative evaluation of the various nonadaptive sampling strategies (Fig. 10), we recommend a blue noise sampling strategy, such as farthest point sampling, for general use in image representation. In particular, farthest point sampling does not perform poorly on any of our seven evaluation criteria. Overall, our qualitative evaluation of nonadaptive sampling strategies is in broad agreement with previous studies, which did not consider quasicrystal sampling. They emphasized the importance of Poisson disk distributions [11] and low discrepancy distributions [48], which are exemplified in our evaluation by farthest point sampling and quasirandom sampling respectively. Hence, the good overall performance of these two techniques should come as no surprise. Farthest point sampling performs better than quasirandom sampling on six out of the seven evaluation criteria. For the majority of our evaluation criteria, quasicrystal sampling performs no better than farthest point sampling and no worse than quasirandom sampling. Nevertheless, from a practical point of view, quasicrystal sampling is significantly simpler to implement and calculate than farthest point sampling, which relies on maintaining complex geometric data structures to keep track of the vertices of a Voronoi diagram. This could be an important consideration for imaging applications on mobile devices that have limited processing and storage capabilities. From a theoretical point of view, the deterministic algebraic construction of quasicrystals renders their sampling patterns particularly well suited to mathematical analysis. Presenting possibilities for future research, the cutandproject method could be adapted for higher dimensional sampling or adaptive sampling applications.
In future work, it would also be interesting to explore the relationship between local symmetry and sampling quality. The cutandproject method can be used to generate nonperiodic point sets with different symmetries, not just the pentagonal and decagonal symmetries associated with the golden ratio, as shown in this work. Just as for periodic sampling it would be interesting to compare the image reconstruction accuracy of square and hexagonal grids, for nonperiodic sampling it would be interesting to compare our decagonal quasicrystal tiling with the dodecagonal Socolar tiling, which was recently proposed for use in sampling applications [39].
4 Conclusion
Cutandproject quasicrystals present new possibilities for image sampling in computer graphics. This nonperiodic sampling approach deterministically generates uniformly spacefilling point sets, ensuring that sample sites are evenly distributed throughout the image. It offers a useful compromise between predictability and randomness, between the standard periodic sampling and the standard Monte Carlo sampling methods. Although blue noise sampling can generate higher quality sampling patterns for photorealistic image reconstruction, quasicrystal sampling may prove much simpler to implement and calculate. In the context of nonphotorealistic image rendering, quasicrystal sampling may prove attractive for its symmetry properties.
Periodic sampling  Quasicrystal sampling 
Farthest point sampling  Jittered sampling 
Quasirandom sampling  Random sampling 
Periodic sampling  Quasicrystal sampling 
Farthest point sampling  Jittered sampling 
Quasirandom sampling  Random sampling 
Periodic sampling (PSNR )  Quasicrystal sampling (PSNR ) 
Farthest point sampling (PSNR )  Jittered sampling (PSNR ) 
Quasirandom sampling (PSNR )  Random sampling (PSNR ) 
Acknowledgements
We are grateful to the referees for valuable remarks that led to improvements of the original manuscript. We also acknowledge the financial support of the Natural Sciences and Engineering Research Council of Canada, le Fonds Québécois de la Recherche sur la Nature et les Technologies, as well as the grants MSM6840770039 and LC06002 of the Ministry of Education of the Czech Republic, and GA201/09/0584 of the Czech Science Foundation. We are also grateful for the support of the MIND Research Institute and the Merck Frosst Company. Mark Grundland further acknowledges the financial assistance of the Celanese Canada Internationalist Fellowship, the British Council, the Cambridge Commonwealth Trust, and the Overseas Research Student Award Scheme. The images were provided by FreeFoto.com and the Waterloo Brag Zone.
Supplementary materials
For a full resolution version of this paper, along with supplementary materials, please visit: http://www.Eyemaginary.com/Portfolio/Publications.html.
[1]Referencesref
References
 [1]
 [2] Amidror I., Scattered data interpolation methods for electronic imaging systems: a survey, J. Electronic Imaging 11 (2002), 157–176.
 [3] Chen L., Moody R.V., Patera J., Noncrystallographic root systems, in Quasicrystals and Discrete Geometry (Toronto, ON, 1995), Fields Inst. Monogr., Vol. 10, Amer. Math. Soc., Providence, RI, 1998, 135–178.
 [4] Cohen M.F., Shade J., Hiller S., Deussen O., Wang tiles for image and texture generation, in Proceedings of SIGGRAPH, 2003, 287–294.
 [5] Cook R.L., Stochastic sampling in computer graphics, ACM Trans. Graphics 5 (1986), 51–72.
 [6] Deussen O., Hiller S., van Overveld C., Strothotte T., Floating points: a method for computing stipple drawings, in Proceedings of EUROGRAPHICS, 2000, 41–50.
 [7] Dippe M.A.Z., Wold E.H., Antialiasing through stochastic sampling, in Proceedings of SIGGRAPH, 1985, 69–78.
 [8] Du Q., Faber V., Gunzburger M., Centroidal Voronoi tessellations: applications and algorithms, SIAM Rev. 41 (1999), 637–676.
 [9] Dunbar D., Humphreys G., A spatial data structure for fast Poissondisk sample generation, in Proceedings of SIGGRAPH, 2006, 503–508.
 [10] Eldar Y., Lindenbaum M., Porat M., Zeevi Y.Y., The farthest point strategy for progressive image sampling, IEEE Trans. Image Process. 6 (1997), 1305–1315.
 [11] Glassner A., Principles of digital image synthesis, Vol. 1, Morgan Kaufmann, 1995.
 [12] Glassner A., Aperiodic tiling, IEEE Computer Graphics and Applications 18 (1998), no. 3, 83–90.
 [13] Glassner A., Penrose tiling, IEEE Computer Graphics and Applications 18 (1998), no. 4, 78–86.
 [14] GoodmanStrauss C., Aperiodic hierarchical tilings, in Foams, Emulsions, and Cellular Materials (Cargèse, 1997), NATO Adv. Sci. Inst. Ser. E Appl. Sci., Vol. 354, Kluwer Acad. Publ., Dordrecht, 1999, 481–496.
 [15] Grunbaum B., Shephard G.C., Tilings and patterns, WH Freeman, 1987.
 [16] Grundland M., Style and content in digital imaging: Reconciling aesthetics with efficiency in image representation, VDM, 2008.
 [17] Grundland M., Gibbs C., Dodgson N.A., Stylized multiresolution image representation, J. Electronic Imaging 17 (2008), 013009, 1–17.
 [18] Hausner A., Simulating decorative mosaics, in Proceedings of SIGGRAPH, 2001, 573–580.
 [19] Hausner A., Pointillist halftoning, in Proceedings of the International Conference on Computer Graphics and Imaging, 2005, 134–139.
 [20] Hiller S., Deussen O., Keller A., Tiled blue noise samples, in Proceedings of Vision, Modeling and Visualization, 2001, 265–271.
 [21] Hiller S., Hellwig H., Deussen O., Beyond stippling – methods for distributing objects on the plane, in Proceedings of EUROGRAPHICS, 2003, 515–522.
 [22] Jones T.R., Efficient generation of Poissondisk sampling patterns, J. Graphics Tools 11 (2006), no. 2, 27–36.
 [23] Klassen R.V., Filtered jitter, Computer Graphics Forum 19 (2000), no. 4, 223–230.
 [24] Kopf J., CohenOr D., Deussen O., Lischinski D., Recursive Wang tiles for realtime blue noise, in Proceedings of SIGGRAPH, 2006, 509–518.
 [25] Lagae A., Dutre P., A procedural object distribution function, ACM Trans. Graphics 24 (2005), 1442–1461.
 [26] Lagae A., Dutre P., An alternative for Wang tiles: colored edges versus colored corners, ACM Trans. Graphics 25 (2006), 1442–1459.
 [27] Lagae A., Dutre P., A comparison of methods for generating Poisson disk distributions, Computer Graphics Forum 27 (2008), no. 1, 114–129.
 [28] Lagae A., Kaplan C.S., Fu C.W., Ostromoukhov V., Deussen O., Tilebased methods for interactive applications, SIGGRAPH 2008 Class Notes, ACM, 2008.
 [29] Lu P.J., Steinhardt P.J., Decagonal and quasicrystalline tilings in medieval Islamic architecture, Science 315 (2007), no. 5815, 1106–1110.
 [30] Masáková Z., Patera J., Zich J., Classification of Voronoi and Delone tiles of quasicrystals. III. Decagonal acceptance window of any size, J. Phys. A: Math. Gen. 38 (2005), 1947–1960.
 [31] McCool M., Fiume E., Hierarchical Poisson disk sampling distributions, in Proceedings of Graphics Interface, 1992, 94–105.
 [32] Meyer Y., Algebraic numbers and harmonic analysis, NorthHolland, 1972.
 [33] Mitchell D.P., Spectrally optimal sampling for distribution ray tracing, in Proceedings of SIGGRAPH, 1991, 157–164.
 [34] Mojsilovic A., Soljanin E., Color quantization and processing by Fibonacci lattices, IEEE Trans. Image Process. 10 (2001), 1712–1725.
 [35] Moody R.V., Patera J., Quasicrystals and icosians, J. Phys. A: Math. Gen. 26 (1993), 2829–2853.
 [36] Moody R.V., Patera J., Dynamical generation of quasicrystals, Lett. Math. Phys. 36 (1996), 291–300.
 [37] Ostromoukhov V., Mathematical tools for computergenerated ornamental patterns, Lecture Notes in Computer Science, Vol. 1375, 1998, 193–223.
 [38] Ostromoukhov V., Donohue C., Jodoin P.M., Fast hierarchical importance sampling with blue noise properties, in Proceedings of SIGGRAPH, 2004, 488–495.
 [39] Ostromoukhov V., Building 2D lowdiscrepancy sequences for hierarchical importance sampling using dodecagonal aperiodic tiling, in Proceedings of GRAPHICON, 2007, 139–142.
 [40] Ostromoukhov V., Sampling with polyominoes, in Proceedings of SIGGRAPH, 2007, 078, 1–6.
 [41] Patera J., Noncrystallographic root systems and quasicrystals. in The Mathematics of LongRange Aperiodic Order (Waterloo, ON, 1995), NATO Adv. Sci. Inst. Ser. C Math. Phys. Sci., Vol. 489, Kluwer Acad. Publ., Dordrecht, 1997, 443–465.
 [42] Press W., Teukolsky S.A., Vetterling W.T., Flannery B.P., Numerical recipes in C, 2nd ed., Cambridge University Press, 1992.
 [43] RangelMondragon J., Abas S.J., Computer generation of Penrose tilings, Computer Graphics Forum 7 (1988), no. 1, 29–37.
 [44] Secord A., Weighted Voronoi stippling, in Proceedings of the Second International Symposium on NonPhotorealistic Animation and Rendering, 2002, 37–43.
 [45] Senechal M., Quasicrystals and geometry, Cambridge University Press, Cambridge, 1995.
 [46] Sharma G., Digital color imaging handbook, CRC Press, 2003.
 [47] Shechtman D., Blech I., Gratias D., Cahn J.W., Metallic phase with longrange orientational order and no translational symmetry, Phys. Rev. Lett. 53 (1984), 1951–1953.
 [48] Shirley P., Discrepancy as a quality measure for sample distributions, in Proceedings of EUROGRAPHICS, 1991, 183–194.
 [49] Stam J., Aperiodic texture mapping, European Research Consortium for Informatics and Mathematics, Technical Report ERCIM01/97R046, 1997.
 [50] Wei L.Y., Tilebased texture mapping on graphics hardware, in Proceedings of the ACM Conference on Graphics Hardware, 2004, 55–63.
 [51] Wei L.Y., Parallel Poisson disk sampling, in Proceedings of SIGGRAPH, 2008, 020, 1–10.
 [52] White K.B., Cline D., Egbert P.K., Poisson disk point sets by hierarchical dart throwing, in Proceedings of the IEEE Symposium on Interactive Ray Tracing, 2007, 129–132.
Comments
There are no comments yet.