stringtranslate.com

Discrete global grid

A discrete global grid (DGG) is a mosaic that covers the entire Earth's surface. Mathematically it is a space partitioning: it consists of a set of non-empty regions that form a partition of the Earth's surface.[1] In a usual grid-modeling strategy, to simplify position calculations, each region is represented by a point, abstracting the grid as a set of region-points. Each region or region-point in the grid is called a cell.

When each cell of a grid is subject to a recursive partition, resulting in a "series of discrete global grids with progressively finer resolution",[2] forming a hierarchical grid, it is called a hierarchical DGG (sometimes "global hierarchical tessellation"[3]or "DGG system").

Discrete global grids are used as the geometric basis for the building of geospatial data structures. Each cell is related with data objects or values, or (in the hierarchical case) may be associated with other cells. DGGs have been proposed for use in a wide range of geospatial applications, including vector and raster location representation, data fusion, and spatial databases.[1]

The most usual grids are for horizontal position representation, using a standard datum, like WGS84. In this context, it is common also to use a specific DGG as foundation for geocoding standardization.

In the context of a spatial index, a DGG can assign unique identifiers to each grid cell, using it for spatial indexing purposes, in geodatabases or for geocoding.

Reference model of the globe

The "globe", in the DGG concept, has no strict semantics, but in geodesy a so-called "grid reference system" is a grid that divides space with precise positions relative to a datum, that is an approximated a "standard model of the Geoid". So, in the role of Geoid, the "globe" covered by a DGG can be any of the following objects:

As a global modeling process, modern DGGs, when including projection process, tend to avoid surfaces like cylinder or a conic solids that result in discontinuities and indexing problems. Regular polyhedra and other topological equivalents of sphere led to the most promising known options to be covered by DGGs,[1] because "spherical projections preserve the correct topology of the Earth – there are no singularities or discontinuities to deal with".[4]

When working with a DGG it is important to specify which of these options was adopted. So, the characterization of the reference model of the globe of a DGG can be summarized by:

NOTE: when the DGG is covering a projection surface, in a context of data provenance, the metadata about reference-Geoid is also important — typically informing its ISO 19111's CRS value, with no confusion with the projection surface.

Types

The main distinguishing feature to classify or compare DGGs is the use or not of hierarchical grid structures:

Other usual criteria to classify a DGG are tile-shape and granularity (grid resolution):

Examples

Non-hierarchical grids

The most common class of discrete global grids are those that place cell center points on longitude/latitude meridians and parallels, or which use the longitude/latitude meridians and parallels to form the boundaries of rectangular cells. Examples of such grids, all based on latitude/longitude:

Hierarchical grids

Successive space partitioning. The grey-and-green grid in the second and third maps are hierarchical.

The right aside illustration show 3 boundary maps of the coast of Great Britain. The first map was covered by a grid-level-0 with 150 km size cells. Only a grey cell in the center, with no need of zoom for detail, remains level-0; all other cells of the second map was partitioned into four-cells-grid (grid-level-1), each with 75 km. In the third map 12 cells level-1 remains as grey, all other was partitioned again, each level-1-cell transformed into a level-2-grid.
Examples of DGGs that use such recursive process, generating hierarchical grids, include:

Standard equal-area hierarchical grids

There is a class of hierarchical DGG's named by the Open Geospatial Consortium (OGC) as "discrete global grid systems" (DGGS), that must to satisfy 18 requirements. Among them, what best distinguishes this class from other hierarchical DGGs, is the Requirement-8, "For each successive level of grid refinement, and for each cell geometry, (...) Cells that are equal area (...) within the specified level of precision".[24]

A DGGS is designed as a framework for information as distinct from conventional coordinate reference systems originally designed for navigation. For a grid-based global spatial information framework to operate effectively as an analytical system it should be constructed using cells that represent the surface of the Earth uniformly.[24] The DGGS standard include in its requirements a set of functions and operations that the framework must to offer.

All DGGS's level-0 cells are equal area faces of a Regular polyhedra...

Regular polyhedra (top) and their corresponding equal area DGG
Regular polyhedra (top) and their corresponding equal area DGG

The DGGS framework

