Computation of Compact Distributions of Discrete Elements

: In our daily lives, many plane patterns can actually be regarded as a compact distribution of a number of elements with certain shapes, like the classic pattern mosaic. In order to synthesize this kind of pattern, the basic problem is, with given graphics elements with certain shapes, to distribute a large number of these elements within a plane region in a possibly random and compact way. It is not easy to achieve this because it not only involves complicated adjacency calculations, but also is closely related to the shape of the elements. This paper attempts to propose an approach that can effectively and quickly synthesize compact distributions of elements of a variety of shapes. The primary idea is that with the seed points and distribution region given as premise, the generation of the Centroidal Voronoi Tesselation (CVT) of this region by iterative relaxation and the CVT will partition the distribution area into small regions of Voronoi, with each region representing the space of an element, to achieve a compact distribution of all the elements. In the generation process of Voronoi diagram, we adopt various distance metrics to control the shape of the generated Voronoi regions, and ﬁnally achieve the compact element distributions of different shapes. Additionally, approaches are introduced to control the sizes and directions of the Voronoi regions to generate element distributions with size and direction variations during the Voronoi diagram generation process to enrich the effect of compact element distributions. Moreover, to increase the synthesis efﬁciency, the time-consuming Voronoi diagram generation process was converted into a graphical rendering process, thus increasing the speed of the synthesis process. This paper is an exploration of elements compact distribution and also carries application value in the ﬁelds like mosaic pattern synthesis.


Introduction
The synthesis of texture patterns has been an important research subject.Traditional texture synthesis technologies are mainly designed for raster textures, which record texture information only in the form of raster images and have no explicit records or expressions of the geometric content expressed in the textures (such as basic units composing the textures and unit distribution information).Gradually, a number of texture patterns have been found to be distributions of numerous discrete elements.If images are described as geometric distributions of many vector elements, a better effect of pattern synthesis will be achieved from the aspect of element distribution.On this basis, many researchers have begun studying the topics of discrete element distribution [1][2][3][4][5][6][7][8][9], among which the synthesis of compact element distributions is the most challenging.Compact element distributions need to arrange one or several given types of elements in a random and possibly compact way (as show in Figure 1, the compact distribution effects of square and circular elements).This process calls for not only the attention to the element location distributions, but also judgment of the crashing and overlapping relations between elements, which are closely connected with the shapes of the elements.
(1) Compact distributions of elements of various shapes were achieved.By adopting different distance metrics of the Voronoi diagram, the shapes of the generated Voronoi regions can be controlled, thus realizing compact distributions of elements of various shapes, such as rectangles, circles, triangles, and ovals.
(2) Our method is highly efficient.By converting the Voronoi generation process into a graphical rendering process, the GPU (Graphics Processing Unit) parallel computing ability can be fully utilized for calculation, and thus the synthesis speed can reach interactive efficiency.
(3) The pattern synthesis can be controlled.By adjusting parameters, users are able to control the size and distribution direction of elements and further control the pattern synthesis to satisfy various synthesis requirements.
This paper not only conducts discussion on geometric topic of compact element distributions, but also achieves the generation of synthesis results with special artistic beauty, which carries practical significance in fields such as pattern design and artistic creation.

Related Work
In the field of computer graphics, a number of researchers have studied compact distributions, and a series of approaches have been proposed.In the following part of this paper, these approaches will be introduced in three categories, namely, physics-based simulation, sample-based synthesis,

Related Work
In the field of computer graphics, a number of researchers have studied compact distributions, and a series of approaches have been proposed.In the following part of this paper, these approaches will be introduced in three categories, namely, physics-based simulation, sample-based synthesis, and process-based modeling.Because the work in this paper is closely related with the Voronoi diagram, a brief introduction of work on the Voronoi diagram is presented in the last part of this chapter.
Physics-based simulation.During physics-based simulation, the compact distribution of the discrete elements will be regarded as a piling process of physical elements.The discrete elements will be arranged tightly in space according to the laws of dynamics, such as collision, contact, and friction, and some manually defined rules.
Cho et al. [15] utilized a physics-based simulation method and managed to accumulate a certain number of randomly initialized elements in a rectangular box by gravity and friction.Researchers such as Guendelman [16] adopted a dual representation to express nonconvex-shaped elements to effectively deal with the collision and contact during element aggregation, and succeeded in generating realistic dynamic distributions of discrete elements.
Using a physics-based simulation method, users can visually control the modeling results.However, they need to be fully aware of the modeling objects and algorithm application and must also be able to repeatedly test the input parameters based on experience.Therefore, this approach is labor-intensive and time-consuming.Moreover, this method calls for a large amount of data calculation, consequently leading to low operational efficiency.
Sample-based synthesis.During sample-based synthesis, given a small region of an element distribution as a sample, users synthesize a larger-scale element distribution that has a similar pattern to the sample.This method develops from sample-based synthesis, mainly on raster textures in the early stage.Starting from synthesizing by pixel [8,9], later synthesizing by patch [17,18] is introduced with better efficiency and performance, followed by many methods on textures synthesis like image analogies [19].Although these methods could synthesize patterns with large-scale elements distribution based on sample image, they merely deal with raster information with no explicit record or expression of the geometric information in the textures (information like basic element of the textures or distributions of elements), thus often resulting in poor performance when synthesizing patterns with compact element distributions.
Afterwards, research began on the synthesis of distribution relationship of geometric elements in texture patterns.The key point is to construct local position relationship between elements and overall distribution rules of these elements based on the samples.Some researchers referred to local neighborhood matching in texture synthesis for the synthesis process [1,2,10,14,20,21] and attempted to guarantee the local distribution characteristics of the synthesized elements.Barla et al. [20] used Delaunay triangulation to construct the neighborhood relationship between elements and then replaced the seed points one by one based on this relationship.Ijiri et al. [10] used the same Delaunay method, but they adopted a local growing method for the distribution synthesis.Xu et al. [2] proposed a new approach to constructing the neighborhood relationship using a histogram to count the element distributions of a certain element type within a circle neighborhood to describe the neighborhood features of this element.Other researchers performed the synthesis based on the statistical features of the entire element distributions and sought to synthesize data with features similar to those of the samples.Barla et al. [3] calculated quantities such as the distance, angle, separation, and superposition between all the elements in the sample pattern and calculated their averages, variances, and Gaussian distributions.During synthesis, distributions of the peak points will gradually accord with the statistical features of the sample data via constant iteration.Other researchers, such as Hurtut [4], calculated the probability density function of elements of all types and synthesized the element distribution based on the calculation results.Some other researchers considered the statistical and local features of the entire set.Ma [1] extended the local neighborhood matching to an overall energy function optimization process, while Liu et al. [5], in addition to proposing a neighborhood coding method to describe the element neighborhood features, considered the overall features, such as the density and minimum distance, of all kinds of elements.
Sample-based methods have proven to be efficient in synthesizing discrete elements.However, when elements are closely arranged, a number of problems will occur due to the neglect of collision and overlapping situations between elements.If we take the collision and overlapping situations into consideration, the amount of calculation will be significantly increased.
Process-based modeling.In the process-based methods for element distribution simulation, finding the element distributions is regarded as a random sampling or an overall optimization process.
Many methods have been adopted to generate different element space distributions, among which the Poisson disk distribution is the most popular.This distribution can ensure that the elements are as evenly arranged as possible, while the distance between any two elements is not less than a user-defined threshold value.The two most classic algorithms to generate Poisson disk distributions are the dart-throwing algorithm [6] and Lloyd's method [14].Lloyd's method, via progressive optimization in a certain region, constructs the Voronoi diagram to create point distributions.Then, the entire region is partitioned into small regions by every Voronoi region.On this basis, if these polygons have similar shapes to those of the discrete elements, the elements can be closely arranged by simply putting them into these polygons.
Voronoi diagram.The Voronoi diagram is a kind of fundamental data structure of space proximity relations [22].On the basis of point set P of a plane, this diagram can divide this plane into Voronoi regions in the shape of a convex polygon.Every point p i in the P corresponds a region V i , making the distance between any point in V i and p i shorter than that between any point in V i and p j (j =i).When generating the Voronoi diagram, various Voronoi region shapes can be achieved by adopting different distance metrics.Normally, the L p -metrics is the most commonly used (also known as Minkowski distance) [23].
When p = 2 (that is, L 2 -metrics), the shapes of Voronoi regions generated will approach a regular hexagon.When p = 1 (also known as Manhattan distance), the shapes will be the square rotated 45 degrees.Hausner et al. [14] attempted to adopt this function for compact distribution of square-shaped elements.
Apart from L p -metrics, other functions like convex polygon-offset distance function [24] and the Moscow metric [25] are also introduced by researchers.There are also related works on the anisotropic Voronoi synthesis [26].Moreover, the definition of the Voronoi diagram has been extended from two dimensions to three or more dimensions.The Voronoi diagram over curved surfaces [27,28] has also been researched.Some researchers even conducted work on the Voronoi diagram in non-Euclidean space metrics (like Riemannian geometry metric) [29,30].
This paper is mainly on how to achieve the compact distributions of elements in various shapes in plane space.L p -metrics, as one of the most commonly used functions in Euclidean space and the most efficient way to generate the Voronoi diagram, has naturally become our key point.Besides L 1 -metric and L 2 -metric, we also try to achieve distributions of square-shaped elements by L ∞ -metric, and distributions of elements in rectangular and oval shape by axial scale.Additionally, the triangle equidistant line method (regarded as a kind of convex polygon-offset distance function) is also adopted to generate compact distributions of triangle-shaped elements.

