# Geological structure measurement by LiDAR

Geological structure measurement by LiDAR technology is a remote sensing method applied in structural geology. It enables monitoring and characterisation of rock bodies. This method's typical use is to acquire high resolution structural and deformational data for identifying geological hazards risk, such as assessing rockfall risks or studying pre-earthquake deformation signs.
Geological structures are the results of tectonic deformations, which control landform distribution patterns. These structures include folds, fault planes, size, persistence, spatial variations, and numbers of the rock discontinuities in a particular region. These discontinuity features significantly impact slope stability, causing slope failures or separating a rock mass into intact rock blocks (rockfall). Some displaced blocks along faults are signs of earthquakes.
Conventionally, geotechnical engineers carried out rock discontinuity studies manually. In post geological hazards studies, such as rockfall, the rockfall source areas are dangerous and are difficult to access, severely hindering the ability to carry out detailed structural measurements and volumetric calculations necessary for hazard assessment. By using LiDAR, geological structures can be evaluated remotely, enabling a 3-D investigation of slopes with virtual outcrops.
LiDAR technology (Light Detection and Ranging) is a remote sensing technique that obtains precise 3-D information and distance. The laser receptor calculates the distance by the travelling time between emitting and receiving laser pulses. LiDAR produces topographic maps, and it is useful for assessing the natural environment.

## Importance of measuring geological structures by LiDAR

Geological structures are responsible for providing distinct physical properties to rock masses. Discontinuous properties and plate tectonic forces may alter rock masses and their geometries. These structures contain joints, fractures, bedding planes, shear zones, mechanical breaks, or any other features ranging from microscopic (<1 cm, foliation development by metamorphism) to macroscopic scale (>100m, mid-oceanic ridges).
Geological structures are typically elongated, their orientations are often described as "strike". If a rock body is extensively tilted, taking its slope resistivity into account, it may have a high potential to cause rockfalls. The use of LiDAR in the structural analysis allows measuring landform features from a single outcrop scale to a terrestrial scale. Some geological structures measurement and their importance are listed below:

### Rock plane orientation measurement and rockfall risk assessment

Rock plane orientations are the natural inclinations that occurred on a rock plane. Some examples of rock planes are bedding planes, fault planes. The planes' orientations are measured by dip and dip direction with a clinometer and compass, where dip represents the maximum inclination of a plane to the horizontal, dip direction is the direction of the intersection line between horizontal and the inclined plane. A stereonet can visualise the distribution of dip and dip directions to analyse the kinematics of a slope.
Kinematics represents the motion of a rock body without external forces that cause them to move. Kinematics analysis concentrates on the possibility of translational failures due to planes sliding although other types of failure modes, such as wedge and toppling failures, can also be recognized.

### Faults behaviours measurement and earthquakes predictions

Faults behaviours can be used to measure the rate of sediment transportation and predict earthquakes. An earthquake can contribute to the formation of faults scraps. One side of a block will be relatively upthrown, causing vertical displacements. Therefore, given the parameters of fault scraps, structural geologists are able to trance the age of it and deduce the time involved to form such features.
Earthquakes are initiated by slow slips. Slips are the displaced blocks along two sides of a fault. However, these slips are undetectable by seismometers (maximum 5mm/day). When the slipping blocks reach a critical rupture velocity, the faults would gradually evolve into a final quake size by linear acceleration along fault planes. The critical displacement of faults is proportional to the initial rupture velocities.
After collecting LiDAR data from pre-earthquake and post-earthquake landforms, by constructing 3-D digital terrain models, the displacement and deformations can be derived. Thus, scientists can predict the final earthquake scale in the future by determining faults and slips' characteristics and areas' size, and short-term earthquake predictions are possible.

### Surface processes and geological mapping

When carrying out geological mapping, interpretations through aerial photographs and satellite imagery are often used, but forest vegetation remained the major challenge for mapping. For example, characterising physical landform features at ridges and valleys are somehow complicated, many of these features are forest-covered. The topographic maps are then constructed by obtaining data manually.
LiDAR provides the full-waveform system, it enables the laser pulse to penetrate canopies and vegetations. This system allows obtaining bare-ground geological data points. Webster et al. have discovered new craters in Northern Canada by harnessing LiDAR data and digital terrain models. A digital terrain model has to be constructed to measure structural parameters (tilt angles, river incision depths). With the precise bedrock and surficial lithology mapping by LiDAR, structural geologists can reconstruct surficial processes involved.

