Assessing of the Road Pavement Roughness by Means of LiDAR Technology

The assessment of the road roughness conditions plays an important role to ensure the required performances related to road safety and ride comfort, furthermore providing a tool for pavement maintenance and rehabilitation planning. In this work, the authors compared the roughness index (International Roughness Index, IRI) derived from high speed inertial profilometer with two other roughness indices, one dynamic and one geometric computed on a digital elevation model (DEM) built by using mobile laser scanner (MLS) data. The MLS data were acquired on an extra-urban road section and interpolated on the nodes of a DEM with a curvilinear abscissa, coinciding with the global navigation satellite system (GNSS) track of the profilometer. To estimate the grid cell elevation, we applied two interpolation methods, ordinary kriging (OK) and inverse distance weighting (IDW), over the same data. The roughness values computed on the surface of the DEM showed a similar trend and a high correlation with those acquired by the profilometer, higher for the dynamic index than for the geometric index. The differences between the IRI values by profilometer and those computed on the DEM were small enough not to significantly affect the judgments on the analyzed sections. Moreover, the road sub-sections derived from profilometer measure that were classified as critical coincided with those derived from light detection and ranging (LiDAR) surveys. The proposed method can be used to perform a network-level analysis. In addition, to evaluate the effects of vibrations on human comfort, we input the DEMs into a dynamic simulation software in order to compute the vertical accelerations, as specified in the UNI ISO 2631 standard. The values obtained were in line and correlated with those inferred from the standard methodology for profilometer measures.


