{"id":499,"date":"2014-01-24T21:12:05","date_gmt":"2014-01-24T21:12:05","guid":{"rendered":"http:\/\/opentextbc.ca\/natureofgeographicinformation\/?post_type=chapter&#038;p=499"},"modified":"2015-05-27T15:45:06","modified_gmt":"2015-05-27T15:45:06","slug":"1-overview-6","status":"publish","type":"chapter","link":"https:\/\/opentextbc.ca\/natureofgeographicinformation\/chapter\/1-overview-6\/","title":{"raw":"National Spatial Data Infrastructure II","rendered":"National Spatial Data Infrastructure II"},"content":{"raw":"<h2>7.1. Overview<\/h2>\r\nChapters 6 and 7 consider the origins and characteristics of the <strong>framework data themes <\/strong>that make up the United States' proposed National Spatial Data Infrastructure (NSDI). Chapter 6 discussed the geodetic control and orthoimagery themes. <strong>This chapter describes the origins, characteristics and current status of the elevation, transportation, hydrography, governmental units and cadastral themes.<\/strong>\r\n<h3>Objectives<\/h3>\r\nStudents who successfully complete Chapter 7 should be able to:\r\n<ol>\r\n\t<li>Given a regular or irregular array of spot elevations, construct a triangulated irregular network, interpolate contour intervals and draw contour lines;<\/li>\r\n\t<li>Compare vector and raster representations of terrain elevation;<\/li>\r\n\t<li>Acquire and view digital elevation data from the National Elevation Dataset;<\/li>\r\n\t<li>Calculate an interpolated spot elevation based on neighboring elevations;<\/li>\r\n\t<li>Contrast the characteristics of three global elevation data products;<\/li>\r\n\t<li>Describe the characteristics and current status of the NSDI hydrography, transportation, and governmental units themes as implemented in USGS' National Map; and<\/li>\r\n\t<li>Interpret the size and relative location of a land parcel designated in terms of the U.S. Public Land Survey System.<\/li>\r\n<\/ol>\r\n<h3>Comments and Questions<\/h3>\r\nRegistered students are welcome to post comments, questions, and replies to questions about the text. Particularly welcome are anecdotes that relate the chapter text to your personal or professional experience. In addition, there are discussion forums available in the ANGEL course management system for comments and questions about topics that you may not wish to share with the whole world.\r\n\r\nTo post a comment, scroll down to the text box under \"Post new comment\" and begin typing in the text box, or you can choose to reply to an existing thread. When you are finished typing, click on either the \"Preview\" or \"Save\" button (Save will actually submit your comment). Once your comment is posted, you will be able to edit or delete it as needed. In addition, you will be able to reply to other posts at any time.\r\n\r\nNote: the first few words of each comment become its \"title\" in the thread.\r\n<h2>7.2. Checklist<\/h2>\r\nThe following checklist is for Penn State students who are registered for classes in which this text, and associated quizzes and projects in the ANGEL course management system, have been assigned. You may find it useful to print this page out first so that you can follow along with the directions.\r\n<table summary=\"Tasks to be completed for the chapter\"><caption>Chapter 7 Checklist (for registered students only)<\/caption>\r\n<thead>\r\n<tr>\r\n<th>Step<\/th>\r\n<th>Activity<\/th>\r\n<th>Access\/Directions<\/th>\r\n<\/tr>\r\n<\/thead>\r\n<tbody>\r\n<tr>\r\n<th>1<\/th>\r\n<td><strong>Read<\/strong>\u00a0Chapter 7<\/td>\r\n<td>This is the second page of the Chapter. Click on the links at the bottom of the page to continue or to return to the previous page, or to go to the top of the chapter. You can also navigate the text via the links in the GEOG 482 menu on the left.<\/td>\r\n<\/tr>\r\n<tr>\r\n<th>2<\/th>\r\n<td>Submit\u00a0<strong>3 practice quizzes<\/strong>including:\r\n<ul>\r\n\t<li>Contouring<\/li>\r\n\t<li>DLGs and DEMs<\/li>\r\n\t<li>Interpolation<\/li>\r\n<\/ul>\r\nPractice quizzes are not graded and may be submitted more than once.<\/td>\r\n<td>Go to ANGEL &gt; [your course section] &gt; Lessons tab &gt; Chapter 7 folder &gt; [quiz]<\/td>\r\n<\/tr>\r\n<tr>\r\n<th>3<\/th>\r\n<td>Perform\u00a0<strong>\u201cTry this\u201d activities<\/strong>including:\r\n<ul>\r\n\t<li>Draw a contour map<\/li>\r\n\t<li>Explore Digital Line Graph hypsography<\/li>\r\n\t<li>Explore a Digital Elevation Model<\/li>\r\n\t<li>Download and view an extract from the National Elevation Dataset<\/li>\r\n<\/ul>\r\n\u201cTry this\u201d activities are not graded.<\/td>\r\n<td>Instructions are provided for each activity.<\/td>\r\n<\/tr>\r\n<tr>\r\n<th>4<\/th>\r\n<td>Submit the<strong>Chapter 7 Graded Quiz<\/strong><\/td>\r\n<td>ANGEL &gt; [your course section] &gt; Lessons tab &gt; Chapter 7 folder &gt; Chapter 7 Graded Quiz. See the Calendar tab in ANGEL for due dates.<\/td>\r\n<\/tr>\r\n<tr>\r\n<th>5<\/th>\r\n<td>\u00a0Read<strong>comments and questions<\/strong>posted by fellow students. Add comments and questions of your own, if any.<\/td>\r\n<td>\u00a0Comments and questions may be posted on any page of the text, or in a Chapter-specific discussion forum in ANGEL.<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n&nbsp;\r\n<h2>7.3. Theme: Elevation<\/h2>\r\nThe NSDI\u00a0<em>Framework Introduction and Guide<\/em>\u00a0(FGDC, 1997, p. 19) points out that \u201celevation data are used in many different applications.\u201d Civilian applications include flood plain delineation, road planning and construction, drainage, runoff, and soil loss calculations, and cell tower placement, among many others. Elevation data are also used to depict the terrain surface by a variety of means, from contours to relief shading and three-dimensional perspective views.\r\n\r\nThe NSDI Framework calls for an \u201celevation matrix\u201d for land surfaces. That is, the terrain is to be represented as a grid of elevation values. The spacing (or resolution) of the elevation grid may vary between areas of high and low relief (i.e., hilly and flat). Specifically, the Framework Introduction states that\r\n\r\nElevation values will be collected at a post-spacing of 2 arc-seconds (approximately 47.4 meters at 40\u00b0 latitude) or finer. In areas of low relief, a spacing of 1\/2 arc-second (approximately 11.8 meters at 40\u00b0 latitude) or finer will be sought (FGDC, 1997, p. 18).\r\n\r\nThe elevation theme also includes bathymetry\u2013depths below water surfaces\u2013for coastal zones and inland water bodies. Specifically,\r\n\r\nFor depths, the framework consists of soundings and a gridded bottom model. Water depth is determined relative to a specific vertical reference surface, usually derived from tidal observations. In the future, this vertical reference may be based on a global model of the geoid or the ellipsoid, which is the reference for expressing height measurements in the Global Positioning System (Ibid).\r\n\r\nUSGS has lead responsibility for the elevation theme. Elevation is also a key component of USGS\u2019 National Map. The next several pages consider how heights and depths are created, how they are represented in digital geographic data, and how they may be depicted cartographically.\r\n<h2>7.4. Vector and Raster Approaches<\/h2>\r\n&nbsp;\r\n\r\nThe terms\u00a0<strong>raster<\/strong>\u00a0and\u00a0<strong>vector<\/strong>\u00a0were introduced back in Chapter 1 to denote two fundamentally different strategies for representing geographic phenomena. Both strategies involve simplifying the infinite complexity of the Earth\u2019s surface. As it relates to elevation data, the raster approach involves\u00a0<em>measuring elevation at a sample of locations<\/em>. The vector approach, on the other hand, involves\u00a0<em>measuring the locations of a sample of elevations<\/em>. I hope that this distinction will be clear to you by the end of this chapter.\r\n\r\n<img alt=\"Diagram comparing contours and elevation grid depicting the same surface\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/vector_raster_surfs.gif\" \/>\r\n\r\nVector and raster representations of the same terrain surface.\r\n\r\nThe illustration above compares how elevation data are represented in vector and raster data. On the left are elevation\u00a0<strong>contours<\/strong>, a vector representation that is familiar with anyone who has used a USGS topographic map. The technical term for an elevation contour is<em>isarithm<\/em>, from the Greek words for \u201csame\u201d and \u201cnumber.\u201d The terms<em>isoline<\/em>,\u00a0<em>isogram<\/em>, and\u00a0<em>isopleth<\/em>\u00a0all mean more or less the same thing. (See any cartography text for the distinctions.)\r\n\r\nAs you will see later in this chapter, when you explore Digital Line Graph hypsography data using Global Mapper or dlgv 32 Pro, elevations in vector data are encoded as attributes of line features. The distribution of elevation points across the quadrangle is therefore irregular. Raster elevation data, by contrast, consist of grids of points at which elevation is encoded at regular intervals. Raster elevation data are what\u2019s called for by the NSDI Framework and the USGS National Map. Digital contours can now be rendered easily from raster data. However, much of the raster elevation data used in the National Map was produced from digital vector contours and hydrography (streams and shorelines). For this reason we\u2019ll consider the vector approach to terrain representation first.\r\n<h2>7.5. Contours<\/h2>\r\n<img alt=\"Perspective view of a terrain surface showing contour lines as traces of parallel horizontal planes\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/raisz_model.gif\" \/>\r\n\r\n&nbsp;\r\n\r\nContour lines trace the elevation of the terrain surface at regularly-spaced intervals (Raisz, 1948. \u00a9 McGraw-Hill, Inc. Used by permission).\r\n\r\nDrawing contour lines is a way to represent a terrain surface with a sample of elevations. Instead of measuring and depicting elevation at every point, you measure only along lines at which a series of imaginary horizontal planes slice through the terrain surface. The more imaginary planes, the more contours, and the more detail is captured.\r\n\r\n<img alt=\"Plan view of contour lines used to depict a terrain surface\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/raisz_plan.gif\" \/>\r\n\r\nContour lines representing the same terrain as in the first figure, but in plan view. (Raisz, 1948. \u00a9 McGraw-Hill, Inc. Used by permission).\r\n\r\nUntil photogrammetric methods came of age in the 1950s, topographers in the field sketched contours on the USGS 15-minute topographic quadrangle series. Since then, contours shown on most of the 7.5-minute quads were compiled from stereoscopic images of the terrain, as described in Chapter 6. Today computer programs draw contours automatically from the spot elevations that photogrammetrists compile stereoscopically.\r\n\r\nAlthough it is uncommon to draw terrain elevation contours by hand these days, it is still worthwhile to know how. In the next few pages you\u2019ll have a chance to practice the technique, which is analogous to the way computers do it.\r\n<h2>7.6. Contouring By Hand<\/h2>\r\nThis page will walk you through a methodical approach to rendering contour lines from an array of spot elevations (Rabenhorst and McDermott, 1989). To get the most from this demonstration, I suggest that you print the\u00a0<a href=\"https:\/\/www.e-education.psu.edu\/files\/geog482\/image\/contouring_lesson.gif\">illustration in the attached image file<\/a>. Find a pencil (preferably one with an eraser!) and straightedge, and duplicate the steps illustrated below. A \u201cTry This!\u201d activity will follow this step-by-step introduction, providing you a chance to go solo.\r\n\r\n&nbsp;\r\n\r\n<img alt=\"Step 1 of contouring demonstration\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/contouring_2.gif\" \/>\r\n\r\nBeginning a triangulated irregular network.\r\n\r\nStarting at the highest elevation, draw straight lines to the nearest neighboring spot elevations. Once you have connected to all of the points that neighbor the highest point, begin again at the second highest elevation. (You will have to make some subjective decisions as to which points are \u201cneighbors\u201d and which are not.) Taking care not to draw triangles across the stream, continue until the surface is completely triangulated.\r\n\r\n<img alt=\"Step 2 of contouring demonstration\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/contouring_3.gif\" \/>\r\n\r\nComplete TIN. Note that the triangle sides must not cross hydrologic features (i.e., the stream) on a terrain surface.\r\n\r\nThe result is a\u00a0<strong>triangulated irregular network (TIN)<\/strong>. A TIN is a vector representation of a continuous surface that consists entirely of triangular facets. The vertices of the triangles are spot elevations that may have been measured in the field by leveling, or in a photogrammetrist\u2019s workshop with a stereoplotter, or by other means. (Spot elevations produced photogrammetrically are called\u00a0<strong>mass points<\/strong>.) A useful characteristic of TINs is that each triangular facet has a single slope degree and direction. With a little imagination and practice, you can visualize the underlying surface from the TIN even without drawing contours.\r\n\r\nWonder why I suggest that you not let triangle sides that make up the TIN cross the stream? Well, if you did, the stream would appear to run along the side of a hill, instead of down a valley as it should. In practice, spot elevations would always be measured at several points along the stream, and along ridges as well. Photogrammetrists refer to spot elevations collected along linear features as\u00a0<strong>breaklines<\/strong>\u00a0(Maune, 2007). I omitted breaklines from this example just to make a point.\r\n\r\nYou may notice that there is more than one correct way to draw the TIN. As you will see, deciding which spot elevations are \u201cnear neighbors\u201d and which are not is subjective in some cases. Related to this element of subjectivity is the fact that the fidelity of a contour map depends in large part on the distribution of spot elevations on which it is based. In general, the density of spot elevations should be greater where terrain elevations vary greatly, and sparser where the terrain varies subtly. Similarly, the smaller the contour interval you intend to use, the more spot elevations you need.\r\n\r\n(There are algorithms for triangulating irregular arrays that produce unique solutions. One approach is called\u00a0<strong>Delaunay Triangulation<\/strong>which, in one of its constrained forms, is useful for representing terrain surfaces. The distinguishing geometric characteristic of a Delaunay triangulation is that a circle surrounding each triangle side does not contain any other vertex.)\r\n\r\n<img alt=\"Step 3 of contouring demonstration\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/contouring_4.gif\" \/>\r\n\r\n&nbsp;\r\n\r\nTick marks drawn where elevation contours cross the edges of each TIN facet.\r\n\r\n&nbsp;\r\n\r\nNow draw ticks to mark the points at which elevation contours intersect each triangle side. For instance, see the triangle side that connects the spot elevations 2360 and 2480 in the lower left corner of the illustration above? One tick mark is drawn on the triangle where a contour representing elevation 2400 intersects. Now find the two spot elevations, 2480 and 2750, in the same lower left corner. Note that three tick marks are placed where contours representing elevations 2500, 2600, and 2700 intersect.\r\n\r\nThis step should remind you of the equal interval classification scheme you read about in Chapter 3. The right choice of contour interval depends on the goal of the mapping project. In general, contour intervals increase in proportion to the variability of the terrain surface. It should be noted that the assumption that elevations increase or decrease at a constant rate is not always correct, of course. We will consider that issue in more detail later.\r\n\r\n<img alt=\"Step 4 of contouring demonstration\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/contouring_5.gif\" \/>\r\n\r\nThreading elevation contours through a TIN.\r\n\r\nFinally, draw your contour lines. Working downslope from the highest elevation, thread contours through ticks of equal value. Move to the next highest elevation when the surface seems ambiguous.\r\n\r\nKeep in mind the following characteristics of contour lines (Rabenhorst and McDermott, 1989):\r\n<ul>\r\n\t<li>Contours should\u00a0<strong>always point upstream<\/strong>\u00a0in valleys<\/li>\r\n\t<li>Contours should\u00a0<strong>always point downridge<\/strong>\u00a0along ridges<\/li>\r\n\t<li>Adjacent contours should<strong>\u00a0always be sequential or equivalent<\/strong><\/li>\r\n\t<li>Contours should\u00a0<strong>never split into two<\/strong><\/li>\r\n\t<li>Contours should\u00a0<strong>never cross or loop<\/strong><\/li>\r\n\t<li>Contours should\u00a0<strong>never spiral<\/strong><\/li>\r\n\t<li>Contours should\u00a0<strong>never stop in the middle of a map<\/strong><\/li>\r\n<\/ul>\r\nHow does your finished map compare with the one I drew below?\r\n\r\n<img alt=\"Outcome of contouring demonstration\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/contouring_6.gif\" \/>\r\n<h3><strong>TRY THIS<\/strong><strong>!<\/strong><\/h3>\r\n<strong>Now try your hand at contouring on your own.<\/strong>\u00a0The purpose of this practice activity is to give you more experience in contouring terrain surfaces.\r\n<ol>\r\n\t<li>First, view\u00a0<a href=\"http:\/\/opentextbc.ca\/natureofgeoinfo\/files\/geog482\/image\/contouring_practice-apr2012.gif\">an image of an irregular array of 16 spot elevations<\/a>.<\/li>\r\n\t<li>Print the image.<\/li>\r\n\t<li>Use the procedure outlined in this lesson to draw contour lines that represent the terrain surface that the spot elevations were sampled from. You may find this to be a moderately challenging task that takes about a half hour to do well. TIP: label the tick marks to make it easier to connect them.<\/li>\r\n\t<li>When finished,\u00a0<a href=\"https:\/\/www.e-education.psu.edu\/files\/geog482\/image\/mt_nittany.jpg\">compare your result to an existing map.<\/a><\/li>\r\n<\/ol>\r\nHere are a couple of somewhat simpler problems and solutions in case you need a little more practice.\r\n<ul>\r\n\t<li><a href=\"https:\/\/www.e-education.psu.edu\/files\/geog482\/image\/cont_practice_will1.gif\">Practice Problem #1<\/a><\/li>\r\n\t<li><a href=\"https:\/\/www.e-education.psu.edu\/files\/geog482\/image\/cont_practice_will6.gif\">Practice Problem #1 Solution<\/a><\/li>\r\n\t<li><a href=\"https:\/\/www.e-education.psu.edu\/files\/geog482\/image\/cont_practice_pitt1.gif\">Practice Problem #2<\/a><\/li>\r\n\t<li><a href=\"https:\/\/www.e-education.psu.edu\/files\/geog482\/image\/cont_practice_pitt6.gif\">Practice Problem #2 Solution<\/a><\/li>\r\n<\/ul>\r\nYou will be asked to demonstrate your contouring ability again in the Lesson 7 Quiz and in the final exam.\r\n\r\nKevin Sabo (personal communication, Winter 2002) remarked that \u201cIf you were unfortunate enough to be hand-contouring data in the 1960\u2032s and 70\u2032s, you may at least have had the aid of a\u00a0<a href=\"http:\/\/www.nzeldes.com\/HOC\/Gerber.htm\">Gerber Variable Scale<\/a>. After hand contouring in Lesson 7, I sure wished I had my Gerber!\u201d\r\n<h3><strong>PRACTICE\u00a0<\/strong><strong>QUIZ<\/strong><\/h3>\r\nRegistered Penn State students should return now to the Chapter 7 folder\u00a0 in ANGEL (via the Resources menu to the left) to take a self-assessment quiz about\u00a0<strong>Contouring<\/strong>. You may take practice quizzes as many times as you wish. They are not scored and do not affect your grade in any way.\r\n<h2>7.7. Digital Line Graph (DLG)<\/h2>\r\n<h3>IDENTIFICATION<\/h3>\r\nDigital Line Graphs (DLGs) are vector representations of most of the features and attributes shown on USGS topographic maps. Individual feature sets (outlined in the table below) are encoded in separate digital files. DLGs exist at three scales: small (1:2,000,000), intermediate (1:100,000) and large (1:24,000). Large-scale DLGs are produced in<strong>tiles<\/strong>\u00a0that correspond to the 7.5-minute topographic quadrangles from which they were derived.\r\n\r\n&nbsp;\r\n<table summary=\"Layers and contents of large-scale Digital Line Graph files\"><caption>Description of Digital Line Graph Layers<\/caption>\r\n<thead>\r\n<tr>\r\n<th><strong>Layer<\/strong><\/th>\r\n<th><strong>Features<\/strong><\/th>\r\n<\/tr>\r\n<\/thead>\r\n<tbody>\r\n<tr>\r\n<td>Public Land Survey System (PLSS)<\/td>\r\n<td>Township, range, and section lines<\/td>\r\n<\/tr>\r\n<tr>\r\n<td>Boundaries<\/td>\r\n<td>State, county, city, and other national and State lands such as forests and parks<\/td>\r\n<\/tr>\r\n<tr>\r\n<td>Transportation<\/td>\r\n<td>Roads and trails, railroads, pipelines and transmission lines<\/td>\r\n<\/tr>\r\n<tr>\r\n<td>Hydrography<\/td>\r\n<td>Flowing water, standing water, and wetlands<\/td>\r\n<\/tr>\r\n<tr>\r\n<td>Hypsography<\/td>\r\n<td>Contours and supplementary spot elevations<\/td>\r\n<\/tr>\r\n<tr>\r\n<td>Non-vegetative features<\/td>\r\n<td>Glacial moraine, lava, sand, and gravel<\/td>\r\n<\/tr>\r\n<tr>\r\n<td>Survey control and markers<\/td>\r\n<td>Horizontal and vertical monuments (third order or better)<\/td>\r\n<\/tr>\r\n<tr>\r\n<td>Man-made features<\/td>\r\n<td>Cultural features, such as building, not collected in other data categories<\/td>\r\n<\/tr>\r\n<tr>\r\n<td>Woods, scrub, orchards, and vineyards<\/td>\r\n<td>Vegetative surface cover<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\nLayers and contents of large-scale Digital Line Graph files. Not all layers available for all quadrangles (USGS, 2006).\r\n\r\n&nbsp;\r\n\r\n&nbsp;\r\n\r\n<img alt=\"Image of Digital Line Graph hyspography, hydrography, and transportation layers viewed in Global Mapper software\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/dlg_bushkill_sw.gif\" \/>\r\n\r\nPortion of three Digital Line Graph (DLG) layers for USGS Bushkill, PA quadrangle; imaged with Global Mapper (dlgv32 Pro) software. Transportation features are arbitrarily colored red, hydrography blue, and hypsography brown. The square symbols are nodes and the triangles represent polygon centroids.\r\n<h4>DATA QUALITY<\/h4>\r\nLike other USGS data products, DLGs conform to National Map Accuracy Standards. In addition, however, DLGs are tested for the logical consistency of the topological relationships among data elements. Similar to the Census Bureau\u2019s TIGER\/Line, line segments in DLGs must begin and end at point features (nodes), and line segments must be bounded on both sides by area features (polygons).\r\n<h4>SPATIAL REFERENCE INFORMATION<\/h4>\r\nDLGs are heterogenous. Some use UTM coordinates, others State Plane Coordinates. Some are based on NAD 27, others on NAD 83. Elevations are referenced either to NGVD 29 or NAVD 88 (USGS, 2006a).\r\n<h4>ENTITIES AND ATTRIBUTES<\/h4>\r\nThe basic elements of DLG files are nodes (positions), line segments that connect two nodes, and areas formed by three or more line segments. Each node, line segment, and area is associated with two-part integer attribute codes. For example, a line segment associated with the attribute code \u201c050 0412\u2033 represents a hydrographic feature (050), specifically, a stream (0412).\r\n<h4>DISTRIBUTION<\/h4>\r\nNot all DLG layers are available for all areas at all three scales. Coverage is complete at 1:2,000,000.\u00a0 At the intermediate scale, 1:100,000 (30 minutes by 60 minutes), all hydrography and transportation files are available for the entire U.S., and complete national coverage is planned. At 1:24,000 (7.5 minutes by 7.5 minutes), coverage remains spotty. The files are in the public domain, and can be used for any purpose without restriction.\r\n\r\nLarge- and Intermediate -scale DLGs are available for download through\u00a0<a href=\"http:\/\/earthexplorer.usgs.gov\/\">EarthExplorer system<\/a>. You can plot 1:2,000,000 DLGs on-line at the\u00a0<a href=\"http:\/\/nationalatlas.gov\/\">USGS\u2019 National Atlas of the United States<\/a>.\r\n<h3><\/h3>\r\n<h3>DIGITAL LINE GRAPH HYPSOGRAPHY<\/h3>\r\nIn one sense, DLGs are as much \u201clegacy\u201d data as the out-of-date topographic maps from which they were produced. Still, DLG data serve as primary or secondary sources for several themes in the USGS National Map, including hydrography, boundaries, and transportation. DLG hypsography data are not included in the National Map, however. It is assumed that GIS users can generate elevation contours as needed from DEMs. DLG hypsography and hydrography layers are the preferred sources from which USGS DEMs are produced, however.\r\n\r\n<img alt=\"Digital Line Graph hypsography and hydrography viewed in Global Mapper software\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/dlg_hypso.gif\" \/>\r\n\r\nPortion of the hypsography and hydrography layers of a large-scale Digital Line Graph (DLG). USGS Bushkill, PA quadrangle; imaged with Global Mapper (dlgv32 Pro) software.\r\n\r\n<strong>Hypsography<\/strong>\u00a0refers to the measurement and depiction of the terrain surface, specifically with contour lines. Several different methods have been used to produce DLG hypsography layers, including:\r\n<ul>\r\n\t<li>Scanning contour lines on photographic film or paper maps, converting the scanned raster data to vectors, then editing and attributing the vector features;<\/li>\r\n\t<li>Manually digitizing and attributing contour lines on photographic film or paper maps; and<\/li>\r\n\t<li>Producing contours by photogrammetric processes.<\/li>\r\n<\/ul>\r\nThe preferred method is to manually digitize contour lines in vector mode, then to key-enter the corresponding elevation attribute data.\r\n\r\n<img alt=\"Attributes of a contour line in a DLG hypsography layer, viewed in Global Mapper software\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/dlg_hypso_attribute.gif\" \/>\r\n\r\nThe highlighted contour line has been selected, and its attributes reported in a Global Mapper window. Notice that the line feature is attributed with a unique Element ID code (LE01, 639) and an elevation (1000 feet).\r\n\r\n&nbsp;\r\n<h3><strong>TRY THIS!<\/strong><\/h3>\r\n<h3><strong>EXPLORING DLGS WITH GLOBAL MAPPER (DLGV32 PRO)<\/strong><\/h3>\r\nNow I\u2019d like you to use Global Mapper (or dlgv32 Pro) software to investigate the characteristics of the hypsography layer of a USGS Digital Line Graph (DLG). The instructions below assume that you have already installed software on your computer. (If you haven\u2019t, return to\u00a0<a href=\"https:\/\/www.e-education.psu.edu\/natureofgeoinfo\/c6_p6.html\">installation instructions<\/a>\u00a0presented earlier in Chapter 6). First you\u2019ll download and a sample DLG file. In a following activity you\u2019ll have a chance to find and download DLG data for your area.\r\n<ol>\r\n\t<li>If you haven\u2019t done so already, create a directory called \u201cUSGS Data\u201d on your hard disk, where you file your course materials.<\/li>\r\n\t<li>Next,\u00a0<a href=\"https:\/\/courseware.e-education.psu.edu\/downloads\/natureofgeoinfo\/DLG.zip\">Download the\u00a0<strong>DLG.zip<\/strong>\u00a0data archive<\/a>. The ZIP archive is 1.2 Mb in size and will take approximately 15 seconds to download via high speed DSL or cable, or about 4 minutes and 15 seconds minutes via 56 Kbps modem.<\/li>\r\n\t<li>Now decompress the archive into a directory on your hard disk.\r\n<ul>\r\n\t<li>Open the archive DLG.zip.<\/li>\r\n\t<li>Create a subdirectory called \u201cDLG\u201d within the directory in which you save data for this class.<\/li>\r\n\t<li>Extract all files in the ZIP archive into your new subdirectory.<\/li>\r\n<\/ul>\r\nThe end result will be five subdirectories, each of which includes the data files that make up a DLG \u201clayer,\u201d along with a master directory.<\/li>\r\n\t<li>Launch Global Mapper or dlgv32 Pro.<\/li>\r\n\t<li>Open a Digital Line Graph by choosing File &gt; Open as New\u2026, then navigate to the directory \u201cDLG\/Hypso.\u201d Open the file \u2018Hp01catd.ddf\u2019 (you can open up to four files at once in the trial version of Global Mapper.) The data correspond with the 7.5 minute quadrangle for Bushkill, PA. The file is encoded in Spatial Data Transfer Standard (SDTS) format. For information about SDTS, see the\u00a0<a href=\"https:\/\/www.e-education.psu.edu\/files\/geog482\/file\/sdts-tutorial.pdf\">SDTS Tutorial<\/a>(PDF format).<\/li>\r\n\t<li>Global Mapper may ask you to direct it to a \u2018Master Data Dictionary\u2019 file. If so, navigate to, and select, the file \u2018Dlg\/MasterDlg\/Dlg3mdir.ddf\u2019<\/li>\r\n\t<li>Experiment with Global Mapper\u2019s tools. Use Zoom and Pan to magnify and scroll across the DLG. The\u00a0<strong>Full View button<\/strong>\u00a0(the one with the house icon) refreshes the initial full view of the data set.<\/li>\r\n\t<li>The<strong>\u00a0Feature Info tool<\/strong>\u00a0allows you to query the attributes of a particular feature. Try clicking a single line segment. Note that you can display the attributes of a feature in the lower left portion of the application window by simply hovering over the feature.<\/li>\r\n\t<li>The\u00a0<strong>Measure tool<\/strong>\u00a0(ruler icon) allows you to not only measure distance as the crow flies, but also to see the area enclosed by a series of line segments drawn by repeated mouse clicks. Note again the location information that is given to you near the bottom of the application window.<\/li>\r\n\t<li>Certain tools, e.g., the\u00a0<strong>3D Path Profile\/Line of Sight tool<\/strong>\u00a0(next to the Feature Info tool) are not functional in the free (unregistered) version of Global Mapper.<\/li>\r\n\t<li>The trial version of Global Mapper allows you to open and view up to four files at once. You might find it interesting to open and compare the Bushkill DLG hypsography file and the corresponding DRG you viewed in Lesson 6. Note that you can turn layers on and off, and even adjust their transparency at Tools &gt; Control Center.\u00a0<strong>How do the contours in the DLG compare with those in the DRG? What explains the difference?<\/strong><\/li>\r\n\t<li>Global Mapper provides the metadata you\u2019ll need to answer questions in a practice quiz.\u00a0<strong>To access the metadata, navigate to Tools &gt; Control Center, then click the Metadata button.<\/strong><\/li>\r\n<\/ol>\r\n<h2>7.8. Digital Elevation Model (DEM)<\/h2>\r\nThe term \u201cDigital Elevation Model\u201d has both generic and specific meanings. In general, a DEM is any raster representation of a terrain surface. Specifically, a DEM is a data product of the U.S. Geological Survey. Here we consider the characteristics of DEMs produced by the USGS Later in this chapter we\u2019ll consider sources of global terrain data.\r\n<h3>IDENTIFICATION<\/h3>\r\nUSGS DEMs are raster grids of elevation values that are arrayed in series of south-north profiles. Like other USGS data, DEMs were produced originally in tiles that correspond to topographic quadrangles. Large scale (7.5-minute and 15-minute), intermediate scale (30 minute), and small scale (1 degree) series were produced for the entire U.S. The resolution of a DEM is a function of the east-west spacing of the profiles and the south-north spacing of elevation points within each profile.\r\n\r\nDEMs corresponding to 7.5-minute quadrangles are available at 10-meter resolution for much, but not all, of the U.S. Coverage is complete at 30-meter resolution. In these large scale DEMs elevation profiles are aligned parallel to the central meridian of the local UTM zone, as shown in the illustration below.\u00a0<strong>See how the DEM tile in the illustration below appears to be tilted? This is because the corner points are defined in unprojected geographic coordinates that correspond to the corner points of a USGS quadrangle<\/strong>. The farther the quadrangle is from the central meridian of the UTM zone, the more it is tilted.\r\n\r\n&nbsp;\r\n\r\n<img alt=\"Schematic illustration of elevation profiles in a USGS digital elevation model\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/dem_structure.png\" \/>\r\n\r\nArrangement of elevation profiles in a large scale USGS Digital Elevation Model (USGS, 1987).\r\n\r\n&nbsp;\r\n\r\nAs shown below, the arrangement of the elevation profiles is different in intermediate- and small-scale DEMs. Like meridians in the northern hemisphere, the profiles in 30-minute and 1-degree DEMs converge toward the north pole. For this reason the resolution of intermediate- and small-scale DEMs (that is to say, the spacing of the elevation values) is expressed differently than for large-scale DEMs.\u00a0<strong>The resolution of 30-minute DEMs is said to be 2 arc seconds and 1-degree DEMs are 3 arc seconds<\/strong>. Since an arc second is 1\/3600 of a degree, elevation values in a 3 arc second DEM are spaced 1\/1200 degree apart, representing a grid cell about 66 meters \u201cwide\u201d by 93 meters \u201ctall\u201d at 45\u00ba latitude.\r\n\r\n<img alt=\"Schematic illustration of elevation profiles in a small-scale USGS digital elevation model\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/dem_structure_smallscale.png\" \/>\r\n\r\nArrangement of elevation profiles in a small scale USGS Digital Elevation Model (USGS, 1987).\r\n\r\n<strong>The preferred method for producing the elevation values that populate DEM profiles is interpolation from DLG hypsography and hydrography layers<\/strong>\u00a0(including the hydrography layer enables analysts to delineate valleys with less uncertainty than hypsography alone). Some older DEMs were produced from elevation contours digitized from paper maps or during photogrammetric processing, then smoothed to filter out errors. Others were produced photogrammtrically from aerial photographs.\r\n<h4>DATA QUALITY<\/h4>\r\nThe vertical accuracy of DEMs is expressed as the root mean square error (RMSE) of a sample of at least 28 elevation points. The target accuracy for large-scale DEMs is seven meters; 15 meters is the maximum error allowed.\r\n<h4>SPATIAL REFERENCE INFORMATION<\/h4>\r\nLike DLGs, USGS DEMs are heterogenous. They are cast on the Universal Transverse Mercator projection used in the local UTM zone. Some DEMs are based upon the North American Datum of 1983, others on NAD 27. Elevations in some DEMs are referenced to either NGVD 29 or NAVD 88.\r\n<h4>ENTITIES AND ATTRIBUTES<\/h4>\r\nEach record in a DEM is a profile of elevation points. Records include the UTM coordinates of the starting point, the number of elevation points that follow in the profile, and the elevation values that make up the profile. Other than the starting point, the positions of the other elevation points need not be encoded, since their spacing is defined. (Later in this lesson you\u2019ll download a sample USGS DEM file. Try opening it in a text editor to see what I\u2019m talking about.)\r\n<h4>DISTRIBUTION<\/h4>\r\nDEM tiles are available for free download through many state and regional clearinghouses. You can find these sources by searching the<a href=\"http:\/\/geo.data.gov\/geoportal\/catalog\/main\/home.page\">GEODATA portion of the Data.Gov site<\/a>, formerly the separate Geospatial One Stop site.\r\n\r\nAs part of its National Map initiative, the USGS has developed a \u201cseamless\u201d\u00a0<a href=\"http:\/\/ned.usgs.gov\/\">National Elevation Dataset<\/a>\u00a0that is derived from DEMs, among other sources. NED data are available at three resolutions: 1 arc second (approximately 30 meters), 1\/3 arc second (approximately 10 meters), and 1\/9 arc second (approximately 3 meters). Coverage ranges from complete at 1 arc second to extremely sparse at 1\/9 arc second. An extensive\u00a0<a href=\"http:\/\/ned.usgs.gov\/Ned\/faq.asp\">FAQ on NED data is published here<\/a>. The second of the two following activities involves downloading NED data and viewing it in Global Mapper.\r\n<h3><strong>TRY THIS!<\/strong><\/h3>\r\n<h3><strong>EXPLORING DEMS WITH GLOBAL MAPPER (DLGV32 PRO)<\/strong><\/h3>\r\nGlobal Mapper time again! This time you\u2019ll investigate the characteristics of a USGS DEM. The instructions below assume that you have already installed the software on your computer. (If you haven\u2019t, return to\u00a0<a href=\"https:\/\/www.e-education.psu.edu\/natureofgeoinfo\/c6_p6.html\">installation instructions<\/a>\u00a0presented earlier in Chapter 6). The instructions will remind you how to open a DEM in dlgv32 Pro. In the practice quiz that follows you\u2019ll be asked questions require you to explore the data for answers.\r\n<ol>\r\n\t<li>First\u00a0<a href=\"https:\/\/courseware.e-education.psu.edu\/downloads\/natureofgeoinfo\/DEM.zip\">Download the\u00a0<strong>DEM.zip<\/strong>\u00a0data archive<\/a>. The ZIP archive is 2.5 Mb in size and will take about 30 seconds to download via high speed DSL or cable, or nearly 9 minutes via 56 Kbps modem. If you can\u2019t download the file, contact my teaching assistant or me right away so we can help you resolve the problem.<\/li>\r\n\t<li>Now decompress the archive into a directory on your hard disk.\r\n<ul>\r\n\t<li>Open the archive<strong>\u00a0DEM.zip.<\/strong><\/li>\r\n\t<li>Create a subdirectory called \u201cDEM\u201d within the directory in which you save class data.<\/li>\r\n\t<li>Extract all files in the ZIP archive into your new subdirectory.<\/li>\r\n<\/ul>\r\nThe end result will be two subdirectories, one of which contains a 30-meter DEM, the other a 10-meter DEM. These datasets are in the earlier distribution format of USGS DEM data \u2014 elevation data in horizontal (pixel) units of meters and representative of the area covered by a 1:24,000 topo map sheet. In the Try This that follows this one you will see that the distribution format options have expanded.<\/li>\r\n\t<li>Launch Global Mapper.<\/li>\r\n\t<li>Open a Digital Elevation Model by choosing\u00a0<strong>File\u00a0<\/strong>&gt;\u00a0<strong>Open Data File(s)\u2026<\/strong>, then navigate to the directory DEM_30m or DEM_10m, then open the file\u00a0<em>bushkill_pa.dem<\/em><\/li>\r\n\t<li>Use the\u00a0<strong>Zoom\u00a0<\/strong>and\u00a0<strong>Pan\u00a0<\/strong>tools to magnify and scroll across the DEM. The<strong>\u00a0Full View<\/strong>\u00a0button (house icon) refreshes the initial full view of the dataset.<\/li>\r\n\t<li>Global Mapper provides access to the metadata you\u2019ll need to answer questions in a practice quiz. To access the metadata, navigate to<strong>\u00a0Tools\u00a0<\/strong>&gt;<strong>\u00a0Control Center<\/strong>, then click the<strong>\u00a0Metadata\u00a0<\/strong>button.<\/li>\r\n<\/ol>\r\nYou can change the appearance of the DEM in the Options section of the Control Center. You can also alter the appearance of the DEM by choosing Tools &gt; Configure, and changing the settings in, especially, Vertical Options and Shader Options. To see the DEM data with(out) hill shading, find the Enable\/Disable Hill Shading button on the Shader toolbar (it has a sunburst in the lower left corner).\r\n<h3><strong>TRY THIS!<\/strong><\/h3>\r\n<h3><strong>DOWNLOAD YOUR OWN NATIONAL ELEVATION DATASET (NED) DATA<\/strong><\/h3>\r\n<ol>\r\n\t<li>Go to the\u00a0<a href=\"http:\/\/nationalmap.gov\/elevation.html\">Elevation page of the USGS National Map site<\/a>.\r\nRead any of the information there that you choose.<\/li>\r\n\t<li>Follow the link to\u00a0<a href=\"http:\/\/viewer.nationalmap.gov\/viewer\/\">The National Map Viewer<\/a>.\r\nIf you wish, also follow the link to the detailed instructions that you see under The National Map Viewer link. The following instructions in this Try This should also suffice, and perhaps elaborate at bit more.You will not see Elevation data listed in the left hand Overlays pane. The USGS is in the process of adding more visualization options when it comes to elevation data. See the Hill Shade button at the upper right of the map area.<\/li>\r\n\t<li>Use the GIS tools, found above the map area, to pan to and zoom in on an area of interest. Then, click the\u00a0<strong>Download Data<\/strong>\u00a0button.\r\nChoose a\u00a0<em>reference area<\/em>\u00a0from the pick list. The default is the index of the areas covered by the 1:24,000 topo map series.\r\nYou could also chose the entire\u00a0<em>current map extent<\/em>, but depending upon your zoom level that may be a huge dataset.\r\nThe instructions found via the Help link mentioned above mention that you can define an area based on creating a\u00a0<em>custom polygon<\/em>, but at this point I do not see how that is done\u2026<\/li>\r\n\t<li>Then click on the map to highlight a specific area of interest. A link to the available data sets will appear in the left hand pane under the<strong>Selection\u00a0<\/strong>tab. (The\u00a0<strong>All Results<\/strong>\u00a0button will list the multiple areas you click on.) Follow the\u00a0<strong>Download\u00a0<\/strong>link of the area you are interested in.<\/li>\r\n\t<li>In the\u00a0<strong>USGS Available Data<\/strong>\u00a0window that opens, select\u00a0<strong>Elevation<\/strong>from the\u00a0<strong>Theme<\/strong>\u00a0column, and choose a file format from the pick list in the\u00a0<strong>Format\u00a0<\/strong>column.\r\nI know that both the GeoTIFF and ArcGRID formats are compatible with Global Mapper. (ArcGRID is the Esri company\u2019s raster format. You\u2019ll be using the Esri software in future courses.)\r\nClick the\u00a0<strong>Next\u00a0<\/strong>button.<\/li>\r\n\t<li>You will be given a list of Elevation\u00a0<em>Product\u00a0<\/em>choices.\r\nThe choices that show\u00a0<em>Dynamic\u00a0<\/em>as the Type will be those that match the Format choice you made in the previous window. The other\u00a0<em>Staged<\/em>datasets are prepackaged and in the formats stated in the Product column, and apparently listed regardless of the Format choice you made.If multiple resolutions are available they will be listed.\r\n\r\nFrom the Product list, check the box for what you wish to download.\r\nClick the\u00a0<strong>Next\u00a0<\/strong>button.\r\nGo to the\u00a0<strong>Cart\u00a0<\/strong>pane on the right. (It may open automatically.)\r\nGo to\u00a0<strong>Checkout<\/strong>, supply your e-mail address, and submit your order via the\u00a0<strong>Place Order<\/strong>\u00a0button.\r\nYou will receive a message telling you that your order has been placed, and that will soon be followed by an e-mail regarding your order. About an hour and a half after I submitted my request I received a second e-mail containing a download link.<\/li>\r\n\t<li>The system will produce a ZIP archive that you can save to your hard disk (e.g., \u201c09647011.zip\u201d)<\/li>\r\n\t<li>Launch Global Mapper and open the ZIP archive. The software can read the data even in its compressed form; you should not need to extract the contents from the .zip file. (It would be a good idea to look at the contents of the .zip archive, though, if only to see the number and type of files included.)<\/li>\r\n\t<li>An image of the DEM data should appear in the Global Mapper window, similar to what you see shown below (even though the image below is from an older version of Global Mapper).<\/li>\r\n\t<li>Again, you can view the\u00a0<em>metadata\u00a0<\/em>associated with the DEM data via the\u00a0<strong>Tools\u00a0<\/strong>&gt;\u00a0<strong>Control Center<\/strong>\u00a0menu. Note the PIXEL dimensions reported in\u00a0<em>arc degrees<\/em>, as opposed to something like meters.<\/li>\r\n<\/ol>\r\n<img alt=\"A portion of the National Elevation Dataset viewed in Global Mapper software\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/ned_global_mapper.png\" \/>\r\n<h3><strong>PRACTICE QUIZ<\/strong><\/h3>\r\nRegistered Penn State students should return now to the Chapter 7 folder\u00a0 in ANGEL (via the Resources menu to the left) to take a self-assessment quiz about\u00a0<strong>DLGs and DEMs<\/strong>. You may take practice quizzes as many times as you wish. They are not scored and do not affect your grade in any way.\r\n<h2>7.9. Interpolation<\/h2>\r\nDEMs are produced by various methods. The method preferred by USGS is to\u00a0<strong>interpolate<\/strong>\u00a0elevations grids from the hypsography and hydrography layers of Digital Line Graphs.\r\n\r\n<img alt=\"Hypsography and hydrography layers of a Digital Line Graph viewed in Global Mapper software\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/dem_dlg.jpg\" \/>\r\n\r\nA USGS 7.5-minute DEM and the DLG hypsography and hydrography layers from which it was produced.\r\n\r\nThe elevation points in DLG hypsography files are not regularly spaced. DEMs need to be regularly spaced to support the slope, gradient, and volume calculations they are often used for. Grid point elevations must be interpolated from neighboring elevation points. In the figure below, for example, the gridded elevations shown in purple were interpolated from the irregularly spaced spot elevations shown in red.\r\n\r\n<img alt=\"A grid of elevation values that were interpolated from an irregularly-spaced array\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/interp_grid.gif\" \/>\r\n\r\nElevation values in DEMs are interpolated from irregular arrays of elevations measured through photogrammetric methods, or derived from existing DLG hypsography and hydrography data.\r\n\r\nHere\u2019s another example of interpolation for mapping. The map below shows how 1995 average surface air temperature differed from the average temperature over a 30-year baseline period (1951-1980). The temperature anomalies are depicted for grid cells that cover 3\u00b0 longitude by 2.5\u00b0 latitude.\r\n\r\n<img alt=\"A map showing gridded temperature data\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/temp_95.jpg\" \/>\r\n\r\n1995 Surface Temperature Anomalies. (National Climatic Data Center, 2005).\r\n\r\nThe gridded data shown above were estimated from the temperature records associated with the very irregular array of 3,467 locations pinpointed in the map below.\u00a0 The irregular array is transformed into a regular array through interpolation. In general,\u00a0<strong>interpolation is the process of estimating an unknown value from neighboring known values.<\/strong>\r\n\r\n<img alt=\"Locations of temperature climate records used to create a gridded temperature map\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/temp_stations.jpg\" \/>\r\n\r\nThe Global Historical Climate Network. (Eischeid et al., 1995).\r\n\r\nElevation data are often not measured at evenly-spaced locations. Photogrammetrists typically take more measurements where the terrain varies the most. They refer to the dense clusters of measurements they take as \u201cmass points.\u201d Topographic maps (and their derivatives, DLGs) are another rich source of elevation data. Elevations can be measured from contour lines, but obviously contours do not form evenly-spaced grids. Both methods give rise to the need for interpolation.\r\n\r\n<img alt=\"Three number lines illustrating how interpolation is affected by assumptions about the underlying distribution\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/linear_interp.gif\" \/>\r\n\r\n&nbsp;\r\n\r\nInterpolating an intermediate value on a number line.\r\n\r\nThe illustration above shows three number lines, each of which ranges in value from 0 to 10. If you were asked to interpolate the value of the tick mark labeled \u201c?\u201d on the top number line, what would you guess? An estimate of \u201c5\u2033 is reasonable, provided that the values between 0 and 10 increase at a constant rate. If the values increase at a geometric rate, the actual value of \u201c?\u201d could be quite different, as illustrated in the bottom number line. The validity of an interpolated value depends, therefore, on the validity of our assumptions about the nature of the underlying surface.\r\n\r\nAs I mentioned in Chapter 1, the surface of the Earth is characterized by a property called\u00a0<strong>spatial dependence<\/strong>. Nearby locations are more likely to have similar elevations than are distant locations. Spatial dependence allows us to assume that it\u2019s valid to estimate elevation values by interpolation.\r\n\r\nMany interpolation algorithms have been developed. One of the simplest and most widely used (although often not the best) is the<strong>inverse distance weighted<\/strong>\u00a0algorithm. Thanks to the property of spatial dependence, we can assume that estimated elevations are more similar to nearby elevations than to distant elevations. The inverse distance weighted algorithm estimates the value z of a point P as a function of the z-values of the nearest n points. The more distant a point, the less it influences the estimate.\r\n\r\n<img alt=\"Diagram and formula explaining inverse distance weighted interpolation\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/inv_dist_interp.gif\" \/>\r\n\r\nThe inverse distance weighted interpolation procedure.\r\n<h3><strong>PRACTICE QUIZ<\/strong><\/h3>\r\nRegistered Penn State students should return now to the Chapter 7 folder\u00a0 in ANGEL (via the Resources menu to the left) to take a self-assessment quiz about\u00a0<strong>Interpolation<\/strong>. You may take practice quizzes as many times as you wish. They are not scored and do not affect your grade in any way.\r\n<h2>7.10. Slope<\/h2>\r\nSlope is a measure of change in elevation. It is a crucial parameter in several well-known predictive models used for environmental management, including the Universal Soil Loss Equation and agricultural non-point source pollution models.\r\n\r\nOne way to express slope is as a percentage. To calculate percent slope, divide the difference between the elevations of two points by the distance between them, then multiply the quotient by 100. The difference in elevation between points is called the rise. The distance between the points is called the run. Thus, percent slope equals (rise \/ run) x 100.\r\n\r\n<img alt=\"Diagram illustrating how slope may be calculated as a percentage\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/slope_percent.gif\" \/>\r\n\r\nCalculating percent slope. A rise of 100 feet over a run of 100 feet yields a 100 percent slope. A 50-foot rise over a 100-foot run yields a 50 percent slope.\r\n\r\n&nbsp;\r\n\r\nAnother way to express slope is as a slope angle, or degree of slope. As shown below, if you visualize rise and run as sides of a right triangle, then the degree of slope is the angle opposite the rise. Since degree of slope is equal to the tangent of the fraction rise\/run, it can be calculated as the arctangent of rise\/run.\r\n\r\n<img alt=\"Illustration showing how slope may be calculated in degrees\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/slope_degree.gif\" \/>\r\n\r\nA rise of 100 feet over a run of 100 feet yields a 45\u00b0 slope angle. A rise of 50 feet over a run of 100 feet yields a 26.6\u00b0 slope angle.\r\n\r\n&nbsp;\r\n\r\nYou can calculate slope on a contour map by analyzing the spacing of the contours. If you have many slope values to calculate, however, you will want to automate the process. It turns out that slope calculations are much easier to calculate for gridded elevation data than for vector data, since elevations are more or less equally spaced in raster grids.\r\n\r\nSeveral algorithms have been developed to calculate percent slope and degree of slope. The simplest and most common is called the<strong>neighborhood method<\/strong>. The neighborhood method calculates the slope at one grid point by comparing the elevations of the eight grid points that surround it.\r\n\r\n<img alt=\"Illustration showing how slope at one point is calculated as a function of the elevations of eight surrounding points\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/slope_grid.gif\" \/>\r\n\r\nThe neighborhood algorithm estimates percent slope in cell 5 by comparing the elevations of neighboring grid cells.\r\n\r\nThe neighborhood algorithm estimates percent slope at grid cell 5 (Z5) as the sum of the absolute values of east-west slope and north-south slope, and multiplying the sum by 100. The diagram below illustrates how east-west slope and north-south slope are calculated. Essentially, east-west slope is estimated as the difference between the sums of the elevations in the first and third columns of the 3 x 3 matrix. Similarly, north-south slope is the difference between the sums of elevations in the first and third rows (note that in each case the middle value is weighted by a factor of two).\r\n\r\n<img alt=\"Algorithm for calculating slope with gridded elevation data\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/slope_algorithm.gif\" \/>\r\n\r\nThe neighborhood algorithm for calculating percent slope.\r\n\r\nThe neighborhood algorithm calculates slope for every cell in an elevation grid by analyzing each 3 x 3 neighborhood. Percent slope can be converted to slope degree later. The result is a grid of slope values suitable for use in various soil loss and hydrologic models.\r\n<h2>7.11. Relief Shading<\/h2>\r\nYou can see individual pixels in the zoomed image of a 7.5-minute DEM below. I used dlgv32 Pro\u2019s \u201cGradient Shader\u201d to produce the image. Each pixel represents one elevation point. The pixels are shaded through 256 levels of gray. Dark pixels represent low elevations, light pixels represent high ones.\r\n\r\n<img alt=\"Image of a DEM with pixels shaded light to dark in proportion to elevation\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/dem.gif\" \/>\r\n\r\nA digital elevation model in which light pixels represent high elevations, and dark pixels represent low elevations.\r\n\r\nIt\u2019s also possible to assign gray values to pixels in ways that make it appear that the DEM is illuminated from above. The image below, which shows the same portion of the Bushkill DEM as the image above, illustrates the effect, which is called\u00a0<strong>terrain shading<\/strong>,\u00a0<strong>hill shading<\/strong>\u00a0or\u00a0<strong>shaded relief.<\/strong>\r\n\r\n<strong><img alt=\"Image of a DEM shaded as though it were illuminated from above\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/dem_shaded.jpg\" \/><\/strong>\r\n\r\nShaded terrain image produced from the same DEM as shown in the above figures, using dlgv32 Pro\u2019s Daylight Shader option, with the Surface Color set to gray.\r\n\r\nThe appearance of a shaded terrain image depends on several parameters, including\u00a0<strong>vertical exaggeration<\/strong>. Click the buttons under the image below to compare the four terrain images of North America shown below, in which elevations are exaggerated 5 times, 10 times, 20 times, and 40 times respectively. (You will need to have the Adobe Flash player installed in order to complete this exercise. If you do not already have the Flash player, you can\u00a0<a href=\"http:\/\/www.adobe.com\/shockwave\/download\/index.cgi?P1_Prod_Version=ShockwaveFlash\">download it for free from Adobe<\/a>.)\r\n\r\nEffects of vertical exaggeration on a shaded terrain image\r\n\r\nAnother influential parameter is the\u00a0<strong>angle of illumination<\/strong>. Click the buttons to compare terrain images that have been illuminated from the northeast, southeast, southwest, and northwest. Does the terrain appear to be inverted in one or more of the images?\u00a0<strong>To minimize the possibility of terrain inversion, it is conventional to illuminate terrain from the northwest.<\/strong>\r\n\r\nEffects of illumination angle on a shaded terrain image.\r\n<h2>7.12. Lidar<\/h2>\r\nFor many applications, 30-meter DEMs whose vertical accuracy is measured in meters are simply not detailed enough. Greater accuracy and higher horizontal resolution can be produced by photogrammetric methods, but precise photogrammetry is often too time-consuming and expensive for extensive areas. Lidar is a digital remote sensing technique that provides an attractive alternative.\r\n\r\nLidar stands for\u00a0<strong>LIght Detection And Ranging<\/strong>. Like radar (RAdio Detecting And Ranging), lidar instruments transmit and receive energy pulses, and enable distance measurement by keeping track of the time elapsed between transmission and reception. Instead of radio waves, however, lidar instruments emit laser light (laser stands for Light Amplifications by Stimulated Emission of Radiation).\r\n\r\nLidar instruments are typically mounted in low altitude aircraft. They emit up to 5,000 laser pulses per second, across a ground swath some 600 meters wide (about 2,000 feet). The ground surface, vegetation canopy, or other obstacles reflect the pulses, and the instrument\u2019s receiver detects some of the backscatter. Lidar mapping missions rely upon GPS to record the position of the aircraft, and upon inertial navigation instruments (gyroscopes that detect an aircraft\u2019s pitch, yaw, and roll) to keep track of the system\u2019s orientation relative to the ground surface.\r\n\r\nIn ideal conditions, lidar can produce DEMs with 15-centimeter vertical accuracy, and horizontal resolution of a few meters. Its cost is prohibitive for small missions, but is justified for larger projects in which detail is essential. For example, lidar has been used successfully to detect subtle changes in the thickness of the Greenland ice sheet that result in a net loss of over 50 cubic kilometers of ice annually.\r\n\r\n<img alt=\"Composite lidar image showing changes in thickness of Greenland ice sheet\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/lidar_greenland.jpg\" \/>\r\n\r\nImage of Greenland, viewed from the south, showing changes in ice thickness measured by airborne lidar. Ice sheet thickness decreasing at 40-60 cm per year in darker blue areas (Goddard Space Flight Center, n.d.).\r\n\r\nTo learn more about the use of lidar in mapping changes in the Greenland ice sheet, visit\u00a0<a href=\"http:\/\/svs.gsfc.nasa.gov\/stories\/greenland\/\">NASA\u2019s Scientific Visualization Studio<\/a>.\r\n<h2>7.13. Global Elevation Data<\/h2>\r\nThis page profiles three data products that include elevation (and, in one case, bathymetry) data for all or most of the Earth\u2019s surface.\r\n<h3>ETOPO1<\/h3>\r\n<img alt=\"World map generated from ETOPO1 global terrain (with ice heights) and bathymetry data\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/etopo1.gif\" \/>\r\nShaded and colored terrain image produced from ETOPO1 data. (National Geophysical Data Center, 2009).\r\n\r\nETOPO1 is a digital elevation model that includes both topography and bathymetry for the entire world. It consists of more than 233 million elevation values which are regularly spaced at 1 minute of latitude and longitude.\u00a0<strong>At the equator, the horizontal resolution of ETOPO1 is approximately 1.85 kilometers.\u00a0<\/strong>Vertical positions are specified in meters, and there are two versions of the dataset: one with elevations at the \u201cIce Surface\u201d of the Greenland and Antarctic ice sheets, and one with elevations at \u201cBedrock\u201d beneath those ice sheets. Horizontal positions are specified in geographic coordinates (decimal degrees). Source data, and thus data quality, vary from region to region.\r\nYou can\u00a0<a href=\"http:\/\/www.ngdc.noaa.gov\/mgg\/global\/global.html\">download ETOPO1 data from the National Geophysical Data Center<\/a>.\r\n<h3><\/h3>\r\n<h3>GTOPO30<\/h3>\r\n<img alt=\"World terrain map generated from GTOPO30 data\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/gtopo30_tiles.gif\" \/>\r\n\r\nShaded and colored terrain image produced from GTOPO30 data. Data are distributed as 33 tiles (USGS, 2006b).\r\n\r\nGTOPO30 is a digital elevation model that extends over the world\u2019s land surfaces (but not under the oceans). GTOPO30 consists of more than 2.5 million elevation values, which are regularly spaced at 30 seconds of latitude and longitude.<strong>\u00a0At the equator, the resolution of GTOPO30 is approximately 0.925 kilometers \u2014 two times greater than ETOPO1<\/strong>. Vertical positions are specified to the nearest meter, and horizontal positions are specified in geographic coordinates. GTOPO30 data are distributed as tiles, most of which are 50\u00b0 in latitude by 40\u00b0 in longitude.\r\n\r\nGTOPO30 tiles are available for download from\u00a0<a href=\"http:\/\/eros.usgs.gov\/#\/Find_Data\/Products_and_Data_Available\/gtopo30_info\">USGS\u2019 EROS Data Center<\/a>. GTOPO60, a resampled and untiled version of GTOPO30, is available through the\u00a0<a href=\"http:\/\/seamless.usgs.gov\/\">USGS\u2019 Seamless Data Distribution Service<\/a>.\r\n<h3><\/h3>\r\n<h3>SHUTTLE RADAR TOPOGRAPHY MISSION (SRTM)<\/h3>\r\nFrom February 11 to February 22, 2000, the space shuttle Endeavor bounced radar waves off the Earth\u2019s surface, and recorded the reflected signals with two receivers spaced 60 meters apart. The mission measured the elevation of land surfaces between 60\u00b0 N and 57\u00b0 S latitude.\u00a0<strong>The highest resolution data products created from the SRTM mission are 30 meters<\/strong>. Access to 30-meter SRTM data for areas outside the U.S. are restricted by the National Geospatial-Intelligence Agency, which sponsored the project along with the National Aeronautics and Space Administration (NASA). A 90-meter SRTM data product is available for free download without restriction (Maune, 2007).\r\n\r\n<img alt=\"Anaglyph stereo image of terrain surface of Fiji mapping by Shuttle Radar Topography Mission\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/shuttle_stereo_fiji.jpg\" \/>\r\n\r\nAnaglyph stereo image derived from Shuttle Radar Topography Mission data (NASA Jet Propulsion Laboratory, 2006).\r\n\r\nThe image above shows Viti Levu, the largest of the some 332 islands that comprise the Sovereign Democratic Republic of the Fiji Islands. Viti Levu\u2019s area is 10,429 square kilometers (about 4000 square miles). Nakauvadra, the rugged mountain range running from north to south, has several peaks rising above 900 meters (about 3000 feet). Mount Tomanivi, in the upper center, is the highest peak at 1324 meters (4341 feet).\r\n\r\nLearn more about the Shuttle Radar Topography Mission at Web sites published by\u00a0<a href=\"http:\/\/www.jpl.nasa.gov\/srtm\">NASA<\/a>\u00a0and\u00a0<a href=\"http:\/\/srtm.usgs.gov\/mission.php\">USGS<\/a>.\r\n<h2>7.14. Bathymetry<\/h2>\r\nThe term bathymetry refers to the process and products of measuring the depth of water bodies. The U.S. Congress authorized the comprehensive mapping of the nation\u2019s coasts in 1807, and directed that the task be carried out by the federal government\u2019s first science agency, the Office of Coast Survey (OCS). That agency is now responsible for mapping some 3.4 million nautical square miles encompassed by the 12-mile territorial sea boundary, as well as the 200-mile Exclusive Economic Zone claimed by the U.S., a responsibility that entails regular revision of about 1,000 nautical charts. The coastal bathymetry data that appears on USGS topographic maps, like the one shown below, is typically compiled from OCS charts.\r\n\r\n<img alt=\"Portion of topographic map showing ocean depths\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/bar_harbor_bathy.gif\" \/>\r\n\r\n\u201dIsobaths\u201d (the technical term for lines of constant depth) shown on a USGS topographic map.\r\n\r\nEarly hydrographic surveys involved sampling water depths by casting overboard ropes weighted with lead and marked with depth intervals called marks and deeps. Such ropes were called leadlines for the weights that caused them to sink to the bottom. Measurements were called soundings. By the late 19th century, piano wire had replaced rope, making it possible to take soundings of thousands rather than just hundreds of fathoms (a fathom is six feet).\r\n\r\n<img alt=\"\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/sounding.jpg\" \/>\r\n\r\nSeaman paying out a sounding line during a hydrographic survey of the East coast of the U.S. in 1916. (NOAA, 2007).\r\n\r\nEcho sounders were introduced for deepwater surveys beginning in the 1920s.\u00a0<strong>Sonar\u00a0<\/strong>(SOund NAvigation and Ranging) technologies have revolutionized oceanography in the same way that aerial photography revolutionized topographic mapping. The seafloor topography revealed by sonar and related shipborne remote sensing techniques provided evidence that supported theories about seafloor spreading and plate tectonics.\r\n\r\nBelow is an artist\u2019s conception of an oceanographic survey vessel operating two types of sonar instruments:\u00a0<strong>multibeam<\/strong>\u00a0and s<strong>ide scan sonar<\/strong>. On the left, a multibeam instrument mounted in the ship\u2019s hull calculates ocean depths by measuring the time elapsed between the sound bursts it emits and the return of echoes from the seafloor. On the right, side scan sonar instruments are mounted on both sides of a submerged \u201ctowfish\u201d tethered to the ship. Unlike multibeam, side scan sonar measures the strength of echoes, not their timing. Instead of depth data, therefore, side scanning produces images that resemble black-and-white photographs of the sea floor.\r\n\r\n<img alt=\"Illustration of sonar in use for bathymetric mapping\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/sonar.jpg\" \/>\r\n\r\nMultibeam and side scan sonar in use for bathymetric mapping. (NOAA, 2002).\r\n\r\nA\u00a0<a href=\"http:\/\/craterlake.wr.usgs.gov\/bathymetry.html\">detailed report of the recent bathymetric survey of Crater Lake, Oregon, USA, is published by the USGS here<\/a>.\r\n<h2>7.15. Statistical Surfaces<\/h2>\r\nStrategies used to represent terrain surfaces can be used for other kinds of surfaces as well. For example, one of my first projects here at Penn State was to work with a distinguished geographer, the late Peter Gould, who was studying the diffusion of the Acquired Immune Deficiency Syndrome (AIDS) virus in the United States. Dr. Gould had recently published the map below.\r\n\r\n<img alt=\"Thematic map depicted HIV\/AIDS as a statistical surface\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/aids_surface.jpg\" \/>\r\n\r\n&nbsp;\r\n\r\nOblique view of contour lines representing distribution of AIDS cases in the U.S. 1988. (Gould, 1989. \u00a9 Association of American Geographers. All rights reserved. Reproduced here for educational purposes only).\r\n\r\n&nbsp;\r\n\r\nGould portrayed the distribution of disease in the same manner as another geographer might portray a terrain surface. The portrayal is faithful to Gould\u2019s conception of the contagion as a continuous phenomenon. It was important to Gould that people understood that there was no location that did not have the potential to be visited by the epidemic. For both the AIDS surface and a terrain surface, a quantitative attribute (<em>z<\/em>) exists for every location (<em>x,y<\/em>). In general, when a\u00a0<strong>continuous phenomenon is conceived as being analogous to the terrain surface, the conception is called a statistical surface<\/strong>.\r\n<h2>7.16. Theme: Hydrography<\/h2>\r\nThe NSDI Framework Introduction and Reference (FGDC, 1997) envisions the hydrography theme in this way:\r\n\r\nFramework hydrography data include surface water features such as lakes and ponds, streams and rivers, canals, oceans, and shorelines. Each of these features has the attributes of a name and feature identification code. Centerlines and polygons encode the positions of these features. For feature identification codes, many federal and state agencies use the Reach schedule developed by the U.S. Environmental Protection Agency (EPA).\r\n\r\n&nbsp;\r\n\r\nMany hydrography data users need complete information about connectivity of the hydrography network and the direction in which the water flows encoded in the data. To meet these needs, additional elements representing flows of water and connections between features may be included in framework data (p. 20).\r\n<h3>IDENTIFICATION<\/h3>\r\nFGDC had the\u00a0<strong>National Hydrography Dataset (NHD)<\/strong>\u00a0in mind when they wrote this description. NHD combines the vector features of<strong>Digital Line Graph (DLG) hydrography<\/strong>\u00a0with the\u00a0<strong>EPA\u2019s Reach files<\/strong>. Reaches are segments of surface water that share similar hydrologic characteristics. Reaches are of three types: transport, coastline, and waterbody. DLG lines features represent the transport and coastline types; polygon features are used to represent waterbodies. Every reach segment in the NHD is assigned a unique reach code, along with a host of other hydrological attributes including stream flow direction (which is encoded in the digitizing order of nodes that make up each segment), network connectivity, and feature names, among others. Because the order of reach codes are sequential from reach to reach, point-source data (such as a pollutant spill) can be geocoded to the affected reach. Used in this way, reaches comprise a<strong>linear referencing system<\/strong>\u00a0comparable to postal addresses along streets\u00a0 (USGS, 2002).\r\n\r\n<img alt=\"Diagram showing how water flow is attributed to reaches\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/nhd_flows.png\" \/>\r\n\r\nHow flow attributes are associated with reaches in the National Hydrography Dataset (USGS, 2000).\r\n\r\nNHD parses the U.S. surface drainage network into four hierarchical categories of units: 21 Regions, 222 Subregions, 352 Accounting units, and 2150 Cataloging units (also called Watersheds). Features can exist at multiple levels of the hierarchy, though they might not be represented in the same way. For example, while it might make the most sense to represent a given stream as a polygon features at the Watershed level, it may be more aptly represented as a line feature at the Region or Subregion level. NHD supports this by allowing multiple features to share the same reach codes. Another distinctive feature of NHD is\u00a0<strong>artificial flowlines<\/strong>\u2013centerline features that represent paths of water flow through polygon features such as standing water bodies. NHD is complex because it is designed to support sophisticated hydrologic modeling tasks, including point-source pollution modeling, flood potential, bridge construction, among others (Ralston, 2004).\r\n\r\n<img alt=\"Diagram illustrating how hydrographic features are represented with points, lines, and polygons\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/nhd_points_lines_areas.png\" \/>\r\n\r\nHow vector features are used to represent various types of reaches in the National Hydrography Dataset (USGS, 2000).\r\n\r\nNHD are available at three levels of detail (scale): medium (1:100,000, which is available for the entire U.S.), high (1:24,000, production of which is underway, \u201caccording to the availability of matching resources from NHD partners\u201d (USGS, 2002, p. 2), and local (larger scales such as 1:5,000), which \u201cis being developed where partners and data exist\u201d for select areas (USGS, 2006c; USGS, 2009; USGS 2013).\r\n<h3>SPATIAL REFERENCE INFORMATION<\/h3>\r\nNHD coordinates are decimal degrees referenced to the NAD 83 horizontal datum.\r\n<h3>DISTRIBUTION<\/h3>\r\n<h3><strong>TRY THIS!<\/strong><\/h3>\r\n<h3>DOWNLOAD AND VIEW AN EXTRACT FROM THE NATIONAL HYDROGRAPHY DATASET<\/h3>\r\n<ol>\r\n\t<li>From the\u00a0<a href=\"http:\/\/nhd.usgs.gov\/\">NHD home page<\/a>\u00a0click the\u00a0<strong>Get Data<\/strong>\u00a0link and then follow the link to the\u00a0<strong>NHD Viewer<\/strong>. (There is a Help button next to the viewer link.)<\/li>\r\n\t<li>Use the GIS tools, found above the map area, to pan to and zoom in on an area of interest. Then, click the\u00a0<strong>Download Data<\/strong>\u00a0button.\r\nChoose a\u00a0<em>reference area<\/em>\u00a0from the pick list. (You could also chose the entire\u00a0<em>current map extent<\/em>, but depending upon your zoom level that may be a huge dataset.)\r\nThen click on the map to highlight a specific area of interest. A link to the available data sets will appear in the left hand pane under the<strong>Selection\u00a0<\/strong>tab. (The\u00a0<strong>All Results<\/strong>\u00a0button will list the multiple areas you click on.) Follow the\u00a0<strong>Download\u00a0<\/strong>link of the area you are interested in.<\/li>\r\n\t<li>In the\u00a0<strong>USGS Available Data<\/strong>\u00a0window that opens, select<strong>Hydrography\u00a0<\/strong>from the\u00a0<strong>Theme<\/strong>\u00a0column, and choose a file format from the pick list in the\u00a0<strong>Format\u00a0<\/strong>column. I chose\u00a0<strong>Shapefile<\/strong>\u00a0format for my extract, because I know it is compatible with Global Mapper \/ dlgv32 Pro. If you were working in ArcGIS you could choose the File Geodatabase option. Ralston (2004, p.187) observes that NHD \u201cis precisely the type of information that could benefit from an integrated data model in an object relational database.\u201d\r\nClick the\u00a0<strong>Next\u00a0<\/strong>button.<\/li>\r\n\t<li>From the list of Hydrography Products check the box for what you wish to download. To be assured of getting the data in Shapefile format select a product referred to as a\u00a0<strong>Dynamic Extract<\/strong>.\r\nClick the\u00a0<strong>Next\u00a0<\/strong>button.\r\nGo to the\u00a0<strong>Cart\u00a0<\/strong>pane on the right. (It may open automatically.)\r\nGo to\u00a0<strong>Checkout<\/strong>, supply your e-mail address, and submit your order via the\u00a0<strong>Place Order<\/strong>\u00a0button.\r\nYou will receive a message telling you that your order has been placed, and that will soon be followed by an e-mail regarding your order. About an hour and a half after I submitted my request I received a second e-mail containing a download link.<\/li>\r\n\t<li>Extract the contents of the .zip file and view the Shapefile data set(s) in Global Mapper.<\/li>\r\n\t<li>Use the Identify pointer tool to reveal attributes of the reaches. In the example below I have highlighted a flowline associated with Cedar Creek in western Michigan.<\/li>\r\n<\/ol>\r\n<img alt=\"Screenshot of the feature information window in Global Mapper\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/NHD_trythis_sample_Aug2011.jpg\" \/>\r\n\r\n&nbsp;\r\n<h2>7.17. Theme: Transportation<\/h2>\r\nTransportation network data are valuable for all sorts of uses, including two we considered in Chapter 4: geocoding and routing. The Federal Geographic Data Committee (1997, p. 19) specified the following vector features and attributes for the transportation framework theme:\r\n<table summary=\"Description of each transportation network feature\"><caption>Transportation Framework Attributes<\/caption>\r\n<thead>\r\n<tr>\r\n<th><strong>Feature<\/strong><\/th>\r\n<th><strong>Attributes<\/strong><\/th>\r\n<\/tr>\r\n<\/thead>\r\n<tbody>\r\n<tr>\r\n<td>Roads<\/td>\r\n<td>Centerlines, feature identification code (using linear referencing systems where available), functional class, name (including route numbers), and street address ranges<\/td>\r\n<\/tr>\r\n<tr>\r\n<td>Trails<\/td>\r\n<td>Centerlines, feature identification code (using linear referencing systems where available), name, and type<\/td>\r\n<\/tr>\r\n<tr>\r\n<td>Railroads<\/td>\r\n<td>Centerlines, feature identification code (using linear referencing systems where available), and type<\/td>\r\n<\/tr>\r\n<tr>\r\n<td>Waterways<\/td>\r\n<td>Centerlines, feature identification code (using linear referencing systems where available), and name<\/td>\r\n<\/tr>\r\n<tr>\r\n<td>Airports and ports<\/td>\r\n<td>Feature identification code and name<\/td>\r\n<\/tr>\r\n<tr>\r\n<td>Bridges and tunnels<\/td>\r\n<td>Feature identification code and name<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n<h3>IDENTIFICATION<\/h3>\r\nAs part of the National Map initiative, USGS and partners are developing a comprehensive national database of vector transportation data. The transportation theme \u201cincludes best available data from Federal partners such as the Census Bureau and the Department of Transportation, State and local agencies\u201d (USGS, 2007).\r\n\r\nAs envisioned by FGDC,\u00a0<strong>centerlines<\/strong>\u00a0are used to represent transportation routes. Like the lines painted down the middle of two-way streets, centerlines are 1-dimensional vector features that approximate the locations of roads, railroads, and navigable waterways. In this sense, road centerlines are analogous to the flowpaths encoded in the National Hydrologic Dataset (see previous page). Also like the NHD (and TIGER), road\u00a0<strong>topology<\/strong>\u00a0must be encoded to facilitate analysis of transportation networks.\r\n\r\nTo get a sense of the complexity of the features and attributes that comprise the transportation theme, see the\u00a0<a href=\"http:\/\/services.nationalmap.gov\/bestpractices\/model\/acrodocs\/Poster_BPTrans_03_01_2006.pdf\">Transportation Data Model<\/a>(This is a 36\u2033 x 48\u2033 poster in a 5.2 Mb PDF file.) [The link to the Transportation Data Model poster recently became disconnected. Instead look at the model diagrams in the\u00a0<a title=\"FGDC Transportation Base\" href=\"http:\/\/www.fgdc.gov\/standards\/projects\/FGDC-standards-projects\/framework-data-standard\/GI_FrameworkDataStandard_Part7_Transportation_Base.pdf\">Part 7: Transportation Base<\/a>of the FGDC Geographic Framework Data Content Standard.]\r\n\r\nIn the U.S. at least, the best road centerline data is that produced by NAVTEQ and Tele Atlas, which license transportation data to routing sites like Google Maps and MapQuest, and to manufacturers of in-car GPS navigation systems. Because these data are proprietary, however, USGS must look elsewhere for data that can be made available for public use. TIGER\/Line data produced by the Census Bureau will likely play an important role after the TIGER\/MAF Modernization project is complete (see Chapter 4).\r\n<h3>DISTRIBUTION<\/h3>\r\n<h3><strong>TRY THIS!<\/strong><\/h3>\r\n<h3>VIEW AND DOWNLOAD NATIONAL MAP TRANSPORTATION DATA<\/h3>\r\n<ol>\r\n\t<li><a href=\"http:\/\/bpgeo.cr.usgs.gov\/\">Access the Viewer here<\/a>.<\/li>\r\n\t<li>Expand the pane containing the layer options by clicking on Overlays at the upper-left.<\/li>\r\n\t<li>Under Base Data Layers, click on Transportation. You can expand the Transportation list and sub-select different layers.<\/li>\r\n\t<li>As you zoom in to larger map scales (using the slider bar at the upper-left of the map), additional transportation layers will become visible.<\/li>\r\n\t<li>If you wish to download an extract from the transportation database, click the Download Data button in the upper-right of the viewer interface and decide how you wish to extract the data. The Transportation data comes down in ESRI\u2019s geodatabase format. Additional information regarding downloading data can be found via the Help button in the upper-right of the viewer interface.<\/li>\r\n<\/ol>\r\n<h2>7.18. Theme: Governmental Units<\/h2>\r\nThe FGDC framework also includes boundaries of governmental units, including:\r\n<ul>\r\n\t<li>Nation<\/li>\r\n\t<li>States and statistically equivalent areas<\/li>\r\n\t<li>Counties and statistically equivalent areas<\/li>\r\n\t<li>Incorporated places and consolidated cities<\/li>\r\n\t<li>Functioning legal minor civil divisions<\/li>\r\n\t<li>Federal- or state-recognized American Indian reservations and trustlands<\/li>\r\n\t<li>Alaska native regional corporations<\/li>\r\n<\/ul>\r\nFGDC specifies that:\r\n\r\nEach of these features includes the attributes of name and the applicable Federal Information Processing Standard (FIPS) code. Features boundaries include information about other features (such as road, railroads, or streams) with which the boundaries are associated and a description of the association (such as coincidence, offset, or corridor. (FGDC, 1997, p. 20-21)\r\n<h3>IDENTIFICATION<\/h3>\r\nThe USGS National Map aspires to include a comprehensive database of boundary data. In addition to the entities outlined above, the National Map also lists congressional districts, school districts, and ZIP Code zones. Sources for these data include \u201cFederal partners such as the U.S. Census Bureau, other Federal agencies, and State and local agencies.\u201d (USGS, 2007).\r\n\r\nTo get a sense of the complexity of the features and attributes that comprise this theme, see the\u00a0<a href=\"http:\/\/services.nationalmap.gov\/bestpractices\/model\/acrodocs\/Poster_BPGovtUnits_03_01_2006.pdf\">Governmental Units Data Model<\/a>\u00a0(This is a 36\u2033 x 48\u2033 poster in a 2.4 Mb PDF file.) [The link to the Governmental Units Data Model poster recently became disconnected. Instead look at the model diagrams in the\u00a0<a title=\"FGDC Governmental Unit Base\" href=\"http:\/\/www.fgdc.gov\/standards\/projects\/FGDC-standards-projects\/framework-data-standard\/GI_FrameworkDataStandard_Part5_GovernmentalUnitBoundaries.pdf\">Part 5: Governemntal unit and other geographic area boundaries<\/a>\u00a0of the FGDC Geographic Framework Data Content Standard.]\r\n<h3>DISTRIBUTION<\/h3>\r\n<h3><strong>TRY THIS!<\/strong><\/h3>\r\n<h3>VIEW AND DOWNLOAD NATIONAL MAP GOVERNMENTAL UNITS DATA<\/h3>\r\n<ol>\r\n\t<li><a href=\"http:\/\/bpgeo.cr.usgs.gov\/\">Access the Viewer here<\/a>.<\/li>\r\n\t<li>Expand the pane containing the layer options by clicking on Overlays at the upper-left.<\/li>\r\n\t<li>Under Base Data Layers, click on Governmental Unit Boundaries. You can expand the this list and sub-select different boundary layers.<\/li>\r\n\t<li>As you zoom in to larger map scales (using the slider bar at the upper-left of the map), additional boundary layers will become visible.<\/li>\r\n\t<li>If you wish to download an extract from the Governmental Unit Boundaries database, click the Download Data button in the upper-right of the viewer interface and decide how you wish to extract the data. The Governmental Unit Boundaries data comes down in ESRI\u2019s geodatabase format. Additional information regarding downloading data can be found via the Help button in the upper-right of the viewer interface.<\/li>\r\n<\/ol>\r\n<h2>7.19. Theme: Cadastral<\/h2>\r\nFGDC (1997, p. 21) points out that:\r\n\r\nCadastral data represent the geographic extent of the past, current, and future rights and interests in real property. The spatial information necessary to describe the geographic extent and the rights and interests includes surveys, legal description reference systems, and parcel-by-parcel surveys and descriptions.\r\n\r\nHowever, no one expects that legal descriptions and survey coordinates of private property boundaries (as depicted schematically in the portion of the\u00a0<strong>plat map<\/strong>\u00a0shown below) will be included in the USGS National Map any time soon. As discussed at the outset of Chapter 6, this is because local governments have authority for land title registration in the U.S., and most of these governments have neither the incentive nor the means to incorporate such data into a publicly-accessible national database.\r\n\r\n<img alt=\"Portion of a plat map showing property boundaries\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/plat.png\" \/>\r\n\r\nPlat maps are supplementary records that depict property parcel boundaries in graphic form. The geometric accuracy of plats is notoriously poor. The investment required to convert plat maps to properly georeferenced digital data is substantial. Many local governments have converted these records to digital form, or are in the process of doing so.\r\n\r\nFGDC\u2019s modest goal for the cadastal theme of the NSDI framework is to include:\r\n\r\n\u2026cadastral reference systems, such as the Public Land Survey System (PLSS) and similar systems not covered by the PLSS \u2026 and publicly administered parcels, such as military reservations, national forests, and state parks. (Ibid, p. 21)\r\n\r\n<a href=\"http:\/\/www.fgdc.gov\/standards\/standards_publications\/\">FGDC\u2019s Cadastral Data Content Standard is published here<\/a>.\r\n\r\nThe colored areas on the map below show the extent of the United States Public Land Surveys, which commenced in 1784 and took nearly a century to complete (Muehrcke and Muehrcke, 1998). The purpose of the surveys was to partition \u201cpublic land\u201d into saleable parcels in order to raise revenues needed to retire war debt, and to promote settlement. A key feature of the system is its nomenclature, which provides concise, unique specifications of the location and extent of any parcel.\r\n\r\n<img alt=\"Map of United States Public Land Survey system\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/USPLS.jpg\" \/>\r\n\r\nExtent of the U.S. Public Land Survey (Thompson, 1988).\r\n\r\nEach Public Land Survey (shown in the colored areas above) commenced from an\u00a0<strong>initial point\u00a0<\/strong>at the precisely surveyed intersection of a\u00a0<strong>base line<\/strong>\u00a0and\u00a0<strong>principal meridian<\/strong>. Surveyed lands were then partitioned into grids of\u00a0<strong>townships\u00a0<\/strong>each approximately six miles square.\r\n\r\n<img alt=\"U.S. Public Land Survey Township grid system\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/uspls_township_range.png\" \/>\r\n\r\n<strong>Townships\u00a0<\/strong>are designated by their locations relative to the base line and principal meridian of a particular survey. For example, the township highlighted in gold above is the second township south of the baseline and the third township west of the principal meridian. The Public Land Survey designation for the highlighted township is \u201cTownship 2 South, Range 3 West.\u201d Because of this nomenclature, the Public Land Survey System is also known as the \u201ctownship and range system.\u201d\u00a0 Township T2S, R3W is shown enlarged below.\r\n\r\n<img alt=\"U.S. Public Land Survey township grid\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/uspls_township.png\" \/>\r\n\r\nTownships are subdivided into grids of 36\u00a0<strong>sections<\/strong>. Each section covers approximately\u00a0<strong>one square mile<\/strong>\u00a0(640 acres). Notice the back-and-forth numbering scheme. Section 14, highlighted in gold above, is shown enlarged below.\r\n\r\n<img alt=\"U.S. Public Land Survey section showing property designations\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/uspls_section.png\" \/>\r\n\r\nInidividual property parcels are designated as shown above. For instance, the\u00a0<strong>NE 1\/4 of Section 14, Township 2 S, Range 3W<\/strong>, is a 160-acre parcel. Public Land Survey designations specify both the location of a parcel and its area.\r\n\r\n<img alt=\"Portion of topographic map showning influence of Public Land Survey on road network in midwest U.S.\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/uspls_topo.png\" \/>\r\n\r\nThe influence of the Public Land Survey grid is evident in the built environment of much of the American Midwest. As Mark Monmonier (1995, p. 114) observes:\r\n\r\nThe result [of the U.S. Public Land Survey] was an \u2018authored landscape\u2019 in which the survey grid had a marked effect on settlement patterns and the shapes of counties and smaller political units. In the typical Midwestern county, roads commonly following section lines, the rural population is dispersed rather than clustered, and the landscape has a pronounced checkerboard appearance.\r\n\r\nFor more information about the\u00a0<a href=\"http:\/\/nationalatlas.gov\/articles\/boundaries\/a_plss.html\">Public Land Survey System, see this article in the in the USGS\u2019 National Atlas<\/a>.\r\n<h2>7.20. Summary<\/h2>\r\nNSDI framework data represent \u201cthe most common data themes [that] users need\u201d (FGDC, 1997, p. 3), including geodetic control, orthoimagery, elevation, hydrography, transportation, governmental unit boundaries, and cadastral reference information. Some themes, like transportation and governmental units, represent things that have well-defined edges. In this sense we can think of things like roads and political boundaries as\u00a0<strong>discrete\u00a0<\/strong>phenomena. The\u00a0<strong>vector\u00a0<\/strong>approach to geographic representation is well suited to digitizing discrete phenomena. Line features do a good job of representing roads, for example, and polygons are useful approximations of boundaries.\r\n\r\nAs you recall from Chapter 1, however, one of the distinguishing properties of the Earth\u2019s surface is that it is\u00a0<strong>continuous<\/strong>. Some phenomena distributed across the surface are continuous too. Terrain elevations, gravity, magnetic declination and surface air temperature can be measured practically everywhere. For many purposes,\u00a0<strong>raster<\/strong>data are best suited to representing continuous phenomena.\r\n\r\nAn implication of continuity is that there is an infinite number of locations at which phenomena can be measured. It is not possible, obviously, to take an infinite number of measurements. Even if it were, the mass of data produced would not be usable. The solution, of course, is to collect a sample of measurements, and to estimate attribute values for locations that are left unmeasured. Chapter 7 also considers how missing elevations in a raster grid can be estimated from existing elevations, using a procedure called\u00a0<strong>interpolation<\/strong>. The inverse distance weighted interpolation procedure relies upon another fundamental property of geographic data,\u00a0<strong>spatial dependence<\/strong>.\r\n\r\nThe chapter concludes by investigating the characteristics and current status of the hydrography, transportation, governmental units, and cadastral themes. You had the opportunity to access, download, and open several of the data themes using viewers provided by USGS as part of its National Map initiative. In general, you should have found that although neither the NSDI or National Map visions have been fully realized, substantial elements of each is in place. Further progress depends on the American public\u2019s continuing commitment to public data, and to the political will of our representatives in government.\r\n<h3><strong>QUIZ<\/strong><\/h3>\r\nRegistered Penn State students should return now to the Chapter 7 folder\u00a0 in ANGEL (via the Resources menu to the left) to access the graded quiz for this chapter. This one counts.\u00a0<strong>You may take graded quizzes only once.<\/strong>\r\n\r\nThe purpose of the quiz is to ensure that you have studied the text closely, that you have mastered the practice activities, and that you have fulfilled the chapter\u2019s learning objectives. You are welcome to review the chapter during the quiz.\r\n\r\nOnce you have submitted the quiz and posted any questions you may have to either our discussion forums or chapter pages, you will have completed Chapter 7.\r\n<h3>COMMENTS AND QUESTIONS<\/h3>\r\nRegistered students are welcome to post comments, questions, and replies to questions about the text. Particularly welcome are anecdotes that relate the chapter text to your personal or professional experience. In addition, there are discussion forums available in the ANGEL course management system for comments and questions about topics that you may not wish to share with the whole world.\r\n\r\nTo post a comment, scroll down to the text box under \u201cPost new comment\u201d and begin typing in the text box, or you can choose to reply to an existing thread. When you are finished typing, click on either the \u201cPreview\u201d or \u201cSave\u201d button (Save will actually submit your comment). Once your comment is posted, you will be able to edit or delete it as needed. In addition, you will be able to reply to other posts at any time.\r\n\r\nNote: the first few words of each comment become its \u201ctitle\u201d in the thread.\r\n<h2>7.21. Bibliography<\/h2>\r\nFederal Geographic Data Committee (1997).\u00a0<em>Framework introduction and guide.\u00a0<\/em>Washington DC: Federal Geographic Data Committee.\r\n\r\nEischeid, J. D., Baker, C. B., Karl, R. R., Diaz, H. F. (1995). The quality control of long-term climatological data using objective data analysis.\u00a0<em>Journal of Applied Meteorology, 34<\/em>, 27-88.\r\n\r\nGould, P. (1989). Geographic dimensions of the AIDS epidemic.<em>Professional Geographer, 41:1<\/em>, 71-77.\r\n\r\nMaune, D. F. (Ed.) (2007).\u00a0<em>Digital elevation model technologies and applications: The DEM users manual, 2nd edition.\u00a0<\/em>Bethesda, MD: American Society for Photogrammetric Engineering and Remote Sensing.\r\n\r\nMonmonier, M. S. (1982).\u00a0<em>Drawing the line: tales of maps and cartocontroversy.<\/em>\u00a0New York, NY: Henry Holt.\r\n\r\nMuehrcke, P. C. and Muehrcke, J. O. (1998)\u00a0<em>Map use,<\/em>\u00a04th Ed. Madison, WI: JP Publications.\r\n\r\nNational Aeronautics and Space Administration, Jet Propulsion Laboratory (2006).\u00a0<em>Shuttle radar topography mission<\/em>. Retrieved May 10, 2006, from\u00a0<a href=\"http:\/\/www.jpl.nasa.gov\/srtm\">http:\/\/www.jpl.nasa.gov\/srtm<\/a>\r\n\r\nGoddard Space Flight Center, National Aeronautics and Space Administration (n.d.). Greenland\u2019s receding ice. Retrieved Feburary 26, 2008, from\u00a0<a href=\"http:\/\/svs.gsfc.nasa.gov\/stories\/greenland\/\">http:\/\/svs.gsfc.nasa.gov\/stories\/greenland\/<\/a>\r\n\r\nNational Geophysical Data Center (2010). ETOPO1 global gridded 1 arc-minute database. Retrieved March 2, 2010, from<a href=\"http:\/\/www.ngdc.noaa.gov\/mgg\/global\/global.html\">http:\/\/www.ngdc.noaa.gov\/mgg\/global\/global.html<\/a>\r\n\r\nNational Oceanic and Atmospheric Administration, National Climatic Data Center (n. d.).\u00a0<em>Merged land-ocean seasonal temperature anomalies<\/em>. Retrieved August 18, 1999, from<a title=\"http:\/\/www.ncdc.noaa.giv\/onlineprod\/landocean\/seasonal\/form.html\" href=\"http:\/\/www.ncdc.noaa.giv\/onlineprod\/landocean\/seasonal\/form.html\">http:\/\/www.ncdc.noaa.giv\/onlineprod\/landocean\/seasonal\/form.html<\/a>(expired)\r\n\r\nNational Oceanic and Atmospheric Administration (2002). Side scan and multibeam sonar. Retrieved February 18, 2008, from<a title=\"http:\/\/www.nauticalcharts.noaa.gov\/hsd\/hydrog.htm\" href=\"http:\/\/www.nauticalcharts.noaa.gov\/hsd\/hydrog.htm\">http:\/\/www.nauticalcharts.noaa.gov\/hsd\/hydrog.htm<\/a>\r\n\r\nNational Oceanic and Atmospheric Administration (2007). NOAA History. Retrieved February 27, 2008, from<a href=\"http:\/\/www.history.noaa.gov\/\">http:\/\/www.history.noaa.gov\/<\/a>\r\n\r\nRabenhorst, T. D. and McDermott, P. D. (1989).\u00a0<em>Applied cartography: source materials for mapmaking.\u00a0<\/em>Columbus, OH: Merrill.\r\n\r\nRaitz, E. (1948).\u00a0<em>General cartography.<\/em>\u00a0New York, NY: McGraw-Hill.\r\n\r\nRalston, B. A. (2004).\u00a0<em>GIS and public data.<\/em>\u00a0Clifton Park NY: Delmar Learning.\r\n\r\nThompson, M. M. (1988) Maps for america, 3rd Ed. Reston, VA: United States Geological Survey.\r\n\r\nUnited States Geological Survey (1987) Digital elevation models. Data users guide 5. Reston, VA: USGS.\r\n\r\nUnited States Geological Survey (1999) The National Hydrography Dataset. Fact Sheet 106-99. Reston, VA: USGS. Retrieved February 19, 2008 from\u00a0<a href=\"http:\/\/erg.usgs.gov\/isb\/pubs\/factsheets\/fs10699.html\">http:\/\/erg.usgs.gov\/isb\/pubs\/factsheets\/fs10699.html<\/a>\r\n\r\nUnited States Geological Survey (2000) The National Hydrography Dataset: Concepts and Contents. Reston, VA: USGS. Retrieved February 19, 2008 from<a href=\"http:\/\/nhd.usgs.gov\/chapter1\/chp1_data_users_guide.pdf\">http:\/\/nhd.usgs.gov\/chapter1\/chp1_data_users_guide.pdf<\/a>\r\n\r\nUnited States Geological Survey (2002) The National Map \u2013 Hydrography. Fact Sheet 060-02. Reston, VA: USGS. Retrieved February 19, 2008 from<a href=\"http:\/\/erg.usgs.gov\/isb\/pubs\/factsheets\/fs06002.html\">http:\/\/erg.usgs.gov\/isb\/pubs\/factsheets\/fs06002.html<\/a>\u00a0Retrieved September 22, 2013 from\u00a0<a href=\"http:\/\/pubs.er.usgs.gov\/publication\/fs06002\">http:\/\/pubs.er.usgs.gov\/publication\/fs06002<\/a>\r\n\r\nUnited States Geological Survey (2006a) Digital Line Graphs (DLG). Reston, VA: USGS. Retrieved February 18, 2008 from<a title=\"http:\/\/edc.usgs.gov\/products\/map\/dlg.html\" href=\"http:\/\/edc.usgs.gov\/products\/map\/dlg.html\">http:\/\/edc.usgs.gov\/products\/map\/dlg.html<\/a>\u00a0(In 2010 the site became<a href=\"http:\/\/eros.usgs.gov\/#\/Find_Data\/Products_and_Data_Available\/DLGs\">http:\/\/eros.usgs.gov\/#\/Find_Data\/Products_and_Data_Available\/DLGs<\/a>)\r\n\r\nUnited States Geological Survey (2006b) GTOPO30. Retrieved February 27, 2008 from<a title=\"http:\/\/edc.usgs.gov\/products\/elevation\/gtopo30\/gtopo30.html\" href=\"http:\/\/edc.usgs.gov\/products\/elevation\/gtopo30\/gtopo30.html\">http:\/\/edc.usgs.gov\/products\/elevation\/gtopo30\/gtopo30.html<\/a>\u00a0since moved to\u00a0<a href=\"http:\/\/www1.gsi.go.jp\/geowww\/globalmap-gsi\/gtopo30\/gtopo30.html\">http:\/\/www1.gsi.go.jp\/geowww\/globalmap-gsi\/gtopo30\/gtopo30.html<\/a>\r\n\r\nUnited States Geological Survey (2006c) National Hydrography Dataset (NHD) \u2013 High-resolution (Metadata). Reston, VA: USGS. Retrieved February 19, 2008 from<a href=\"http:\/\/nhdgeo.usgs.gov\/metadata\/nhd_high.htm\">http:\/\/nhdgeo.usgs.gov\/metadata\/nhd_high.htm<\/a>\r\n\r\nUnited States Geological Survey (2007). Vector data theme development of<em>\u00a0The National Map.\u00a0<\/em>Retrieved 24 February 2008 from<a title=\"http:\/\/bpgeo.cr.usgs.gov\/model\/\" href=\"http:\/\/bpgeo.cr.usgs.gov\/model\/\">http:\/\/bpgeo.cr.usgs.gov\/model\/<\/a>\u00a0(expired or moved)\r\n\r\nUnited States Geological Survey (2009) The National Map \u2013 Hydrography Dataset. Reston, VA: USGS. Retrieved September 22, 2013 from\u00a0<a href=\"http:\/\/pubs.usgs.gov\/fs\/2009\/3054\/pdf\/FS2009-3054.pdf\">http:\/\/pubs.usgs.gov\/fs\/2009\/3054\/pdf\/FS2009-3054.pdf<\/a>\r\n\r\nUnited States Geological Survey (2013) National Hydrography Dataset (NHD) \u2013 Get NDH Data. Reston, VA: USGS. Retrieved September 22, 2013 from\u00a0<a href=\"http:\/\/nhd.usgs.gov\/data.html\">http:\/\/nhd.usgs.gov\/data.html<\/a>","rendered":"<h2>7.1. Overview<\/h2>\n<p>Chapters 6 and 7 consider the origins and characteristics of the <strong>framework data themes <\/strong>that make up the United States&#8217; proposed National Spatial Data Infrastructure (NSDI). Chapter 6 discussed the geodetic control and orthoimagery themes. <strong>This chapter describes the origins, characteristics and current status of the elevation, transportation, hydrography, governmental units and cadastral themes.<\/strong><\/p>\n<h3>Objectives<\/h3>\n<p>Students who successfully complete Chapter 7 should be able to:<\/p>\n<ol>\n<li>Given a regular or irregular array of spot elevations, construct a triangulated irregular network, interpolate contour intervals and draw contour lines;<\/li>\n<li>Compare vector and raster representations of terrain elevation;<\/li>\n<li>Acquire and view digital elevation data from the National Elevation Dataset;<\/li>\n<li>Calculate an interpolated spot elevation based on neighboring elevations;<\/li>\n<li>Contrast the characteristics of three global elevation data products;<\/li>\n<li>Describe the characteristics and current status of the NSDI hydrography, transportation, and governmental units themes as implemented in USGS&#8217; National Map; and<\/li>\n<li>Interpret the size and relative location of a land parcel designated in terms of the U.S. Public Land Survey System.<\/li>\n<\/ol>\n<h3>Comments and Questions<\/h3>\n<p>Registered students are welcome to post comments, questions, and replies to questions about the text. Particularly welcome are anecdotes that relate the chapter text to your personal or professional experience. In addition, there are discussion forums available in the ANGEL course management system for comments and questions about topics that you may not wish to share with the whole world.<\/p>\n<p>To post a comment, scroll down to the text box under &#8220;Post new comment&#8221; and begin typing in the text box, or you can choose to reply to an existing thread. When you are finished typing, click on either the &#8220;Preview&#8221; or &#8220;Save&#8221; button (Save will actually submit your comment). Once your comment is posted, you will be able to edit or delete it as needed. In addition, you will be able to reply to other posts at any time.<\/p>\n<p>Note: the first few words of each comment become its &#8220;title&#8221; in the thread.<\/p>\n<h2>7.2. Checklist<\/h2>\n<p>The following checklist is for Penn State students who are registered for classes in which this text, and associated quizzes and projects in the ANGEL course management system, have been assigned. You may find it useful to print this page out first so that you can follow along with the directions.<\/p>\n<table summary=\"Tasks to be completed for the chapter\">\n<caption>Chapter 7 Checklist (for registered students only)<\/caption>\n<thead>\n<tr>\n<th>Step<\/th>\n<th>Activity<\/th>\n<th>Access\/Directions<\/th>\n<\/tr>\n<\/thead>\n<tbody>\n<tr>\n<th>1<\/th>\n<td><strong>Read<\/strong>\u00a0Chapter 7<\/td>\n<td>This is the second page of the Chapter. Click on the links at the bottom of the page to continue or to return to the previous page, or to go to the top of the chapter. You can also navigate the text via the links in the GEOG 482 menu on the left.<\/td>\n<\/tr>\n<tr>\n<th>2<\/th>\n<td>Submit\u00a0<strong>3 practice quizzes<\/strong>including:<\/p>\n<ul>\n<li>Contouring<\/li>\n<li>DLGs and DEMs<\/li>\n<li>Interpolation<\/li>\n<\/ul>\n<p>Practice quizzes are not graded and may be submitted more than once.<\/td>\n<td>Go to ANGEL &gt; [your course section] &gt; Lessons tab &gt; Chapter 7 folder &gt; [quiz]<\/td>\n<\/tr>\n<tr>\n<th>3<\/th>\n<td>Perform\u00a0<strong>\u201cTry this\u201d activities<\/strong>including:<\/p>\n<ul>\n<li>Draw a contour map<\/li>\n<li>Explore Digital Line Graph hypsography<\/li>\n<li>Explore a Digital Elevation Model<\/li>\n<li>Download and view an extract from the National Elevation Dataset<\/li>\n<\/ul>\n<p>\u201cTry this\u201d activities are not graded.<\/td>\n<td>Instructions are provided for each activity.<\/td>\n<\/tr>\n<tr>\n<th>4<\/th>\n<td>Submit the<strong>Chapter 7 Graded Quiz<\/strong><\/td>\n<td>ANGEL &gt; [your course section] &gt; Lessons tab &gt; Chapter 7 folder &gt; Chapter 7 Graded Quiz. See the Calendar tab in ANGEL for due dates.<\/td>\n<\/tr>\n<tr>\n<th>5<\/th>\n<td>\u00a0Read<strong>comments and questions<\/strong>posted by fellow students. Add comments and questions of your own, if any.<\/td>\n<td>\u00a0Comments and questions may be posted on any page of the text, or in a Chapter-specific discussion forum in ANGEL.<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>&nbsp;<\/p>\n<h2>7.3. Theme: Elevation<\/h2>\n<p>The NSDI\u00a0<em>Framework Introduction and Guide<\/em>\u00a0(FGDC, 1997, p. 19) points out that \u201celevation data are used in many different applications.\u201d Civilian applications include flood plain delineation, road planning and construction, drainage, runoff, and soil loss calculations, and cell tower placement, among many others. Elevation data are also used to depict the terrain surface by a variety of means, from contours to relief shading and three-dimensional perspective views.<\/p>\n<p>The NSDI Framework calls for an \u201celevation matrix\u201d for land surfaces. That is, the terrain is to be represented as a grid of elevation values. The spacing (or resolution) of the elevation grid may vary between areas of high and low relief (i.e., hilly and flat). Specifically, the Framework Introduction states that<\/p>\n<p>Elevation values will be collected at a post-spacing of 2 arc-seconds (approximately 47.4 meters at 40\u00b0 latitude) or finer. In areas of low relief, a spacing of 1\/2 arc-second (approximately 11.8 meters at 40\u00b0 latitude) or finer will be sought (FGDC, 1997, p. 18).<\/p>\n<p>The elevation theme also includes bathymetry\u2013depths below water surfaces\u2013for coastal zones and inland water bodies. Specifically,<\/p>\n<p>For depths, the framework consists of soundings and a gridded bottom model. Water depth is determined relative to a specific vertical reference surface, usually derived from tidal observations. In the future, this vertical reference may be based on a global model of the geoid or the ellipsoid, which is the reference for expressing height measurements in the Global Positioning System (Ibid).<\/p>\n<p>USGS has lead responsibility for the elevation theme. Elevation is also a key component of USGS\u2019 National Map. The next several pages consider how heights and depths are created, how they are represented in digital geographic data, and how they may be depicted cartographically.<\/p>\n<h2>7.4. Vector and Raster Approaches<\/h2>\n<p>&nbsp;<\/p>\n<p>The terms\u00a0<strong>raster<\/strong>\u00a0and\u00a0<strong>vector<\/strong>\u00a0were introduced back in Chapter 1 to denote two fundamentally different strategies for representing geographic phenomena. Both strategies involve simplifying the infinite complexity of the Earth\u2019s surface. As it relates to elevation data, the raster approach involves\u00a0<em>measuring elevation at a sample of locations<\/em>. The vector approach, on the other hand, involves\u00a0<em>measuring the locations of a sample of elevations<\/em>. I hope that this distinction will be clear to you by the end of this chapter.<\/p>\n<p><img decoding=\"async\" alt=\"Diagram comparing contours and elevation grid depicting the same surface\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/vector_raster_surfs.gif\" \/><\/p>\n<p>Vector and raster representations of the same terrain surface.<\/p>\n<p>The illustration above compares how elevation data are represented in vector and raster data. On the left are elevation\u00a0<strong>contours<\/strong>, a vector representation that is familiar with anyone who has used a USGS topographic map. The technical term for an elevation contour is<em>isarithm<\/em>, from the Greek words for \u201csame\u201d and \u201cnumber.\u201d The terms<em>isoline<\/em>,\u00a0<em>isogram<\/em>, and\u00a0<em>isopleth<\/em>\u00a0all mean more or less the same thing. (See any cartography text for the distinctions.)<\/p>\n<p>As you will see later in this chapter, when you explore Digital Line Graph hypsography data using Global Mapper or dlgv 32 Pro, elevations in vector data are encoded as attributes of line features. The distribution of elevation points across the quadrangle is therefore irregular. Raster elevation data, by contrast, consist of grids of points at which elevation is encoded at regular intervals. Raster elevation data are what\u2019s called for by the NSDI Framework and the USGS National Map. Digital contours can now be rendered easily from raster data. However, much of the raster elevation data used in the National Map was produced from digital vector contours and hydrography (streams and shorelines). For this reason we\u2019ll consider the vector approach to terrain representation first.<\/p>\n<h2>7.5. Contours<\/h2>\n<p><img decoding=\"async\" alt=\"Perspective view of a terrain surface showing contour lines as traces of parallel horizontal planes\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/raisz_model.gif\" \/><\/p>\n<p>&nbsp;<\/p>\n<p>Contour lines trace the elevation of the terrain surface at regularly-spaced intervals (Raisz, 1948. \u00a9 McGraw-Hill, Inc. Used by permission).<\/p>\n<p>Drawing contour lines is a way to represent a terrain surface with a sample of elevations. Instead of measuring and depicting elevation at every point, you measure only along lines at which a series of imaginary horizontal planes slice through the terrain surface. The more imaginary planes, the more contours, and the more detail is captured.<\/p>\n<p><img decoding=\"async\" alt=\"Plan view of contour lines used to depict a terrain surface\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/raisz_plan.gif\" \/><\/p>\n<p>Contour lines representing the same terrain as in the first figure, but in plan view. (Raisz, 1948. \u00a9 McGraw-Hill, Inc. Used by permission).<\/p>\n<p>Until photogrammetric methods came of age in the 1950s, topographers in the field sketched contours on the USGS 15-minute topographic quadrangle series. Since then, contours shown on most of the 7.5-minute quads were compiled from stereoscopic images of the terrain, as described in Chapter 6. Today computer programs draw contours automatically from the spot elevations that photogrammetrists compile stereoscopically.<\/p>\n<p>Although it is uncommon to draw terrain elevation contours by hand these days, it is still worthwhile to know how. In the next few pages you\u2019ll have a chance to practice the technique, which is analogous to the way computers do it.<\/p>\n<h2>7.6. Contouring By Hand<\/h2>\n<p>This page will walk you through a methodical approach to rendering contour lines from an array of spot elevations (Rabenhorst and McDermott, 1989). To get the most from this demonstration, I suggest that you print the\u00a0<a href=\"https:\/\/www.e-education.psu.edu\/files\/geog482\/image\/contouring_lesson.gif\">illustration in the attached image file<\/a>. Find a pencil (preferably one with an eraser!) and straightedge, and duplicate the steps illustrated below. A \u201cTry This!\u201d activity will follow this step-by-step introduction, providing you a chance to go solo.<\/p>\n<p>&nbsp;<\/p>\n<p><img decoding=\"async\" alt=\"Step 1 of contouring demonstration\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/contouring_2.gif\" \/><\/p>\n<p>Beginning a triangulated irregular network.<\/p>\n<p>Starting at the highest elevation, draw straight lines to the nearest neighboring spot elevations. Once you have connected to all of the points that neighbor the highest point, begin again at the second highest elevation. (You will have to make some subjective decisions as to which points are \u201cneighbors\u201d and which are not.) Taking care not to draw triangles across the stream, continue until the surface is completely triangulated.<\/p>\n<p><img decoding=\"async\" alt=\"Step 2 of contouring demonstration\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/contouring_3.gif\" \/><\/p>\n<p>Complete TIN. Note that the triangle sides must not cross hydrologic features (i.e., the stream) on a terrain surface.<\/p>\n<p>The result is a\u00a0<strong>triangulated irregular network (TIN)<\/strong>. A TIN is a vector representation of a continuous surface that consists entirely of triangular facets. The vertices of the triangles are spot elevations that may have been measured in the field by leveling, or in a photogrammetrist\u2019s workshop with a stereoplotter, or by other means. (Spot elevations produced photogrammetrically are called\u00a0<strong>mass points<\/strong>.) A useful characteristic of TINs is that each triangular facet has a single slope degree and direction. With a little imagination and practice, you can visualize the underlying surface from the TIN even without drawing contours.<\/p>\n<p>Wonder why I suggest that you not let triangle sides that make up the TIN cross the stream? Well, if you did, the stream would appear to run along the side of a hill, instead of down a valley as it should. In practice, spot elevations would always be measured at several points along the stream, and along ridges as well. Photogrammetrists refer to spot elevations collected along linear features as\u00a0<strong>breaklines<\/strong>\u00a0(Maune, 2007). I omitted breaklines from this example just to make a point.<\/p>\n<p>You may notice that there is more than one correct way to draw the TIN. As you will see, deciding which spot elevations are \u201cnear neighbors\u201d and which are not is subjective in some cases. Related to this element of subjectivity is the fact that the fidelity of a contour map depends in large part on the distribution of spot elevations on which it is based. In general, the density of spot elevations should be greater where terrain elevations vary greatly, and sparser where the terrain varies subtly. Similarly, the smaller the contour interval you intend to use, the more spot elevations you need.<\/p>\n<p>(There are algorithms for triangulating irregular arrays that produce unique solutions. One approach is called\u00a0<strong>Delaunay Triangulation<\/strong>which, in one of its constrained forms, is useful for representing terrain surfaces. The distinguishing geometric characteristic of a Delaunay triangulation is that a circle surrounding each triangle side does not contain any other vertex.)<\/p>\n<p><img decoding=\"async\" alt=\"Step 3 of contouring demonstration\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/contouring_4.gif\" \/><\/p>\n<p>&nbsp;<\/p>\n<p>Tick marks drawn where elevation contours cross the edges of each TIN facet.<\/p>\n<p>&nbsp;<\/p>\n<p>Now draw ticks to mark the points at which elevation contours intersect each triangle side. For instance, see the triangle side that connects the spot elevations 2360 and 2480 in the lower left corner of the illustration above? One tick mark is drawn on the triangle where a contour representing elevation 2400 intersects. Now find the two spot elevations, 2480 and 2750, in the same lower left corner. Note that three tick marks are placed where contours representing elevations 2500, 2600, and 2700 intersect.<\/p>\n<p>This step should remind you of the equal interval classification scheme you read about in Chapter 3. The right choice of contour interval depends on the goal of the mapping project. In general, contour intervals increase in proportion to the variability of the terrain surface. It should be noted that the assumption that elevations increase or decrease at a constant rate is not always correct, of course. We will consider that issue in more detail later.<\/p>\n<p><img decoding=\"async\" alt=\"Step 4 of contouring demonstration\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/contouring_5.gif\" \/><\/p>\n<p>Threading elevation contours through a TIN.<\/p>\n<p>Finally, draw your contour lines. Working downslope from the highest elevation, thread contours through ticks of equal value. Move to the next highest elevation when the surface seems ambiguous.<\/p>\n<p>Keep in mind the following characteristics of contour lines (Rabenhorst and McDermott, 1989):<\/p>\n<ul>\n<li>Contours should\u00a0<strong>always point upstream<\/strong>\u00a0in valleys<\/li>\n<li>Contours should\u00a0<strong>always point downridge<\/strong>\u00a0along ridges<\/li>\n<li>Adjacent contours should<strong>\u00a0always be sequential or equivalent<\/strong><\/li>\n<li>Contours should\u00a0<strong>never split into two<\/strong><\/li>\n<li>Contours should\u00a0<strong>never cross or loop<\/strong><\/li>\n<li>Contours should\u00a0<strong>never spiral<\/strong><\/li>\n<li>Contours should\u00a0<strong>never stop in the middle of a map<\/strong><\/li>\n<\/ul>\n<p>How does your finished map compare with the one I drew below?<\/p>\n<p><img decoding=\"async\" alt=\"Outcome of contouring demonstration\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/contouring_6.gif\" \/><\/p>\n<h3><strong>TRY THIS<\/strong><strong>!<\/strong><\/h3>\n<p><strong>Now try your hand at contouring on your own.<\/strong>\u00a0The purpose of this practice activity is to give you more experience in contouring terrain surfaces.<\/p>\n<ol>\n<li>First, view\u00a0<a href=\"http:\/\/opentextbc.ca\/natureofgeoinfo\/files\/geog482\/image\/contouring_practice-apr2012.gif\">an image of an irregular array of 16 spot elevations<\/a>.<\/li>\n<li>Print the image.<\/li>\n<li>Use the procedure outlined in this lesson to draw contour lines that represent the terrain surface that the spot elevations were sampled from. You may find this to be a moderately challenging task that takes about a half hour to do well. TIP: label the tick marks to make it easier to connect them.<\/li>\n<li>When finished,\u00a0<a href=\"https:\/\/www.e-education.psu.edu\/files\/geog482\/image\/mt_nittany.jpg\">compare your result to an existing map.<\/a><\/li>\n<\/ol>\n<p>Here are a couple of somewhat simpler problems and solutions in case you need a little more practice.<\/p>\n<ul>\n<li><a href=\"https:\/\/www.e-education.psu.edu\/files\/geog482\/image\/cont_practice_will1.gif\">Practice Problem #1<\/a><\/li>\n<li><a href=\"https:\/\/www.e-education.psu.edu\/files\/geog482\/image\/cont_practice_will6.gif\">Practice Problem #1 Solution<\/a><\/li>\n<li><a href=\"https:\/\/www.e-education.psu.edu\/files\/geog482\/image\/cont_practice_pitt1.gif\">Practice Problem #2<\/a><\/li>\n<li><a href=\"https:\/\/www.e-education.psu.edu\/files\/geog482\/image\/cont_practice_pitt6.gif\">Practice Problem #2 Solution<\/a><\/li>\n<\/ul>\n<p>You will be asked to demonstrate your contouring ability again in the Lesson 7 Quiz and in the final exam.<\/p>\n<p>Kevin Sabo (personal communication, Winter 2002) remarked that \u201cIf you were unfortunate enough to be hand-contouring data in the 1960\u2032s and 70\u2032s, you may at least have had the aid of a\u00a0<a href=\"http:\/\/www.nzeldes.com\/HOC\/Gerber.htm\">Gerber Variable Scale<\/a>. After hand contouring in Lesson 7, I sure wished I had my Gerber!\u201d<\/p>\n<h3><strong>PRACTICE\u00a0<\/strong><strong>QUIZ<\/strong><\/h3>\n<p>Registered Penn State students should return now to the Chapter 7 folder\u00a0 in ANGEL (via the Resources menu to the left) to take a self-assessment quiz about\u00a0<strong>Contouring<\/strong>. You may take practice quizzes as many times as you wish. They are not scored and do not affect your grade in any way.<\/p>\n<h2>7.7. Digital Line Graph (DLG)<\/h2>\n<h3>IDENTIFICATION<\/h3>\n<p>Digital Line Graphs (DLGs) are vector representations of most of the features and attributes shown on USGS topographic maps. Individual feature sets (outlined in the table below) are encoded in separate digital files. DLGs exist at three scales: small (1:2,000,000), intermediate (1:100,000) and large (1:24,000). Large-scale DLGs are produced in<strong>tiles<\/strong>\u00a0that correspond to the 7.5-minute topographic quadrangles from which they were derived.<\/p>\n<p>&nbsp;<\/p>\n<table summary=\"Layers and contents of large-scale Digital Line Graph files\">\n<caption>Description of Digital Line Graph Layers<\/caption>\n<thead>\n<tr>\n<th><strong>Layer<\/strong><\/th>\n<th><strong>Features<\/strong><\/th>\n<\/tr>\n<\/thead>\n<tbody>\n<tr>\n<td>Public Land Survey System (PLSS)<\/td>\n<td>Township, range, and section lines<\/td>\n<\/tr>\n<tr>\n<td>Boundaries<\/td>\n<td>State, county, city, and other national and State lands such as forests and parks<\/td>\n<\/tr>\n<tr>\n<td>Transportation<\/td>\n<td>Roads and trails, railroads, pipelines and transmission lines<\/td>\n<\/tr>\n<tr>\n<td>Hydrography<\/td>\n<td>Flowing water, standing water, and wetlands<\/td>\n<\/tr>\n<tr>\n<td>Hypsography<\/td>\n<td>Contours and supplementary spot elevations<\/td>\n<\/tr>\n<tr>\n<td>Non-vegetative features<\/td>\n<td>Glacial moraine, lava, sand, and gravel<\/td>\n<\/tr>\n<tr>\n<td>Survey control and markers<\/td>\n<td>Horizontal and vertical monuments (third order or better)<\/td>\n<\/tr>\n<tr>\n<td>Man-made features<\/td>\n<td>Cultural features, such as building, not collected in other data categories<\/td>\n<\/tr>\n<tr>\n<td>Woods, scrub, orchards, and vineyards<\/td>\n<td>Vegetative surface cover<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p>Layers and contents of large-scale Digital Line Graph files. Not all layers available for all quadrangles (USGS, 2006).<\/p>\n<p>&nbsp;<\/p>\n<p>&nbsp;<\/p>\n<p><img decoding=\"async\" alt=\"Image of Digital Line Graph hyspography, hydrography, and transportation layers viewed in Global Mapper software\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/dlg_bushkill_sw.gif\" \/><\/p>\n<p>Portion of three Digital Line Graph (DLG) layers for USGS Bushkill, PA quadrangle; imaged with Global Mapper (dlgv32 Pro) software. Transportation features are arbitrarily colored red, hydrography blue, and hypsography brown. The square symbols are nodes and the triangles represent polygon centroids.<\/p>\n<h4>DATA QUALITY<\/h4>\n<p>Like other USGS data products, DLGs conform to National Map Accuracy Standards. In addition, however, DLGs are tested for the logical consistency of the topological relationships among data elements. Similar to the Census Bureau\u2019s TIGER\/Line, line segments in DLGs must begin and end at point features (nodes), and line segments must be bounded on both sides by area features (polygons).<\/p>\n<h4>SPATIAL REFERENCE INFORMATION<\/h4>\n<p>DLGs are heterogenous. Some use UTM coordinates, others State Plane Coordinates. Some are based on NAD 27, others on NAD 83. Elevations are referenced either to NGVD 29 or NAVD 88 (USGS, 2006a).<\/p>\n<h4>ENTITIES AND ATTRIBUTES<\/h4>\n<p>The basic elements of DLG files are nodes (positions), line segments that connect two nodes, and areas formed by three or more line segments. Each node, line segment, and area is associated with two-part integer attribute codes. For example, a line segment associated with the attribute code \u201c050 0412\u2033 represents a hydrographic feature (050), specifically, a stream (0412).<\/p>\n<h4>DISTRIBUTION<\/h4>\n<p>Not all DLG layers are available for all areas at all three scales. Coverage is complete at 1:2,000,000.\u00a0 At the intermediate scale, 1:100,000 (30 minutes by 60 minutes), all hydrography and transportation files are available for the entire U.S., and complete national coverage is planned. At 1:24,000 (7.5 minutes by 7.5 minutes), coverage remains spotty. The files are in the public domain, and can be used for any purpose without restriction.<\/p>\n<p>Large- and Intermediate -scale DLGs are available for download through\u00a0<a href=\"http:\/\/earthexplorer.usgs.gov\/\">EarthExplorer system<\/a>. You can plot 1:2,000,000 DLGs on-line at the\u00a0<a href=\"http:\/\/nationalatlas.gov\/\">USGS\u2019 National Atlas of the United States<\/a>.<\/p>\n<h3><\/h3>\n<h3>DIGITAL LINE GRAPH HYPSOGRAPHY<\/h3>\n<p>In one sense, DLGs are as much \u201clegacy\u201d data as the out-of-date topographic maps from which they were produced. Still, DLG data serve as primary or secondary sources for several themes in the USGS National Map, including hydrography, boundaries, and transportation. DLG hypsography data are not included in the National Map, however. It is assumed that GIS users can generate elevation contours as needed from DEMs. DLG hypsography and hydrography layers are the preferred sources from which USGS DEMs are produced, however.<\/p>\n<p><img decoding=\"async\" alt=\"Digital Line Graph hypsography and hydrography viewed in Global Mapper software\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/dlg_hypso.gif\" \/><\/p>\n<p>Portion of the hypsography and hydrography layers of a large-scale Digital Line Graph (DLG). USGS Bushkill, PA quadrangle; imaged with Global Mapper (dlgv32 Pro) software.<\/p>\n<p><strong>Hypsography<\/strong>\u00a0refers to the measurement and depiction of the terrain surface, specifically with contour lines. Several different methods have been used to produce DLG hypsography layers, including:<\/p>\n<ul>\n<li>Scanning contour lines on photographic film or paper maps, converting the scanned raster data to vectors, then editing and attributing the vector features;<\/li>\n<li>Manually digitizing and attributing contour lines on photographic film or paper maps; and<\/li>\n<li>Producing contours by photogrammetric processes.<\/li>\n<\/ul>\n<p>The preferred method is to manually digitize contour lines in vector mode, then to key-enter the corresponding elevation attribute data.<\/p>\n<p><img decoding=\"async\" alt=\"Attributes of a contour line in a DLG hypsography layer, viewed in Global Mapper software\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/dlg_hypso_attribute.gif\" \/><\/p>\n<p>The highlighted contour line has been selected, and its attributes reported in a Global Mapper window. Notice that the line feature is attributed with a unique Element ID code (LE01, 639) and an elevation (1000 feet).<\/p>\n<p>&nbsp;<\/p>\n<h3><strong>TRY THIS!<\/strong><\/h3>\n<h3><strong>EXPLORING DLGS WITH GLOBAL MAPPER (DLGV32 PRO)<\/strong><\/h3>\n<p>Now I\u2019d like you to use Global Mapper (or dlgv32 Pro) software to investigate the characteristics of the hypsography layer of a USGS Digital Line Graph (DLG). The instructions below assume that you have already installed software on your computer. (If you haven\u2019t, return to\u00a0<a href=\"https:\/\/www.e-education.psu.edu\/natureofgeoinfo\/c6_p6.html\">installation instructions<\/a>\u00a0presented earlier in Chapter 6). First you\u2019ll download and a sample DLG file. In a following activity you\u2019ll have a chance to find and download DLG data for your area.<\/p>\n<ol>\n<li>If you haven\u2019t done so already, create a directory called \u201cUSGS Data\u201d on your hard disk, where you file your course materials.<\/li>\n<li>Next,\u00a0<a href=\"https:\/\/courseware.e-education.psu.edu\/downloads\/natureofgeoinfo\/DLG.zip\">Download the\u00a0<strong>DLG.zip<\/strong>\u00a0data archive<\/a>. The ZIP archive is 1.2 Mb in size and will take approximately 15 seconds to download via high speed DSL or cable, or about 4 minutes and 15 seconds minutes via 56 Kbps modem.<\/li>\n<li>Now decompress the archive into a directory on your hard disk.\n<ul>\n<li>Open the archive DLG.zip.<\/li>\n<li>Create a subdirectory called \u201cDLG\u201d within the directory in which you save data for this class.<\/li>\n<li>Extract all files in the ZIP archive into your new subdirectory.<\/li>\n<\/ul>\n<p>The end result will be five subdirectories, each of which includes the data files that make up a DLG \u201clayer,\u201d along with a master directory.<\/li>\n<li>Launch Global Mapper or dlgv32 Pro.<\/li>\n<li>Open a Digital Line Graph by choosing File &gt; Open as New\u2026, then navigate to the directory \u201cDLG\/Hypso.\u201d Open the file \u2018Hp01catd.ddf\u2019 (you can open up to four files at once in the trial version of Global Mapper.) The data correspond with the 7.5 minute quadrangle for Bushkill, PA. The file is encoded in Spatial Data Transfer Standard (SDTS) format. For information about SDTS, see the\u00a0<a href=\"https:\/\/www.e-education.psu.edu\/files\/geog482\/file\/sdts-tutorial.pdf\">SDTS Tutorial<\/a>(PDF format).<\/li>\n<li>Global Mapper may ask you to direct it to a \u2018Master Data Dictionary\u2019 file. If so, navigate to, and select, the file \u2018Dlg\/MasterDlg\/Dlg3mdir.ddf\u2019<\/li>\n<li>Experiment with Global Mapper\u2019s tools. Use Zoom and Pan to magnify and scroll across the DLG. The\u00a0<strong>Full View button<\/strong>\u00a0(the one with the house icon) refreshes the initial full view of the data set.<\/li>\n<li>The<strong>\u00a0Feature Info tool<\/strong>\u00a0allows you to query the attributes of a particular feature. Try clicking a single line segment. Note that you can display the attributes of a feature in the lower left portion of the application window by simply hovering over the feature.<\/li>\n<li>The\u00a0<strong>Measure tool<\/strong>\u00a0(ruler icon) allows you to not only measure distance as the crow flies, but also to see the area enclosed by a series of line segments drawn by repeated mouse clicks. Note again the location information that is given to you near the bottom of the application window.<\/li>\n<li>Certain tools, e.g., the\u00a0<strong>3D Path Profile\/Line of Sight tool<\/strong>\u00a0(next to the Feature Info tool) are not functional in the free (unregistered) version of Global Mapper.<\/li>\n<li>The trial version of Global Mapper allows you to open and view up to four files at once. You might find it interesting to open and compare the Bushkill DLG hypsography file and the corresponding DRG you viewed in Lesson 6. Note that you can turn layers on and off, and even adjust their transparency at Tools &gt; Control Center.\u00a0<strong>How do the contours in the DLG compare with those in the DRG? What explains the difference?<\/strong><\/li>\n<li>Global Mapper provides the metadata you\u2019ll need to answer questions in a practice quiz.\u00a0<strong>To access the metadata, navigate to Tools &gt; Control Center, then click the Metadata button.<\/strong><\/li>\n<\/ol>\n<h2>7.8. Digital Elevation Model (DEM)<\/h2>\n<p>The term \u201cDigital Elevation Model\u201d has both generic and specific meanings. In general, a DEM is any raster representation of a terrain surface. Specifically, a DEM is a data product of the U.S. Geological Survey. Here we consider the characteristics of DEMs produced by the USGS Later in this chapter we\u2019ll consider sources of global terrain data.<\/p>\n<h3>IDENTIFICATION<\/h3>\n<p>USGS DEMs are raster grids of elevation values that are arrayed in series of south-north profiles. Like other USGS data, DEMs were produced originally in tiles that correspond to topographic quadrangles. Large scale (7.5-minute and 15-minute), intermediate scale (30 minute), and small scale (1 degree) series were produced for the entire U.S. The resolution of a DEM is a function of the east-west spacing of the profiles and the south-north spacing of elevation points within each profile.<\/p>\n<p>DEMs corresponding to 7.5-minute quadrangles are available at 10-meter resolution for much, but not all, of the U.S. Coverage is complete at 30-meter resolution. In these large scale DEMs elevation profiles are aligned parallel to the central meridian of the local UTM zone, as shown in the illustration below.\u00a0<strong>See how the DEM tile in the illustration below appears to be tilted? This is because the corner points are defined in unprojected geographic coordinates that correspond to the corner points of a USGS quadrangle<\/strong>. The farther the quadrangle is from the central meridian of the UTM zone, the more it is tilted.<\/p>\n<p>&nbsp;<\/p>\n<p><img decoding=\"async\" alt=\"Schematic illustration of elevation profiles in a USGS digital elevation model\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/dem_structure.png\" \/><\/p>\n<p>Arrangement of elevation profiles in a large scale USGS Digital Elevation Model (USGS, 1987).<\/p>\n<p>&nbsp;<\/p>\n<p>As shown below, the arrangement of the elevation profiles is different in intermediate- and small-scale DEMs. Like meridians in the northern hemisphere, the profiles in 30-minute and 1-degree DEMs converge toward the north pole. For this reason the resolution of intermediate- and small-scale DEMs (that is to say, the spacing of the elevation values) is expressed differently than for large-scale DEMs.\u00a0<strong>The resolution of 30-minute DEMs is said to be 2 arc seconds and 1-degree DEMs are 3 arc seconds<\/strong>. Since an arc second is 1\/3600 of a degree, elevation values in a 3 arc second DEM are spaced 1\/1200 degree apart, representing a grid cell about 66 meters \u201cwide\u201d by 93 meters \u201ctall\u201d at 45\u00ba latitude.<\/p>\n<p><img decoding=\"async\" alt=\"Schematic illustration of elevation profiles in a small-scale USGS digital elevation model\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/dem_structure_smallscale.png\" \/><\/p>\n<p>Arrangement of elevation profiles in a small scale USGS Digital Elevation Model (USGS, 1987).<\/p>\n<p><strong>The preferred method for producing the elevation values that populate DEM profiles is interpolation from DLG hypsography and hydrography layers<\/strong>\u00a0(including the hydrography layer enables analysts to delineate valleys with less uncertainty than hypsography alone). Some older DEMs were produced from elevation contours digitized from paper maps or during photogrammetric processing, then smoothed to filter out errors. Others were produced photogrammtrically from aerial photographs.<\/p>\n<h4>DATA QUALITY<\/h4>\n<p>The vertical accuracy of DEMs is expressed as the root mean square error (RMSE) of a sample of at least 28 elevation points. The target accuracy for large-scale DEMs is seven meters; 15 meters is the maximum error allowed.<\/p>\n<h4>SPATIAL REFERENCE INFORMATION<\/h4>\n<p>Like DLGs, USGS DEMs are heterogenous. They are cast on the Universal Transverse Mercator projection used in the local UTM zone. Some DEMs are based upon the North American Datum of 1983, others on NAD 27. Elevations in some DEMs are referenced to either NGVD 29 or NAVD 88.<\/p>\n<h4>ENTITIES AND ATTRIBUTES<\/h4>\n<p>Each record in a DEM is a profile of elevation points. Records include the UTM coordinates of the starting point, the number of elevation points that follow in the profile, and the elevation values that make up the profile. Other than the starting point, the positions of the other elevation points need not be encoded, since their spacing is defined. (Later in this lesson you\u2019ll download a sample USGS DEM file. Try opening it in a text editor to see what I\u2019m talking about.)<\/p>\n<h4>DISTRIBUTION<\/h4>\n<p>DEM tiles are available for free download through many state and regional clearinghouses. You can find these sources by searching the<a href=\"http:\/\/geo.data.gov\/geoportal\/catalog\/main\/home.page\">GEODATA portion of the Data.Gov site<\/a>, formerly the separate Geospatial One Stop site.<\/p>\n<p>As part of its National Map initiative, the USGS has developed a \u201cseamless\u201d\u00a0<a href=\"http:\/\/ned.usgs.gov\/\">National Elevation Dataset<\/a>\u00a0that is derived from DEMs, among other sources. NED data are available at three resolutions: 1 arc second (approximately 30 meters), 1\/3 arc second (approximately 10 meters), and 1\/9 arc second (approximately 3 meters). Coverage ranges from complete at 1 arc second to extremely sparse at 1\/9 arc second. An extensive\u00a0<a href=\"http:\/\/ned.usgs.gov\/Ned\/faq.asp\">FAQ on NED data is published here<\/a>. The second of the two following activities involves downloading NED data and viewing it in Global Mapper.<\/p>\n<h3><strong>TRY THIS!<\/strong><\/h3>\n<h3><strong>EXPLORING DEMS WITH GLOBAL MAPPER (DLGV32 PRO)<\/strong><\/h3>\n<p>Global Mapper time again! This time you\u2019ll investigate the characteristics of a USGS DEM. The instructions below assume that you have already installed the software on your computer. (If you haven\u2019t, return to\u00a0<a href=\"https:\/\/www.e-education.psu.edu\/natureofgeoinfo\/c6_p6.html\">installation instructions<\/a>\u00a0presented earlier in Chapter 6). The instructions will remind you how to open a DEM in dlgv32 Pro. In the practice quiz that follows you\u2019ll be asked questions require you to explore the data for answers.<\/p>\n<ol>\n<li>First\u00a0<a href=\"https:\/\/courseware.e-education.psu.edu\/downloads\/natureofgeoinfo\/DEM.zip\">Download the\u00a0<strong>DEM.zip<\/strong>\u00a0data archive<\/a>. The ZIP archive is 2.5 Mb in size and will take about 30 seconds to download via high speed DSL or cable, or nearly 9 minutes via 56 Kbps modem. If you can\u2019t download the file, contact my teaching assistant or me right away so we can help you resolve the problem.<\/li>\n<li>Now decompress the archive into a directory on your hard disk.\n<ul>\n<li>Open the archive<strong>\u00a0DEM.zip.<\/strong><\/li>\n<li>Create a subdirectory called \u201cDEM\u201d within the directory in which you save class data.<\/li>\n<li>Extract all files in the ZIP archive into your new subdirectory.<\/li>\n<\/ul>\n<p>The end result will be two subdirectories, one of which contains a 30-meter DEM, the other a 10-meter DEM. These datasets are in the earlier distribution format of USGS DEM data \u2014 elevation data in horizontal (pixel) units of meters and representative of the area covered by a 1:24,000 topo map sheet. In the Try This that follows this one you will see that the distribution format options have expanded.<\/li>\n<li>Launch Global Mapper.<\/li>\n<li>Open a Digital Elevation Model by choosing\u00a0<strong>File\u00a0<\/strong>&gt;\u00a0<strong>Open Data File(s)\u2026<\/strong>, then navigate to the directory DEM_30m or DEM_10m, then open the file\u00a0<em>bushkill_pa.dem<\/em><\/li>\n<li>Use the\u00a0<strong>Zoom\u00a0<\/strong>and\u00a0<strong>Pan\u00a0<\/strong>tools to magnify and scroll across the DEM. The<strong>\u00a0Full View<\/strong>\u00a0button (house icon) refreshes the initial full view of the dataset.<\/li>\n<li>Global Mapper provides access to the metadata you\u2019ll need to answer questions in a practice quiz. To access the metadata, navigate to<strong>\u00a0Tools\u00a0<\/strong>&gt;<strong>\u00a0Control Center<\/strong>, then click the<strong>\u00a0Metadata\u00a0<\/strong>button.<\/li>\n<\/ol>\n<p>You can change the appearance of the DEM in the Options section of the Control Center. You can also alter the appearance of the DEM by choosing Tools &gt; Configure, and changing the settings in, especially, Vertical Options and Shader Options. To see the DEM data with(out) hill shading, find the Enable\/Disable Hill Shading button on the Shader toolbar (it has a sunburst in the lower left corner).<\/p>\n<h3><strong>TRY THIS!<\/strong><\/h3>\n<h3><strong>DOWNLOAD YOUR OWN NATIONAL ELEVATION DATASET (NED) DATA<\/strong><\/h3>\n<ol>\n<li>Go to the\u00a0<a href=\"http:\/\/nationalmap.gov\/elevation.html\">Elevation page of the USGS National Map site<\/a>.<br \/>\nRead any of the information there that you choose.<\/li>\n<li>Follow the link to\u00a0<a href=\"http:\/\/viewer.nationalmap.gov\/viewer\/\">The National Map Viewer<\/a>.<br \/>\nIf you wish, also follow the link to the detailed instructions that you see under The National Map Viewer link. The following instructions in this Try This should also suffice, and perhaps elaborate at bit more.You will not see Elevation data listed in the left hand Overlays pane. The USGS is in the process of adding more visualization options when it comes to elevation data. See the Hill Shade button at the upper right of the map area.<\/li>\n<li>Use the GIS tools, found above the map area, to pan to and zoom in on an area of interest. Then, click the\u00a0<strong>Download Data<\/strong>\u00a0button.<br \/>\nChoose a\u00a0<em>reference area<\/em>\u00a0from the pick list. The default is the index of the areas covered by the 1:24,000 topo map series.<br \/>\nYou could also chose the entire\u00a0<em>current map extent<\/em>, but depending upon your zoom level that may be a huge dataset.<br \/>\nThe instructions found via the Help link mentioned above mention that you can define an area based on creating a\u00a0<em>custom polygon<\/em>, but at this point I do not see how that is done\u2026<\/li>\n<li>Then click on the map to highlight a specific area of interest. A link to the available data sets will appear in the left hand pane under the<strong>Selection\u00a0<\/strong>tab. (The\u00a0<strong>All Results<\/strong>\u00a0button will list the multiple areas you click on.) Follow the\u00a0<strong>Download\u00a0<\/strong>link of the area you are interested in.<\/li>\n<li>In the\u00a0<strong>USGS Available Data<\/strong>\u00a0window that opens, select\u00a0<strong>Elevation<\/strong>from the\u00a0<strong>Theme<\/strong>\u00a0column, and choose a file format from the pick list in the\u00a0<strong>Format\u00a0<\/strong>column.<br \/>\nI know that both the GeoTIFF and ArcGRID formats are compatible with Global Mapper. (ArcGRID is the Esri company\u2019s raster format. You\u2019ll be using the Esri software in future courses.)<br \/>\nClick the\u00a0<strong>Next\u00a0<\/strong>button.<\/li>\n<li>You will be given a list of Elevation\u00a0<em>Product\u00a0<\/em>choices.<br \/>\nThe choices that show\u00a0<em>Dynamic\u00a0<\/em>as the Type will be those that match the Format choice you made in the previous window. The other\u00a0<em>Staged<\/em>datasets are prepackaged and in the formats stated in the Product column, and apparently listed regardless of the Format choice you made.If multiple resolutions are available they will be listed.<\/p>\n<p>From the Product list, check the box for what you wish to download.<br \/>\nClick the\u00a0<strong>Next\u00a0<\/strong>button.<br \/>\nGo to the\u00a0<strong>Cart\u00a0<\/strong>pane on the right. (It may open automatically.)<br \/>\nGo to\u00a0<strong>Checkout<\/strong>, supply your e-mail address, and submit your order via the\u00a0<strong>Place Order<\/strong>\u00a0button.<br \/>\nYou will receive a message telling you that your order has been placed, and that will soon be followed by an e-mail regarding your order. About an hour and a half after I submitted my request I received a second e-mail containing a download link.<\/li>\n<li>The system will produce a ZIP archive that you can save to your hard disk (e.g., \u201c09647011.zip\u201d)<\/li>\n<li>Launch Global Mapper and open the ZIP archive. The software can read the data even in its compressed form; you should not need to extract the contents from the .zip file. (It would be a good idea to look at the contents of the .zip archive, though, if only to see the number and type of files included.)<\/li>\n<li>An image of the DEM data should appear in the Global Mapper window, similar to what you see shown below (even though the image below is from an older version of Global Mapper).<\/li>\n<li>Again, you can view the\u00a0<em>metadata\u00a0<\/em>associated with the DEM data via the\u00a0<strong>Tools\u00a0<\/strong>&gt;\u00a0<strong>Control Center<\/strong>\u00a0menu. Note the PIXEL dimensions reported in\u00a0<em>arc degrees<\/em>, as opposed to something like meters.<\/li>\n<\/ol>\n<p><img decoding=\"async\" alt=\"A portion of the National Elevation Dataset viewed in Global Mapper software\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/ned_global_mapper.png\" \/><\/p>\n<h3><strong>PRACTICE QUIZ<\/strong><\/h3>\n<p>Registered Penn State students should return now to the Chapter 7 folder\u00a0 in ANGEL (via the Resources menu to the left) to take a self-assessment quiz about\u00a0<strong>DLGs and DEMs<\/strong>. You may take practice quizzes as many times as you wish. They are not scored and do not affect your grade in any way.<\/p>\n<h2>7.9. Interpolation<\/h2>\n<p>DEMs are produced by various methods. The method preferred by USGS is to\u00a0<strong>interpolate<\/strong>\u00a0elevations grids from the hypsography and hydrography layers of Digital Line Graphs.<\/p>\n<p><img decoding=\"async\" alt=\"Hypsography and hydrography layers of a Digital Line Graph viewed in Global Mapper software\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/dem_dlg.jpg\" \/><\/p>\n<p>A USGS 7.5-minute DEM and the DLG hypsography and hydrography layers from which it was produced.<\/p>\n<p>The elevation points in DLG hypsography files are not regularly spaced. DEMs need to be regularly spaced to support the slope, gradient, and volume calculations they are often used for. Grid point elevations must be interpolated from neighboring elevation points. In the figure below, for example, the gridded elevations shown in purple were interpolated from the irregularly spaced spot elevations shown in red.<\/p>\n<p><img decoding=\"async\" alt=\"A grid of elevation values that were interpolated from an irregularly-spaced array\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/interp_grid.gif\" \/><\/p>\n<p>Elevation values in DEMs are interpolated from irregular arrays of elevations measured through photogrammetric methods, or derived from existing DLG hypsography and hydrography data.<\/p>\n<p>Here\u2019s another example of interpolation for mapping. The map below shows how 1995 average surface air temperature differed from the average temperature over a 30-year baseline period (1951-1980). The temperature anomalies are depicted for grid cells that cover 3\u00b0 longitude by 2.5\u00b0 latitude.<\/p>\n<p><img decoding=\"async\" alt=\"A map showing gridded temperature data\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/temp_95.jpg\" \/><\/p>\n<p>1995 Surface Temperature Anomalies. (National Climatic Data Center, 2005).<\/p>\n<p>The gridded data shown above were estimated from the temperature records associated with the very irregular array of 3,467 locations pinpointed in the map below.\u00a0 The irregular array is transformed into a regular array through interpolation. In general,\u00a0<strong>interpolation is the process of estimating an unknown value from neighboring known values.<\/strong><\/p>\n<p><img decoding=\"async\" alt=\"Locations of temperature climate records used to create a gridded temperature map\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/temp_stations.jpg\" \/><\/p>\n<p>The Global Historical Climate Network. (Eischeid et al., 1995).<\/p>\n<p>Elevation data are often not measured at evenly-spaced locations. Photogrammetrists typically take more measurements where the terrain varies the most. They refer to the dense clusters of measurements they take as \u201cmass points.\u201d Topographic maps (and their derivatives, DLGs) are another rich source of elevation data. Elevations can be measured from contour lines, but obviously contours do not form evenly-spaced grids. Both methods give rise to the need for interpolation.<\/p>\n<p><img decoding=\"async\" alt=\"Three number lines illustrating how interpolation is affected by assumptions about the underlying distribution\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/linear_interp.gif\" \/><\/p>\n<p>&nbsp;<\/p>\n<p>Interpolating an intermediate value on a number line.<\/p>\n<p>The illustration above shows three number lines, each of which ranges in value from 0 to 10. If you were asked to interpolate the value of the tick mark labeled \u201c?\u201d on the top number line, what would you guess? An estimate of \u201c5\u2033 is reasonable, provided that the values between 0 and 10 increase at a constant rate. If the values increase at a geometric rate, the actual value of \u201c?\u201d could be quite different, as illustrated in the bottom number line. The validity of an interpolated value depends, therefore, on the validity of our assumptions about the nature of the underlying surface.<\/p>\n<p>As I mentioned in Chapter 1, the surface of the Earth is characterized by a property called\u00a0<strong>spatial dependence<\/strong>. Nearby locations are more likely to have similar elevations than are distant locations. Spatial dependence allows us to assume that it\u2019s valid to estimate elevation values by interpolation.<\/p>\n<p>Many interpolation algorithms have been developed. One of the simplest and most widely used (although often not the best) is the<strong>inverse distance weighted<\/strong>\u00a0algorithm. Thanks to the property of spatial dependence, we can assume that estimated elevations are more similar to nearby elevations than to distant elevations. The inverse distance weighted algorithm estimates the value z of a point P as a function of the z-values of the nearest n points. The more distant a point, the less it influences the estimate.<\/p>\n<p><img decoding=\"async\" alt=\"Diagram and formula explaining inverse distance weighted interpolation\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/inv_dist_interp.gif\" \/><\/p>\n<p>The inverse distance weighted interpolation procedure.<\/p>\n<h3><strong>PRACTICE QUIZ<\/strong><\/h3>\n<p>Registered Penn State students should return now to the Chapter 7 folder\u00a0 in ANGEL (via the Resources menu to the left) to take a self-assessment quiz about\u00a0<strong>Interpolation<\/strong>. You may take practice quizzes as many times as you wish. They are not scored and do not affect your grade in any way.<\/p>\n<h2>7.10. Slope<\/h2>\n<p>Slope is a measure of change in elevation. It is a crucial parameter in several well-known predictive models used for environmental management, including the Universal Soil Loss Equation and agricultural non-point source pollution models.<\/p>\n<p>One way to express slope is as a percentage. To calculate percent slope, divide the difference between the elevations of two points by the distance between them, then multiply the quotient by 100. The difference in elevation between points is called the rise. The distance between the points is called the run. Thus, percent slope equals (rise \/ run) x 100.<\/p>\n<p><img decoding=\"async\" alt=\"Diagram illustrating how slope may be calculated as a percentage\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/slope_percent.gif\" \/><\/p>\n<p>Calculating percent slope. A rise of 100 feet over a run of 100 feet yields a 100 percent slope. A 50-foot rise over a 100-foot run yields a 50 percent slope.<\/p>\n<p>&nbsp;<\/p>\n<p>Another way to express slope is as a slope angle, or degree of slope. As shown below, if you visualize rise and run as sides of a right triangle, then the degree of slope is the angle opposite the rise. Since degree of slope is equal to the tangent of the fraction rise\/run, it can be calculated as the arctangent of rise\/run.<\/p>\n<p><img decoding=\"async\" alt=\"Illustration showing how slope may be calculated in degrees\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/slope_degree.gif\" \/><\/p>\n<p>A rise of 100 feet over a run of 100 feet yields a 45\u00b0 slope angle. A rise of 50 feet over a run of 100 feet yields a 26.6\u00b0 slope angle.<\/p>\n<p>&nbsp;<\/p>\n<p>You can calculate slope on a contour map by analyzing the spacing of the contours. If you have many slope values to calculate, however, you will want to automate the process. It turns out that slope calculations are much easier to calculate for gridded elevation data than for vector data, since elevations are more or less equally spaced in raster grids.<\/p>\n<p>Several algorithms have been developed to calculate percent slope and degree of slope. The simplest and most common is called the<strong>neighborhood method<\/strong>. The neighborhood method calculates the slope at one grid point by comparing the elevations of the eight grid points that surround it.<\/p>\n<p><img decoding=\"async\" alt=\"Illustration showing how slope at one point is calculated as a function of the elevations of eight surrounding points\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/slope_grid.gif\" \/><\/p>\n<p>The neighborhood algorithm estimates percent slope in cell 5 by comparing the elevations of neighboring grid cells.<\/p>\n<p>The neighborhood algorithm estimates percent slope at grid cell 5 (Z5) as the sum of the absolute values of east-west slope and north-south slope, and multiplying the sum by 100. The diagram below illustrates how east-west slope and north-south slope are calculated. Essentially, east-west slope is estimated as the difference between the sums of the elevations in the first and third columns of the 3 x 3 matrix. Similarly, north-south slope is the difference between the sums of elevations in the first and third rows (note that in each case the middle value is weighted by a factor of two).<\/p>\n<p><img decoding=\"async\" alt=\"Algorithm for calculating slope with gridded elevation data\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/slope_algorithm.gif\" \/><\/p>\n<p>The neighborhood algorithm for calculating percent slope.<\/p>\n<p>The neighborhood algorithm calculates slope for every cell in an elevation grid by analyzing each 3 x 3 neighborhood. Percent slope can be converted to slope degree later. The result is a grid of slope values suitable for use in various soil loss and hydrologic models.<\/p>\n<h2>7.11. Relief Shading<\/h2>\n<p>You can see individual pixels in the zoomed image of a 7.5-minute DEM below. I used dlgv32 Pro\u2019s \u201cGradient Shader\u201d to produce the image. Each pixel represents one elevation point. The pixels are shaded through 256 levels of gray. Dark pixels represent low elevations, light pixels represent high ones.<\/p>\n<p><img decoding=\"async\" alt=\"Image of a DEM with pixels shaded light to dark in proportion to elevation\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/dem.gif\" \/><\/p>\n<p>A digital elevation model in which light pixels represent high elevations, and dark pixels represent low elevations.<\/p>\n<p>It\u2019s also possible to assign gray values to pixels in ways that make it appear that the DEM is illuminated from above. The image below, which shows the same portion of the Bushkill DEM as the image above, illustrates the effect, which is called\u00a0<strong>terrain shading<\/strong>,\u00a0<strong>hill shading<\/strong>\u00a0or\u00a0<strong>shaded relief.<\/strong><\/p>\n<p><strong><img decoding=\"async\" alt=\"Image of a DEM shaded as though it were illuminated from above\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/dem_shaded.jpg\" \/><\/strong><\/p>\n<p>Shaded terrain image produced from the same DEM as shown in the above figures, using dlgv32 Pro\u2019s Daylight Shader option, with the Surface Color set to gray.<\/p>\n<p>The appearance of a shaded terrain image depends on several parameters, including\u00a0<strong>vertical exaggeration<\/strong>. Click the buttons under the image below to compare the four terrain images of North America shown below, in which elevations are exaggerated 5 times, 10 times, 20 times, and 40 times respectively. (You will need to have the Adobe Flash player installed in order to complete this exercise. If you do not already have the Flash player, you can\u00a0<a href=\"http:\/\/www.adobe.com\/shockwave\/download\/index.cgi?P1_Prod_Version=ShockwaveFlash\">download it for free from Adobe<\/a>.)<\/p>\n<p>Effects of vertical exaggeration on a shaded terrain image<\/p>\n<p>Another influential parameter is the\u00a0<strong>angle of illumination<\/strong>. Click the buttons to compare terrain images that have been illuminated from the northeast, southeast, southwest, and northwest. Does the terrain appear to be inverted in one or more of the images?\u00a0<strong>To minimize the possibility of terrain inversion, it is conventional to illuminate terrain from the northwest.<\/strong><\/p>\n<p>Effects of illumination angle on a shaded terrain image.<\/p>\n<h2>7.12. Lidar<\/h2>\n<p>For many applications, 30-meter DEMs whose vertical accuracy is measured in meters are simply not detailed enough. Greater accuracy and higher horizontal resolution can be produced by photogrammetric methods, but precise photogrammetry is often too time-consuming and expensive for extensive areas. Lidar is a digital remote sensing technique that provides an attractive alternative.<\/p>\n<p>Lidar stands for\u00a0<strong>LIght Detection And Ranging<\/strong>. Like radar (RAdio Detecting And Ranging), lidar instruments transmit and receive energy pulses, and enable distance measurement by keeping track of the time elapsed between transmission and reception. Instead of radio waves, however, lidar instruments emit laser light (laser stands for Light Amplifications by Stimulated Emission of Radiation).<\/p>\n<p>Lidar instruments are typically mounted in low altitude aircraft. They emit up to 5,000 laser pulses per second, across a ground swath some 600 meters wide (about 2,000 feet). The ground surface, vegetation canopy, or other obstacles reflect the pulses, and the instrument\u2019s receiver detects some of the backscatter. Lidar mapping missions rely upon GPS to record the position of the aircraft, and upon inertial navigation instruments (gyroscopes that detect an aircraft\u2019s pitch, yaw, and roll) to keep track of the system\u2019s orientation relative to the ground surface.<\/p>\n<p>In ideal conditions, lidar can produce DEMs with 15-centimeter vertical accuracy, and horizontal resolution of a few meters. Its cost is prohibitive for small missions, but is justified for larger projects in which detail is essential. For example, lidar has been used successfully to detect subtle changes in the thickness of the Greenland ice sheet that result in a net loss of over 50 cubic kilometers of ice annually.<\/p>\n<p><img decoding=\"async\" alt=\"Composite lidar image showing changes in thickness of Greenland ice sheet\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/lidar_greenland.jpg\" \/><\/p>\n<p>Image of Greenland, viewed from the south, showing changes in ice thickness measured by airborne lidar. Ice sheet thickness decreasing at 40-60 cm per year in darker blue areas (Goddard Space Flight Center, n.d.).<\/p>\n<p>To learn more about the use of lidar in mapping changes in the Greenland ice sheet, visit\u00a0<a href=\"http:\/\/svs.gsfc.nasa.gov\/stories\/greenland\/\">NASA\u2019s Scientific Visualization Studio<\/a>.<\/p>\n<h2>7.13. Global Elevation Data<\/h2>\n<p>This page profiles three data products that include elevation (and, in one case, bathymetry) data for all or most of the Earth\u2019s surface.<\/p>\n<h3>ETOPO1<\/h3>\n<p><img decoding=\"async\" alt=\"World map generated from ETOPO1 global terrain (with ice heights) and bathymetry data\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/etopo1.gif\" \/><br \/>\nShaded and colored terrain image produced from ETOPO1 data. (National Geophysical Data Center, 2009).<\/p>\n<p>ETOPO1 is a digital elevation model that includes both topography and bathymetry for the entire world. It consists of more than 233 million elevation values which are regularly spaced at 1 minute of latitude and longitude.\u00a0<strong>At the equator, the horizontal resolution of ETOPO1 is approximately 1.85 kilometers.\u00a0<\/strong>Vertical positions are specified in meters, and there are two versions of the dataset: one with elevations at the \u201cIce Surface\u201d of the Greenland and Antarctic ice sheets, and one with elevations at \u201cBedrock\u201d beneath those ice sheets. Horizontal positions are specified in geographic coordinates (decimal degrees). Source data, and thus data quality, vary from region to region.<br \/>\nYou can\u00a0<a href=\"http:\/\/www.ngdc.noaa.gov\/mgg\/global\/global.html\">download ETOPO1 data from the National Geophysical Data Center<\/a>.<\/p>\n<h3><\/h3>\n<h3>GTOPO30<\/h3>\n<p><img decoding=\"async\" alt=\"World terrain map generated from GTOPO30 data\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/gtopo30_tiles.gif\" \/><\/p>\n<p>Shaded and colored terrain image produced from GTOPO30 data. Data are distributed as 33 tiles (USGS, 2006b).<\/p>\n<p>GTOPO30 is a digital elevation model that extends over the world\u2019s land surfaces (but not under the oceans). GTOPO30 consists of more than 2.5 million elevation values, which are regularly spaced at 30 seconds of latitude and longitude.<strong>\u00a0At the equator, the resolution of GTOPO30 is approximately 0.925 kilometers \u2014 two times greater than ETOPO1<\/strong>. Vertical positions are specified to the nearest meter, and horizontal positions are specified in geographic coordinates. GTOPO30 data are distributed as tiles, most of which are 50\u00b0 in latitude by 40\u00b0 in longitude.<\/p>\n<p>GTOPO30 tiles are available for download from\u00a0<a href=\"http:\/\/eros.usgs.gov\/#\/Find_Data\/Products_and_Data_Available\/gtopo30_info\">USGS\u2019 EROS Data Center<\/a>. GTOPO60, a resampled and untiled version of GTOPO30, is available through the\u00a0<a href=\"http:\/\/seamless.usgs.gov\/\">USGS\u2019 Seamless Data Distribution Service<\/a>.<\/p>\n<h3><\/h3>\n<h3>SHUTTLE RADAR TOPOGRAPHY MISSION (SRTM)<\/h3>\n<p>From February 11 to February 22, 2000, the space shuttle Endeavor bounced radar waves off the Earth\u2019s surface, and recorded the reflected signals with two receivers spaced 60 meters apart. The mission measured the elevation of land surfaces between 60\u00b0 N and 57\u00b0 S latitude.\u00a0<strong>The highest resolution data products created from the SRTM mission are 30 meters<\/strong>. Access to 30-meter SRTM data for areas outside the U.S. are restricted by the National Geospatial-Intelligence Agency, which sponsored the project along with the National Aeronautics and Space Administration (NASA). A 90-meter SRTM data product is available for free download without restriction (Maune, 2007).<\/p>\n<p><img decoding=\"async\" alt=\"Anaglyph stereo image of terrain surface of Fiji mapping by Shuttle Radar Topography Mission\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/shuttle_stereo_fiji.jpg\" \/><\/p>\n<p>Anaglyph stereo image derived from Shuttle Radar Topography Mission data (NASA Jet Propulsion Laboratory, 2006).<\/p>\n<p>The image above shows Viti Levu, the largest of the some 332 islands that comprise the Sovereign Democratic Republic of the Fiji Islands. Viti Levu\u2019s area is 10,429 square kilometers (about 4000 square miles). Nakauvadra, the rugged mountain range running from north to south, has several peaks rising above 900 meters (about 3000 feet). Mount Tomanivi, in the upper center, is the highest peak at 1324 meters (4341 feet).<\/p>\n<p>Learn more about the Shuttle Radar Topography Mission at Web sites published by\u00a0<a href=\"http:\/\/www.jpl.nasa.gov\/srtm\">NASA<\/a>\u00a0and\u00a0<a href=\"http:\/\/srtm.usgs.gov\/mission.php\">USGS<\/a>.<\/p>\n<h2>7.14. Bathymetry<\/h2>\n<p>The term bathymetry refers to the process and products of measuring the depth of water bodies. The U.S. Congress authorized the comprehensive mapping of the nation\u2019s coasts in 1807, and directed that the task be carried out by the federal government\u2019s first science agency, the Office of Coast Survey (OCS). That agency is now responsible for mapping some 3.4 million nautical square miles encompassed by the 12-mile territorial sea boundary, as well as the 200-mile Exclusive Economic Zone claimed by the U.S., a responsibility that entails regular revision of about 1,000 nautical charts. The coastal bathymetry data that appears on USGS topographic maps, like the one shown below, is typically compiled from OCS charts.<\/p>\n<p><img decoding=\"async\" alt=\"Portion of topographic map showing ocean depths\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/bar_harbor_bathy.gif\" \/><\/p>\n<p>\u201dIsobaths\u201d (the technical term for lines of constant depth) shown on a USGS topographic map.<\/p>\n<p>Early hydrographic surveys involved sampling water depths by casting overboard ropes weighted with lead and marked with depth intervals called marks and deeps. Such ropes were called leadlines for the weights that caused them to sink to the bottom. Measurements were called soundings. By the late 19th century, piano wire had replaced rope, making it possible to take soundings of thousands rather than just hundreds of fathoms (a fathom is six feet).<\/p>\n<p><img decoding=\"async\" alt=\"\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/sounding.jpg\" \/><\/p>\n<p>Seaman paying out a sounding line during a hydrographic survey of the East coast of the U.S. in 1916. (NOAA, 2007).<\/p>\n<p>Echo sounders were introduced for deepwater surveys beginning in the 1920s.\u00a0<strong>Sonar\u00a0<\/strong>(SOund NAvigation and Ranging) technologies have revolutionized oceanography in the same way that aerial photography revolutionized topographic mapping. The seafloor topography revealed by sonar and related shipborne remote sensing techniques provided evidence that supported theories about seafloor spreading and plate tectonics.<\/p>\n<p>Below is an artist\u2019s conception of an oceanographic survey vessel operating two types of sonar instruments:\u00a0<strong>multibeam<\/strong>\u00a0and s<strong>ide scan sonar<\/strong>. On the left, a multibeam instrument mounted in the ship\u2019s hull calculates ocean depths by measuring the time elapsed between the sound bursts it emits and the return of echoes from the seafloor. On the right, side scan sonar instruments are mounted on both sides of a submerged \u201ctowfish\u201d tethered to the ship. Unlike multibeam, side scan sonar measures the strength of echoes, not their timing. Instead of depth data, therefore, side scanning produces images that resemble black-and-white photographs of the sea floor.<\/p>\n<p><img decoding=\"async\" alt=\"Illustration of sonar in use for bathymetric mapping\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/sonar.jpg\" \/><\/p>\n<p>Multibeam and side scan sonar in use for bathymetric mapping. (NOAA, 2002).<\/p>\n<p>A\u00a0<a href=\"http:\/\/craterlake.wr.usgs.gov\/bathymetry.html\">detailed report of the recent bathymetric survey of Crater Lake, Oregon, USA, is published by the USGS here<\/a>.<\/p>\n<h2>7.15. Statistical Surfaces<\/h2>\n<p>Strategies used to represent terrain surfaces can be used for other kinds of surfaces as well. For example, one of my first projects here at Penn State was to work with a distinguished geographer, the late Peter Gould, who was studying the diffusion of the Acquired Immune Deficiency Syndrome (AIDS) virus in the United States. Dr. Gould had recently published the map below.<\/p>\n<p><img decoding=\"async\" alt=\"Thematic map depicted HIV\/AIDS as a statistical surface\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/aids_surface.jpg\" \/><\/p>\n<p>&nbsp;<\/p>\n<p>Oblique view of contour lines representing distribution of AIDS cases in the U.S. 1988. (Gould, 1989. \u00a9 Association of American Geographers. All rights reserved. Reproduced here for educational purposes only).<\/p>\n<p>&nbsp;<\/p>\n<p>Gould portrayed the distribution of disease in the same manner as another geographer might portray a terrain surface. The portrayal is faithful to Gould\u2019s conception of the contagion as a continuous phenomenon. It was important to Gould that people understood that there was no location that did not have the potential to be visited by the epidemic. For both the AIDS surface and a terrain surface, a quantitative attribute (<em>z<\/em>) exists for every location (<em>x,y<\/em>). In general, when a\u00a0<strong>continuous phenomenon is conceived as being analogous to the terrain surface, the conception is called a statistical surface<\/strong>.<\/p>\n<h2>7.16. Theme: Hydrography<\/h2>\n<p>The NSDI Framework Introduction and Reference (FGDC, 1997) envisions the hydrography theme in this way:<\/p>\n<p>Framework hydrography data include surface water features such as lakes and ponds, streams and rivers, canals, oceans, and shorelines. Each of these features has the attributes of a name and feature identification code. Centerlines and polygons encode the positions of these features. For feature identification codes, many federal and state agencies use the Reach schedule developed by the U.S. Environmental Protection Agency (EPA).<\/p>\n<p>&nbsp;<\/p>\n<p>Many hydrography data users need complete information about connectivity of the hydrography network and the direction in which the water flows encoded in the data. To meet these needs, additional elements representing flows of water and connections between features may be included in framework data (p. 20).<\/p>\n<h3>IDENTIFICATION<\/h3>\n<p>FGDC had the\u00a0<strong>National Hydrography Dataset (NHD)<\/strong>\u00a0in mind when they wrote this description. NHD combines the vector features of<strong>Digital Line Graph (DLG) hydrography<\/strong>\u00a0with the\u00a0<strong>EPA\u2019s Reach files<\/strong>. Reaches are segments of surface water that share similar hydrologic characteristics. Reaches are of three types: transport, coastline, and waterbody. DLG lines features represent the transport and coastline types; polygon features are used to represent waterbodies. Every reach segment in the NHD is assigned a unique reach code, along with a host of other hydrological attributes including stream flow direction (which is encoded in the digitizing order of nodes that make up each segment), network connectivity, and feature names, among others. Because the order of reach codes are sequential from reach to reach, point-source data (such as a pollutant spill) can be geocoded to the affected reach. Used in this way, reaches comprise a<strong>linear referencing system<\/strong>\u00a0comparable to postal addresses along streets\u00a0 (USGS, 2002).<\/p>\n<p><img decoding=\"async\" alt=\"Diagram showing how water flow is attributed to reaches\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/nhd_flows.png\" \/><\/p>\n<p>How flow attributes are associated with reaches in the National Hydrography Dataset (USGS, 2000).<\/p>\n<p>NHD parses the U.S. surface drainage network into four hierarchical categories of units: 21 Regions, 222 Subregions, 352 Accounting units, and 2150 Cataloging units (also called Watersheds). Features can exist at multiple levels of the hierarchy, though they might not be represented in the same way. For example, while it might make the most sense to represent a given stream as a polygon features at the Watershed level, it may be more aptly represented as a line feature at the Region or Subregion level. NHD supports this by allowing multiple features to share the same reach codes. Another distinctive feature of NHD is\u00a0<strong>artificial flowlines<\/strong>\u2013centerline features that represent paths of water flow through polygon features such as standing water bodies. NHD is complex because it is designed to support sophisticated hydrologic modeling tasks, including point-source pollution modeling, flood potential, bridge construction, among others (Ralston, 2004).<\/p>\n<p><img decoding=\"async\" alt=\"Diagram illustrating how hydrographic features are represented with points, lines, and polygons\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/nhd_points_lines_areas.png\" \/><\/p>\n<p>How vector features are used to represent various types of reaches in the National Hydrography Dataset (USGS, 2000).<\/p>\n<p>NHD are available at three levels of detail (scale): medium (1:100,000, which is available for the entire U.S.), high (1:24,000, production of which is underway, \u201caccording to the availability of matching resources from NHD partners\u201d (USGS, 2002, p. 2), and local (larger scales such as 1:5,000), which \u201cis being developed where partners and data exist\u201d for select areas (USGS, 2006c; USGS, 2009; USGS 2013).<\/p>\n<h3>SPATIAL REFERENCE INFORMATION<\/h3>\n<p>NHD coordinates are decimal degrees referenced to the NAD 83 horizontal datum.<\/p>\n<h3>DISTRIBUTION<\/h3>\n<h3><strong>TRY THIS!<\/strong><\/h3>\n<h3>DOWNLOAD AND VIEW AN EXTRACT FROM THE NATIONAL HYDROGRAPHY DATASET<\/h3>\n<ol>\n<li>From the\u00a0<a href=\"http:\/\/nhd.usgs.gov\/\">NHD home page<\/a>\u00a0click the\u00a0<strong>Get Data<\/strong>\u00a0link and then follow the link to the\u00a0<strong>NHD Viewer<\/strong>. (There is a Help button next to the viewer link.)<\/li>\n<li>Use the GIS tools, found above the map area, to pan to and zoom in on an area of interest. Then, click the\u00a0<strong>Download Data<\/strong>\u00a0button.<br \/>\nChoose a\u00a0<em>reference area<\/em>\u00a0from the pick list. (You could also chose the entire\u00a0<em>current map extent<\/em>, but depending upon your zoom level that may be a huge dataset.)<br \/>\nThen click on the map to highlight a specific area of interest. A link to the available data sets will appear in the left hand pane under the<strong>Selection\u00a0<\/strong>tab. (The\u00a0<strong>All Results<\/strong>\u00a0button will list the multiple areas you click on.) Follow the\u00a0<strong>Download\u00a0<\/strong>link of the area you are interested in.<\/li>\n<li>In the\u00a0<strong>USGS Available Data<\/strong>\u00a0window that opens, select<strong>Hydrography\u00a0<\/strong>from the\u00a0<strong>Theme<\/strong>\u00a0column, and choose a file format from the pick list in the\u00a0<strong>Format\u00a0<\/strong>column. I chose\u00a0<strong>Shapefile<\/strong>\u00a0format for my extract, because I know it is compatible with Global Mapper \/ dlgv32 Pro. If you were working in ArcGIS you could choose the File Geodatabase option. Ralston (2004, p.187) observes that NHD \u201cis precisely the type of information that could benefit from an integrated data model in an object relational database.\u201d<br \/>\nClick the\u00a0<strong>Next\u00a0<\/strong>button.<\/li>\n<li>From the list of Hydrography Products check the box for what you wish to download. To be assured of getting the data in Shapefile format select a product referred to as a\u00a0<strong>Dynamic Extract<\/strong>.<br \/>\nClick the\u00a0<strong>Next\u00a0<\/strong>button.<br \/>\nGo to the\u00a0<strong>Cart\u00a0<\/strong>pane on the right. (It may open automatically.)<br \/>\nGo to\u00a0<strong>Checkout<\/strong>, supply your e-mail address, and submit your order via the\u00a0<strong>Place Order<\/strong>\u00a0button.<br \/>\nYou will receive a message telling you that your order has been placed, and that will soon be followed by an e-mail regarding your order. About an hour and a half after I submitted my request I received a second e-mail containing a download link.<\/li>\n<li>Extract the contents of the .zip file and view the Shapefile data set(s) in Global Mapper.<\/li>\n<li>Use the Identify pointer tool to reveal attributes of the reaches. In the example below I have highlighted a flowline associated with Cedar Creek in western Michigan.<\/li>\n<\/ol>\n<p><img decoding=\"async\" alt=\"Screenshot of the feature information window in Global Mapper\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/NHD_trythis_sample_Aug2011.jpg\" \/><\/p>\n<p>&nbsp;<\/p>\n<h2>7.17. Theme: Transportation<\/h2>\n<p>Transportation network data are valuable for all sorts of uses, including two we considered in Chapter 4: geocoding and routing. The Federal Geographic Data Committee (1997, p. 19) specified the following vector features and attributes for the transportation framework theme:<\/p>\n<table summary=\"Description of each transportation network feature\">\n<caption>Transportation Framework Attributes<\/caption>\n<thead>\n<tr>\n<th><strong>Feature<\/strong><\/th>\n<th><strong>Attributes<\/strong><\/th>\n<\/tr>\n<\/thead>\n<tbody>\n<tr>\n<td>Roads<\/td>\n<td>Centerlines, feature identification code (using linear referencing systems where available), functional class, name (including route numbers), and street address ranges<\/td>\n<\/tr>\n<tr>\n<td>Trails<\/td>\n<td>Centerlines, feature identification code (using linear referencing systems where available), name, and type<\/td>\n<\/tr>\n<tr>\n<td>Railroads<\/td>\n<td>Centerlines, feature identification code (using linear referencing systems where available), and type<\/td>\n<\/tr>\n<tr>\n<td>Waterways<\/td>\n<td>Centerlines, feature identification code (using linear referencing systems where available), and name<\/td>\n<\/tr>\n<tr>\n<td>Airports and ports<\/td>\n<td>Feature identification code and name<\/td>\n<\/tr>\n<tr>\n<td>Bridges and tunnels<\/td>\n<td>Feature identification code and name<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<h3>IDENTIFICATION<\/h3>\n<p>As part of the National Map initiative, USGS and partners are developing a comprehensive national database of vector transportation data. The transportation theme \u201cincludes best available data from Federal partners such as the Census Bureau and the Department of Transportation, State and local agencies\u201d (USGS, 2007).<\/p>\n<p>As envisioned by FGDC,\u00a0<strong>centerlines<\/strong>\u00a0are used to represent transportation routes. Like the lines painted down the middle of two-way streets, centerlines are 1-dimensional vector features that approximate the locations of roads, railroads, and navigable waterways. In this sense, road centerlines are analogous to the flowpaths encoded in the National Hydrologic Dataset (see previous page). Also like the NHD (and TIGER), road\u00a0<strong>topology<\/strong>\u00a0must be encoded to facilitate analysis of transportation networks.<\/p>\n<p>To get a sense of the complexity of the features and attributes that comprise the transportation theme, see the\u00a0<a href=\"http:\/\/services.nationalmap.gov\/bestpractices\/model\/acrodocs\/Poster_BPTrans_03_01_2006.pdf\">Transportation Data Model<\/a>(This is a 36\u2033 x 48\u2033 poster in a 5.2 Mb PDF file.) [The link to the Transportation Data Model poster recently became disconnected. Instead look at the model diagrams in the\u00a0<a title=\"FGDC Transportation Base\" href=\"http:\/\/www.fgdc.gov\/standards\/projects\/FGDC-standards-projects\/framework-data-standard\/GI_FrameworkDataStandard_Part7_Transportation_Base.pdf\">Part 7: Transportation Base<\/a>of the FGDC Geographic Framework Data Content Standard.]<\/p>\n<p>In the U.S. at least, the best road centerline data is that produced by NAVTEQ and Tele Atlas, which license transportation data to routing sites like Google Maps and MapQuest, and to manufacturers of in-car GPS navigation systems. Because these data are proprietary, however, USGS must look elsewhere for data that can be made available for public use. TIGER\/Line data produced by the Census Bureau will likely play an important role after the TIGER\/MAF Modernization project is complete (see Chapter 4).<\/p>\n<h3>DISTRIBUTION<\/h3>\n<h3><strong>TRY THIS!<\/strong><\/h3>\n<h3>VIEW AND DOWNLOAD NATIONAL MAP TRANSPORTATION DATA<\/h3>\n<ol>\n<li><a href=\"http:\/\/bpgeo.cr.usgs.gov\/\">Access the Viewer here<\/a>.<\/li>\n<li>Expand the pane containing the layer options by clicking on Overlays at the upper-left.<\/li>\n<li>Under Base Data Layers, click on Transportation. You can expand the Transportation list and sub-select different layers.<\/li>\n<li>As you zoom in to larger map scales (using the slider bar at the upper-left of the map), additional transportation layers will become visible.<\/li>\n<li>If you wish to download an extract from the transportation database, click the Download Data button in the upper-right of the viewer interface and decide how you wish to extract the data. The Transportation data comes down in ESRI\u2019s geodatabase format. Additional information regarding downloading data can be found via the Help button in the upper-right of the viewer interface.<\/li>\n<\/ol>\n<h2>7.18. Theme: Governmental Units<\/h2>\n<p>The FGDC framework also includes boundaries of governmental units, including:<\/p>\n<ul>\n<li>Nation<\/li>\n<li>States and statistically equivalent areas<\/li>\n<li>Counties and statistically equivalent areas<\/li>\n<li>Incorporated places and consolidated cities<\/li>\n<li>Functioning legal minor civil divisions<\/li>\n<li>Federal- or state-recognized American Indian reservations and trustlands<\/li>\n<li>Alaska native regional corporations<\/li>\n<\/ul>\n<p>FGDC specifies that:<\/p>\n<p>Each of these features includes the attributes of name and the applicable Federal Information Processing Standard (FIPS) code. Features boundaries include information about other features (such as road, railroads, or streams) with which the boundaries are associated and a description of the association (such as coincidence, offset, or corridor. (FGDC, 1997, p. 20-21)<\/p>\n<h3>IDENTIFICATION<\/h3>\n<p>The USGS National Map aspires to include a comprehensive database of boundary data. In addition to the entities outlined above, the National Map also lists congressional districts, school districts, and ZIP Code zones. Sources for these data include \u201cFederal partners such as the U.S. Census Bureau, other Federal agencies, and State and local agencies.\u201d (USGS, 2007).<\/p>\n<p>To get a sense of the complexity of the features and attributes that comprise this theme, see the\u00a0<a href=\"http:\/\/services.nationalmap.gov\/bestpractices\/model\/acrodocs\/Poster_BPGovtUnits_03_01_2006.pdf\">Governmental Units Data Model<\/a>\u00a0(This is a 36\u2033 x 48\u2033 poster in a 2.4 Mb PDF file.) [The link to the Governmental Units Data Model poster recently became disconnected. Instead look at the model diagrams in the\u00a0<a title=\"FGDC Governmental Unit Base\" href=\"http:\/\/www.fgdc.gov\/standards\/projects\/FGDC-standards-projects\/framework-data-standard\/GI_FrameworkDataStandard_Part5_GovernmentalUnitBoundaries.pdf\">Part 5: Governemntal unit and other geographic area boundaries<\/a>\u00a0of the FGDC Geographic Framework Data Content Standard.]<\/p>\n<h3>DISTRIBUTION<\/h3>\n<h3><strong>TRY THIS!<\/strong><\/h3>\n<h3>VIEW AND DOWNLOAD NATIONAL MAP GOVERNMENTAL UNITS DATA<\/h3>\n<ol>\n<li><a href=\"http:\/\/bpgeo.cr.usgs.gov\/\">Access the Viewer here<\/a>.<\/li>\n<li>Expand the pane containing the layer options by clicking on Overlays at the upper-left.<\/li>\n<li>Under Base Data Layers, click on Governmental Unit Boundaries. You can expand the this list and sub-select different boundary layers.<\/li>\n<li>As you zoom in to larger map scales (using the slider bar at the upper-left of the map), additional boundary layers will become visible.<\/li>\n<li>If you wish to download an extract from the Governmental Unit Boundaries database, click the Download Data button in the upper-right of the viewer interface and decide how you wish to extract the data. The Governmental Unit Boundaries data comes down in ESRI\u2019s geodatabase format. Additional information regarding downloading data can be found via the Help button in the upper-right of the viewer interface.<\/li>\n<\/ol>\n<h2>7.19. Theme: Cadastral<\/h2>\n<p>FGDC (1997, p. 21) points out that:<\/p>\n<p>Cadastral data represent the geographic extent of the past, current, and future rights and interests in real property. The spatial information necessary to describe the geographic extent and the rights and interests includes surveys, legal description reference systems, and parcel-by-parcel surveys and descriptions.<\/p>\n<p>However, no one expects that legal descriptions and survey coordinates of private property boundaries (as depicted schematically in the portion of the\u00a0<strong>plat map<\/strong>\u00a0shown below) will be included in the USGS National Map any time soon. As discussed at the outset of Chapter 6, this is because local governments have authority for land title registration in the U.S., and most of these governments have neither the incentive nor the means to incorporate such data into a publicly-accessible national database.<\/p>\n<p><img decoding=\"async\" alt=\"Portion of a plat map showing property boundaries\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/plat.png\" \/><\/p>\n<p>Plat maps are supplementary records that depict property parcel boundaries in graphic form. The geometric accuracy of plats is notoriously poor. The investment required to convert plat maps to properly georeferenced digital data is substantial. Many local governments have converted these records to digital form, or are in the process of doing so.<\/p>\n<p>FGDC\u2019s modest goal for the cadastal theme of the NSDI framework is to include:<\/p>\n<p>\u2026cadastral reference systems, such as the Public Land Survey System (PLSS) and similar systems not covered by the PLSS \u2026 and publicly administered parcels, such as military reservations, national forests, and state parks. (Ibid, p. 21)<\/p>\n<p><a href=\"http:\/\/www.fgdc.gov\/standards\/standards_publications\/\">FGDC\u2019s Cadastral Data Content Standard is published here<\/a>.<\/p>\n<p>The colored areas on the map below show the extent of the United States Public Land Surveys, which commenced in 1784 and took nearly a century to complete (Muehrcke and Muehrcke, 1998). The purpose of the surveys was to partition \u201cpublic land\u201d into saleable parcels in order to raise revenues needed to retire war debt, and to promote settlement. A key feature of the system is its nomenclature, which provides concise, unique specifications of the location and extent of any parcel.<\/p>\n<p><img decoding=\"async\" alt=\"Map of United States Public Land Survey system\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/USPLS.jpg\" \/><\/p>\n<p>Extent of the U.S. Public Land Survey (Thompson, 1988).<\/p>\n<p>Each Public Land Survey (shown in the colored areas above) commenced from an\u00a0<strong>initial point\u00a0<\/strong>at the precisely surveyed intersection of a\u00a0<strong>base line<\/strong>\u00a0and\u00a0<strong>principal meridian<\/strong>. Surveyed lands were then partitioned into grids of\u00a0<strong>townships\u00a0<\/strong>each approximately six miles square.<\/p>\n<p><img decoding=\"async\" alt=\"U.S. Public Land Survey Township grid system\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/uspls_township_range.png\" \/><\/p>\n<p><strong>Townships\u00a0<\/strong>are designated by their locations relative to the base line and principal meridian of a particular survey. For example, the township highlighted in gold above is the second township south of the baseline and the third township west of the principal meridian. The Public Land Survey designation for the highlighted township is \u201cTownship 2 South, Range 3 West.\u201d Because of this nomenclature, the Public Land Survey System is also known as the \u201ctownship and range system.\u201d\u00a0 Township T2S, R3W is shown enlarged below.<\/p>\n<p><img decoding=\"async\" alt=\"U.S. Public Land Survey township grid\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/uspls_township.png\" \/><\/p>\n<p>Townships are subdivided into grids of 36\u00a0<strong>sections<\/strong>. Each section covers approximately\u00a0<strong>one square mile<\/strong>\u00a0(640 acres). Notice the back-and-forth numbering scheme. Section 14, highlighted in gold above, is shown enlarged below.<\/p>\n<p><img decoding=\"async\" alt=\"U.S. Public Land Survey section showing property designations\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/uspls_section.png\" \/><\/p>\n<p>Inidividual property parcels are designated as shown above. For instance, the\u00a0<strong>NE 1\/4 of Section 14, Township 2 S, Range 3W<\/strong>, is a 160-acre parcel. Public Land Survey designations specify both the location of a parcel and its area.<\/p>\n<p><img decoding=\"async\" alt=\"Portion of topographic map showning influence of Public Land Survey on road network in midwest U.S.\" src=\"http:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-content\/uploads\/sites\/19\/2014\/01\/uspls_topo.png\" \/><\/p>\n<p>The influence of the Public Land Survey grid is evident in the built environment of much of the American Midwest. As Mark Monmonier (1995, p. 114) observes:<\/p>\n<p>The result [of the U.S. Public Land Survey] was an \u2018authored landscape\u2019 in which the survey grid had a marked effect on settlement patterns and the shapes of counties and smaller political units. In the typical Midwestern county, roads commonly following section lines, the rural population is dispersed rather than clustered, and the landscape has a pronounced checkerboard appearance.<\/p>\n<p>For more information about the\u00a0<a href=\"http:\/\/nationalatlas.gov\/articles\/boundaries\/a_plss.html\">Public Land Survey System, see this article in the in the USGS\u2019 National Atlas<\/a>.<\/p>\n<h2>7.20. Summary<\/h2>\n<p>NSDI framework data represent \u201cthe most common data themes [that] users need\u201d (FGDC, 1997, p. 3), including geodetic control, orthoimagery, elevation, hydrography, transportation, governmental unit boundaries, and cadastral reference information. Some themes, like transportation and governmental units, represent things that have well-defined edges. In this sense we can think of things like roads and political boundaries as\u00a0<strong>discrete\u00a0<\/strong>phenomena. The\u00a0<strong>vector\u00a0<\/strong>approach to geographic representation is well suited to digitizing discrete phenomena. Line features do a good job of representing roads, for example, and polygons are useful approximations of boundaries.<\/p>\n<p>As you recall from Chapter 1, however, one of the distinguishing properties of the Earth\u2019s surface is that it is\u00a0<strong>continuous<\/strong>. Some phenomena distributed across the surface are continuous too. Terrain elevations, gravity, magnetic declination and surface air temperature can be measured practically everywhere. For many purposes,\u00a0<strong>raster<\/strong>data are best suited to representing continuous phenomena.<\/p>\n<p>An implication of continuity is that there is an infinite number of locations at which phenomena can be measured. It is not possible, obviously, to take an infinite number of measurements. Even if it were, the mass of data produced would not be usable. The solution, of course, is to collect a sample of measurements, and to estimate attribute values for locations that are left unmeasured. Chapter 7 also considers how missing elevations in a raster grid can be estimated from existing elevations, using a procedure called\u00a0<strong>interpolation<\/strong>. The inverse distance weighted interpolation procedure relies upon another fundamental property of geographic data,\u00a0<strong>spatial dependence<\/strong>.<\/p>\n<p>The chapter concludes by investigating the characteristics and current status of the hydrography, transportation, governmental units, and cadastral themes. You had the opportunity to access, download, and open several of the data themes using viewers provided by USGS as part of its National Map initiative. In general, you should have found that although neither the NSDI or National Map visions have been fully realized, substantial elements of each is in place. Further progress depends on the American public\u2019s continuing commitment to public data, and to the political will of our representatives in government.<\/p>\n<h3><strong>QUIZ<\/strong><\/h3>\n<p>Registered Penn State students should return now to the Chapter 7 folder\u00a0 in ANGEL (via the Resources menu to the left) to access the graded quiz for this chapter. This one counts.\u00a0<strong>You may take graded quizzes only once.<\/strong><\/p>\n<p>The purpose of the quiz is to ensure that you have studied the text closely, that you have mastered the practice activities, and that you have fulfilled the chapter\u2019s learning objectives. You are welcome to review the chapter during the quiz.<\/p>\n<p>Once you have submitted the quiz and posted any questions you may have to either our discussion forums or chapter pages, you will have completed Chapter 7.<\/p>\n<h3>COMMENTS AND QUESTIONS<\/h3>\n<p>Registered students are welcome to post comments, questions, and replies to questions about the text. Particularly welcome are anecdotes that relate the chapter text to your personal or professional experience. In addition, there are discussion forums available in the ANGEL course management system for comments and questions about topics that you may not wish to share with the whole world.<\/p>\n<p>To post a comment, scroll down to the text box under \u201cPost new comment\u201d and begin typing in the text box, or you can choose to reply to an existing thread. When you are finished typing, click on either the \u201cPreview\u201d or \u201cSave\u201d button (Save will actually submit your comment). Once your comment is posted, you will be able to edit or delete it as needed. In addition, you will be able to reply to other posts at any time.<\/p>\n<p>Note: the first few words of each comment become its \u201ctitle\u201d in the thread.<\/p>\n<h2>7.21. Bibliography<\/h2>\n<p>Federal Geographic Data Committee (1997).\u00a0<em>Framework introduction and guide.\u00a0<\/em>Washington DC: Federal Geographic Data Committee.<\/p>\n<p>Eischeid, J. D., Baker, C. B., Karl, R. R., Diaz, H. F. (1995). The quality control of long-term climatological data using objective data analysis.\u00a0<em>Journal of Applied Meteorology, 34<\/em>, 27-88.<\/p>\n<p>Gould, P. (1989). Geographic dimensions of the AIDS epidemic.<em>Professional Geographer, 41:1<\/em>, 71-77.<\/p>\n<p>Maune, D. F. (Ed.) (2007).\u00a0<em>Digital elevation model technologies and applications: The DEM users manual, 2nd edition.\u00a0<\/em>Bethesda, MD: American Society for Photogrammetric Engineering and Remote Sensing.<\/p>\n<p>Monmonier, M. S. (1982).\u00a0<em>Drawing the line: tales of maps and cartocontroversy.<\/em>\u00a0New York, NY: Henry Holt.<\/p>\n<p>Muehrcke, P. C. and Muehrcke, J. O. (1998)\u00a0<em>Map use,<\/em>\u00a04th Ed. Madison, WI: JP Publications.<\/p>\n<p>National Aeronautics and Space Administration, Jet Propulsion Laboratory (2006).\u00a0<em>Shuttle radar topography mission<\/em>. Retrieved May 10, 2006, from\u00a0<a href=\"http:\/\/www.jpl.nasa.gov\/srtm\">http:\/\/www.jpl.nasa.gov\/srtm<\/a><\/p>\n<p>Goddard Space Flight Center, National Aeronautics and Space Administration (n.d.). Greenland\u2019s receding ice. Retrieved Feburary 26, 2008, from\u00a0<a href=\"http:\/\/svs.gsfc.nasa.gov\/stories\/greenland\/\">http:\/\/svs.gsfc.nasa.gov\/stories\/greenland\/<\/a><\/p>\n<p>National Geophysical Data Center (2010). ETOPO1 global gridded 1 arc-minute database. Retrieved March 2, 2010, from<a href=\"http:\/\/www.ngdc.noaa.gov\/mgg\/global\/global.html\">http:\/\/www.ngdc.noaa.gov\/mgg\/global\/global.html<\/a><\/p>\n<p>National Oceanic and Atmospheric Administration, National Climatic Data Center (n. d.).\u00a0<em>Merged land-ocean seasonal temperature anomalies<\/em>. Retrieved August 18, 1999, from<a title=\"http:\/\/www.ncdc.noaa.giv\/onlineprod\/landocean\/seasonal\/form.html\" href=\"http:\/\/www.ncdc.noaa.giv\/onlineprod\/landocean\/seasonal\/form.html\">http:\/\/www.ncdc.noaa.giv\/onlineprod\/landocean\/seasonal\/form.html<\/a>(expired)<\/p>\n<p>National Oceanic and Atmospheric Administration (2002). Side scan and multibeam sonar. Retrieved February 18, 2008, from<a title=\"http:\/\/www.nauticalcharts.noaa.gov\/hsd\/hydrog.htm\" href=\"http:\/\/www.nauticalcharts.noaa.gov\/hsd\/hydrog.htm\">http:\/\/www.nauticalcharts.noaa.gov\/hsd\/hydrog.htm<\/a><\/p>\n<p>National Oceanic and Atmospheric Administration (2007). NOAA History. Retrieved February 27, 2008, from<a href=\"http:\/\/www.history.noaa.gov\/\">http:\/\/www.history.noaa.gov\/<\/a><\/p>\n<p>Rabenhorst, T. D. and McDermott, P. D. (1989).\u00a0<em>Applied cartography: source materials for mapmaking.\u00a0<\/em>Columbus, OH: Merrill.<\/p>\n<p>Raitz, E. (1948).\u00a0<em>General cartography.<\/em>\u00a0New York, NY: McGraw-Hill.<\/p>\n<p>Ralston, B. A. (2004).\u00a0<em>GIS and public data.<\/em>\u00a0Clifton Park NY: Delmar Learning.<\/p>\n<p>Thompson, M. M. (1988) Maps for america, 3rd Ed. Reston, VA: United States Geological Survey.<\/p>\n<p>United States Geological Survey (1987) Digital elevation models. Data users guide 5. Reston, VA: USGS.<\/p>\n<p>United States Geological Survey (1999) The National Hydrography Dataset. Fact Sheet 106-99. Reston, VA: USGS. Retrieved February 19, 2008 from\u00a0<a href=\"http:\/\/erg.usgs.gov\/isb\/pubs\/factsheets\/fs10699.html\">http:\/\/erg.usgs.gov\/isb\/pubs\/factsheets\/fs10699.html<\/a><\/p>\n<p>United States Geological Survey (2000) The National Hydrography Dataset: Concepts and Contents. Reston, VA: USGS. Retrieved February 19, 2008 from<a href=\"http:\/\/nhd.usgs.gov\/chapter1\/chp1_data_users_guide.pdf\">http:\/\/nhd.usgs.gov\/chapter1\/chp1_data_users_guide.pdf<\/a><\/p>\n<p>United States Geological Survey (2002) The National Map \u2013 Hydrography. Fact Sheet 060-02. Reston, VA: USGS. Retrieved February 19, 2008 from<a href=\"http:\/\/erg.usgs.gov\/isb\/pubs\/factsheets\/fs06002.html\">http:\/\/erg.usgs.gov\/isb\/pubs\/factsheets\/fs06002.html<\/a>\u00a0Retrieved September 22, 2013 from\u00a0<a href=\"http:\/\/pubs.er.usgs.gov\/publication\/fs06002\">http:\/\/pubs.er.usgs.gov\/publication\/fs06002<\/a><\/p>\n<p>United States Geological Survey (2006a) Digital Line Graphs (DLG). Reston, VA: USGS. Retrieved February 18, 2008 from<a title=\"http:\/\/edc.usgs.gov\/products\/map\/dlg.html\" href=\"http:\/\/edc.usgs.gov\/products\/map\/dlg.html\">http:\/\/edc.usgs.gov\/products\/map\/dlg.html<\/a>\u00a0(In 2010 the site became<a href=\"http:\/\/eros.usgs.gov\/#\/Find_Data\/Products_and_Data_Available\/DLGs\">http:\/\/eros.usgs.gov\/#\/Find_Data\/Products_and_Data_Available\/DLGs<\/a>)<\/p>\n<p>United States Geological Survey (2006b) GTOPO30. Retrieved February 27, 2008 from<a title=\"http:\/\/edc.usgs.gov\/products\/elevation\/gtopo30\/gtopo30.html\" href=\"http:\/\/edc.usgs.gov\/products\/elevation\/gtopo30\/gtopo30.html\">http:\/\/edc.usgs.gov\/products\/elevation\/gtopo30\/gtopo30.html<\/a>\u00a0since moved to\u00a0<a href=\"http:\/\/www1.gsi.go.jp\/geowww\/globalmap-gsi\/gtopo30\/gtopo30.html\">http:\/\/www1.gsi.go.jp\/geowww\/globalmap-gsi\/gtopo30\/gtopo30.html<\/a><\/p>\n<p>United States Geological Survey (2006c) National Hydrography Dataset (NHD) \u2013 High-resolution (Metadata). Reston, VA: USGS. Retrieved February 19, 2008 from<a href=\"http:\/\/nhdgeo.usgs.gov\/metadata\/nhd_high.htm\">http:\/\/nhdgeo.usgs.gov\/metadata\/nhd_high.htm<\/a><\/p>\n<p>United States Geological Survey (2007). Vector data theme development of<em>\u00a0The National Map.\u00a0<\/em>Retrieved 24 February 2008 from<a title=\"http:\/\/bpgeo.cr.usgs.gov\/model\/\" href=\"http:\/\/bpgeo.cr.usgs.gov\/model\/\">http:\/\/bpgeo.cr.usgs.gov\/model\/<\/a>\u00a0(expired or moved)<\/p>\n<p>United States Geological Survey (2009) The National Map \u2013 Hydrography Dataset. Reston, VA: USGS. Retrieved September 22, 2013 from\u00a0<a href=\"http:\/\/pubs.usgs.gov\/fs\/2009\/3054\/pdf\/FS2009-3054.pdf\">http:\/\/pubs.usgs.gov\/fs\/2009\/3054\/pdf\/FS2009-3054.pdf<\/a><\/p>\n<p>United States Geological Survey (2013) National Hydrography Dataset (NHD) \u2013 Get NDH Data. Reston, VA: USGS. Retrieved September 22, 2013 from\u00a0<a href=\"http:\/\/nhd.usgs.gov\/data.html\">http:\/\/nhd.usgs.gov\/data.html<\/a><\/p>\n","protected":false},"author":1,"menu_order":1,"template":"","meta":{"pb_show_title":"on","pb_short_title":"","pb_subtitle":"","pb_authors":["david-dibiase"],"pb_section_license":""},"chapter-type":[],"contributor":[47],"license":[],"class_list":["post-499","chapter","type-chapter","status-publish","hentry","contributor-david-dibiase"],"part":85,"_links":{"self":[{"href":"https:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-json\/pressbooks\/v2\/chapters\/499","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-json\/pressbooks\/v2\/chapters"}],"about":[{"href":"https:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-json\/wp\/v2\/types\/chapter"}],"author":[{"embeddable":true,"href":"https:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-json\/wp\/v2\/users\/1"}],"version-history":[{"count":4,"href":"https:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-json\/pressbooks\/v2\/chapters\/499\/revisions"}],"predecessor-version":[{"id":886,"href":"https:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-json\/pressbooks\/v2\/chapters\/499\/revisions\/886"}],"part":[{"href":"https:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-json\/pressbooks\/v2\/parts\/85"}],"metadata":[{"href":"https:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-json\/pressbooks\/v2\/chapters\/499\/metadata\/"}],"wp:attachment":[{"href":"https:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-json\/wp\/v2\/media?parent=499"}],"wp:term":[{"taxonomy":"chapter-type","embeddable":true,"href":"https:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-json\/pressbooks\/v2\/chapter-type?post=499"},{"taxonomy":"contributor","embeddable":true,"href":"https:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-json\/wp\/v2\/contributor?post=499"},{"taxonomy":"license","embeddable":true,"href":"https:\/\/opentextbc.ca\/natureofgeographicinformation\/wp-json\/wp\/v2\/license?post=499"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}