Generate Point Cloud (Data Management)

Summary

Generates 3D points from stereo pairs and outputs a point cloud as a set of LAS files.

The tiling of the LAS files is based on 1,000 by 1,000 ground spacing. The points in each LAS tile are computed by selecting pairs, based on user-defined criteria, and filter points from the selected pairs. The input of this tool is a mosaic dataset that contains a stereo model. The output of this tool can be used to generate a digital terrain model (DTM) or a digital surface model (DSM).

Usage

  • If this tool is run multiple times with the same input parameters, the output may be slightly different due to random sampling.

  • The order of the stereo pairs in the Number of Image Pairs parameter is first based on the values defined for the Adjustment Quality Threshold, GSD Difference Threshold, and Omega/Phi Difference Threshold parameters. Each of these thresholds will give a pair of images a score of 0 if it did not meet the threshold, or a score of 1 for each threshold that was met, for a maximum of 3. The highest scores will be ordered at the top of the priority list. Next, the Overlap Area Threshold value is taken into account for any pairs with the same scores, and the intersection angle between the image pairs will be used as the last criteria to order the list, when the larger intersection angle is ordered higher.

  • If you want a specific pair to be used for point cloud generation, set a high value for this pair in the Use field of the stereo table. To open the stereo table, right-click the mosaic layer in the Contents pane and click Open > Stereo.

Parameters

LabelExplanationData Type
Input Mosaic Dataset

The input mosaic dataset, which must have completed the block adjustment process and have a stereo model.

To block adjust the mosaic dataset, use the Apply Block Adjustment tool. To build a stereo model on the mosaic dataset, use the Build Stereo Model tool.

Mosaic Dataset; Mosaic Layer
Matching Method

Specifies the method that will be used to generate 3D points.

References:

  • Heiko Hirschmuller et al., "Memory Efficient Semi-Global Matching," ISPRS Annals of the Photogrammetry, Remote Sensing and Spatial Information Sciences, Volume 1–3, (2012): 371–376.
  • Hirschmuller, H. "Stereo Processing by Semiglobal Matching and Mutual Information." Pattern Analysis and Machine Intelligence, (2008).

  • Extended terrain matchingExtended terrain matching, a feature-based stereo matching in which the Harris operator is used in detecting feature points, will be used as the matching method. Since fewer feature points are extracted, this method is fast and can be used for data with fewer terrain variations and less detail.
  • Semiglobal matchingSemi-Global Matching, which produces points that are denser and have more detailed terrain information, will be used as the matching method. It can be used for images of urban areas. This method is more computational intensive than ETM.
  • Multi-view image matchingMulti-view image matching, which is based on the SGM matching method followed by a fusion step in which the redundant depth estimations across a single stereo model are merged, will be used as the matching method. It produces dense 3D points and is computationally efficient.
String
Output LAS Folder

The folder used to store the output LAS files, including cloud storage.

If this tool is run multiple times with the same input parameters, the output may be slightly different due to random sampling.

Folder
Output LAS Base Name

A string used as a prefix to formulate the output LAS file names. For example, if name is used as the base, the output files will be named name1.las, name2.las, and so on.

String
Maximum Object Size (in meter)
(Optional)

A search radius within which surface objects, such as buildings or trees, will be identified. It is the linear size in map units.

Double
DSM Ground Spacing (in meter)
(Optional)

The ground spacing, in meters, at which the 3D points will be generated.

The default is five times the source image pixel size.

Double
Number of Image Pairs
(Optional)

The maximum number of image pairs that an image can contribute to generate 3D points. The default value is a minimum of 8 image pairs.

If the image is involved in more image pairs than specified, those image pairs will not be considered when constructing the 3D points. In this case, the tool will order the pairs based on the various threshold parameters specified in the tool. The pairs with the highest scores will be used to generate the points.

Double
Overlap Area Threshold
(Optional)

The minimum overlap threshold area that is acceptable, which is a percentage of overlap between a pair of images. Image pairs with overlap areas smaller than this threshold will receive a score of 0 for this criteria and will descend in the ordered list. The range of values for the threshold is 0 to 1. The default threshold value is 0.5, which is equal to 50 percent.

Double
Adjustment Quality Threshold
(Optional)

The minimum adjustment quality that is acceptable. The threshold value will be compared to the adjustment quality value that is stored in 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 range of values for the threshold is 0 to 1. The default value is 0.2, which is equal to 20 percent.

Double
GSD Difference Threshold
(Optional)

The maximum allowable threshold for the ground sample distance (GSD) between two images in a pair. The resolution ratio between the two images will be compared to the threshold value. Image pairs with a ground sample ratio greater than this threshold will receive a score of 0 for this criteria and will descend in the ordered list. The default threshold ratio is 2.