Introduction
Roughness is a major parameter for the pavement surface. It depends on surface irregularities and therefore affects drive quality and road safety [1]. Road maintenance is usually based on the assessment of roughness [2], as well as on vehicle delay costs and fuel consumption.
The regularity of the pavement affects vehicle vibration (produces accelerations that can be felt by the driver of the vehicle, consequently reducing ride comfort), driving speed, and tire wear. The conditions of a road pavement are defined by specific or global performance indices, which summarize its regularity by numerical values and evaluation scales. In particular, they can be geometric or dynamic and are computed along the longitudinal and transverse surface profiles of the road surface, while others are based on the punctual measurement of the various degradations.
The International Organization for Standardization ISO 13473-2:2002 [3] defines surface regularity as the deviation of the road surface from a true planar surface, and in part 3 (ISO 13473-3:2002 [4]), requirements for profilometer are specified. The American Society Such configuration, besides producing a partial survey of the road surface, provides a result that is strongly influenced by the trajectory.
Sayers and Karamihas [21] proved how a little variation in the vehicle's trajectory can produce a high variation (20%) on the IRI value computed in a 300 m long road segment. Another major issue is not being able to determine univocally the "start location and lateral position" [9], a prerequisite to ensure the repeatability of the measurement.
Transport agencies in recent years needed to produce better results/profits using fewer and fewer resources. One way they can increase productivity is through the use of new technologies. These include LiDAR (light detection and ranging) technology, which is currently one of the most interesting remote sensing techniques for infrastructure applications [22][23][24][25]. Laser scanners (LS), which are based on this technology, provide very accurate and high-resolution 3D data through dense point clouds.
Recent advances in sensor electronics and data processing make these technologies worth using. Assessments of the accuracy of an LS survey carried out on roads showed that this 3D data output is well suited for the analysis of both the regularity and geometry of roads where high accuracy in surveying is required [26].
The LS technique makes it possible to determine a "true" profile more effectively than traditional static methods. It also offers a further advantage because it provides data on the whole road surface and therefore gives the possibility to characterize the regularity of the road along the whole carriageway, instead of only along some trajectories [27][28][29]. A test showed that the spectral density of the profiles derived from LS is in line with that derived from the main traditional techniques; this highlights the reliability of the technique since it shows the same wavelength sensitivity obtained with the standard methods [27].
Profilometer measurements performed using multi-laser profilometers combined with high precision level have a 99% correlation with the data acquired by static LS [30]. A very useful device is the mobile laser scanner (MLS), a laser scanning system that allows for the acquisition of 3D data by means of one or more laser scanners mounted on a mobile platform [31].
MLS systems are significantly more efficient than static systems, terrestrial laser scanner (TLS). A freeway section about 80 km long can be surveyed in 3 h using an MLS, whereas using a TLS, the time required can exceed 120 working days [32]. Mobile LiDAR, a technology widely used since the end of the 20st century, when the first mobile LiDAR system became commercially available [33], has gained attention on applications mainly related to transports [34].
The point cloud obtained from MLS makes it possible a complete inventory of all the elements that make up the infrastructure, from the characterization of road geometry to that of signage and artwork [31]. Safety, efficiency, and cost benefits from MLS technology are highlighted [35].
Mendenhall [36] analyzed the efficiency of the MLS technique in terms of cost and acquisition time. The study was conducted on a 15-mile stretch of an urban street in the city of San Francisco; compared to standardized techniques, the cost savings was estimated to be between USD 200,000 and 300,000, with the time of acquisition reduced by 6 to 8 weeks, further reducing the management time by 4 weeks. The measurement was carried out by a team of operators on a vehicle equipped with MLS travelling on the road at operating speed. MLS scans repeated over time helped the initial design, estimation of the percentage of completion, project compliance, and as-built project drawings, significantly reducing costs and acquisition time compared to traditional methods [37].
There are many reasons that lead the authors to propose LiDAR technology as a support/alternative to the profilometer technique for the measurement of road surface roughness; for example, unlike most profilometers, LiDAR technology does not need to be limited in speed range for correct operating (for profilometers usually from 20 to 100 km/h). Some profilometers are not able to measure on very badly damaged roads or stone and dirt roads, and they are also very sensitive to humidity/temperature of the pavement. Another key aspect is that, in most cases, the equipped vehicle is able to measure the whole surface in one-way pass; these data are georeferenced with respect to an external reference system, allowing the repeatability of the measurement and allowing the unambiguous positioning of the points belonging to the surface to be compared with measurements made later for the study of the damaging trend.
In more recent times, several applications in the road field have been carried out to validate the use of LiDAR technology. Alhasan [38] acquired point clouds by TLS on different types of concrete Portland cement (PCC) and asphalt concrete (AC) pavements for the assessment of road roughness on unpaved roads; correlations with California bearing ratio (CBR) and dynamic cone penetration (DCP) have also been performed.
Alhasan [39] used the laser technique to evaluate roughness using a variety of methodologies: IRI, fast Fourier transform (FFT), and statistical analysis for the measurement of heights.
Kumar et al. [40] proposed an algorithm that allows users to assess roughness from a point cloud acquired through MLS. The aim of the authors was to identify a methodology able to provide fast, inexpensive, and complete information about regularity by determining the standard deviation of elevation values with respect to an interpolated planar surface.
A fast and automatic identification of localized surface defects is possible even using only the intensity of the reflected beam; morphological functions and special filters allow their identification [41].
Tran and Taweep [42] developed a segmentation algorithm, based on point cloud voxelization, able to estimate the roughness of road surfaces. The method has given excellent results, demonstrating the validity of operating with geometric methods. In addition, a few studies have been performed on potential improvements in inertial profilometers, adding new sensors recently made available such as non-contact wide-footprint and multi-point height sensors [43,44].
Over the years, interdisciplinarity has been a key factor for research, given the complexity of the various phenomena that occur in the characterization of the road/vehicle/user interaction. Liu and Herman [45] introduced a theoretical and detailed description of the vehicle-road interaction together with the analysis of the drive comfort. By using the quarter-car model and applying spectral density analysis, they studied the analytical expressions of dynamic indices, including the acceleration rate (jerk), the rectified mean speed (ARV), the ARS, and the IRI index. One year later, the same authors pointed out that the jerk index is one of the most important parameters for assessing the accessibility and regularity of roads in terms of comfort [46].
A few studies have recently been published to add to the computation of the IRI on the basis of riding comfort. Fuentes et al. [47] analyzed the subjective perception of the user about the level of comfort of a road pavement and the vertical accelerations through the analysis of deterministic and probabilistic models for the estimation of the pavement maintenance, specific for urban roads at low operating speed. Zhang et al. [48] introduced new IRI threshold values on the basis of riding comfort; three different evaluation techniques were used: a questionnaire, a vibration acceleration method, and a psychological and physiological index root mean square of the successive differences (RMSSD) method.
It is then easy to understand how there is already a need to introduce more complex methodologies that are able to simulate the dynamic response of a vehicle as faithfully as possible in order to evaluate its effects. The limitations were mainly due to measurement techniques that did not make it possible to address the problem in its real complexity.
The possibility of obtaining 3D data of the road surface on which to simulate "full car" models has only recently been applied [49] due to the rapid development and reliability of commercial dynamic simulation software. Studies have confirmed the good correspondence between the dynamic responses obtained from simulation and those measured by the instrumented vehicle. There are different studies on this subject; starting from the estimated stresses in the simulation environment, it is possible to deduce the effects of vibrations on the user's well-being without the aid of instrumented vehicles [50][51][52].
The aim of our work was to quantify the regularity of the whole road surface on the basis of a MLS survey, as opposed to what is usually done with the standardized method (profilometer), for which the judgment of regularity is based only on the measurement of linear profiles.
The indices of regularity resulting from the application of the two methodologies on a stretch of road more than 3.5 km long were compared, as well as the judgments of regularity that resulted.
The chosen parameters for the creation of the digital model of the pavement were analyzed and justified. In addition, an estimate of the drive comfort was made on the generated surface and the judgement was compared with that deduced from the IRI and the geometric index.

