GEOG 482
The Nature of Geographic Information

9. Interpolation


DEMs are produced by various methods. The method preferred by USGS is to interpolate elevations grids from the hypsography and hydrography layers of Digital Line Graphs.

Hypsography and hydrography layers of a Digital Line Graph viewed in Global Mapper software
Figure 7.9.1 A USGS 7.5-minute DEM and the DLG hypsography and hydrography layers from which it was produced.

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 Figure 7.9.2 for example, the gridded elevations shown in purple were interpolated from the irregularly spaced spot elevations shown in red.

A grid of elevation values that were interpolated from an irregularly-spaced array
Figure 7.9.2 Elevation values in DEMs are interpolated from irregular arrays of elevations measured through photogrammetric methods, or derived from existing DLG hypsography and hydrography data.

Here's another example of interpolation for mapping. The map belowin Figure 7.9.3 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° longitude by 2.5° latitude.

A map showing gridded temperature data
Figure 7.9.3 1995 Surface Temperature Anomalies.
National Climatic Data Center, 2005

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. The irregular array is transformed into a regular array through interpolation. In general, interpolation is the process of estimating an unknown value from neighboring known values.

Locations of temperature climate records used to create a gridded temperature map
Figure 7.9.4 The Global Historical Climate Network.
Eischeid et al., 1995

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 "mass points." 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.

Three number lines illustrating how interpolation is affected by assumptions about the underlying distribution
Figure 7.9.5 Interpolating an intermediate value on a number line.

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 "?" on the top number line, what would you guess? An estimate of "5" 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 "?" 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.

As I mentioned in Chapter 1, the surface of the Earth is characterized by a property called spatial dependence. Nearby locations are more likely to have similar elevations than are distant locations. Spatial dependence allows us to assume that it's valid to estimate elevation values by interpolation.

Many interpolation algorithms have been developed. One of the simplest and most widely used (although often not the best) is the inverse distance weighted algorithm. 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.

Diagram and formula explaining inverse distance weighted interpolation
Figure 7.9.6 The inverse distance weighted interpolation procedure.

Practice Quiz

Registered Penn State students should return now to the Chapter 7 section of the Modules page in Canvas to take a self-assessment quiz about Interpolation. You may take practice quizzes as many times as you wish. They are not scored and do not affect your grade in any way.

Penn State logo
Students who register for this Penn State course gain access to assignments and instructor feedback, and earn academic credit. Information about Penn State's Online Geospatial Education programs is available at the Geospatial Education Program Office.