Module Normals
See also
opals::INormals

Aim of module

Estimates local surface normal vectors for each point of a specific point cloud (ODM).

General description

Estimation of local plane information on a point basis is useful (or even necessary) for different tasks. Parameters that are derived by opalsNormal can be used to compute local slope and exposure values, roughness measures, values describing the point distribution and are input for other OPALS modules (e.g. opalsEchoRatio).

The module attaches the computed information to the input OPALS Data Manager (ODM). Different neighbour selection strategies, different plane estimation methods but also a very flexible way to process only a subset of input points are provided. Additionally it is possible to define the level of meta information that are stored with each point. Depending on the subsequent task appropriate parameter have to be set.

Fig. 1: Normals of buildings and vegetation visualised with ParaView

Normals estimator

Currently three different normal estimator (see parameter normalsAlg ) are provided which can be mainly differentiate by their degree of robustness. simplePlane is a standard least-squares plane adjustment offering no blunder detection (fastest method). Hence, it should be used if the normal quality is of minor concern or the focus is set on a simple roughness detection strategy (non-planar region will lead to a higher sigma0 of the adjustment). robustPlane detects blunders using observation weighting based on their residual. After convergence all blunder are eliminated and a final adjustment, using equally weighted observations, is performed. The robust fmcd (Fast Minimum Covariance Determinant) estimator uses an additional inverse distance weighting strategy for better approximation of the local point environment. This complex estimation method is appropriate for very dense and very inhomogeneous point distributions.

If a robust estimator is selected, the number of neighbours (see parameter neighbours) should be clearly increased (>10). To reduce extrapolation effects it is possible to make a quadrant (or octant) orientated neighbbour selection. For more details on the quadrant selection strategie see section data selection strategies. Additionally it should be mentoined that the searchMode and selMode parameters are independently and can be unrestricted combined. Hence, it is possible to select the nearest neigbhours in 3d and apply a quadrant oriented selection (and vice versa).

Stored result (level of meta information)

The normal estimation deliver more interesting data than just the normal vector itself. Beside some statistical values (sigma0, number of used observations, ...) specially the eigenvalues of the covariance matrix can be interesting for subsequent processes. Since the computed data accrue for each point, storing all values can dramatically increase the size of the used ODM files (twice the original size or even more). Therefore, the module provides the parameter storeMetaInfo which allows setting the level of meta data that are stored with each point.

  • minimum: stores NormalX, NormalY, NormalZ, NormalSigma0 and the estimation method (NormalEstimationMethod)
  • medium: stores all minimum values + the three eigenvalues (NormalEigenvalue1-3) of the covariance matrix, plane offset, number of given points (NormalPtsGiven) and number of used points (NormalPtsUsed)
  • maximum: stores all medium values + the corresponding eigenvectors

For additional details of the normal attributes see the predefined attributes of the ODM.

The orientation of the normal vector can be controlled by the parameter direction. As default the normals are oriented upwards towards the positive z-axis (NormalZ >= 0). This can lead to flipping orientations on vertical (or near vertical) surfaces (see opals::NormalsDirection for more options).

For completeness it should be mentioned that NormalSigma0 is only set if the redundancy of the adjustment is greater than zero (otherwise it is set to null). If this sigma is required by subsequent processes (like Module EchoRatio) the number of neighbours should be carefully chosen, especially, if robust estimators are used.

Process only a subset of data (filtering)

Within Module Normals the powerful OPALS filters can be used in two different ways. The first filter defines in which points a normal should be estimated. The second filter entry selects the points the are used for the neighbour search. If only one filter is provided, it is used for both filtering stages. Hence, it is possible to estimate normals for single echoes only but all last echoes are used as neighbouring information. Situations are imagineable where different estimators should be used for different subsets. Or a second (or third) run of Module Normals with adapted parameters is performed to compute normals for points where previous runs have failed (without changing existing results). All these issues can be handled by setting appropriate filter strings. For details see the Examples section.

Parameter description