The standard defines the requirements of a hierarchical DGG, including how to operate the grid. Any DGG that satisfies these requirements can be named DGGS. "A DGGS specification SHALL include a DGGS Reference Frame and the associated Functional Algorithms as defined by the DGGS Core Conceptual Data Model".[25]

For an Earth grid system to be compliant with this Abstract Specification it must define a hierarchical tessellation of equal area cells that both partition the entire Earth at multiple levels of granularity and provide a global spatial reference frame. The system must also include encoding methods to: address each cell; assign quantized data to cells; and perform algebraic operations on the cells and the data assigned to them. Main concepts of the DGGS Core Conceptual Data Model:
  1. reference frame elements, and,
  2. functional algorithm elements; comprising:
    1. quantization operations,
    2. algebraic operations, and
    3. interoperability operations.

Database modeling

In all DGG databases the grid is a composition of its cells. The region and centralPoint are illustrated as typical properties or subclasses. The cell identifier (cell ID) is also an important property, used as internal index and/or as public label of the cell (instead the point-coordinates) in geocoding applications. Sometimes, as in the MGRS grid, the coordinates make the role of ID.

There are many DGGs because there are many representational, optimization and modeling alternatives. All DGG grid is a composition of its cells, and, in the Hierarchical DGG each cell uses a new grid over its local region.

The illustration is not adequate to TIN DEM cases and similar "raw data" structures, where the database not use the cell concept (that geometrically is the triangular region), but nodes and edges: each node is an elevation and each edge is the distance between two nodes.

In general, each cell of the DGG is identified by the coordinates of its region-point (illustrated as the centralPoint of a database representation). It is also possible, with loss of functionality, to use a "free identifier", that is, any unique number or unique symbolic label per cell, the cell ID. The ID is usually used as spatial index (such as internal Quadtree or k-d tree), but is also possible to transform ID into a human-readable label for geocoding applications.

Modern databases (e.g. using S2 grid) use also multiple representations for the same data, offering both, a grid (or cell region) based in the Geoid and a grid-based in the projection.

History

Discrete global grids with cell regions defined by parallels and meridians of latitude/longitude have been used since the earliest days of global geospatial computing. Before it, the discretization of continuous coordinates for practical purposes, with paper maps, occurred only with low granularity. Perhaps the most representative and main example of DGG of this pre-digital era was the 1940s military UTM DGGs, with finer granulated cell identification for geocoding purposes. Similarly some hierarchical grid exists before geospatial computing, but only in coarse granulation.

A global surface is not required for use on daily geographical maps, and the memory was very expensive before the 2000s, to put all planetary data into the same computer. The first digital global grids were used for data processing of the satellite images and global (climatic and oceanographic) fluid dynamics modeling.

The first published references to hierarchical geodesic DGG systems are to systems developed for atmospheric modeling and published in 1968. These systems have hexagonal cell regions created on the surface of a spherical icosahedron.[26][27]

The spatial hierarchical grids were subject to more intensive studies in the 1980s,[28] when main structures, as Quadtree, were adapted in image indexing and databases.

While specific instances of these grids have been in use for decades, the term discrete global grids was coined by researchers at Oregon State University in 1997[2] to describe the class of all such entities.

... OGC standardization in 2017...

Comparison and evolution

Comparing grid-cell identifier schemas of two different curves, Morton and Hilbert. The Hilbert curve was adopted in DGGs like S2-geometry, Morton curve in DGGs like Geohash. The adoption of Hilbert curve was an evolution because have less "jumps", preserving nearest cells as neighbours.

The evaluation discrete global grid consists of many aspects, including area, shape, compactness, etc. Evaluation methods for map projection, such as Tissot's indicatrix, are also suitable for evaluating map projection-based discrete global grid.

In addition, averaged ratio between complementary profiles (AveRaComp) [29] gives a good evaluation of shape distortions for quadrilateral-shaped discrete global grid.

Database development-choices and adaptations are oriented by practical demands for greater performance, reliability or precision. The best choices are being selected and adapted to necessities, propitiating the evolution of the DGG architectures. Examples of this evolution process: from non-hierarchical to hierarchical DGGs; from the use of Z-curve indexes (a naive algorithm based in digits-interlacing), used by Geohash, to Hilbert-curve indexes, used in modern optimizations, like S2.

