Next Article in Journal
AIS and VBD Data Fusion for Marine Fishing Intensity Mapping and Analysis in the Northern Part of the South China Sea
Next Article in Special Issue
Methodological Proposal for Automated Detection of the Wildland–Urban Interface: Application to the Metropolitan Regions of Madrid and Barcelona
Previous Article in Journal
3D Visibility Analysis for Evaluating the Attractiveness of Tourism Routes Computed from Social Media Photos
Previous Article in Special Issue
Maximizing Impacts of Remote Sensing Surveys in Slope Stability—A Novel Method to Incorporate Discontinuities into Machine Learning Landslide Prediction
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exploiting a Semi-Automatic Point Cloud Segmentation Method to Improve the Quality of Rock-Mass Characterization. The Cima Grappa Conservative Restoration Case Study

1
Department of Civil and Environmental Engineering, University of Florence, Via di S. Marta 3, 50139 Firenze, Italy
2
Geoapp s.r.l., Viale Spartaco Lavagnini 70/72, 50129 Firenze, Italy
*
Author to whom correspondence should be addressed.
ISPRS Int. J. Geo-Inf. 2021, 10(5), 276; https://doi.org/10.3390/ijgi10050276
Submission received: 10 February 2021 / Revised: 13 April 2021 / Accepted: 26 April 2021 / Published: 28 April 2021
(This article belongs to the Special Issue Advancements in Remote Sensing Derived Point Cloud Processing)

Abstract

:
This paper presents results from applying semi-automatic point cloud segmentation methods in the underground tunnels within the Military Shrine’s conservative restoration project in Cima Grappa (Italy). The studied area, which has a predominant underground development distributed in a network of tunnels, is characterized by diffuse rock collapsing. In such a context, carrying out surveys and other technical operations are dangerous activities. Considering safety restrictions and unreachable impervious tunnels, having approached the study area with the scan-line survey technique resulted in only partial rock mass characterization. Hence, the geo-mechanical dataset was integrated, applying a semi-automatic segmentation method to the point clouds acquired through terrestrial laser scanning (TLS). The combined approach allowed for remote performance of detailed rock mass characterization, even remotely, in a short time and with a limited operators presence on site. Moreover, it permitted extending assessing tunnels’ stability and state of conservation to the inaccessible areas.

1. Introduction

The Cima Grappa Military Shrine restoration project has been promoted and financed by the Italian Ministry of Defence of the Italian Presidency of the Council of the Ministers. It concerns a mountain area on the top of Monte Grappa (also called “Cima del Grappa”) in the Venetian Prealps context (Figure 1).
The Cima Grappa Military Shrine is the largest Italian military ossuary of the First World War. During the First World War, many Italian soldiers died. Most of them were buried in mass graves close to battlefields. Between 1920 and 1930, remains of the fallen were unearthed and reburied in a large war memorial. Politicizing the memory of the fallen was one of the primary purposes. The war memorials, designed by architects close to the regime, were constituted near former front lines and kept under the patronage of a special commission from the former Ministry of War [1].
The hereafter presented activities that insisted on the underground connection area between the Milan Barracks and the former underground infirmary, focus on securing the former infirmary rooms and tunnels. These areas were dug in a fractured limestone rock mass. The spaces are cramped, and the accessible areas are tunnels that span from 1.2 to 2.2 m in height and 0.9 to 2.4 m in width. The overall study area’s extension is 300 m2 with a tunnel’s lengths of about 100 linear meters (Figure 2). During past years, several collapses have occurred in the underground infirmary and the access tunnels. Thanks to the restoration project, the site will become a new tourist area. A detailed geo-mechanical characterization of the former infirmary’s underground area has been necessary to identify and design the optimal stabilization and restoration interventions [2,3]. Discontinuities play a crucial role in determining rock mass stability [4,5], and we need to investigate extensive rock-mass portions [6,7,8,9,10] to perform sound statistical analyses. The manual surveying of fractures and planar rock mass surfaces is one of the most fundamental but time-consuming activities performed by surveyors.
The classical geo-mechanical survey [11] on its own would not have brought sound results [12]. Primarily because the tunnels’ roof was too high to be reached with classic instrumentation; and secondly, because the amount of data acquired through the classical scan-line survey would not have a good statistical significance, considering the geometrical complexity of the tunnels. Finally, the relatively poor conservative status of the underground area and some inaccessible spaces would not have allowed the technicians to access all the spaces. Furthermore, the restoration project has a strong conservative purpose and using non-destructive investigation methods is a compulsory request. The laser scanning technique was applied to collect high-density point clouds to identify the discontinuities main sets [13,14,15,16,17,18,19]. The geo-mechanics surveying community have been extracting fractures and planar surfaces orientations from 3D outcrop point clouds through semi-automatic and automatic methods developed in the last 15 years.
Using random sample consensus (RANSAC) algorithm-based methods to segment point clouds into subsets [20,21,22], deriving triangulated irregular networks from point cloud and grouping neighbouring polygons with a similar orientation to obtain planar features [23,24], and some other methods based on k-means clustering [25], moving sampling cube, point attributes, neighbouring points coplanarity testing, and principal component analysis (PCA), were developed and can be applied to gather rock mass geometrical settings instead of performing a classical geo-mechanic survey. Considering the narrow and complicated geometrical setting of the study area, a method designed with supervision capability was chosen and applied within the presented paper. Finally, kinematic analyses were carried out by applying 3D methodology [26,27,28,29].

2. Materials and Methods

In correspondence to inaccessible areas and areas prone to collapsing events, the classic geo-mechanical survey (Figure 3) scan-line detection method [30] of the rock mass was substituted with measures obtained from laser scanning datasets [31,32]. The geo-mechanical analysis results were combined with the analysis of the point clouds acquired by the laser scanner. In particular, a semi-automatic extraction of rock mass structural data from high-resolution laser point clouds was adopted [33,34,35]. This section first introduces the classical geo-mechanical scan-line method (Figure 3) and then uses segmentation procedures to extract structural data from point clouds.

