A New Coordinate System for Constructing Spherical Grid Systems

: In astronomy, physics, climate modeling, geoscience, planetary science, and many other disciplines, the mass of data often comes from spherical sampling. Therefore, establishing an efﬁcient and distortion-free representation of spherical data is essential. This paper introduces a novel spherical (global) coordinate system that is free of singularity. Contrary to classical coordinates, such as Cartesian or spherical polar systems, the proposed coordinate system is naturally deﬁned on the spherical surface. The basic idea of this coordinate system originated from the classical planar barycentric coordinates that describe the positions of points on a plane concerning the vertices of a given planar triangle; analogously, spherical area coordinates (SACs) describe the positions of points on a sphere concerning the vertices of a given spherical triangle. In particular, the global coordinate system is obtained by decomposing the globe into several identical triangular regions, constructing local coordinates for each region, and then combining them. Once the SACs have been established, the coordinate isolines form a new class of global grid systems. This kind of grid system has some useful properties: the grid cells exhaustively cover the globe without overlapping and have the same shape, and the grid system has a congruent hierarchical structure and simple relationship with traditional coordinates. These beneﬁcial characteristics are suitable for organizing, representing, and analyzing spatial data.


Introduction
Astronomy, geoscience, computer graphics, and many other disciplines rely on research that uses spherical data.A fundamental problem lies in determining the optimal means of organizing, representing, or approximating these data.The most natural representation of spherical data is the spherical function f : S 2 → X or samples of the function f .Because of the intrinsic geometric properties of the sphere, the choice of an appropriate coordinate system is crucial for formulating an efficient and accurate numerical method for solving problems in spherical geometry.Spherical polar coordinates, although arguably the most natural coordinates for representing a point on a sphere p ∈ S 2 , present several disadvantages (collectively known as pole problems) when used for numerical computations over the entire spherical surface.The other approach to representing data is to map the spherical domain onto a plane M : S 2 → R 2 , such as the method used in traditional cartography.Then, points on a plane can be employed to process spherical signals.Unfortunately, in mathematics, the spherical surface is non-developable.That is, it is impossible to map a portion of the sphere onto the plane without introducing distortion (i.e., 'stretching' or 'compressing').For these reasons, a new class of spherical coordinates, spherical area coordinates (SACs), is presented in this paper.The SAC system is free of any singularity and naturally defined on the sphere surface, thus avoiding planar mapping distortion.
Research in science and engineering often needs a method for partitioning the sphere or some subset of it into sampling or analysis units (sphere tessellation or discretization).In geographical information systems (GISs), researchers have presented different kinds of Discrete Global Grid Systems (DGGSs) for organizing global data [1][2][3].The most commonly used regular DGGSs are those based on the geographic (latitude-longitude) coordinate system.Raster global data sets often employ cell regions with the edges defined by arcs of equal-angle increments of latitude and longitude (for example, the 2.5 • × 2.5 • and 10 • × 10 • NASA Earth Radiation Budget Experiment (ERBE) grids [4] and National Geophysical Data Center (NGDC) 5-minute Gridded Earth Topography Data [5]).Grid systems based on the geographic coordinate system have numerous practical advantages.The geographic coordinate system is in line with people's habits, and a wide array of existing data sets are based on it, so square grid regions are by far the most familiar to users.However, such grids have limitations: the cells of the grid become increasingly distorted in area, shape, and interpoint spacing as one moves north and south from the equator (Figure 1).The north and south poles, both points on the surface of the globe, map to lines on the longitude × latitude plane.The top and bottom row of grid cells are, in fact, triangles, not squares as they appear in the other row.These polar singularities have forced applications, such as global climate modeling to make use of specific grids for polar regions.To solve these problems, researchers have proposed another kind of discrete global grid system called the Geodesic Global Grid System (Geodesic DGGS) [6], which is based on the geodesic dome invented by [7].A geodesic dome is a spherical or partially spherical shell structure or lattice shell network on the surface of a sphere (Figure 2).These geodesic DGGSs are not a coordinate isoline grid but rather constructed by some specialized subdivision schemes [8].Furthermore, some of these grid systems have useful features, such as an equal-area hierarchical structure and regularity in shape.Rather than using great circle arcs for subdivision, Song et al. [9] proposed a small circle arc method that was optimized to achieve equal-area cell regions.However, the algorithms of some DGGSs are very complicated and have no analytical solution [1,3,9].(In our approach, we have straightforward subdivision schemes and clear transformation laws between a spherical point's Cartesian coordinate (x, y, z) and the SAC (λ 1 , λ 2 , λ 3 ).) Global Grid Systems are widely used in different areas of science and engineering.The cubed sphere, a well-known Cartesian-to-spherical mapping, was proposed by Ronchi et al. [10] and has been used to solve partial differential equations on the sphere.On the basis of the same grid system, Simons et al. [11] proposed a class of spherical wavelet bases for the analysis of geophysical models and the tomographic inversion of global seismic data.Furthermore, an axis-free grid named the Yin-Yang grid was used by Wongwathanarat et al. [12] to simulate 3D self-gravitating flows.In computer graphics, researchers have developed some discrete spherical wavelet bases that rely on geodesic sphere construction [13][14][15][16][17].They have been used successfully for spherical data visualization texturing and the bidirectional reflectance distribution function (BRDF).
In our approach, we adopt different strategies.First, we partition the globe into a series of spherical triangular cells without overlapping.Second, the local SACs on each triangular region are established, and then the global SACs can be constructed by gluing these local coordinate charts together along boundaries.Finally, we consider coordinate isolines to be grid cell edges and the intersection of isolines to be grid cell vertices, analogous to the Cartesian coordinate grid on a plane.The main properties of

