DITM: “tf1732_c015” — 2004/10/25 — 12:48 — page 297 — #1 CHAPTER 15 Beyond Digital Terrain Modeling Chapter 14 discusses the more traditional development and applications of terrain models. Here, we look at some extensions of these models for specific problems. 15.1 DIGITAL TERRAIN MODELING WITH COMPLEX CONSTRUCTION 15.1.1 Manual Addition of Constructions on Terrain Surface For simplicity it is usually assumed that terrain models are monotonic in X and Y — there is only one possible Z for each XY location. This is often true in the real world, but not always — occasionally there are caves, tunnels, overhanging cliffs, bridges, and overpasses. In the work by Tse and Gold (2002), the standard TIN model is extended by merging some aspects of terrain modeling (TINs), com- putational geometry (the Quad-Edge data structure), and computer aided design or CAD (Euler operators, which guarantee to preserve the connectivity of the surface after they are applied). They found it easy to combine them to give the usual oper- ations on a 2D triangulation — as well as add an operator that generates a hole between any two nonadjacent triangles (which is really the same thing as adding a bridge or handle to the surface). Figure 15.1, Figure 15.2, and Figure 15.3 give simple examples, and Figure 15.4 and Figure 15.5 show part of a Hong Kong city model. Thus, a simple modification of the basic triangulation algorithm allows one to interactively modify the terrain model to add complex features that are other- wise unavailable. Because one is still forming a connected surface, a variety of topological operations, such as neighborhood selection and flow modeling, may be performed. Clearly another, even higher, layer of operations would permit one to add predesigned features such as buildings, dams, tunnels, etc. to our terrain model. 297 © 2005 by CRC Press DITM: “tf1732_c015” — 2004/10/25 — 12:48 — page 298 — #2 298 DIGITAL TERRAIN MODELING: PRINCIPLES AND METHODOLOGY (a) (b) Figure 15.1 TIN model of a surface with a tunnel: (a) a tunnel on a TIN model and (b) the enlarged tunnel. (a) (b) Figure 15.2 TIN model of a surface with bridges: (a) two bridges on a TIN model and (b) the enlarged bridge. (a) (b) Figure 15.3 TIN model of a ground surface with buildings and bridge: (a) a building on a TIN model and (b) a bridge connecting two buildings. 15.1.2 Semiautomated Modification of the Terrain Surface Section 15.1.1 showed simple terrain modification based on modifying the TIN by adding and deleting individual points with (X, Y ,Z) coordinates. This is effective but slow to do by hand. An alternative approach is to “cut” the triangulated surface with a “knife” in order to sculpt it to the form desired. One first sets the knife size, location, and orientation and then performs the cut (or intersect) operation. One may © 2005 by CRC Press DITM: “tf1732_c015” — 2004/10/25 — 12:48 — page 299 — #3 BEYOND DIGITAL TERRAIN MODELING 299 Figure 15.4 A partial view of Hong Kong harbor. Figure 15.5 Part of the Hong Kong city model. either lower the terrain surface to the knife position, for example, cutting into the side of a hill, or else raise the surface to the knife position, creating an embankment or dam. More points are added to the triangulation to form the intersection lines between the knife and the original terrain, and one assumes a maximum slope (less than vertical) for the edges of the cut or embankment. Figure 15.6(a) shows a simple TIN model with the knife in place. Figure 15.6(b) shows the result after the surface is lowered to the knife (with a 45 ◦ embankment specified). Figure 15.6(c) shows the knife positioned across a valley, and Figure 15.6(d) shows the result of raising the terrain surface to the knife, forming a dam structure across the valley. © 2005 by CRC Press DITM: “tf1732_c015” — 2004/10/25 — 12:48 — page 300 — #4 300 DIGITAL TERRAIN MODELING: PRINCIPLES AND METHODOLOGY (a) (c) (b) (d) Figure 15.6 Terrain modification by “cutting” the triangulated surface with a “knife.” (a) The knife positioned on the terrain. (b) Modified terrain after lowering the surface to the knife blade. (c) The knife positioned across a valley. (d) The dam created after raising the surface to the knife blade. 15.2 DIGITAL TERRAIN MODELING ON THE SPHERE With the introduction of the concept digital earth, global modeling of the Earth’s surface (Gold and Mostafavi 2000) has become a hot topic. The digital terrain modeling techniques described in this book can also be extended to spherical terrain modeling. 15.2.1 Generation of TIN and Voronoi Diagram on Sphere In planimetric terrain modeling, as discussed in Chapter 4, grids and TINs have been widely used to tessellate the terrain. On the sphere, a similar tessellation model needs to be used. The concept of spherical surface tessellation was presented by Fuller, a German cartographer, for map projection in the 1940s (Dutton 1996). Since then, many researchers have approached this problem to project, analyze, and index global data. Many methods are based on inscribed polyhedrons, such as the tetra- hedron, the cube (Snyder 1992), the octahedron (Dutton 1989, 1996; Goodchild et al. 1991; Goodchild and Yang 1992; Otoo and Zhu 1993; Clarke and Mulcahy 1995), the dodecahedron (Wickman and Elvers 1974), and the icosahedron (Fekete 1990; White et al. 1992; Lee and Samet 2000), as shown in Figure 15.7. The edges of the polyhedron are projected to the spherical surface and form the edges of spherical triangles. The octahedron-based tessellation is a regular triangular mesh on the sphere, called the octahedral quaternary triangular mesh (O-QTM). Figure 15.8 shows an © 2005 by CRC Press DITM: “tf1732_c015” — 2004/10/25 — 12:48 — page 301 — #5 BEYOND DIGITAL TERRAIN MODELING 301 (a) (b) Figure 15.7 Spherical surface tessellation based on inscribed polyhedra (Reprinted with per- mission fromWhite et al. 1992): (a) five polyhedra and (b) projectedto the spherical surface. (a) (b) (c) Figure 15.8 Hierarchical tessellation of the spherical facet based on octahedron (Dutton 1996): (a) level 1; (b) level 2; and (c) level 3. example of O-QTM at three difference levels (Dutton 1996). Terrain modeling can then be applied to the QTM. The QTM can also be used as a coordinate system on the sphere, just like the regular grid or triangular network on a 2D plane. In the QTM, a point is represented by a triangle, an arc is represented by a series of neighbor triangles, and a region is represented by a series of neighbor triangles on and within its boundary trace. From the QTM, a TIN can then be constructed. Alternatively, spherical TINs can also be derived from spherical Voronoi diagrams (Augenbaum 1985; Robert 1997; Chen et al. 2003). Figure 15.9 shows an example of spherical Voronoi diagram and its dual — the spherical TIN. 15.2.2 Voronoi Diagram for Modeling Changes in Sea Level on Sphere Mostafavi and Gold (2004) used the dynamic Voronoi diagram on the sphere to model the continually changing height of the sea, rather than of terrain. Figure 15.10(a) shows an initial set of cells, each representing a fixed mass of water, and uses the free Lagrange method to simulate flow under lunar gravitational influence, and hence the sea height. Coastlines were modeled by a double line of fixed Voronoi cell generators. Figure 15.10(b)shows theresult after simulationstarted: high water (HW) isindicated by smaller, and therefore higher, cells, while low water (LW) is shown by larger, © 2005 by CRC Press DITM: “tf1732_c015” — 2004/10/25 — 12:48 — page 302 — #6 302 DIGITAL TERRAIN MODELING: PRINCIPLES AND METHODOLOGY (a) (b) Figure 15.9 Spherical TIN formation: (a) Voronoi diagram and (b) Delaunay triangulation. HW (a) (b) LW Figure 15.10 Dynamic Voronoi diagram on sphere to model the continually changing heights of sea: (a) initial configuration of Voronoi cells and (b) Voronoi cell configuration indicating lunar tides. lower cells. Figure 15.11 is a Mercator projection showing the flow directions and velocities of each cell. 15.3 THREE-DIMENSIONAL VOLUMETRIC MODELING Two dimensions are required in terrain modeling to generate the underlying triangulation or grid. Once elevations are added as an attribute, the result is usually known as “2.5D” modeling, although the data structures remain 2D. Once the topol- ogy (or connectedness) can no longer be represented on the plane (as in 3D objects in CAD or games), a surface representation, composed usually of triangles, is often used, as in Section 15.1. However, for some applications a surface model is inappropriate, and a full 3D volumetric model is needed. Examples include geological, atmospheric, and oceano- graphic models, where attributes need to be assigned to arbitrary locations in 3D space. In some cases a 3D grid may be used, or an octree where nodes repre- sent volumes. A more flexible approach is to replace the 2D triangulation structure © 2005 by CRC Press DITM: “tf1732_c015” — 2004/10/25 — 12:48 — page 303 — #7 BEYOND DIGITAL TERRAIN MODELING 303 Figure 15.11 Mercator projection showing cell velocities and directions. Figure 15.12 Delaunay and Voronoi cells in three dimensions. © 2005 by CRC Press DITM: “tf1732_c015” — 2004/10/25 — 12:48 — page 304 — #8 304 DIGITAL TERRAIN MODELING: PRINCIPLES AND METHODOLOGY Figure 15.13 A three-dimensional isosurface. The color plate can be viewed at http://www. crcpress.com/e_products/downloads/download.asp?cat_no=TF1732. with a 3D Delaunay tetrahedral model, thus allowing the connection of arbitrarily located observations, andthen toadd the3D Voronoicells. Whileconceptually simple, the implementation is often difficult, due to the large number of degenerate (coplanar, cocircular) cases. Figure 15.12 shows a 3D data set with one Delaunay cell and its neighbors shown on the left and a Voronoi cell on the right. Various interpolation techniques, such as the equivalent of the 2D Sibson or natural neighbor interpolation, may be used, and these behave well for the aniso- tropic distributions of data that are often found in three dimensions. For visualization purposes the individual tetrahedral may be sliced, based on the values at the corner vertices, to give the 3D equivalent of 2D contours. Figure 15.13 shows a single 3D isosurface constructed in this way. © 2005 by CRC Press DITM: “tf1732_c016” — 2004/10/20 — 15:45 — page 305 — #1 Epilogue It was natural that we felt relieved and excited somehow after having completed the final draft of this book and having uploaded the materials onto the ftp site of the publisher. However, soon we started to feel obliged to write this epilogue because there are a few issues confronting us. We thought it is really a pity that no authored book in this discipline had been made after over 40 years of development although there are two edited works, Terrain Modelling inSurveying andCivil Engineeringby Petrie andKennie(1990) andDigital Elevation Model Techniques and Applications: The DEM User Manual by Maune (2001). Our aim was to write a book systematically covering a wide range of topics in digital terrain modeling so as to fill in the gap in this area. While writing, we were faced with a number of challenges. The first challenge was related to the selective omission of materials. It was difficult to make decisions. This is because the term “digital terrain modeling” would mean different things to different groups of terrainspecialists and practitioners. To the producers (including photogrammetrists and surveyors), data acquisition and terrain surface modeling are of most concern; to geographers, terrain analysis and appli- cations are the most important; to geologists, interpolation techniques seem to be critical; It is really hard to satisfy all these groups. In the end, we decided that those topics are simplified if they have rich bodiesof literature available. For example, we did notinclude many algorithmsand techniques for interpolationand triangulation as there is a huge body of literature (e.g., Su 1989; Chin 1995; Sakhnovich 1997; de Berg 2000; Phillips 2003) in these areas covering the techniques developed in computational geometry and geosciences. Contouring is a traditional topic in digital terrain modeling but is only briefly discussed in this book because a book authored by Watson (1992) has been dedicated to this topic. Similarly, DTM-based terrain analysis is briefly discussed because of a recent book edited by Wilson and Gallant (2000). The second challenge was related to the depth of discussion. We may disappoint those readers who are interested in mathematics because we present neither mathe- matical proofs nor technical details. Indeed, it is the main aim of this book to present a systematic accounting of stories in digital terrain modeling at the level of principles and methodology, as the title of the book suggests. 305 © 2005 by CRC Press DITM: “tf1732_c016” — 2004/10/20 — 15:45 — page 306 — #2 306 DIGITAL TERRAIN MODELING: PRINCIPLES AND METHODOLOGY The third challenge was related to the boundary of the discipline. Because the terrain surface is of concern to all geosciences and a huge body of related model- ing methodology and applications is available, we have to cut down the contents somewhere. Therefore, we did not cover much on Voronoi diagrams (e.g., Davies 2000; Okabe et al. 2000) although we are very interested in this topic. Similarly, we did not cover much on geostatistics (e.g., Olea 1999) and even omitted the famous Kriging technique. We only simply mentioned the surface modeling on sphere and with construction in Chapter 15. All in all, we are pleased with the compilation of some materials presented to you, but also feel guilty about the imperfection. Your comments are appreciated so that we could make improvements in another edition, if possible. © 2005 by CRC Press [...]... 312 DIGITAL TERRAIN MODELING: PRINCIPLES AND METHODOLOGY Kraak, M.J and Brown, A 2001 Web Cartography: Developments and Prospectus Taylor & Francis, London Kraus, K and Pfeifer, N 1998 Determination of terrain models in wooded areas with airborne laser scanner data ISPRS Journal of Photogrammetry and Remote Sensing, 53(4):193–203 Kubik, K and Botman, A 1976 Interpolation accuracy for topographic and. .. T (Eds.), Terrain Modelling in Surveying and Civil Engineering, Whittles Publishing, Caithness, England, pp 112–137 Petrie, G 1990b Terrain data acquisition and modelling from existing maps In: Petrie, G and Kennie, T (Eds.), Terrain Modelling in Surveying and Civil Engineering Whittles Publishing, Caithness, England, pp 85–111 Petrie, G and Kennie T (Eds.) 1987 An introduction to terrain modeling:. .. 155 :636–638 Mandelbrot, B 1981 The Fractal Geometry of Nature W.H Freeman and Company, San Francisco, CA Mark, D 1975 Geomorphometric parameters: a review and evaluation Geografiska Annaler, 57A:165–177 Masry, S.E 1974 Digital correlation principles Photogrammetric Engineering, 40(3): 303–308 © 2005 by CRC Press DITM: “tf1732_c017” — 2004/10/22 — 16:38 — page 313 — #7 314 DIGITAL TERRAIN MODELING: PRINCIPLES. .. of Speech and Visual Form, MIT Press, Cambridge, MA, pp 153 –171 Borgefors, G 1986 Distance transformations in digital images Computer Vision, Graphics and Image Processing, 34:344–371 Bowyer, A 1981 Computing Dirichlet tessellations The Computer Journal, 24(2):162–166 307 © 2005 by CRC Press DITM: “tf1732_c017” — 2004/10/22 — 16:38 — page 307 — #1 308 DIGITAL TERRAIN MODELING: PRINCIPLES AND METHODOLOGY. .. Engineering and Remote Sensing, 48(3):375–381 Serra, J 1982 Image Processing and Mathematical Morphology Academic Press, New York © 2005 by CRC Press DITM: “tf1732_c017” — 2004/10/22 — 16:38 — page 315 — #9 316 DIGITAL TERRAIN MODELING: PRINCIPLES AND METHODOLOGY Sharpnack, D.A and Akin, G 1969 An algorithm for computing slope and aspect from elevations Photogrammetric Survey, 35:247–248 Sibson, R 1980... REFERENCES 315 Ostman, A 1986 A graphic editor for digital elevation models Geo-processing, 3: 143 154 Otoo, E and Zhu, H 1993 Indexing on spherical surfaces using semi-quadcodes In: Abel, J and Beng, C.O (Eds.), Advances in Spatial Databases, Lecture Notes in Computer Science 692, Springer, Singapore, pp 509–529 Petrie, G 1990a Modelling, interpolation and contouring procedures In: Petrie, G and Kennie,... Geometry: Algorithms and Applications, 2nd ed Springer, Berlin de Berg, M and Dobrindt, K 1998 On levels of detail in terrains Graphical Models and Image Processing, 60(1):1–12 de Floriani, L 1989 A pyramidal data structure for triangle-based surface description IEEE Computer Graphics and Applications, 9(2):67–78 De Floriani, L and Magillo, P 1994 Visibility algorithms on triangulated digital terrain models... Phong, B.-T 1975 Illumination for computer-generated pictures Communication of ACM, 18(6):311–317 Polidori, L 1991 Digital terrain models from radar images: a review In: Guyenne, T.D and Hunt, J.J (Eds.), Proceedings of the International Symposium on Radars and Lidars in Earth and Planetary Science, France, Cannes, pp 141–146 Quattrochi, D.A and Goodchild, M.F (Eds.) 1997 Scale in Remote Sensing and GIS... Photogrammetry and Remote Sensing, 49:29–33 Flood, M 2001 Laser altimetry: from science to commercial Lidar mapping Photogrammetric Engineering & Remote Sensing, 67(11):1209–1218 © 2005 by CRC Press DITM: “tf1732_c017” — 2004/10/22 — 16:38 — page 309 — #3 310 DIGITAL TERRAIN MODELING: PRINCIPLES AND METHODOLOGY Forture, S 1975 A sweep line method for Voronoi diagram Methodica, 2 :153 –174 Frederiksen, P 1980 Terrain. .. Accuracy of Digital Elevation Models Science Press, Beijing Tang, L 1989 Surface modelling and visualization based upon digital image processing techniques In: Grun, A and Kahmen, H.K (Eds.), Optical 3-D Measurement Techniques, Wichmann Verlag, Karlsruhe, Germany, pp 317–325 Tempfli, K 1980 Spectral analysis of terrain relief for the accuracy estimation of digital terrain models ITC Journal, 198 0-3 :487–510 . etc. to our terrain model. 297 © 2005 by CRC Press DITM: “tf1732_c 015 — 2004/10/25 — 12:48 — page 298 — #2 298 DIGITAL TERRAIN MODELING: PRINCIPLES AND METHODOLOGY (a) (b) Figure 15. 1 TIN model. Press DITM: “tf1732_c 015 — 2004/10/25 — 12:48 — page 302 — #6 302 DIGITAL TERRAIN MODELING: PRINCIPLES AND METHODOLOGY (a) (b) Figure 15. 9 Spherical TIN formation: (a) Voronoi diagram and (b) Delaunay. — #4 310 DIGITAL TERRAIN MODELING: PRINCIPLES AND METHODOLOGY Forture, S. 1975. A sweep line method for Voronoi diagram. Methodica, 2 :153 –174. Frederiksen, P. 1980. Terrain analysis and accuracy