The rapid growth in Earth’s global geospatial data necessitates an efficient system for organizing the data, facilitating data fusion from diverse sources, and promoting interoperability. Mapping the spheroidal surface of the planet presents significant challenges as it involves balancing distortion and splitting the surface into multiple partitions. The distortion decreases as the number of partitions increases, but, at the same time, the complexity of data processing increases since each partition represents a separate dataset and is defined in its own local coordinate system.
1. Introduction
Organizing and referencing geospatial data pose increasingly complex challenges due to the sheer volume and rapid growth of the data collected. For example, the daily influx of high-resolution satellite imagery alone amounts to terabytes of data. It is important to store these data in a format that allows easy access, referencing, sharing, and analysis without frequent re-projections to maintain accuracy.
Addressing these challenges requires the development of a spatial reference frame capable of fusing data from diverse sources into a global mosaic at multiple resolutions. Discrete Global Grid Systems (DGGSs) have emerged as a promising class of such reference frames. They use hierarchical tessellation to partition and address the entire planet without gaps or overlaps. The development of DGGS began in the mid 20th century
[1] but only became popular at the turn of the century
[2][3][4]. A prominent subclass of DGGSs are Geodesic DGGSs (GDGGSs)
[4], which project the surface of the planet onto the faces of regular or semi-regular circumscribed polyhedra. Commonly used base polyhedra are the five Platonic solids, in particular the hexahedron (cube), while the most common semi-regular polyhedron is the truncated icosahedron. The appearance of unfolded regular polyhedra, the number and shape of faces, and the way they partition Earth’s surface can be observed in
Figure 1.
Figure 1. Unfolded regular polyhedra used in GDDGSs: (a) tetrahedron; (b) hexahedron; (c) octahedron; (d) dodecahedron; (e) icosahedron.
While increasing the number of faces in a polyhedron improves its approximation of a spherical surface, it introduces complications when merging adjacent partitions due to separate coordinate systems and datasets. Moreover, memory and CPU consumption increase with the number of partitions in systems handling streaming geospatial data, like out-of-core planet-scale terrain rendering applications
[5]. Therefore, minimizing the number of partitions becomes desirable.
2. Development of Discrete Global Grid Systems (DGGSs)
The need to develop a system for organizing global geospatial data is not new. One of the first studies on the feasibility of implementing the Earth Data Base System
[6] was conducted in the early 1970s for the needs of the US Navy. The system was based on Quadrilateralized Spherical Cube (QSC), one of the first hexahedral projections implemented on digital computers. The proposed system was soon modified
[7], and, in the following years, it was also used as part of the Cosmic Background Explorer (COBE) project at NASA.
Due to the regular and uniform structure of the grid consisting of square cells, consistency with the Cartesian coordinate system, ease of interpolation and extrapolation, and straightforward visualization, hexahedral projections have gained wide popularity and are used in many different fields.
Hierarchical Equal-Area isoLatitude Pixelization (HEALPix) is a class of spherical projections with the property of distributing 12N
2 points as uniformly as possible over the surface of the unit sphere
[8]. These hybrid projections combine the Lambert cylindrical equal-area projection for the equatorial region with the interrupted Collignon projection for the polar regions. Of this infinite class of projections, only the projection with three base resolution pixel layers between the north and south poles and four equatorial base resolution pixels can be rearranged to a hexahedral projection.
Rotated HEALPix (rHEALPix)
[9] is an extension of the HEALPix scheme that introduces rotation capabilities, is better adapted to standards, and inherently combines polar triangles into quadratic partitions. Further, rHEALPix has found wide application in organizing global geospatial data
[9][10][11]. All the previously mentioned projections are equal-area, but with significant angular distortions and even discontinuities. Due to their simpler implementation and relatively good balance between angular and areal distortion, many hexahedral projections also find application in computer graphics
[12] and, in particular, in planet-scale terrain visualization
[13]. Among the best known are Adjusted Spherical Cube (ASC)
[14], Continuous Cube Mapping (CCM)
[15], and Cartesian Spherical Cube (CSC)
[16].
In addition to polyhedral projections, the UTM and Equi7 grid systems and equidistant cylindrical projection, also known as plate carree, are commonly used today. Equidistant cylindrical projection, despite its considerable distortion, is still widely used to organize global data because it is simple, has a regular grid structure, and can cover the entire world in a single rectangular map. The Universal Transverse Mercator (UTM) grid system, on the other hand, divides Earth into 60 zones, which reduces distortion but limits its application to medium-scale mapping. Equi7
[17] is an example of a compromise that combines seven continental grids based on equidistant azimuthal projections.
The expansion of the system for organizing and referencing global geospatial data occurs at the turn of the century, when the first classifications appear and the term Discrete Global Grid System (DGGS) is introduced for a very significant class of such systems. The importance of DGGSs is also reflected in the fact that the Open Geospatial Consortium (OGC) established the DGGS Standard and Domain Working Groups to support the standardization of these systems. In 2017, the OGC published the first version of the DGGS Abstract Specification
[18]. The standardization process continued, resulting in a formal specification defined by ISO 19170-1:2021 standard
[19] and a revised version of the OGC Abstract Specification
[20]. The popularity of DGGS has also increased due to numerous open source implementations
[21].
A recent trend in geospatial data processing involves the implementation of datacubes based on DGGSs
[22], enabling efficient management of big data workflows. DGGSs, as a standardized representation of Earth, provide the foundational platform for Digital Earth
[23]. Digital Earth is a concept that aims to create an interactive digital replica of the entire planet, fostering a shared understanding of the relationships between the physical and natural environment and society
[24].
The design of GDGGSs is characterized by five fundamental elements
[4]:
- A regular base polyhedron;
- The orientation of the base polyhedron with respect to the planet;
- A hierarchical spatial partitioning of the polyhedron faces;
- The mapping of a spherical or ellipsoidal surface to polyhedral faces and vice versa;
- Methods for indexing and addressing cells.