### Traditional structural measurement

Traditional structural orientations can only be assessed on reachable exposed rock mass manually. Conventionally, engineering geologists investigated rock discontinuities studies with only a limited number at a time. The discontinuities may not represent the whole outcrop. Thus, the traditional rock plane orientation measurement may contain bias.
Geotechnical studies also investigate other geomechanical parameters, such as persistence, block size and rock joint spacing.

## LiDAR technology

LiDAR (Light Detection and Ranging) is a rapid surveying process that emits and receives laser pluses to acquire 3-D information. By illuminating lights with different wavelengths to the object of interest, LiDAR can be used to create precise topographic maps, with applications in: geology, geomorphology, surveying and other applications. Topographic maps are possible because of the Inertial Measuring Unit and Global Positioning System. Furthermore, it is a technology that can carry out studies on a steep slope and rock-cliffs.
Accurately measured data is necessary for the LiDAR data to be geo-referenced such as locating the data in a local or global coordinates system. Therefore, the produced LiDAR can be overlaid onto the aerial photographs collected previously to observe the topography changes over time.

### Principle of LiDAR

LiDAR system emits pulsed and continuous-wave lasers to acquire 3-D information. The laser scanner is the main component of LiDAR. Lasers with a wavelength of 550-600 nm are used on a ground-based system (handheld laser scanning and terrestrial laser scanning), whereas airborne systems use lasers with 1000-1600 nm wavelength.
Laser calculates the reachable range by the following formula:
$R={1 \over 2}\cdot c\cdot t_{s}$
* R represents Range in meters
* c represents speed of light in m/s
* t is the travelling time of a laser pulser in seconds
LiDAR receives information by discrete and full-waveform return. Full-waveform (multi-return) is often used for forest analysis by Airborne LiDAR, while discrete return (single return) is used by a ground-based laser scanning method. A laser is reflected whenever it reaches any surfaces. The full-waveform return is able to penetrate down into canopies and return vegetation information at different heights. Discrete return can only return the superficial materials. Thus, airborne LiDAR is often used for forestry studies.

### Data representation and data format

LiDAR data is mainly stored in a point cloud format(.las). The captured point cloud data store X, Y Z geometric data. Each data point is obtained from a single laser scan and represents a local geo-referenced spatial datum. It can represent realistic and three-dimensional rock faces in a remote and inaccessible natural terrain.
The LiDAR data have the following parameters:
* Each data point contains Co-ordinates (X, Y, Z co-ordinates in a local reference system)
* Each point has its associated color information(R, G, B format)
* The reflected intensity of the surface
* Possible automatic superposition of photographs onto point clouds (captured during scanning process)
* Any other hyperspectral data
These data help rock body features analysis, they include recorded geometrical or radiometric information of natural, excavated, or blasted rock slopes.
Previously, the LiDAR data are in the form of the American Standard Code for Information Interchange format(ASCII), which has several problems :
1) Low reading and interpreting speed of ASCII files
2) Useful data loss during data processing
3) ASCII is unstandardised
After 2003, the American Society for Photogrammetry and Remote Sensing (ASPRS) has standardized the LiDAR data in sequential binary laser(.las file) containing LiDAR or other point cloud data records.

### Geo-referencing

Geo-referencing means the co-ordinates on an aerial photograph/digital terrain model can be referenced on a global/regional geographic system. Thus, the users are able to determine and locate every point of collected data on the Earth's surface. The typical global positioning system uses the World Geodetic System of 1984(WGS84) datum, and stores the geo-referenced data in GeoTIFF/GeoPDF format. In addition, the users may require orthometric elevation (elevation above sea-level or geoid model) in different scenarios. For example, analysing sea-level change by hydrological data.
Geo-referencing can be done by adding control points at the base of the slope. By shifting the point cloud data onto the known correct co-ordinates with at least 3 points, the point cloud data can be repositioned into an accurate co-ordinate system. This information is useful for calculating distances, volume and areas.

## Types of LiDAR