Test Case
The test area consisted of a section of a secondary suburban road, about 4 km long, from km 10 + 350 to km 13 + 870.
It consisted of a single carriageway with two lanes, one in each driving direction, 3.60 m wide. It was located near a large industrial district, which was why there were different levels of irregularities produced by industrial vehicles on it. The analyses were carried out on the north carriageway ( Figure 1). possible in order to evaluate its effects. The limitations were mainly due to measurement techniques that did not make it possible to address the problem in its real complexity.
The possibility of obtaining 3D data of the road surface on which to simulate "full car" models has only recently been applied [49] due to the rapid development and reliability of commercial dynamic simulation software. Studies have confirmed the good correspondence between the dynamic responses obtained from simulation and those measured by the instrumented vehicle. There are different studies on this subject; starting from the estimated stresses in the simulation environment, it is possible to deduce the effects of vibrations on the user's well-being without the aid of instrumented vehicles [50][51][52].
The aim of our work was to quantify the regularity of the whole road surface on the basis of a MLS survey, as opposed to what is usually done with the standardized method (profilometer), for which the judgment of regularity is based only on the measurement of linear profiles.
The indices of regularity resulting from the application of the two methodologies on a stretch of road more than 3.5 km long were compared, as well as the judgments of regularity that resulted.
The chosen parameters for the creation of the digital model of the pavement were analyzed and justified. In addition, an estimate of the drive comfort was made on the generated surface and the judgement was compared with that deduced from the IRI and the geometric index.

Test Case
The test area consisted of a section of a secondary suburban road, about 4 km long, from km 10 + 350 to km 13 + 870.
It consisted of a single carriageway with two lanes, one in each driving direction, 3.60 m wide. It was located near a large industrial district, which was why there were different levels of irregularities produced by industrial vehicles on it. The analyses were carried out on the north carriageway ( Figure 1). The LiDAR survey was carried out with an MLS Riegl VMX-450 mounted on the roof of a car. The MLS system was equipped with inertial and global navigation satellite system (GNSS) sensors, housed under an aerodynamic protection dome, in order to be able to perform moving scans and to georeference them in an external reference system (ETRF00). The system consists of two Riegl VQ-450 laser scanners inclined by 35°, each with a 360° field of view (type of configuration called "Butterfly"). Table 1 provides some details about the data acquisition and the point cloud obtained. The LiDAR survey was carried out with an MLS Riegl VMX-450 mounted on the roof of a car. The MLS system was equipped with inertial and global navigation satellite system (GNSS) sensors, housed under an aerodynamic protection dome, in order to be able to perform moving scans and to georeference them in an external reference system (ETRF00). The system consists of two Riegl VQ-450 laser scanners inclined by 35 • , each with a 360 • field of view (type of configuration called "Butterfly"). Table 1 provides some details about the data acquisition and the point cloud obtained. For comparing the roughness indices computed on the surface modelled using MLS data with the standardized index, we carried out a regularity measurement with a highspeed inertial profilometer. The profilometer bar, installed on the front of the instrumented vehicle, consisted of two lasers about 2 m apart; the elevation values were acquired at longitudinal intervals of 100 mm, and the set frequency was 62.5 kHz.
The profilometer was equipped with a GNSS positioning system to measure its trajectory. The profilometer measurement was made for a stretch of road between km 10 + 350 and km 13 + 870. IRI values provided by the profilometer were given for 10 m long segments, which led to 352 values being obtained. Table 2 shows the main parameters of measurement made with profilometer.

Methods
The road surface roughness analysis method we developed was based on the computation of two main roughness indices:

•
The International Roughness Index, computed using data by mobile laser scanner (IRI MLS ) and coming from the inertial profilometer (IRI PROF ).

•
The standard deviation of longitudinal roughness (σ), computed using data by mobile laser scanner (σ MLS ).
Data processing was carried out in three main steps, as follows: • The first step focused on the construction of a numerical model of the road surface built by interpolating the data acquired with mobile laser scanner. The built digital elevation model (DEM) was one especially designed for road pavements. Two different interpolation methods were implemented to estimate the elevation on the grid nodes in order to compare the results and choose the most efficient in terms of adherence to the surface and computational time.