Voronoi Diagram Based on Several Distance Metrics
In this paper, the calculation process of the centroidal Voronoi diagram is adopted for the Voronoi diagram generation with the following calculation process.First, generate the Voronoi diagram of N non-overlapping seed points in a plane, which partitions this plane into N regions with all the points in each region closer to the corresponding seed points than to other seed points, as shown in Figure 2a, Then, move the seed point of each region to the centroid of the corresponding region and generate the Voronoi diagram of the seed points again, and repeat the above process.According to the centroidal Voronoi diagram, the distribution of seed points will become increasingly balanced.Normally, the Voronoi region of a set of evenly distributed points will be similar to a relatively regular shape (as shown in Figure 2b).In the most classical Voronoi diagram, Lp-metric is used to measure the distance between the p0( ,  ) and p1( ,  ) and among all the Lp-metrics, L2-metric is the most commonly used.

d(p0, p1) = |𝑥 − 𝑥 | + |𝑦 − 𝑦 |
The polygons of the generated Voronoi diagram were close to square.To obtain Voronoi regions with more shapes, we attempt to adopt more metrics.Four additional distance metrics are introduced, which are Chebyshev distance, Chebyshev distance with axial scale, L2-metric distance with axial scale, and equidistant line distance, and finally, distributions of squares, rectangles, ovals, and triangles are achieved, respectively.In the following part of this paper, the different metrics will be discussed in detail.In the most classical Voronoi diagram, L p -metric is used to measure the distance between the p 0 (x 0 , y 0 ) and p 1 (x 1 , y 1 ) and among all the L p -metrics, L 2 -metric is the most commonly used.
With this metric, the generated Voronoi region is close to a regular hexagon (as shown in Figure 3a).In the most classical Voronoi diagram, Lp-metric is used to measure the distance between the p0( ,  ) and p1( ,  ) and among all the Lp-metrics, L2-metric is the most commonly used.

d(p0, p1) = |𝑥 − 𝑥 | + |𝑦 − 𝑦 |
The polygons of the generated Voronoi diagram were close to square.To obtain Voronoi regions with more shapes, we attempt to adopt more metrics.Four additional distance metrics are introduced, which are Chebyshev distance, Chebyshev distance with axial scale, L2-metric distance with axial scale, and equidistant line distance, and finally, distributions of squares, rectangles, ovals, and triangles are achieved, respectively.In the following part of this paper, the different metrics will be discussed in detail.The polygons of the generated Voronoi diagram were close to square.To obtain Voronoi regions with more shapes, we attempt to adopt more metrics.Four additional distance metrics are introduced, which are Chebyshev distance, Chebyshev distance with axial scale, L 2 -metric distance with axial scale, and equidistant line distance, and finally, distributions of squares, rectangles, ovals, and triangles are achieved, respectively.In the following part of this paper, the different metrics will be discussed in detail.

Chebyshev Distance
The Chebyshev distance defines the distance between two points (x 1 , y 1 ) and (x 0 , y 0 ) as max(|x 1 − x 0 |, |y 1 − y 0 |).It is not difficult to prove that it is equivalent to the L p -metric with p→∞.The Voronoi region generated with this distance is the same shape, namely square, as the one generated with Manhattan metrics (as shown in Figure 3c).However, the differences between these two metrics are the directions of the coordinate axes.The axes of the Manhattan distance are inclined 45 degrees, while the axes of the Chebyshev distance are horizontal and vertical (as shown in Figure 4).In this paper, we use the Chebyshev distance because its axial directions conform with normal image space, which is easier to scale and control.
As shown in Figure 4, the side of the square is defined as L a , while the diagonal line is defined as L b ; the axis xy corresponding to the Manhattan distance is shown in Figure 4a.The distance between any point A(x 1 , y 1 ) on the square side and the origin O(x 0 , y 0 ) is |x 1 − x 0 | + |y 1 − y 0 |, which is the sum of the red segment and blue segment in Figure 4a, that is, L b /2.The axis xy corresponding to the Chebyshev distance is shown in Figure 4b.The distance between any point A(x 1 , y 1 ) on the square side and the origin O(x 0 , y 0 ) is max(|x 1 − x 0 |, |y 1 − y 0 |), which is the length of the longer segment among the red and blue segments in the picture, La /2. with axial scale; and (f) equidistant line distance.The partitioned regions are presented in different colors.