•
The construction process is very simple; • Since our grid system has a slower growth rate between levels, it is more flexible and can be used for many applications;

•
The grid cells have a low distortion in shape and area;

•
The grid cells exhaustively cover the globe without overlapping; • The grid system has a congruent hierarchical structure; • The grid system has a simple relationship with conventional coordinates.
Our main contribution in this paper is the proposal of a new method to construct spherical area coordinates and a new kind of DGGS that is based on these coordinate systems.

A Brief Review of the Geodesic DGGS
The following introductory paragraphs are intended to acquaint the reader with the concepts associated with the Geodesic DGGS.
A Discrete Global Grid (DGG) consists of a set of regions that form a partition of the sphere S 2 , where each region has a single point contained in the associated region.Each combination of a region and point is called a cell .A series of increasingly finer-resolution DGGs composes the DGGS.If the grids are defined consistently using regular planar polygons on a surrogate surface, the aperture of a DGGS is defined using the ratio of the areas of a planar polygon cell at resolutions k and k + 1 [18].The aperture of the system reflects the relationship between consecutive levels of the grids.Regular hierarchical relationships between DGGS resolutions play an important role in creating efficient data structures [19,20].Two types of hierarchical relationships, congruent and aligned, are commonly used [6].
Consider Cr k to be a cell region and Cp k to be a cell point of a cell at resolution k.
• A DGGS is congruent if and only if If a DGGS does not have these properties, the system is defined as incongruent or unaligned.In Sahr's study [6], the construction of the Geodesic DGGS was treated as a series of design choices that were, for the most part, independent.The following five design choices fully specify a Geodesic DGGS: 1.
Choice of cell points.

Base Regular Polyhedron
The first step in the construction of a Geodesic DGGS is to establish an initial partition on the globe.Since multi-resolution discrete grids are obtained by subdividing the initial grid, the cells of the initial grid should preferably be the same regular spherical polygons.A geometrically rigorous solution to the problem of covering a spherical surface with tiles shows that there are only five classes of inscribed polyhedrons: the tetrahedron, the cube, the octahedron, the dodecahedron, and the icosahedron.All these polyhedrons produce a covering that is both regular and homogeneous.The covering can be obtained by a central projection of the sides and corners of the inscribed polyhedron on the spherical surface (Figure 3).These five platonic solids represent the only ways in which the sphere can be partitioned into cells, each consisting of the same regular spherical polygon, with the same number of polygons meeting at each vertex.Platonic solids have thus been commonly used to construct Geodesic DGGSs.In general, platonic solids with smaller faces reduce the distortion introduced when transforming between a face of the polyhedron and the corresponding spherical surface.The icosahedron is thus the most commonly chosen base regular polyhedron.

