Generalization for Generalization for The The National Map National Map with emphasis on with emphasis on the National Hydrography the National Hydrography Dataset Dataset Lawrence Stanislawski, Science Applications International Lawrence Stanislawski, Science Applications International Corporation (SAIC) Corporation (SAIC) Michael Starbuck, U.S. Geological Survey Michael Starbuck, U.S. Geological Survey Michael Finn, U.S. Geological Survey Michael Finn, U.S. Geological Survey E. Lynn Usery, U.S. Geological Survey E. Lynn Usery, U.S. Geological Survey Pat Turley, U.S. Geological Survey Pat Turley, U.S. Geological Survey Josh Wilkerson, U.S. Geological Survey Josh Wilkerson, U.S. Geological Survey 2005 ESRI International User Conference, July 25-29. San Diego, CA.
54
Embed
Generalization for The National Map with emphasis on the National Hydrography Dataset Lawrence Stanislawski, Science Applications International Corporation.
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Generalization for Generalization for The National Map The National Map with emphasis onwith emphasis on
the National Hydrography Datasetthe National Hydrography Dataset
Lawrence Stanislawski, Science Applications International Corporation (SAIC)Lawrence Stanislawski, Science Applications International Corporation (SAIC)Michael Starbuck, U.S. Geological SurveyMichael Starbuck, U.S. Geological SurveyMichael Finn, U.S. Geological SurveyMichael Finn, U.S. Geological SurveyE. Lynn Usery, U.S. Geological SurveyE. Lynn Usery, U.S. Geological SurveyPat Turley, U.S. Geological SurveyPat Turley, U.S. Geological SurveyJosh Wilkerson, U.S. Geological SurveyJosh Wilkerson, U.S. Geological Survey
2005 ESRI International User Conference, July 25-29. San Diego, CA.
Overview of The National MapOverview of The National Map
Geometric characteristics of hydrographic networkGeometric characteristics of hydrographic network
Test DataTest Data
ResultsResults
2005 ESRI International User Conference, July 25-29. San Diego, CA.
The National MapThe National Map
•Consistent framework for geographic knowledge needed by the Consistent framework for geographic knowledge needed by the Nation Nation •Provides public access to high-quality, geospatial data andProvides public access to high-quality, geospatial data andinformation from multiple partners to help support decision-information from multiple partners to help support decision-making by resource managers and the public. making by resource managers and the public. •Product of a consortium of Federal, State, and local partnersProduct of a consortium of Federal, State, and local partnerswho provide geospatial data to enhance America's ability towho provide geospatial data to enhance America's ability toaccess, integrate, and apply geospatial data at global, national,access, integrate, and apply geospatial data at global, national,and local scales. and local scales. •Includes eight primary data themes:Includes eight primary data themes:
Orthoimagery,Orthoimagery,Elevation,Elevation,Hydrography,Hydrography,Transportation,Transportation,Structures,Structures,Boundaries, Boundaries, Geographic names, Geographic names, Land use/land cover, Land use/land cover,
•Most themes are populated, and content will continue to beMost themes are populated, and content will continue to beadded and updatedadded and updated
2005 ESRI International User Conference, July 25-29. San Diego, CA.
The National MapThe National Map
•Created from many data sources and Web Map Services Created from many data sources and Web Map Services (WMS) (WMS)
owned by many organizations.owned by many organizations.•Data sources are electronically stored at various scales and Data sources are electronically stored at various scales and resolutions. Consequently, when data are extracted from several resolutions. Consequently, when data are extracted from several WMS and horizontally integrated in geospatial data applications, WMS and horizontally integrated in geospatial data applications, the differences between data sources become apparent. the differences between data sources become apparent. •The ability to render one or more scales of The ability to render one or more scales of The National MapThe National Map data into an appropriate, or functionally equivalent, data into an appropriate, or functionally equivalent,
representation at a user-specified scale could substantially representation at a user-specified scale could substantially enhance analysis capabilities of these data.enhance analysis capabilities of these data.
Research Topic: development of automated generalization approach that renders a functionally equivalent dataset at a user-specified scale, focusing on the National Hydrography Dataset (NHD) layer of The National Map.
2005 ESRI International User Conference, July 25-29. San Diego, CA.
The National MapThe National MapNational Hydrography Dataset (NHD)National Hydrography Dataset (NHD)
2005 ESRI International User Conference, July 25-29. San Diego, CA.
National Hydrography Dataset
• Vector data layer of The National Map representing surface waters of the United States.
• Includes a set of surface water reaches
Reach: significant segment of surface water having similar hydrologic characteristics, such as a stretch of river between two confluences, a lake, or a pond.
A unique address, called a reach code, is assigned to each reach, which enables linking of ancillary data to specific features and locations on the NHD.
Reach code from Lower Mississippi subbasin
08010100000413
region-subregion-accounting unit-subbasin-reach number
08010100000696
08010100000413
The National MapThe National MapNational Hydrography Dataset (NHD)National Hydrography Dataset (NHD)
2005 ESRI International User Conference, July 25-29. San Diego, CA.
National Hydrography Dataset
• Divided and distributed at watershed basin and subbasin boundaries.
• Stored in an ArcGIS geographic database (geodatabase) model.
• Three levels of detail (resolutions)
Medium (1:100,000-scale source)
only complete layer
High (1:24,000-scale source)
three-quarters complete
Local (1:12,000 or larger source)Subregions along northern shore of Gulf of Mexico
ObjectivesObjectives
2005 ESRI International User Conference, July 25-29. San Diego, CA.
Develop a generalization strategy that can be implemented on subsets of the NHD.
Strategy should produce a dataset in the NHD model format that maintains:• feature definitions, • reach delineations, • feature relationships, and • flow connections between remaining generalized features.
Extracted dataset should function with NHD applications • less detail • faster processing speed
Extracted level of detail:• user-specified, or • defined by display scale
Generalization ApproachGeneralization Approach
2005 ESRI International User Conference, July 25-29. San Diego, CA.
Base data: user specified, or highest resolution that covers desired area.
Feature pruning – removal of feature that are too small for desired output scale.• select a subset of network features• select a subset of area features• remove point features associated with pruned line or area features
Feature simplification• removal of vertices• aggregation, amalgamation, merging, etc.
Focus so far has been on network pruning
National Hydrography Dataset (NHD)National Hydrography Dataset (NHD)Flowline NetworkFlowline Network
2005 ESRI International User Conference, July 25-29. San Diego, CA.
Drainage network: flow-oriented digital lines that represent the connectivity and flow between surface water features.
Network features types:• Stream/river• Canal/ditch• Pipeline• Connector• Artificial path
Connector represents a path where surface flow is known to exist but was not included in the source material.
Artificial path represents a flow path through an areal water feature that is connected to the drainage network. Areal features that may include artificial paths are:
area of complex channel, canal/ditch, estuary, ice mass, lake/pond, playa, reservoir, stream/river, swamp/marsh, and wash
NHD FeaturesNHD Features
Areal Stream/River
Areal Lake/Pond
Areal Lake/Pond
Linear Streams
Linear Canal/Ditch
Linear Connector
Artificial Paths
National Hydrography Dataset (NHD)National Hydrography Dataset (NHD)Flowline NetworkFlowline Network
2005 ESRI International User Conference, July 25-29. San Diego, CA.
2005 ESRI International User Conference, July 25-29. San Diego, CA.
Quantities used to compare geometric characteristics of a stream drainage network:
Catchment area – surface runoff in this area flows into the associated network segment.Upstream arbolate sum – sum of the lengths of all upstream network features that flow
into a point of a network.Drainage proportion – proportion of the arbolate sum at a point on the network to the
total length of the drainage network.Upstream drainage area – sum of all upstream catchment areas up to a point in the
network.Upstream drainage density – ratio of the arbolate sum to the upstream drainage area at
a point in the network.Basin drainage density – the ratio of length of all network features to the sum of all
catchment areas in the network. The sum of all catchment areas in the network may also be referred to as the basin area.
Number of confluence-to-confluence segments – the number of segments between network feature confluences (intersections) that exist in the network.
Stream order – method of numbering network segments within a drainage system, in which the smallest unbranched mapped tributary is called first order. Segments where two or more first order tributaries intersect becomes a second order segment and so on (Horton 1945).
Bifurcation ratio – the ratio of the number of segments of any order to the number of segments of the next highest order.
Basin bifurcation ratio – the average of the bifurcation ratios of each stream order in the drainage network.
MethodsMethods
2005 ESRI International User Conference, July 25-29. San Diego, CA.
Network drainage characteristics summarized by NHD subbasin
Arc Macro Language Programs
Python geoprocessing scripts
Rapid method that uses Thiessen polygons to the estimate catchment area for each network drainage feature was developed for these computations
MethodsMethods
2005 ESRI International User Conference, July 25-29. San Diego, CA.
Flowline Network Pruning
• Extract the most prominent network features.
• Eliminate network features having an upstream
drainage area that is less than a
threshold. Features that drain larger areas are more likely to be
maintained than features that drain smaller areas. Minimizes the effect of disparate data collection. Produces a uniform distribution of network features.
• Iteratively increase upstream drainage
area threshold until desired output drainage
density is achieved.Requires scale-dependency for ease of use.
MethodsMethodsFlowline Network Pruning
2005 ESRI International User Conference, July 25-29. San Diego, CA.
Variation in density of mapped drainage features along 7.5-minute quadrangle boundaries for medium resolution NHD. Background colors delineate subbasins.
2005 ESRI International User Conference, July 25-29. San Diego, CA.
Modeled relationship between • Elevation-derived drainage density• Minimum drainage area for stream formation
Elevation data extracted from the National Elevation Dataset (NED, USGS 2005) and resampled to a 30 meter cell size.
Approximate stream channels were derived from the elevation data through commonly used methods.
Varied minimum drainage area for stream formation on the flow accumulation grid between 0.135 km2 (150 cells) to 90 km2
(100,000 cells)
Completed for two study areas.
MethodsMethods
2005 ESRI International User Conference, July 25-29. San Diego, CA.
Gasconade-Osage subregion (1029) falls in the Interior Plains and Interior Highlands physiographic divisions.
MethodsMethods
2005 ESRI International User Conference, July 25-29. San Diego, CA.
Colorado Headwaters and Gunnison subregions (1401 and 1402) fall in Intermontane Plateaus and Rocky Mountain System physiographic divisions.
MethodsMethods
2005 ESRI International User Conference, July 25-29. San Diego, CA.
Determined that an approximate linear relationship should exist between the square root of map scale and drainage density derived from mapped network features.
Approximate drainage densities determined for the two study areas at four scales of hydrographic vector data:
2005 ESRI International User Conference, July 25-29. San Diego, CA.
Sample Dataset from Gasconade-Osage subregion:•Pomme De Terre (PDT) subbasin (NHD 10290107).•Medium level of drainage density difference between the high and medium resolution NHD subbasins for the PDT (difference about 0.6421 km/km2).
Relative magnitude of drainage density differences between the high resolution layer and the medium resolution layer of the region. Darker red indicates a higher difference in density.
2005 ESRI International User Conference, July 25-29. San Diego, CA.
Pomme De Terre network generalization
Graduated symbols of arbolate sum values by network reach for Pomme De Terre subbasin. High resolution (green), Medium resolution (purple), Generalized (pink).
2005 ESRI International User Conference, July 25-29. San Diego, CA.
Pomme De Terre network generalization
0
50
100
150
200
250
300
350
400
450
500
0 5 10 15 20
Catchment Drainage Density (km/sq km)
Fre
qu
en
cy
Drainage density distribution of reach catchments for Pomme De Terre subbasin of high-resolution NHD network generalized to 0.942 km/km2 basin drainage density. Maximum frequency at 1.32 km/km2
2005 ESRI International User Conference, July 25-29. San Diego, CA.
Pomme De Terre network generalization
Catchment area and upstream drainage density relationships with arbolate sum for high-resolution Pomme De Terre subbasin generalized to 0.942 km/sq km basin drainage density
2005 ESRI International User Conference, July 25-29. San Diego, CA.
Pomme De Terre network generalization
Summary of stream frequency by stream order for high resolution, medium resolution, and generalized high resolution Pomme De Terre subbasin.
High Resolution
0
500
1000
1500
2000
0 3 6 9 12 15
Stream Order
Fre
qu
en
cy
Generalized High Resolution
0
500
1000
1500
2000
0 3 6 9 12 15
Stream Order
Fre
qu
en
cy
Medium Resolution
0
500
1000
1500
2000
0 3 6 9 12 15
Stream Order
Fre
qu
en
cy
SummarySummary
2005 ESRI International User Conference, July 25-29. San Diego, CA.
Identified and automated computations for quantities that characterize an NHD flowline network.
Reach and basin statistics, along with visual inspections assist analysis and comparison of the network drainage.
Identified a linear relationship to estimate watershed drainage density of elevation-derived streams from the minimum upstream drainage area for stream formation.
Identifed a linear relationship to estimate drainage density appropriate for a user-specified map scales, to be used for network pruning.
Developed and tested a network pruning strategy that provides a user-specified network drainage density based on upstream drainage area, which also maintains reaches and network connectivity.
Network pruning process furnishes a network that is structurally similar to a target dataset. However, when compared to the target dataset, a larger number of shorter, confluence-to-confluence segments and reaches were maintained in the generalized network because that was a pre-existing condition of the source dataset.
SummarySummary
2005 ESRI International User Conference, July 25-29. San Diego, CA.
Upstream drainage area quickly rises to overwhelming sizes. We may need to implement a numeric scaling technique to the pruning strategy.
Additional work is required to tailor the pruning strategy to be affected by terrain or topographic variations.
Next steps in developing an automated generalization include defining a pruning strategy for area features and rules for feature simplification, and also implementing a line simplification process.
Aside from enhancing our understanding of the relationship between terrain and surface hydrography, it is possible that similar automated pruning strategies can be developed for other vector layers of The National Map.
ReferencesReferences
2005 ESRI International User Conference, July 25-29. San Diego, CA.
Gold, C. M. and D. Thibault. 2001. Map Generalization by Skeleton Retraction, In Proceedings, International Cartographic Association, Beijing,China, pp.2072-2081.
Horn, R.C., L. McKay, and, S.A. Hanson. 1994. " History of the U.S. EPA's River Reach File: A National Hydrographic Database Available for ArcInfo Applications ," in Proceedings of the Fourteenth Annual ESRI User Conference, Environmental Systems Research Institute, Redlands, CA.
Horton, R.E. 1945. Erosional development of streams and their drainage basins: Hydrophysical approach to quantitative morphology. Geological Society of America Bulletin, 56, 275-370.
Maidment, D. R. 1996. GIS and hydrologic modeling – an assessment of progress. Third International Conference on GIS and Environmental Modeling. January 22-26, 1996. Sante Fe, New Mexico.
Northcott, W. 2002. Lecture notes to Biosystems Engineering 481: Agricultural and Small Watershed Hydrology. College of Engineering, Michigan State University.
Pidwirny, M. 2004. Stream Morphometry. Accessed June 9, 2005, at URL: http://www.physicalgeography.net/fundamentals/10ab.html
Sanchez, R.R. 2002. GIS-Based Upland Erosion Modeling, Geovisualization, and Grid Size Effects on Erosion Simulations with CASC2D-SED. Ph.D. Dissertation, Civil Engineering Dept., Colorado State University.
2005 ESRI International User Conference, July 25-29. San Diego, CA.
Smith, P.N.H, and D.R. Maidment. 1995. Hydrologic Data Development System. Report 95-1 Center for Research in Water Resources. University of Texas at Austin.
Tarbotom, D. G. and D.P. Ames. 2001. Advances in the mapping of flow networks from digital elevation data. World Water and Environmental Resources Congress, May 20-24,2001, Orlando, Florida.
U.S. Geological Survey. 2005. National Elevation Dataset. Accessed June 7, 2005, at URL: http://ned.usgs.gov/.
U.S. Geological Survey. 2004. Physiographic divisions of conterminous U.S. Accessed June 14, 2005, at URL: http://aa179.cr.usgs.gov/metadata/wrdmeta/physio.htm.
U.S. Geological Survey. 2003. Streams and Waterbodies of the United States. Accessed June 14, 2005, at URL: http://nationalatlas.gov/atlasftp.html.
U.S. Geological Survey. 2000. The National Hydrography Dataset: Concepts and Contents. Accessed June 8, 2005, at URL http://nhd.usgs.gov/chapter1/index.html.
U.S. Geological Survey. 1955. Map Publication Scales. Book 1, Part B, Chapter 1 of the Geological Survey Topographic Instructions. March 1955. 13pp.
Verdin, K. L. 1997. A system for topologically coding global drainage basins and stream networks. 1997 ESRI International GIS User Conference Proceedings.
Generalization for Generalization for The National Map The National Map wwith emphasis onith emphasis onthe National Hydrography Datasetthe National Hydrography Dataset
Lawrence Stanislawski, Lawrence Stanislawski, [email protected]@usgs.govMichael Starbuck, U.S. Geological SurveyMichael Starbuck, U.S. Geological SurveyMichael Finn, U.S. Geological SurveyMichael Finn, U.S. Geological SurveyE. Lynn Usery, U.S. Geological SurveyE. Lynn Usery, U.S. Geological SurveyPat Turley, U.S. Geological SurveyPat Turley, U.S. Geological SurveyJosh Wilkerson, U.S. Geological SurveyJosh Wilkerson, U.S. Geological Survey
2005 ESRI International User Conference, July 25-29. San Diego, CA.
Pruning test on Gasconade-Osage subregionPruning test on Gasconade-Osage subregion
Green: 1:100,000Green: 1:100,000
Blue: 1:500,000Blue: 1:500,000
Red: 1:2,000,000Red: 1:2,000,000
2005 ESRI International User Conference, July 25-29. San Diego, CA.
Difference in Network Difference in Network DistributionDistribution
Green: original high Green: original high resolution resolution
data.data.
Blue: pruned data Blue: pruned data
(1:100,000)(1:100,000)
2005 ESRI International User Conference, July 25-29. San Diego, CA.
Difference in Difference in Network DistributionNetwork Distribution
Green: original high Green: original high resolution resolution
data.data.Blue: pruned data Blue: pruned data
(1:500,000) (1:500,000)
2005 ESRI International User Conference, July 25-29. San Diego, CA.