1. Introduction
With the rapid advancement of space science and technology, coupled with the continuous deepening of global economic integration, the scope of human research has progressively expanded from localized regions to encompass the entire globe. The evolution of modern surveying techniques, including photogrammetry, remote sensing, and GPS, has enabled the acquisition of extensive, multi-scale, and multi-temporal data on a global scale. To effectively store, manage, and analyze spatial data across the globe, there is a necessity for the development of a global, continuous, and multi-scale dynamic data model. Spherical Voronoi diagrams, as a significant data model endowed with beneficial attributes such as natural proximity and dynamic stability, have been widely utilized in critical domains including oceanography, meteorology, airspace management, and satellite planning [
1,
2,
3,
4,
5]. They stand out as one of the most promising approaches for the dynamic management and analysis of global spatial information.
Existing algorithms for generating spherical Voronoi diagrams are often extensions of their planar counterparts. They can be categorized based on the foundational models they employ: algorithms based on vector models and those based on raster models. Both categories of algorithms possesses distinct characteristics, making them suitable for a variety of applications.
Vector algorithms [
6,
7,
8,
9,
10,
11] generally perform well with sets of points on the sphere, but they encounter complex calculations and produce intricate data structures when dealing with sets of lines and surfaces on the sphere. Moreover, the results of vector algorithms are not conducive to the multi-scale representation of spatial data. Additionally, constructing vector Voronoi diagrams on an ellipsoidal surface requires iterative calculations based on approximate coordinates of an auxiliary sphere [
2]. The vector Voronoi diagrams struggle to fundamentally address the challenges of multi-scale integrated processing and the comprehensive analysis of massive global data [
12], which is one of the bottleneck issues limiting their wider application on a larger scale.
In Earth science applications, the generation of Voronoi diagrams based on raster models is of significant importance. First, raster Voronoi diagrams allow for the multi-scale visualization of complex shapes through a set of grids, reducing the complexity associated with constructing Voronoi diagrams for spherical line and surface sets and simplifying the generation process for these geometric elements. Second, the algorithms for generating these diagrams are more easily extended to ellipsoidal surfaces, eliminating the need for additional computational steps. Lastly, grids are pivotal within the digital Earth system, serving as the foundational structure for contemporary digital and intelligent management applications [
13]. A wide variety of data, including remote sensing imagery, are stored and represented in raster format. Raster-based Voronoi diagrams utilize this grid data directly, thus avoiding the complex conversion process from raster to vector formats [
14].
Many scholars have conducted research on the generation of raster-based spherical Voronoi diagrams. In 2003, Zhao et al. [
12] proposed a classic dilation algorithm for constructing spherical raster Voronoi diagrams, which partitions the sphere into a Quaternary Triangular Mesh (QTM) grid and then uses dilation from digital morphology to progressively expand generation points throughout the sphere. However, utilizing the QTM grid distance as the benchmark for dilation results in an approximation of the exact Voronoi diagram on Earth, with errors that are challenging to control [
14]. Furthermore, the dilation process necessitates the manipulation of all neighboring grids, including the repeated selection and subsequent elimination of duplicates, which introduces significant computational overhead and substantially reduces algorithmic efficiency. In 2014, Hu et al. [
14] introduced an ellipsoid-based raster Voronoi diagram generation algorithm that employs an ellipsoidal distance metric during the raster scan distance transformation process. In 2022, Kazemi et al. [
15] leveraged the multi-resolution properties of Discrete Global Grid Systems, employing a hierarchical and progressive approach to implement an efficient spherical distance transformation algorithm. The algorithms developed by Hu and Kazemi significantly improved the construction efficiency of raster Voronoi diagrams and were effective in controlling errors within the confines of a single grid.
However, the algorithms by Hu and Kazemi encounter significant challenges in dynamically constructing Voronoi diagrams. In the myriad of practical applications for Voronoi diagrams, the underlying environments frequently experience dynamic changes. These changes necessitate real-time updates and ongoing maintenance to ensure the accuracy and effectiveness of the Voronoi diagram. Given the dynamic stability characteristic of Voronoi diagrams, only the parts corresponding to the local regions affected by environmental changes undergo alteration, while the majority of the overall diagram’s structure remains unaffected. Utilizing its local mechanism, this targeted process, which addresses changes solely within the affected Voronoi regions and precludes the need for a comprehensive reconstruction of the entire diagram, is referred to as a Voronoi diagram update. In most instances, this approach can significantly reduce computational costs. The algorithm by Hu and Kazemi does not possess a strategy capable of accurately pinpointing the regions that necessitate updates, thereby rendering it unsuitable for direct application in Voronoi diagram updates. As illustrated in
Figure 1,
Figure 1a depicts the Voronoi diagram constructed with generation points A, B, C, and D. Upon the insertion of a new generation point E within the Voronoi region of point D, the resulting updated Voronoi diagram is shown in
Figure 1b. The gray area in this figure represents the Voronoi region corresponding to the new generation point E and is also the region where the Voronoi diagram has undergone change. It can be observed from
Figure 1b that this area overlaps with parts of the Voronoi regions originally associated with generation points A, B, C, and D. This region is neither enclosed by the region formed by the neighboring generation points of point D nor covered by the Voronoi regions corresponding to the neighboring generation points of point D. Consequently, it is not feasible to determine the area needing an update based solely on the topological adjacency of the inserted point.
Canonical ordering, a strategy for the rearrangement of matrices, is utilized to reorganize complex topological structures into a simplified, ordered list or sequence. This method reduces the complexity of processing, facilitating more convenient analysis and visualization. It has been widely applied in various fields, including optimized image processing, real-time pathfinding on regular grids, and distance transformations on planar lattices [
16,
17,
18,
19]. However, due to significant topological differences between spherical and planar spaces, and considering the singularities at the Earth’s poles and the distortions caused by projection [
15], the rules of canonical ordering applicable to planar lattices are not directly transferable to spherical spaces.
This paper introduces a canonical ordering strategy and presents a dynamic construction scheme for spherical Voronoi diagrams based on the Quaternary Triangular Mesh (QTM) system. This scheme is tailored for the generation, dynamic updating, and multi-scale visualization of spherical Voronoi diagrams. Initially, the Earth’s surface is partitioned according to the QTM methodology. Subsequently, a three-level data structure, composed of dictionaries and queues, is designed. Thereafter, an ordered dilation algorithm is devised. By incorporating the canonical ordering strategy and considering the properties of the spherical triangular grid, a procedural sequence for the ordered dilation of the spherical grids is designed. Special treatment methods are also developed to address changes in the topological structure of neighborhoods across the basic faces. A unified and standardized sorting framework is constructed for the grid dilation process, aiming to enhance overall efficiency and performance. Ultimately, the spherical Voronoi diagram is constructed and updated through the ordered dilation of points, lines, and surfaces on the sphere. The delineation of the update region’s boundary is determined by the grid dilation. Additionally, the multi-resolution attribute of the grid enables the multi-scale visualization of the global Voronoi diagram.
This paper is organized into five sections, beginning with an introduction.
Section 2 details the specific steps of the ordered dilation algorithm proposed for the generation and rapid local updating of spherical Voronoi diagrams.
Section 3 presents a comparative analysis of the algorithm’s feasibility and efficiency, complemented by an application demonstration.
Section 4 offers discussions on the findings.
Section 5 concludes the paper.
3. Experiment
The programming language utilized for the algorithm presented in this paper is Python (3.9.0). The three-dimensional visualization is accomplished using the Matplotlib toolkit. Calculations for spherical and ellipsoidal distances are facilitated by the PYPROJ library. The hardware environment consists of an 11th Gen Intel(R) Core(TM) i5-1135G7 @ 2.40 GHz, with 24.0 GB of RAM.
3.1. Feasibility and Completeness of the Algorithm
This paper conducts an experimental validation of the feasibility of the proposed method.
Figure 10a–c present illustrations of the results generated by the ordered expansion algorithm under various grid levels and different numbers of generation points.
Figure 10d,e showcase illustrations of the spherical Voronoi diagrams constructed by the algorithm, utilizing the major landmasses as the generation surface.
An experiment is designed to measure the error by directly applying the definition of Voronoi, which involves randomly sampling points on the boundary lines and calculating the absolute value of the distance difference between the sampled points and the two sides of the Voronoi generator. As shown in
Figure 11, the boundary points, represented by red dots, are randomly located on the shared edges of adjacent grid cells. The yellow dots represent the Voronoi generators. The spherical distances between these points and their respective Voronoi generators, indicated by black lines, are used to calculate the Voronoi error, which is the absolute difference between these two distances.
Equation (1) is used to calculate the algorithm error results. In the equation, the function
is utilized to compute the spherical distance between the generator and the boundary point,
and
represent the two sides of the Voronoi generator,
means the randomly sample point on the boundary line, and
denotes the maximum edge length of the grid at the current hierarchical level.
Due to the close relationship between error magnitude and grid cell size,
is defined as the ratio of the error to the edge length of the grid cell. Since the grid cells at the same level have slight variations in size and edge length, we select the maximum edge length of the grid cells as the evaluation standard. It should be noted that the maximum grid edge length at each level is approximately 1.5 times the average edge length. The experimental results of the algorithm’s error are depicted in
Table 1, with the 6th, 7th, 8th, 10th, and 14th levels of QTM serving as the foundations.
As shown in the results, the error for most boundary random points falls within the range of 0 to 0.3. The error of all boundary points is confined within one grid cell of the corresponding hierarchical level, i.e., the error is less than 280,006 m at level 6 and less than 18,086 m at level 10. By increasing the subdivision level locally, according to the specific needs, the error can be controlled more precisely.
The algorithm by Kazemi et al. [
15] has demonstrated its accuracy through experimental validation. The spherical distance transformation produced by their method exhibits only minor errors when compared to the results generated by ArcGIS Pro, which can be disregarded. This study reproduces their algorithm and employs the same spherical distance calculation tool used in our algorithm. The spherical results of their algorithm are utilized as a benchmark for the accuracy assessment of our algorithm’s results.
This paper employs Equation (2) to calculate the error of the algorithm’s results. In the equation, the function
Dist(
a,
b) is utilized to compute the spherical distance between the site and the grid, and
and
represent the closest growth points calculated using the algorithm presented in this paper and Kazemi’s algorithm, respectively.
L denotes the average edge length of the grid at the current hierarchical level. The experimental results of the algorithm’s error are depicted in
Table 2. The experiments are based on a set of 100 sites that are randomly distributed across the surface of a sphere, with the 7th, 8th, 9th, and 10th levels of QTM serving as the foundations.
As can be observed, using Kazemi’s algorithm as a benchmark, the Voronoi diagrams generated by the algorithm presented in this paper have controlled the error within one grid. Given that raster data inherently possess an error equivalent to one grid cell in expressing geospatial entities, the algorithm proposed is considered to have achieved a comparable level of accuracy to Kazemi’s algorithm. This affirms the completeness of our algorithm.
3.2. Efficiency Evaluation
This paper reproduces the algorithms of Kazemi [
15], Zhao [
12], and Hu [
14] and conducts efficiency comparison experiments with them. The reproduced algorithms of others adopted the same distance metric as this paper and utilized the same spherical toolbox. The sites consist of a set of 100 randomly distributed spherical points. The experimental results of the overall construction of the spherical Voronoi diagram are presented in
Table 3 and illustrated in
Figure 12.
Figure 12 represents the time required by the algorithm, while the horizontal axis delineates the levels of the QTM. It should be noted that Hu’s algorithm is predicated on a planar grid; hence, the time recorded for Hu’s algorithm is under an equivalent number of grids. The experimental results demonstrate that, in the context of constructing a complete spherical Voronoi diagram, the efficiency of the ordered dilation algorithm proposed in this paper is significantly higher than that of Zhao’s algorithm and slightly superior to the algorithms of Kazemi and Hu. Overall, the efficiency of our algorithm is broadly considered to be on par with the algorithms of Kazemi and Hu in terms of overall efficiency.
Similarly, for scenarios where changes occur in local environments, the algorithm presented in this paper is capable of selectively updating the affected regions without the need for a complete reconstruction of the Voronoi diagram. An efficiency comparison experiment was conducted, based on a nine-level QTM, consisting of 2,097,152 grids, with 100 random spherical points serving as the sites. This paper has designed four distinct scenarios: In the first scenario, all sites undergo changes, necessitating the construction of an entirely new spherical Voronoi diagram; in the second scenario, the coordinates of 20% of the sites are altered, which includes operations such as insertions, movements, or deletions; in the third scenario, the coordinates of 10% of the sites are changed; and in the fourth scenario, the change affects 1% of the sites.
Figure 13 presents the time required by the algorithms of Hu and Kazemi and the algorithm proposed in this paper for the first scenario, as well as the time consumed by our algorithm for the second, third, and fourth scenarios. In the second, third, and fourth scenarios, the algorithms of Hu and Kazemi are not readily amenable to localized updates, resulting in an invariant expenditure of time. The experimental results substantiate that our proposed update algorithm can effectively enhance efficiency when changes occur in localized environments.
From the perspective of time complexity, when reconstructing the entire spherical Voronoi diagram, the algorithm presented in this paper, as well as the algorithms by Hu and Kazemi, exhibit a time complexity of O(n), where n is the number of spherical grids. In the context of dynamic updates, our algorithm’s time complexity is O(m), where m is the number of grids requiring updates. Should the area experiencing dynamic changes be considerably smaller than the entire domain, that is, when m is much less than n ((m << n), the efficiency can be significantly enhanced. The current best vector algorithms have time complexity O(n), where n depends on the number and shape of spatial objects. Given the distinct logical constructs inherent in vector and raster algorithms, a straightforward comparison based on time complexity alone is not feasible.
3.3. Application Examples of Dynamic Navigation and Management for Ocean Vessels
Voronoi diagrams are extensively utilized in the fields of maritime navigation and vessel management [
28,
29,
30]. By employing Voronoi diagrams, a series of distance-based boundary lines can be generated in maritime areas. This approach facilitates the analysis of ship distribution and aids in the planning and optimization of shipping routes. It ensures safe distances between vessels and coordinates their direction of travel to prevent congestion and navigate around obstacles, such as land, reefs, and other ships. Furthermore, the analysis of Voronoi diagrams enables rapid emergency response and rescue operations and allows for the more effective management of maritime resources, including the identification of the nearest rescue points, oil platforms, and fishing resources.
Taking into account the advancing automation in maritime navigation and the imperatives for safety and real-time capabilities within practical applications, this application calls for an algorithm that is not only highly precise but also efficient. To achieve this, it is essential to integrate a variety of data types and sources, including weather and vessel information, within a cohesive framework. Furthermore, the application necessitates the multi-scale construction of Voronoi diagrams and the real-time, dynamic updating of Voronoi diagrams.
3.3.1. Dynamic Construction of Spherical Voronoi Diagrams and Ship Routing
This paper employs a multi-level Quadrilateralized Triangular Mesh (QTM) as a framework for global data integration. Utilizing the sixth level of the QTM, with an approximate edge length of 200 km, the study incorporates authentic ocean storm wave height data sourced from the National Oceanic and Atmospheric Administration of the United States. Oceanic regions with wave heights exceeding 5 m are delineated by orange grids, while areas with waves over 9 m are represented by red grids, as illustrated in
Figure 14a.
Subsequently, the spherical Voronoi diagram is constructed using the land and storm regions as the generating surfaces, as depicted in
Figure 14b. Thereafter, the A* algorithm is employed to derive safe routes that circumvent both land and storm areas, followed by the application of a polynomial interpolation method to smooth these paths, with the results presented in
Figure 14c. Maritime weather conditions are subject to change, with significant alterations in the oceanic wind and wave hazard areas, as shown in
Figure 14d, necessitating adjustments to the navigation paths of ocean vessels. Utilizing the spherical Voronoi diagram update algorithm presented in this paper, the Voronoi diagram is updated, with the outcomes displayed in
Figure 14e. Subsequently, the ship’s path is adjusted accordingly, as depicted in
Figure 14f.
Figure 15 presents a comparative analysis of the number of grids accessed and the time required for the updates of spherical Voronoi diagrams over time between the algorithm introduced in this paper and Kazemi’s algorithm. Kazemi’s algorithm involves a complete reconstruction of the Voronoi diagram, resulting in an invariant number of grid accesses and a nearly constant time requirement. In contrast, our algorithm accesses a significantly smaller number of grids within the updated regions. This indicates that our algorithm offers improved efficiency under dynamic conditions.
3.3.2. Global Vessel Management Based on Multi-Scale Voronoi Diagrams
This paper employs the 10th level of the QTM, with edges roughly 10 km in length, to integrate over 9000 global fishing vessel coordinates from Global Fishing Watch, constructing a global Voronoi diagram, as shown in
Figure 16a. The Voronoi diagram delineates fishing resource areas for different vessels. Over time, the coordinates of approximately 500 vessels shifted, prompting updates to the Voronoi diagram. The outcomes of these updates are presented in
Figure 16b, with the Voronoi regions corresponding to the altered vessels highlighted in green.
Within the context of a localized area, this study selects the Bohai Bay as the experimental region, based on the 14th level of the QTM, where the edge length of the grids is approximately 0.5 km.
Figure 17 displays a grid-based Voronoi diagram constructed on a local scale within a globally unified framework, utilizing vessels and peninsulas as the generating elements. The positions of ships at sea change over time, and the updates to the Voronoi diagram are reflected in the results presented in
Figure 17. An enlarged view on the right side of the figure highlights the changes in the grid-based Voronoi diagram following these updates.
In our global maritime dynamics management experiment, we constructed Voronoi diagrams at various scales. The environment incorporated static elements like landmasses and some ships, alongside dynamic factors such as storms and moving vessels, causing localized Voronoi diagram alterations over time. Our algorithm performs these updates locally rather than resorting to full reconstruction.
Table 4 contrasts our approach with Kazemi’s algorithm in terms of the average update time and grid numbers across scales, demonstrating the superior efficiency of our method in the dynamic experiment.
It is noteworthy that beyond the application to ocean vessel management as covered in the experiments of this paper, spherical Voronoi diagram structures have found extensive application in significant domains such as airspace management and satellite planning. The algorithm presented in this paper has the potential to be extended to these fields, demonstrating a broad spectrum of application prospects.
4. Discussion
4.1. Spherical and Ellipsoidal Spatial Reference
Due to the distortions caused by projection and significant topological differences between the spherical and planar spaces [
15], existing algorithms for updating planar space Voronoi diagrams [
19,
31,
32,
33,
34,
35] are not readily applicable to spherical space. Kazemi et al. [
15] have experimentally demonstrated the superior accuracy of DGGS in performing spherical distance transformations compared to planar algorithms. Additionally, planar algorithms face challenges in processing at the Earth’s poles, as the singularity of the topological relationship at the poles also makes it difficult for planar algorithms to adapt. The Discrete Global Grid System (DGGS) is a global grid model based on the true geospatial space, which fits the Earth’s surface by recursively subdividing the spherical surface, thereby avoiding the various deformations and data fractures caused by planar projection. This paper, grounded in a spherical spatial reference, presents the design of Voronoi diagram generation and updating algorithms that are tailored to the characteristics of spherical grids, effectively resolving the aforementioned issues.
In the realms of geography and Earth science, to more effectively cater to the geospatial demands of the real world, numerous applications favor the use of a geodetic ellipsoid model—such as the WGS84 ellipsoid—over a basic spherical model as their foundational framework. These applications employ geodetic ellipsoidal distances as the standard measure, rather than spherical distances. The construction of vector Voronoi diagrams on an ellipsoidal surface is a complex task that necessitates iterative calculations based on the approximate coordinates derived from an auxiliary sphere [
2]. For instance, Kastrisios et al. [
6] constructed vector Voronoi diagrams on an ellipsoidal surface, where 17% of the vertices required iteration in excess of a hundred times.
The algorithm presented in this paper is applicable to both spherical and ellipsoidal spaces, contingent upon the underlying QTM framework employed. It may be predicated on a generic spherical model or aligned with an ellipsoidal model, with the distance metric selected being either spherical or ellipsoidal. The foundational triangles of the spherical QTM are constituted by spherical triangles with their vertices situated on the spherical surface. Conversely, the foundational triangles of the ellipsoidal QTM consist of ellipsoidal triangles with vertices positioned on the surface of the reference ellipsoid. This paper delineates the spherical distance benchmark and the ellipsoidal distance benchmark as the great-circle arc distance between the centroids of the spherical triangle grids or the ellipsoidal distance, respectively. For the computation of ellipsoidal distances, this paper adopts the WGS84 reference ellipsoid as the basis for the ellipsoidal model. The Python PYPROJ package, renowned for its robust ellipsoidal computation functions [
36], is utilized, which is predicated on the classical Jacobi ellipsoidal functions; the ellipsoid is defined as
the ellipsoidal distance is defined as
Herein, (
β,
ω) represents the ellipsoidal latitude and longitude, and
represents the path azimuth angle:
This paper leverages the PYPROJ toolkit to test the efficiency of distance calculations on spherical and ellipsoidal surfaces. The left axis of
Figure 18 shows the time needed to compute distances 100,000 times. Subsequently, based on the eighth-level QTM tessellation, which includes approximately over half a million cells, Voronoi diagrams are constructed using both spherical and ellipsoidal distance metrics. The results are presented on the right axis of
Figure 18. The experimental results indicate that the construction time for the ellipsoidal Voronoi diagram is slightly higher than that for the spherical Voronoi diagram, but it is still within an acceptable range.
This paper defines the distance metric as the spherical or ellipsoidal distance between the centroids of spherical triangular grids. The spherical distance is calculated based on the great-circle arc length, while the ellipsoidal distance is measured by the shortest path length over the reference ellipsoidal surface. The distance measurement within the grid space is not entirely equivalent to the point-to-point distance in vector space, and this process inevitably introduces some error. The results of the algorithm presented in this paper limit the error to within a single grid cell. The grid of QTM possesses attributes of multi-resolution and infinite subdivision. The average unit area and average edge length of the grids at different levels of the QTM, based on ellipsoidal and spherical benchmarks, are shown in
Table 5. In practical applications, one can select the grid scale based on specific requirements, balancing the needs for efficiency and accuracy to control the error. For instance, in the experiments of
Section 3.3, a global-scale route is approximated for ocean-going vessels, utilizing the sixth-level grid to establish a 20 km navigational interval for vessel transit. When allocating fishing resources for marine fishing boats, the 14th-level grid is chosen, with a grid edge length of approximately 500 m and a boundary error within 500 m. For detailed navigation and obstacle avoidance operations of the vessels, a higher-resolution grid can be delineated in the local area. If the requirement is to control the error within 10 m, the 20th-level grid can be selected.
4.2. Voronoi Diagram Construction Based on Multi-Resolution Grids
In global spatial applications, there is a significant disparity in the scales of various applications and data. Spherical vector Voronoi diagrams face certain challenges when applied across multiple scales. This is due to the fact that vector data are typically stored as a collection of geometric elements such as points, lines, and polygons, each with specific positional and attribute information. When attempting to represent spatial information at different scales, the vector data structure may need to adapt by either subdividing or simplifying geometric elements, which can lead to an increase in data complexity and storage volume. This also augments the complexity of data processing and analysis.
In comparison, the QTM allows for infinite recursive subdivision to closely fit the Earth’s surface, with grids of varying resolutions nested within one another. The Voronoi diagrams constructed in this paper based on the QTM are inherently capable of adapting to geographic data from different sources and with different resolutions. The multi-level nature of the QTM provides a unified framework for Voronoi diagrams at various scales, enabling the integration and interaction of heterogeneous data and Voronoi diagrams with different resolutions within the same framework. This approach enhances the consistency and operability of the algorithms.
Taking
Section 3.3 on ocean-going vessel navigation and management as an example, our approach balances efficiency and accuracy to generate Voronoi diagrams of varying resolutions based on application context. At a macro scale, we utilize continents and storm areas as generators to create high-resolution Voronoi diagrams for overall direction and route planning. In contrast, micro-scale local areas are addressed with finer resolution for more detailed boundary delineation.
In the research presented in this paper, the transformation between Voronoi diagrams of different resolutions is achieved through the subdivision and merging of the Voronoi diagram’s boundary grids, as depicted in
Figure 19. Specifically, when transitioning from a high-resolution Voronoi diagram to a lower-resolution one, the process begins with the subdivision of the boundary grids, yielding four smaller triangles. The subsequent goal is to determine the closest feature point to each of these four smaller triangles individually. Given that grids within the same Voronoi region are generally continuous, the four smaller triangles can each select the closest feature point from the nearest feature points of their parent grid and the neighboring grids of the parent grid. After calculating the distances to each of these points individually, the new boundary of the Voronoi region is established. Subsequently, the aforementioned process is repeated until the grid resolution of the Voronoi region boundary meets the desired target resolution. Conversely, the conversion from a low-resolution grid Voronoi diagram to a high-resolution grid is relatively straightforward. This is because the nearest feature point for the high-resolution grid can be derived from the feature point of the most central small triangle. This operation can be accomplished based on the hierarchical indexing of the grids.
This study has achieved interconversion between Voronoi diagrams of different resolutions by the hierarchical characteristics of the QTM. However, the research has not fully utilized the multi-level features within the process of the Voronoi diagram update algorithm. Specifically, in this study, when dealing with high-resolution and low-resolution Voronoi diagrams, updates were performed independently for each, without fully integrating the update process for multi-scale Voronoi diagrams at the algorithmic level. This is due to the fact that we must consider the update frequency and the update areas of Voronoi diagrams at different scales. If we aim to establish a unified update process for multi-resolution Voronoi diagrams, we should fully consider the relationships between Voronoi diagrams of various scales, including the selection of the scale and direction for the dilation. For example, during an ordered dilation of the grid, one might choose to dilate by expanding 10 units of the 6th-level grid, 2 units of the 8th-level grid, and 1 unit of the 10th-level grid in the northeast direction. This approach helps to reduce computational redundancy and the need for redundant judgment when converting between different resolution Voronoi diagrams, thereby achieving consistency and efficiency in the global update of Voronoi diagrams. It allows the algorithm to be more flexible in adapting to the needs of different scales, which may potentially enhance the efficiency of the algorithm. We express a keen interest in further exploring this aspect and plan to focus on resolving this issue in our subsequent research.
4.3. Interoperability between Different Grids in Voronoi Diagrams
Global discrete grids are available in various fundamental shapes, with the shape of grid cells determining the storage, indexing, and relational characteristics of the grid system. Each type of partitioned grid possesses distinct structural advantages. A triangular grid offers high rendering efficiency and more accurate boundary fitting. A diamond-shaped grid can leverage the rich algorithms of planar grids. A hexagonal grid benefits from adjacency consistency. Since different grids are suitable for different application scenarios, to meet the application needs of various fields and to achieve the optimal expression, processing, and analysis of various types of global spatiotemporal data, there is a requirement for transformation between different grid models.
The global discrete grid system is constituted by “grid points”, “grid edges,” and “grid cells”. Grid points act as the central control elements of the DGGS, usually arranged at certain intervals across the Earth’s surface. Their density is contingent upon the selected grid system and its corresponding level of partitioning. Grid edges represent the boundaries or connectors between adjacent grid points, delineating the regional boundaries and the structural framework of the grid system. Grid cells are polygons formed by the grid edges surrounding adjacent grid points, with each grid cell signifying a spatial unit within the global discrete grid system. Grid points that share the same spatial distribution can create various grid cell systems through distinct arrangements of grid edges. There is potential to investigate a unified expression pattern from grid points to different grid cells.
Based on this concept, we can develop a set of spherical Voronoi diagram construction and update models centered around grid points, independent of the shape of the grid cells. Specifically, using the centroids of octahedral QTM as grid points, we can construct isomorphic diamond-shaped and hexagonal grids. Each triangular grid encompasses a single grid point, while a diamond-shaped grid, which is composed of two triangular grids, encompasses two grid points. Furthermore, a hexagonal grid, consisting of six triangular grids, encompasses six grid points, as depicted in
Figure 20a. Utilizing this unified model centered on grid points enables facile transformation between grid cells of diverse shapes.
For raster spherical Voronoi diagrams, the Voronoi region attribute of each grid is determined by the proximity of the central point to the nearest feature point. The algorithm presented in this paper is capable of determining the nearest feature point for each grid point within a triangular grid. In diamond-shaped Voronoi diagrams, the grids are defined by the two internal grid points, whereas in hexagonal Voronoi diagrams, the grids are determined by the six internal grid points, as depicted in
Figure 20b. Utilizing this method, the transformation between grids of different shapes based on the grid point system can be easily accomplished, thus meeting the needs of various application scenarios. It should be noted that this conversion process may lead to some information loss.
In conclusion, the construction of spherical Voronoi diagrams based on grid points offers the flexibility to adapt to various grid shapes and application backgrounds, along with high scalability and universality. Moreover, the mobility of grid points allows for a better fit to the data, thereby enhancing the accuracy of the model. As illustrated in
Figure 21, grid points can be aggregated into areas of interest without altering the topological structure of the grid, which does not impact the application of the Voronoi diagram algorithm presented in this paper. This aggregation results in a denser Voronoi diagram in specific regions, effectively increasing the local resolution. Consequently, this approach enables a clearer acquisition of information from focal areas, facilitating a more in-depth analysis of the characteristics and variations within those regions.
The perturbation of grid points enables the construction of spherical Voronoi diagrams with variable resolutions. The unaltered topological relationships of the grid ensure the usability of the algorithm presented in this paper. However, the acceptable degree of perturbation requires further research. Perturbation can alter the distribution of grid points to a certain extent, and it is essential to ensure that such alterations do not adversely affect the accuracy and stability of the Voronoi diagram algorithm to an unacceptable degree. Therefore, additional experimental testing is necessary to determine the acceptable range of perturbation. Moreover, the multi-level of the grid point system is also worth exploring, implementing hierarchical subdivision based on the perturbed grid points and constructing spherical Voronoi diagrams on this foundation. This provides new perspectives and methodologies for the modeling and updating of Voronoi diagrams based on a global discrete grid.