Chebyshev Distance
The Chebyshev distance defines the distance between two points ( ,  ) and ( ,  ) as max(| −  |, | −  |).It is not difficult to prove that it is equivalent to the Lp-metric with p→∞.The Voronoi region generated with this distance is the same shape, namely square, as the one generated with Manhattan metrics (as shown in Figure 3c).However, the differences between these two metrics are the directions of the coordinate axes.The axes of the Manhattan distance are inclined 45 degrees, while the axes of the Chebyshev distance are horizontal and vertical (as shown in Figure 4).In this paper, we use the Chebyshev distance because its axial directions conform with normal image space, which is easier to scale and control.
As shown in Figure 4, the side of the square is defined as La, while the diagonal line is defined as Lb; the axis xy corresponding to the Manhattan distance is shown in Figure 4a

Chebyshev Distance with Axial Scale
The Chebyshev distance between p0( ,  ) and p1( ,  ) with axial scale is as follows: With this metric, a rectangle-shaped Voronoi diagram is generated (as shown in Figure 3d).To obtain the Chebyshev distance with axial scale, the original Chebyshev space is scaled along two axial directions; as a result, the unit distance has been changed and the original square-shaped Voronoi diagram becomes a rectangle-shaped one.
As shown in Figure 5, let  be 1/2 and the length of each square side be La before the axial scaling.For the rectangle after the scaling, the length is 2 × La, while the width is 1/2La.The distance from any point ( ,  ) on the rectangle to the origin point ( ,  ) is (| −  |, (1/)| −  |).To measure this distance, a correspondence between points A and point  ( , (1/) ) on the square is defined, and the distance between point  and the origin point is calculated based on the original Chebyshev space.Hence, the distance from any point ( ,  ) on the rectangle to the origin point ( ,  ) is La/2.