LiDAR data can be collected on ground-based, airborne and mobile platforms. For example, Airborne LiDAR Scanning (ALS), Unmanned Aerial Vehicles (UAV), Terrestrial Laser Scanning (TLS) and Handheld Laser Scanning (HLS). The tables below will compare the mentioned data collection platforms:
1) The data collection methods
2) Geo-referencing technique (how to acquire the exact coordinate of a point)
3) Advantages and disadvantages
Sample images of different data acquisition methods :
| Airborne LiDAR Scanning (ALS) | Unmanned Aerial Vehicle (UAV) | Terrestrial Laser Scanning | Handheld Laser Scanning |
| An airborne LiDAR scanning system would scan the objects along its flight path with a zigzag pattern | An Unmanned Aerial Vehicle can be controlled by a small scale remote control | A terrestrial laser scanning method that stands on a flat and stable terrain | A handheld laser scanning device with a small size, it is highly portable |
A comparison of different laser scanning methods :
| LiDAR scanning methods | Data collection methods | Geo-referencing technique | Advantages | Disadvantages |
| Airborne LiDAR Scanning (ALS) | An airborne LiDAR system that emits lasers with a wavelength of 1000-1600 nm Data returned from a single pulse with multiple returns for classifying different objects according to canopy structure With a maximum 167 MHz laser repetition rate The distance can be calculated by LiDAR 2 way return time Point density can be more than 8 points/m$^2$ Multi-return with full-waveform 0.1m vertical accuracy 0.3m horizontal accuracy | Composed of an inertial measurement unit (IMU), geospatial global navigation satellite systems(GNSS) for positioning and orientations purposes | Large area coverage (whole territory-wide) The returned pulse can help to classify between the ground and non-ground objects effectively | The flight is limited by air-traffic condition & airspace sharing problem¨ Require take-off and landing sites, an experienced pilot Data loss may occur due to atmospheric absorption Commission error (classifying non-ground points as the ground point) and omission error (rejecting ground points) may occur |
| Unmanned Aerial Vehicle (UAV) | A drone that is integrated with cameras and sensors on its platform for taking ortho-images Point clouds are generated through photogrammetrical workflows, such as structure from motion (SfM) and Multi-view(MVS) techniques by extracting and matching images features Multi-return with full-waveform | A Real-Time Kinematic(RTK) calibration is needed before flying the drone Using geospatial global navigation satellite systems(GNSS) for navigation and geo-referencing Flight direction can be acquired by communicating between web-based map servers and GNSS receivers | Low cost and ready to use within a short period of time No aircraft, airport and experienced workers needed Able to study unreachable areas by air instead of carrying out traditional ground-based surveying technique | The air-traffic condition limits the flight The size of detecting sensors are determined by the payload/ weight capacity of the UAV Commission error and omission error may occur Its battery capacity limits flight duration |
| Terrestrial Laser Scanning (TLS) | A ground-based method that emits 500-600 nm infra-red pulses in the known direction The distance can be calculated by LiDAR 2 way return time About 8mm accuracy | Through setting ground control points with known coordinates | Point cloud data are returned with colours (Red, Green &Blue) | Having the possibility of occlusion due to line of sight problem(Shadow areas issue) Multiple TLS devices have to be positioned on site Only suitable for a site with stable and gentle slope variations |
| Handheld Laser Scanning (HLS) | A ground-based method that can measure micro to mesoscale structure The distance can be calculated by LiDAR 2 way return time 30mm accuracy Single return | Through setting ground control points with known coordinates | High mobility with small device volume Effective in recording site data on an unstable surface, such as a fresh landslide site | Only greyscale coloured points are returned, the point cloud data is difficult for manual filtering Only receive 1 return at a time |

## Digital terrain modelling

A digital terrain model (DTM) is a mathematical model that represents Earth's visible terrain surfaces. A DTM transforms discernible LiDAR data points into a continuous 3-D surface, the model connects discrete points with distinct height values to form planes. Thus, structural geologists are able to derive structural orientations from these 3-D planes. This modelling technique is also used to create digital planetary surfaces and there are more other applications.

### Principle of Digital Terrain Model

The DTMs are classified by the basic geometric units, such as triangles, squares. 3 major mathematic functions are used, as listed below :
* Point-based
* Triangular-based
* Grid-based
Moreover, triangle-based and grid-based are the most used.

#### Point-based

Point-based surface modelling reconstructs a surface by forming a series of small contiguous discontinuous surfaces. Each of the surfaces is planar and formed by connecting the individual data point. This function can form a surface with regular and irregular patterns by considering the regional boundaries of influence of each point, where a Voronoi diagram defines the regional boundaries. However, regular patterns are most widely used for simpler computations, such as hexagons, squared patterns.
The mathematical expression of the formation for each horizontal planar-surface is:
$Z_{i}=H_{i}$
where Z is the height of the planar surface, and H is the ith point's height.