Geocode variants

In general each cell of the grid is identified by the coordinates of its region-point, but it is also possible to simplify the coordinate syntax and semantics, to obtain an identifier, as in a classic alphanumeric grids — and find the coordinates of a region-point from its identifier. Small and fast coordinate representations is a goal in the cell-ID implementations, for any DGG solutions.

There is no loss of functionality when using a "free identifier" instead of a coordinate, that is, any unique number (or unique symbolic label) per region-point, the cell ID. So, to transform a coordinate into a human-readable label, and/or compressing the length of the label, is an additional step in the grid representation. This representation is named geocode.

Some popular "global place codes" as ISO 3166-1 alpha-2 for administrative regions or Longhurst code for ecological regions of the globe, are partial in globe's coverage. By other hand, any set of cell-identifiers of a specific DGG can be used as "full-coverage place codes". Each different set of IDs, when used as a standard for data interchange purposes, are named "geocoding system".

There are many ways to represent the value of a cell identifier (cell-ID) of a grid: structured or monolithic, binary or not, human-readable or not. Supposing a map feature, like the Singapore's Merlion fountaine (~5m scale feature), represented by its minimum bounding cell or a center-point-cell, the cell ID will be:

All these geocodes represents the same position in the globe, with similar precision, but differ in string-length, separators-use and alphabet (non-separator characters). In some cases the "original DGG" representation can be used. The variants are minor changes, affecting only final representation, for example the base of the numeric representation, or interlacing parts of the structured into only one number or code representation. The most popular variants are used for geocoding applications.

Alphanumeric global grids

DGGs and its variants, with human-readable cell-identifiers, has been used as de facto standard for alphanumeric grids. It is not limited to alphanumeric symbols, but "alphanumeric" is the most usual term.

Geocodes are notations for locations, and in a DGG context, notations to express grid cell IDs. There are a continuous evolution in digital standards and DGGs, so a continuous change in the popularity of each geocoding convention in the last years. Broader adoption also depends on country's government adoption, use in popular mapping platforms, and many other factors.

Examples used in the following list are about "minor grid cell" containing the Washington obelisk, 38° 53 22.11″ N, 77° 2 6.88″ W.

Other documented systems:

See also