Modelling of the Road Surface
The pavement surface modelling required the interpolation of the MLS data on the nodes of a grid. Since the road planimetry usually includes straight and curved lines and also the vehicle equipped with the profilometer had in turn variations with respect to the

Modelling of the Road Surface
The pavement surface modelling required the interpolation of the MLS data on the nodes of a grid. Since the road planimetry usually includes straight and curved lines and also the vehicle equipped with the profilometer had in turn variations with respect to the layout, we believe that a grid DEM with a north-south direction was not suitable for building a reliable model.
To overcome this issue, we designed and implemented a specific algorithm able to generate a curvilinear abscissa grid, hereafter called DEMc [24]. The direction of the abscissa can follow the axis of the road or the GNSS track of the profilometer. The curvilinear grid was organized in a matrix at different levels, composed by n rows equal to the number of nodes constituting the generic longitudinal profile and by c columns equal to the number of nodes constituting the generic cross-section. To each node was assigned an elevation value derived from interpolation.
To choose the grid step to be set, we used the formula proposed by Hengl [53], which relates the data density to the resolution of the regular grid. The minimum resolution value ρ can be estimated with the following formula: where D is the average density of the point cloud (number of points/dm 2 ). The elevation value to be attributed to the node was estimated using two different interpolators in order to evaluate the effect of the interpolation method on the regularity analysis: ordinary kriging (OK) and inverse distance weighted (IDW) [24]. The latter is a less rigorous method than kriging and consumes less CPU (central processing unit) time.
Kriging is a geo-statistical interpolation method that uses a variogram that depends on the spatial distribution of data rather than on actual values [54]. Ordinary kriging is the most widely used kriging method. IDW is a local deterministic interpolation technique that computes the value to be attributed to the node as a weighted average of the distance of the sample points in a given neighborhood [55]. It considers that the points closer to the node will have more influence and weights the sample points with the inverse of their distance from the node. In the IDW method, the main factor influencing the accuracy of the interpolator is the power value.
The best results are obtained using a power of 2 [56]. Powers greater than 4 produce slight differences in the resulting surface while increasing computational time, whereas a power equal to 1 produces a smoothed output [57]. Both interpolation methods were implemented in MATLAB to build the DEMc. The quality of the built DEMc was evaluated by estimating the "cloud to cloud distance", implemented in the open source software CloudCompare [58].
The default way to compute distances between two point clouds is the "nearest neighbor distance"-for each point of the compared cloud, CloudCompare searches the nearest point in the reference cloud and computes their (Euclidean) distance. When the nearest point in the reference cloud has been determined, the reference cloud (underlying) surface is locally modelled by fitting a mathematical model on the "nearest" point and several of its neighbors.
The distance from each point of the compared cloud to its nearest point in the reference cloud was replaced by the distance to this model. This was statistically more precise and less dependent on the cloud sampling. The local model used was a plan.
The criterion used is the Bayesian information criterion (BIC), a criterion for selecting the model from a finite set of models; the model with the lowest BIC was chosen. When adapting models, it is possible to increase the probability by adding parameters, but this can lead to over-adaptation. The BIC tries to solve this issue by introducing a penalty term for the number of parameters in the model. The whole process was implemented in MATLAB.