Polyhedron Orientation
When a base polyhedron is chosen, a fixed orientation relative to the actual surface of the sphere must be determined.The best choice for the orientation of a base polyhedron mostly depends on the background.In geoscience, Fekete and Treinish [21] and Thuburn [22] placed a vertex of the icosahedron at each of the poles and then aligned one of the edges emanating from the vertex at the north pole with the prime meridian; this is the most common orientation.Fortunately, regardless of the chosen orientation, the geometric properties of the grid are not changed.

Subdivision Scheme
In the next step, we need to choose a method for subdividing this polyhedron to create multiple resolution grids.When the subdivision scheme has been chosen, one can partition a single face of the polyhedron or a set of faces that constitute a unit that tiles the polyhedron.A number of approaches to subdividing the surface of a polyhedron have been suggested, and the approaches can be classified according to different criteria.There are four types of subdivision schemes, and they have different partition shapes: triangles, squares, diamonds, and hexagons.Depending on the aperture of the grids, the approaches can be classified into two classes, namely, 4-fold and 9-fold.According to the kind of line that we use to connect the subdivision points, the approaches can be classified into two categories: Class I (connecting the subdivision points with lines parallel to the edges of the polyhedron face) and Class II (connecting the subdivision points with lines perpendicular to the edges of the polyhedron face).Figure 4 shows the most commonly chosen subdivision scheme for DGGSs.There are, however, other partition shapes.In this study, we investigated partitions that are based on the triangular face, but other partition shapes can be chosen.

Transformation
When choosing the partitioning method, we must select a method for creating a similar shape on the corresponding spherical surface.There are two types of approaches: Direct Spherical Subdivision (DSS) and Map Projection [20].DSS approaches involve the creation of a partition directly on the spherical surface that maps to the corresponding partition on the planar face(s).The simplest direct subdividing method is using great circle arcs that correspond to the cell edges on the planar face(s).It should be noted that the aperture 4 Class I DSS is actually the geodesic bisector subdivision.This partitioning method has been used extensively to construct spherical wavelets [13][14][15][16].Rather than using great circle arcs, Song et al. [9] proposed using small circle arcs that were optimized to achieve equal-area cell regions.Map projection approaches use an inverse map projection to transform a partition defined on the planar face(s) to the sphere.The most popular projections are the Gnomonic azimuthal projection, the Fuller Dymaxion projection [3,7], the Quaternary Triangular Mesh (QTM) [2,23], and the Snyder equal-area polyhedral projection.

Choice of Cell Points
Once the spherical grids have been specified, the points associated with the grid cells are usually chosen to be the center points of the cell regions.In map projection approaches, it is often convenient to choose the center points of the planar cell regions.In DSS approaches, the choice of the center points may be complicated.

