Generate elevation data using the DEMs wizard

Available with Advanced license.

Elevation data can be derived from stereo image pairs using photogrammetry algorithms. A stereo pair is comprised of two images of the same geolocation taken from different perspectives.

Stereo image pairs collected from different perspectives

The stereo image pairs of an image collection are used to generate a point cloud (3D points) from which elevation data can be derived. The derived elevation data will be used to orthorectify the image collection in the ortho mapping workspace.

Note:

Nadir (vertical) viewing, high-resolution satellite imagery is not affected much by distortion that is inherent in aerial imagery due to the large distance between sensor and ground, long sensor focal length (on the order of 10 meters), and small field of view. These factors, together with accurate orientation information in the form of rational polynomial coefficients (RPCs), result in the condition that DEM accuracy and dense postings are less important in producing accurate orthoimages, as long as the adjusted exterior orientation and control points are adequate. Thus, the DEM generation step is often not used, and existing USGS NED DEMs or SRTM DEMs, together with accurate GCPs, can produce Class I or Class II orthoimages at a scale of 1:5,000 or smaller.

The 3D points generated from stereo pairs in the DEMs wizard are classified into two categories:

  • Digital terrain model (DTM)—Digital elevation of the earth, not including the elevation of any objects on it. This is also referred to as bare-earth elevation. The bare earth DTM dataset is used to produce the orthoimage and orthomosaics.
  • Digital surface model (DSM)—Digital elevation of the earth, including the elevation of objects on it such as trees and buildings. The DSM is a valuable analytical dataset used for classifying features in orthoimages, such as discriminating asphalt pavement and asphalt roofs. It should not be used for image orthorectification unless the source imagery is nadir looking, with no building or feature lean, to produce true orthoimages.

DTM and hillshaded DSM

Note:

If a forest area is heavily wooded, or has other dense vegetation cover, it will not be possible to derive a DTM ground surface because the ground is not visible. The most appropriate elevation surface product for densely forested land cover is a DSM, which specifically creates a surface depicting the top of the tree canopy.

Elevation can be derived when the image collection has a good amount of overlap to form the stereo pairs. Typical image overlap to produce point clouds is 80 percent forward overlap along a flight line and 60 percent overlap between flight lines, so that every location on the ground is covered by multiple images. It is often used in processing drone images, digital aerial images, or certain satellite images captured specifically for stereo applications.

This wizard provides two preconfigured steps to produce the outputs:

  1. Construct stereo pairs from the image collection to compute the point cloud.
  2. Interpolate a raster at a user-defined resolution from the point cloud.

You can modify the default processing parameters, but you cannot remove a step. If you wish to omit or perform a specific step, you can use the Custom wizard.

Point Cloud Settings page

First, specifying the parameters for the Point Cloud Settings page.

Parameters for Point Cloud Settings

Parameter nameDescription

Matching Method

There are three matching methods for generating point cloud:

  • ETM (Extended Terrain Matching)—is a feature-based stereo matching technique in which the Harris operator is used in detecting feature points. Since fewer feature points are extracted, this method is fast and can be used for data with less terrain variation and detail. This is the default.1
  • SGM (Semi-Global Matching)—produces points that are denser and provide more detailed terrain information. It can be used for images of urban areas. This is more computational intensive than ETM.2
  • MVM (Multi-View Matching)— is based on the SGM matching method followed by a fusion step in which the redundant elevation estimations across single stereo model are merged. It produces dense 3D points and is computationally efficient.3

Maximum Object Size (in meters)

Defines a search radius that is used to filter out objects above ground. Objects smaller than the threshold will be filtered as ground, otherwise objects will be treated as above-ground features such as buildings, bridges, or trees. The default object size is 10 meters.

Point Ground Spacing

Defines the spacing, in meters, at which the 3D points are generated.

The suggested spacing is five times the source image pixel size.

Minimum Intersection Angle (in degrees)

The point cloud is generated from stereo pairs. This value, in degrees, defines the minimum angle the stereo pair must meet. The default is 5 degrees.

Note:

A stereo pair with too small of an intersection angle will produce unstable results when triangulating 3D points.

Maximum Intersection Angle (in degrees)

The point cloud is generated from stereo pairs. This value, in degrees, defines the maximum angle the stereo pair must meet. The default is 70 degrees.

Note:

A stereo pair with too large of an intersection angle will produce few or no match points.

Minimum Area Overlap

The percentage of overlapping area over the whole image. The default is 0.6.