IRI Evaluation
In order to compute the IRI MLS , we implemented an algorithm based on quarter-car simulation (QCS) [7] in MATLAB. The mathematics of IRI MLS computation is where L is the length of the road section; t is the time of the simulation; and z s and z u are the vertical speeds of the sprung and un-sprung mass, respectively. IRI MLS was computed on each longitudinal profile of the road surface.
The σ MLS formula is as follows: where d i is the deviation of the ith elevation from a simple linear regression for an assigned base-length (3.5 m), and n is the number of elevation values in the base-length [13]. When implementing the algorithm, the weight of the σ depends on the length of the fixed interval (mm/m). Both indices were computed on any profile extracted from the DEMc of the road surface measured with MLS, having a curvilinear abscissa corresponding to the GNSS track of the profilometer. These indices were compared with the values derived from standardized profilometer measurements.
In order to take the uncertainty of the GNSS positioning system mounted on the profilometer into account, we implemented an algorithm that identifies the values numerically closest to the reference values (profilometer) in correspondence with the same milestone.
Such values are considered to be those with the lowest residue, computed as the absolute value of the differences between IRI PROF and IRI MLS . At the same time, the planimetric deviation in transverse direction between the track composed by these "minimum distance" points and the theoretical track given by the laser profiles on the interpolated surface is computed in order to verify if the deviation between the two tracks is in accordance with the accuracy of GNSS.
The uncertainty on the IRI computation is also related to the sampling step of the longitudinal profile; the maximum error peak is around 167 mm, the maximum allowable sampling interval is around 300 mm [59]. Between two consecutive points, the profile trend can be schematized as a linear function, a quadratic function, or a zero slope function (horizontal line). For intervals less than or equal to 50 mm, the difference between local models is not significant, it is significant for larger intervals [7]. A few studies show that the function that best approximates the trend is the linear function (the straight line connecting the points of the profile), as long as the inter-distance between the points does not exceed 300 mm for accurate measurements and 600 mm for less accurate measurements [7].
Chin [27] suggested an optimal sampling step for MLS data by demonstrating that a sampling range between 100 and 300 mm provides IRI values comparable (with a tolerance of ±5%) with those derived from standardized high-precision techniques (rod and level) and inertial profilometer. Following their indications, we used a linear model between two profile points and a spacing in accordance with the data density.
The analysis was carried out by computing the roughness values according to two different methods: IRI MLS and σ MLS on the DEMc surface, distinguishing the values of the right track from those of the left track. Each of these groups of values was compared with the recorded values on the corresponding traces.
The values of IRI MLS and σ MLS were computed over 10 m intervals in order to make a reliable comparison with the corresponding values derived by the profilometer (IRI PROF ). Since the road surface measured was more affected on the right side than on the left side, it was advisable for us to choose a small sampling interval in the IRI calculation (10 m) in order to highlight the localized effects. In particular, the sections were analyzed separately, dividing into two groups according to the value of the IRI (ASTM-E1926 [60])-the first group for IRI < 3.5 m/km and the second group for IRI > 3.5 m/km. The regularity computed for these sub-segments was compared with the regularity conditions of the corresponding segments on the DEMc, measured with the dynamic and geometric method.

IRI Evaluation of Ride Comfort
An assessment of the effects of vibrations on human comfort can be made by analyzing the vertical acceleration; in particular, according to the ISO 2631-1 [61], the frequencyweighted root mean square (RMS) accelerations (a wz ) is where w k,i is the value of the weight factor, and a RMS iz is the value of the RMS acceleration, computed using the power spectral density.
The vertical accelerations were computed using the dynamic simulation software package CarSim@ by Mechanical Simulation Corporation on the DEMc of the road surface.
Finally, the evaluation of comfort levels was assessed according to the values shown in Table 3.

Modelling the Road Surface
In order to build the DEMc using the kriging interpolation algorithm, we estimated the experimental semi-variogram; the lag was chosen so as not to have loss of detail and excessive smoothing of the variogram. A Gaussian' variogram model with nugget and absence of anisotropy was adapted to the experimental variogram model (Figure 3). Greater than 2 m/s 2 Extremely uncomfortable

Modelling the Road Surface
In order to build the DEMc using the kriging interpolation algorithm, we estimated the experimental semi-variogram; the lag was chosen so as not to have loss of detail and excessive smoothing of the variogram. A Gaussian' variogram model with nugget and absence of anisotropy was adapted to the experimental variogram model (Figure 3).   The grid step of the DEMc was set at 100 mm, on the basis of the density of the point clouds, and equally spaced with the data acquired by the profilometer. The residuals were computed as the distance between the local plane fitting the point cloud and the built DEMs (with kriging and IDW); the local plane model was interpolated using a number of points equal to the default number (6 points).
We used the BIC method (see Section 3.1) to find the probability distribution functions that best fit the residuals. For the DEM interpolated with ordinary kriging, the function that best interpolated the residuals was Weibull, while for IDW, it was the beta distribution ( Figure 4). The grid step of the DEMc was set at 100 mm, on the basis of the density of the point clouds, and equally spaced with the data acquired by the profilometer. The residuals were computed as the distance between the local plane fitting the point cloud and the built DEMs (with kriging and IDW); the local plane model was interpolated using a number of points equal to the default number (6 points).
We used the BIC method (see Section 3.1) to find the probability distribution functions that best fit the residuals. For the DEM interpolated with ordinary kriging, the function that best interpolated the residuals was Weibull, while for IDW, it was the beta distribution ( Figure 4).  Looking at the statistical results, it is possible to observe how the surfaces generated with ordinary kriging were similar to those generated with IDW. The processing times were very different, around 30 min for IDW and 10 h for kriging. Data processing was carried out with a Dell precision tower 3620 workstation, 32 Gb Ram, 500 Gb SSD storage; an Intel Xeon E3 1200 V5 processor; and a Radeon Pro WX5100 GPU. The choice of the interpolator to be used was related to the more or less homogeneous distribution of the data and was a function of the characteristics that are to be highlighted. The data distribution was therefore a key element in the choice of the interpolation method, regardless of its robustness. Tan and Xu [62], following tests on a set of experimental data, deduced that, to model homogeneously distributed data like our own, IDW and kriging give results, in terms of root mean square error and mean absolute error, of the same order of magnitude. In this work, we chose to use kriging, even if its computational cost is higher, as it is recognized in the literature as a more robust interpolator.