Previous Work
Numerous spherical surface tessellations have been proposed to establish DGGSs.The most commonly used DGGS is the traditional latitude-longitude graticule, in which they are termed equal-angle grids.Another approach is the quadrilateral cells, which are called "constant-area" grids [24].Constant-area tessellations begin with an arbitrarily sized quadrilateral cell at the equator and then define the parallel and meridian cell boundaries across the globe to achieve approximately equal-area cells.Fekete's quad-tree data structure [21] and Song's small circle subdivision [9] are both based on the recursive subdivision of spherical triangles obtained by projecting the faces of a regular polyhedron onto a sphere, and the edges of the cells are all great and/or small circle.The Hipparchus System [25] allows for the creation of arbitrarily regular DGGs by generating Voronoi cells on the surface of a sphere from a specified set of points.However, for many applications, it is desirable to have cells that consist of highly regular regions.The approaches discussed above are all a DSS, but there are some popular map projection methods.In the common Gnomonic azimuthal projection, all straight lines in the plane correspond to great circle arcs on the sphere but exhibit relatively large area and shape distortions and the geodesic bisector method give the same results.The Snyder equal-area polyhedral projection [1] is defined on all regular polyhedra, but it has greater shape distortion.The Dymaxion projection [7] given in [3] has less shape distortion than Snyder's projection, and it is not an equal-area projection.Dutton's Quaternary Triangular Mesh (QTM) [26] is based on the 4-fold recursive subdivision of the octahedron, and two projections have been proposed for transforming the QTM from the plane to the sphere: Goodchild and Yang's Plate Carree projection [2] and Dutton's Zenithial OrthoTriangular (ZOT) projection [23].The difference between these two projections is that they use different planar surfaces.Dutton's method uses an isosceles right triangle, and Goodchild and Yang's method uses an equilateral triangle; however, both projections give identical geographic coordinates when taken to the sphere.

Spherical Area Coordinates (SACs)
Let us now more intuitively and visually introduce the basic concept that underlies the construction of the SAC system.The fundamental idea of our approach is inspired directly by F. Möbius's work: barycentric coordinates [27].

Planar Area Coordinates
Let P = (v j ) j=1...n be a polygon with vertices v j .The barycentric coordinates λ i (v; P) of a point v are continuous functions that satisfy the following properties for all points v inside the polygon.
Figure 5 illustrates the idea of barycentric coordinates for the case of planar triangles.Any point v on the plane E 2 can be expressed by a linear combination of the vertices of the reference triangle ).The point v outside the reference triangle can lead to negative coordinates, but this situation is not considered in this paper.Now, the question is, 'What are the continuous functions λ i (v; P) of the barycentric coordinates of point v?' In the context of a triangle, barycentric coordinates are also known as area coordinates because the coordinates of v with respect to the triangle v 1 v 2 v 3 are proportional to the areas of S vv 2 v 3 ,S v 1 vv 3 and S v 1 v 2 v .The equations of area coordinates are written as So, according to Equation (1), we can easily calculate the area coordinates of points v 1 , v 2 , and v 3 .The area coordinates of v 1 are Since the triangles v 1 v 1 v 3 and v 1 v 2 v 1 are degenerate, the area of these two triangles is 0. As a result, the area coordinates of v 1 is v 1 (λ 1 , λ 2 , λ 3 ) = (1, 0, 0).What is more, the area coordinates of the verties v 2 and v 3 are v 2 (λ 1 , λ 2 , λ 3 ) = (0, 1, 0), and v 3 (λ 1 , λ 2 , λ 3 ) = (0, 0, 1).

Local SACs
Spherical barycentric coordinates were first studied by Möbius [28] and introduced to computer graphics by Alfeld et al. [29].In these papers, the spherical barycentric coordinates express the location of a point v on a sphere to the vertices v i of a given spherical triangle P. In fact, Bronw and Worsey [30] showed that such coordinates that satisfy properties 1, 2, 3, just as planar barycentric coodinates, do not exist.Langer et al. [31] chose to relax property 2 to ∑ i λ i ≥ 1 and to preserve property 3 because the main issue that they aimed to address was defining curves and surfaces over the sphere.Since our approach focuses on the problem of spherical surface tessellation, we chose different strategies: relax property 3 and preserve property 2.

