This is machine translation

Translated by Microsoft
Mouse over text to see original. Click the button below to return to the English verison of the page.


Read USGS 1-degree (3-arc-second) Digital Elevation Model


[Z,refvec] = usgsdem(filename,scalefactor)
[Z,refvec] = usgsdem(filename,scalefactor,latlim,lonlim)


[Z,refvec] = usgsdem(filename,scalefactor) reads the specified file and returns the data in a regular data grid along with referencing vector refvec, a 1-by-3 vector having elements [cells/degree north-latitude west-longitude] with latitude and longitude limits specified in degrees. The data can be read at full resolution (scalefactor = 1), or can be downsampled by the scalefactor. A scalefactor of 3 returns every third point, giving 1/3 of the full resolution.

[Z,refvec] = usgsdem(filename,scalefactor,latlim,lonlim) reads data within the latitude and longitude limits. These limits are two-element vectors with the minimum and maximum values specified in units of degrees.


The U.S. Geological Survey has made available a set of digital elevation maps of 1-degree quadrangles covering the contiguous United States, Hawaii, and limited portions of Alaska. The data is on a regular grid with a spacing of 30 arc-seconds (or about 100-meter resolution). 1-degree DEMs are also referred to as 3-arc-second or 1:250,000 scale DEM data.

The data is derived from the U.S. Defense Mapping Agency's DTED-1 digital elevation model, which itself was derived from cartographic and photographic sources. The cartographic sources were maps from the 7.5-minute through 1-degree series (1:24,000 scale through 1:250,000 scale).


Read every fifth point in the file containing part of Rhode Island and Cape Cod:

[Z,refvec] = usgsdem('providence-e',5);

Read the elevation data for Martha's Vineyard at full resolution:

[Z,refvec] = usgsdem('providence-e',1,...
	[41.2952 41.4826],[-70.8429 -70.4392]);
whos Z

  Name        Size         Bytes  Class

  Z         226x485       876880  double array

More About

collapse all


The grid for the digital elevation maps is based on the 1984 World Geodetic System (WGS84). Older DEMs were based on WGS72. Elevations are in meters relative to National Geodetic Vertical Datum of 1929 (NGVD 29) in the continental U.S. and local mean sea level in Hawaii.

The absolute horizontal accuracy of the DEMs is 130 meters, while the absolute vertical accuracy is ±30 meters. The relative horizontal and vertical accuracy is not specified, but is probably much better than the absolute accuracy.

These DEMs have a grid spacing of 3 arc-seconds in both the latitude and longitude directions. The exception is DEM data in Alaska, where latitudes between 50 and 70 degrees North have grid spacings of 6 arc-seconds, and latitudes greater than 70 degrees North have grid spacings of 9 arc-seconds.

Statistical data in the files is not returned.

You can obtain the data files from the U.S. Geological Survey and from commercial vendors. Other agencies have made some local area data available online.

    Note   For details on locating map data for download over the Internet, see the following documentation at the MathWorks Web site: Finding Geospatial Data.

Introduced before R2006a

Was this topic helpful?