IRI Evaluation
The regularity analysis was carried out by computing the IRI values on a number of segments 10 m long; over the total length analyzed (3520 m), 352 values were obtained for each wheel track.
For each segment, the IRI MLS and σ MLS values on the DEMc were computed and compared with the corresponding IRI values derived from the inertial profilometer (IRI PROF ). The results for right and left tracks were analyzed separately. The boxplot of IRI differences (IRI PROF − IRI LSM and IRI PROF − σ LSM ) and the associated standard deviations (SD) are shown in the box in Figure 5. each wheel track.
For each segment, the IRIMLS and σMLS values on the DEMc were computed and compared with the corresponding IRI values derived from the inertial profilometer (IRIPROF). The results for right and left tracks were analyzed separately. The boxplot of IRI differences (IRIPROF − IRILSM and IRIPROF − σLSM) and the associated standard deviations (SD) are shown in the box in Figure 5. The correlation between the σMLS and IRIPROF was lower, given that the Pearson coefficients were 0.67 (left track) and 0.59 (right track) ( Figure 6). The diagrams in Figure 7 show the trend of the values computed on the left wheel track and the differences (in absolute value) between both IRIMLS and σMLS with IRIPROF The red and blue dots at the top of the figures show the segments (10 m long) with IRI greater than 3.5 m/km. The black line in the IRI diagrams highlights the values within the regularity threshold.
The diagrams in Figure 8 show results for the right wheel track. The correlation between the σ MLS and IRI PROF was lower, given that the Pearson coefficients were 0.67 (left track) and 0.59 (right track) ( Figure 6).
The diagrams in Figure 7 show the trend of the values computed on the left wheel track and the differences (in absolute value) between both IRI MLS and σ MLS with IRI PROF . The red and blue dots at the top of the figures show the segments (10 m long) with IRI greater than 3.5 m/km. The black line in the IRI diagrams highlights the values within the regularity threshold.
The diagrams in Figure 8 show results for the right wheel track.
Observing the box plots ( Figure 5), the correlation diagrams (Figure 6), and the IRI diagrams (Figures 7 and 8), one notices a large spread of the measured values, which highlights the occurrence of outliers in the dataset. Their occurrence may be due to

•
The vehicle with the profilometer bar being subjected to transverse oscillations during the motion, since the operator was not able to fully maintain the intended trajectory; skidding in the direction transverse to the motion significantly afflicted the values, resulting from the measurement of the profiles, mainly when the road segment was highly distressed (right trace, Figure 8). In detail, the IRI computed on 10 m long sec-tions was more affected by localized defects [20]; hence, the probability of highlighting one due to the skidding of the profilometric bar was quite high. • The differences observed were probably due to the dynamic response of the profilometer. The acquisition speed was variable (from 30 to 60 km/h), and thus large variations in acceleration probably led to errors in the measurement of the profiles and therefore in the computation of IRI [16].

•
The poor accuracy in trajectory measurement due to positioning errors of the profilometer's single-frequency GNSS receiver together with the imperfect knowledge of the offsets between trajectory and profilometer bar did not allow an accurate determination of laser profilometer traces [63].
Although we implemented an algorithm capable of accounting for transverse oscillation, the effects produced and highlighted in the second and third bullet points can certainly contribute to the generation of outliers, mainly concentrated in very distressed sections.
Observing the red and blue dots (Figures 7 and 8), one can notice how these tend to form homologous clusters only in correspondence of sections with adjacent dots (contiguous and greater than 10 m). The diagrams in Figure 7 show the trend of the values computed on the left wheel track and the differences (in absolute value) between both IRIMLS and σMLS with IRIPROF. The red and blue dots at the top of the figures show the segments (10 m long) with IRI greater than 3.5 m/km. The black line in the IRI diagrams highlights the values within the regularity threshold.
The diagrams in Figure 8 show results for the right wheel track.  Observing the box plots ( Figure 5), the correlation diagrams ( Figure 6), and the IRI diagrams (Figures 7 and 8), one notices a large spread of the measured values, which highlights the occurrence of outliers in the dataset. Their occurrence may be due to