2.1. Classic Manually Performed Geo-Mechanical Scan-Line Survey

A geo-mechanical survey of the calcareous rock mass surfacing in the infirmary rooms was performed. It consisted of the following steps, (a) geo-mechanical characterization of intact rock, (b) geo-mechanical characterization of discontinuities, and (c) geo-mechanical characterization of the rock mass. Skipping a detailed description of a geo-mechanical survey [36] beyond the present paper’s scope, the essential parameters to describe joints, cracks, and discontinuities are reported. In particular, cracks orientation, number of families of discontinuity, spacing, persistence, opening, roughness, surface resistance, degradation, filling, and filtration [37].

2.2. Remote Geo-Mechanical Survey

Using terrestrial laser scanning (TLS) for rock mass discontinuity characterization [38,39,40,41] was reported in several scientific works [42,43]. The DiAna3D algorithm is applied to extend the rock-mass geo-mechanical assessment to the less accessible areas. The method is more rapid than the standard scan line survey, and it can minimize the operator’s exposure to rock collapse risk. Thanks to the laser scanning data high spatial resolution, it is possible to use point clouds to integrate the traditional survey data and obtain information on the rock mass geometry remotely if compared to the manual magnetic compass method.
In the specific case, a calculation algorithm was applied to the point clouds to extract accurate geometric information, valid for the rock mass reconstruction and characterization. In particular, through DiAna, the discontinuity plans were identified semi-automatically, and their spatial orientations were calculated.
The analysis was carried out in correspondence with the limestone rock mass’s main outcrop areas and extracted 2761 discontinuous surfaces. However, a surface of discontinuity that is not perfectly planar and has a high persistence may be fragmented into a certain number of "sub-surfaces". The number of discontinuities extracted is overestimated.
The laser scanning technique was used to generate point clouds and measure geometrical parameters as persistence, spacing, roughness, and discontinuities orientation.

2.3. Acquisition Campaign

The laser scanning acquisition campaign involved the underground rooms of the infirmary behind the Milan barracks. The used laser scanner model is a Riegl LMS-Z420i ground-based (Table 1).
A series of high-resolution laser scans were carried out from 13 different scanning positions starting from the square in front of the Milan barracks to mapping the underground hospital area. The obtained point clouds were then georeferenced by aligning them to that acquired externally. The scans were carried out mainly with the vertical instrument setup; a horizontal configuration was also used in correspondence with the major collapses to allow the complete acquisition of the tunnel roof. For each scanning position, low-resolution framing scans and detail scans were acquired. A total of over 40 million points were thus acquired, with the average angular resolution of 0.056°, a minimum acquisition distance of 1.5 m, and a maximum acquisition distance of 20 meters. One point/cm2 is the average point cloud density of the processed point clouds after optimization and unification. Before proceeding with the laser scan, some targets/reflectors were positioned within the scene (Figure 4).
Up to 15 cylindrical reflectors (retroreflector targets) with a pre-determined size of 50 × 50 mm were deployed, performing each scanning acquisitions. Reflectors’ use is linked to the possibility of the acquisition software recognizing them in the various scans and combining them into a single point cloud, minimizing the problems related to shaded areas. Thanks to the reflector’s high reflectivity, their positions were measured with the highest TLS resolution.
The TLS scans registration accuracy was improved by applying the multi-station adjustment algorithm (MSA). It was developed by the laser manufacturer and based on the recognition of homologous plans materialized for each scan (which therefore require a scan overlap between them) and their iterative alignment, up to minimizing the position differences (Figure 5).
The overlapping value varied from 35 to 60%. Further, 60% could be considered a relatively high value. It is the partial consequence of the very articulated tunnels’ morphology, characterized by many bumps, depressions, and irregularities of the rock mass surface. Indeed, the consequent massive number of shaded areas led to complicated scanning geometries and made designing scan positions effectively a problematic task.
Cross-sections of the infirmary point cloud are reported in Figure 6 to better understanding the former infirmary’s spaces and shapes.

2.4. Point Clouds Segmentation for Remote Geo-Mechanical Survey