Definition of SACs on the Spherical Triangle
In the context of SACs, the reference triangle is the spherical triangle T. The reference triangle can be any non-degenerate spherical triangle.Any area coordinates of the point in this reference triangle can be calculated through the vertices of this triangle.If there is a point p inside T, then we can divide T into three sub-triangles by connecting p and the vertices of the reference triangle p i , i = 1, 2, 3 with great circular arcs C i , i = 1, 2, 3 (see Figure 6a).Definition 3.1.A system of spherical area coordinates over a geodesic triangle T on S 2 is a one-to-one correspondence between points p ∈ T and ordered triples of real numbers u = (λ 1 , λ 2 , λ 3 ) such that ∑ i λ i = 1; 3.
where T 1 = (p, p 2 , p 3 ),T 2 = (p 1 , p, p 3 ),T 3 = (p 1 , p 2 , p).It is obvious that the sub-triangles are all spherical triangles.Until this step, the coordinates of point p can be written in the planar form (Equation ( 1)), but the issue lies in the method for calculating the areas of these four spherical triangles.Let a spherical triangle be drawn on the surface of a sphere of radius R, centered at a point O = (0, 0, 0) (in Cartesian coordinates), with vertices P 1 , P 2 , and P 3 .Therefore, the vector from the center of the sphere to the vertices is given by p 1 = OP 1 ,p 2 = OP 2 , and p 3 = OP 3 .We assume that the area of the spherical triangle is S, with R = 1.According to Girard's theorem in spherical trigonometry, the area of a spherical triangle on a sphere of radius R is where α, β, and γ denote the angles between the sides (dihedral angles) in radians (shown in Figure 6b), and the angle α is opposite the side that subtends the greater circle arc C 1 , and so forth.Hence, the SACs can be written as When the spherical reference triangle T is determined, the area of T is a constant, and S i , i = 1, 2, 3, is a continuous function about the spherical point p.According to cosine laws of spherical trigonometry, The definitions of the dot product and cross-product are Note that p = p 1 = p 2 = p 2 = 1.Then, we obtain the formula below.
and so forth.

Transformation Laws between SACs and Cartesian Coordinates
Because Cartesian coordinates have been widely used in many areas of science and engineering, the transformation between Cartesian coordinates and SACs must be established.However, Equation ( 3) is too complicated for obtaining the transformation; thus, we need a simpler relationship between the spherical triangle area and its vertices.In spherical trigonometry, the following identity can be derived from sine and cosine laws [32]: From Equation ( 3), ( 6), we obtain If the vectors in Equation ( 7) are given in the Cartesian coordinate system, then it determines the transformation rules from Cartesian coordinates to SAC.
Given a point P at u = (λ 1 , λ 2 , λ 3 ), if the Cartesian coordinates of the vertices p 1 , p 2 , p 3 are known, then Equation ( 7) can be used to obtain Equation ( 8) can be written in the matrix form: where p = (x, y, z), and matrix A and vector k are given by Finally, we obtain transformation laws between the Cartesian coordinates and SACs.

The Global SAC and Grid Systems
From the discussion above, the SACs are defined only in a portion of the sphere (inside the reference triangle).In other words, they are local coordinates.In this section, we introduce the extension of the SACs to the entire sphere.The first step is to partition the globe into several cells (the initial global grid), each consisting of the same regular spherical triangle.Further, the cells constitute a complete tiling of the globe, exhaustively covering the globe without overlapping.This criterion ensures that the sum of cell areas equals the surface area of the globe and that the coordinates of all shared cell edges are identical.Since cells of the initial global grid are spherical triangles, we can treat them as reference triangles and establish local SACs for each cell.When putting these patches together, we finally create a global coordinate system over the sphere.We expect that the initial global grid is regular.Since we require that the grid cells be a triangular cell, we now have only three options: the spherical tetrahedron, octahedron, and icosahedron versions.After the global SAC system has been developed, the coordinate grid naturally spans the globe.

