Beyond Mapping IV

Topic 10 – Future Directions and Trends (Further Reading)

 

 

GIS Modeling book

 

 

From a Map Pancake to a Soufflé — continues the discussion of concepts and configuration of a 3D GIS (December 2009)

 

<Click here> for a printer-friendly version of this topic (.pdf).

(Back to the Table of Contents)
______________________________

From a Map Pancake to a Soufflé

(GeoWorld, December 2009)  

(return to top of Topic)

 

As the Time Traveller noted in H. G. Wells’ classic “The Time Machine,” the real world has three geometric dimensions not simply the two we commonly use in mapping.  In fact, he further suggested that “...any real body must have extension in four directions: Length, Breadth, Thickness—and Duration (time)” …but that’s a whole other story.  

 

Recall from earlier discussion of 3D GIS that Geodetic Referencing (geographic position) used in identifying an “areal extent” in two-dimensions on the earth’s surface can be extended to a Database Referencing system (matrix location) effectively defining a 3-dimensional “project block” (see the left side of figure 1).  The key is the use of Geodetic Height above and below the earth’s ellipsoid as measured along the perpendicular from the ellipsoid to provide the vertical (Z) axis for any location in 3-dimensional space. 

 

The result is a coordinate system of columns (X), rows (Y), and verticals (Z) defining an imaginary matrix of grid elements, or “voxels,” that are a direct conceptual extension of the “pixels” in a 2D raster image.  For example, the top-right inset in figure 1 shows a 3D map of a cave system using ArcGIS 3D Analyst software.  The X, Y and Z positioning forms a 3-dimensional display of the network of interconnecting subterranean passages.  The lower-right inset shows an analogous network of blood vessels for the human body except at much different scale.  The important point is that both renderings are 3D visualizations and not a 3D GIS as they are unable to perform volumetric analyzes, such as directional flows along the passageways.

 

The distinction between 3D visualization and analytical systems arises from differences in their data structures.  A 3D visualization system stores just three values—X and Y for “where” and Z for “what (elevation).”   A 3-dimensional mapping system stores at least four values—X, Y and Z for “where” and an attribute value for “what” describing the characteristic/condition at each location within a project block. 

 

T27_9a

 

Figure 1.  Storage of a vertical (Z) coordinate extends traditional 2D mapping to 3D volumetric representation.

 

Figure 2 illustrates two ways of storing 3-dimensional grid data.  A Flat File stores a single map value for each grid element in a map block.  The individual records can explicitly identify each grid element (grayed columns—“where”) along with the attribute (black column—“what”).  Or, much more efficiently, the information can be implicitly organized as a header line containing the grid block configuration/size/referencing followed by a long string of numbers with each value’s position in the string determining its location in the block through standard nested programming loops.  This shortened format provides for advanced compression techniques similar to those used in image files to greatly reduce file size.  

 

An alternative strategy, termed an Interleaved File, stores a series of map attributes as separate fields for each record that in turn represents each grid element, either implicitly or explicitly organized into a table.  Note that in the interleaved file in figure 2, the map values for Elevation, %Slope and Cover type identify surface characteristics with a “null value (---)” assigned to grid elements both above and below the surface.  Soil type, on the other hand, contains values for the grid elements on and immediately below the surface with null values only assigned to locations above ground.  This format reduces the number of files in a data set but complicates compression and has high table maintenance overhead for adding and deleting maps.

 

T27_9b

 

Figure 2.  A 3-dimensional matrix structure can be used to organize volumetric mapped data.

 

Figure 3 outlines some broader issues and future directions in 3D GIS data storage and processing.  The top portion suggests that the inconsistent geometry of the traditional Cube results in differing distances and facet adjacency relationship to the surrounding twenty six neighbors, thereby making a cube a poor grid element for 3D data storage.  A Dodecahedron, on the other hand, aligns with a consistent set of 12 equidistant pentagonal faces that “nest” without gaps ...an important condition in spatial analysis of movements, flows and proximal conditions.

 

The lower portion of figure 3 illustrates the knurly reality of geographic referencing in 3-dimensions—things change as distance from the center of the earth or bounding ellipsoid changes.   Nicely nesting grid elements of a fixed size separate as distance increases (diverge); overlap as distance decreases (converge).  To maintain a “close-packing” arrangement either the size of the grid element needs to adjust or the progressive errors of fixed size zones are tolerated.

 

Similar historical changes in mapping paradigms and procedures occurred when we moved from a flat earth perspective to a round earth one that generated a lot of room for rethinking.  There are likely some soon-to-be-famous mathematicians and geographers who will match the likes of Claudius Ptolemy (90-170), Gerardus Mercator (1512-1594) and Rene Descartes (1596-1650)— I wonder who among us will take us beyond mapping as we know it?     

 

T27_9c

 

Figure 3.  Alternative grid element shapes and new procedures for dealing with radial divergence form the basis for continued 3D GIS research and development.

 

_____________________________

 

Author’s Notes:  A good discussion of polyhedral and other 3-dimensional coordinate systems is in Topic 12, “Modeling locational uncertainty via hierarchical tessellation,” by Geoffrey Dutton in Accuracy of Spatial Databases edited by Goodchild and Gopal.  In his discussion he notes that “One common objection to polyhedral data models for GIS is that computations on the sphere are quite cumbersome … and for many applications the spherical/geographical coordinates … must be converted to and from Cartesian coordinates for input and output.”  

 

 (return to top of Topic)

 

(Back to the Table of Contents)