Maximum Omega / Phi Difference (in degrees)

The maximum threshold for the Omega/Phi difference between the two image pairs. The Omega values and Phi values for the image pairs are compared. If the difference between either the two Omega or the two Phi values is above the threshold, the pairs will not be formatted as a stereo pair.

Maximum GSD Difference

The threshold for the maximum ground sample distance (GSD) between two images in a stereo pair. If the resolution ratio between the two images is greater than the threshold value, the pairs will not be built as a stereo pair. The default is 2.

Number of Image Pairs

The number of pairs used to generate 3D points. For a project that has dense overlaps and many stereo pairs, increasing this number means more computation time. The suggested value is 4.

Sometimes a location may be covered with many image pairs. In this case, the tool will order the pairs based on the various threshold parameters specified in this tool. The pairs with the highest scores will be used to generate the points.

This parameter limits one pair from being used too many times. The parameters that affect the order of the stereo pair, besides Minimum Intersection Angle, Maximum Intersection Angle, Minimum Area Overlap, may also include Maximum Omega / Phi Difference, Maximum GSD Difference, and Adjustment Quality Threshold.

Adjustment Quality Threshold

Specify the minimum adjustment quality that is acceptable. The threshold value will be compared to the adjustment quality value that is stored within the stereo model. Image pairs with an adjustment quality less than the specified threshold will receive a score of 0 for this criteria and will descend in the ordered list. The values for the threshold ranges from 0 to 1. The suggested value is 0.2.

References

  1. Harris, Christopher G., and Mike Stephens. "A combined corner and edge detector." Alvey vision conference, vol. 15, no. 50, pp. 10-5244. 1988.
  2. Hirschmuller, Heiko, Maximilian Buder, and Ines Ernst. "Memory Efficient Semi-Global Matching." ISPRS Annals of the Photogrammetry, Remote Sensing and Spatial Information Sciences, Volume 1-3, (2012): 371-376.
  3. Hirschmuller, Heiko. "Stereo Processing by Semiglobal Matching and Mutual Information." IEEE Transactions on pattern analysis and machine intelligence, Volume 30, no. 2 (2007): 328-341.

DEM Interpolation Settings page

Next, specify the parameters for the DEM Interpolation Settings page.

Parameters for DEM Interpolation Settings

Parameter nameDescription

Surface Type

Create a digital terrain model or a digital surface model.

  • DTM—Create a digital terrain model by interpolating only the raster surface using ground only points.
  • DSM—Create a digital surface model by interpolating a raster using all the points; both ground and above around points.

Cellsize

The cell size of the output raster dataset.

Format

The format for the output raster dataset:

  • Cloud Raster Format. This is the default.
  • TIFF Format

Compression

The method used to compress the output raster dataset.

  • None—Do not compress the output raster dataset. This is the default.
  • LERC—Use LERC compression to compress the output raster dataset. LERC is a compression method that works on any data type, such as byte, integer, real, and double precision. The efficiency of the compression algorithm increases with the pixel depth.

Max Error

The maximum error allowed for LERC compression. The maximum error is a tolerance value that is applicable per pixel (not an average for the image).

Set the Compression to LERC to use this option.

Interpolation Method

The method used to interpolate the output raster dataset from the point cloud.

  • TIN Linear Interpolation—Also known as triangulated irregular network (TIN), linear interpolation is designed for irregularly distributed sparse points, such as solution points from block adjustment computation. This is the default.
  • TIN Natural Neighbor Interpolation—This is similar to triangulation, but generates a smoother surface and is more computationally intensive.
  • Inverse Distance Weighted Average Interpolation—This is used for regularly distributed dense points, such as point cloud LAS files from the Generate Point Cloud tool. The IDW search radius is automatically computed based on average point density.

Smoothing Method

Select a filter to smooth the output raster dataset.

  • Gaussian 3 by 3—A Gaussian filter with a 3 by 3 window.
  • Gaussian 5 by 5—A Gaussian filter with a 5 by 5 window. This is the default.
  • Gaussian 7 by 7—A Gaussian filter with a 7 by 7 window.
  • Gaussian 9 by 9—A Gaussian filter with a 9 by 9 window.
  • No smoothing—No smoothing filter is applied.

Fill Missing Pixels Using

Specify a DEM input that is used to fill NoData areas.

Note:

Areas of NoData may exist where stereo overlap is insufficient, or no matching points are found in the area during point cloud generation.

Related topics