#### Triangular-based

A triangular-based function can form a more tilted or irregular digital surface model. This approach is treated as the primary way to construct a complex DTM. Triangles have great flexibility, any polygons (e.g. a square, rectangle) can be decomposed into other smaller triangles. A linked triangular network can incorporate break lines for plane fitting, it facilitates the formation of a curved facet/surface. The minimum requirement for forming a triangle requires 3 data points, where the nearest 3 points are grouped to form a triangle by Delaunay triangulation without overlapping.
As the data points are unevenly distributed, the triangular-based method effectively constructs DTM because this method can create surface variations. Even some data points of the point cloud are removed or added, a local triangle can be reformed without complete DTM reconstruction.
A few parameters control the surface formation process:
1) Density of Mesh
2) Maximum angle of neighbouring triangles
3) Minimum Patch size

#### Grid-based

The grid-based model has less relevance for constructing broken terrain or sharp terrain discontinuities with steep slopes. 4 data points are the minimum requirement for forming a grid-based model. The resulting surface is named a bi-linear surface that quadrilaterals of any shapes are linked to create a DTM, such as parallelograms, squares, rectangles or other irregular polygons. This method has the advantage of data handling that the data is in the form of squared grids, evenly distributing data points. In this case, some software has provided a "Random-to-Grid" operation to ensure that the data has the grid form.

## Approach for digital terrain modelling

As the data collected are in the form of a point cloud, there is a need to transform these 3-D co-ordinates of laser points into a digital terrain model with 2 major procedures:
1) Point cloud classification and ground filtering
2) Reconstruction of the ground surface from discrete laser point cloud data by interpolation.

### Data processing and classification (Airborne LiDAR scanning only)

Data classification and noise cleaning are the processes of obtaining a non-biased slope surface. When ALS collects the with multiple returns, this principle can classify the objects into different categories. Classification algorithms can be performed by TerraScan & TerraModel, computer software for classifying point cloud data automatically developed by TerraSolid. However, some manual adjustments and validations are needed to ensure the data points are classified correctly.
Algorithms can identify pre-dominate landform features, these algorithms assume the surfaces with significant variations are non-ground features. For creating a surface model of rock slope, only classified ground data is needed. The classifications of the ground and non-ground surfaces are in the following categories:
Major LiDAR Data Classification Categories for Airborne LiDAR Scanning :
| Classifications | Examples |
| Flat, ground or open terrain | Open terrain, sand, rock, flat rooftops, golf courses |
| Low-vegetation area | Tall crops and weeds <2m |
| Medium-vegetation area | Low trees or brushlands between 2 and 10m |
| High-vegetation area | Forested areas >10m |

### Data filtering

Data filtering enables the extraction of bare-Earth's surface by removing unnecessary data or data noise. For rock orientation studies, outcrop mapping and topographical studies, they only acquire information about the rock slope bodies. Thus, data filtering is implemented to separate the point clouds as ground and non-ground features. Bugs, vegetations or other artificial infrastructures belong to non-ground features.
In order to filter 3-D laser points, few methods from open source software can be applied to filter ground points out of the whole area of interests :
A comparison of different data filtering methods :
| Author | Technique | Description of the method | Accuracy |
| Axelsson et al., 2000 | Triangular-irregular network (TIN) | Multiple points are divided into grids. The lower elevation within each grid is treated as ground points. Users can set the grid size on their own. The points are added into a triangle facet when they are close enough. | Satisfactory |
| Zeybek, 2019 | Cloth-simulation method | A simulated cloth is put on the up-side-down point cloud, points are not filtered when the cloth can touch the data points. The softer the cloth, the more intersection points can be preserved. The intersected points are then classified as a ground. | 93% |
| Zeybek, 2019 | Surface-based | This method uses a linear algorithm to remove the non-ground points | 81% |
| Zeybek, 2019 | Curvature based | Fitting is applied to the nearest neighbour, the points with large elevation variations(Z component) are filtered out. | 88% |

### Ground verification survey

The purpose of the ground verification survey is to test the accuracy of LiDAR. LiDAR data is acquired by sending laser pulses at different angles or receiving the returning signal, these signals may include some error induced by atmosphere absorption of wavelength. Therefore, a ground-truth survey is needed to ensure the collected data's co-ordinates match the local coordinate system. For example, horizontal accuracy will be tested by comparing the collected data through different data collection techniques. Also, data can be corrected by setting multiple control points with known coordinates.