The Tetrahedron Version of SACs
For the first example, we show the development of the tetrahedron version of the SAC system.To acquire an elegant form of the transformation, we place the vertices of the tetrahedron at The initial global grid consists of four equal-area spherical equilateral triangles, T i , i = 1, 2, 3, 4 (shown in Figure 7), and the area is π.Now, we use the method in Section 3.2 to establish the coordinates on each cell.The tetrahedron version of the SACs is an ordered tetrad of real numbers (λ i ), i = 1, 2, 3, 4, and we specify that Then, substituting the coordinates of points A, B, C, D into Equation ( 7) and rearranging it, we have p ∈ T 1 : p ∈ T 2 : ; ; . p ∈ T 3 : p ∈ T 4 : ; ; Note that if the point lies on the edges of T i , then two of the coordinate components will be zero: Finally, every point on the sphere corresponds to unique SACs (λ i ).To obtain the inverse transformations, we substitute A, B, C, D into Equation ( 9) and solve the equation.

The Octahedron Version of the Global Grid System
We can also define SACs on the octahedron.First, we need to orient the coordinate sphere triangles.The octahedron is oriented in the following way: the center of the sphere is at the (x, y, z) coordinates (0, 0, 0), and the six vertices are at and V 6 = −V 3 (shown in Figure 11).Then, we divide the sphere into eight spherical triangles.We denote these regions as N i (north) and S i (south), i = 1 . . . 4. Therefore, any point on the sphere surface must be in one of these eight cells, and the cells can serve as the reference triangle.Because the derivation is the same as the tetrahedron version, the forward transformation laws have the following forms: Different from the tetrahedron version of SACs, the octahedron version of SACs have three components (λ 1 , λ 2 , λ 3 ).In Equation ( 13), λ 1 , λ 2 , and λ 3 are all positive.The SACs in different regions have different sign patterns and are determined by the regions containing the point P. The sign patterns are the same as the octants in the 3D Cartesian system (illustrated in Figure 12).systems.This subsection will give you some visual results that use SACs to create hexagonal based global grid.Figure 16 illustrates three levels of hexagonal based DGGSs that are based on SACs.Hexagonal based grid has some advantages.They are more compact than triangles and squares and also have the smallest average error.Because hexagons are more symmetry, they provide the highest angular resolution.Unlike square and triangle grids, hexagon grids have uniform adjacency; each hexagon cell has six neighbors.But the hexagonal based grids also have some disadvantages.The grid's cells are not all hexagon, some pentagons exist in the initial vertices of Platonic solid.A simple data structure, like quadtrees, octrees cannot use to index the grid systems.In this paper, we focus on triangular based grid systems.

DGGS Evaluation Criteria
When there is a new global grid, a common question arises: What is the difference between this grid and existing ones?What are the advantages of this grid?Some researchers have proposed a series of standards for evaluating global spherical grids.These standards include the following: 1.
The grid cells complete tiling the globe without overlapping in any resolution; 2.
The grid cells at one resolution have equal areas; 3.
The grid cells have the same topology; 4.
The grid cells have the same shape; 5.
The grid system has a simple relationship with the latitude and longitude graticule; 6.
The grid system contains grids of any defined spatial resolution; 7.
The grid cells are compact.
Table 1 gives the result between different methods based on the criteria.

Performance Analysis
We compared the performance of different DGGSs by measuring two general properties: area and shape distortions.We used the perimeter and area as the two direct measurements of partition cells.Since the edges of the partition cells are not always great circle arcs in all methods, we used an approximation method to compute the perimeter and area of the partition cells.At all levels of subdivision, we calculated the resolution of representation as not instantly level but level 2 (16−subdivision level) of the partition cells.A large number of flat triangular facets approximate a spherical triangle.We calculated the perimeters as the sum of the lengths of straight lines, and we calculated the area as the sum of the area of small planar triangles.For each comparison, we started with a set of partition cells on the sphere that were generated from one of the three polyhedrons, one of the two density ratios, one of four methods (including our approach) for mapping to the sphere, and one of the eight or five levels of recursion, depending on the density ratio.
For a more natural interpretation of the comparisons, we used standardized measures.For area measurements, we converted the actual areas to regulated areas for each level of recursion by dividing the exact areas by the values listed in Table 2 to obtain regions relative to a unit value of 1.For shape measurements, we used a standardized perimeter-to-area ratio and called this the 'Zone Standardized Compactness' (ZSC) [20] using the equation below.
where a is the cell area, p is the cell perimeter, and r is the radius of the sphere.The resulting values were then dimensionless numbers between 0 and 1.
In this subsection, we compare the distortion performance of our approach with that of some classical DGGSs.We implemented two recursive partitioning methods (4-fold, 9-fold) and three methods of mapping the partition from the polyhedron to the surface of the sphere.These three methods are the Snyder projection method, the geo bisector (geodesic bisector) method, and the Dymaxion projection.All these global partitions are based on the tetrahedron, octahedron, and icosahedron.The standardized surface area was derived by dividing the surface area of the sphere with radius = 1 by the number of cells.