The vehicle with the profilometer bar being subjected to transverse oscillations during the motion, since the operator was not able to fully maintain the intended trajectory; skidding in the direction transverse to the motion significantly afflicted the values, resulting from the measurement of the profiles, mainly when the road segment was highly distressed (right trace, Figure 8). In detail, the IRI computed on 10 m long sections was more affected by localized defects [20]; hence, the probability of highlighting one due to the skidding of the profilometric bar was quite high.  The differences observed were probably due to the dynamic response of the profilometer. The acquisition speed was variable (from 30 to 60 km/h), and thus large variations in acceleration probably led to errors in the measurement of the profiles and therefore in the computation of IRI [16].  The poor accuracy in trajectory measurement due to positioning errors of the profilometer's single-frequency GNSS receiver together with the imperfect knowledge of the offsets between trajectory and profilometer bar did not allow an accurate determination of laser profilometer traces [63].
Although we implemented an algorithm capable of accounting for transverse oscillation, the effects produced and highlighted in the second and third bullet points can The analysis of these clusters was used to subdivide the entire section in more subsections, in particular: • Sub-sections with grouped dots (grey background, 2L, 4L, and 6L for the left track and 2R, 4R, 6R, 8R, 10R, and 12R for the right track); these sub-sections identified sections with repetitive irregularities, not localized; • Sub-sections with isolated dots (1L, 3L, and 5L for the left track and 1R, 3R, 5R, 7R, 9R, and 11R for the right track); isolated dots represent sections that are probably characterized by the presence of localized distress (isolated defects) and therefore not very significant for network analysis.
Eighteen sub-sections with a length ranging between 60 and 1320 m were identified, 6 for the left wheel track and 12 for the right one. Among all 18, 9 were classified as having a good regularity with an average IRI PROF value of 1.79 m/km and 9 sections were classified as having low regularity with an average IRI value of 4.56 m/km (Tables 4 and 5, showing the profilometer values).  The IRI values shown in Figures 7 and 8 were computed on 10 m segments; the higher differences occurred in correspondence with very damaged surfaces, i.e., where isolated defects were also recorded.
The roughness value on each homogeneous section of the development between 60 and 1320 m was computed as an average of the values recorded on the relative 10 m segments; in this case, the differences of IRIs in the different methods were decreased, exceeding 2 m/km only on two sections, i.e., on a total length of 230 m, equal to 6.5% of the total development.
For each of these sections, the average IRI MLS and σ MLS values were computed and compared with the corresponding profilometer values (Table 6). In most cases, the IRI MLS and σ MLS values were lower than the profilometric values (IRI PROF ), and this sometimes led to a change in class. This happened for a total section length of 1280 m (18.2% of total length) for the dynamic index and 1700 m (24.1% of total length) for the σ MLS . It should be noted that the change of class mainly concerned the sections characterized by a few points above the threshold; especially in the case of the dynamic index, there was a significant reduction, greater than one unit, only for three sections (4.6% of the total length of the sample road). In the case of the geometric index, there was a more significant change in class for eight sections (16.1% of the total length of the road). It should be noted that the most troublesome sub-sections (IRI PROF > 3.5) were well identifiable also with the indices derived from laser scanner measurements; in most cases, they were the same identified with the standard method (IRI PROF ).
It should be noted that the results depended strongly on the development of each section analyzed-the major differences in regularity were found on the shorter sections with the most damaged surfaces ( Figure 5). These differences were, in fact, mainly present for developments of less than 250m (Table 6), and this can be explained by the effect of overlapping between localized distress and regularity. This points to the need for a standardization of the development of the sections to be measured, in accordance with the results of the experimental studies conducted by Sayers et al. [6].
However, the analysis of the correlation of the regularity values computed with the two indices (IRI MLS and σ MLS ) showed in these sections a very good correlation with the values obtained by the profilometer. Figure 9 shows the correlation diagrams together with the corresponding correlation and determination coefficients.

