Topic 5: Terrain Analysis

1. Topic 5: Terrain Analysis 第五讲 地 形 分 析

2. Chapter 12 12.1 Introduction 12.2 Data for Terrain Mapping and Analysis 12.2.1 DEM 12.2.2 TIN 12.3 Terrain Mapping 12.3.1 Contouring 12.3.2 Vertical Profiling 12.3.3 Hill Shading 12.3.4 Hypsometric Tinting 12.3.5 Perspective View

3. 12.4 Terrain Analysis 12.4.1 Slope and Aspect 12.4.1.1 Computing Algorithms for Slope and Aspect Using Grid 12.4.1.2 Computing Algorithms for Slope and Aspect Using TIN 12.4.1.3 Factors Influencing Slope and Aspect Measures 12.4.2 Surface Curvature 12.4.3 Viewshed Analysis 12.4.4 Watershed Analysis 12.5 Grid versus TIN

4. Applications: Terrain Mapping and Analysis Task 1: Use DEM for Terrain Mapping Task 2: Perform Viewshed Analysis Task 3: Create a New Lookout Shapefile for Viewshed Analysis Task 4: Delineate Watersheds Task 5: Build and Display a TIN in ArcView

5. Data for Terrain Mapping and Analysis DEM TIN

6. DEMs The USGS (U.S. Geological Survey) classifies the quality of 7.5-minute DEMS into three levels, with Level 1 having the poorest quality. Level-1 accuracy has a RMSE (root mean square error) target of 7 meters and a maximum RMSE of 15 meters. Level-2 accuracy has a maximum RMSE of one half the contour interval. And Level-3 accuracy has a maximum RMSE of one-third of a contour interval—not to exceed 7 meters. Most USGS DEMs in use are either Level 1 or Level 2.