Measures of Performance
Figure 17 provides the visual results of the four DGGSs based on different methods.Our DGGSs have, on average, the most compact cells and smaller area distortions, except for the Snyder equal-area DGGSs.Since equivalent surface areas for the cells are maintained, Snyder's DGGSs exhibit more shape distortions, especially in the cells along the radial lines that extend from a polyhedron face center to its vertices.We used five statistics of shape and area measurements: the range of the different area (standardized area), the standard deviation of the area, the range of compactness, the standard deviation of compactness, and the mean compactness.Figure 18-20 illustrate the performance of DGGSs by the level of resolution for each of the five response measures.18a, 19a, 20a).The greatest area distortion is found with the Geo Bisector method, followed by (in descending order) the Dymaxion projection method, SACs, and Snyder.The Snyder projection has an area range of zero, of course, since it is an equal-area mapping method.
For the compactness measures, the Snyder method has the highest scores both in range and standard deviation (Figures 18c, 19c, 20c).These results indicate that the Snyder equal-area projection has the most shape distortion.The compactness performance of our method is almost the same as that of the Dymaxion method.As the mean of the compactness increases, it approaches the most compact shape possible.The results of the experiments for this measure show that the Geo Bisector has the best performance (highest values).
The histograms of the distributions of the area bias and compactness values (in Figures 21 and 22) show distinct differences between methods.The Snyder projection histogram for the area has only one bar for the single equal-area value (no bias).Our method has a continuous, unimodal histogram.It also has a better shape distortion performance than the Dymaxion and Geo Bisector methods on the tetrahedron (Figure 21a).For the icosahedron, no matter what fold method is used the histogram for our method has a smoother distribution (Figures 22a and23a).The maps of the spatial variation in the area and ZSC values have interesting patterns (in Figure 24-26).For the 4-fold density ratio, the spatial variation in the area (Figures 24a and 25a) has one of two patterns.The Snyder method result is homogeneous with an equal area.The other three methods created partitions with a decrease in area per cell radially from the center.For compactness, each method has a unique pattern (in Figures 24b and 25b).For the 9-fold density ratio, all methods have patterns of area and compactness variation that are very similar to those for the 4-fold ratio (in Figure 26).

Discussion
This paper presents several complementary views of the distortion performance of different methods for partitioning the triangular faces of polyhedrons.These views show that our method has a generally smoother distribution in both area and shape measurements, especially for the tetrahedron.
Another key advantage of our method is not only the association of points but also the continuous functions with the cells.In classical DGGSs, each cell can only associate with one point.If more than one point needs to be represented, a higher-resolution DGGS is required.

Conclusions and Further Work
In this paper, we detail a method for defining SACs on the globe.Further, we compare our approach with some classical DGGS methods.We also compare the performance between SAC DGGSs and three classical DGGSs.The analysis results illustrate that the DGGSs based on SACs have an average performance.This means that they have small area and shape distortions simultaneously.Another advantage is the explicit transformation equation between the Cartesian coordinates and SACs.
We show that the SAC system is a global coordinate system and can be used to represent different kinds of spherical signals: continuous function (Figure 27) and discrete sampling data (Figure 28).On the other hand, contrary to many other approaches, such as the methods of [1,9], the transformations between Cartesian coordinates and SACs are analytical and compact.Because the global grid system constructed by our method is triangular, techniques that have been proposed for processing the triangular grid [33] can also be used to deal with our system.
It would be interesting to establish a general theory of SACs for spherical signal processing.Could the SACs be used in spherical CAGD (Computer Aided Geometric Design) or solve partial differential equations in spherical geometry?Another future work is how to describe the surface with a single binary number.Using the binary description it would be much more computationally efficient to calculate gradients, etc.These topics need exploration in future studies.