Evaluation of Ride Comfort
The evaluation of the ride comfort was carried out using the CarSim software package. To run the simulation, we used a Generic C-Class vehicle hatchback with Strut front and 5-Link rear suspension. C-Class vehicles include the following car models: Audi A3, Chevrolet Cruze, Ford Focus, and Opel Astra.
The DEM of the road surface was provided as input for the simulator in order to compute the weighted root mean square of the vertical accelerations acting on the vehicle. The computation of the accelerations was carried out at different speeds (40,60,80, and 100 km/h). The results are shown in the diagram ( Figure 10). It is possible to notice that there was a good correlation between the IRIPROF values and the vertical accelerations, meaning with the ride comfort.
At segments characterized by an IRI value close to or greater than 3.5 m/km, the vertical acceleration was greater than 0.8 m/s 2 , thus falling within the class defined as "uncomfortable". If a threshold value for the perception of discomfort equal to 0.4 m/s 2 [64] was set, the judgments were even more consistent with those derived from the IRIPROF values. This further highlights the potential of the methodology we proposed.
The influence of driving speed generally does not change the comfort class according

Evaluation of Ride Comfort
The evaluation of the ride comfort was carried out using the CarSim software package. To run the simulation, we used a Generic C-Class vehicle hatchback with Strut front and 5-Link rear suspension. C-Class vehicles include the following car models: Audi A3, Chevrolet Cruze, Ford Focus, and Opel Astra.
The DEM of the road surface was provided as input for the simulator in order to compute the weighted root mean square of the vertical accelerations acting on the vehicle. The computation of the accelerations was carried out at different speeds (40,60,80, and 100 km/h). The results are shown in the diagram ( Figure 10).

Evaluation of Ride Comfort
The evaluation of the ride comfort was carried out using the CarSim software package. To run the simulation, we used a Generic C-Class vehicle hatchback with Strut front and 5-Link rear suspension. C-Class vehicles include the following car models: Audi A3, Chevrolet Cruze, Ford Focus, and Opel Astra.
The DEM of the road surface was provided as input for the simulator in order to compute the weighted root mean square of the vertical accelerations acting on the vehicle. The computation of the accelerations was carried out at different speeds (40,60,80, and 100 km/h). The results are shown in the diagram ( Figure 10). It is possible to notice that there was a good correlation between the IRIPROF values and the vertical accelerations, meaning with the ride comfort.
At segments characterized by an IRI value close to or greater than 3.5 m/km, the vertical acceleration was greater than 0.8 m/s 2 , thus falling within the class defined as "uncomfortable". If a threshold value for the perception of discomfort equal to 0.4 m/s 2 [64] was set, the judgments were even more consistent with those derived from the IRIPROF It is possible to notice that there was a good correlation between the IRI PROF values and the vertical accelerations, meaning with the ride comfort.
At segments characterized by an IRI value close to or greater than 3.5 m/km, the vertical acceleration was greater than 0.8 m/s 2 , thus falling within the class defined as "uncomfortable". If a threshold value for the perception of discomfort equal to 0.4 m/s 2 [64] was set, the judgments were even more consistent with those derived from the IRI PROF values. This further highlights the potential of the methodology we proposed.
The influence of driving speed generally does not change the comfort class according to UNI ISO 2631-1 [61], although at higher speeds, the influence on comfort is stronger if the road is more damaged [64].
A comparison of the vertical accelerations recorded at 40 km/h with those at 100 km/h shows that the average increase was 23.1%. The same average increase reached 40% along the most damaged segments.

Conclusions
The availability of new technologies is focusing researchers' attention back on geometric road surveys. In the recent past, the dynamic analysis of the response of an instrumented vehicle has represented a technological alternative to overcome the drawbacks of geometric surveys (low productivity, high impact on traffic flow, etc.). The aim of our study was to compare regularity measurements derived by standardized techniques, namely, by profilometer, with those by MLS.
On the basis of the experimental tests we performed and the results we obtained, we can conclude that

•
The DEM with a curvilinear abscissa that was specifically designed and implemented for the road pavements allows us to model road sections several tens of kilometers long in order to run on the entire surface regularity analysis, without limitations due to the computational power, an issue experienced on grid models with northsouth direction.

•
The comparison between the indices derived by the profilometer and those computed on the DEM produced positive results; the correlations were good and proved the potential to quantify the regularity of the whole surface against the measurement on linear profiles only, but were unsuitable in terms of representing the surface in its real complexity. • The differences between the IRI values derived from profilometer and the values derived from DEM were small enough not to change and substantially affect the evaluation of the single sections, which is useful for planning maintenance works.

•
The critical sub-sections identified using the standardized method (profilometer) coincided, except for those where localized defects were present, with the critical sub-sections identified by the indices derived from MLS measurements. This meant that the method could be used for a quick analysis at the network level.

•
The values obtained, in terms of comfort, were in line and were strongly correlated with those deduced from standard methodology. This showed a double use of the implemented DEM; in addition to being used for an estimation of the regularity extended to the entire paved surface, it also allowed for the estimation of driving comfort through dynamic simulation.