1 Introduction
Low isometric distortion parameterization plays a fundamental role in many computer graphics and geometric processing tasks, such as atlas generation [1, 2], remeshing [3, 4], and bijective surface correspondence [5]. Fig. 1 shows three applications. Before being parameterized to the plane, a closed mesh needs to be cut to a disk topology. The distortion of a parameterization is highly related to not only the parameterization method, but also the position of the cut. Since a surface is generally not developable, the isometric parameterization distortion is usually high, unless the cut passes through some vertices, where the distortion is usually concentrated. We called these vertices distortion points.
A good distortion point detection algorithm usually satisfies the following properties: (1) the algorithm is automatic and does not require manual intervention; (2) the parameterizations exhibit low isometric distortion if the cut path passes through the detected distortion points; (3) the number of distortion points is small; (4) and the algorithm is efficient, otherwise manual selection is preferable. On the one hand, more distortion points tend to reduce isometric distortion, but too many distortion points are very likely to have a negative impact on subsequent applications. For example, more distortion points usually lead to longer cut, which may worsen the rendering performance [6] and increase the texturing artifacts [1]. On the other hand, fewer distortion points usually produce higher isometric distortion. Therefore, there is a tradeoff between properties (2) and (3).
[width=0.90]highdistbound
[width=0.95]curdist
From manual picking to automatic generation, there have been many attempts to detect distortion points [8, 7, 10]. Since isometric distortion is the key to our problem, a possible strategy for detecting distortion points is to directly use isometric distortion as a core measure. For instance, some methods iteratively parameterize the mesh to the plane and add the point of greatest isometric distortion into the distortion point set [7, 8]. However, the point of greatest distortion is often located at a boundary, thereby leading to early termination and causing high isometric distortion (see Fig. 2). Instead of using planar parameterizations, Chai [10] cluster high isometric distortion points using an asconformalaspossible (ACAP) spherical parameterization of the input mesh [11]. Although the clustering technique used in this method is very effective for detecting distortion points, it remains a challenge to efficiently compute bijective spherical parameterizations with theoretical guarantees.
In this paper, we propose a novel method to detect distortion points. Due to recent developments in planar parameterization techniques [12, 13], a flipfree planar parameterization is more computationally efficient and robust than a spherical parameterization. Thus, we hope to use planar parameterizations instead of spherical parameterizations to find distortion points efficiently and robustly. However, it is difficult due to the following two reasons. First, planar parameterizations requires that the input closed mesh is cut to a disk topology. However, since different cuts result in different parameterizations, the cut location has a significant effect on the detection result. Hence how to prepare a cut and handle large distortion near boundaries is very challenging. Second, robustly identifying distortion points based on planar parameterizations is nontrivial.
In order to eliminate the effects of different cuts, we present a simple and effective voting strategy to detect distortion points. To generate candidate distortion points to vote on, we repeatedly perform the following threestep procedure: first cut the input mesh randomly, then compute a low conformal distortion planar parameterization, and finally detect a set of vertices as candidates. To consistently seek candidates with strong practical robustness, we customize an effective clustering strategy, shown in [10], which uses a parameter to control the balance between the number of distortion points and the final parameterization distortion. Note that the generated candidates may include some vertices near boundaries, but they are usually detected only once since the random cut paths have almost no intersections. Therefore, in the voting phase, we only recognize vertices that appear more than once as the final candidates.
Our approach is simple and performs better than current stateoftheart methods. We conduct various experiments and comparisons that demonstrate the feasibility and effectiveness of our method. Furthermore, we successfully employ the detected distortion points in three applications, including lowdistortion parameterizations, semiautomatic corresponding landmark detection for bijective surface mappings, and isotropic remeshing (Fig. 1).
2 Related Work
Distortion point detection. Since vertices with high curvatures often induce high isometric distortion, some algorithms try to detect distortion points by using the curvature information as a guide [15, 9]. However, as curvaturebased methods do not consider distortion directly, some distortion points may be ignored (see Fig. 3). Distortion measures are used to develop an iterative algorithm that alternately flattens the mesh and puts the maximum distortion vertex into the distortion point set [7]. However, as mentioned in the previous section and in Fig. 2, it is difficult to avoid early termination due to the maximum distortion vertex’s location on the boundary. To overcome this issue, spherical parameterizations are used [10]; however, efficiently computing bijective spherical parameterizations for complex objects is a nontrivial task (see the comparison in Fig. 17).
Conformal cone singularities. In the context of conformal flattening, finding cone singularities is an effective way to mitigate isometric distortion. Kharevych [16] manually place cone singularities and adjust cone angles to reduce distortion. Similar to [7], an iterative algorithm is used in [8] to compute cone configurations, including the number, placement, and size. Curvatures are also used to find the cone singularities [17, 18, 19]; however, curvatures do not always provide useful information about how cones should be arranged, as observed in [20]. Thus, starting from a different perspective, Soliman [20] study the problem of finding optimal cones as an approximation problem. However, their optimal cone configurations are affected by a parameter that has no intuitive geometric interpretation. When using default values, some important vertices are missed (see Fig. 15 in [20]). In our method, the tradeoff between the number of distortion points and the level of isometric distortion is controlled by a parameter that represents the range of influence of distortion points.
Singularities of fields. Some recent quadrilateral remeshing methods [4, 21] cut the input mesh to a disk topology by adding the singularities of the designed fields into the cut paths. The field design procedure [22] and parameterization computation step are separate, so the connections between the field singularities and parameterization distortion are not close and indirect. As shown in Fig. 17 in [20], the highly regular fields generated by [23] do not immediately result in a low isometric distortion parameterization. Thus, the isometric distortion of parameterizations may still be large, even if the singularities are from highly smooth fields.
Cut construction and atlas generation. Instead of defining cuts by connecting vertices, some methods directly construct cut paths to generate low isometric distortion parameterizations. Mesh segmentation approaches [24, 25, 26, 27, 28, 29] partition an input mesh into multiple charts, and each chart is parameterized with very low isometric distortion. Generating a quad layout [30] is an another way to segment an input mesh into multiple charts, such as the domain simplicity [31], aligned global parameterizations [32], perfect matching [33, 34], reliable quad meshing [35], and filed tracing [36, 37]. Techniques [38, 1] that simultaneously cut and flatten input meshes have also been developed. Li [39] propose an algorithm to jointly optimize the parameterization and the cutpreserving bijection of the mapping. In contrast to these methods, our goal is to automatically detect distortion points, and use them in applications beyond parameterization generation. Liu [40] generate cuts for peeling art.
3 Voting for distortion points
3.1 Overview
We study distortion point detection on 3D surface models. Given a triangular mesh containing vertices and triangles , our goal is to find a set of distortion points . To achieve this goal, we propose a votingbased algorithm. In Section 3.2, we present a threestep method to generate candidate distortion points. After performing the threestep procedure many times, we develop a voting strategy to determine the resulting distortion points, as described in Section 3.3. Some implementation details and discussions are provided in Section 3.4. Fig. 4 illustrates the pipeline of our method using a bear model.
3.2 Candidate distortion points
Instead of spherical parameterizations [10], we use distortion measures on planar parameterizations to detect distortion points. Thus, we first cut as a disk topology mesh (Section 3.2.1), then parameterize onto the plane (Section 3.2.2), and finally find the distortion points (Section 3.2.3).
3.2.1 Constructing cuts
Genuszero surfaces. We generate cuts for genuszero surfaces as follows:

Randomly select a vertex ;

Find the farthest vertex from , i.e., satisfies , for all , where refers to the Euclidean norm;

Compute the shortest path between and along the mesh edges as a cut.
Highgenus surfaces. For highgenus surfaces, we first convert them into genuszero surfaces, and then find the cut using the above procedure for genuszero surfaces. The conversion from high genus to zero genus is as follows:

Compute the handles of the highgenus surface using [41];

Randomly perturb the generated handles by adding a random offset in terms of the geodesic distance;

Cut the highgenus surface along these perturbed handles to create some holes;

Fill in the holes to generate a genuszero surface.
Fig. 5 uses an example to show the distortion detection process for highgenus models. In general, vertices that are located on handles are usually identified as distortion points using our detection algorithm. To distinguish them from true distortion points, we randomly perturb the generated handles of [41] so that the vertices on the handles become different. After the voting process, the vertices caused by the handles are filtered out since they only appear once.
3.2.2 Parameterization of
Inspired by methods [7, 8] that iteratively compute a conformal parameterization and place a distortion point at the vertex with the highest isometric distortion, our method also generates a parameterization that is as conformal as possible (ACAP) with no flips for detecting distortion points. Numerous techniques for ACAP parameterizations have been developed (cf. the surveys in [42, 43, 44]). To achieve a flipfree and ACAP conformal parameterization, we initialize the parameterization using Tutte’s embedding method and optimize the AMIPS energy [45] to reduce the conformal distortion.
Formulation. A parameterization of is constituted by piecewise constant affine maps defined on triangles . After defining a coordinate system on each triangle , each affine map has a linear form , where refers to the Jacobian of . The MIPS energy [46] is defined as follows:
(1) 
where and
are singular values of
, and is the Frobenius norm. reaches its minimum when represents a similarity transformation. The AMIPS energy [45] is defined as follows:(2)  
Optimizing penalizes the maximum distortion and generates evenly and smoothly distributed distortion.
Optimization. Unfortunately, conformal parameterizations are not unique. For example, a new parameterization , where is an arbitrary Möbius transformation, contains the same conformal distortion as . However, the area distortion of and are likely to be very different. Thus, if the generated ACAP parameterizations are not appropriate, the vertices with high isometric distortion may be located incorrectly. These conditions would affect our technique’s ability to robustly detect distortion points.
In our method, is initialized by Tutte’s embedding method [12]. As observed in our experiments, if we try to minimize the MIPS and AMIPS energy while fixing the boundary vertices, the vertices with high isometric distortion often become our desired distortion points. We also test free boundary cases and observe that the distortion distribution is similar except for the distortion near boundaries. (see the comparison in Fig. 6). However, the running time of free boundary parameterization is faster. Thus, the boundary vertices are always free during the conformal distortion optimization in our experiments, which use the KPNewton method [47] to optimize .
3.2.3 Detecting distortion points
Isometric distortion. After achieving an ACAP parameterization , we select the vertices with high isometric distortion as our distortion points. The isometric distortion [45] is defined as follows:
(3)  
If is a rotation matrix, reaches its minimum.
Clustering method. Since the clustering method proposed in [10] is practically robust for detecting distortion points , we first review it and then modify it for our ACAP planar parameterizations.
In [10], the distortion triangles each have local maximum distortion over a region. These triangles are first detected, and then one vertex from each distortion triangle is added to . The pseudocode of the detection algorithm for is listed in Alg. 1. The threshold indicates the distortion point scale (see a more detailed discussion in Section 3.4).
In the first iteration, since contains all of the triangles, the algorithm adds the triangle with maximum distortion into and filters out half of the triangles. For some models with ACAP planar parameterizations, some distortion triangles may be filtered out by the median (see Fig. 7 left). To avoid this issue, we modify the clustering method by choosing a different strategy for the first iteration. Since our goal is to detect the triangles at the local maxima, the isometric distortion in our target triangles is at least greater than the small value . In the first iteration, we filter out the triangles with isometric distortions that are smaller than . Then, we group the remaining triangles into new isolated regions to initialize the queue . Finally, we run the “while” loop of Alg. 1 to compute using the median filter. With this simple modification, our method catches obvious distortion triangles, resulting in parameterizations with lower isometric distortion (see Fig. 7 right). In our experiments, we set in all of the models.
[width=0.99]bird3
3.3 Voting
We present a simple voting strategy to generate final distortion points. This strategy is based on the fact that no matter how cuts change, some vertices are almost always detected as distortion points. The voting strategy contains two steps: (1) perform the candidate distortion point generation procedure many times (ten in our experiments); and (2) vote for the candidates that appear multiple times (more than twice in our experiments) to serve as the resulting distortion points. Fig. 4 (b) shows the voting results from the candidates shown in Fig. 4 (a1–a4).
Postfiltering. Due to numerical issues, some distortion points are very close to each other. In fact, grouping these close distortion points into one point does not affect the resulting isometric distortion (see the comparison in Fig. 8). If one distortion point is within the ring neighbors of another distortion point , we keep the one that has more votes and discard the other one. In the experiments, we observe that the very close points always lie in onering or tworing neighborhoods, so we set as 5 to filter out close points effectively.
3.4 Implementation details
Mesh resolution. Since we compute ACAP parametrizations many times, this process is very timeconsuming for largescale models. In fact, if a lowresolution mesh approximates the highresolution model well and and are cut to disk topologies by similar cut paths, the distortion distributions of their planar parameterizations are similar as well. As a result, the positions of the distortion points detected by are similar to the positions detected by . Fig. 9 shows an example. In our experiments, if , our detection algorithm contains three steps:

Simplify using QEM [48] to contain vertices;

Detect the distortion points of the simplified mesh;

For each detected distortion point, find the nearest vertex on , which becomes the resulting distortion point for .
We set in all examples. If the threshold is too small, some features may be missed. Fig. 11left shows an example of . Note that the top of the octopus is missing, leading to high distortion.
We consider the simplification as an optional operation that aims to accelerate the computation. In our experiments, the simplification is not applied everywhere. Specifically, if the mesh has a small number of vertices, there is no need to simplify the mesh. Besides, if we do not simplify largescale meshes, our method is still able to detect distortion points but takes more time. Our experiments show that the simplification does not affect the robustness in practice.
[width=0.99]simplification
[width=0.99]diffn
Choice of . The parameter controls the balance between the number of distortion points and the isometric distortion. A large value indicates that a distortion point has a large range of influence, so the number of distortion points is often small, resulting in high isometric distortion. A small value usually leads to a large number of distortion points and small isometric distortion. By default, we set if ; otherwise, . Fig. 10 shows a comparison between different values for .
[width=0.95]parameters
Parameter selection. There are some parameters in the voting strategy. We run ten candidate generation processes to generate candidates and select the candidate points that get at least three votes as distortion points. The reasons for these settings are based on three cases of candidate points. First, the candidate distortion points that are near the random cut have at most two votes (e.g. the points near a cut in Fig. 4 (a2) top). These points are caused by the cut and should not be involved in the final distortion points. Second, when some candidate distortion points are located very close together (as shown in the red box in Fig. 8), they are not detected at the same time in one candidate generation process. We observe that there are usually two or three candidate distortion points that are closely located in our experiments. This case should be distinguished from the first case. Third, the isolated candidate distortion points get almost all of the votes. Therefore, in ten candidate generation processes, the candidate distortion points typically get one or two votes for the first case, four or five votes for the second case, and 8 to ten votes for the third case. On one hand, if there are fewer candidate generation processes, the points in the first and second situations cannot be distinguished. On the other hand, it is feasible to conduct more candidate generation processes, but the entire algorithm requires more time to yield the final result. We show a result with five candidates in the middle of Fig. 11. In this example, the algorithm does not find the point at the top of the model, and thus producing a highdistortion result.
[width=0.90]persistence
[width=0.95]tessellation
4 Experiments and evaluations
We apply our distortion point detection method to various models. In Section 4.1, we first evaluate the results and components of our algorithm and then demonstrate the practical robustness of our algorithm on different tessellations and complex models. We select the Seamster [9], Geometry Images method [7], and a spherebased clustering method [10] as the competitors (Section 4.2). The experiments were performed on a desktop PC with an Intel Core i74790 processor and 8GB memory.
Qualitative and quantitative evaluation methods. Since the main purpose of finding distortion points is to decrease the distortion in parameterizations, a simple quantitative evaluation method is to conduct parameterizations and then compute the resulting distortion distribution. Specifically, we adopt the isometric distortion metric defined in [45]
to evaluate the quality of the parameterizations. Then we report the distortion distribution by offering the maximum (worst case), average, and standard deviations for all triangles, denoted as
, , and , respectively. Meanwhile, we visualize the distortion distribution by shading the triangles according to the isometric distortion, with white being optimal (See Fig. 4 (b) top). In addition, another popular visualization method we adopt is rendering mappings using a checkerboard texture with letters (See Fig. 4 (b) bottom).4.1 Evaluations
Other clustering methods. There are other existing clustering methods that detect points according to the distortion of a parameterization. We select one typical persistencebased method [49] to compare with the effectiveness of our method. As shown in Fig. 12, the distortion points detected by the persistencebased clustering method are inappropriate, i.e., some points are superfluous, but some important regions are not detected. In fact, the persistencebased method is regarded as a topdown approach, which finds the most significant features according to the topology persistence of the distortion. However, our clustering method is a bottomup approach, which filters out low distortion areas hierarchically. As seen in the comparison, our bottomup approach is more suitable for finding points that cause high distortion.
Tessellations. In Fig. 13, distortion point detection is performed for six types of tessellations that represent the same shape. The distortion points detected by our method are in very similar places, and the isometric distortions of the parameterizations are all at a low level. Since the final distortion points are obtained by voting from several parameterizations, our method exhibits reliable results when using different tessellations as inputs.
Practical robustness. Our method successfully detects the desired distortion points, resulting in parameterizations with low isometric distortion in all of the models. We show eleven complex models in Fig. 18, including five genuszero models and six highgenus models. The C++ implementation for detecting distortion points for genuszero surfaces is provided in the supplementary material.
Timings. The running time for distortion point detection on an ant model with 10k vertices is 8.64 seconds. For a dragon model with 50k vertices, the mesh simplification and the distortion point detection take 3.44 seconds and 14.06 seconds, respectively. Thanks to the use of mesh simplification, our method takes about 15 seconds on average. The more structurally complex the input model is, the longer our method will take, since more iterations are required to compute ACAP parameterizations.
[width=0.90]seamster
[width=0.90]geoimg
[width=0.99]hcc
4.2 Comparisons
4.2.1 Comparison to Seamster [9]
The Seamster method uses curvatures as guides for detecting distortion points, so the vertices with low curvatures are not recogized. However, some low curvature vertices are important for reducing isometric distortion. In Fig. 3 and Fig. 15, we compare our method to the Seamster method. Due to some missing distortion points with low curvatures, the parameterizations generated by the Seamster method often have greater distortion than ours.
4.2.2 Comparison to Geometry Images [7]
In Fig. 2 and Fig. 16, comparisons between our method and the Geometry Images method are conducted. Due to early termination issues, some important points are missed, e.g., the points at the top of the human head (Fig. 2) and the gorilla head (Fig. 16). Thus, the parameterization distortion caused by the Geometry Images method is usually higher than ours.
4.2.3 Comparison to the spherebased method [10]
In Fig. 17, we compare our method with the spherebased method [10]. Note that since we use a different method to compute the final parameterizations, the distortion values are different from those in [10]. The AQP method [51] and the SA method [52] are used in [10], while the progressive parameterization method [14] is used in this paper. Thus, we recompute a parameterization using the same method and measure the distortion for fair comparison. The detected distortion points are similar, so the parameterization distortion is also similar. However, as the bijective spherical parameterization takes 42.09 seconds, the spherebased method takes more time than our method.
4.2.4 Extensive comparisons
Data set construction. To verify the robustness of our method and compare it with other methods, we construct a data set containing 20,000 random genuszero meshes. The meshes in the data set are constructed from a primitive shape, such as a sphere, a cube, a pyramid, a dodecahedron, etc. Then, we conduct random shape manipulations, including rotation, scaling, random affine transformation, stretching, bending, twisting, random vertex perturbation, fractal displacement [53], CSG manipulations (merge, intersect, and subtract with other models), remeshing, smoothing, and simplification [54]. After each manipulation, if the mesh is not a genus zero manifold, we discard this manipulation and roll back to the previous shape. Each mesh is finally simplified or subdivided to 13,000 vertices. We show five examples in Fig. 18.
Results. We have run the Geometry Images method [7], the Seamster method [9], the spherebased method [10], and our method on the constructed data set. The Seamster method, the Geometry Images method, and our method succeed in all models, while the spherebased method fails in 508 models. The spherebased method depends on a spherical parameterization. If the input model has a long and thin neck, the spherebased method [10] usually fails due to its spherical parameterization method [11]. Specifically, such a long and thin neck provides no enough space to insert vertices, so the spherical parameterization gets stuck in the vertex insertion phase. Since the possible vertex positions are enumerated discretely by sampling points in a triangle, if flips occur for all vertex positions, the algorithm terminates and does not generate a flipfree spherical parameterization. Fig. 18 shows its five failed examples.
5 Applications
We employ our detected distortion points in three applications.
[width=0.99]autocut
[width=0.99]optcuts
5.1 Parameterizations
Computing parameterizations is a straightforward application of our method. Given a closed mesh, we generate a cut by simply constructing a minimal spanning tree along the mesh edges that passes through all the distortion points detected by our method. After cutting the mesh into a disk topology, any of the recent parameterization methods are able to generate a lowdistortion UV map, which is used for texture mapping. In our experiments, we use the progressive parameterizations method [14] to generate UV maps and color the triangles according to the isometric distortion with red indicating high distortion. In addition, we also use a checkerboard texture with letters to show the levels of distortion.
In Fig. 20 and Fig. 21, we give a brief comparison with the recent Autocuts method [1] and OptCuts method [39]. We only adjust the parameters and in the Autocuts method to obtain results with similar distortion levels to ours. For the OptCuts method, we turn off the bijective option and also adjust the parameters to obtain similar distortion results. Then, we compare the cut lengths with these two methods. The cut lengths of our result are shorter than theirs. Moreover, our method is fully automatic and does not require any user intervention.
5.2 Semiautomatic landmark correspondence
Given two meshes, we first generate the distortion points for each mesh and use these points as corresponding landmarks, where the correspondence is manually specified. Then, we find cuts that pass through these landmarks and use the lifted bijection method [5] to compute a bijective surface mapping. We give an example of this application in Fig. 1 center, and a comparison in Fig. 22. From the comparison, we find that the locations of distortion points (center) and manually selected landmarks (top) are similar, and the distortions of the resulting intersurface mappings are also very similar. However, if some of the distortion points are missing (bottom), the distortion becomes dramatically larger. Therefore, our method has the ability to automatically detect the most important features in models and reduce user interactions when generating intersurface mappings.
[width=0.95]correspondence
5.3 Isotropic remeshing
The incremental remeshing method [3] is an efficient and effective algorithm for generating isotropic triangle meshes. Given a target edge length, this method repeatedly carries out four steps: (1) splitting long edges, (2) collapsing short edges, (3) equalizing valences using edge flipping, and (4) relocating vertices. However, certain features, such as thin parts and narrow ditches, are often collapsed (see Fig. 23center). When a cut passes through the distortion points, we usually get a parameterization with very low isometric distortion. Since the mapping from the original surface to the parameterization domain is able to preserve isotropy, we do an isotropic remeshing on the parameterization mesh. First, we split and collapse edges on the cut before the parameterization, and fix the boundary when remeshing the 2D parameterization mesh. Then, the four aforementioned steps in the incremental remeshing method are carried out to create a 2D isotropic mesh, which is mapped back onto the original mesh. Fig. 23 shows a comparison of three dog models. The thin features collapse in the incremental remeshing method [3], while our method preserves these features successfully. The reason is that the edge lengths in [3] are computed using the 3D Euclidean distance, which yields short lengths at the thin parts; thus, the collapsing operations destroy these thin features. However, the edge lengths in our method are computed using 2D Euclidean distances that approximate the geodesic distances of an input mesh, which is able to preserves thin features. This characteristic is very useful in the errorbounded remeshing [56].
6 Conclusion
In this paper, we present a novel approach for detecting distortion points on 3D triangle meshes. This method can be further applied to planar parameterization, semiautomatic landmark correspondence, and isotropic remeshing.
Our method has some limitations. First, there is no theoretical analysis regarding the relationships between the detected distortion points and the resulting parameterization distortion. Actually, our method is empirical, based on the observation that the isometric distortion will be low if the cuts pass through the distortion points. Although cutting through the distortion points tends to generate lowdistortion parameterizations, this condition is not always guaranteed. For example, as shown in Fig. 13, the noisy model has higher distortion than the others, even if the cut passes through the distortion points. The second limitation is that for a sphereshaped mesh, our method cannot find any distortion points. The reason is that the vertices are almost all symmetrically equivalent, and the random cuts generate random candidate distortion points, each of which has only one vote.
References
 [1] R. Poranne, M. Tarini, S. Huber, D. Panozzo, and O. SorkineHornung, “Autocuts: Simultaneous Distortion and Cut Optimization for UV Mapping,” ACM Trans. Graph., vol. 36, no. 6, pp. 215:1–215:11, 2017.
 [2] H.Y. Liu, X.M. Fu, C. Ye, S. Chai, and L. Liu, “Atlas refinement with bounded packing efficiency,” ACM Trans. Graph., vol. 38, no. 4, pp. 33:1–33:13, 2019.
 [3] M. Botsch and L. Kobbelt, “A remeshing approach to multiresolution modeling,” in Proceedings of the 2004 Eurographics/ACM SIGGRAPH Symposium on Geometry processing, 2004, pp. 185–192.
 [4] D. Bommes, H. Zimmer, and L. Kobbelt, “Mixedinteger quadrangulation,” ACM Trans. Graph., vol. 28, no. 3, pp. 77:1–77:10, 2009.
 [5] N. Aigerman, R. Poranne, and Y. Lipman, “Lifted bijections for low distortion surface mappings,” ACM Trans. Graph., vol. 33, no. 4, pp. 69:1–69:12, 2014.
 [6] Z. S. Hakura and A. Gupta, “The Design and Analysis of a Cache Architecture for Texture Mapping,” SIGARCH Comput. Archit. News, vol. 25, no. 2, pp. 108–120, 1997.
 [7] X. Gu, S. J. Gortler, and H. Hoppe, “Geometry Images,” ACM Trans. Graph., vol. 21, no. 3, pp. 355–361, 2002.
 [8] B. Springborn, P. Schröder, and U. Pinkall, “Conformal equivalence of triangle meshes,” ACM Trans. Graph., vol. 27, no. 3, pp. 77:1–77:11, 2008.
 [9] A. Sheffer and J. C. Hart, “Seamster: inconspicuous lowdistortion texture seam layout,” in Proceedings of the conference on Visualization’02, 2002, pp. 291–298.
 [10] S. Chai, X.M. Fu, X. Hu, Y. Yang, and L. Liu, “Spherebased Cut Construction for Planar Parameterizations,” Computer & Graphics (SMI 2018), vol. 74, pp. 66–75, 2018.
 [11] X. Hu, X.M. Fu, and L. Liu, “Advanced Hierarchical Spherical Parameterizations,” IEEE. T. Vis. Comput. Gr., vol. 24, no. 6, pp. 1930–1941, 2018.
 [12] W. T. Tutte, “How to draw a graph,” in Proceedings of the London Mathematical Society, vol. 13, 1963, pp. 747–767.
 [13] H. Shen, Z. Jiang, D. Zorin, and D. Panozzo, “Progressive Embedding,” ACM Trans. Graph., vol. 38, no. 4, pp. 32:1–32:13, 2019.
 [14] L. Liu, C. Ye, R. Ni, and X.M. Fu, “Progressive Parameterizations,” ACM Trans. Graph., vol. 37, no. 4, pp. 41:1–41:12, 2018.
 [15] A. Sheffer, “Spanning tree seams for reducing parameterization distortion of triangulated surfaces,” in Shape Modeling International, 2002, pp. 61–66.
 [16] L. Kharevych, B. Springborn, and P. Schröder, “Discrete conformal mappings via circle patterns,” ACM Trans. Graph., vol. 25, no. 2, pp. 412–438, 2006.
 [17] M. BenChen, C. Gotsman, and G. Bunin, “Conformal flattening by curvature prescription and metric scaling,” in Comput. Graph. Forum, vol. 27, no. 2, 2008, pp. 449–458.
 [18] A. Myles and D. Zorin, “Global Parametrization by Incremental Flattening,” ACM Trans. Graph., vol. 31, no. 4, pp. 109:1–109:11, 2012.
 [19] A.M. Vintescu, F. Dupont, G. Lavoué, P. Memari, and J. Tierny, “Conformal factor persistence for fast hierarchical cone extraction,” in Eurographics 2017  Short Papers, 2017, pp. 57–60.
 [20] Y. Soliman, D. Slepčev, and K. Crane, “Optimal Cone Singularities for Conformal Flattening,” ACM Trans. Graph., vol. 37, no. 4, pp. 105:1–105:17, 2018.
 [21] D. Bommes, B. Lévy, N. Pietroni, E. Puppo, C. Silva, M. Tarini, and D. Zorin, “Quadmesh generation and processing: a survey,” Comput. Graph. Forum, vol. 32, no. 6, pp. 51–76, 2013.
 [22] A. Vaxman, M. Campen, O. Diamanti, D. Panozzo, D. Bommes, K. Hildebrandt, and M. BenChen, “Directional field synthesis, design, and processing,” Comput. Graph. Forum, vol. 35, no. 2, pp. 545–572, 2016.
 [23] F. Knöppel, K. Crane, U. Pinkall, and P. Schröder, “Globally optimal direction fields,” ACM Trans. Graph., vol. 32, no. 4, pp. 59:1–59:10, 2013.
 [24] B. Lévy, S. Petitjean, N. Ray, and J. Maillot, “Least squares conformal maps for automatic texture atlas generation,” ACM Trans. Graph., vol. 21, no. 3, pp. 362–371, 2002.
 [25] P. V. Sander, S. J. Gortler, J. Snyder, and H. Hoppe, “Signalspecialized Parametrization,” in Proceedings of the 13th Eurographics Workshop on Rendering, 2002, pp. 87–98.
 [26] P. V. Sander, Z. J. Wood, S. J. Gortler, J. Snyder, and H. Hoppe, “Multichart Geometry Images,” in Proceedings of the 2003 Eurographics/ACM SIGGRAPH Symposium on Geometry Processing, 2003, pp. 146–155.
 [27] K. Zhou, J. Synder, B. Guo, and H.Y. Shum, “Isocharts: Stretchdriven Mesh Parameterization Using Spectral Analysis,” in Proceedings of the 2004 Eurographics/ACM SIGGRAPH Symposium on Geometry Processing, 2004, pp. 45–54.
 [28] D. Julius, V. Kraevoy, and A. Sheffer, “DCharts: QuasiDevelopable Mesh Segmentation,” in Comput. Graph. Forum, vol. 24, no. 3, 2005, pp. 581–590.
 [29] E. Zhang, K. Mischaikow, and G. Turk, “Featurebased surface parameterization and texture mapping,” ACM Trans. Graph., vol. 24, no. 1, pp. 1–27, 2005.
 [30] M. Campen, “Partitioning surfaces into quadrilateral patches: a survey,” Comput. Graph. Forum, vol. 36, no. 8, pp. 567–588, 2017.
 [31] M. Tarini, E. Puppo, D. Panozzo, N. Pietroni, and P. Cignoni, “Simple Quad Domains for Field Aligned Mesh Parametrization,” ACM Trans. Graph., vol. 30, no. 6, pp. 142:1–142:12, 2011.
 [32] M. Campen and L. Kobbelt, “Quad Layout Embedding via Aligned Parameterization,” Comput. Graph. Forum, vol. 33, no. 8, pp. 69–81, 2014.
 [33] F. H. Razaflndrazaka, U. Reitebuch, and K. Polthier, “Perfect Matching Quad Layouts for Manifold Meshes,” Comput. Graph. Forum, vol. 34, no. 5, pp. 219–228, 2015.
 [34] F. H. Razafindrazaka and K. Polthier, “Optimal base complexes for quadrilateral meshes,” Comput. Aided Geom. Des., vol. 5253, pp. 63 – 74, 2017.
 [35] D. Bommes, M. Campen, H.C. Ebke, P. Alliez, and L. Kobbelt, “Integergrid maps for reliable quad meshing,” ACM Trans. Graph., vol. 32, no. 4, pp. 98:1–98:12, 2013.
 [36] M. Campen, D. Bommes, and L. Kobbelt, “Dual loops meshing: quality quad layouts on manifolds,” ACM Trans. Graph., vol. 31, no. 4, pp. 110:1–110:11, 2012.
 [37] N. Pietroni, E. Puppo, G. Marcias, R. Scopigno, and P. Cignoni, “Tracing fieldcoherent quad layouts,” Comput. Graph. Forum, vol. 35, no. 7, pp. 485–496, 2016.
 [38] O. Sorkine, D. CohenOr, R. Goldenthal, and D. Lischinski, “Boundeddistortion piecewise mesh parameterization,” in Proceedings of the Conference on Visualization ’02, 2002, pp. 355–362.
 [39] M. Li, D. M. Kaufman, V. G. Kim, J. Solomon, and A. Sheffer, “OptCuts: Joint Optimization of Surface Cuts and Parameterization,” ACM Trans. Graph., vol. 37, no. 6, pp. 247:1–247:13, 2018.
 [40] H. Liu, X.T. Zhang, X.M. Fu, Z.C. Dong, and L. Liu, “Computational peeling art design,” ACM Trans. Graph., vol. 38, no. 4, pp. 64:1–64:12, 2019.
 [41] T. K. Dey, F. Fan, and Y. Wang, “An Efficient Computation of Handle and Tunnel Loops via Reeb Graphs,” ACM Trans. Graph., vol. 32, no. 4, pp. 32:1–32:10, 2013.
 [42] M. S. Floater and K. Hormann, “Surface parameterization: a tutorial and survey,” in In Advances in Multiresolution for Geometric Modelling. Springer, 2005, pp. 157–186.
 [43] A. Sheffer, E. Praun, and K. Rose, “Mesh parameterization methods and their applications,” Found. Trends. Comput. Graph. Vis., vol. 2, no. 2, pp. 105–171, 2006.
 [44] K. Hormann, B. Lévy, and A. Sheffer, “Mesh Parameterization: Theory and Practice,” in ACM SIGGRAPH 2007 Courses, ser. SIGGRAPH ’07, 2007.
 [45] X.M. Fu, Y. Liu, and B. Guo, “Computing locally injective mappings by advanced MIPS,” ACM Trans. Graph., vol. 34, no. 4, pp. 71:1–71:12, 2015.
 [46] K. Hormann and G. Greiner, “MIPS: An efficient global parametrization method,” in Curve and Surface Design: SaintMalo 1999. Vanderbilt University Press, 2000, pp. 153–162.
 [47] B. Golla, H.P. Seidel, and R. Chen, “Piecewise linear mapping optimization based on the complex view,” Comput. Graph. Forum, vol. 37, no. 7, pp. 233–243, 2018.
 [48] M. Garland and P. S. Heckbert, “Surface Simplification Using Quadric Error Metrics,” in Proceedings of the 24th Annual Conference on Computer Graphics and Interactive Techniques, ser. SIGGRAPH ’97, 1997, pp. 209–216.
 [49] F. Chazal, L. J. Guibas, S. Y. Oudot, and P. Skraba, “Persistencebased clustering in riemannian manifolds,” Journal of the ACM (JACM), vol. 60, no. 6, p. 41, 2013.
 [50] X.M. Fu, Y. Liu, J. Snyder, and B. Guo, “Anisotropic simplicial meshing using local convex functions,” ACM Trans. Graph., vol. 33, no. 6, pp. 182:1–182:11, 2014.
 [51] S. Z. Kovalsky, M. Galun, and Y. Lipman, “Accelerated Quadratic Proxy for Geometric Optimization,” ACM Trans. Graph., vol. 35, no. 4, pp. 134:1–134:11, 2016.
 [52] X.M. Fu and Y. Liu, “Computing InversionFree Mappings by Simplex Assembly,” ACM Trans. Graph., vol. 35, no. 6, 2016.
 [53] D. S. Ebert and F. K. Musgrave, Texturing & modeling: a procedural approach. Morgan Kaufmann, 2003.
 [54] M. Botsch, L. Kobbelt, M. Pauly, P. Alliez, and B. Lévy, Polygon mesh processing. AK Peters/CRC Press, 2010.
 [55] Y. Yang, X.M. Fu, S. Chai, S.W. Xiao, and L. Liu, “Volumeenhanced compatible remeshing of 3d models,” IEEE. T. Vis. Comput. Gr., vol. 25, no. 10, pp. 2999–3010, 2019.
 [56] X.X. Cheng, X.M. Fu, C. Zhang, and S. Chai, “Practical errorbounded remeshing by adaptive refinement,” Computers & Graphics, vol. 82, pp. 163–173, 2019.