Figure 1 .
Figure 1.The distortions of the geographical coordinate grid.The top and bottom row of grid cells are, in fact, triangles, not squares as they appear in the other row.

Figure 3 .
Figure 3. Spherical tessellations of the five regular polyhedra.

Figure 4 .
Figure 4. Four-and 9-fold triangular and hexagonal partitioning of a single triangular face.

Figure 5 .
Figure 5. Barycentric coordinates (in the case of planar triangles).
(a) The basis of spherical area coordinates (SACs): a point p in the spherical reference triangle.(b) The notation of a spherical triangle.

Figure 7 .
Figure 7.The tetrahedron version of the initial global grid.

Figure 11 .
Figure 11.The octahedron version of the initial global grid.

Figure 17 .
Figure 17.The four polyhedral DGGSs.Each figure shows a 4-fold, level-5 triangle subdivision on polyhedrons.The surfaces are shown in an orthographic projection to accentuate the distortion patterns present in each method.

Figure 18 .
Figure18.Graphs of the five statistics of the area and Zone Standardized Compactness (ZSC) response measurements for each method as a sequence by recursion level (level 1 to level 8).These results are for the tetrahedron using the 4-fold density ratio.

Figure 19 .
Figure19.Graphs of the five statistics of the area and ZSC response measurements for each method as a sequence by recursion level (level 1 to level 8).These results are for the icosahedron using the 4-fold density ratio.

Figure 20 .
Figure20.Graphs of the five statistics of the area and ZSC response measurements for each method as a sequence by recursion level (level 1 to level 5).These results are for the icosahedron using the 9-fold density ratio.Differences in the range of the area distortion in the different methods are clearly distinguished as the resolution increases (shown in Figures18a, 19a, 20a).The greatest area distortion is found with the Geo Bisector method, followed by (in descending order) the Dymaxion projection method, SACs, and Snyder.The Snyder projection has an area range of zero, of course, since it is an equal-area mapping method.For the compactness measures, the Snyder method has the highest scores both in range and standard deviation (Figures18c, 19c, 20c).These results indicate that the Snyder equal-area projection has the most shape distortion.The compactness performance of our method is almost the same as that of the Dymaxion method.As the mean of the compactness increases, it approaches the most compact

Figure 21 .
Figure 21.The distribution histograms of area and compactness for different methods.These values are for the tetrahedron partitions.

Figure 22 .
Figure 22.The distribution histograms of area and compactness for different methods.These values are for the icosahedron partitions.

Figure 23 .
Figure 23.The distribution histograms of area and compactness for different methods.These values are for the icosahedron partitions.

Figure 24 .
Figure 24.Maps of the spatial variation in the area and ZSC measurements for the Geo Bisector, Snyder, Dymaxion, and SAC methods for the tetrahedron and 4-fold density ratio.

Figure 25 .
Figure 25.Maps of the spatial variation in the area and ZSC measurements for the Geo Bisector, Snyder, Dymaxion, and SAC methods for the icosahedron and 4-fold density ratio.

Figure 26 .
Figure 26.Maps of the spatial variation in the area and ZSC measurements for the Geo Bisector, Snyder, Dymaxion, and SAC methods for the icosahedron and 9-fold density ratio.

Figure 27 .
Figure 27.The shape of the squirrel is composed of 10 parametric curves (equations).In the planar triangle, there is the planar version of the curves that use the same form of expressions as the spherical version.

Table 1 .
Criteria comparison between different methods.

Table 2 .
Cell counts and standardized surface areas for DGGSs based on a single face of the spherical tetrahedron, octahedron, and icosahedron.