Chebyshev Distance with Axial Scale
The Chebyshev distance between p 0 (x 0 , y 0 ) and p 1 (x 1 , y 1 ) with axial scale is as follows: With this metric, a rectangle-shaped Voronoi diagram is generated (as shown in Figure 3d).To obtain the Chebyshev distance with axial scale, the original Chebyshev space is scaled along two axial directions; as a result, the unit distance has been changed and the original square-shaped Voronoi diagram becomes a rectangle-shaped one.
As shown in Figure 5, let β be 1/2 and the length of each square side be L a before the axial scaling.For the rectangle after the scaling, the length is 2 × L a , while the width is 1/2L a .The distance from any point A(x 1 , y 1 ) on the rectangle to the origin point O(x 0 , To measure this distance, a correspondence between points A and point A 1 (βx 1 , (1/β)y 1 ) on the square is defined, and the distance between point A 1 and the origin point is calculated based on the original Chebyshev space.Hence, the distance from any point A(x 1 , y 1 ) on the rectangle to the origin point O(x 0 , y 0 ) is L a /2.

L2-Metric Distance with Axial Scale
L2-metric with axial scale can be expressed as follows: With this metric, the generated Voronoi region is close to an oval shape (as shown in Figure 6).In the same manner as the Chebyshev distance with axial scale mentioned above, for the L2-metric distance with axial scale, the original Euclidean space is scaled along two axial directions; as a result, the unit distance is changed and the original circle-shaped Voronoi diagram becomes an oval-shaped one.
As shown in Figure 6, let  be 1/2 and the radius of the circle be .With the transverse length of the oval defined as 4×r after the axial scaling and the minor-axis length of the oval as r, the distance between any point ( ,  ) on the oval and the origin point ( ,  ) equals that between the point  and origin point ( ,  ). is the corresponding point of  on the circle, whose coordinates are  and  multiplied by the factors  and ( ) respectively.This means that the distance between any point on the oval and the origin point ( ,  ) is r.

Equidistant Line Distance
The equidistant line distance is a kind of metric adopted in the paper to generate regular triangular regions.On the basis of this metric, the distance between two points is measured by the equidistant line of a regular triangle.The direction of the regular triangle in the space is defined as shown in Figure 7.Given a regular triangle with a side length of La, with the origin position defined as the center of gravity, draw a similar regular triangle within the original regular triangle with the centers of gravity of all the triangles at the origin point, and then the distance between any point on the triangle side and the origin point is defined as the distance between the center of gravity and the triangle vertex.When measuring the distance between ( ,  ) and ( ,  ) with this metric, we need to define the origin point as P, construct regular triangles with P as the gravity center, and let the point Q be located on the side of the regular triangle.Hence, the distance between P and the vertex of the triangle is the same as that between P and Q.
The distance between any point ( ,  ) on the triangle and the origin point ( ,  ) is √3/3.However, instead of being a strict triangle, the final Voronoi region obtained with this With this metric, the generated Voronoi region is close to an oval shape (as shown in Figure 6).In the same manner as the Chebyshev distance with axial scale mentioned above, for the L 2 -metric distance with axial scale, the original Euclidean space is scaled along two axial directions; as a result, the unit distance is changed and the original circle-shaped Voronoi diagram becomes an oval-shaped one.With this metric, the generated Voronoi region is close to an oval shape (as shown in Figure 6).In the same manner as the Chebyshev distance with axial scale mentioned above, for the L2-metric distance with axial scale, the original Euclidean space is scaled along two axial directions; as a result, the unit distance is changed and the original circle-shaped Voronoi diagram becomes an oval-shaped one.
As shown in Figure 6, let  be 1/2 and the radius of the circle be .With the transverse length of the oval defined as 4×r after the axial scaling and the minor-axis length of the oval as r, the distance between any point ( ,  ) on the oval and the origin point ( ,  ) equals that between the point  and origin point ( ,  ). is the corresponding point of  on the circle, whose coordinates are  and  multiplied by the factors  and ( ) respectively.This means that the distance between any point on the oval and the origin point ( ,  ) is r.

Equidistant Line Distance
The equidistant line distance is a kind of metric adopted in the paper to generate regular triangular regions.On the basis of this metric, the distance between two points is measured by the equidistant line of a regular triangle.The direction of the regular triangle in the space is defined as shown in Figure 7.Given a regular triangle with a side length of La, with the origin position defined as the center of gravity, draw a similar regular triangle within the original regular triangle with the centers of gravity of all the triangles at the origin point, and then the distance between any point on the triangle side and the origin point is defined as the distance between the center of gravity and the triangle vertex.When measuring the distance between ( ,  ) and ( ,  ) with this metric, we need to define the origin point as P, construct regular triangles with P as the gravity center, and let the point Q be located on the side of the regular triangle.Hence, the distance between P and the vertex of the triangle is the same as that between P and Q.
The distance between any point ( ,  ) on the triangle and the origin point ( ,  ) is √3/3.However, instead of being a strict triangle, the final Voronoi region obtained with this As shown in Figure 6, let β be 1/2 and the radius of the circle be r.With the transverse length of the oval defined as 4 × r after the axial scaling and the minor-axis length of the oval as r, the distance between any point A(x 1 , y 1 ) on the oval and the origin point O(x 0 , y 0 ) equals that between the point A 1 and origin point O(x 0 , y 0 ).A 1 is the corresponding point of A on the circle, whose coordinates are x 1 and y 1 multiplied by the factors β and 1 β respectively.This means that the distance between any point on the oval and the origin point O(x 0 , y 0 ) is r.

Equidistant Line Distance
The equidistant line distance is a kind of metric adopted in the paper to generate regular triangular regions.On the basis of this metric, the distance between two points is measured by the equidistant line of a regular triangle.The direction of the regular triangle in the space is defined as shown in Figure 7.Given a regular triangle with a side length of L a , with the origin position defined as the center of gravity, draw a similar regular triangle within the original regular triangle with the centers of gravity of all the triangles at the origin point, and then the distance between any point on the triangle side and the origin point is defined as the distance between the center of gravity and the triangle vertex.When measuring the distance between P(x 1 , y 1 ) and Q(x 0 , y 0 ) with this metric, we need to define the origin point as P, construct regular triangles with P as the gravity center, and let the point Q be located on the side of the regular triangle.Hence, the distance between P and the vertex of the triangle is the same as that between P and Q.
The distance between any point A(x 1 , y 1 ) on the triangle and the origin point O(x 0 , y 0 ) is √ 3/3La.However, instead of being a strict triangle, the final Voronoi region obtained with this metric is a dodecagon that is close to a triangle (as shown in Figure 7, right).The reason is that triangles with consistent directions cannot fill the entire plane with no gaps or overlaps.Therefore, during the repeating iteration process, the triangle is turned into the shape shown in Figure 7, left.Finding an approach to obtain a more accurate triangular region partition will be our future research direction.metric is a dodecagon that is close to a triangle (as shown in Figure 7, right).The reason is that triangles with consistent directions cannot fill the entire plane with no gaps or overlaps.Therefore, during the repeating iteration process, the triangle is turned into the shape shown in Figure 7, left.
Finding an approach to obtain a more accurate triangular region partition will be our future research direction.

Voronoi Diagram Generation Based on GPU Acceleration
To generate a Voronoi diagram, the iteration process must be repeated and every iteration process needs to loop through all the distribution points, which is quite time-consuming.Hoff et al. [30] proposed a method of converting the Voronoi generation process into a graphical rendering process.First, the distribution of the initial points is determined, and then a circular cone with height equal to the radius is constructed at each initial point location.The vertices of the cones are used as the initial seed points.The cones are the same in terms of all parameters except color.In this way, the height of the side part of the cone can be adopted to measure the distance from the corresponding points on the bottom surface to the center; thus, the correct Voronoi diagram can be obtained by rendering all the cones with visibility detection via the z-buffer algorithm.The Voronoi generation process is successfully transformed into a graphical rendering process, and the GPU parallel computing ability can be fully utilized to increase the speed.This paper, taking this as a reference, implemented a GPU rendering process to achieve quick generation of Voronoi regions with many shapes.In the following part of this paper, the Voronoi regions generated by GPU rendering are described-they are of five shapes; namely, circle, square, rectangle, oval, and triangle.
The key to accomplishing the Voronoi generation via geometric rendering and visibility detection with the z-buffer algorithm is to ensure that the height of any point on the side of the geometric shape is equal to the distance between the corresponding points of the random point on the bottom surface and the center.Taking the cone shown in Figure 8 9, by rendering a mass of three-dimensional cones, we can obtain the Voronoi diagram.In the rendering, the orthographic projection is adopted, and the z-buffer method is used for visibility detection.

Voronoi Diagram Generation Based on GPU Acceleration
To generate a Voronoi diagram, the iteration process must be repeated and every iteration process needs to loop through all the distribution points, which is quite time-consuming.Hoff et al. [30] proposed a method of converting the Voronoi generation process into a graphical rendering process.First, the distribution of the initial points is determined, and then a circular cone with height equal to the radius is constructed at each initial point location.The vertices of the cones are used as the initial seed points.The cones are the same in terms of all parameters except color.In this way, the height of the side part of the cone can be adopted to measure the distance from the corresponding points on the bottom surface to the center; thus, the correct Voronoi diagram can be obtained by rendering all the cones with visibility detection via the z-buffer algorithm.The Voronoi generation process is successfully transformed into a graphical rendering process, and the GPU parallel computing ability can be fully utilized to increase the speed.This paper, taking this as a reference, implemented a GPU rendering process to achieve quick generation of Voronoi regions with many shapes.In the following part of this paper, the Voronoi regions generated by GPU rendering are described-they are of five shapes; namely, circle, square, rectangle, oval, and triangle.
The key to accomplishing the Voronoi generation via geometric rendering and visibility detection with the z-buffer algorithm is to ensure that the height of any point on the side of the geometric shape is equal to the distance between the corresponding points of the random point on the bottom surface and the center.Taking the cone shown in Figure 8 as an example, the height of the cone and the radius of the bottom surface are both L a , and the corresponding point of any random point A on the side is A 0 .Given that the distance between A 0 and center is the yellow segment |OA 0 | and the height of A is the blue segment |AB|, based on |TB| = |AB| and |TB| = |OA 0 |, we conclude that |AB| = |OA 0 |.As shown in Figure 9, by rendering a mass of three-dimensional cones, we can obtain the Voronoi diagram.In the rendering, the orthographic projection is adopted, and the z-buffer method is used for visibility detection.
the bottom surface and the center.Taking the cone shown in Figure 8 9, by rendering a mass of three-dimensional cones, we can obtain the Voronoi diagram.In the rendering, the orthographic projection is adopted, and the z-buffer method is used for visibility detection.To generate Voronoi diagrams in many shapes, we adopt many kinds of geometries that conform to the above distance conditions; the detailed information is shown in Table 1.To generate Voronoi diagrams in many shapes, we adopt many kinds of geometries that conform to the above distance conditions; the detailed information is shown in Table 1.Step 2. Show geometries (circular cone, square pyramid, pyramid, elliptic cone, and triangular pyramid) in different colors at every p(x, y) location; different colors can make it easier to distinguish different Voronoi regions.
Step 3. Rotate the geometry ψ(x, y) according to the direction field of p(x, y).
Step 4. Perform the depth measurement, and read the current drawing results; calculate the center of gravity of each Voronoi region; the results are the new distribution point set p.
Step 5. Repeat steps 2 to 4 for the new point set p until the specified number of iterations is reached.
Step 6. Place the elements at the locations of the final centers of gravity.

Parametric Control of Element Distribution
To make the synthesis output more variable and controllable, this paper introduces several approaches to control the size and direction of the generated Voronoi regions, thus achieving To generate Voronoi diagrams in many shapes, we adopt many kinds of geometries that conform to the above distance conditions; the detailed information is shown in Table 1.Step 2. Show geometries (circular cone, square pyramid, pyramid, elliptic cone, and triangular pyramid) in different colors at every p(x, y) location; different colors can make it easier to distinguish different Voronoi regions.
Step 3. Rotate the geometry ψ(x, y) according to the direction field of p(x, y).
Step 4. Perform the depth measurement, and read the current drawing results; calculate the center of gravity of each Voronoi region; the results are the new distribution point set p.
Step 5. Repeat steps 2 to 4 for the new point set p until the specified number of iterations is reached.
Step 6. Place the elements at the locations of the final centers of gravity.

Parametric Control of Element Distribution
To make the synthesis output more variable and controllable, this paper introduces several approaches to control the size and direction of the generated Voronoi regions, thus achieving

Manhattan distance
Side length = La To generate Voronoi diagrams in many shapes, we adopt many kinds of geometries that conform to the above distance conditions; the detailed information is shown in Table 1.Step 2. Show geometries (circular cone, square pyramid, pyramid, elliptic cone, and triangular pyramid) in different colors at every p(x, y) location; different colors can make it easier to distinguish different Voronoi regions.
Step 3. Rotate the geometry ψ(x, y) according to the direction field of p(x, y).
Step 4. Perform the depth measurement, and read the current drawing results; calculate the center of gravity of each Voronoi region; the results are the new distribution point set p.
Step 5. Repeat steps 2 to 4 for the new point set p until the specified number of iterations is reached.
Step 6. Place the elements at the locations of the final centers of gravity.

Parametric Control of Element Distribution
To make the synthesis output more variable and controllable, this paper introduces several To generate Voronoi diagrams in many shapes, we adopt many kinds of geometries that conform to the above distance conditions; the detailed information is shown in Table 1.Step 2. Show geometries (circular cone, square pyramid, pyramid, elliptic cone, and triangular pyramid) in different colors at every p(x, y) location; different colors can make it easier to distinguish different Voronoi regions.
Step 3. Rotate the geometry ψ(x, y) according to the direction field of p(x, y).
Step 4. Perform the depth measurement, and read the current drawing results; calculate the center of gravity of each Voronoi region; the results are the new distribution point set p.
Step 5. Repeat steps 2 to 4 for the new point set p until the specified number of iterations is reached.
Step 6. Place the elements at the locations of the final centers of gravity.

Parametric Control of Element Distribution
To make the synthesis output more variable and controllable, this paper introduces several To generate Voronoi diagrams in many shapes, we adopt many kinds of geometries that conform to the above distance conditions; the detailed information is shown in Table 1.Step 2. Show geometries (circular cone, square pyramid, pyramid, elliptic cone, and triangular pyramid) in different colors at every p(x, y) location; different colors can make it easier to distinguish different Voronoi regions.
Step 3. Rotate the geometry ψ(x, y) according to the direction field of p(x, y).
Step 4. Perform the depth measurement, and read the current drawing results; calculate the center of gravity of each Voronoi region; the results are the new distribution point set p.
Step 5. Repeat steps 2 to 4 for the new point set p until the specified number of iterations is reached.
Step 6. Place the elements at the locations of the final centers of gravity.

Parametric Control of Element Distribution
To make the synthesis output more variable and controllable, this paper introduces several To generate Voronoi diagrams in many shapes, we adopt many kinds of geometries that conform to the above distance conditions; the detailed information is shown in Table 1.

Calculation process:
Step 1. Randomly generate initial distribution point set p.
Step 2. Show geometries (circular cone, square pyramid, pyramid, elliptic cone, and triangular pyramid) in different colors at every p(x, y) location; different colors can make it easier to distinguish different Voronoi regions.
Step 3. Rotate the geometry ψ(x, y) according to the direction field of p(x, y).
Step 4. Perform the depth measurement, and read the current drawing results; calculate the center of gravity of each Voronoi region; the results are the new distribution point set p.
Step 5. Repeat steps 2 to 4 for the new point set p until the specified number of iterations is reached.
Step 6. Place the elements at the locations of the final centers of gravity.

Parametric Control of Element Distribution
To make the synthesis output more variable and controllable, this paper introduces several approaches to control the size and direction of the generated Voronoi regions, thus achieving Calculation process: Step 1. Randomly generate initial distribution point set p.
Step 2. Show geometries (circular cone, square pyramid, pyramid, elliptic cone, and triangular pyramid) in different colors at every p(x, y) location; different colors can make it easier to distinguish different Voronoi regions.
Step 3. Rotate the geometry ψ(x, y) according to the direction field of p(x, y).
Step 4. Perform the depth measurement, and read the current drawing results; calculate the center of gravity of each Voronoi region; the results are the new distribution point set p.
Step 5. Repeat steps 2 to 4 for the new point set p until the specified number of iterations is reached.
Step 6. Place the elements at the locations of the final centers of gravity.To generate Voronoi diagrams in many shapes, we adopt many kinds of geometries that conform to the above distance conditions; the detailed information is shown in Table 1.

Calculation process:
Step 1. Randomly generate initial distribution point set p.
Step 2. Show geometries (circular cone, square pyramid, pyramid, elliptic cone, and triangular pyramid) in different colors at every p(x, y) location; different colors can make it easier to distinguish different Voronoi regions.
Step 3. Rotate the geometry ψ(x, y) according to the direction field of p(x, y).
Step 4. Perform the depth measurement, and read the current drawing results; calculate the center of gravity of each Voronoi region; the results are the new distribution point set p.
Step 5. Repeat steps 2 to 4 for the new point set p until the specified number of iterations is reached.
Step 6. Place the elements at the locations of the final centers of gravity.

Parametric Control of Element Distribution
To make the synthesis output more variable and controllable, this paper introduces several approaches to control the size and direction of the generated Voronoi regions, thus achieving compact element distributions of changeable size and direction.

Parametric Control of Element Distribution
To make the synthesis output more variable and controllable, this paper introduces several approaches to control the size and direction of the generated Voronoi regions, thus achieving compact element distributions of changeable size and direction.

Changing of Element Sizes
Our research results indicate that the distribution of the initial point set has a significant effect on the synthesis results.A region with a more compact distribution of initial points usually eventually results in a smaller Voronoi region.Therefore, we propose to control the sizes of the Voronoi polygons by adjusting the parameters while generating the initial points to control the density.However, the experimental results show that if we only use the method described above to control the size, there will be situations in which the size gradually varies in the junctions of large and small polygons.This will create irregular polygons at the boundary (as shown in Figure 11a), which is not desirable.For this reason, in addition to controlling the density of the initial points, we utilize the method of changing the size of the geometry that was adopted during the synthesis process for better results.
Let α be the slope of the geometry used during the synthesis process, and take the circular cone in Figure 10 as an example.While generating the Voronoi region, use geometry with a smaller size (geometry with a larger α) when a smaller region needs to be obtained and geometry with a larger size (geometry with a smaller α) in other regions.Together with the initial point distribution in this region, better results can be achieved by changing the sizes of the elements (as shown in Figure 11b).results in a smaller Voronoi region.Therefore, we propose to control the sizes of the Voronoi polygons by adjusting the parameters while generating the initial points to control the density.However, the experimental results show that if we only use the method described above to control the size, there will be situations in which the size gradually varies in the junctions of large and small polygons.This will create irregular polygons at the boundary (as shown in 11a), which is not desirable.For this reason, in addition to controlling the density of the initial points, we utilize the method of changing the size of the geometry that was adopted during the synthesis process for better results.
Let  be the slope of the geometry used during the synthesis process, and take the circular cone in Figure 10 as an example.While generating the Voronoi region, use geometry with a smaller size (geometry with a larger ) when a smaller region needs to be obtained and geometry with a larger size (geometry with a smaller ) in other regions.Together with the initial point distribution in this region, better results can be achieved by changing the sizes of the elements (as shown in 11b).

Changing of the Vector Field
During the iterative optimization process of generating the centroidal Voronoi diagram, the method in this paper can adjust the direction of the generated Voronoi region, which is shown in Step 3 of the calculation process in Chapter 3.2.Users can define a field ψ for the distribution region.Let the direction of the vector field of seed point p(x,y) be ψ(x,y).By rotating the geometry at p(x,y) to the direction ψ(x,y), we can obtain element distribution results by conforming to the vector field.Figure 12a is the vector field and Figure 12b is the generated Voronoi diagram based on the Chebyshev distance metric.results in a smaller Voronoi region.Therefore, we propose to control the sizes of the Voronoi polygons by adjusting the parameters while generating the initial points to control the density.However, the experimental results show that if we only use the method described above to control the size, there will be situations in which the size gradually varies in the junctions of large and small polygons.This will create irregular polygons at the boundary (as shown in 11a), which is not desirable.For this reason, in addition to controlling the density of the initial points, we utilize the method of changing the size of the geometry that was adopted during the synthesis process for better results.
Let  be the slope of the geometry used during the synthesis process, and take the circular cone in Figure 10 as an example.While generating the Voronoi region, use geometry with a smaller size (geometry with a larger ) when a smaller region needs to be obtained and geometry with a larger size (geometry with a smaller ) in other regions.Together with the initial point distribution in this region, better results can be achieved by changing the sizes of the elements (as shown in 11b).

Changing of the Vector Field
During the iterative optimization process of generating the centroidal Voronoi diagram, the method in this paper can adjust the direction of the generated Voronoi region, which is shown in Step 3 of the calculation process in Chapter 3.2.Users can define a vector field ψ for the distribution region.Let the direction of the vector field of seed point p(x,y) be ψ(x,y).By rotating the geometry at p(x,y) to the direction ψ(x,y), we can obtain element distribution results by conforming to the vector field.Figure 12a is the vector field and Figure 12b is the generated Voronoi diagram based on the Chebyshev distance metric.

Changing of the Vector Field
During the iterative optimization process of generating the centroidal Voronoi diagram, the method in this paper can adjust the direction of the generated Voronoi region, which is shown in Step 3 of the calculation process in Chapter 3.2.Users can define a vector field ψ for the distribution region.Let the direction of the vector field of seed point p(x,y) be ψ(x,y).By rotating the geometry at p(x,y) to the direction ψ(x,y), we can obtain element distribution results by conforming to the vector field.Figure 12a is the vector field and Figure 12b

Placement of the Elements
Through the above methods, compact distributions of Voronoi regions with different shapes can be obtained.In order to realize the compact distributions of the discrete elements, we need to place the elements into every Voronoi region.However, because the Voronoi regions have different sizes and directions, the sizes and directions of the elements need to be adjusted while putting them into the regions to obtain the distribution results.Taking the size adjustment as an example, as shown in Figure 13, we need to extract the contour of every Voronoi region and adjust the size and direction based on the contour information.

Results
We use C++ programming to realize the algorithm and, based on OpenGL, the Voronoi diagram is successfully generated through three-dimensional graphical rendering.

Experimental Results
The pictures shown in Figure 14 are synthesis results that were obtained using the methods proposed in this paper.Figures 14a-d show distribution effects of squares, rectangles, ellipses, and triangles, respectively.Figures 14e,f show the compact distribution effects with elements varying from sizes and directions.

Placement of the Elements
Through the above methods, compact distributions of Voronoi regions with different shapes can be obtained.In order to realize the compact distributions of the discrete elements, we need to place the elements into every Voronoi region.However, because the Voronoi regions have different sizes and directions, the sizes and directions of the elements need to be adjusted while putting them into the regions to obtain the distribution results.Taking the size adjustment as an example, as shown in Figure 13, we need to extract the contour of every Voronoi region and adjust the size and direction based on the contour information.

Placement of the Elements
Through the above methods, compact distributions of Voronoi regions with different shapes can be obtained.In order to realize the compact distributions of the discrete elements, we need to place the elements into every Voronoi region.However, because the Voronoi regions have different sizes and directions, the sizes and directions of the elements need to be adjusted while putting them into the regions to obtain the distribution results.Taking the size adjustment as an example, as shown in Figure 13, we need to extract the contour of every Voronoi region and adjust the size and direction based on the contour information.

Results
We use C++ programming to realize the algorithm and, based on OpenGL, the Voronoi diagram is successfully generated through three-dimensional graphical rendering.

Experimental Results
The pictures shown in Figure 14 are synthesis results that were obtained using the methods proposed in this paper.Figures 14a-d show distribution effects of squares, rectangles, ellipses, and triangles, respectively.Figures 14e,f show the compact distribution effects with elements varying from sizes and directions.

Results
We use C++ programming to realize the algorithm and, based on OpenGL, the Voronoi diagram is successfully generated through three-dimensional graphical rendering.

Experimental Results
The pictures shown in Figure 14 are synthesis results that were obtained using the methods proposed in this paper.Figure 14a-d show distribution effects of squares, rectangles, ellipses, and triangles, respectively.Figure 14e,f show the compact distribution effects with elements varying from sizes and directions.In particular, if random points are distributed throughout an image and the image information is used as a constraint, our algorithm can simulate the mosaic effect based on the image features, which means for any picture specified by the user, our algorithm can automatically convert the image into a mosaic-style image and obtain more artistic mosaic effects.As shown in Figures 16-19

Synthesis Effect Comparison
This paper conducts research on the compact distributions of discrete elements.The more compact the distribution is in a certain region, the higher the regional coverage scale.Therefore, the synthesis effect could be evaluated by the regional coverage scale.Methods of FritzscheL, P. [12], Smith, K. [26] and Hu, W. [13] are adopted to compare the synthesis effects from the following three aspects: synthesis result, coverage scale, and time consumption.The synthesis results of the method introduced in this paper are closely related to the number of iterations used to generate the centroidal Voronoi diagram.Table 2 shows the results after 20 iterations.On the basis of the data in Table 2, the method of this paper is characterized by a larger coverage scale, lower time consumption, and a better synthesis result.

Synthesis Effect Comparison
This paper conducts research on the compact distributions of discrete elements.The more compact the distribution is in a certain region, the higher the regional coverage scale.Therefore, the synthesis effect could be evaluated by the regional coverage scale.Methods of FritzscheL, P. [12], Smith, K. [26] and Hu, W. [13] are adopted to compare the synthesis effects from the following three aspects: synthesis result, coverage scale, and time consumption.The synthesis results of the method introduced in this paper are closely related to the number of iterations used to generate the centroidal Voronoi diagram.Table 2 shows the results after 20 iterations.On the basis of the data in Table 2, the method of this paper is characterized by a larger coverage scale, lower time consumption, and a better synthesis result.

Synthesis Effect Comparison
This paper conducts research on the compact distributions of discrete elements.The more compact the distribution is in a certain region, the higher the regional coverage scale.Therefore, the synthesis effect could be evaluated by the regional coverage scale.Methods of FritzscheL, P. [12], Smith, K. [26] and Hu, W. [13] are adopted to compare the synthesis effects from the following three aspects: synthesis result, coverage scale, and time consumption.The synthesis results of the method introduced in this paper are closely related to the number of iterations used to generate the centroidal Voronoi diagram.Table 2 shows the results after 20 iterations.On the basis of the data in Table 2, the method of this paper is characterized by a larger coverage scale, lower time consumption, and a better synthesis result.

Synthesis Effect Comparison
This paper conducts research on the compact distributions of discrete elements.The more compact the distribution is in a certain region, the higher the regional coverage scale.Therefore, the synthesis effect could be evaluated by the regional coverage scale.Methods of FritzscheL, P. [12], Smith, K. [26] and Hu, W. [13] are adopted to compare the synthesis effects from the following three aspects: synthesis result, coverage scale, and time consumption.The synthesis results of the method introduced in this paper are closely related to the number of iterations used to generate the centroidal Voronoi diagram.Table 2 shows the results after 20 iterations.On the basis of the data in Table 2, the method of this paper is characterized by a larger coverage scale, lower time consumption, and a better synthesis result.

Synthesis Effect Comparison
This paper conducts research on the compact distributions of discrete elements.The more compact the distribution is in a certain region, the higher the regional coverage scale.Therefore, the synthesis effect could be evaluated by the regional coverage scale.Methods of FritzscheL, P. [12], Smith, K. [26] and Hu, W. [13] are adopted to compare the synthesis effects from the following three aspects: synthesis result, coverage scale, and time consumption.The synthesis results of the method introduced in this paper are closely related to the number of iterations used to generate the centroidal Voronoi diagram.Table 2 shows the results after 20 iterations.On the basis of the data in Table 2, the method of this paper is characterized by a larger coverage scale, lower time consumption, and a better synthesis result.

Synthesis Effect Comparison
This paper conducts research on the compact distributions of discrete elements.The more compact the distribution is in a certain region, the higher the regional coverage scale.Therefore, the synthesis effect could be evaluated by the regional coverage scale.Methods of FritzscheL, P. [12], Smith, K. [26] and Hu, W. [13] are adopted to compare the synthesis effects from the following three aspects: synthesis result, coverage scale, and time consumption.The synthesis results of the method introduced in this paper are closely related to the number of iterations used to generate the centroidal Voronoi diagram.Table 2 shows the results after 20 iterations.On the basis of the data in Table 2, the method of this paper is characterized by a larger coverage scale, lower time consumption, and a better synthesis result.

Figure 1 .
Figure 1.Synthesis effect of the work in this paper: (a) square elements; (b) circular elements.

Figure 1 .
Figure 1.Synthesis effect of the work in this paper: (a) square elements; (b) circular elements.

Figure 3 .
Figure 3. Voronoi diagrams generated with different distance function: (a) L 2 -metric distance; (b) Manhattan distance; (c) Chebyshev distance; (d) Chebyshev distance with axial scale; (e) L 2 -metric with axial scale; and (f) equidistant line distance.The partitioned regions are presented in different colors.

Figure 4 .
Figure 4.The axial directions of the Manhattan distance and Chebyshev distance: (a) Manhattan distance; (b) Chebyshev distance.

Figure 4 .
Figure 4.The axial directions of the Manhattan distance and Chebyshev distance: (a) Manhattan distance; (b) Chebyshev distance.

Algorithms 2019 , 16 Figure 5 .
Figure 5. Different appearances of the square in Chebyshev space and Chebyshev space with axial scale.

Figure 6 .
Figure 6.Different appearances of the circle in L2-metric distance and L2-metric distance with axial scale.

Figure 5 .
Figure 5. Different appearances of the square in Chebyshev space and Chebyshev space with axial scale.3.1.3.L 2 -Metric Distance with Axial Scale L 2 -metric with axial scale can be expressed as follows:

16 Figure 5 .
Figure 5. Different appearances of the square in Chebyshev space and Chebyshev space with axial scale.

Figure 6 .
Figure 6.Different appearances of the circle in L2-metric distance and L2-metric distance with axial scale.

Figure 6 .
Figure 6.Different appearances of the circle in L 2 -metric distance and L 2 -metric distance with axial scale.

Figure 7 .
Figure 7. Left: equidistant line distance.Right: the Voronoi diagram and the Voronoi region shape (lower right corner) by the equidistant line distance.
as an example, the height of the cone and the radius of the bottom surface are both La, and the corresponding point of any random point A on the side is  .Given that the distance between  and center is the yellow segment | | and the height of A is the blue segment ||, based on || = || and || = | |, we conclude that || = | |.As shown in Figure

Figure 8 .
Figure 8.The size of the cones.

Figure 7 .
Figure 7. Left: equidistant line distance.Right: the Voronoi diagram and the Voronoi region shape (lower right corner) by the equidistant line distance.
as an example, the height of the cone and the radius of the bottom surface are both La, and the corresponding point of any random point A on the side is  .Given that the distance between  and center is the yellow segment | | and the height of A is the blue segment ||, based on || = || and || = | |, we conclude that || = | |.As shown in Figure

Figure 8 .
Figure 8.The size of the cones.

Figure 8 .
Figure 8.The size of the cones.

Figure 11 .
Figure 11.Results of controlling the sizes of the Voronoi regions (a) by simply controlling the initial point density and (b) by controlling the initial point density and geometric slopes.

Figure 11 .
Figure 11.Results of controlling the sizes of the Voronoi regions (a) by simply controlling the initial point density and (b) by controlling the initial point density and geometric slopes.

Figure 11 .
Figure 11.Results of controlling the sizes of the Voronoi regions (a) by simply controlling the initial point density and (b) by controlling the initial point density and geometric slopes.
is the generated Voronoi diagram based on the Chebyshev distance metric.

Figure 12 .
Figure 12.Centroidal Voronoi diagram with vector field: (a) vector field; (b) centroidal Voronoi diagram with the vector field shown in (a) based on the Chebyshev distance metric.

Figure 13 .
Figure 13.Putting the elements into Voronoi regions after adjusting the sizes.

Figure 12 .
Figure 12.Centroidal Voronoi diagram with vector field: (a) vector field; (b) centroidal Voronoi diagram with the vector field shown in (a) based on the Chebyshev distance metric.

Figure 12 .
Figure 12.Centroidal Voronoi diagram with vector field: (a) vector field; (b) centroidal Voronoi diagram with the vector field shown in (a) based on the Chebyshev distance metric.

Figure 13 .
Figure 13.Putting the elements into Voronoi regions after adjusting the sizes.

Figure 13 .
Figure 13.Putting the elements into Voronoi regions after adjusting the sizes.

Figure 14 .
Figure 14.Demonstration of synthesis results: (a) distribution effects of squares; (b) rectangles; (c) ellipses; (d) triangles; (e) the distribution effects of elements with varying sizes; (f) the distribution effects of elements with varying directions.Furthermore, the methods introduced in this paper can simulate the artistic effects of mosaic works, as shown in Figure 15.

Figure 15 .
Figure 15.Synthesis results in specified regions.

Figure 14 . 16 Figure 14 .
Figure 14.Demonstration of synthesis results: (a) distribution effects of squares; (b) rectangles; (c) ellipses; (d) triangles; (e) the distribution effects of elements with varying sizes; (f) the distribution effects of elements with varying directions.Furthermore, the methods introduced in this paper can simulate the artistic effects of mosaic works, as shown in Figure 15.

Figure 15 .
Figure 15.Synthesis results in specified regions.

Figure 15 .
Figure 15.Synthesis results in specified regions.
below, the mosaic effects are synthesized based on the image features.

Figure 16 .
Figure 16.Converting the image into mosaic effects.The vector field controls the directions of the elements; the density graph controls the sizes of the elements.

Figure 17 .
Figure 17.Synthesis results based on image feature analysis.

Figure 18 .
Figure 18.Synthesis results without (middle) and with (right) the vector field.(The element's directions in the right image are controlled by the vector field extracted from the left image.).

Figure 16 . 16 Figure 16 .
Figure 16.Converting the image into mosaic effects.The vector field controls the directions of the elements; the density graph controls the sizes of the elements.

Figure 17 .
Figure 17.Synthesis results based on image feature analysis.

Figure 18 .
Figure 18.Synthesis results without (middle) and with (right) the vector field.(The element's directions in the right image are controlled by the vector field extracted from the left image.).

Figure 17 .
Figure 17.Synthesis results based on image feature analysis.

Figure 18 .
Figure 18.Synthesis results without (middle) and with (right) the vector field.(The element's directions in the right image are controlled by the vector field extracted from the left image.).

Figure 18 .
Figure 18.Synthesis results without (middle) and with (right) the vector field.(The element's directions in the right image are controlled by the vector field extracted from the left image.).

Figure 19 .
Figure 19.Synthesis results without (middle) and with (right) size variation of the elements (the size variation of the right image is controlled by the density graph extracted from the left image).

Figure 19 .
Figure 19.Synthesis results without (middle) and with (right) size variation of the elements (the size variation of the right image is controlled by the density graph extracted from the left image).

Figure 19 .
Figure 19.Synthesis results without (middle) and with (right) size variation of the elements (the size variation of the right image is controlled by the density graph extracted from the left image).

Figure 19 .
Figure 19.Synthesis results without (middle) and with (right) size variation of the elements (the size variation of the right image is controlled by the density graph extracted from the left image).

Figure 19 .
Figure 19.Synthesis results without (middle) and with (right) size variation of the elements (the size variation of the right image is controlled by the density graph extracted from the left image).

Figure 19 .
Figure 19.Synthesis results without (middle) and with (right) size variation of the elements (the size variation of the right image is controlled by the density graph extracted from the left image).

Table 1 .
Shapes and sizes of geometries corresponding to different distance metrics.

Table 1 .
Shapes and sizes of geometries corresponding to different distance metrics.
Calculation process:Step 1. Randomly generate initial distribution point set p.

Table 1 .
Shapes and sizes of geometries corresponding to different distance metrics.
Calculation process:Step 1. Randomly generate initial distribution point set p.

Table 1 .
Shapes and sizes of geometries corresponding to different distance metrics.
Calculation process:Step 1. Randomly generate initial distribution point set p.

Table 1 .
Shapes and sizes of geometries corresponding to different distance metrics.
Calculation process:Step 1. Randomly generate initial distribution point set p.

Table 1 .
Shapes and sizes of geometries corresponding to different distance metrics.
Calculation process:Step 1. Randomly generate initial distribution point set p.

Table 1 .
Shapes and sizes of geometries corresponding to different distance metrics.

Table 1 .
Shapes and sizes of geometries corresponding to different distance metrics.

Table 2 .
Comparison among different methods.

Table 2 .
Comparison among different methods.

Table 2 .
Comparison among different methods.

picture Average coverage scale Average time
Algorithms 2019, 12, x FOR PEER REVIEW 14 of 16

Table 2 .
Comparison among different methods.

picture Average coverage scale Average time
Algorithms 2019, 12, x FOR PEER REVIEW 14 of 16

Table 2 .
Comparison among different methods.

picture Average coverage scale Average time
Algorithms 2019, 12, x FOR PEER REVIEW 14 of 16

Table 2 .
Comparison among different methods.