As mentioned in the section above, a semi-automatic segmentation method, DiAna3D, was used to extract geometric characteristics of discontinuities and relative planes. By processing point clouds, some of the rock-mass parameters [44,45] are block size, the number of sets, persistence, orientation (Figure 7), spacing/frequency (and derived RQD), and scale-dependent roughness can be determined.
The three-dimensional approach of the segmentation method is hereafter shortly reported. The first step is aimed at identifying the principal planes on the cloud. A searching cube contains a selection of a point cloud subset. The best-fitting plane is identified for the searching cube volume by applying the least square method. The plane is described with the general equation (Equation (1)) end its direction cosines characterize it l , m , n :
a x + b y + c z + d = 0
The algorithm computes the selected points’ standard deviation for the best fitting plane. Simultaneously, the operator can designate a threshold value (st) depending on the point cloud resolution, roughness of the rock slope, and searching cube dimension. A cluster is formed with the points extracted from the point cloud when s < st. A box analysis can be performed; in this case, the searching cube is moved along the geographic axes according to a regular pattern to process the entire point cloud. A discrete analysis can also be processed; in this case, the searching cube is centred at every point of the point cloud. For each valid determined cluster, the correlate plane orientation is calculated by applying Equations (2) and (3):
α = a r c t a n m l + Q
β = a r c t a n n l 2 + m 2
where α and β are the dip and dip/direction of the plan; l , m , n are the direction cosines of the plane; Q is a constant which depend on l , m .
Although the method is quite far from moving with a compass and clinometer on the rock face, taking measurements on discontinuity planes allows local face orientation measurements.
The operator can choose a variation on cube size, and it returns details that can be observed on the point cloud.
The extraction process of the discontinuity planes consisted of 3 main steps: In the first one, organized point clouds were analyzed (Figure 8), looking for the points around which the local geometry is sufficiently planar (standard deviation of the points concerning the interpolation plane, below a user-defined threshold (blue areas in Figure 9a; starting from these points, step 2 consisted of a progressive germination algorithm that allowed to gradually add all the neighbouring points, co-planar with the initial nucleus (Figure 9b), and to define its external 3D polygon (Figure 9c,d). The last step allowed to combine the results obtained from the individual scan positions by eliminating redundancies in the overlapping areas. For each identified plane, its persistence was calculated by measuring the maximum linear distance between two points of the cluster of points that compose it. The 3D polygons that delimit the discontinuity planes with persistence greater than 10 cm; thus, identified are shown in Figure 10a,b.
Once cluster orientations are determined (Figure 11), they are plotted on stereographic projections (Figure 12).
In this way, defining a discontinuity set as a series of discontinuities, for which orientation is broadly the same, discontinuity sets can be easily identified and extracted from the point cloud. A common ID is assigned to the clusters belonging to the same set. This last-mentioned process is manually driven to have more substantial supervision of the set recognition procedure and guarantees a real-time integration with possible validated in-field data.
By family of discontinuities, we mean the set of discontinuities with a similar orientation (sub-parallels), generally of the same type and, presumably, originating with the same genetic mechanism.
Four families of main discontinuities were recognized (Figure 12), whose respective modal planes’ orientations (Table 2).

2.5. 3D Kinematic Stability Analysis of the Rock Mass

After doing the geo-mechanical characterization of the intact rock, the discontinuity system, and the rock mass, the data obtained were used to carry out a 3D kinematic analysis of the rock mass [46]. The stability of the blocks composing a rock mass is essentially controlled by the discontinuities’ geometric and geo-mechanical properties. The term “kinematic analysis” refers to analyzing blocks of the rock’s behaviour under their weight and with purely attractive sliding resistance along the discontinuity surfaces [47]. The analysis allows investigating the following instability mechanisms (Figure 13): plane failure (PF), wedge failure (WF), block toppling (BT), flexural toppling (FT), and free fall (FF).
The kinematic analysis of slope stability is a valuable tool for understanding rock mass behaviour. It represents a valid connection point between the technical geological survey and the design phase. The analysis is generally carried out by graphically verifying certain geometric conditions on the lower hemispheric projection. Equiangular projection is usually used as it allows to maintain the relations between the angles.
The basic assumptions for kinematic analysis are:
  • the rock mass is divided into blocks by flat and infinitely persistent discontinuities;
  • the shear resistance along the discontinuity planes is purely due to friction;
  • the block system is subject to weight force only.
The term “kinematic analysis” refers to analyzing individual blocks of the rock’s behaviour under its weight and with purely attractive sliding resistance on discontinuous surfaces.
This type of analysis applied to the slopes constitutes a valuable tool for understanding rock mass behaviour. They allow for a rapid assessment of the propensity to failure for the various breaking mechanisms in the sectors of the slopes analyzed.
The kinematic analysis of the different failure mechanisms can be performed using geometric relationships amongst discontinuity’s planes, intersection lines between different discontinuity planes, and the excavation face plane.
These can be expressed as vector relations between the vector units that identify these planes, expressed in polar coordinates using the azimuth and the inclination on the horizontal β .
In particular, the planes of discontinuity are identified by the versors of the lines of maximum slope p   α p ; β p , or their normal n α n ; β n or poles;
The intersection lines between pairs of plans are identified by versors i   α i ; β i .
While the excavation line is identified by the versor of the maximum slope line f   α f ; β f .
By convention, all versors units are considered with the extreme in the lower half-space of the horizontal plane.
A quantitative approach can be pursued by defining a kinematic hazard index for each kinematic mechanism, which expresses the relative probability that a given mechanism failure may occur based on the structural and geo-mechanical conditions of the rock mass. Such indexes can be calculated by counting the following quantities:
  • Npf = number of poles of the discontinuities that satisfy the conditions for planar sliding;
  • Iwf = number of the intersections of the discontinuities that satisfy the conditions for the sliding of wedges;
  • Nbt = number of the poles of the discontinuities that satisfy the conditions for direct overturning;
  • Ibt = number of intersection lines that satisfy the conditions for direct overturning;
  • Nft = number of the poles of the discontinuities that satisfy the conditions for bending overturning.
Indicating with N the total number of discontinuities sampled in the domain and with I = 0.5 (N2-N), the number of all possible intersections, it is possible to calculate the hazard kinematics indexes reported in Table 3.

2.6. Extraction of Detailed Meshes

The overall point cloud, obtained by combining the information acquired from each scanning position, was triangulated, getting a continuous surface from discrete data of a high point density.
The resulting 3D model is shown in Figure 14.
From the detailed mesh obtained by applying the Poisson surface reconstruction method [49], it was then possible to extract contour lines with an equidistance of 10 cm, as shown in Figure 15.

3. Results

Thanks to the performed survey campaigns, a considerable amount of the 3D data of the underground spaces was obtained. A series of elaborations were processed, capitalizing on the high-resolution point clouds acquired by laser scanning.
Works were explicitly aimed at:
(a)
Extracting a digital 3D model of the limestone rock-mass;
(b)
Integrating the traditional geo-mechanical survey by identifying the main discontinuities affecting the rock mass and dividing them into families;
(c)
Providing detailed geometric reconstruction of the collapse events that occurred;
(d)
Identifying main failure mechanisms of the rock mass in correspondence of the tunnels and underground spaces (kinematic analysis).
The high-resolution point clouds acquired with the laser scanner were also used to reconstruct the position and volumes of the main collapse phenomena that have affected the underground rooms under investigation over the years. In particular, six main events were identified, the location of which is shown in Figure 16.
The corresponding volumes were calculated using laser scanner data for each collapse, considering the detachment niches and accumulations (Figure 17).
In particular, collapse 1 (0.4 m3), collapse 2 (1.2 m3), collapse 3 (1.7 m3), collapse 4 (6.0 m3), collapse 5 (1.5 m3), collapse 6 (7.0 m3).
Figure 18 shows the results of the graphic processing relating to the quantitative kinematic analysis of the rock mass understudy, as described in Section 2.3. The most probable failure mechanisms are wedge failure and, subordinately, planar sliding by observing above mentioned figures. The maximum values of the hazard indices reach 35% for the individual mechanisms. These data testify to a high propensity for instability of the investigated tunnels and allow us to identify which sectors of the rock mass are most likely to generate phenomena of kinematic instability.

4. Discussions

By applying point clouds segmentation, remotely extracted orientations of planes were measured. Orientations of joint families modal plans were identified (BG: 13°/245°, JN1: 76°/069°, JN2: 76°/001°, JN3: 75°/125°) and a stereographic projection was populated as presented in Figure 12.
From the comparison between the stereographic projections (Figure 14 and Table 2), a good correspondence between the joint families measured by scan line survey and the ones extracted applying the semi-automatic method is noticeable. The discrepancy stays within 10% for the studied areas. Furthermore, a good correspondence between the main collapse and main joint families is evident, as Figure 16 shows.
The definition of pre-assigned size value to the searching cube of the DiAna3D algorithm to efficiently perform the semi-automatic family extraction at a different scale and with a different desired detail could be an alternative solution to implement.
There are many remote sensing applications to support restoration projects of cultural heritage sites, including laser scanning and photogrammetry [50,51,52,53,54,55]. Performing considerable quantities of measurements, including inaccessible areas with fewer efforts than manual measurements, is one of the main advantages. It allows for drastically reducing the time it takes to perform the survey and enables a broader and more detailed view of the context. Moreover, the reduction of operators’ time exposure to risk-prone areas and the capability to acquire valuable data from inaccessible spaces completes the benefits list.
However, to survey the whole area, a group of 3 operators spent more than three days. Considering the relatively small extension, the long time spent appears to have resulted from the highly complex surface’s morphology and the cramped spaces. Although static TLS is an appropriate technique for this kind of environment, a mobile and presumably remoted operating system, e.g., a TLS mounted on a robotized platform, would have performed the same in less time. Investigating results obtained by applying other segmentation methods [56] to the same point cloud dataset could be further researched to compare different performances and identify the most suitable method.

5. Conclusions

A laser scanning survey was effectively exploited to remotely perform rock mass characterization in the presented context, a conservative restoration project for a high-value cultural heritage underground area.
As the studied area has inaccessible portions, two different types of survey were carried out, the traditional manual geo-mechanical survey and the remote geo-mechanical survey.
Main joint families were extracted by applying the semi-automatic classification algorithm DiAna3D. A dataset of the whole investigation area’s joints was defined by integrating results from the traditional geo-mechanical survey and semi-automatic classification.
The acquired high-resolution point clouds have made it possible to extract the discontinuities that pervade the rock mass dividing it into small blocks. This process, combining the results of traditional geo-mechanical surveys, has allowed us to characterize the rock mass from a quantitative perspective.
The traditional geo-mechanical survey effectiveness (scan-line survey) is minimal since the investigated area is inaccessible or unreachable. As in the proposed case, it is very articulated and intensely fractured.
Performing a remote geo-mechanical survey was allowed to carry out a high-quality geo-mechanical rock-mass characterization quickly compared to the traditional method and maintain the survey activities under an acceptable level of risk exposure.
Comparing results from traditional scan line surveys and the remote one, in terms of joint family orientations, points out the remote method’s good reliability, even in narrow and unconventional environments as relatively small underground tunnels.

Author Contributions

Conceptualization, Francesco Mugnai; methodology, Francesco Mugnai; software, Francesco Mugnai and Paolo Farina; validation, Francesco Mugnai, Paolo Farina and Grazia Tucci; formal analysis, Francesco Mugnai; investigation, Francesco Mugnai and Paolo Farina; resources, Paolo Farina; data curation, Francesco Mugnai; writing—original draft preparation, Francesco Mugnai; writing—review and editing, Francesco Mugnai, Paolo Farina and Grazia Tucci; visualization, Francesco Mugnai; supervision, Francesco Mugnai, Grazia Tucci; project administration, Paolo Farina; funding acquisition, Paolo Farina. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

3rd Party Data. Restrictions apply to the availability of these data. Data was obtained from Politecnica Ingegneria e Architettura, and are available from the authors upon permission given by Politecnica Ingegneria e Architettura.

Acknowledgments

The authors want to thank: (1) Italian Ministry of Defence and the Ministry Council’s Presidency for having founded this study, (2) Geol. Pietro Accolti of Politecnica and Architettura, the leading promoting company of the temporary grouping of professionals (consortium) for designing, together with: (3) ACPV studio, (4) Parallelo and Bit Movies, (5) Architect Gianfranco Tedeschi as Project Manager, (6) Engineer Fabio Camorani as Structural and Geotechnical design, and (7) engineer Francesco Frassineti, (8) Massimo Cavazzuti and (9) Claudio Pongolini for the technical project harmonization. (10) GEOAPP s.r.l. for technical support and outstanding specialized experience in managing such a sensitive and complex survey project.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Malone, H. Fascist Italy’s Ossuaries of the First World War: Objects or symbols? RIHA J. 2017, 0166. Available online: http://www.riha-journal.org/articles/2017/0150-0176-special-issue-war-graves/0166-malone (accessed on 10 February 2020).
  2. Barton, N.R.; Lien, R.; Lunde, J. Engineering classification of rock masses for the design of tunnel support. Rock Mech. 1974, 6, 189–239. [Google Scholar] [CrossRef]
  3. Bieniawski, Z.T. Engineering Rock Mass Classification; John Wiley and Sons: New York, NY, USA, 1989; 251p. [Google Scholar]
  4. Jaboyedoff, M.; Metzger, R.; Oppikofer, T.; Couture, R.; Derron, M.-H.; Locat, J.; Turmel, D. New insight techniques to analyze rock-slope relief using DEM and 3Dimaging cloud points: COLTOP-3D. In Proceedings of the 1st Canada-U.S. Rock Mechanics Symposium, American Rock Mechanics Association (ARMA), Vancouver, BC, Canada, 27–31 May 2007. [Google Scholar] [CrossRef]
  5. Li, X.; Chen, Z.; Chen, J.; Zhu, H. Automatic characterization of rock mass discontinuities using 3D point clouds. Eng. Geol. 2019, 259, 105131. [Google Scholar] [CrossRef]
  6. Riquelme, A.; Tomás, R.; Cano, M.; Pastor, J.L.; Abellán, A. Automatic mapping of discontinuity persistence on rock masses using 3D point clouds. Rock Mech. Rock Eng. 2018, 51, 3005–3028. [Google Scholar] [CrossRef]
  7. Bolkas, D.; Vazaios, I.; Peidou, A.; Vlachopoulos, N. Detection of rock discontinuity traces using terrestrial LiDAR data and space-frequency transforms. Geotech. Geol. Eng. 2018, 36, 1745–1765. [Google Scholar] [CrossRef]
  8. Feng, Q.; Röshoff, K. A survey of 3D laser scanning techniques for application to rock mechanics and rock engineering. In The ISRM Suggested Methods for Rock Characterization, Testing and Monitoring: 2007–2014; Springer: Cham, Switzerland, 2015; pp. 265–293. [Google Scholar]
  9. Gigli, G.; Casagli, N. Semi-automatic extraction of rock mass structural data from high-resolution LIDAR point clouds. Int. J. Rock Mech. Min. Sci. 2011, 48, 187–198. [Google Scholar] [CrossRef]
  10. Gomes, R.K.; De Oliveira, L.P.L.; Gonzaga, L.; Tognoli, F.M.W.; Veronez, M.R.; De Souza, M.K. An algorithm for automatic detection and orientation estimation of planar structures in LiDAR-scanned outcrops. Comput. Geosci. 2016, 90, 170–178. [Google Scholar] [CrossRef]
  11. Admassu, Y.; Shakoor, A. Cut slope design recommendations for sub-horizontal sedimentary rock units in Ohio, USA. Geotech. Geol. Eng. 2013, 31, 1207–1219. [Google Scholar] [CrossRef]
  12. Battulwar, R.; Zare-Naghadehi, M.; Emami, E.; Sattarvand, J. A state-of-the-art review of automated extraction of rock mass discontinuity characteristics using three-dimensional surface models. J. Rock Mech. Geotech. Eng. 2021, 1674–7755. [Google Scholar] [CrossRef]
  13. Guo, J.; Wu, L.; Zhang, M.; Liu, S.; Sun, X. Towards automatic discontinuity trace extraction from rock mass point cloud without triangulation. Int. J. Rock Mech. Min. Sci. 2018, 112, 226–237. [Google Scholar] [CrossRef]
  14. Henk, A.; Laux, D. Terrestrial laser scanning and fracture network characterization-perspectives for a (semi-) automatic analysis of point cloud data from outcrops. Z. Dt. Ges. Geowiss. Ger. J. Geol. 2015, 166, 99–118. [Google Scholar]
  15. Hu, L.; Xiao, J.; Wang, Y. Efficient and automatic plane detection approach for 3-D rock mass point clouds. Multimed. Tools Appl. 2019, 79, 839–864. [Google Scholar] [CrossRef]
  16. Kong, D.; Wu, F.; Saroglou, C. Automatic identification and characterization of discontinuities in rock masses from 3D point clouds. Eng. Geol. 2020, 265, 105442. [Google Scholar] [CrossRef]
  17. Zhang, P.; Li, J.; Yang, X.; Zhu, H. Semi-automatic extraction of rock discontinuities from point clouds using the ISODATA clustering algorithm and deviation from mean elevation. Int. J. Rock Mech. Min. Sci. 2018, 110, 76–87. [Google Scholar] [CrossRef]
  18. Riquelme, A.J.; Abellán, A.; Tomás, R.; Jaboyedoff, M. A new approach for semi-automatic rock mass joints recognition from 3D point clouds. Comput. Geosci. 2014, 68, 38–52. [Google Scholar] [CrossRef] [Green Version]
  19. Cao, T.; Xiao, A.; Wu, L.; Mao, L. Automatic fracture detection based on terrestrial laser scanning data: A new method and case study. Comput. Geosci. 2017, 106, 209–216. [Google Scholar] [CrossRef]
  20. Roncella, R.; Forlani, G. Extraction of planar patches from point clouds to retrieve dip and dip direction of rock discontinuities. In Proceedings of the ISPRS WG III/3, III/4, V/3 Workshop “Laser Scanning 2005”, Enschede, The Netherlands, 12–14 September 2005; pp. 162–167. [Google Scholar]
  21. Voyat, I.; Roncella, R.; Forlani, G.; Ferrero, A.M. Advanced techniques for geo structural surveys in modelling fractured rock masses: Application to two Alpine sites. In Proceedings of the Golden Rocks 2006, the 41st US Symposium on Rock Mechanics (USRMS), American Rock Mechanics Association, Golden, CO, USA, 17–21 June 2006. [Google Scholar]
  22. Ferrero, A.M.; Forlani, G.; Roncella, R.; Voyat, H.I. Advanced geostructural survey methods applied to rock mass characterization. Rock. Mech. Rock. Eng. 2009, 42, 631–665. [Google Scholar] [CrossRef]
  23. Lato, M.; Diederichs, M.S.; Hutchinson, D.J.; Harrap, R. Optimisation of LiDAR scanning and processing for automated structural evaluation of discontinuities in rockmasses. Int. J. Rock Mech. Min. Sci. 2009, 46, 194–199. [Google Scholar] [CrossRef]
  24. Slob, S.; Hack, R. 3D terrestrial laser scanning as a new field measurement and monitoring technique. In Engineering Geology for Infrastructure Planning in Europe. Lecture Notes in Earth Sciences; Springer: Berlin/Heidelberg, Germany, 2004; pp. 179–189. [Google Scholar]
  25. García-Sellés, D.; Falivene, O.; Arbués, P.; Gratacos, O.; Tavani, S.; Muñoz, J.A. Supervised identification and reconstruction of near-planar geological surfaces from terrestrial laser scanning. Comput. Geosci. 2011, 37, 1584–1594. [Google Scholar] [CrossRef]
  26. Gigli, G.; Casagli, N. Extraction of rock mass structural data from high resolution laser scanning products. In Landslide Science and Practice: Spatial Analysis and Modelling; Springer: Berlin/Heidelberg, Germany, 2013; pp. 89–94. [Google Scholar]
  27. Van Knapen, B.; Slob, S. Identification and characterization of rock mass discontinuity sets using 3D laser scanning. Procedia Eng. 2006, 191, 838–845. [Google Scholar]
  28. Weidner, L.; Walton, G.; Kromer, R. Classification methods for point clouds in rock slope monitoring: A novel machine learning approach and comparative analysis. Eng. Geol. 2019, 263, 105326. [Google Scholar] [CrossRef]
  29. Hoek, E.; Bray, J.W. Rock Slope Engineering, 3rd ed.; The Institution of Mining and Metallurgy: London, UK, 1981; pp. 341–351. [Google Scholar]
  30. Priest, S.D. Discontinuity Analysis for Rock Engineering; Chapman and Hall: London, UK, 1993. [Google Scholar]
  31. Olariu, M.I.; Ferguson, J.F.; Aiken, C.L.V.; Xu, X. Outcrop fracture characterization using terrestrial laser scanners: Deep-water Jackfork sandstone at big rock quarry. Ark. Geosphere 2008, 4, 247–259. [Google Scholar] [CrossRef]
  32. Abellán, A.; Vilaplana, J.M.; Martínez, J. Application of a long-range terrestrial laser scanner to a detailed rockfall study at Vall de Núria (Eastern Pyrenees, Spain). Eng. Geol. 2006, 88, 136–148. [Google Scholar] [CrossRef]
  33. García-Cortés, S.; Galán, C.O.; Argüelles-Fraga, R.; Díaz, A.M. Automatic detection of discontinuities from 3D point clouds for the stability analysis of jointed rock masses. In Proceedings of the 18th International Conference on Virtual Systems and Multimedia, VSMM 2012 Virtual Systems in the Information Society, Milan, Italy, 2–5 September 2012; pp. 595–598. [Google Scholar]
  34. Ge, Y.; Tang, H.; Xia, D.; Wang, L.; Zhao, B.; Teaway, J.W.; Chen, H.; Zhou, T. Automated measurements of discontinuity geometric properties from a 3D-point cloud based on a modified region growing algorithm. Eng. Geol. 2018, 242, 44–54. [Google Scholar] [CrossRef]
  35. Lai, P.; Samson, C.; Bose, P. Surface roughness of rock faces through the curvature of triangulated meshes. Comput. Geosci. 2014, 70, 229–237. [Google Scholar] [CrossRef]
  36. Priest, S.D.; Hudson, J.A. Estimation of discontinuity spacing and trace length using scan-line surveys. Int. J. Rock Mech. Min. Sci. Geomech. Abstr. 1981, 18, 183–197. [Google Scholar] [CrossRef]
  37. Commission on The Standardization of Laboratory and Field Test (ISRM). Suggested methods for the quantitative description of discontinuities in rock masses. Int. J. Rock Mech. Min. Sci. Geomech. Abs. 1978, 15, 319–368. [Google Scholar]
  38. Slob, S.; van Knapen, B.; Hack, R.; Turner, K.; Kemeny, J. Method for automated discontinuity analysis of rock slopes with three-dimensional laser scanning. Transp. Res. Rec. 2005, 1913, 187–194. [Google Scholar] [CrossRef]
  39. Mah, J.; Samson, C.; McKinnon, S.D. 3D laser imaging for joint orientation analysis. Int. J. Rock Mech. Min. Sci. 2011, 48, 932–941. [Google Scholar] [CrossRef]
  40. Reid, T.R.; Harrison, J.P. A semi-automated methodology for discontinuity trace detection in digital images of rock mass exposures. Int. J. Rock Mech. Min. Sci. 2000, 37, 1073–1089. [Google Scholar] [CrossRef]
  41. Thiel, K.H.; Wehr, A. Performance capabilities of laser-scanners—An overview and measurement principle analysis. Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2004, 36, 14–18. [Google Scholar]
  42. Vöge, M.; Lato, M.J.; Diederichs, M.S. Automated rockmass discontinuity mapping from 3-dimensional surface data. Eng. Geol. 2013, 164, 155–162. [Google Scholar] [CrossRef]
  43. Mavrouli, O.; Corominas, J.; Jaboyedoff, M. Size distribution for potentially unstable rock masses and in situ rock blocks using LIDAR-generated digital elevation models. Rock Mech. Rock Eng. 2015, 48, 1589–1604. [Google Scholar] [CrossRef] [Green Version]
  44. Battulwar, R.; Emami, E.; Zare Naghadehi, M.; Sattarvand, J. Automatic extraction of joint orientations in rock mass using PointNet and DBSCAN. In Proceedings of the International Symposium on Visual Computing, San Diego, CA, USA, 5–7 October 2020; Springer: Cham, Switzerland, 2020; pp. 1–10. [Google Scholar]
  45. Casagli, N.; Pini, G. Analisi cinematica della stabilità in versanti naturali e fronti di scavo in roccia. Geol. Appl. Idrogeol. 1993, 28, 223–232. [Google Scholar]
  46. Sturzenegger, M.; Stead, D. Quantifying discontinuity orientation and persistence on high mountain rock slopes and large landslides using terrestrial remote sensing techniques. Nat. Hazards Earth Syst. Sci. 2009, 9, 267–287. [Google Scholar] [CrossRef]
  47. Goodman, R.E.; Bray, J.W. Toppling of rock slopes. Proc. Spec. Conf. Rock Eng. Found. Slopes. 1976, 2, 201–234. [Google Scholar]
  48. Farina, P.; Bardi, F.; Lombardi, L.; Gigli, G. Combining structural data with monitoring data in open pit mines to interpret the failure mechanism and calibrate radar alarms. In Proceedings of the 2020 International Symposium on Slope Stability in Open Pit Mining and Civil Engineering; Dight, P., Ed.; Australian Centre for Geomechanics: Perth, Australia, 2020; pp. 523–534. [Google Scholar] [CrossRef]
  49. Kazhdan, M.; Bolitho, M.; Hoppe, H. Poisson surface reconstruction. In Proceedings of the Fourth Eurographics Symposium on Geometry Processing, Goslar, Germany, 26–28 June 2006; pp. 61–70. [Google Scholar]
  50. Castagnetti, C.; Bertacchini, E.; Capra, A.; Dubbini, A. Terrestrial laser scanning for preserving cultural heritage: Analysis of geometric anomalies for ancient structures. In Proceedings of the FIG Working Week 2012—Knowing to Manage the Territory, Protect the Environment, Evaluate the Cultural Heritage, Rome, Italy, 6–10 May 2012. [Google Scholar]
  51. Chías, P.; Abad, T.; Echeverría, E.; De Miguel, M.; Llorente, P. Underground surveying: 16th-century cellar vaults in the Gallería de Convalecienties, monastery of San Lorenzo del Ecorial. In Proceedings of the International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences, Nafplio, Greece, 1–3 March 2017; Volume XLII2/W3, pp. 179–186. [Google Scholar] [CrossRef] [Green Version]
  52. Grussenmeyer, P.; Alby, E.; Landes, T.; Koehl, M.; Guillemin, S.; Hullo, J.F.; Assali, P.; Smigiel, E. Recording approach of heritage sites based on merging point clouds from high resolution photogrammetry and terrestrial laser scanning. Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2012, XXXIXB5, 553–558. [Google Scholar] [CrossRef] [Green Version]
  53. Olsen, M.; Kuester, M.; Chang, B.; Hutchinson, T. Terrestrial laser scanning-base structural damage assessment. J. Comput. Civ. Eng. 2010, 24, 264–272. [Google Scholar] [CrossRef]
  54. Assali, P.; Grussenmeyer, P.; Villemin, T.; Pollet, N.; Viguier, F. Surveying and modeling of rock discontinuities by terrestrial laser scanning and photogrammetry: Semi-automatic approaches for linear outcrop inspection. J. Struct. Geol. 2014, 66, 102–114. [Google Scholar] [CrossRef]
  55. Gigli, G.; Mugnai, F.; Leoni, L.; Casagli, N. Brief communication “Analysis of deformations in historic urban areas using terrestrial laser scanning”. Nat. Hazards Earth Syst. Sci. 2014, 9, 1759–1761. [Google Scholar] [CrossRef]
  56. Xin, W.; Lejun, Z.; Xiaohua, S.; Yupeng, R.; Yi, Q. A region-growing approach for automatic outcrop fracture extraction from a three-dimensional point cloud. Comput. Geosci. 2017, 99, 100–106. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Map of the Cima Grappa Military Shrine. Ex-infirmary and underground tunnels network.
Figure 1. Map of the Cima Grappa Military Shrine. Ex-infirmary and underground tunnels network.
Ijgi 10 00276 g001
Figure 2. View of (a) vault collapse within the tunnel; (b) limestone outcrop; (c) vault protected by metal liner; (d) limestone outcrop + concrete coating.
Figure 2. View of (a) vault collapse within the tunnel; (b) limestone outcrop; (c) vault protected by metal liner; (d) limestone outcrop + concrete coating.
Ijgi 10 00276 g002
Figure 3. Geo-mechanical survey with the scan-line method (S4).
Figure 3. Geo-mechanical survey with the scan-line method (S4).
Ijgi 10 00276 g003
Figure 4. Laser scanner working in correspondence of a vault collapse.
Figure 4. Laser scanner working in correspondence of a vault collapse.
Ijgi 10 00276 g004
Figure 5. Top view of the infirmary from 3D point cloud. S from 1 to 4 are the scanning positions.
Figure 5. Top view of the infirmary from 3D point cloud. S from 1 to 4 are the scanning positions.
Ijgi 10 00276 g005
Figure 6. Top view of the infirmary’s point cloud (A) and cross-sections (BD). Dimensions are expressed in meters.
Figure 6. Top view of the infirmary’s point cloud (A) and cross-sections (BD). Dimensions are expressed in meters.
Ijgi 10 00276 g006
Figure 7. Identified plains on an outcrop’s portion point cloud.
Figure 7. Identified plains on an outcrop’s portion point cloud.
Ijgi 10 00276 g007
Figure 8. Point cloud acquired from scan position no. 4 used to search for discontinuity plans.
Figure 8. Point cloud acquired from scan position no. 4 used to search for discontinuity plans.
Ijgi 10 00276 g008
Figure 9. Extraction of discontinuity planes from point cloud TLS. (a) Step1: identify the flat sectors (blue areas) and germination of discontinuities on an organized point cloud. (b) Step2: progressive germination by adding co-planar points with the initial nucleus. (c,d) Definition of the 3D polygons that delimit the discontinuity planes identified.
Figure 9. Extraction of discontinuity planes from point cloud TLS. (a) Step1: identify the flat sectors (blue areas) and germination of discontinuities on an organized point cloud. (b) Step2: progressive germination by adding co-planar points with the initial nucleus. (c,d) Definition of the 3D polygons that delimit the discontinuity planes identified.
Ijgi 10 00276 g009
Figure 10. Extraction of discontinuity planes: step 3. (a) Union of the extracted polygons for the different scanning positions and eliminate redundancies. (b) Union of the extracted polygons for the different scanning positions and elimination of redundancies. The polygons are projected onto the coloured 3D model according to the local immersion direction.
Figure 10. Extraction of discontinuity planes: step 3. (a) Union of the extracted polygons for the different scanning positions and eliminate redundancies. (b) Union of the extracted polygons for the different scanning positions and elimination of redundancies. The polygons are projected onto the coloured 3D model according to the local immersion direction.
Ijgi 10 00276 g010
Figure 11. 3D representation of the discontinuity plans remotely extracted and divided by family. White: no family; red: BG; yellow: JN1; green: JN2; blue: JN3.
Figure 11. 3D representation of the discontinuity plans remotely extracted and divided by family. White: no family; red: BG; yellow: JN1; green: JN2; blue: JN3.
Ijgi 10 00276 g011
Figure 12. Stereographic projections of the poles and discontinuity planes’ concentration.
Figure 12. Stereographic projections of the poles and discontinuity planes’ concentration.
Ijgi 10 00276 g012
Figure 13. Main mechanisms of rock instability and kinematic conditions for structurally controlled failures occurrence (αd = discontinuity dip direction, αs = slope dip direction, βd = discontinuity dip angle, βs = slope dip angle, φ = friction angle) [48].
Figure 13. Main mechanisms of rock instability and kinematic conditions for structurally controlled failures occurrence (αd = discontinuity dip direction, αs = slope dip direction, βd = discontinuity dip angle, βs = slope dip angle, φ = friction angle) [48].
Ijgi 10 00276 g013
Figure 14. Detail mesh (in purple) of the underground portions where the rock mass is outcropping.
Figure 14. Detail mesh (in purple) of the underground portions where the rock mass is outcropping.
Ijgi 10 00276 g014
Figure 15. Contours with equidistance 10 cm relative to the underground portions where the rock- mass is outcropping.
Figure 15. Contours with equidistance 10 cm relative to the underground portions where the rock- mass is outcropping.
Ijgi 10 00276 g015
Figure 16. High-resolution cloud of points relating to underground infirmary rooms. The white circles indicate the position of the larger collapses. The green and blue planes indicate the position and orientation of two bands of weakness, parallel respectively to the families of discontinuities JN2 and JN1 identified.
Figure 16. High-resolution cloud of points relating to underground infirmary rooms. The white circles indicate the position of the larger collapses. The green and blue planes indicate the position and orientation of two bands of weakness, parallel respectively to the families of discontinuities JN2 and JN1 identified.
Ijgi 10 00276 g016
Figure 17. Images of collapsing zones: (a) picture of collapsing niche 3, (b) point clouds of collapsing niche 3.
Figure 17. Images of collapsing zones: (a) picture of collapsing niche 3, (b) point clouds of collapsing niche 3.
Ijgi 10 00276 g017
Figure 18. Kinematic analysis: (a) planar failure mechanism; (b) wedge failure; (c) toppling failure (d) free fall.
Figure 18. Kinematic analysis: (a) planar failure mechanism; (b) wedge failure; (c) toppling failure (d) free fall.
Ijgi 10 00276 g018
Table 1. Laser scanner hardware parameters.
Table 1. Laser scanner hardware parameters.
ParametersValue
Measurement Range0.2–1000 m
Accuracy10 mm
Laser class1
Minimum angle step width0.004°
Table 2. Values (dip and dip direction) of joint families recorded by operators using classic scan-line survey and values of joint families extracted by the remote geo-mechanical survey.
Table 2. Values (dip and dip direction) of joint families recorded by operators using classic scan-line survey and values of joint families extracted by the remote geo-mechanical survey.
Scan-Line SurveyRemote
JN1 059/78°JN1: 65°/76°
JN2 002/87°JN2: 001°/82°
JN3 115/75°JN3: 121°/75°
BG 250/09°BG: 245°/13°
Table 3. Hazard kinematic indexes [49], N = number of poles satisfying a specific failure condition; I = number of intersections satisfying a specific failure condition.
Table 3. Hazard kinematic indexes [49], N = number of poles satisfying a specific failure condition; I = number of intersections satisfying a specific failure condition.
Instability MechanismHazard Index
Planar failure C p f = N p f N   %
Wedge failure C w f = I w f I   %
Flexural toppling C f t = N f t N   %
Block toppling C b t = N b t N · I b t I   %
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Mugnai, F.; Farina, P.; Tucci, G. Exploiting a Semi-Automatic Point Cloud Segmentation Method to Improve the Quality of Rock-Mass Characterization. The Cima Grappa Conservative Restoration Case Study. ISPRS Int. J. Geo-Inf. 2021, 10, 276. https://doi.org/10.3390/ijgi10050276

AMA Style

Mugnai F, Farina P, Tucci G. Exploiting a Semi-Automatic Point Cloud Segmentation Method to Improve the Quality of Rock-Mass Characterization. The Cima Grappa Conservative Restoration Case Study. ISPRS International Journal of Geo-Information. 2021; 10(5):276. https://doi.org/10.3390/ijgi10050276

Chicago/Turabian Style

Mugnai, Francesco, Paolo Farina, and Grazia Tucci. 2021. "Exploiting a Semi-Automatic Point Cloud Segmentation Method to Improve the Quality of Rock-Mass Characterization. The Cima Grappa Conservative Restoration Case Study" ISPRS International Journal of Geo-Information 10, no. 5: 276. https://doi.org/10.3390/ijgi10050276

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop