Next Article in Journal
Game Theoretic Honeypot Deployment in Smart Grid
Next Article in Special Issue
Multi-Constellation Software-Defined Receiver for Doppler Positioning with LEO Satellites
Previous Article in Journal
High Precision Wide Bandwidth DC Current Transducer Based on the Platiše Flux Sensor
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Pavement Crack Detection from Mobile Laser Scanning Point Clouds Using a Time Grid

1
College of Geological Engineering and Geomatics, Chang’an University, Xi’an 710054, China
2
State Key Laboratory of Resources and Environment Information System, Institute of Geographic Sciences and Natural Resources Research, CAS, Beijing 100101, China
*
Author to whom correspondence should be addressed.
Sensors 2020, 20(15), 4198; https://doi.org/10.3390/s20154198
Submission received: 1 July 2020 / Revised: 15 July 2020 / Accepted: 26 July 2020 / Published: 28 July 2020

Abstract

:
This paper presents a novel algorithm for detecting pavement cracks from mobile laser scanning (MLS) data. The algorithm losslessly transforms MLS data into a regular grid structure to adopt the proven image-based methods of crack extraction. To address the problem of lacking topology, this study assigns a two-dimensional index for each laser point depending on its scanning angle or acquisition time. Next, crack candidates are identified by integrating the differential intensity and height changes from their neighbors. Then, morphology filtering, a thinning algorithm, and the Freeman codes serve for the extraction of the edge and skeleton of the crack curves. Further than the other studies, this work quantitatively evaluates crack shape parameters: crack direction, width, length, and area, from the extracted crack points. The F1 scores of the quantity of the transverse, longitudinal, and oblique cracks correctly extracted from the test data reached 96.55%, 87.09%, and 81.48%, respectively. In addition, the average accuracy of the crack width and length exceeded 0.812 and 0.897. Experimental results demonstrate that the proposed approach is robust for detecting pavement cracks in a complex road surface status. The proposed method is also promising in serving the extraction of other on-road objects.

1. Introduction

Pavement cracking serves as a vital quantitative reference index for road maintenance [1,2]. Currently, large-scale efforts of manual inspection of cracks are rare owing to their low efficiency, high cost, and exposure to vehicle traffic [3,4]. Automatic detection and quantification of pavement cracks from road surface images and videos acquired using a variety of special sensors is becoming more common [5,6,7]. Various methods of crack recognition have been developed based on image processing and computer vision [8]. Excellent performances have been achieved in recent years since the successful application of neural networks and deep learning technology for image segmentation [9,10,11,12,13,14,15,16,17,18,19]. However, road surface images and videos are often obscured by light, shadows, stains, rust, and noise, making crack recognition challenging. Moreover, two-dimensional (2D) detection technology is limited regarding the illustration of the terrain and details of road damage, which reduces the performance of crack evaluation. Researchers have been exploring the methods for three-dimensional (3D) pavement crack detection and quantification. For example, the application of a line–structured light scanner and 3D camera can collect the detailed 3D terrain textures (with an accuracy of 1 mm or better) of pavement profiles [20,21,22,23,24,25,26]. However, the advanced 3D measurement systems that integrate numerous technologies are expensive, which hinders their large-scale application in pavement inspection [27].
In recent years, mobile laser scanning (MLS) systems have been gradually applied to fast, high-precision road and pavement condition surveys, and some researchers have investigated the applicability of MLS point clouds to crack extraction. Two kinds of algorithms were summarized in Reference [28]: (1) Projecting dense 3D point clouds into 2D images and extracting cracks using image-based algorithms [29], and (2) identifying cracks directly from a discrete point cloud. For example, Guan et al. [30,31] interpolated the point cloud into geo-referenced feature (GRF) images with ground resolution of 2 cm by applying a modified inverse distance-weighted (IDW) algorithm. Then, they picked out five test areas that contained a variety of cracks, ranging in width from a few centimeters up to 10 cm. The low-intensity pixels were selected as possible cracks using the maximum entropy threshold algorithm and enhanced by tensor voting to highlight the crack curves. The authors compared their results with that from pavement images. Their method was verified as effective for pavement cracks with low contrast and bad continuity and demonstrated much better performance compared to road surface images and videos. However, its computational time is considerable and mostly consumed by the iterative tensor voting process. Moreover, points with low intensities may be due to wet pavements, low-reflectance repair materials, stains, or other noises. In their study, Chen and Li rasterized all the ground points into a 4 cm digital terrain model (DTM) by selecting the minimum height within the grid as the pixel value [32]. Next, they used a high-pass filter to convolute the DTM to detect local elevation change. Finally, pavement cracks were delineated using a two-step matched filter. This method was computationally efficient and simple when compared to Guan’s. Convoluting a 2500 × 760 DTM only took a few seconds. Another advantage of the proposed method was that it successfully detected some cracks of only a few millimeters width. However, slight cracks are difficult to distinguish from the normal pavement height. In addition, it is difficult to select the grid size to retain the features, such as the crack width, needed for crack quantification. Without dimension reduction, Yu et al. [33] extracted 3D crack skeletons directly from mobile LiDAR point clouds. They normalized the intensities of the point data and identified crack candidate points with lower intensities using Otsu’s thresholding method. Next, they grouped crack-lines using a Euclidean distance clustering method. Finally, they extracted 3D crack skeletons by applying a L1-medial skeleton extraction method. This method executed quickly because it avoided time-consuming neighborhood searches in dense and discrete point clouds. However, crack points only account for a very small proportion of MLS data, causing Otsu’s threshold to be unreliable regarding segment cracks.
From all the research above, cracks often exhibit low-intensity value and can be detected by the intensity-threshold segmentation method. What is more, the method of intensity-threshold segmentation is independent of neighborhood contexts. Therefore, it can be applied to both images and point clouds. However, distinguishing cracks using an intensity threshold is not always reliable when encountering certain interference factors, such as uneven intensity, road markings, stains, low-reflectance mending materials, and sensor noise. The methods that interpolate dense 3D point clouds into 2D images can provide the neighbors information to delineate cracks by elevation changes. However, some cracks are too shallow to be detected by MLS data with a coordinate resolution of the mm-level. The completeness and robustness of extracted cracks only by low-intensity values or elevation changes may be insufficient. Also, high-resolution rasterization of massive point clouds consumes time and resources. More importantly, using a fixed size grid to interpolate point cloud weakens the quantitative evaluation performance of crack width.
This study attempts to propose a novel algorithm to detect 3D pavement cracks directly from MLS data and to quantitatively evaluate their shape parameters. To adapt to a noisy background, the method selects crack candidates based on the integration of the intensity difference and height changes from their neighbors. To address the problem of lacking topology in MLS data, this paper introduced a time grid (Tgrid) for the data collected by the 2D scanner with the same incremental angles. Each point is assigned a row and column number based on either its scan angle or acquisition time. Thus, through a point’s row and column numbers, it is easy to find its neighbor points. Moreover, with the regular grid indexes, MLS data can be processed like images, and then crack detection can be performed using proven image-based methods. From this, the Tgrid can also serve the detection of other on-road objects from MLS data, such as curbs, road markings, and manholes. Moreover, it offers promising potential for application in convolutional neural networks and deep learning with MLS point clouds. The proposed method was tested in road conditions over a 220 m stretch of road with hundreds of cracks. The remainder of this study is organized as follows: Section 2 introduces the materials and details the proposed method; Section 3 presents the results and discussion; and Section 4 presents our conclusions.

2. Materials and Methods

2.1. Materials

The test dataset shown in Figure 1a was collected using the SSW-IV vehicle laser modeling and measurement system, including an approximately 220 m long highway with numerous cracks. A 2D scanner with a 360° field of view (FOV) was assembled at the rear of an SUV roof 2.2 m above the road surface and oriented perpendicular to the moving direction. The 2D scanning plane was set at an angle of 45 degrees inclined from the horizontal plane. The scanner provided 500,000 measurements and 200 scan lines per second, featuring a detection range of up to 300 m. The concrete road surface was measured at a speed of 30 km/h, resulting in the spacing between adjacent scan lines of approximately 4 cm and distance between adjacent laser points of 0.7–3 cm in the road. Figure 1b shows the common phenomenon of divergent point intensity in two lanes. The point intensities of the lane on the other side were so low (black box shown in Figure 1b) that some cracks were difficult to distinguish with a fixed intensity threshold. In addition, the road surface contained road markings, ruts, stains, black spots, and other intensity-altering factors, which would have an unfavorable effect on algorithms that are based on a global intensity threshold. This situation indicates that the use of point intensity differences and altitude changes to the background may be more effective than a fixed intensity threshold alone.

2.2. Methods

As a linear feature, cracks should be continuous; however, discrete MLS points have no connectivity. This study attempts to identify crack points from their neighbors, which are defined by an especially topological structure within which the discrete MLS points are grid-adjacent. The introduced structure is built based on the imaging mechanism of the MLS system. Figure 2a sketches the arrangement of the MLS point cloud collected over time using a 2D scanner with equal incremental angles. The laser beams are projected onto different scan lines as the vehicle moves. The neighbors of the crack points on the same scan line can be simply identified according to the collection sequence of the point data, whereas neighbors across multiple scanning lines can be deduced using similar scanning angles. Based on the above assumptions, we created a 2D grid index for MLS point clouds, as shown in Figure 2b. We named the structure Time grid (Tgrid) because it records the arrangement of the MLS point cloud over time.
The proposed method comprises three steps: (i) Construction of the Tgrid for the MLS point clouds, (ii) detection of road surface point and crack candidates, and (iii) generation of the crack skeleton and calculation of the crack shape parameters.

2.3. Construction of a Tgrid for MLS Point Clouds

Let P(X, Y, Z, T, I) represent the test datasets with N points, where (X, Y, Z) denotes the spatial coordinates of a point, and T and I express its acquisition time and intensity. The Tgrid records the 2D arrangement of MLS points along and across the scan lines. Each point is assigned a grid index (i, j), from which it can search for its neighbors. A point with the index (i, j) is located on the ith scanning line, and its scanning angle θ is j times the angular resolution ∆θ. The values of j can be positive or negative.
This study provides two approaches to construct a Tgrid using scanning angle (θ) or acquisition time (T). Figure 3 shows a flow chart of the two approaches.
As shown in Figure 4a, the scan angle, θ, increases or decreases monotonously with T along the same scan line. Thus, we divided the point clouds to their scanning lines by detecting the sign change of the angle difference between the adjacent laser points. When the neighboring point satisfies Sgn(θk+1 − θk) ≠ 0, k = 1, 2, …, N, the kth point is a piecewise point. Let m represent the number of k detected. The numerical array [1, k1, k2, …, km, N] then forms (m + 1) non-overlapping intervals, which correspond to (m + 1) scan lines. On each scanning line, the multiples of scan angles, θ, with respect to the angular resolution, ∆θ, is calculated to obtain the j value of each point, j = θ/∆θ.
When θ is not available or the MLS data provides a scan angle level with insufficient accuracy, the above method is ineffective. In such a situation, we propose a different approach of constructing Tgrid based on T. In our previous study [34], we deduced that the time series of the scanner ground track, i.e., the laser points that the scanner projects directly below, should be in accordance with an isochromatic sequence. As shown in Figure 4b, let TG be the time series of the scanner ground track, TG = {TG1, TG2, TG3, …, TGm+1}, and TGi+1TGi = Δt, where Δt is the elapsed time on each scan line. The time range [TGi − Δt/2, TGi + Δt/2] divides the laser points to the ith scanning line. The j indexes of points on ith scanning line can be derived as j = T T G i / Δ T + 0.5 , where ΔT denotes the emission time interval of the laser beams. If ∆T is unknown, it can be estimated using the time difference (Ti+1Ti) between the adjacent laser points. Note that (Ti+1Ti) may change over time. We use the symbol ⌊  ⌋ to reduce the estimation error (for more detail, please refer to Reference [34]). In the Tgrid, each point (i, j) has a maximum of four side neighbors: (I − 1, j), (i + 1, j), (i, j − 1), and (i, j + 1), and four diagonal neighbors: (i − 1, j − 1), (i + 1, j − 1), (i − 1, j + 1), and (i + 1, j + 1). The spatial position of these neighbors is free. In addition, eight neighbors may not always exist. When there are no laser beam returns and no neighbors, they are shown as empty nodes (as shown in (i − 1, j − 1) in Figure 2b).
For easy understanding, we added two columns (xt, yt) before P(X, Y, Z, T, I) to the expressed MLS point clouds as a Tgrid map, G(xt, yt, X, Y, Z, T, I), where (xt, yt) represents the row and column of a point in the Tgrid map. Similar to a multi-band image, (X, Y, Z, T, I) represent the multi-band characteristics. The cost of building and maintaining such a 2D grid is quite low.

2.4. Detection of Road Surface Points and Crack Candidates

2.4.1. Detection of Road Surface Points

Various algorithms have been developed to extract roads from MLS point clouds. The advantage of using the Tgrid is that the road boundaries can be defined as a continuous curve. Figure 5a shows the cross-section of the studied road with natural slopes on both sides. To identify road surface points quickly, we used the method in Reference [35], which roughly approximated the local road patch as a 3D best-fitting plane. Next, this method took the trajectory data as the initial road seed to search the road boundary area along and perpendicular to the scan line using a statistical hypothesis testing method based on the point altitudes. Figure 6 shows a flow chart of the method of detection of road surface points. First, the entire Tgrid is divided into non-overlapping sub-blocks, the size of which is Wxt × Wyt. The sub-block located by the scanner ground track serves as the initialized searching sub-block, P0 (shown in red in Figure 5b). Next, both the road surface and road boundary areas are searched along xt and yt directions through a statistical hypothesis test based on the point altitudes in the sub-block. Then, the possible road points are extracted from the sub-blocks of the road surface and road boundary and finally, identified by finding the largest connected region in the Tgrid.
Assume the point altitude, ZS, in a road patch conforms to the Gaussian distribution of N (μ, σ) and that μ is unknown. Let the number of points within the sub-block be Nb. A statistical test problem for the population variance with unknown mean is as follows:
H 0 :   σ 2     σ 0 2 ,   H 1 : σ 2   <   σ 0 2
Z ¯ = 1 N b i = 1 N b Z s i
s 2 = 1 N b 1 i = 1 N b ( Z s i Z ¯ ) 2
where Z ¯ denotes the mean value and s2 is the standard deviation. The statistical variable and its region of rejection are:
χ 2 = N b 1 σ 0 2 s 2
χ 2 χ 2 1 α
where α is the significance level or the confidence; generally, α = 0.05. If the statistical variable of the sub-block is located inside of the region of rejection, H1 is accepted (shown in green in Figure 4b) and H0 is denied. This means that no road boundary area is found; otherwise, H0 is accepted (shown in blue in Figure 4b). The associated sub-block is selected for the road boundary area. Figure 4b illustrates the process of detecting this area. The proposed method can be summarized as follows:
(1)
Set an initialized searching sub-block, P0, where trajectory data are located, and calculate the statistical variance, σ02, of the point altitudes within it.
(2)
Search forward to find the road boundary area along the direction of xt, yt, and the opposite xt and opposite yt, until a sub-block (shown in blue) that passes the test is found. Pause searching.
(3)
Iterate for all trajectory points until all road boundary areas are beside each other.
(4)
Filter the non-road points in the road boundary area (shown in blue) and its closed road surface neighbors (the closed sub-block shown in green) using a height threshold, hth, to the local road plane.
(5)
Extract the complete road surface by finding the largest connected region in the Tgrid.
The size of the sub-block (Wxt × Wyt) can be determined through the geometric dimension, WD. Although the distance of the adjacent points is inconsistent in yt directions, we unified the Wyt for convenience. The mean spacing between the adjacent scanner ground track, d s (in xt direction), and the largest distance between the adjacent road grid points, d a (in yt direction), are used to determine Wxt and Wyt. This means that Wxt =   W D / d s + 0.5 and Wyt = W D / d a + 0.5 . Due to the fixed size of the sub-blocks, it is also possible that the road boundary lies in the road sub-block adjacent to the boundary area. Thus, we derive the local road plane by fitting the 3D coordinates of the points in a road boundary area (shown in blue) and its closed road neighbors (the closed sub-block shown in green) using the robust least-squares method. Hence, more road points are participating in the fitting of the local plane. The conformity of the fitting plane to the road surface is thereby ensured. The height of points involved in the fitting operation to the local plane is taken as the evaluation standard of whether this is a non-road point. The height threshold, hth, can be given according to the roughness of the road surface.

2.4.2. Detection of Crack Candidates

The crack points are detected by two approaches: detection (1) of the valley point that is lower than local pavement and (2) of the low-intensity point with height no higher than the road surface.
(1) Detection of the valley point
The valley point has two typical characteristics: a height difference exceeding the road surface roughness and a nearby gradient change. We defined a circular neighbor region with radius, LR, that is larger than the common maximum width of the crack. An M estimator of Andrew’s type was applied to derive the robust altitude, ZM, of the road surface and the standard interquartile interval, ZNIQR, is used to measure the degree of the pavement roughness. The two conditions are combined as follows:
G i , j : i f ( Z i , j > Z M h d ) p e n e t r a t i o n   p o i n t i f ( Z i , j > Z M h d   &   Z i , j Z M Z NIQR   &   grad f i , j v t h ) v a l l e y   p o i n t o t h e r w i s e , i r r e l e v a n t   p o i n t
where hd expresses the maximum possible depth of the crack, which is designed to eliminate the penetration laser point caused by loose pavement facilities, such as rainwater grates, well covers, etc. Next, gradf(i, j) represents the average gradient value centered on (i, j), which is calculated through the geometric coordinates (X, Y, Z) of eight neighbors using Equation (7). Note that only the existing neighbors are calculated. Following this, vth defines the minimum gradient change requirement of the crack candidates. Figure 7 depicts the characteristics of the valley points. The extracted valley points were recorded to dataset GV:
grad f ( i ,   j ) = 1 k 1 k gradient i , j , k gradient i , j , k   = Z i , j , k Z i , j ( X i , j , k X i , j ) 2 + ( Y i , j , k Y i , j ) 2
where k is the number of the existing neighbors centered on (i, j).
(2) Detection of low-intensity points with height no higher than the road surface
Some slight cracks have only a few points with the valley effect, and further candidates are points with low intensities. However, the intensity of the road surface points varies greatly owing to many factors. The low-intensity point clouds do not necessarily belong to cracks, and vice versa. In addition, the proportion of the crack points is so small that it is difficult to segment through a threshold automatically provided by the commonly used method. Our approach is to first detect the low-intensity points and then filter out those that are higher than the road surface, thereby retaining the points with low-intensity that are no higher than the road surface.
This study detects the lower intensity points through a background-differential map, within which the crack proportion limitation (TP) is preset to reduce the unfavorable effects of any generous non-crack points. Figure 8 shows the flow of the detection method.
We first generated the background map, GB, using a median filter from the input intensity Tgrid map, GI, which derives from the grid index (i, j) and its intensity. Next, the background-differential map, GF, is obtained using a minus operation as follows: GF = GI − GB. In Figure 8b, the proportion limitation of the crack points (CP) is given and the relevant truncation value, FB, is calculated on the histogram of GF. The differential truncation map, GFC, is generated in Figure 8c by transforming the value of GF(i,j) that is higher than FB to FB. The maximum entropy method is then adopted to determine the segment threshold, TC, of the possible crack points (recorded in GP) and non-crack points. Finally, GP is filtered through the valley effect to remove the nearest non-crack point on both sides of the GP point along the tx and ty directions to compare their altitude. As shown in Figure 8f, two rectangles frame the crack points in the directions of tx and ty respectively, with non-crack points intercepted at two ends when centered on a candidate. If the altitude of the candidate is no greater than the largest of the four non-crack neighbors, the candidate is reserved; otherwise, it is filtered out. The remaining crack candidates are saved to the binary map, GC.
Some crack candidates obtained by these two approaches may overlap. Thus, the union of GV ∪ GC is taken as the candidate crack points.

2.5. Generation of the Crack Skeleton and Calculation of Crack-Shape Parameters

2.5.1. Generation of the Crack Skeleton

The extracted crack candidates may be noisy and discontinuous. Thus, a series of morphological operations are implemented to extract the continuous crack skeleton. First, to eliminate small holes and connect the nearest candidates, the morphology closing algorithm is adopted to cluster the crack regions. Next, the one-point-thick continuous crack curves are extracted using the thinning algorithm proposed in Reference [36]. Following this, the burrs and branches are separated according to the connectivity of the points on the skeleton.
Morphological image processing techniques are often used to connect crack fragments into crack curves and for complicated crack shapes in image-based methods [2,37]. Unlike the geometric grid, the Euclidean distance of the adjacent points on the Tgrid is constantly changing. The size of structural elements of the closing algorithm is not fixed but is defined by a geometric dimension. As shown in Figure 9, we set a double-layer search structure, the first layer of which is a rectangle with a fixed size, and the second layer of which is a preset geometric search circle with the radius D2. Given that the width and the height of the rectangle are RW and RH, they can be determined by the distance between the densest point cloud under the vehicle in the direction of xt and yt. We assumed these to be dx and dy, respectively. Thus, RW = 2 × D 2 / d x + 0.5 + 1, RH = 2 × D 2 / d y + 0.5 + 1. The points within the second layer form the structural elements of the morphology algorithm. Note that some empty nodes without geometric coordinates may be expanded during this process. To mark the connected regions of the cracks using the 8-connected neighborhood, the 3D coordinates of the expanded empty nodes are linearly interpolated using the closest existing neighbors at both sides along the ty direction.
To remove the isolated noise, we extract the maximum width (represented by the number of points in the tx direction) and height (the number of points in the ty direction) of the connected regions and recover them to the geometrical space to calculate their geometrical dimensions. If the longest of two sides is larger than the set length threshold, Lth1, the connected region indicates a curve region. Otherwise, it is marked as noise to be removed. The remaining regions are thinned to extract the continuous crack skeletons.
To remove the burrs, the connectivity of points in the skeletons are distinguished as endpoints, general skeleton points, and intersections, which correspond to points that are single-connected, dual-connected, and connected by no less than three. The burrs are recognized by the space distance threshold, Lth1, between the endpoint and the intersections that were first encountered. Each skeleton only leaves one longest main curve, and the rest are defined as branches, i.e., another skeleton. The skeleton curve is traced by the Freeman chain code and the coordinates of crack points are recorded to a table to calculate the crack-shape parameters.

2.5.2. Calculation of Crack-Shape Parameters

The crack-shape parameters here include crack direction, width, length, and area.
(1) The crack direction
To determine the crack direction, the crack points and skeleton are recovered to the geometric space. We use the length/width ratio RLW of the smallest circumscribed rectangle of the skeleton points and its long axis direction to determine the crack direction. When RLW exceeds the set threshold, RLWth, the direction of the crack must be evaluated in sections. Otherwise, the crack is oriented by the angle α between the long axis direction of the minimum circumscribed rectangle and the vehicle path. The crack direction is classified as transverse, vertical, or oblique with the related inclination ranges of the α value [0, π/6], (π/6, π/3), and [π/3, π/2]. In Figure 10, crack 2# and crack 3# are classified as transverse and vertical due to the qualified RLW and the scope of the values of α. Crack 1# has an RLW exceeding the specified threshold and is further separated into two sections, the piecewise points of which are given at the farthest point from the straight-line connecting head and tail. The RLW ratio of the smallest external rectangle of the cracks is recalculated and split again until RLWRLWth, and the angle α between the long axis direction and the path line is taken to identify the trends of the segmented cracks.
(2) The crack width
The minimum distance method is often used to calculate the crack width in image-based approaches; that is, to take the shortest distance from a point on one edge to the other as the crack width at this edge point. Equation (8) expresses the calculation, where wi corresponds to the width of the ith point of the left edge. Next, j = 1, 2, 3, …, NR, where NR denotes the number of points on the right edge. However, limited by the geometric resolution of the MLS data, the extracted cracks contain many single-point-thickness points, particularly the transverse cracks along the scan line:
w = min   ( X i X j 2 + Y i Y j 2 )
This study selects the single-point-thickness edge point, from which a search radius, RS, is given to fit the edge points within the search range to a virtual edge line. The maximum distances wiL and wiR from the edge points to the line on both sides are calculated, and wi = max (wiL, wiR) is taken as the related width of the point. The procedures of the crack width calculations are as follows:
  • In the Tgrid, the Freeman chain code is used to track the edge of the crack connection area. The closed edge curve is disconnected from the endpoint of the skeleton to split the edge into the left and right borders (Figure 11a).
  • The edge points shared by the left and right borders are searched and marked as single-point-thickness edge points. In the above method, wi = max (wiL, wiR) is adopted to calculate the related crack width (Figure 11b), while other edge points use Equation (8).
  • Output the average width of all edge points to measure the severity of the cracks. The maximum width of the crack and its corresponding location, Pm(i, j, X, Y, Z, T, I), serve as the supplementary information.
(3) The crack length
The crack length is calculated by accumulating the distance of the adjacent points from the crack skeleton using Equation (9). Here, the complex network of cracks is evaluated separately from the main skeleton curves and branches:
L D = k = 1 N L 1 X k + 1 X k 2 + Y k + 1 Y k 2
where (Xk, Yk) and (Xk+1, Yk+1) express the plane coordinates of two adjacent crack skeleton points, assuming there are NL skeleton points in the crack.
(4) The crack area
The crack area is obtained by adding up the total area of the triangulated irregular network (TIN) produced by the crack edge points. For the single-point-thickness region, the TIN is flattened. The single-point-thickness region can be roughly estimated using the length of the skeleton multiplied by the average width of the corresponding point.

3. Results and Discussion

The pavement cracks in the test data were automatically extracted on a computer with 16 GB RAM and 4 CPU cores (each @3.6 GHz). Table 1 lists the optional parameters used in the process. The parameters WD, α, and hth in Section 2.4.1 represent the geometric dimension of the initialized searching sub-block P0, the significance level or the confidence of the statistical test problem, and the roughness of the road surface, respectively. The parameters hd and vth express the maximum possible depth and the minimum gradient change requirement of the valley point. The parameter CP in Section 2.4.2 presets the proportion of the crack points in the number of all point clouds. The parameter D2 is the preset searching radius of the second layer of the structural elements in the morphological closing algorithm. The parameter Lth1 denotes the minimum length threshold of the connected region of crack points. The parameter Rlwth in Section 2.5.2 is the threshold of the length/width ratio to determine the crack direction. The parameter Rs illustrates a search radius to define a virtual edge line when dealing with the single-point-thickness edge point.
Figure 12 shows the intensity Tgrid map of the road surface, of which each point with the regular numbers of rows and columns illustrates the point adjacency. As shown, there are many transverse and longitudinal cracks in the pavement, and the width of the longitudinal cracks is expanded near the path but compressed far from the path due to the change of geometric resolution along the scanning line. Moreover, the transverse cracks in the Tgrid map are extremely narrow and often appear to be of single-pixel-thickness.
Figure 13a shows the extracted crack candidates, which contain significant noise. Figure 13b expresses the remaining crack points after morphological filtering and Figure 13c overlies the crack curves to the road surface. To enhance the display effect, the crack points are magnified in Figure 11.
For a quantitative assessment of the results, we obtained the ground truth by means of manual measurement and field investigation. Due to the geometric resolution of the MLS data, we compared only the extracted cracks where the width was above the minimum calculated value (approximately 1 cm). The number of crack curves detected, those not detected, and the error detection were counted to calculate the precision (P) and recall (R). Precision (P) represents the proportion of all “correctly detected items” to all “actually detected”. Recall (R) refers to the proportion of all “correctly detected items” to all “should be detected”. To evaluate the quality of the proposed method, we used a common evaluation standard, the F1-Measure (also known as F1 score), that combines the results of precision (P) and recall (R). The quotas of P, R, and F1 are calculated using Equations (10) and (11), respectively. Table 2 shows the quantitative results of the proposed method on these three types.
P = T P ( T P + F P ) R = T P ( T P + F N ) ,  
F 1 = 2 P × R P + R ,
where TP is the true positive, i.e., the correctly detected crack curves, while FP expresses the false positive, i.e., the number of cracks detected incorrectly. FN denotes false negative, i.e., the number of crack curves not detected.
As shown in Table 2, the performance of detection for the three types of cracks varies greatly. Taking the F1 score as the evaluation standard, the scores of transverse crack, longitudinal crack, and oblique crack are 96.55%, 87.09%, and 81.48%, respectively. Among them, the transverse crack has the highest detection rate (R), while the R is relatively poor for the longitudinal and oblique cracks (less than 80%). This means that the cracks distributed along the scanning line were detected successfully (R = 98.00%), and the cracks across multiple scanning lines went undetected. In addition, the false detection rate of transverse and longitudinal cracks was very low. This score was followed by that of the oblique cracks. This result may be related to the aggregation effect of crack point distribution on the Tgrid map, in which the points of transverse and longitudinal cracks were basically clustered vertically and horizontally while the points of oblique crack were distributed in a relatively scattered fashion. These distributions resulted in poor performance with the lowest detection rate and the highest error detection rate. Finally, three potholes were mistakenly identified as cracks.
To quantify the accuracy of the crack extraction, 60 curves were selected to test the width and length conformance of three types of cracks with 20 curves for each type. Equation (12) was used to calculate the width and length conformance, represented by CW and CL:
C w = 1 T w E w / T w   C L = 1 T L E L / T L    
where TW and EW represent the average width of the actual crack and extracted crack respectively, and TL and EL are the total length of the actual crack and extracted crack. Figure 14 shows the histogram of the width and length compliance of 60 extracted cracks at 0.05 intervals, in piles of three types of cracks.
The width compliance of the transverse cracks was relatively scattered and that of the longitudinal cracks was the most concentrated. For the length compliance, the oblique cracks were dispersed and the other two were relatively centralized. Table 3 shows the quantitative results of the 60 samples after averaging the compliance. The variance of the width /length compliance illustrates the above trends. Specifically, an incomplete detection of a longitudinal crack, as an abnormal case, increased the variance of this type.
From the mean value of CW and CL, the estimated width of the longitudinal crack was the closest to the actual width, followed by that of the oblique crack. In addition, the width compliance of the transverse curves was the lowest at only 0.812. This may be attributable to the fact that most of the transverse cracks are single-point-width cracks. With only one laser point on the width direction of the crack curve, it is difficult to obtain high-precision calculation results. The width accuracy of the oblique cracks was second, which may owe to its scattered distribution across multiple scan lines, bringing challenges to its identification and width measurement. The longitudinal crack points gather in the tx direction with high geometric resolution. With multiple points in the crack width direction, the longitudinal crack achieved the highest accuracy of the width measurement, i.e., up to 91.05%. It is worth noting that spot checks were also made on the widest section and related position of the extracted cracks. However, no satisfactory results were obtained. Considering the positioning error rate of the MLS system, a position error within ±5 cm is regarded as an accurate detection. Nonetheless, the correct detection rate of the widest section was less than 60%. It may be inferred that the geometric resolution requiring improvement is still the bottleneck of the crack detection from the MLS data.
The length error mainly manifested in the lack of extracted length. Some of the transverse, longitudinal, and inclined cracks suffered a portion of the length loss at the endpoint. There were some cases of one longitudinal or oblique crack splitting into two pieces with one section undetected. We regarded such cases as correct detections with insufficient length. Finally, no large gaps were found in the length compliance of the three types: all were approximately 90%.
To further illustrate the merits of integrating the low-intensity and altitude changes for the detection of crack candidates, we conducted comparative tests. For comparison, we used Guan’s method that adopts the maximum entropy sum thresholding to detect possible crack pixels, and with Yu’s method that used Otsu thresholding. Unlike our proposed method, these other two methods only use intensity values. To have a unified basis for comparison, only the methods used to determine the threshold value could be compared. This excludes the use of the subsampling process to generate interpolation images. The results are shown in Figure 15.
It is evident that using Otsu’s threshold or the maximum entropy threshold extracted over-sized crack candidates. Although they are excellent methods for automatically determining the global threshold, they cannot adaptively adjust to the extraction of cracks with very small proportions in large-scale data ranges. Moreover, the proposed background difference algorithm improves the negative effects of inconsistent intensity due to road roughness, road markings, different pavement materials, etc. In addition, the ruts expressed as raised ridges are not distinguished by low intensities (the black box in Figure 15a,c,d). This study embeds the analysis of the altitude changes near the crack points: the ridged ruts have been removed and the robustness of the method has improved. When compared with the aforementioned threshold-segmentation methods, the proposed method is advantageous in its adaption to more complex environments and global applicability; thus, it effectively segments cracks from large-scale MLS point clouds.
Finally, the sensitivity of the CP with TC was also tested. Figure 16a shows the changing trends of the value of TC with CP; where CP increases, TC increases. When CP reaches 6%, TC approaches saturation. We did not test for cases greater than 10%, as crack points rarely account for more than 10% of the surface points on large roadways.
Although valuable (but not universal), interference-related noise caused by road markings, etc., that accounted for a larger proportion of the point cloud than cracks, is greatly reduced after setting CP (not exceeding 10%). It is also evident that when CP exceeds 4% (Figure 16b), the extracted crack curve is almost complete. As CP continues to increase, the noise increases. It is preferable to infer that the real proportion of cracks of the test data was approximately 4%. The results show that presetting the crack ratio is effective for eliminating interference.
The proposed method uses the regular grid indexes to transform discrete MLS point clouds into structures that are similar to multispectral images, from which the altitude and intensity information are integrated to detect the robust crack candidates. Image processing methods, such as morphological algorithms, can be used for post-processing. The execution time of building the Tgrid, detecting the valley points, and selecting points with low-intensity that are no higher than the road surface was 13.267 s, 8.9 s, and 78.21 s, respectively. The time cost of the proposed method is only one-thousandth that of Guan’s method of projecting 3D point clouds to 2D to construct geometric grid topology and using the tensor voting algorithm. Another benefit of the proposed method is that it does not lose the original accuracy of the point cloud and connects crack points using adjacency relation, both of which offer a sound basis for detecting cracks.
To summarize, the proposed method can successfully detect cracks with a width of not less than 1 cm using MLS data. As the geometric resolution of the state-of-the-art MLS system improves (the distance between adjacent points tx and ty direction can be compressed to 0.1–1 cm), detecting cracks from MLS data would have more potential and research value. In addition, the proposed structure makes point cloud processing easier and can introduce mature image processing methods into point cloud processing. More efforts are expected to be devoted to the detection of other on-road facilities, such as curbs, road markings, and manholes, from MLS data using the Tgrid method. Furthermore, the potential of Tgrid for improving the current deep learning point cloud frameworks is promising.

4. Conclusions

In previous studies, MLS data may not have been the first choice to evaluate pavement cracks owing to its geometric resolution and discrete distribution without connectivity. However, as the geometric resolution of the MLS system has been greatly improved, detecting cracks from MLS data has gained potential and significant research value. Based on the scanning angle or acquisition time of MLS point cloud, this study has losslessly transformed discrete MLS data into a “grid map” (Tgrid) through assigning each point cloud a grid index. Thus, numerous proven image-based methods of crack extraction can be implemented in point cloud processing. The approach was tested on a 220 m stretch of road with hundreds of cracks. This work has successfully constructed the “Tgrid map” for the test data and used the adjacency relation provided by this map to detect road surfaces, extract continuous road boundaries, and allow for image-based methods of identifying and quantifying pavement cracks in MLS point clouds. The results prove the feasibility and effectiveness of the proposed method on transverse, longitudinal, and oblique cracks, respectively. 3D crack skeletons and their shape parameters were extracted. The F1 scores of the quantity of correctly extracted transverse, longitudinal, and oblique cracks reached 96.55%, 87.09%, and 81.48%, respectively. Moreover, the average accuracies of the crack width and length were above 0.812 and 0.897. Experiments have proven that using the presented method, the MLS point cloud can be processed similarly to images. This finding is valuable for extraction of many on-road features, such as road points, curbs, road markings, and manholes, from MLS data. Furthermore, the proposed structure (Tgrid) can avoid the loss of information due to dimension reduction and has the potential to improve the existing deep learning framework for point clouds.

Author Contributions

Conceptualization, L.S. and Z.W.; methodology, M.Z. and L.S.; software, M.Z.; validation, Z.W.; investigation, M.Z. and L.S.; resources, L.S.; writing—original draft preparation, M.Z.; writing—review and editing, M.Z., D.H., and Z.W.; supervision, Z.W.; funding acquisition, L.S. and Z.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, Grant Nos. 41890854; 41372330. This research was also funded by the Transportation Research Project of Shaanxi Provincial Transport Department (SPTD) Nos. 16-01K, 18-06K and the Fundamental Research Funds for the Central Universities, Nos. CHD 300102269304.

Acknowledgments

Thanks are due to Tengfei Zhang, Chuanshuai Zhang for his data collection and technological help, and to Xue Wang for assistance with the experiments. We are grateful to anonymous reviewers for their thoughtful comments and suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Koch, C.; Georgieva, K.; Kasireddy, V.; Akinci, B.; Fieguth, P. A review on computer vision based defect detection and condition assessment of concrete and asphalt civil infrastructure. Adv. Eng. Inform. 2015, 29, 196–210. [Google Scholar] [CrossRef] [Green Version]
  2. Weng, X.; Huang, Y.; Wang, W. Segment-based pavement crack quantification. Autom. Constr. 2019, 105, 102819. [Google Scholar] [CrossRef]
  3. Hoang, N.D.; Nguyen, Q.L.; Tran, V.D. Automatic recognition of asphalt pavement cracks using metaheuristic optimized edge detection algorithms and convolution neural network. Autom. Constr. 2018, 94, 203–213. [Google Scholar]
  4. Yousaf, M.H.; Azhar, K.; Murtaza, F.; Hussain, F. Visual analysis of asphalt pavement for detection and localization of potholes. Adv. Eng. Inform. 2018, 38, 527–537. [Google Scholar] [CrossRef]
  5. Kim, J.; Lee, H.D. Development of New Automated Crack Measurement Algorithm to Analyze Laser Images of Pavement Surface. Diss. Theses Gradworks 2009, 112, 1231–1236. [Google Scholar]
  6. Jahanshahi, M.R.; Jazizadeh, F.; Masri, S.F.; Becerik-Gerber, B. Unsupervised Approach for Autonomous Pavement-Defect Detection and Quantification Using an Inexpensive Depth Sensor. J. Comput. Civ. Eng. 2013, 27, 743–754. [Google Scholar] [CrossRef]
  7. Du, Y.; Zhang, X.; Li, F.; Sun, L. Detection of Crack Growth in Asphalt Pavement Through Use of Infrared Imaging. Transp. Res. Rec. J. Transp. Res. Board 2017, 2645, 24–31. [Google Scholar] [CrossRef]
  8. Radopoulou, S.C.; Brilakis, I. Automated Detection of Multiple Pavement Defects. J. Comput. Civ. Eng. 2017, 31, 04016057. [Google Scholar] [CrossRef] [Green Version]
  9. Dorafshan, S.; Thomas, R.J.; Maguire, M. Comparison of deep convolutional neural networks and edge detectors for image-based crack detection in concrete. Constr. Build. Mater. 2018, 186, 1031–1045. [Google Scholar] [CrossRef]
  10. Fan, Z.; Li, C.; Chen, Y.; di Mascio, P.; Chen, X.; Zhu, G.; Loprencipe, G. Ensemble of Deep Convolutional Neural Networks for Automatic Pavement Crack Detection and Measurement. Coatings 2020, 10, 152. [Google Scholar] [CrossRef] [Green Version]
  11. Jeong, J.-H.; Jo, H.; Ditzler, G. Convolutional neural networks for pavement roughness assessment using calibration-free vehicle dynamics. Comput. Aided Civ. Infrastruct. Eng. 2020. [Google Scholar] [CrossRef]
  12. Ju, H.; Li, W.; Tighe, S.; Xu, Z.; Zhai, J. CrackU-net: A novel deep convolutional neural network for pixelwise pavement crack detection. Struct. Control. Heal. Monit. 2020, e2551. [Google Scholar] [CrossRef]
  13. Kalfarisi, R.; Wu, Z.Y.; Soh, K. Crack Detection and Segmentation Using Deep Learning with 3D Reality Mesh Model for Quantitative Assessment and Integrated Visualization. J. Comput. Civ. Eng. 2020, 34, 04020010. [Google Scholar] [CrossRef]
  14. Mei, Q.; Gul, M.; Azim, M.R. Densely connected deep neural network considering connectivity of pixels for automatic crack detection. Autom. Constr. 2020, 110, 103018. [Google Scholar] [CrossRef]
  15. Yang, F.; Zhang, L.; Yu, S.; Prokhorov, D.; Mei, X.; Ling, H. Feature Pyramid and Hierarchical Boosting Network for Pavement Crack Detection. IEEE Trans. Intell. Transp. Syst. 2020, 21, 1525–1535. [Google Scholar] [CrossRef] [Green Version]
  16. Hadjidemetriou, G.M.; Christodoulou, S.E. Vision-and Entropy-Based Detection of Distressed Areas for Integrated Pavement Condition Assessment. J. Comput. Civ. Eng. 2019, 33, 04019020. [Google Scholar] [CrossRef]
  17. Cha, Y.-J.; Choi, W.; Suh, G.; Mahmoudkhani, S.; Büyükztürk, O. Autonomous Structural Visual Inspection Using Region-Based Deep Learning for Detecting Multiple Damage Types. Comput. Aided Civ. Infrastruct. Eng. 2017. [Google Scholar] [CrossRef]
  18. Gopalakrishnan, K.; Khaitan, S.K.; Choudhary, A.; Agrawal, A. Deep Convolutional Neural Networks with transfer learning for computer vision-based data-driven pavement distress detection. Constr. Build. Mater. 2017, 157, 322–330. [Google Scholar] [CrossRef]
  19. Cha, Y.J.; Choi, W.; Buyukozturk, O. Deep Learning-Based Crack Damage Detection Using Convolutional Neural Networks. Comput. Civ. Infrastruct. Eng. 2017, 32, 361–378. [Google Scholar] [CrossRef]
  20. Zhang, D.; Zou, Q.; Lin, H.; Xu, X.; He, L.; Gui, R.; Li, Q. Automatic pavement defect detection using 3D laser profiling technology. Autom. Constr. 2018, 96, 350–365. [Google Scholar] [CrossRef]
  21. Tsai, Y.-C.; Wu, Y.-C.; Price, G. A Cost-Effective and Objective Full-Depth Patching Identification Method using 3D Sensing Technology with Automated Crack Detection and Classification. Transp. Res. Rec. J. Transp. Res. Board 2018, 2672, 50–58. [Google Scholar] [CrossRef]
  22. Pavement crack image acquisition methods and crack extraction algorithms: A review. J. Traffic Transp. Eng. 2019, 6, 535–556.
  23. Gui, R.; Xu, X.; Zhang, D.; Pu, F. Object-Based Crack Detection and Attribute Extraction From Laser-Scanning 3D Profile Data. IEEE Access 2019, 7, 172728–172743. [Google Scholar] [CrossRef]
  24. Xu, Z.-G.; Chen, Y.-L.; Li, J.-L.; Zhao, X.-M.; Pan, Y.; Wang, Z.-R.; Wei, N.; Song, H.-X. Research progress on automatic image processing technology for pavement distress. J. Traffic Transp. Eng. 2019, 19, 172–190. [Google Scholar]
  25. Li, L.; Wang, K.C.P. Bounding Box-Based Technique for Pavement Crack Classification and Measurement Using 1 mm 3D Laser Data. J. Comput. Civ. Eng. 2016, 30, 04016011. [Google Scholar] [CrossRef]
  26. Tsai, Y.C.J.; Li, F. Critical Assessment of Detecting Asphalt Pavement Cracks under Different Lighting and Low Intensity Contrast Conditions Using Emerging 3D Laser Technology. J. Transp. Eng. 2012, 138, 649–656. [Google Scholar] [CrossRef]
  27. Woo, S.; Yeo, H. Optimization of Pavement Inspection Schedule with Traffic Demand Prediction. Procedia Soc. Behav. Sci. 2016, 218, 95–103. [Google Scholar] [CrossRef] [Green Version]
  28. Chen, K.; Cheng, M.; Zhou, M.; Chen, X.; Chen, Y.; Jonathan, L.; Nie, H. Automated Object Extraction from MLS Data: A Survey. In Proceedings of the 2015 10th International Conference on Intelligent Systems and Knowledge Engineering, Taipei, Taiwan, 24–27 November 2015; pp. 331–334. [Google Scholar]
  29. Zhang, Z.; Cheng, M.; Chen, X.; Zhou, M.; Chen, Y.; Li, J.; Nie, H. Turning Mobile Laser Scanning Points Into 2d/3d On-Road Object Models: Current Status. In Proceedings of the 2015 IEEE International Geoscience and Remote Sensing Symposium, Milan, Italy, 26–31 July 2015; pp. 3524–3527. [Google Scholar]
  30. Guan, H.; Li, J.; Yu, Y.; Chapman, M.; Wang, H.; Wang, C.; Zhai, R. Iterative Tensor Voting for Pavement Crack Extraction Using Mobile Laser Scanning Data. IEEE Trans. Geosci. Remote Sens. 2015, 53, 1527–1537. [Google Scholar] [CrossRef]
  31. Guan, H.; Li, J.; Yu, Y.; Chapman, M.; Wang, C. Automated Road Information Extraction From Mobile Laser Scanning Data. IEEE Trans. Intell. Transp. Syst. 2015, 16, 194–205. [Google Scholar] [CrossRef]
  32. Chen, X.; Li, J. A Feasibility Study On Use Of Generic Mobile Laser Scanning System for Detecting Asphalt Pavement Cracks. In Proceedings of the XXIII Isprs Congress, Prague, Czech Republic, 12–19 July 2016; pp. 545–549. [Google Scholar]
  33. Yongtao, Y.; Li, J.; Haiyan, G.; Cheng, W. 3D crack skeleton extraction from mobile LiDAR point clouds. In Proceedings of the 2014 IEEE Geoscience and Remote Sensing Symposium, Québec City, QC, Canada, 13–18 July 2014. [Google Scholar]
  34. Zhong, M.; Sui, L.; Wang, Z.; Yang, X.; Zhang, C.; Chen, N. Recovering Missing Trajectory Data for Mobile Laser Scanning Systems. Remote Sens. 2020, 12, 899. [Google Scholar] [CrossRef] [Green Version]
  35. Wang, H.; Cai, Z.; Luo, H.; Cheng, W.; Li, J. Automatic road extraction from mobile laser scanning data. In Proceedings of the International Conference on Computer Vision in Remote Sensing, Xiamen, China, 16–18 December 2012. [Google Scholar]
  36. Wei, C.; Lichun, S.; Zhengchao, X.; Yu, L. Improved Zhang-Suen thinning algorithm in binary line drawing applications. In Proceedings of the 2012 International Conference on Systems and Informatics (ICSAI 2012), Yantai, China, 19–21 May 2012; pp. 1947–1950. [Google Scholar] [CrossRef]
  37. Yun, H.B.; Mokhtari, S.; Wu, L.L. Crack Recognition and Segmentation Using Morphological Image-Processing Techniques for Flexible Pavements. Transp. Res. Rec. J. Transp. Res. Board 2015, 2523, 115–124. [Google Scholar] [CrossRef]
Figure 1. The test data. (a) The overall view of test data, and (b) the divergent point intensity in two lanes.
Figure 1. The test data. (a) The overall view of test data, and (b) the divergent point intensity in two lanes.
Sensors 20 04198 g001
Figure 2. The connectivity provided by the arrangement of mobile laser scanning (MLS) point clouds. (a) The arrangement of MLS point clouds collected by a two-dimensional (2D) scanner with equal incremental angle. (b) Create a Time grid (Tgrid) for MLS point cloud by using the scanning line and scan angle information.
Figure 2. The connectivity provided by the arrangement of mobile laser scanning (MLS) point clouds. (a) The arrangement of MLS point clouds collected by a two-dimensional (2D) scanner with equal incremental angle. (b) Create a Time grid (Tgrid) for MLS point cloud by using the scanning line and scan angle information.
Sensors 20 04198 g002
Figure 3. The flow chart of construction of a Tgrid for MLS data.
Figure 3. The flow chart of construction of a Tgrid for MLS data.
Sensors 20 04198 g003
Figure 4. Construction of a Tgrid for MLS data through (a) the scan angle (θ), and (b) the acquisition time (T).
Figure 4. Construction of a Tgrid for MLS data through (a) the scan angle (θ), and (b) the acquisition time (T).
Sensors 20 04198 g004
Figure 5. Detection of road surface and boundary area in Tgrid. (a) The cross-section of a road segment, and (b) the searching process of a road boundary area.
Figure 5. Detection of road surface and boundary area in Tgrid. (a) The cross-section of a road segment, and (b) the searching process of a road boundary area.
Sensors 20 04198 g005
Figure 6. The flow chart of the method of detection of road surface points.
Figure 6. The flow chart of the method of detection of road surface points.
Sensors 20 04198 g006
Figure 7. Characteristics of valley points. (a) Height feature, and (b) gradient effect near crack candidates.
Figure 7. Characteristics of valley points. (a) Height feature, and (b) gradient effect near crack candidates.
Sensors 20 04198 g007
Figure 8. The flow of detection of the low-intensity points with height no higher than the road surface. (a) The background-differential map (GF), (b) the relevant truncation value (FB) by the given proportion of cracks’ point (TP), (c) the differential truncation map (GFC), (d) the segment threshold, TC, of the possible cracks, (e) the binary map of possible crack candidates (GP), and (f) the filter about valley effect.
Figure 8. The flow of detection of the low-intensity points with height no higher than the road surface. (a) The background-differential map (GF), (b) the relevant truncation value (FB) by the given proportion of cracks’ point (TP), (c) the differential truncation map (GFC), (d) the segment threshold, TC, of the possible cracks, (e) the binary map of possible crack candidates (GP), and (f) the filter about valley effect.
Sensors 20 04198 g008aSensors 20 04198 g008b
Figure 9. The structural elements of the morphological closing algorithm.
Figure 9. The structural elements of the morphological closing algorithm.
Sensors 20 04198 g009
Figure 10. The judgment of the crack direction.
Figure 10. The judgment of the crack direction.
Sensors 20 04198 g010
Figure 11. Edge points classification to calculate crack width, (a) in Tgrid, and (b) in geometrical space.
Figure 11. Edge points classification to calculate crack width, (a) in Tgrid, and (b) in geometrical space.
Sensors 20 04198 g011
Figure 12. The extracted point cloud of the road surface and time grid map rendered with intensity. (a) The extracted point cloud of the road surface, and (b) the intensity-based Tgrid map of the road surface.
Figure 12. The extracted point cloud of the road surface and time grid map rendered with intensity. (a) The extracted point cloud of the road surface, and (b) the intensity-based Tgrid map of the road surface.
Sensors 20 04198 g012
Figure 13. The extracted crack candidates in the Tgrid map and the results of crack craves. (a) The extracted crack candidates (in red) on the intensity-based Tgrid map of the road surface, (b) the remaining crack candidates after Morphological filtering, and (c) overlay of crack craves (in magenta) and road surface point cloud.
Figure 13. The extracted crack candidates in the Tgrid map and the results of crack craves. (a) The extracted crack candidates (in red) on the intensity-based Tgrid map of the road surface, (b) the remaining crack candidates after Morphological filtering, and (c) overlay of crack craves (in magenta) and road surface point cloud.
Sensors 20 04198 g013
Figure 14. The histogram of length and width compliance of cracks in 60 samples. (a) Compliance of the crack width, and (b) compliance of the crack length.
Figure 14. The histogram of length and width compliance of cracks in 60 samples. (a) Compliance of the crack width, and (b) compliance of the crack length.
Sensors 20 04198 g014
Figure 15. The detected crack candidates by using different methods. (a) Original point cloud, (b) the detected crack candidates by using the proposed method, (c) by using the maximum entropy sum method, and (d) by using Otsu’s thresholding method.
Figure 15. The detected crack candidates by using different methods. (a) Original point cloud, (b) the detected crack candidates by using the proposed method, (c) by using the maximum entropy sum method, and (d) by using Otsu’s thresholding method.
Sensors 20 04198 g015
Figure 16. The sensitivity test of the CP with TC. (a) The change of the value of TC with CP, and (b) crack candidates corresponding to the selected CP.
Figure 16. The sensitivity test of the CP with TC. (a) The change of the value of TC with CP, and (b) crack candidates corresponding to the selected CP.
Sensors 20 04198 g016aSensors 20 04198 g016b
Table 1. Parameters used in the process.
Table 1. Parameters used in the process.
ParametersValuesParametersValuesParametersValues
WD (m)0.50vth (°)10Rlwth3
α0.05CP (%)5Rs (m)0.20
hth (m)0.03D2 (m)0.04
hd (m)0.20Lth1 (m)0.25
Table 2. Quantitative assessment of the proposed method on three types of cracks.
Table 2. Quantitative assessment of the proposed method on three types of cracks.
Crack TypeActual QuantityPrecision(P)Recall (R)F1-Measure
Transverse crack9895.15%98.00%96.55%
Longitudinal crack4797.91%78.43%87.09%
Oblique crack4284.62%78.57%81.48%
Table 3. Quantitative evaluation of the width/length compliance of 3 types cracks.
Table 3. Quantitative evaluation of the width/length compliance of 3 types cracks.
Crack TypeNumber of SamplesMean (CW)Variance (VW)Mean (CL)Variance (VL)
Transverse crack200.8120.1490.9210.066
Longitudinal crack200.9100.0710.9350.089
Oblique crack200.8740.1310.8970.076

Share and Cite

MDPI and ACS Style

Zhong, M.; Sui, L.; Wang, Z.; Hu, D. Pavement Crack Detection from Mobile Laser Scanning Point Clouds Using a Time Grid. Sensors 2020, 20, 4198. https://doi.org/10.3390/s20154198

AMA Style

Zhong M, Sui L, Wang Z, Hu D. Pavement Crack Detection from Mobile Laser Scanning Point Clouds Using a Time Grid. Sensors. 2020; 20(15):4198. https://doi.org/10.3390/s20154198

Chicago/Turabian Style

Zhong, Mianqing, Lichun Sui, Zhihua Wang, and Dongming Hu. 2020. "Pavement Crack Detection from Mobile Laser Scanning Point Clouds Using a Time Grid" Sensors 20, no. 15: 4198. https://doi.org/10.3390/s20154198

APA Style

Zhong, M., Sui, L., Wang, Z., & Hu, D. (2020). Pavement Crack Detection from Mobile Laser Scanning Point Clouds Using a Time Grid. Sensors, 20(15), 4198. https://doi.org/10.3390/s20154198

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