These data sets report elevation in meters above the WGS84 ellipsoid. The fill value for pixels containing no data is -9999.
High Mountain Asia 8-meter DEMs Derived from Along-track Optical Imagery, Version 1
This data set contains 8-meter Digital Elevation Models (DEMs) of high mountain Asia glacier and snow regions generated from very-high-resolution commercial stereoscopic satellite imagery.
This is the most recent version of these data.
Initial release
STANDARD Level of Service
Data: Data integrity and usability verified
Documentation: Key metadata and user guide available
User Support: Assistance with data access and usage; guidance on use of data in tools
|
Geographic Coverage |
Other Access Options
Close
Once you have logged in, you will be able to click and download files via a Web browser. There are also options for downloading via a command line or client. For more detailed instructions, please see Options Available for Bulk Downloading Data from HTTPS with Earthdata Login. Earthdata Search: This application allows you to search, visualize, and access data across thousands of Earth science data sets. Additional customization services are available for select data sets, including subsetting, reformatting, and reprojection. |
As a condition of using these data, you must cite the use of this data set using the following citation. For more information, see our Use and Copyright Web page.
Shean, D. 2017. High Mountain Asia 8-meter DEMs Derived from Along-track Optical Imagery, Version 1. [Indicate subset used]. Boulder, Colorado USA. NASA National Snow and Ice Data Center Distributed Active Archive Center. doi: https://doi.org/10.5067/GSACB044M4PK. [Date Accessed].Detailed Data Description

Sample 8 m DEM strip generated from along-track stereo imagery. Produced by D. Shean for the NASA High Mountain Asia project.
The Digital Elevation Models (DEMs) provided in these data sets were generated from very-high-resolution (VHR) along- and cross-track stereoscopic imagery from DigitalGlobe satellites.
The DEM strips have unique timestamps, allowing users to perform change detection analysis and compare observations of topographic data acquired in different seasons and/or years.
DEM strips range from approximately 13.1 km to 17.4 km wide by 111 km long. Figure 1 shows a sample shaded relief map.
DEMs are provided as LZW-compressed* GeoTIFFs (.tif) with elevations stored as single precision, floating point numbers (Float32).
Shaded relief browse images in PNG (.png) format accompany each DEM file. These images were created using a standard azimuth setting of 315 degrees and an elevation of 45 degrees.
An associated Extensible Markup Language (XML) metadata file is also provided for each data file.
*Refers to lossless Lempel–Ziv–Welch (LZW) compression.
Along-track (AT) and cross-track (CT) data are available via HTTPS at:
- https://n5eil01u.ecs.nsidc.org/HMA/HMA_DEM8m_AT.001/
- http://n5eil01u.ecs.nsidc.org/HMA/HMA_DEM8m_CT.001/
Within these directories, files are organized by the center date of the stereo-pair images, such as 2002.01.01 (in YYYY.MM.DD format).
Note: Subdirectory names (e.g. 2002.01.01) conform to NASA archiving requirements for these data, but in this case do not reflect the temporal resolution/timestamps of the data.
File names vary slightly as shown in the following examples. Variables used in the file names are defined in Table 1.
Example File Names
- HMA_DEM8m_[AT/CT]_20141226_0534_1050410011D3AB00_1050410011D3AD00.tif
- HMA_DEM8m_[AT/CT]_20141226_0534_1050410011D3AB00_1050410011D3AD00_hs_az315.png
Naming Conventions
- HMA_DEM8m_[AT/CT]_[YYYYmmdd_HHMM]_[CATID1_CATID2].tif
- HMA_DEM8m_[AT/CT]_[YYYYmmdd_HHMM]_[CATID1_CATID2]_hs_az315.png
Variable | Definitions |
---|---|
HMA_DEM8m_[AT/CT] | Data set IDs [abbreviations for High Mountain Asia (HMA) DEM 8 m Along-track or Cross-track Product] |
YYYYmmdd_HHMM | Stereo-pair center acquisition date/time [4-digit year, 2-digit month, 2-digit day, 2-digit hour, 2-digit minute), Coordinated Universal Time (UTC); indicates mean of the two image timestamps] |
[CATID1_CATID2] |
Catalog IDs of first and second DigitalGlobe images in the stereo pair; the first 3 digits indicate the DigitalGlobe sensor, including:
|
hs_az315 | Shaded relief/hillshade (hs) browse image generated at 315 azimuth (az) |
File sizes vary significantly. Individual 8 m DEMs range from approximately 7 MB to 51 MB and browse images range from approximately 1 KB to 46 KB.
Northernmost Latitude: 46° N
Southernmost Latitude: 26° N
Easternmost Longitude: 103° E
Westernmost Longitude: 67° E
Coverage extends to portions of the following eleven countries:
Afghanistan
Bhutan
China
India
Kazakhstan
Kyrgyzstan
Myanmar
Nepal
Pakistan
Tajikistan
Uzbekistan
To view the coverage of a particular DEM, refer to the Knowledge Base FAQ: What is the spatial extent of individual HMA DEMs?
Note: Initial archive request and processing was limited to images that intersected one or more of the Randolph Glacier Inventory (RGI) glacier polygons. Future releases will include expanded HMA coverage.
The DEMs have 8 m posting (output grid pixel size).
Projection and Grid
The along- and cross-track DEMs are provided in Universal Transverse Mercator (UTM) projection with a standard parallel of 0° N and central longitudes ranging from 69° N to 105° N. The output UTM zone is automatically determined based on the longitude of the stereo intersection geometry. Projected coordinates are meters, with the origin at 0° N. Elevation values are in meters above the WGS84 ellipsoid. Refer to Tables 2 and 3 for details.
Projection | Transverse Mercator, UTM Zones 12N-48N (details vary by zone; see Table 3) |
---|---|
Grid Spacing | 8 m |
Latitudes of True Scale (Standard Parallels) | 0° N |
Datum | WGS84 |
Ellipsoid | WGS84 |
Units | Meters |
UTM Zone | Central Longitude (Meridan) | EPSG Code | PROJ4 String |
---|---|---|---|
42N | 69° E | 32642 | +proj=utm +zone=42 +ellps=WGS84 +datum=WGS84 +units=m +no_defs |
43N | 75° E | 32643 | +proj=utm +zone=43 +ellps=WGS84 +datum=WGS84 +units=m +no_defs |
44N | 81° E | 32644 | +proj=utm +zone=44 +ellps=WGS84 +datum=WGS84 +units=m +no_defs |
45N | 87° E | 32645 | +proj=utm +zone=45 +ellps=WGS84 +datum=WGS84 +units=m +no_defs |
46N | 93° E/td> | 32646 | +proj=utm +zone=46 +ellps=WGS84 +datum=WGS84 +units=m +no_defs |
47N | 99° E | 32647 | +proj=utm +zone=47 +ellps=WGS84 +datum=WGS84 +units=m +no_defs |
48N | 105° E | 32648 | +proj=utm +zone=48 +ellps=WGS84 +datum=WGS84 +units=m +no_defs |
The optical images used for the DEMs were obtained between 2002 and 2017, with the majority of coverage from 2013 to 2016. For HMA stereo:
- The WorldView-1 (WV-1) record began in late 2008.
- The WorldView-2 (WV-2) record began in mid 2010.
- The WorldView-3 (WV-3) record began in late 2014.
- The GeoEye-1 (GE-1) record began in late 2009.
- The QuickBird-2 (QB-2) record began in 2002.
The temporal resolution for this data set varies. Time intervals between overlapping DEM strips can be hours to decades, depending on the location, with typical values of approximately one to three years.
Figure 2 shows a histogram of the annual DEM count.

Software and Tools
Data Acquisition and Processing
For detailed information on the various satellites and sensors used to produce this data set, refer to the DigitalGlobe, Inc. Web site. Satellite and sensor details are also provided in Shean et al. (2016).
The DEM strips were generated from Level-1B (L1B) panchromatic (450-800 nm) imagery acquired by the following DigitalGlobe satellites:
- GeoEye-1
- WorldView-1, WorldView-2, and WorldView-3
- QuickBird-2
The pushbroom swath imagery was collected in strips, ranging from approximately 13.1 km to 17.4 km wide and up to 111 km long.
All DigitalGlobe imagery was made available through the National Geospatial-Intelligence Agency's NextView License. The original L1B imagery used to generate the DEMs is archived and distributed by the Polar Geospatial Center at the University of Minnesota.
All DEMs were generated using the NASA Ames Stereo Pipeline (ASP) (v2.5.3 Build ID: 0be5052 to v2.6.0 Build ID: 4521300) and methodology described in Shean et al. (2016).. The data were processed using the NASA High-Performance Computing (HPC) resources, primarily the NASA Pleiades Supercomputer at the NASA Ames Research Center.
The DEMs in these data sets were created from two sets of images. One set was generated from along-track (also referred to as in-track) stereoscopic imagery in which both images were collected minutes apart along the same orbital pass. The other set was generated from cross-track imagery collected at different times on different orbits, and potentially from different spacecraft.
All subscenes for each 11-bit panchromatic image were mosaicked, and new rational polynomial coefficients (RPCs) were derived for each image strip. The mosaicked images were orthorectified using the 1-arcsec (30 m) SRTM-GL1, the void-filled Global Shuttle Radar Topography Mission (SRTM) DEM. Orthorectification using a low-resolution DEM removes most of the feature offsets due to terrain and viewing geometry, providing a significant performance improvement during stereo correlation.
The ASP stereo utility was run for the input orthoimages with the following parameters:
- Integer correlation kernel size of 21x21 pixels
- Parabolic sub-pixel refinement (using the ASP sub-pixel-mode 1 utility) with kernel size 11x11 pixels
- Conservative filtering to remove isolated "islands" up to 1024 pixels
- The ASP point2dem utility was used to produce gridded DEMs from the triangulated point cloud.
Along-track Stereo Strips
These are generated from in-track/along-track stereo pairs, involving images acquired on the same orbit, typically 60-90 seconds apart.
Cross-track Stereo Strips
Cross-track (also referred to as coincident mono) DEMs were generated from pairs of independent monoscopic images with geometry suitable for stereo reconstruction. Candidate pairs in the DigitalGlobe image archive were identified based on the criteria in Table 4. Some of these cross-track pairs were acquired on the same orbit, while others were acquired on different orbits or by different spacecraft (e.g., WorldView-1 and WorldView-2).
Parameter | Minimum Value | Maximum Value |
---|---|---|
Convergence angle (degrees) | 15 | 70 |
Pair time offset (days) | - | 7 |
Intersection area (km2) | 40 | - |
Intersection width (km) | 6 | - |
Intersection height (km) | 6 | - |
Cloud cover (%) | 0 | 75 |
The cross-track DEMs potentially have increased error due to residual horizontal displacement errors (i.e., errors due to ice flow between image acquisitions), non-ideal stereo geometry (e.g., smaller convergence angles) and the fact that some errors in ephemeris data for the two images are independent (as opposed to highly correlated errors for along-track pairs). In practice, these issues can result in increased DEM vertical/horizontal bias and increased relative error (e.g., DEM “tilt”). Despite potentially increased error, we include these cross-track DEMs to fill critical gaps in along-track stereo coverage.
Processing Steps
In summary, two images are collected for the same ground extent from different angles. Cross-correlation is used to measure offsets between features in the images, which are related to the different viewing geometry and topography. These offsets are combined with spacecraft attitude/ephemeris data to identify precise sensor locations, and stereo triangulation is used to determine the three dimensional coordinates of matched features. The resulting point cloud is interpolated to create the high-resolution DEM.
See Shean et al. (2016) for details of the full processing workflow. For information regarding ASP, refer to the official documentation (PDF, 21 MB).
The absolute geolocation accuracy of the DEMs depends on the sensor. In general, the WorldView and GeoEye along-track pairs will have horizontal/vertical accuracy of <5 m CE90/LE90. Quickbird-2 pairs have reduced horizontal (CE90 23 m) and vertical accuracy.
Relative error within an individual DEM strip depends on stereo-pair geometry and surface slope, but generally should be <1-2 m. Individual DEMs can contain "blunders" or local elevation errors due to image correlation errors that are not removed by subsequent filters.
Data voids can occur due to clouds, occlusions (surface not visible from both images), saturated images (e.g., overexposure for surfaces with fresh snow on illuminated slopes), errors in the low-resolution DEM used for initial orthorectification, and/or correlation failures. The investigators elected to preserve only valid DEM pixels, and do not perform any additional interpolation or filling to remove voids. Future releases will include void-filled products.
For more information on errors and limitations, refer to Shean et al. (2016).
Before comparing two DEMs, the user must remove horizontal and vertical offsets over control surfaces via co-registration (e.g., Berthier et al., 2007; Nuth and Kääb, 2011; Shean et al., 2016). This can be accomplished with the ASP pc_align utility (see the demcoreg tools https://github.com/dshean/demcoreg). Future data releases will include DEM products co-registered to satellite laser altimetry data for improved absolute accuracy.
References and Related Publications
Contacts and Acknowledgments
David Shean
University of Washington
Polar Science Center, Applied Physics Laboratory
Department of Civil and Environmental Engineering
Seattle, Washington 98195 USA
http://psc.apl.uw.edu/people/post-docs/david-shean/
Technical Contact
NSIDC User Services
National Snow and Ice Data Center
CIRES, 449 UCB
University of Colorado
Boulder, CO 80309-0449 USA
phone: +1 303.492.6199
fax: +1 303.492.2468
e-mail: nsidc@nsidc.org
These data were generated through funding from the NASA Cryosphere program and the NASA HiMAT project (NNX16AQ88G). Support for the Polar Geospatial Center was provided by the National Science Foundation (ANT-1043681). Resources supporting this work were provided by the NASA High-End Computing (HEC) Program through the NASA Advanced Supercomputing (NAS) Division at Ames Research Center.
Document Information
Document Creation Date
December 2017
Document Revision Date
N/A