## Structural orientation analysis

The DTMs are capable of identifying the structural parameters of geological features(tilt angle of a fold limb). For example, dip and dip directions of rock planes can be derived from DTM. The general methodology has the following steps:
1) Obtain LiDAR data of a targeted slope
2) Noise cleaning and filtering of point data
3) Transform point data into DTM through mathematical functions
4) By using programs, such as Coltop 3D, the structural geometry of the joint/slope planes can be calculated automatically.
5) The rock planes are determined by color-coding or statistical method for evaluating whether the points are grouped to form a single or different planes.
6) Output the structural orientations data into suitable computer software files, and plotting the distribution of dip and dip directions of different planes on a stereonet.

### Data partition

Data partitioning helps redistribute the unevenly distributed point data by dividing data into cubes. Returned point cloud data has different point density, due to the variation of data collection devices and parameters. The point cloud's density determines the cube size, no fewer than 4 points form each cube. The points that are considered as non-conformers will be removed from the point cloud. Afterwards, the normal vector of each cube will be calculated.
Octree Partitioning in open source softwares, including CloudCompare and Geomagic, can achieve data partitioning. Considering the rock masses in different terrains, they have variations in the roughness. Thus, it is necessary for users to set the cube size manually to obtain the best cube. The average number of points in each cube ranges from 15 to 30 points, by setting the range of the point spacing between 4mm and 7mm.

### Cluster analysis by normal vectors

The purpose of rock discontinuities clustering is to group sub-planes of a slope into the same discontinuity set they belong to. The discontinuities demonstrate waviness, roughness or an undulating surface. The points within the same joint set show similar orientation habits. The algorithm determines whether the points are located within the main orientation of bedding planes or sub-parallel planes with each other. The clustering technique can classify different joint sets based on assigned normal vectors on each face within approximately the same orientation.

## Challenges for LiDAR Technology

Although LiDAR has high efficiency in collecting data in a large area within a short time, some issues remain challenging for the overall data processing and generating the expected result.

### Data filtering and surface generation parameters

Vegetation may not be entirely removed during data filtering. It may affect the smoothness and the clustering of rock surfaces. In most cases, 3-D models are generated through the triangulation method, spikes are often formed. The spikes may affect the clustering of rock surfaces. They affect the smoothness of the 3-D surface, inducing errors in calculating the rock plane orientations. Spikes are formed when the data points have similar X and Y coordinates but very different Z coordinates, the needle-like shape triangles are produced as a result of non-smooth surfaces.
The filtering parameters for surface clustering often depend on users' past experiences. For example, the octree parameter can be adjusted with different point densities. The result of DTM formation has to be determined by users. Hence, a repeated testing procedure is needed for generating a satisfactory surface.

### Point density

Point cloud density, known as the spacing between each data point in an obtained LiDAR datasets. This parameter affects the accuracy in measuring rock slopes. It is one of the characteristics that need to be considered during data processing. For example, data filtering, classification, feature extraction and object recognition. Point cloud density depends on various factors:
* Sight overlapping of TLS
* Reflective properties of the object of interest
* Scanning frequency of ALS

### Cost

Although LiDAR has been an effective data acquisition method on a rock slope, its high cost caused the technique to be not practical. When the area of interest is of a very small scale, it limits the usability of TLS & ALS, given carrying out an ALS survey requires aircraft, an experienced pilot, a designed flying path and height that got the approval of the local aviation department.
By employing UAV, the problem mentioned above can be resolved. It can collect data in inaccessible areas that enable data collection at a low cost at a small scale, since the device is portable and lightweight.

## Other LiDAR applications

With the high accuracy and high-efficiency data acquisition method of laser scanning, it can potentially apply to different areas other than structural measurements:
* Emergency inspections and landslide mapping
* Project management (with the need of 3-D models / topographic data)
* Site progress, the volume of excavation or filling
* Volumetric calculations
* Topographic mapping
* Building and highway construction
* Natural disasters assessments
* Unreachable area data collection
* Precision farming
* As-built construction records
* Unauthorised works (e.g. capturing of field situation)
* Monitoring surface deformation with continuous and repeated LiDAR scanning
* Archaeology and site reconstructions
* Seismology
* Forestry
* Atmospheric physics
* Extraterrestrial explorations