Double
Omega/Phi Difference Threshold
(Optional)

The maximum threshold for the difference between the Omega values and Phi values for the two image pairs. The Omega values and Phi values for the image pairs are compared. Image pairs with an Omega or a Phi difference greater than this threshold will receive a score of 0 for this criteria and will descend in the ordered list. The default threshold difference for each comparison is 8.

Double

arcpy.management.GeneratePointCloud(in_mosaic_dataset, matching_method, out_folder, out_base_name, {object_size}, {ground_spacing}, {minimum_pairs}, {minimum_area}, {minimum_adjustment_quality}, {maximum_diff_gsd}, {maximum_diff_OP})
NameExplanationData Type
in_mosaic_dataset

The input mosaic dataset, which must have completed the block adjustment process and have a stereo model.

To block adjust the mosaic dataset, use the Apply Block Adjustment tool. To build a stereo model on the mosaic dataset, use the Build Stereo Model tool.

Mosaic Dataset; Mosaic Layer
matching_method

Specifies the method that will be used to generate 3D points.

  • ETMExtended terrain matching, a feature-based stereo matching in which the Harris operator is used in detecting feature points, will be used as the matching method. Since fewer feature points are extracted, this method is fast and can be used for data with fewer terrain variations and less detail.
  • SGMSemi-Global Matching, which produces points that are denser and have more detailed terrain information, will be used as the matching method. It can be used for images of urban areas. This method is more computational intensive than ETM.
  • MVMMulti-view image matching, which is based on the SGM matching method followed by a fusion step in which the redundant depth estimations across a single stereo model are merged, will be used as the matching method. It produces dense 3D points and is computationally efficient.

References:

  • Heiko Hirschmuller et al., "Memory Efficient Semi-Global Matching," ISPRS Annals of the Photogrammetry, Remote Sensing and Spatial Information Sciences, Volume 1–3, (2012): 371–376.
  • Hirschmuller, H. "Stereo Processing by Semiglobal Matching and Mutual Information." Pattern Analysis and Machine Intelligence, (2008).

String
out_folder

The folder used to store the output LAS files, including cloud storage.

If this tool is run multiple times with the same input parameters, the output may be slightly different due to random sampling.

Folder
out_base_name

A string used as a prefix to formulate the output LAS file names. For example, if name is used as the base, the output files will be named name1.las, name2.las, and so on.

String
object_size
(Optional)

A search radius within which surface objects, such as buildings or trees, will be identified. It is the linear size in map units.

Double
ground_spacing
(Optional)

The ground spacing, in meters, at which the 3D points will be generated.

The default is five times the source image pixel size.

Double
minimum_pairs
(Optional)

The maximum number of image pairs that an image can contribute to generate 3D points. The default value is a minimum of 8 image pairs.

If the image is involved in more image pairs than specified, those image pairs will not be considered when constructing the 3D points. In this case, the tool will order the pairs based on the various threshold parameters specified in the tool. The pairs with the highest scores will be used to generate the points.

Double
minimum_area
(Optional)

The minimum overlap threshold area that is acceptable, which is a percentage of overlap between a pair of images. Image pairs with overlap areas smaller than this threshold will receive a score of 0 for this criteria and will descend in the ordered list. The range of values for the threshold is 0 to 1. The default threshold value is 0.5, which is equal to 50 percent.

Double
minimum_adjustment_quality
(Optional)

The minimum adjustment quality that is acceptable. The threshold value will be compared to the adjustment quality value that is stored in 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 range of values for the threshold is 0 to 1. The default value is 0.2, which is equal to 20 percent.

Double
maximum_diff_gsd
(Optional)

The maximum allowable threshold for the ground sample distance (GSD) between two images in a pair. The resolution ratio between the two images will be compared to the threshold value. Image pairs with a ground sample ratio greater than this threshold will receive a score of 0 for this criteria and will descend in the ordered list. The default threshold ratio is 2.

Double
maximum_diff_OP
(Optional)

The maximum threshold for the difference between the Omega values and Phi values for the two image pairs. The Omega values and Phi values for the image pairs are compared. Image pairs with an Omega or a Phi difference greater than this threshold will receive a score of 0 for this criteria and will descend in the ordered list. The default threshold difference for each comparison is 8.

Double

Code sample

GeneratePointCloud example 1 (Python window)

This is a Python sample for the GeneratePointCloud function.

import arcpy
arcpy.management.GeneratePointCloud('c:/data/BD.gdb/SpringMD', 'ETM',
                                    'c:/data/output', 'SpringLAS', '10')

Related topics