-inFileodm input file
Type: opals::Path
Remarks: mandatory
Specifies the opals datamanager input file
-filterfilter string
Type: opals::Vector<opals::String>
Remarks: optional
The first filter specifies for which points a normal should be estimated (If not defined all echos are used). The second filter defines which points are used for the plane fit (If not set, the same filter as value one is used)
-normalsAlgnormals algorithm
Type: opals::NormalsAlgorithm
Remarks: default=simplePlane
-neighboursnr of nearest neighbours used for normals estimation
Type: int32
Remarks: default=4
-searchRadiusmaximum search radius for point selection
Type: float
Remarks: optional
Estimation rule: If not set, the search radius is estimated using the average point density and the number of neighbours
-searchModedimension of nearest neighbor search
Type: opals::SearchMode
Remarks: default=d3
Possible values:  
  d2 ... Search based on 2D coordinates (x and y) only 
  d3 ... Search based on full 3D coordinates (x,y and z)
Defines if 2D or 3D distances should be used for the nearest neighbor search
-selModedata selection mode
Type: opals::SelectionMode
Remarks: default=nearest
Please note that quadrant or octant wise data selection on the one hand reduces extrapolation effects, but on the other hand decreases the processing performance
-storeMetaInfolevel of meta info that are stored for each point
Type: opals::NormalsMetaInfo
Remarks: default=minimum
-directionuniquely define the direction of the resulting normals
Type: opals::NormalsDirection
Remarks: default=upwards
Normal directions are generally only defined up to +/- pi.
-debugOutFilenormals debug output as winput file
Type: opals::Path
Remarks: optional
Outputs a winput file containing the normals as structure lines and points where the estimation failed as bulk data
-debugNormalsLenlength of the normal vectors in the debug output file
Type: float
Remarks: optional
This parameter only effects the length of the normal vectors in the debug output file. The stored attribute values are unchanged.

Examples

The data used in the following example are located in the $OPALS_ROOT/demo/ directory. The example shows how to use opalsNormals including a short analysis of the results.

As a prerequisite, the ALS point cloud data must be imported into the ODM. To achieve that, change to the demo directory and type:

opalsImport -infile strip21.laz

Now, run the following command


opalsNormals -inFile strip21.odm -neighb 8 -searchMode d3 -selmode quadrant

to perform the local plane estimation selecting 8 neighbours in an 3d neighbourhood. To avoid extrapolations the selection mode is set to quadrant.

One possiblity to check the computed normals is a visualization of the vector output, e.g. using the debugOutFile parameter. On a larger scale, however, opalsHisto is a valuable tool for similar verifications. Executing the following commands


opalsHisto -infile strip21.odm -attribute NormalSigma0 -sampleRange 0 0.5
opalsHisto -infile strip21.odm -attribute NormalY

produces two histograms as shown in Figure 2. The data set covers an urban area (see Figure 3) with little degree of vegetation. This accords to the textual output from this first opalsHisto run. It can be seen that 75% of all sigma0 values are below 6 [cm]. In the histogram of the the y component of the normal vectors (right diagram in Figure 2) two interesting bumps (at position +/- 0.6) are cleary visible. Those result from the significant amount of saddle roofs with an east-west ridge.

Fig. 2: Histogram of sigma0 of plane estimations (left: NormalSigma0) and the y component of normal vectors (right: NormalY)
Fig. 3: Z coding of strip21.laz

References

Dorninger, P. and Nothegger, C., 2007. 3D Segmentation of Unstructured Point Clouds for Building Modelling. In: IAPRS XXXVI Part 3/W49A, pp. 1919996.

Author
jo
Date
23.06.2016
@ quadrant
quadrant-wise nn selection, ie. nn per quadrant, then 2nd nn per quadrant, ...
opalsNormals is the executable file of Module Normals
Definition: ModuleExecutables.hpp:148
@ NormalSigma0
Sigma0 of normal estimation.
opalsImport is the executable file of Module Import
Definition: ModuleExecutables.hpp:113
@ sampleRange
sample (attribute) range (e.g. opalsHisto)
opalsHisto is the executable file of Module Histo
Definition: ModuleExecutables.hpp:103
@ searchMode
dimension of nearest neighbor search (opalsNormals)
@ NormalY
y value of normal (unit) vector
@ d3
Search based on full 3D coordinates (x,y and z)