7. Global DEMs DEMs at different resolutions are now available on the global scale. ETOPO5 (Earth Topography-5 Minute) data cover both the land surface and ocean floor of the Earth, with a grid spacing of 5 minutes of latitude by 5 minutes of longitude (http://edcwww.cr.usgs.gov/glis/hyper/guide/etopo5). Both GTOPO30 (http://edcdaac.usgs.gov/gtopo30/gtopo30.html/) and GLOBE (http://www.ngdc.noaa.gov/seg/topo/globe.shtml/) offer global DEMs with a horizontal grid spacing of 30 arc seconds or approximately 1 kilometer.

8. LIDAR LIDAR data have the advantage over USGS DEMs in providing high-quality DEMs with a spatial resolution of 0.5 to 2 meters and a vertical accuracy of about 15 centimeters. a – Canopy elevation b – Ground elevation c – Estimated tree height

9. TIN A TIN approximates the land surface with a series of non-overlapping triangles. Elevation values (z values) along with x-, y-coordinates are stored at nodes that make up the triangles. In contrast to DEMs, TINs are based on an irregular distribution of elevation points.

10. Data Sources for TINs GIS users typically compile TINs using DEMs as the primary data source in a process sometimes referred to as conversion of DEM to TIN. Additional point data may include surveyed elevation points, GPS (Global Positioning System), and LIDAR data. Line data may include contour lines and breaklines. Breaklines are line features that represent changes of the land surface such as streams, shorelines, ridges, and roads. And area data may include lakes and reservoirs.

11. DEM to TIN Conversion • VIP (Very Important Points) used by ArcInfo Workstation • Maximum z-tolerance used by 3D Analyst and ArcInfo Workstation.

12. a A B C H P D G F E Ph s d b Pe G P C VIP evaluates the significance of an elevation point by measuring how well its value can be estimated from the neighboring point values. Figure 12.1b shows the case of using elevations at G and C to estimate the elevation at P. Pe is the estimated elevation, Ph is the actual elevation, and d represents the offset between Pe and Ph. Rather than using d as the measure of significance, VIP uses s.

13. Maximum Z-tolerance The maximum z-tolerance algorithm selects points from an elevation grid to construct a TIN such that, for every point in the elevation grid, the difference between the original elevation and the estimated elevation from the TIN is within the specified maximum z-tolerance. The algorithm uses an iterative process.

14. TIN to DEM Conversion The process requires each elevation point of the DEM to be estimated (interpolated) from its neighboring nodes that make up the TIN. Each of these nodes has its x-, y-coordinates as well as a Z (elevation) value. The default method used by ArcGIS to convert a TIN to a DEM is local first-order polynomial interpolation.

15. Terrain Mapping • Contouring • Vertical profiling • Hill shading • Hypsometric tinting • Perspective view

16. 950 900 1000 850 850 800 The contour line of 900 connects points that are interpolated to have the value of 900 along the triangle edges.

17. A vertical profile showing changes in elevation along a stream tributary. The profile has a vertical exaggeration factor of 1.0 (i.e., or no vertical exaggeration).

18. An example of hill shading, with the sun’s azimuth at 3150 (NW) and the sun’s altitude at 450

20. An example of a 3-D perspective view

21. An example of 3-D draping. In this case, streams and shorelines are draped on a 3-D surface

22. Terrain Analysis • Slope and Aspect • Surface curvature • Viewshed analysis • Watershed analysis

23. a θ b Percent slope in the diagram is 100 x (a/b). a is the vertical distance or rise, and b is the horizontal distance or run. Degree slope can be calculated by arctan (a/b).

24. N W E S N NE NW W E SE SW S Aspect is a directional measure in degrees. Aspect measures are often grouped into 4 principal directions (top) or 8 principal directions (bottom).

25. Z North Y Slope Aspect X South The normal vector to the cell is the directed line perpendicular to the cell. The quantity and direction of tilt of the normal vector determine the slope and aspect of the cell. (Redrawn from Hodgson, 1998, CaGIS vol. 25, no. 3, pp. 173–185; reprinted with the permission of the American Congress on Surveying and Mapping.)

26. e2 e1 C0 e3 e4 Slope Aspect Ritter’s algorithm for computing slope and aspect at C0 uses the four immediate neighbors of C0.

27. e1 e2 e3 e4 C0 e5 e6 e7 e8 Slope Aspect Horn’s algorithm for computing slope and aspect at C0 uses the eight neighboring cells of C0. The algorithm also applies a weight of 2 to e2, e4, e5, and e7, and a weight of 1 to e1, e6, e3, and e8.

28. A (X1,Y1,Z1) B (X2,Y2,Z2) C (X3,Y3,Z3) nx = (y2 - y1) (z3 - z1) - (y3 - y1) (z2 - z1) ny = (z2 - z1) (x3 - x1) - (z3 - z1) (x2 - x1) nz = (x2 - x1) (y3 - y1) - (x3 - x1) (y2 - y1) The algorithm for computing slope and aspect of a triangle in a TIN uses the x, y, and z values at the three nodes of the triangle.

29. Surface Curvature GIS applications in hydrological studies often require computation of surface curvature to determine if the surface at a cell location is upwardly convex or concave.

30. Viewshed analysis divides the study area into (1) not visible from observation points and (2) visible from observation points.

31. Viewshed analysis is based on the line of sight operation.

32. 1014 1011 1004 +11 +4 +1 1019 1015 1007 -4 +8 1021 1012 1025 -10 -6 +3 (b) (a) (c) The flow direction of the center cell in (a) is determined by first calculating the distance-weighted gradient to each of its eight neighbors. For the four immediate neighbors, the gradient is calculated by dividing the elevation difference between the center cell and the neighbor by 1. For the four corner neighbors, the gradient is calculated by dividing the elevation difference by 1.414. The results in (b) show that the steepest gradient, and therefore the flow direction, is from the center cell to the right cell (+8).

33. 996 1004 1011 1014 1007 1015 1019 999 (a) 1003 1021 1012 1025 1029 1020 1003 1033 (b) 0 0 1 2 0 0 2 6 (c) 0 0 2 3 0 3 1 2 The illustration shows a filled elevation grid (a), a flow direction grid (b), and a flow accumulation grid (c). Both shaded cells in (c) have the same flow accumulation value of 2. The top cell receives its flow from its left and lower left cells. The bottom cell receives its flow from its lower left cell, which already has a flow accumulation value of 1.

34. A flow accumulation grid is shown in (a). The darkness of the symbol corresponds to the flow accumulation value. A drainage network derived from the flow accumulation grid using a threshold value of 500 cells is shown in (b). The delineation of watersheds is shown in (c).

35. A watershed (shaded area) is derived for each of the three pour points.

36. Grid or TIN? Advantages of using TIN A main advantage of using a TIN lies in the flexibility with input data sources. One can construct a TIN using inputs from DEM, contour lines, GPS data, LIDAR data, and survey data. The user can add elevation points to a TIN at their precise locations and add breaklines, such as streams, roads, ridgelines, and shorelines, to define surface discontinuities. Besides data flexibility, TIN is also an excellent data model for terrain mapping and 3-D display. The triangular facets of a TIN better define the land surface than an elevation grid and create a sharper image. Most GIS users seem to prefer the look of a map based on a TIN rather than an elevation grid.

37. Grid or TIN? Advantages of using Grid Computational efficiency is the main advantage of using grids for terrain analysis. The simple data structure makes it relatively easy to perform neighborhood operations on an elevation grid.

38. Thelin and Pike’s (1991) digital shaded-relief map of the United States http://www.usgs.gov/reports/misc/Misc._Investigations_Series_Maps_(I_Series)/I_2206/usa_dem.gif USGS’ National Hydrography Dataset http://nhd.usgs.gov/

39. Thank You ! Thank You !