References

  1. ^ a b c d Sahr, Kevin; White, Denis; Kimerling, A.J. (2003). "Geodesic discrete global grid systems" (PDF). Cartography and Geographic Information Science. 30 (2): 121–134. doi:10.1559/152304003100011090. S2CID 16549179.
  2. ^ a b Sahr, Kevin; White, Denis; Kimerling, A.J. (18 March 1997), "A Proposed Criteria for Evaluating Discrete Global Grids", Draft Technical Report, Corvallis, Oregon: Oregon State University, archived from the original on 3 March 2016, retrieved 10 September 2014
  3. ^ Geoffrey Dutton. "What's the big deal about global hierarchical tessellation?". quote: "a few prototype systems that are either hierarchical tessellations, global tessellations, or both".
  4. ^ "Overview".
  5. ^ "Global 30 Arc-Second Elevation (GTOPO30)". USGS. Archived from the original on July 10, 2017. Retrieved October 8, 2015.
  6. ^ "Global Multi-resolution Terrain Elevation Data 2010 (GMTED2010)". USGS. Retrieved October 8, 2015.
  7. ^ Arakawa, Akio; Lamb, Vivian R. (1977). "Computational Design of the Basic Dynamical Processes of the UCLA General Circulation Model". General Circulation Models of the Atmosphere. Methods in Computational Physics: Advances in Research and Applications. Vol. 17. pp. 173–265. doi:10.1016/B978-0-12-460817-7.50009-4. ISBN 978-0-12-460817-7.
  8. ^ "Research Institute for World Grid Squares". Yokohama City University. Retrieved November 21, 2017.
  9. ^ Snyder, J.P. (1992). "An equal-area map projection for polyhedral globes". Cartographica. 29 (1): 10–21. doi:10.3138/27h7-8k88-4882-1752.
  10. ^ Barnes, Richard (2019). "Optimal orientations of discrete global grids and the Poles of Inaccessibility". International Journal of Digital Earth. 13 (7): 803–816. doi:10.1080/17538947.2019.1576786. S2CID 134622203.
  11. ^ Martin Suess, Pedro Matos (2004). "Processing of SMOS level 1C data onto a discrete global grid". IEEE International IEEE International IEEE International Geoscience and Remote Sensing Symposium, 2004. IGARSS '04. Proceedings. 2004. Vol. 3. pp. 1914–1917. doi:10.1109/IGARSS.2004.1370716. ISBN 978-0-7803-8742-3. S2CID 7267254.
  12. ^ "Global Grid Systems Insight". Global Grid Systems. Archived from the original on November 16, 2020. Retrieved October 8, 2015.
  13. ^ "LAMBDA – COBE Quadrilateralized Spherical Cube".
  14. ^ Dutton, Geoffrey H. (1999). A Hierarchical Coordinate System for Geoprocessing and Cartography. Lecture Notes in Earth Sciences. Vol. 79. doi:10.1007/BFb0011617. ISBN 978-3-540-64980-9.[page needed]
  15. ^ Snyder, Van (10 August 2016). "Quaternary triangular mesh: a hierarchical coordinate system for geoprocessing and cartography". hdl:2014/46916. {{cite journal}}: Cite journal requires |journal= (help)
  16. ^ "HEALPix Background Purpose". NASA Jet Propulsion Laboratory. Retrieved October 8, 2015..
  17. ^ "S2 Cells". Archived from the original on 2020-10-20. Retrieved 2018-09-06.
  18. ^ "HTM Overview".
  19. ^ "ADASS 2003 Conference Proceedings".
  20. ^ "S2 Geometry".
  21. ^ "Overview". s2geometry.io. Retrieved 2018-05-11.
  22. ^ Kreiss, Sven (2016-07-27). "S2 cells and space-filling curves: Keys to building better digital map tools for cities". Medium. Retrieved 2018-05-11.
  23. ^ "S2 Cell Statistics".
  24. ^ a b Open Geospatial Consortium (2017), "Topic 21: Discrete Global Grid Systems Abstract Specification". Document 15-104r5 version 1.0.
  25. ^ Section 6.1, "DGGS Core Data Model Overview", of the DGGS standard
  26. ^ Sadourny, R.; Arakawa, A.; Mintz, Y. (1968). "Integration of the nondivergent barotropic vorticity equation with an icosahedral-hexagonal grid for the sphere". Monthly Weather Review. 96 (6): 351–356. Bibcode:1968MWRv...96..351S. CiteSeerX 10.1.1.395.2717. doi:10.1175/1520-0493(1968)096<0351:iotnbv>2.0.co;2.
  27. ^ Williamson, D.L. (1968). "Integration of the barotropic vorticity equation on a spherical geodesic grid". Tellus. 20 (4): 642–653. doi:10.1111/j.2153-3490.1968.tb00406.x.
  28. ^ Kleiner, Andreas; Brassel, Kurt E. (1986). "Hierarchical grid structures for static geographic data bases" (PDF). In Blakemore, Michael (ed.). Auto Carto London: proceedings International Conference on the Acquisition, Management and Presentation of Spatial Data, 14–19 September 1986. Auto Carto London. pp. 485–496. ISBN 978-0-85406-312-3. OCLC 898826374.
  29. ^ Yan, Jin; Song, Xiao; Gong, Guanghong (2016). "Averaged ratio between complementary profiles for evaluating shape distortions of map projections and spherical hierarchical tessellations". Computers & Geosciences. 87: 41–55. Bibcode:2016CG.....87...41Y. doi:10.1016/j.cageo.2015.11.009.
  30. ^ "What3words: Find and share very precise locations via Google Maps with just 3 words". 2013-07-02. Retrieved 8 July 2014.
  31. ^ "Open Location Code is a library to generate short codes that can be used like street addresses, for places where street addresses don't exist.: Google/open-location-code". GitHub. 2019-02-18.
  32. ^ "Computational geometry and spatial indexing on the sphere: Google/s2geometry". GitHub. 2019-02-18.

External links