Ellipsoidal Density-Equalizing Map for Genus-0 Closed Surfaces
Abstract
Surface parameterization is a fundamental task in geometry processing and plays an important role in many science and engineering applications. In recent years, the density-equalizing map, a shape deformation technique based on the physical principle of density diffusion, has been utilized for the parameterization of simply connected and multiply connected open surfaces. More recently, a spherical density-equalizing mapping method has been developed for the parameterization of genus-0 closed surfaces. However, for genus-0 closed surfaces with extreme geometry, using a spherical domain for the parameterization may induce large geometric distortion. In this work, we develop a novel method for computing density-equalizing maps of genus-0 closed surfaces onto an ellipsoidal domain. This allows us to achieve ellipsoidal area-preserving parameterizations and ellipsoidal parameterizations with controlled area change. We further propose an energy minimization approach that combines density-equalizing maps and quasi-conformal maps, which allows us to produce ellipsoidal density-equalizing quasi-conformal maps for achieving a balance between density-equalization and quasi-conformality. Using our proposed methods, we can significantly improve the performance of surface remeshing for genus-0 closed surfaces. Experimental results on a large variety of genus-0 closed surfaces are presented to demonstrate the effectiveness of our proposed methods.
1 Introduction
Surface parameterization is the process of mapping a complicated surface onto a simpler domain. It is closely related to many tasks in geometry processing including surface registration, texture mapping, surface remeshing, and shape analysis. In recent decades, numerous efforts have been devoted to the development of efficient surface parameterization algorithms with different desired properties [1, 2, 3, 4] for various science and engineering applications.
As genus-0 closed surfaces are topologically equivalent to a sphere, most prior works have focused on the problem of parameterizing genus-0 closed surfaces onto the unit sphere. For conformal parameterization, a vast number of spherical conformal parameterization methods have been developed over the past few decades based on linearization [5], harmonic energy minimization [6], Ricci flow [7], conformal curvature flow [8], quasi-conformal theory [9], partial welding [10] and Dirichlet energy minimization [11]. For area-preserving parameterization, the theory and computation of optimal mass transportation (OMT) have been extensively studied [12, 13]. Based on OMT, some area-preserving parameterization methods for genus-0 closed surfaces have been developed [14, 15, 16, 17]. Besides the above-mentioned conformal and area-preserving methods, some other prior works have also developed spherical parameterization algorithms with different mapping criteria [18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28]. While there is a topological equivalence between genus-0 closed surfaces and the sphere, the geometric difference between them may be large. For instance, mapping a genus-0 closed surface with an elongated shape to a sphere may induce large geometric distortion. Therefore, a few recent studies have focused on ellipsoidal parameterizations [29, 30, 31], in which an ellipsoid is used as the parameter domain for parameterizing genus-0 closed surfaces.
Density-equalizing maps [32] are a class of mappings with shape deformation produced based on prescribed density distribution. More specifically, given a planar domain with a density distribution prescribed in it, a density-equalizing map will deform the domain with the high-density regions enlarged and the low-density regions shrunk. Over the past few decades, density-equalizing maps have been widely used in cartogram creation and data visualization [33]. Also, different variants and improved algorithms have been developed [34, 35, 36]. Choi and Rycroft [37] developed a novel method for computing surface density-equalizing maps and demonstrated its effectiveness in surface parameterization for simply connected open surfaces. Later, the method was extended for other mapping problems for simply connected open surfaces [38, 39, 40] and volumetric domains [41] with applications to medical visualization and shape modeling. Lyu et al. [42] developed a surface parameterization method for multiply connected open surfaces by combining density-equalizing maps and quasi-conformal maps. More recently, they have developed spherical density-equalizing mapping methods for genus-0 closed surfaces [43]. Using their methods, spherical area-preserving parameterizations and spherical parameterizations with controlled area change can be achieved. However, as mentioned above, mapping genus-0 closed surfaces with extreme geometry onto a spherical domain may induce large geometric distortion. For instance, the spherical area-preserving parameterization of an elongated surface may possess low area distortion but extremely high angular distortion, which is undesirable for many practical applications.
Motivated by the above works, here we develop a novel method for computing ellipsoidal density-equalizing maps for genus-0 closed surfaces in . Specifically, given a genus-0 closed surface, we compute density-equalizing maps onto ellipsoidal domains with different prescribed radii based on a density distribution encoding the desired mapping effect (Fig. 1(a)). This allows us to easily achieve ellipsoidal area-preserving parameterizations and ellipsoidal parameterizations with controlled area change. We then further propose an algorithm for computing ellipsoidal density-equalizing quasi-conformal maps, which combines density-equalizing maps and quasi-conformal maps to achieve a balance between area and angle distortions. Moreover, throughout the algorithm, we can optimize the shape of the target ellipsoidal domain to further reduce overall geometric distortion (Fig. 1(b)). This provides us with an effective and automatic way of representing any genus-0 closed surface by an optimal ellipsoid. We apply our proposed methods for surface remeshing of genus-0 closed surfaces and demonstrate the improvement over prior approaches via various examples.
The organization of this paper is as follows. In Section 2, we introduce the mathematical background of our work. In Section 3, we describe our proposed methods for computing ellipsoidal density-equalizing maps (EDEM) and ellipsoidal density-equalizing quasi-conformal maps (EDEQ). In Section 4, we present experimental results of our proposed methods on various genus-0 closed surfaces. In Section 5, we describe the application of our proposed methods to genus-0 surface remeshing. Finally, in Section 6, we conclude this paper and discuss possible future works.
2 Mathematical background
Our proposed methods are primarily based on the theory of density-equalizing maps and quasi-conformal geometry. The relevant concepts are introduced in this section.
2.1 Density-equalizing maps
Gastner and Newmann [32] proposed a method for computing the density-equalizing maps for 2D planar domains based on the principle of density diffusion. Given a positive density function defined on the planar domain, the method produces shape deformation following the density gradient. More specifically, note that the advection equation is given by
(1) |
where is the density flux by Fick’s law. This gives the diffusion equation
(2) |
Now, since the flux can be expressed as , where is the velocity field, we have
(3) |
Therefore, the position of any tracer particle at time can be traced by:
(4) |
where is the initial position. Taking the time , the density will be fully equalized on the entire domain, and the resulting shape deformation is a density-equalizing map. In particular, as the shape deformation is induced by the density diffusion process, it is easy to see that high-density regions will be enlarged throughout the process and low-density regions will be shrunk (see Fig. 2 for an illustration).
The above density-equalizing mapping method has been widely applied to cartogram creation and sociological data visualization. Specifically, the input density is commonly defined as some prescribed quantity (known as the “population”) per unit area, where the “population” can either be the actual population of a certain region on the geographical map, the average income of a region, or any other sociological data to be visualized.
In recent years, the density-equalizing mapping method has been further introduced to the field of surface parameterization. In particular, Choi et al. [37, 38] proposed methods for computing the density-equalizing maps for simply connected open surfaces in . Lyu et al. [42] developed parameterization methods for multiply connected open surfaces based on density-equalizing maps. More recently, Lyu et al. [43] developed a method for computing spherical density-equalizing maps for genus-0 closed surfaces.
2.2 Quasi-conformal theory
It is well-known that conformal maps [44] are angle-preserving. Intuitively, they map infinitesimal circles to infinitesimal circles. Quasi-conformal maps [45, 46] are a generalization of conformal maps taking infinitesimal circles to infinitesimal ellipses with bounded eccentricity. Mathematically, a quasi-conformal map satisfies the Beltrami equation:
(5) |
for some complex-valued function satisfying . Here, is called the Beltrami coefficient, which encodes important information about the conformal distortion of the map . Specifically, is a conformal map if and only if , as Eq. (5) becomes the Cauchy–Riemann equation. Intuitively, around a point , the first order approximation of can be expressed as:
(6) |
The above approximation suggests that maps an infinitesimal circle centered at to an infinitesimal ellipse centered at . Additionally, we can determine the angles and scaling factors of the maximal magnification and maximal shrinkage using (see Fig. 3). More specifically, the angle of maximal magnification is given by , with the magnifying factor . Also, the angle of maximal shrinkage is given by , with the shrinking factor . The maximal dilation of is .
Besides encoding important geometric information about the associated quasi-conformal map, the Beltrami coefficient is also related to the bijectivity of the mapping [45, 46]:
2.1.
If is a mapping satisfying , then is bijective.
Moreover, the Beltrami coefficient of a composition of two quasi-conformal maps can be expressed in the following way. Let be two quasi-conformal maps with Beltrami coefficient and , respectively. The Beltrami coefficient of the composition map is given by
(7) |
where . In particular, if is a conformal map, then . In other words, composing a conformal map with a given quasi-conformal map will not change its Beltrami coefficient.
Lui et al. [47] proposed an algorithm called Linear Beltrami Solver (LBS) for efficiently reconstructing a quasi-conformal mapping based on a given Beltrami coefficient . Specifically, for any given , the corresponding quasi-conformal map can be computed by solving the following equations:
(8) |
where with
(9) |
In the discrete case, Eq. (8) becomes two sparse symmetric positive definite linear systems and can be easily solved numerically. In the following, we denote the above method for reconstructing a quasi-conformal map from a Beltrami coefficient by .
3 Proposed methods
In this section, we first propose a novel algorithm for computing ellipsoidal density-equalizing maps (abbreviated as EDEM) of genus-0 closed surfaces onto a prescribed ellipsoid. Next, we propose another algorithm for computing ellipsoidal density-equalizing quasi-conformal maps (abbreviated as EDEQ), which allows us to simultaneously optimize the shape of the ellipsoidal domain and the mapping onto it to achieve minimal geometric distortions.
3.1 Ellipsoidal density-equalizing map (EDEM)
Consider a given genus-0 closed surface in discretized as a triangle mesh , where is the set of vertices, is the set of edges, and is the set of triangle faces. Our goal is to compute a bijective ellipsoidal density-equalizing map , where is a prescribed ellipsoid with elliptic radii :
(10) |
Here, the density involved in the density-equalizing mapping method will be defined by some prescribed positive quantity (denoted as the “population”) per unit area. By changing the “population”, we can easily control the shape deformation and achieve different desired parameterization results.
3.1.1 Initial ellipsoidal parameterization
First, we compute an initial ellipsoidal conformal parameterization .
To accomplish this, we apply the Fast Ellipsoidal Conformal Mapping (FECM) method [30], which is an efficient algorithm for computing an ellipsoidal conformal parameterization by leveraging quasi-conformal theory. Initially, the input genus-0 closed surface is conformally mapped to the unit sphere . Then, the FECM method maps the sphere onto the complex plane using the stereographic projection, followed by a Möbius transformation step that maps two desired polar points of to and , respectively. However, it is important to note that the stereographic projection may lead to an uneven distribution of points. To address this issue, the FECM method incorporates an additional rescaling step to enhance the distribution of points.
Next, the FECM method looks for an inverse projection that maps the planar domain to the target ellipsoid conformally. It is worth noting that while one can extend the ordinally stereographic projection and its inverse to the ellipsoidal case, the inverse ellipsoidal stereographic projection is not inherently conformal. Therefore, by utilizing the idea of quasi-conformal composition, the FECM method identifies a suitable additional quasi-conformal map on the plane and composes it with the inverse ellipsoidal stereographic projection to form a conformal map from the plane to the ellipsoid .
Finally, by combining all of the above-mentioned mappings, the FECM method produces an ellipsoidal conformal parameterization . Readers are referred to [30] for more details.
3.1.2 Density-equalizing map on ellipsoid
After obtaining an initial ellipsoidal parameterization, we handle the subsequent shape deformation problem on the ellipsoidal domain .
We define the initial density as the “population” per unit area, where the “population” is some prescribed positive quantity defined everywhere on the ellipsoid. In the following, we propose an iterative scheme for computing ellipsoidal density-equalizing maps on , with the bijectivity of the mappings enforced. As the density diffusion process is time-dependent, here we denote the position of a vertex on at time as . The initial position of is represented as .
Analogous to the original density-equalizing mapping formulation [32], the diffusion equation of on the ellipsoid is given by:
(11) |
where is the Laplace–Beltrami operator. The velocity field induced by the density gradient on the ellipsoid is given by
(12) |
Then, the position of any vertex on the ellipsoid can be computed using the following equation:
(13) |
As time progresses towards infinity (), the density is equalized and we obtain the desired ellipsoidal density-equalizing map.
Note that in the discrete case, the initial density needs to be discretized on the ellipsoidal triangle mesh. Here, we first define the initial density on every triangular face on the ellipsoid. The subsequent iterative process involves updating the positions of all vertices , where and at different time points , where with the time step size . In addition to the face density , we also denote the vertex density as . Note that can be obtained from by a simple formula [37]:
(14) |
where is a face-to-vertex conversion matrix such that:
(15) |
Here, is the 1-ring neighborhood of the vertex .
Now, to solve Eq. (11) in the discrete case, note that the Laplace-Beltrami operator at the -th iteration can be discretized as:
(16) |
Here, is a diagonal matrix (known as the lumped mass matrix) of size given by:
(17) |
where is a triangle in . It is worth noting that every element of the lumped mass matrix is the area sum of the 1-ring neighborhood of the corresponding vertex at . The matrix is a symmetric sparse matrix with cotangent weights [48]:
(18) |
where and are the two angles opposite to the edge . Now, Eq. (11) can be solved using the semi-discrete backward Euler method:
(19) |
from which we have
(20) |
After solving Eq. (11), the velocity field can be updated using Eq. (12). More specifically, we first compute the discretized density gradient on the triangle mesh. For any triangle element , using the Whitney 0-forms, can be interpolated at any point on the triangle by:
(21) |
where , , are hat functions and , , are the vertex densities at . By the property (see Ref. [37]), the density gradient on the triangle is given by:
(22) |
where , , are the three directed edges , , , and is the outward unit normal vector of the triangle face . Then, using Eq. (23), we can obtain on every vertex by:
(23) |
Therefore, the velocity field is given by:
(24) |
However, it is noteworthy that in the discrete case, the velocity in Eq. (24) may not be located in the tangent space. To resolve this issue, we add a projection step to project the velocity onto the admissible space. Since the equation of the ellipsoid is
(25) |
it is easy to see that the outward unit normal vector of is
(26) |
Therefore, the projected velocity field can be obtained by
(27) |
where is the outward unit normal vector at the vertex . We can then update the vertex positions based on Eq. (13):
(28) |
Note that there may still be small numerical errors that cause the vertices to move outside the ellipsoidal surface. To ensure that all vertices remain on the ellipsoid exactly, we further divide the coordinates of each vertex by , where . This ensures that all updated vertices satisfy the parametric equation in Eq. (10) and hence lie on the ellipsoid .
3.1.3 Ensuring the bijectivity of the ellipsoidal mapping throughout the iterative process
As mentioned in [42, 43], conventional density-equalizing mapping methods do not provide any guarantee of bijectivity throughout the diffusion process. Specifically, mesh overlaps may occur under the mapping update in Eq. (28) if the density gradient or the time step size is too large. To resolve this issue, we follow the idea in [43] and propose a correction scheme that ensures the bijectivity of the ellipsoidal mapping at each iteration.
We first introduce a rescaling transformation defined by:
(29) |
where , , and are the radii of the ellipsoid. Using the rescaling transformation , the initial ellipsoidal map and the ellipsoidal map after the -th iterative update can be mapped onto two unit spheres and , respectively. It is easy to see that is a bijection and hence will not have any effect on the presence or absence of local mesh fold-overs in the ellipsoids.
Now, we can apply the spherical overlap correction scheme proposed by the SDEM work [43], which is designed for enforcing the bijectivity of a spherical mapping. Specifically, the scheme involves projecting the two spheres and onto the extended complex plane via the north-pole stereographic projection and the south-pole stereographic projection and correcting the local mesh fold-overs on the plane using the LBS method [47]. Readers are referred to [43] for the detailed description of the spherical overlap correction scheme. We denote the process of the spherical overlap correction scheme as a map , where the mapping result is folding-free.
Finally, we apply the inverse rescaling transformation to rescale the spherical mapping result to the ellipsoid. Again, it is easy to see that is a bijection and will not have any effect on the presence or absence of mesh fold-overs. Therefore, the composition is a mapping from to that can correct any local mesh fold-overs in . This completes the ellipsoidal overlap correction scheme.
3.1.4 Re-coupling the deformation and density
Recall that the proposed ellipsoidal density-equalizing mapping scheme creates shape deformations on a prescribed ellipsoid based on density diffusion. Specifically, the velocity field at each iteration is determined by the density and its gradient at the mesh vertices. As the computational procedure involves multiple discretization schemes, numerical errors may accumulate throughout the iterations. Also, the above-mentioned overlap correction scheme may alter the vertex positions when correcting the local mesh fold-overs, thereby leading to a discrepancy between the vertex positions and the actual density flow. To address this issue, we follow the idea in [43] and introduce an extra step to re-couple the density and deformation at the end of every iteration.
More specifically, at the end of the -th iteration, we do not directly use the density obtained by solving the diffusion equation in Eq. (20) to continue the next iteration. Instead, we re-define the density on the triangle element using the updated vertex positions on the ellipsoid:
(30) |
Once the updated density is obtained, we calculate using Eq. (14) and proceed with the next iteration. Altogether, this additional re-coupling step allows us to reduce the accumulation of numerical errors as well as the discrepancy between the shape deformation and the density.
3.1.5 Summary
By integrating the initial ellipsoidal conformal parameterization (Section 3.1.1) and the diffusion iteration scheme (Section 3.1.2) with the overlap correction scheme (Section 3.1.3) and the re-coupling scheme (Section 3.1.4), we have the proposed EDEM algorithm for computing bijective ellipsoidal density-equalizing maps for genus-0 closed surfaces. The proposed algorithm is summarized in Algorithm 1. In practice, we set the step size , the stopping parameter , and the maximum number of iterations allowed .
3.2 Ellipsoidal density-equalizing quasi-conformal map (EDEQ)
While our EDEM algorithm can produce shape deformations and achieve prescribed area changes on the given ellipsoid , other geometric distortions such as conformal distortions may be uncontrolled. Also, note that the EDEM method can be applied to ellipsoids with any prescribed elliptic radii. It is therefore natural to ask whether one can further reduce the overall geometric distortions by optimizing both the shape of the ellipsoid and the associated ellipsoidal mappings.
As described in the DEQ method [42], the density-equalizing process can be considered as an energy minimization problem involving the prescribed density. Analogously, as our EDEM method aims to create a shape deformation following the density diffusion process on the ellipsoid, the density-equalizing effect of a map can be assessed using the following energy:
(31) |
where is the density associated with the map . Also, note that by quasi-conformal theory, the conformal distortion of a map can be represented using the norm of its Beltrami coefficient , where a smaller implies a smaller eccentricity of the local ellipses and hence a smaller angle distortion. Therefore, we can assess the conformal distortion of a mapping by considering the following energy:
(32) |
In this section, our goal is to develop an ellipsoidal density-equalizing quasi-conformal mapping algorithm, which we abbreviate as EDEQ, to reduce both of the above energies and . Now, it is noteworthy that the shape of the ellipsoid , characterized by the elliptic radii , can effectively provide us with extra degrees of freedom in reducing the energies. Therefore, we consider minimizing the following combined energy with the radii and the map being the variables:
(33) |
where is a nonnegative weighting parameter for balancing the density-equalizing error and the conformal distortion. The existence of the minimizer of can be proved by following the arguments in [42, 19].
3.2.1 Decoupling the combined energy
To simplify the optimization process, we first decouple the minimization problem of the combined energy into two subproblems of minimizing and below:
(34) |
(35) |
More specifically, the energy in Eq. (34) aims to equalize the density distribution across a fixed ellipsoidal domain with radii , promoting a more uniform density throughout the iterations. The energy in Eq. (35) focuses on optimizing the shape of the ellipsoid to minimize the combined geometric distortion. We can then utilize our developed computational scheme in the EDEM method for the subproblem and focus on updating the ellipsoidal geometry in the subproblem .
3.2.2 Subproblem
We first consider the descent direction of , with the elliptic radii fixed. As in the formulation of the EDEM method in Section 3.1, the velocity field can be computed by . Also, to preserve the given ellipsoidal shape, we remove the normal component , where is the outward normal unit vector. The descent direction of is then given by:
(36) |
Hence, we have the following iterative scheme for the energy :
(37) |
where is the time step size, is the map at the -th iteration, and is the descent direction of at the -th iteration.
3.2.3 Subproblem
To minimize the energy , we first consider the Beltrami coefficient term . Recall that the Beltrami coefficient is a complex-valued function defined on the complex plane. To extend it for assessing the conformal distortion of an ellipsoidal map, we use a triangle-based approach as follows. Let be a triangle element on the initial ellipsoidal conformal parameterization and be the corresponding triangle at the -th iteration. We rigidly map both and onto the complex plane. Specifically, this can be done by fixing one of the vertices of each triangle at the origin, one of the triangle edges containing this vertex onto the real axis, and then computing the position of the remaining vertex based on the edge lengths and angles of the triangle. We denote the resulting triangles on the complex plane as and . Then, we can consider the mapping between and as a quasi-conformal map . The Beltrami coefficient of can be obtained by:
(38) |
From the above formula, we can get the norm of the Beltrami coefficient and use it to represent the conformal distortion between and . We can then repeat the above procedure for all triangular faces of the surface to get the term .
Next, we consider optimizing the shape of the ellipsoid (i.e. the radii ) based on the energy . To simplify the formulation, we keep unchanged and focus on optimizing the two other elliptic radii and using an iterative procedure. We first define the radius step sizes along and as and , respectively. Now, there are nine possible combinations of the radii in the form of with . We can then consider the following energies associated with the nine combinations, including the original one:
(39) |
the ones with only one among and changed:
(40) |
and the ones with both and changed:
(41) |
Here, the terms and are computed based on the updated ellipsoid . By selecting the combination with the smallest energy value, we can update the shape of the ellipsoid and reduce . The shape update procedure is summarized in Algorithm 2.
3.2.4 Summary
In Sections 3.2.2 and 3.2.3, we discussed how to solve the subproblems and . Combining these approaches, we have our proposed ellipsoidal density-equalizing quasi-conformal map (EDEQ) algorithm.
Specifically, analogous to the EDEM algorithm, here we start by computing an initial ellipsoidal conformal parameterization (Section 3.1.1) with radii . The initial density and the initial combined energy can be obtained using the initial ellipsoidal parameterization. Then, to solve the subproblem , we deform the initial ellipsoid based on the descent direction (Section 3.2.2) iteratively. The overlap correction scheme (Section 3.1.3) and the re-coupling scheme (Section 3.1.4) are also applied at each iteration. We repeat the above process for a certain number of iterations (set by a prescribed parameter ) so that the energy defined based on the initial elliptic radii is sufficiently reduced. We then include an additional shape update step (Section 3.2.3) after every iteration to modify the ellipsoidal radii by solving the subproblem . With the elliptic radii updated, we solve the subproblem as described above again and repeat the above process. Also, note that a finer adjustment of the elliptic radii will be needed as the iterations continue. Therefore, instead of using fixed radius step sizes in the shape update algorithm, we add a scaling factor to the prescribed step sizes at the -th shape update step. We repeat the iterations until for some stopping parameter . The final EDEQ map is obtained by , from which we also get the final optimal elliptic radii .
The proposed EDEQ algorithm is summarized in Algorithm 3. In practice, the initial radii can be chosen arbitrarily, the time step size is set to be , the initial radius step sizes are set to be and , the number of iterations for each set of fixed radii is set to be , the error threshold is set to be , and the maximum number of iterations is set to be .
4 Experiments
In this section, we present experimental results to demonstrate the effectiveness of our proposed EDEM and EDEQ algorithms. The algorithms are implemented using MATLAB R2021a on the Windows platform. All experiments are conducted on a computer with an Intel(R) Core(TM) i9-12900 2.40 GHz processor and 32GB memory. The surface meshes are from online mesh repositories [49]. All surfaces are discretized in the form of triangular meshes.
4.1 Ellipsoidal density-equalizing map
To test our proposed EDEM algorithm, we first consider mapping an ellipsoidal surface with two different density distributions, including a discontinuous density distribution (Fig. 4(a)) and a continuous density distribution (Fig. 4(b)). In both examples, it can be observed that the input densities are highly non-uniform. By applying the proposed EDEM method, we obtain the corresponding ellipsoidal density-equalizing maps. It can be observed that the domains with a high initial density are enlarged while the domains with a low initial density are shrunk, which shows that our method can accurately produce shape deformations on the ellipsoid based on the prescribed density. Also, from the histograms of the initial and final densities, it can be observed that the density is effectively equalized.
Next, we consider some more complicated ellipsoidal mapping examples. Fig. 5(a) shows an ellipsoidal mesh with non-uniform distributed triangle elements. Specifically, the mesh is denser at the top part of the ellipsoid and much coarser at the bottom part. We define different populations on the ellipsoid such that the initial density of all the triangular regions is three times the density of all the pentagonal regions. From our EDEM mapping result, it can be observed that different regions are enlarged or shrunk accordingly. This shows that our method is capable of computing ellipsoidal density-equalizing maps on a non-uniform mesh. One may also wonder whether our method can be applied to more extreme ellipsoidal geometries and more extreme density distributions. In Fig. 5(b), we consider a more elongated ellipsoid with a more extreme density distribution prescribed on it, where the maximum and minimum density values differ by 10 times. It can be observed from our EDEM mapping result that different regions are effectively enlarged or shrunk based on the input density. This suggests that our method is capable of handling the extreme ellipsoidal geometry and extreme input density. Fig. 5(c) shows another ellipsoid with non-uniform mesh elements and a complicated continuous density distribution. Again, we can see that our proposed method produces an ellipsoidal density-equalizing map successfully. Lastly, in Fig. 5(d) we test our EDEM algorithm on a spherical surface (i.e. an ellipsoid with all three elliptic radii being identical). It can be observed from the mapping result that our method also works well for the spherical case.
For a more quantitative analysis, we record the computational time, radii of the ellipsoid, variance of the initial density, variance of the final density, and number of overlaps for all the above ellipsoidal mapping examples in Table 1. It can be observed that our EDEM method can be effectively applied to a wide range of ellipsoidal shapes with different elliptic radii. Specifically, for ellipsoids with different radii and input densities, the variances of the final density in the EDEM result are all close to 0. This shows that all EDEM results are highly density-equalizing. Also, it can be observed that all mapping results are folding-free. Overall, the experiments have demonstrated the versatility of the EDEM method.
Surface | # Faces | Time (s) | # Overlaps | |||
---|---|---|---|---|---|---|
Ellipsoid 1 (Fig. 4(a)) | 7808 | 1.8171 | (1,2,4) | 0.1176 | 0.0032 | 0 |
Ellipsoid 2 (Fig. 4(b)) | 7808 | 4.6973 | (1,2,4) | 0.2084 | 0.0028 | 0 |
Ellipsoid 3 (Fig. 5(a)) | 20480 | 1.4869 | (1,1,1.5) | 0.3333 | 0.0137 | 0 |
Ellipsoid 4 (Fig. 5(b)) | 18904 | 0.3811 | (1,0.6,2) | 0.2135 | 0.0055 | 0 |
Ellipsoid 5 (Fig. 5(c)) | 20480 | 2.6060 | (1,1,1.5) | 0.2580 | 0.0010 | 0 |
Ellipsoid 6 (Fig. 5(d)) | 18904 | 0.6203 | (1,1,1) | 0.2940 | 0.0068 | 0 |
Next, we consider computing the ellipsoidal density-equalizing maps for more general genus-0 closed surfaces. In Fig. 6, we set the population as the face area of the original mesh and apply our EDEM algorithm to achieve ellipsoidal area-preserving parameterizations. For the initial parameterization step, we map each surface conformally onto the ellipsoid using the FECM method [30]. Then, we apply the EDEM method to compute the ellipsoidal density-equalizing parameterizations. To assess the area-preserving property of the ellipsoidal parameterizations, we consider the area distortion for each triangular face as follows:
(42) |
Here denotes the ellipsoidal parameterization and is the set of all triangular faces. The normalization factors in the numerator and denominator ensure that a perfectly area-preserving mapping would result in . In Fig. 6, we present various genus-0 surface models, the initial ellipsoidal parameterizations, the ellipsoidal parameterizations obtained by our EDEM method, and the initial and final area distortions. It can be observed that the EDEM method performs very well for genus-0 closed surfaces with complex structures. Specifically, by visually comparing the triangle element distributions in the initial ellipsoidal parameterizations and the EDEM parameterization results, it is easy to see that the EDEM method achieves a more uniform distribution on the prescribed ellipsoid. By comparing the initial and final area distortion histograms, we can also see that the area distortion is significantly reduced by our EDEM method. For a more quantitative analysis, Table 2 provides detailed statistics on the performance of our EDEM algorithm for ellipsoidal area-preserving parameterization. We can see that the EDEM method can efficiently reduce the area distortion by over 95% on average when compared to the initial parameterization, and the bijectivity of the parameterization is well-preserved in all examples. Altogether, the experiments show that our EDEM method is capable of computing bijective ellipsoidal area-preserving parameterizations for a large variety of genus-0 closed surfaces.
Surface | # Faces | Time (s) | # Overlaps | |||||
Mean | SD | Mean | SD | |||||
David (Fig. 6(a)) | 21338 | 2.6982 | (1,1.2,1.4) | 0.7078 | 0.5248 | 0.0132 | 0.0136 | 0 |
Buddha (Fig. 6(b)) | 50002 | 10.5964 | (1,1,1.4) | 0.4365 | 0.3907 | 0.0366 | 0.0648 | 0 |
Hippocampus (Fig. 6(c)) | 12000 | 3.1002 | (1,1,2.2) | 1.3635 | 1.4740 | 0.0515 | 0.0973 | 0 |
Twisted Ball (Fig. 6(d)) | 38620 | 7.5060 | (1,0.8,1.4) | 0.8998 | 0.8195 | 0.0239 | 0.0289 | 0 |
Chinese Lion 1 (Fig. 1(a)) | 10000 | 10.1311 | (1,1,0.8) | 1.8306 | 1.8851 | 0.1407 | 0.2081 | 0 |
Chinese Lion 2 (Fig. 1(a)) | 10000 | 9.4919 | (1,1,1.1) | 1.6529 | 1.7923 | 0.1180 | 0.1759 | 0 |
Chinese Lion 3 (Fig. 1(a)) | 10000 | 7.9418 | (1,1,1.4) | 1.5394 | 1.8125 | 0.1137 | 0.1892 | 0 |
Mean | # Overlaps | |||||
1 | 1 | 1 | 0.0895 | 0.1684 | 0.4443 | 0 |
0.7 | 1 | 0.1349 | 0.3481 | 0.4620 | 0 | |
1.2 | 1 | 0.1042 | 0.2063 | 0.4446 | 0 | |
1 | 0.8 | 0.0996 | 0.1838 | 0.4709 | 0 | |
1 | 1.3 | 0.0816 | 0.1518 | 0.4017 | 0 | |
1 | 1.5 | 0.0730 | 0.1322 | 0.3730 | 0 | |
1 | 2.2 | 0.0515 | 0.0973 | 0.2862 | 0 | |
1.3 | 1.6 | 0.1484 | 0.3474 | 0.4026 | 0 | |
1.2 | 1.8 | 0.1340 | 0.3179 | 0.3942 | 0 | |
1.1 | 2 | 0.0975 | 0.2178 | 0.3718 | 0 |
It is natural to ask whether the performance of our EDEM method is affected by the shape of the target ellipsoid. Here, we map the Hippocampus model to various target ellipsoids with different elliptic radii using our EDEM method and analyze the area and angle distortions of the resulting mappings. In Table 3, we report the mean and variance of the area distortion, as well as the mean value of the norm of the Beltrami coefficient . It can be observed that for all combinations of the elliptic radii, all mapping results are folding-free and possess low area distortion. This shows that our EDEM method is capable of computing area-preserving parameterizations of genus-0 closed surfaces onto different prescribed ellipsoids. However, if we also take the angle distortion into consideration, then some differences can be observed. For instance, if we increase the value of , we can see that the angle distortion is further reduced. This experiment suggests that the overall ellipsoidal shape plays an important role and motivates the need for the proposed EDEQ method, which optimizes the ellipsoidal geometry and achieves ellipsoidal density-equalizing quasi-conformal maps.
4.2 Ellipsoidal density-equalizing quasi-conformal map
Next, we test our EDEQ algorithm for computing the ellipsoidal density-equalizing quasi-conformal mappings. In Fig. 7(a), we define different populations at different regions on the ellipsoid to obtain a discontinuous initial density. Note that the maximum density is three times larger than the minimum one. It can be observed in the EDEQ mapping result that the shape of the ellipsoid is changed under the mapping process, while the higher-density domain expands and the lower-density domain contracts. Comparing the initial and final density histograms, it can be observed that the density is highly equalized. Fig. 7(b) shows another example with a continuous initial density defined on an ellipsoid. Under the EDEQ algorithm, the ellipsoidal shape is stretched along the z-direction and becomes more elongated, and different regions are expanded or shrunk according to the given density. From the density histograms, it can again be observed that the mapping result is highly density-equalizing.
In Fig. 8, we further apply the EDEQ method to the same set of ellipsoidal examples in Fig. 5. It can be observed that our EDEQ method produces ellipsoidal density-equalizing mapping results with both the ellipsoidal geometry and the vertex positions optimized. For a more quantitative analysis, we record the initial and final elliptic radii, the variance of the initial and final densities, the mean value of the Beltrami coefficient , and the number of overlaps for all the above-mentioned examples in Table 4. From the table, we can see that the ellipsoidal shapes are effectively optimized under the EDEQ method. Furthermore, in addition to possessing a highly density-equalizing effect as reflected in the variance of the final density, the mapping results also achieve a low conformal distortion. Specifically, note that the norm of the Beltrami coefficient normally ranges from 0 to 1, where indicates that the mapping is conformal. As the value of is small for all examples, we can see that the EDEQ method not only achieves density equalization but also reduces the conformal distortion of the ellipsoidal mappings.
Surface | # Faces | mean() | # Overlaps | ||||
---|---|---|---|---|---|---|---|
Ellipsoid 1 (Fig. 7(a)) | 7808 | (1,1,2) | (1,1.8973,2.8163) | 0.2503 | 0.0049 | 0.1434 | 0 |
Ellipsoid 2 (Fig. 7(b)) | 7808 | (1,2,4) | (1,1.9888,4.3110) | 0.2084 | 0.0029 | 0.1953 | 0 |
Ellipsoid 3 (Fig. 8(a)) | 20480 | (1,1,1.5) | (1,1,1.4250) | 0.333 | 0.0153 | 0.1852 | 0 |
Ellipsoid 4 (Fig. 8(b)) | 18904 | (1,0.6,2) | (1,0.9821,2.5162) | 0.2135 | 0.0082 | 0.2285 | 0 |
Ellipsoid 5 (Fig. 8(c)) | 20480 | (1,1,1.5) | (1,1.0310,1.0783) | 0.2580 | 0.0011 | 0.1679 | 0 |
Ellipsoid 6 (Fig. 8(d)) | 18904 | (1,1,1) | (1,1.3095,0.9100) | 0.2940 | 0.0072 | 0.2267 | 0 |
Next, we consider using the EDEQ method to compute the ellipsoidal area-preserving parameterizations for genus-0 closed surfaces. For each surface in Fig. 9, we first conformally map it onto an ellipsoid with some arbitrary elliptic radii using the FECM method [30]. Then, we apply the EDEQ algorithm with the population set as the original face area to obtain the ellipsoidal area-preserving parameterizations. It can be observed that the ellipsoidal shapes are changed significantly under the EDEQ method. Comparing the initial and final area distortion histograms, we can see that the mapping results are highly area-preserving.
As the choice of the initial elliptic radii is arbitrary, we can even start with a unit sphere and run the EDEQ method to eventually obtain an optimal ellipsoidal parameterization. In Fig. 10, we consider several genus-0 closed surface models and compute an initial spherical conformal parameterization. Then, we apply the proposed EDEQ algorithm to compute an optimal ellipsoidal area-preserving parameterization. It can be observed that throughout the EDEQ iterative process, both the ellipsoidal shapes and the vertex positions on the ellipsoids are changed gradually. For a more quantitative analysis, Table 5 summarizes the performance of our EDEQ algorithm in terms of the final radii, the initial and final area distortions, the conformal distortion, and the number of overlaps. Comparing the initial and final area distortions, it can be observed that the area distortion is significantly reduced in all examples. The mean of the norm of the Beltrami coefficient is also small, indicating that the conformal distortion is low. Besides, the number of overlaps is 0 for all examples, which shows that our mapping results are all folding-free. In Table 6, We further compare the performance of the EDEQ method and the EDEM method in terms of the area and conformal distortions, where both methods use the same initial elliptic radii. It can be observed that the EDEQ method achieves a significantly lower conformal distortion when compared to the EDEM method, while the area distortion remains comparable. Altogether, the above experiments demonstrate the effectiveness of our proposed EDEQ method for ellipsoidal density-equalizing quasi-conformal maps.
Surface | # Faces | mean() | # Overlaps | |||||
Mean | SD | Mean | SD | |||||
Duck (Fig. 10(a)) | 20000 | (1,1.8797,1.8521) | 1.0561 | 1.2190 | 0.0241 | 0.0379 | 0.2216 | 0 |
Hippocampus (Fig. 10(b)) | 12000 | (1,0.1886,1.8818) | 2.5651 | 1.7272 | 0.1028 | 0.1990 | 0.2602 | 0 |
Lion-Vase (Fig. 10(c)) | 38620 | (1,0.3833,1.6136) | 0.8264 | 0.7195 | 0.1032 | 0.1749 | 0.2306 | 0 |
David (Fig. 1(b)) | 21338 | (1,1.0062,1.5126) | 0.5652 | 0.4349 | 0.0125 | 0.0133 | 0.1279 | 0 |
Surface | EDEM | EDEQ | ||||||
---|---|---|---|---|---|---|---|---|
Mean of | SD of | Mean | # Overlaps | Mean of | SD of | Mean | # Overlaps | |
Duck | 0.0268 | 0.0507 | 0.2467 | 0 | 0.0241 | 0.0379 | 0.2216 | 0 |
Hippocampus | 0.0895 | 0.1684 | 0.4443 | 0 | 0.1028 | 0.1990 | 0.2602 | 0 |
Lion-Vase | 0.0909 | 0.1604 | 0.3137 | 0 | 0.1032 | 0.1749 | 0.2306 | 0 |
5 Application to genus-0 surface remeshing
Note that one important application of surface parameterization is surface remeshing [50]. Specifically, after computing the parameterization of a given surface, one can generate a high-quality mesh on the parameter domain and use the inverse of the parameterization mapping to map the mesh onto the input surface, thereby yielding a remeshed surface. For instance, Sheffer and de Sturler [51] utilized the angle-based flattening method to compute parameterizations for surface remeshing. Praun and Hoppe [52] computed spherical parameterization for surface remeshing. Later, different point cloud parameterization methods [53, 54] have been developed and applied to surface meshing. Analogously, our proposed EDEM and EDEQ methods can be easily applied to genus-0 surface remeshing.
In Figs. 11, 12, 13, and 14, we consider various genus-0 surface models and compare the remeshing results obtained by different parameterization methods for genus-0 closed surfaces, including spherical conformal mapping (SCM) [9], spherical density-equalizing mapping (SDEM) [43], ellipsoidal conformal mapping (ECM) [30], our ellipsoidal density-equalizing mapping (EDEM) method, and our ellipsoidal density-equalizing quasi-conformal mapping (EDEQ) method. For each genus-0 surface model in and each parameterization method, denote the parameterization mapping by , where the parameter domain is either the unit sphere (for SCM and SDEM) or an ellipsoid (for ECM, EDEM and EDEQ). Then we construct a regular triangular mesh on using the DistMesh method [55]. Finally, using the inverse mappings , the regular mesh can be mapped back onto , yielding a remeshed surface . For a fair comparison, we enforce the number of vertices of the triangular mesh generated by DistMesh (and hence the number of vertices of the final remeshing result) for each method to be approximately 8500 (with variation ).
From the remeshing results in Figs. 11, 12, 13, and 14, it can be observed that the two conformal approaches (SCM and ECM) lead to significantly non-uniform remeshing results. Moreover, the shape features near the sharp regions are lost in the SCM and ECM results. This can be explained by the fact that the conformal mappings preserve angles but may yield large area distortions. By contrast, the remeshing results obtained by the density-equalizing approaches (SDEM, EDEM, and EDEQ) are much more uniform and have a higher mesh quality. The features located near the sharp regions are well-preserved.
To further compare the remeshing results obtained using the density-equalizing approaches, Fig. 15 shows the zoom-in images of the SDEM, EDEM and EDEQ remeshing results for the Pig, Bear, Hippocampus, and Lion-Vase models. From the zoom-in images, we can see that the EDEQ remeshing results are more uniform than the ones by the other two methods, especially at the regions corresponding to certain sharp features in the surface models. Moreover, comparing the mesh qualities in terms of the angles of the triangles, it can be observed that the SDEM remeshing results contain more skinny and irregular triangle elements when compared to the other two methods. This can be explained by the fact that the SDEM method maps all surfaces onto the unit sphere regardless of their geometries, which may lead to large geometric distortions. Also, while the EDEM method allows us to map the surfaces onto a prescribed ellipsoid, we may need to carefully choose a suitable ellipsoidal shape in order to reduce both the area and angle distortions of the mappings. By contrast, the EDEQ method optimizes both the ellipsoidal shape and the mapping onto it, giving the highest flexibility among the three methods. Therefore, it gives the most uniformly distributed and regular triangle elements in the remeshing results.
For a more quantitative analysis, we evaluate the surface remeshing performance of different approaches by considering the shape and size variations of the triangle elements in the remeshed surfaces. Ideally, the triangle elements on the remeshed surface should be as uniform as possible in terms of their size (i.e. the area of the triangles) and shape (i.e. the regularity of the triangles). To assess the size variation of the remeshed surface, we define the size variation measure as follows:
(43) |
where is the maximum triangle area and is the minimum triangle area in the remeshed surface. It is easy to see that for any remeshed surface, and the equality holds if and only if , i.e. all triangle elements of the remeshed surface are equal in size. This shows that can effectively capture the size variation in the remeshed surface.
To assess the shape variation, we first define the face regularity for each face of the remeshed surface as follows:
(44) |
where represent the length of the three edges of . Note that if and only if is an equilateral triangle. Then, we define the shape variation measure of the entire remeshed surface as
(45) |
It is easy to see that for any remeshed surface, and the equality holds if and only if all triangles in the remeshed surface are equilateral. Therefore, can effectively quantify the overall shape variation of the remeshed surface.
Table 7 records the values of and for all remeshing results shown in Figs. 11, 12, 13, and 14. For the two conformal approaches (SCM and ECM), it can be observed that the values of are small but the values of are extremely large, which means the remeshing results obtained by those two methods are highly non-uniform in size. By contrast, the values of for the three density-equalizing approaches (SDEM, EDEM, and EDEQ) are small. This demonstrates the advantage of utilizing density-equalizing maps for surface remeshing. Furthermore, by considering both and for the three density-equalizing methods, it can be observed that the EDEQ method gives the best overall remeshing performance among the three methods, which also matches our observations in Fig. 15.
Surface | Method | Shape variation | Size variation |
---|---|---|---|
Pig (Fig. 11) | SCM | 0.0499 | 4.5846 |
ECM | 0.0505 | 4.7123 | |
SDEM | 0.1562 | 3.3251 | |
EDEM | 0.1493 | 2.4706 | |
EDEQ | 0.1478 | 1.4218 | |
Bear (Fig. 12) | SCM | 0.0535 | 3.7168 |
ECM | 0.0526 | 3.9181 | |
SDEM | 0.1296 | 1.9349 | |
EDEM | 0.1282 | 1.4401 | |
EDEQ | 0.1263 | 1.3977 | |
Hippocampus (Fig. 13) | SCM | 0.0546 | 7.7204 |
ECM | 0.0549 | 7.0796 | |
SDEM | 0.2593 | 4.1060 | |
EDEM | 0.2428 | 1.3482 | |
EDEQ | 0.1568 | 1.2508 | |
Lion-Vase (Fig. 14) | SCM | 0.0587 | 5.5245 |
ECM | 0.0576 | 5.4088 | |
SDEM | 0.1744 | 2.9125 | |
EDEM | 0.1583 | 1.7410 | |
EDEQ | 0.1249 | 1.7081 |
6 Discussion
In this work, we have proposed a novel method for computing bijective ellipsoidal density-equalizing maps (EDEM) for genus-0 closed surfaces. Using this approach, we can efficiently obtain parameterizations with different desired mapping effects, including ellipsoidal area-preserving parameterizations and ellipsoidal parameterizations with controlled area change. Also, by considering a combined energy involving both a density-equalizing term and a quasi-conformal term, we can achieve ellipsoidal density-equalizing quasi-conformal maps (EDEQ) balancing the angle and area distortions. Experimental results on various ellipsoidal surfaces and genus-0 closed surface models have demonstrated the effectiveness of our proposed methods. Our proposed methods can be easily applied to surface remeshing of genus-0 closed surfaces, with the remeshing quality significantly improved when compared to other prior approaches.
We remark that while our proposed methods have considered both the density-equalizing property and the quasi-conformality, landmark-matching constraints have not been incorporated into our formulations. In the future, we plan to further extend our methods to facilitate landmark-matching surface parameterization and surface registration. Our models may also be combined with curvature-based methods to achieve curvature-based surface registrations with controlled area change for surfaces with extreme geometries.
References
- [1] M. S. Floater and K. Hormann, “Surface parameterization: a tutorial and survey,” in Advances in Multiresolution for Geometric Modelling, pp. 157–186, Springer, 2005.
- [2] 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, 2007.
- [3] K. Hormann, K. Polthier, and A. Sheffer, “Mesh parameterization: theory and practice,” in ACM SIGGRAPH ASIA 2008 courses, pp. 1–87, Association for Computing Machinery, 2008.
- [4] G. P. T. Choi and L. M. Lui, “Recent developments of surface parameterization methods using quasi-conformal geometry,” Handbook of Mathematical Models and Algorithms in Computer Vision and Imaging: Mathematical Imaging and Vision, pp. 1483–1523, 2023.
- [5] S. Haker, S. Angenent, A. Tannenbaum, R. Kikinis, G. Sapiro, and M. Halle, “Conformal surface parameterization for texture mapping,” IEEE Trans. Vis. Comput. Graph., vol. 6, no. 2, pp. 181–189, 2000.
- [6] X. Gu, Y. Wang, T. F. Chan, P. M. Thompson, and S.-T. Yau, “Genus zero surface conformal mapping and its application to brain surface mapping,” IEEE Trans. Med. Imaging, vol. 23, no. 8, pp. 949–958, 2004.
- [7] X. Chen, H. He, G. Zou, X. Zhang, X. Gu, and J. Hua, “Ricci flow-based spherical parameterization and surface registration,” Comput. Vis. Image Underst., vol. 117, no. 9, pp. 1107–1118, 2013.
- [8] K. Crane, U. Pinkall, and P. Schröder, “Robust fairing via conformal curvature flow,” ACM Trans. Graph., vol. 32, no. 4, pp. 1–10, 2013.
- [9] P. T. Choi, K. C. Lam, and L. M. Lui, “FLASH: Fast landmark aligned spherical harmonic parameterization for genus-0 closed brain surfaces,” SIAM J. Imaging Sci., vol. 8, no. 1, pp. 67–94, 2015.
- [10] G. P. T. Choi, Y. Leung-Liu, X. Gu, and L. M. Lui, “Parallelizable global conformal parameterization of simply-connected surfaces via partial welding,” SIAM J. Imaging Sci., vol. 13, no. 3, pp. 1049–1083, 2020.
- [11] W.-H. Liao, T.-M. Huang, W.-W. Lin, and M.-H. Yueh, “Convergence analysis of Dirichlet energy minimization for spherical conformal parameterizations,” J. Sci. Comput., vol. 98, no. 29, pp. 1–28, 2024.
- [12] Z. Su, W. Zeng, R. Shi, Y. Wang, J. Sun, and X. Gu, “Area preserving brain mapping,” in Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition, pp. 2235–2242, 2013.
- [13] X. Gu, F. Luo, J. Sun, and S. T. Yau, “Variational principles for minkowski type problems, discrete optimal transport, and discrete monge-ampere equations,” Asian J. Math., vol. 20, no. 2, pp. 383–398, 2016.
- [14] L. Cui, X. Qi, C. Wen, N. Lei, X. Li, M. Zhang, and X. Gu, “Spherical optimal transportation,” Comput. Aided Des., vol. 115, pp. 181–193, 2019.
- [15] A. Pumarola, J. Sanchez-Riera, G. P. T. Choi, A. Sanfeliu, and F. Moreno-Noguer, “3DPeople: Modeling the geometry of dressed humans,” in Proceedings of the IEEE International Conference on Computer Vision, pp. 2242–2251, 2019.
- [16] A. Giri, G. P. T. Choi, and L. Kumar, “Open and closed anatomical surface description via hemispherical area-preserving map,” Signal Process., vol. 180, p. 107867, 2021.
- [17] G. P. T. Choi, A. Giri, and L. Kumar, “Adaptive area-preserving parameterization of open and closed anatomical surfaces,” Comput. Biol. Med., vol. 148, p. 105715, 2022.
- [18] C. Gotsman, X. Gu, and A. Sheffer, “Fundamentals of spherical parameterization for 3D meshes,” ACM Trans. Graph., vol. 22, no. 3, pp. 358–363, 2003.
- [19] L. M. Lui, Y. Wang, T. F. Chan, and P. Thompson, “Landmark constrained genus zero surface conformal mapping and its application to brain mapping research,” Appl. Numer. Math., vol. 57, no. 5-7, pp. 847–858, 2007.
- [20] J. Lefèvre and G. Auzias, “Spherical parameterization for genus zero surfaces using laplace-beltrami eigenfunctions,” in International Conference on Geometric Science of Information, pp. 121–129, Springer, 2015.
- [21] S. Nadeem, Z. Su, W. Zeng, A. Kaufman, and X. Gu, “Spherical parameterization balancing angle and area distortions,” IEEE Trans. Vis. Comput. Graph., vol. 23, no. 6, pp. 1663–1676, 2016.
- [22] C. Wang, X. Hu, X. Fu, and L. Liu, “Bijective spherical parametrization with low distortion,” Comput. Graph., vol. 58, pp. 161–171, 2016.
- [23] G. P.-T. Choi, M. H.-Y. Man, and L. M. Lui, “Fast spherical quasiconformal parameterization of genus-0 closed surfaces with application to adaptive remeshing,” Geom. Imaging Comput., vol. 3, no. 1–2, pp. 1–29, 2016.
- [24] X. Hu, X.-M. Fu, and L. Liu, “Advanced hierarchical spherical parameterizations,” IEEE Trans. Vis. Comput. Graph., vol. 24, no. 6, pp. 1930–1941, 2017.
- [25] Z. Wang, Z. Luo, J. Zhang, and E. Saucan, “A novel local/global approach to spherical parameterization,” J. Comput. Appl. Math., vol. 329, pp. 294–306, 2018.
- [26] M.-H. Yueh, T.-M. Huang, T. Li, W.-W. Lin, and S.-T. Yau, “Projected gradient method combined with homotopy techniques for volume-measure-preserving optimal mass transportation problems,” J. Sci. Comput., vol. 88, pp. 1–24, 2021.
- [27] Z. Lyu, Q. Chen, and L. M. Lui, “A two-stage algorithm for combined quasiconformal and optimal mass transportation spherical parameterization,” Math. Comput. Geom. Data, vol. 3, no. 1, pp. 29–57, 2023.
- [28] T.-M. Huang, W.-H. Liao, and W.-W. Lin, “Fundamental theory and R-linear convergence of stretch energy minimization for spherical equiareal parameterization,” J. Numer. Math., vol. 32, no. 1, pp. 1–25, 2024.
- [29] J.-W. Lin, T. Li, W.-W. Lin, and T.-M. Huang, “Ellipsoidal conformal and area-/volume-preserving parameterizations and associated optimal mass transportations,” Adv. Comput. Math., vol. 49, no. 4, p. 50, 2023.
- [30] G. P. T. Choi, “Fast ellipsoidal conformal and quasi-conformal parameterization of genus-0 closed surfaces,” J. Comput. Appl. Math., p. 115888, 2024.
- [31] M. Shaqfa and W. M. van Rees, “Spheroidal harmonics for generalizing the morphological decomposition of closed parametric surfaces,” arXiv preprint arXiv:2407.03350, 2024.
- [32] M. T. Gastner and M. E. J. Newman, “Diffusion-based method for producing density-equalizing maps,” Proc. Natl. Acad. Sci., vol. 101, no. 20, pp. 7499–7504, 2004.
- [33] D. Dorling, M. Newman, and A. Barford, The atlas of the real world. Thames & Hudson., 2008.
- [34] M. T. Gastner, Spatial distributions: Density-equalizing map projections, facility location, and two-dimensional networks. Ph.D. Thesis, University of Michigan, 2005.
- [35] M. T. Gastner, V. Seguy, and P. More, “Fast flow-based algorithm for creating density-equalizing map projections,” Proc. Natl. Acad. Sci., vol. 115, no. 10, pp. E2156–E2164, 2018.
- [36] Z. Li and S. Aryana, “Diffusion-based cartogram on spheres,” Cartogr. Geogr. Inf. Sci., vol. 45, no. 5, pp. 464–475, 2018.
- [37] G. P. T. Choi and C. H. Rycroft, “Density-equalizing maps for simply connected open surfaces,” SIAM J. Imaging Sci., vol. 11, no. 2, pp. 1134–1178, 2018.
- [38] G. P. T. Choi, B. Chiu, and C. H. Rycroft, “Area-preserving mapping of 3D carotid ultrasound images using density-equalizing reference map,” IEEE Trans. Biomed. Eng., vol. 67, no. 9, pp. 1507–1517, 2020.
- [39] M. Shaqfa, G. P. T. Choi, G. Anciaux, and K. Beyer, “Disk harmonics for analysing curved and flat self-affine rough surfaces and the topological reconstruction of open surfaces,” arXiv preprint arXiv:2403.07001, 2024.
- [40] G. P. T. Choi and M. Shaqfa, “Hemispheroidal parameterization and harmonic decomposition of simply connected open surfaces,” arXiv preprint arXiv:2407.15417, 2024.
- [41] G. P. T. Choi and C. H. Rycroft, “Volumetric density-equalizing reference maps with applications,” J. Sci. Comput., vol. 86, no. 3, p. 41, 2021.
- [42] Z. Lyu, G. P. T. Choi, and L. M. Lui, “Bijective density-equalizing quasiconformal map for multiply connected open surfaces,” SIAM J. Imaging Sci., vol. 17, no. 1, pp. 706––755, 2024.
- [43] Z. Lyu, L. M. Lui, and G. P. T. Choi, “Spherical density-equalizing map for genus-0 closed surfaces,” SIAM J. Imaging Sci., to appear.
- [44] Z. Nehari, Conformal mapping. Courier Corporation, 2012.
- [45] O. Lehto and K. I. Virtanen, Quasiconformal mappings in the plane, vol. 126. Springer Berlin, Heidelberg, 1973.
- [46] L. V. Ahlfors, Lectures on quasiconformal mappings, vol. 38. American Mathematical Society, Providence, RI, 2006.
- [47] L. M. Lui, K. C. Lam, T. W. Wong, and X. Gu, “Texture map and video compression using Beltrami representation,” SIAM J. Imaging Sci., vol. 6, no. 4, pp. 1880–1902, 2013.
- [48] U. Pinkall and K. Polthier, “Computing discrete minimal surfaces and their conjugates,” Exp. Math., vol. 2, no. 1, pp. 15–36, 1993.
- [49] A. Jacobson, “Common 3D test models.” https://github.com/alecjacobson/common-3d-test-models, 2023. Accessed: 2023-10-03.
- [50] P. Alliez, G. Ucelli, C. Gotsman, and M. Attene, “Recent advances in remeshing of surfaces,” in Shape Analysis and Structuring, pp. 53–82, Springer, 2008.
- [51] A. Sheffer and E. de Sturler, “Parameterization of faceted surfaces for meshing using angle-based flattening,” Eng. Comput., vol. 17, pp. 326–337, 2001.
- [52] E. Praun and H. Hoppe, “Spherical parametrization and remeshing,” ACM Trans. Graph., vol. 22, no. 3, pp. 340–349, 2003.
- [53] M. Zwicker and C. Gotsman, “Meshing point clouds using spherical parameterization,” in Proceedings of the First Eurographics conference on Point-Based Graphics, pp. 173–180, 2004.
- [54] G. P.-T. Choi, K. T. Ho, and L. M. Lui, “Spherical conformal parameterization of genus-0 point clouds for meshing,” SIAM J. Imaging Sci., vol. 9, no. 4, pp. 1582–1618, 2016.
- [55] P.-O. Persson and G. Strang, “A simple mesh generator in MATLAB,” SIAM Review, vol. 46, no. 2, pp. 329–345, 2004.