Next Article in Journal
An Enhanced Hyper-Parameter Optimization of a Convolutional Neural Network Model for Leukemia Cancer Diagnosis in a Smart Healthcare System
Next Article in Special Issue
Differences in Driver Behavior between Manual and Automatic Turning of an Inverted Pendulum Vehicle
Previous Article in Journal
Improvement of Visible Photodetection of Chemical Vapor Deposition-Grown MoS2 Devices via Graphene/Au Contacts
Previous Article in Special Issue
Bimodal Extended Kalman Filter-Based Pedestrian Trajectory Prediction
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Urban Road Surface Discrimination by Tire-Road Noise Analysis and Data Clustering

by
Carlos Ramos-Romero
1,*,
César Asensio
2,*,
Ricardo Moreno
3 and
Guillermo de Arcas
2
1
Acoustics Research Centre, University of Salford, Manchester M5 4WT, UK
2
ETSI Sistemas de Telecomunicación, Departamento de Ingeniería Audiovisual y Comunicaciones, Grupo de Investigación en Instrumentación y Acústica Aplicada (I2A2), Universidad Politécnica de Madrid, 28031 Madrid, Spain
3
Institute for Chemical-Physical Processes of the Italian Research Council (CNR-IPCF), Via Giuseppe Moruzzi 1, 56124 Pisa, Italy
*
Authors to whom correspondence should be addressed.
Sensors 2022, 22(24), 9686; https://doi.org/10.3390/s22249686
Submission received: 14 November 2022 / Revised: 8 December 2022 / Accepted: 9 December 2022 / Published: 10 December 2022
(This article belongs to the Special Issue The Intelligent Sensing Technology of Transportation System)

Abstract

:
The surface condition of roadways has direct consequences on a wide range of processes related to the transportation technology, quality of road facilities, road safety, and traffic noise emissions. Methods developed for detection of road surface condition are crucial for maintenance and rehabilitation plans, also relevant for driving environment detection for autonomous transportation systems and e-mobility solutions. In this paper, the clustering of the tire-road noise emission features is proposed to detect the condition of the wheel tracks regions during naturalistic driving events. This acoustic-based methodology was applied in urban areas under nonstop real-life traffic conditions. Using the proposed method, it was possible to identify at least two groups of surface status on the inspected routes over the wheel-path interaction zone. The detection rate on urban zone reaches 75% for renewed lanes and 72% for distressed lanes.

1. Introduction

Due to continuous improvements in technologies applied to transportation systems, from autonomous driving to e-mobility, it is more common to find sensors embedded in vehicles that allow continuous scanning of the road environment. One of the crucial elements of this environment is the superficial pavement condition which is closely related to traffic safety and rolling noise emissions [1,2].
The Tire-Pavement Interaction Noise (TPIN) is a complex phenomenon which depends on a number of parameters, such as tire characteristics (e.g., tread pattern, inflation), driver influence (e.g., speed, acceleration), environmental conditions, surface contamination (e.g., wet, dry, dusty), and the parameter on which the present work will focus, the superficial asphalt condition [3,4,5,6].
Even though noise reduction properties of low-noise pavements decrease over time, resulting in heterogeneous tire-road noise generation along the route [7], the road inspection activities are mainly focused on roughness, skid resistance, and distress regarding functional and asphalt serviceability [8]. These conditions of the rolling surface are closely related to the noise generation from the wheel-path interaction zone, including both the megatexture and the macrotexture [9,10].
Traditionally, the identification and evaluation of asphalt defects has been performed by labour-consuming visual inspections. These activities have some subjective assessments due to experience and the assessor’s judgment of the catalogued defects. Other approaches, such as coring tests, are also available as mechanical techniques [11]. These methods are time-consuming, require further laboratory analysis, and are not easy to apply to large road sections.
Conversely, new remote sensing methodologies can address some of these shortcomings; for example, analysing large road sections is possible while handling large amounts of data. These methods include tools such as ground penetration radar, infrared thermography [12], laser scanning [13], image-based [14,15], vibration-based [16,17,18], and acoustic-based [19] methods. As these techniques are not mutually exclusive, more than one technique can be used simultaneously [10,20].
In particular, the acoustic-based methodology exploits the information from the rolling noise measurements, (i.e., Tire Pavement Interaction Noise (TPIN) [21], which is dominant over powertrain noise above the “crossover speed”, usually 35   km / h for combustion engine cars and lower speeds for electric vehicles [22].
The objective of this paper is to present an acoustic-based approach for unattended discrimination of the changes of the road surface conditions over the wheel-track path through the analysis of the tire-road noise and unsupervised machine learning (UL). Thus, the resulting clusters are presented accordingly the road sections of similar asphalt status. Although the obtained noise measurement data would also gather information about other aspects present during rolling, such as acoustical comfort inside the car, tire-pavement adhesion mechanism, vibrations, and aerodynamic phenomena, their identification is out of the scope of the present work. This study addresses the influences of driving conditions into the analysis, such as speed and acceleration of the vehicle, and several types and degrees of damage could be found on a typical route. However, the design of a dataset with all distress types for supervised classification could be a difficult and time-consuming task. Hence, UL algorithms will be included in superficial pavement monitoring of three road circuits. Table 1 shows a summary of the state-of-the-art research focused on road materiality studies based on the acoustic information contained in TPIN signals.
Although some standards for rolling noise measurement in near field have been implemented, such as CPX and OBSI, the data analysis from several sensors needed for their configuration can be time and computationally consuming tasks for other purposes added to policy noise level reports. Behind-The-Tire (BTT) measurement setup is a practical option, which allows to analyse the relative influence of pavement texture on tire-road noise through a simpler electroacoustic setup configuration [21,23]. The BTT technique has provided good performance in data acquisition stage for further identification of the actual condition of roads [5,24,25,26].
Table 1. Literature review: Road condition identification based on vibroacoustic data.
Table 1. Literature review: Road condition identification based on vibroacoustic data.
AuthorSensorsMethodClassificationFeature Extraction *Learning Approach
SoundVibrationRoughnessAsphalt TypeWet/dryAnomaliesMacro-mega
Texture
Road LabellingClassifier Algorithm *SupervisedUnsupervised
Masino et al. [27]x Tire cavity sound measurement x PSDSVM, ANNx
Ambrosini et al. [28] x Multichannel arrayx MFCCCNNx
Doğan [29]x Behind the tire (BTT) x MFCC, PSC, LPCANNx
Paulo et al. [30]x Close-Proximity method (CPX) x 1/3 OCTBayesianx
Kongrattanaprasert et al. [31] x Coast-by method x Power spectrumANNx
Alonso et al. [5] x BTT x 1/3 OCTSVMx
Abdic et al. [32]x BTT x Mel-Frequency scale.SVM, RNN-LSTMx
Pepe et al. [33]x Internal and external microphones x MFCCCNNx
Kalliris et al. [25]x On board single microphone x 1/1 OCTSVMx
Gueta and Sato [17] xOn board smartphone monitoring x x Peaks on the signal. WTSVM, LDA, KNNx
Ramos-Romero et al. [26]x BTT x x 1/3 OCTKNN. x
Safont et al. [34]xx10 channel sensor system x 256 featuresPCA, LDA, SVM, RFCx
David et al. [35]x Microphone pointed to the wheelxx 1/3 OCT—speedClustering x
Ganji et al. [36]x CPX x MFCCSVMx
Zhang et al. [37]x BTT x PCAStatistic modelx
Van Hauwermeiren et al. [38]xxOpportunistic sensing box xSound levels and acceleration. DAEDAEx
Del Pizzo et al. [9]x Tyre Cavity Microphone x x 1/3 OCTStatistic modelx
* Power spectral density (PSD). Power spectrum coefficients (PSC). Mel-frequency cepstrum coefficients (MFCC). Linear predictive coefficients (LPC). Principal components analysis (PCA). Wavelet transform (WT). Octave frequency spectrum (OCT). Support vector machine (SVM). Linear discriminant analysis (LDA). Convolutional neural networks (CNN). Artificial neural networks (ANN). Denoising autoencoder (DAE). Recurrent neural networks (RNN). Long-short term memory (LSTM). K-nearest neighbour (KNN). Principal component analysis (PCA). Random forest classifier (RFC).
Most of the literature in this field, exploiting the advantages of supervised machine learning in surface-asphalt detection using the frequency features of TPIN measurements, reports high accuracy on classification tasks, for instance the surface-road type, materiality class or wet-dry condition. The supervised method has been effective for these tasks, although it requires a labelled database with enough examples of all possible deterioration classes to identify them on the actual route. However, the road upper layer to be inspected might have many forms and degrees of deterioration, and several of them can even be unknown; especially, when the inspection is over long distances. In addition, the difficulty of collecting this complete database is the reason for the decision to use the unsupervised approach for monitoring the superficial asphalt condition along the wheel-path interaction zone.
This paper is organized as follows: Section 2 presents a detailed structure of the methodology, including data acquisition, signal processing, the UL algorithm, and geoprocessing. Then, in Section 3, the obtained results applying the methodology on urban roads are presented. Section 4 discusses the experimental results, the benefits of the application of this experimental approach, and the detected limitations. Finally, main conclusions of this work are detailed in Section 5.

2. Materials and Methods

The implemented method consists of four main steps, as depicted in Figure 1.
Firstly, the data are collected from three different types of input, such as the recordings of the TPIN signal, speed, and acceleration from the vehicle’s electronic control unit via the On-Board Diagnostics (OBD) port, and coordinates of the surveyed routes from a smartphone GPS. Afterwards, the dataset construction is performed through sound signal processing by frequency domain transformations, feature extraction, and driving condition dependences. The conformed dataset can be sliced by geographic areas or specific streets to deepen distress road identification.
Next, the similarities between the acoustic footprint data are clustered using UL algorithms. Finally, the available data from GPS were used to visually represent the detected differences in road conditions. Since the GPS of the mobile phone has an accuracy of 4.9 m, the geolocation data is not intended to determine the location of a specific type of deterioration, but rather a generalized inspection of the road surface.

2.1. Data Collecting

As far as possible during the data recording rounds, some constraints of the experiment were conserved to reduce variability factors and facilitate the subsequent interpretation of results. Only one class of tire tread pattern was installed on the experimental vehicle, which is depicted in Figure 2. This tread pattern is recommended for urban mobility in all seasons by the manufacturer (Pirelli-Cinturato P1™ Verde). The tire inflation pressure was checked before each experiment and kept constant at 2.2   bar . A single driver drove the routes to minimize variability in driving behaviour. In addition, the own vehicle’s mass was checked to ensure that it remained the same as much as possible. For this purpose, only the driver, regulatory equipment and data acquisition hardware were carried in the vehicle during the experiments. In addition, the measurement runs were driven with the fuel tank above its half-capacity whenever possible.
The TPIN data acquisition campaign was carried out during fall season under stable weather conditions, i.e., average air temperature maintained between 22–25 °C and humidity of 40–50%. Only dry asphalt was considered in the experiments.
Sound pressure signals were collected by using a unique instrumented diesel passenger car, equipped with two 1/2” type-1 microphones and BSWA-MA231 preamplifiers, according to the BTT technique. Although there is not an international standard for this simple instrumentation technique, it offers certain advantages for inspections. For applications in opportunistic scenarios, the sensor configuration is performed with the minimum disturbance to the vehicle structure, no additional trailer is required, and the microphones remain hidden during inspections.
The microphones were connected to an NI-9234 Data Acquisition (DAQ) device with sample rate of 51,200 Hz. The DAQ software is controlled by a portable PC on board. Although the measurement tasks were performed with two microphones, one on each rear wheel, the data captured by each microphone were processed independently. This means that the instrumentation setup simulates two vehicles of the same type making the same run with only one BTT microphone.
Simultaneously, an ELM327 interface connected at the car OBD port sends the driving parameter data to a smartphone via Bluetooth at 1 Hz. Because the PC was connected to the mobile network through the Wi-Fi smartphone portable hotspot, their clocks were matched to enable recordings using a single clock.

2.2. Dataset Design

The method seeks to recognize the differences between the sounds coming from different road conditions when the tire rolls over them. In this context, the information that makes the sound comparison between the signals is extracted from the TPIN registers and then constitutes a data set. The selected features are then expected to contain the relevant information of the input data.
The first step is to pass the audio data through a cut-off high-pass filter from 20 Hz, which removes the unrelated signal to rolling phenomenon at lower frequencies, such as sources of mechanical vibration and nonaudible acoustical data. Then, the audio data were split into 1-s nonoverlapped frames. This frame size allows matching between the OBD and GPS data sampled at 1 Hz with the corresponding TPIN frame trough date-time data. A Hamming window was applied to reduce the discontinuity effects at the boundaries of each chunk before the transformation in the frequency domain [39]. Each audio frame is processed by the Discrete Fourier Transform (DFT). The power amplitude is computed as the absolute value from the first half of the coefficients of DFT. Subsequently, the data conformed by DFT spectrum are further processed in the filter bank.
In general, literature reports that the sound radiation of TPIN mechanisms is prominent below 4000 Hz, and it is closed related with parameters of both the pavement (e.g., wavelength texture based on the aggregate particles in the mixture, the road condition and temperature) and the wheel (e.g., speed, torque, load, and inflation) [3,40]. Moreover, the rolling noise amplitude, at frequencies below 1000 Hz, increases with texture amplitude within the texture wavelength range of 10–500 mm. Indeed, above 1000 Hz, the noise amplitude decreases with texture amplitude within the texture wavelength range of 0.5–10 mm. [3]. Additionally, the changes in asphalt macrotexture are closely related to sound levels in the frequency range between 40 and 400 Hz at a specific speed [41]. Likewise, the sound intensity level is correlated with pavement surface texture at different frequency bands, from 315 to 2500 Hz [42].
Although TPIN noise characterisation is most commonly performed in the frequency domain through the 1/n-octave bands [5,25,26,29,35], there are alternatives for noise representation that are capable of handling the subjective impression of frequency, such as Mel’s triangular filter bank.
Indeed, the triangular filter bank was typically selected for audio processing in tasks related to the non-linear perception of sound by humans [39,43]. The reason for this type of filtering approach is that passengers inside the vehicle can perceive changes in the road surface by both auditory and vibrational stimuli, so this bank filtering is generally associated with acoustic comfort in the vehicle cabin [44].
In accordance with the frequency behaviour of the TPIN described above, the triangular filter bank was designed to extract the data at frequencies below 4000 Hz.
On the whole, each signal frame passes through a set of triangular filters ( T f = 50 ). The first 15 central frequencies ( f c ) are linearly spaced by 15 Hz ( f c = 65 Hz to f c = 260 Hz), and the next 35 central frequencies ( f c = 278.5 Hz to f c = 3089.4 Hz) are logarithmic spaced [45]. Figure 3 depicts the triangular filter bank representation, and the frequency bands are described in the Appendix A.
The OBD data allow the selection of only those signal frames that were recorded at more than 35 km/h, i.e., the data that continue the process are above the crossover speed. The detailed length of the studied roadways and corresponding observations [ n ] are listed in Table 2.

2.2.1. Corrections by Driving Conditions

To minimize the influences of driving conditions (speed and acceleration) on the frequency extracted sound features, a linear multivariate model was implemented. The data for the model were recorded under a wide speed range (35 to 60 km/h) on urban roads with the same vehicle, set of tires, and driver in order to minimize the data variability (see Table 2). The data used to derive the linear model was excluded from the data sets detailed in the results.
The speed s [km/h] and acceleration a [m/s2] of the vehicle registered during the trips are considered for corrections of noise level on each frequency band [35]. It is important to note that if speed and acceleration are taken as independent features, the clustering algorithm could return an incorrect detection of asphalt zones because it could group the zones by driving performance or roads constraints rather than by rolling noise and pavement conditions.
Consequently, every triangular band with sound pressure level in dB of a 1-s element [n] L T [ n ] f on the dataset is adjusted by s and a influence by the Equation (1), where L T [ n ] ( f ) is the corrected level and B s T f and B a T f are the coefficients for linear regression for speed and acceleration, respectively, at the reference speed s r e f = 70   km / h [35]. Although this s r e f value was applied to tire-road noise experiments related to road roughness [35] other values could also be considered [46].
L T [ n ] ( f ) = L T [ n ] f B s T f l o g 10 s [ n ] s r e f B a T f a [ n ]
Then, the feature space is made up of frequency bands with the coefficient of determination R 2 75 % and p v a l u e 0.005 . The bands with R 2 < 75 % or p v a l u e > 0.005 are rejected. Hence, only the bands from f c = 392.8   Hz to f c = 3089.4   Hz will be kept for the following stages of the experiment. These resulting 31 bands correspond to the range of useful frequencies indicated in the literature for tire/noise and pavement condition studies [42]. The selected triangular bands and the coefficients for linear regression are presented in Appendix A.
Furthermore, deterioration in certain areas of the pavement has been observed due to several punctual defects distributed over the surface. This surface irregularities generate impact noise events. To represent these noise dynamics during rolling, three additional overall-level features have been included for each observation [n] of the data. These are the peak level L P e a k [ n ] , equivalent continuous level L e q [ n ] , and the difference between them C F [ n ] = L P e a k [ n ] L e q [ n ] as crest factor. Finally, the main dataset is constructed with the featured TPIN by each element [n] as x = { L T [ n ] ( f ) , L P e a k [ n ] , L e q [ n ] , C F [ n ] } : x 34 .

2.2.2. Trip Segmentation

Once completed, the dataset can be filtered by geographic regions for further local and specific analysis. The available GPS data tracking makes it possible to select the amount of data for the next tasks by a geospatial query. Each zone has a very specific pavement structural capacity and corresponds to the same administration. The similarity in the initial asphalt mixture type within each geographical section has been assumed. This consideration would facilitate further analysis since the variability of the noise signal can reach 10 dB due to differences in pavement texture [7]. Therefore, data segmentation allows the selection of roads for a more specific analysis purpose, e.g., roads with similar usage and traffic density.

2.2.3. Pre-Processing and Feature Space Reduction

After data segmentation, the available data were pre-processed to ensure that each feature x contributed equally to the estimation of the parameters of the unsupervised model. Thus, the standardized feature x s t d was computed by Equation (2), with mean μ x and standard deviation σ x . These new scaled features are centered with mean 0 and standard deviation 1 [47].
x s t d = x μ x σ x
Next, the transformation of high-dimensional data into a meaningful representation of reduced dimensionality is included through dimensionality reduction, also called the feature reduction task. As a result, feature reduction facilitates and improves the cluster discrimination, classification, visualization and compression of high-dimensional data [48].
In this regard, unsupervised nonlinear dimensionality reduction was applied by the t-Distributed Stochastic Neighbor Embedding algorithm (t-SNE) introduced by Maaten and Hinton [49]. t-SNE attempts to maintain the local neighbourhood structure of input data points X = { x 1 , x 2 , , x n } d in low-dimensional space Y = { y 1 , y 2 , , y n } s , where s d , usually s = 2   o r   3 .
The above algorithm allows us to preserve the local structure of the data through pairwise similarity based on the Euclidean distance, while preserving much of the global structure of the data [49,50,51]. This distance-based feature reduction approach has provided better results for clustering tasks than other feature reduction methods, such as principal component analysis. Although this method was suggested for the graphical representation of features in a reduced dimensional space, it can also be employed with the clustering process, as it has been recently applied in fault detections and monitoring experiments [52]. Finally, the dataset reduced space results in y 2 , by t-SNE components equal to 2.

2.3. Unsupervised Learning: Cluster Model and Validity

Because there are no a priori labelled classes about the current condition of the inspected streets, UL may be able to explore the similarity and separability criteria among the observations in the data set. Then, two types of clustering techniques with different approaches were applied, but similar results were obtained: hierarchical clustering and probabilistic clustering. In many cases, the first attempt at grouping data set results in a clustering that may not be the most effective, so multiple clustering configurations must be studied [53].
The hierarchical approach does not apply a random initialization. Particularly, the data pooling algorithm starts with each observation as a cluster itself. These smaller clusters merge into larger ones by a series of successive fusions of the observations by minimizing the distance between clusters criteria [54]. The result of this agglomerative algorithm is a tree of clusters with distance relations. The number of clusters k or disjointed groups of data is obtained by cutting trees or dendrograms at a desired level [55].
Alternatively, the probabilistic Gaussian Mixture Models (GMMs) technique form ellipsoidal-shaped clusters based on the iterative expectation-maximization algorithm. It provides the basis for the Bayesian Gaussian Mixture Model (BGMM) [47,55,56].
The approximation of the minimum number of clusters could be established by the “elbow diagram”. The elbow method helps to determine the optimal number of clusters K by means of the sum of squared errors or inertia ( S S E ) function minimization [47]. The inflexion point “elbow” of the S S E vs. K plot shows the optimal number of clusters.
An example of the selection of the number of clusters k = 3 by both described clustering methods is depicted in Figure 4. Clustering algorithm results in a unique cluster label C k : k { 1 , 2 , , K } for each element of the dataset. This label is derived from the order in which the clusters have emerged.

2.4. Geo-Procesing of Results

The representation of all the passes over the inspected route is carried out by a basic geo-processing stage. In this process, the cluster type assigned to each element of the dataset was plotted on a map by GPS tracking data. However, after the cluster assignment task, the elements per cluster class might appear during normal driving with a certain level of randomness because the area of the asphalt on which the wheel contacts is not always the same. Moreover, defects may not be reached by the wheel in a single pass, especially in deteriorated areas with very local defects; in that case, the deteriorated area of the road upper layer is not registered just because the defect has not been passed over. Therefore, the possibility of detecting defects increases with the number of records belonging to the same street segment.
To address this drawback, the travelled road is first segmented, and then each element of the clustered dataset is related to the nearest piece of the road line. Road segments were set every 20 m for urban roads. This segment length allows us to include enough cluster elements in each segment when velocities are not constant. Finally, the segment label is assigned by the mode M o of the cluster type events C k | s e g m e n t associated with it, according to Equation (3). This information reveals a continuous report of the asphalt status along the inspected route. The basic geo-processing procedure was also applied in a previous work [26].
L a b e l s e g m e n t = M o { C k | s e g m e n t }
When a segment of route is classified as multimodal, the label legend not assigned will be present in the mapping report as a “n/a”. In the same way, the segments without available data (due to stop-car or low-speed events) will be present as “n/d”, from legend no data.

3. Results

The performance of the introduced acoustic-based methodology for surveying road condition was tested on urban roads (see Table 2).

3.1. Reference Route

The first experiment was carried out in the Reference-route. The benefit of applying the proposed methodology on this route is to evaluate the clustering performance on a controlled and known road environment. The road has two different types of superficial conditions (Figure 5) which includes the last renewed pavement (6 months before the data acquisition) conforming the main ring-shaped route, and the straight access to the ring where the distress sections are presented, such as longitudinal, transversal and alligator cracking.
For the feature reduction task, the t-SNE algorithm was fitted to respond closer to the known condition of the reference surface road, i.e., t-SNE components (n = 2) and (perplexity = 30). This last one is a parameter that means (loosely) how to balance attention between local and global aspects of the data [57]. Figure 5 shows the correspondence between the feature reduction step and the resulting clusters for pavement-condition zones. Figure 6 shows the clustered original features for the two resulting types of asphalt conditions.
Frequency features are showed in Figure 6 as the mean spectrum of the dataset for each cluster with its standard deviation. Temporal features are shown as box and whisker plots.
The cluster C 2 has the highest readings of noise levels per band, mainly between 896.0 Hz and 1553.4 Hz bands. Since the initial materiality of the route is considered the same, and the driving speed and acceleration are quiet constant into the selected data it can be assumed that the section assigned to cluster C 2 would correspond to the most deteriorated surface condition. Likewise, the differences in the means of the global noise level indicators between clusters C 2 and C 1 : Δ L P e a k   10   dB , Δ L r m s   7 dB and Δ C F   4   dB ; reinforce the premise that the route sections with distressed conditions is located on the route assigned to cluster C 2 .
Subsequently, the geographic processing of the instances assigned to each cluster was carried out through the stage described in Section 2.4. It was possible to visualise the sectors of the surveyed route according to the similar behaviour of the rolling noise footprint obtained from the clustering stage in Figure 5.

3.2. Urban Avenue

The second group of data refers to a single avenue in urban environment. It was known that maintenance work had been done previously, and a short section of this road had been completely repaved. Thus, the present experiment was motivated by the evaluation of the model’s ability to identify both distressed and renewed areas. For this experiment, the resulting number of clusters was k = 3 , as is depicted in Figure 7.
Subsequently, the instances assigned to each cluster could be displayed in their original domain, i.e., according to noise levels by frequency bands and by overall noise levels. Although the clustered noise spectrum does not show more than 3 dB difference between the average levels of each cluster (Figure 8) and the main values of the L P e a k and L e q values of C 1 and C 2 are closer to each other <3 dB, the C 3 is 6 dB lower than the other two clusters. The cluster C 3 corresponds to the data in driving deacceleration. The resulting clusters could be reported by the surveyed route line illustrated in Figure 7. In the mapped route report two main segments can be identified the one corresponding to cluster C 1 and a shorter one corresponding to cluster C 2 . The third cluster appears located in the street intersections, roundabout approximation zones and pedestrian crossings which could explain the data obtained during the car’s deacceleration. Cluster C 3 combines the lowest TPIN levels on feature plots. Nevertheless, this is not related completely with the superficial condition of the wheel-path interaction zone but with the driving conditions during the experiments.

3.3. Urban Street Circuit

The reason of this last experiment was to test the identification of the condition of the wheel-path interaction zone among several neighbouring streets. These streets, with similar characteristics, such as vehicle density flow and travel speed, do not necessarily have the same asphalt condition, but in general, the initial installed materiality could. Figure 9 shows the obtained results correspond to cluster number k = 3 .
Figure 10 shows the averages of the characteristics in frequency and time domain, as well as their standard deviation for each cluster assigned to the studied data. It is evident that instances assigned to cluster C 1 have the highest amplitudes, followed by cluster C 2 and cluster C 3 with lower amplitudes. In general, a difference of about 5 dB is observed between the band spectrum means and consecutive clusters. This slight difference in the frequency characteristics of clusters C 1 and C 2 could indicate that the track section assigned to cluster C 1 is in worse condition than the section assigned to cluster C 2 .
Besides, cluster C 3 shows the lowest amplitudes, and it could be associated with a quieter routing condition than the other clusters found. Its behaviour is like the third cluster of the previous experiment’s data set, i.e., the urban avenue. Hence, we can deduct that cluster C 3 is not related to the superficial condition of the wheel-path interaction zone.
Mapping report depicted in Figure 9, a marked zone assigned to cluster C 2 is dominant in almost a whole street of the inspected group. Because of the features on cluster C 2 are lower than C 1 , a better quality of asphalt is assumed. The greater section of the inspected group of streets belongs to C 1 , with the noisiest acoustic footprint. Finally, cluster C 3 describes a similar effect as the previous case of groups deaccelerating zones.
After, a further visual inspection of the studied street group was carried out. The area assigned to C 2 certainly showed signs of recent repaving, and zones assigned to   C 1 correspond to the old pavement.
A straightforward interpretation of the cluster correspondence of a road-surface condition is possible if the road surfaces along each route have similar materiality at the time of installation. This condition is assumed for the roads considered separately for each case study. From this point of view, the cluster with higher amplitude of acoustic footprint can be linked to the old areas, while new areas would belong to a lower amplitude cluster. Details of the visual inspection and the comparison with cluster estimation are presented at the end of the section in Table 3.
A comparison is made between interpretation of the resulting clusters and visual inspection. The term “estimation” is used since there are no previously assigned class labels to compare the detections. Nonetheless, it was possible to obtain the length of the road sections of each circuit l R S , which were visually labeled as “renewed” or “distressed”. Then, a ratio between l R S and the length assigned to each cluster l C is calculated. The length assigned to each cluster l C was obtained by counting the total length of segments assigned to each cluster in the geoprocessing step. In the same way, the total length of the “no-assigned” sections l n a and the total length of “no-data” sections l n d were counted. Each section of the route has a dominant cluster C k , which covers the greatest length of that section. Using this cluster, the condition of each section was rated.
The dominant clusters on each route section agree with the discriminated wheel-path condition obtained from the cluster interpretation and with the assigned condition by means of visual inspection. Then, the correct identification of each route section is estimated by the ratio between the dominant cluster length l C and the total distance of the route section, excluding the part that has not been considered on the datasets l R S l n d , in percentage. The percentage of the route that was not clearly identified and the percentage of the route that did not contribute data to the study are also presented, using l n a / l R S and l n d / l R S , respectively.
The results show that the surface discrimination for two pavement condition categories was obtained from the clusters within the reduced feature space. In the roadway of the University campus, the detection of both renewed and distressed sections shows more separable resulting clusters in the feature space (e.g., ΔLeq = 7.2 dB), and the a priori knowledge of the superficial roadway status allowed a high estimation of the real condition of the road. In particular, the renewed section length was detected up to 96% and the distressed section length up to 76%.
On the other hand, the system is less effective in road environments with greater variability. These are the cases of the urban avenue and the urban circuit, where the presence of a third cluster is related with external urban factors instead rolling noise. Therefore, this cluster was not considered for the estimation of the road surface condition. In addition, the classes proposed by the number of clusters are closer to each other into the feature space (i.e., ΔLeq < 3 dB). As a result, the discrimination rate of distressed areas between 60% and 72% was achieved for urban roads environment.

4. Discussion

The results of the automatic discrimination of the asphalt surface condition over the wheel-path interaction zone have been compared with the conventional visual inspection technique, which is performed on reference points of the road. In addition, local road evaluations were generalized to make them comparable with the results of the large length inspection attained from the proposed method. As a result, the discrimination rate presented in Table 3 can be underestimated
For all urban-roads study case, the corresponding C 3 appears in localized deceleration zones. It suggests that several low-noise data observations do not contain information purely related to tire-road noise. The speed limits and the car used for this traffic environment would make it difficult to record signals without contamination, such as ambient noise or phenomena related to sound reflections. This effect would be present in all clusters due to data acquired at speeds very close to the crossover speed (35 km/h) but mostly assigned to C 3 .
Due to the location of the microphones, BTT method depends on the wheel path, which is the contact area between the tire and the road during normal traffic. Therefore, the assigned group is not necessarily the same for all passes neither for both tires at the same time. However, with increasing passes, the detections show homogeneity across road segments, and the subsequent geoprocessing step allows smoothing the cluster assignment by a majority vote summary.
This acquiring system proposed allows to process signals from two microphone positions independently. This suggests that two cars with similar engine and tire characteristics and with only one instrumented tire, would provide similar results. However, the influence of driver and vehicle variability must be considered in the clustering interpretation. This can lead to collecting information collaboratively in the horizon.
  • This method, based on the clustering of the acoustic features of rolling noise, provides an unsupervised alternative for the discrimination of the asphalt surface status along the trajectory followed by a vehicle wheel. Notwithstanding, certain possible improvements, such as the selection of the vehicle and the placement of the microphone for data acquisition, should be considered in future works.
  • The placement of the microphone in the wheel housing produces signals that are not “purely” from the interaction between the tire and the road. In fact, there could be other types of sounds, both from the outside (when driving on busy roads) and from the vehicle itself (noise from the exhaust pipe or from the ventilation of the wheels themselves). The presence of these other sounds could negatively affect the classification process by increasing the background noise and masking the signal containing the contact surface information. However, most of these components are not sensitive to the type of road surface the vehicle is driving on. Therefore, the filtering, feature extraction, and dimensionality reduction processes have allowed to minimize their influence, as their acoustic fingerprints are separable. These non-overlapping classes were also observed in experiments with supervised classifiers with similar microphone placement [5,25,32,58]. Furthermore, the impact of these spurious sound phenomena is minimised by processing multiple observations at the same location using the geoprocessing step.
  • The selection of the speed of reference in the linear model could be improved in future works. Although the literature reports the linear relationship of the TPIN levels in dB and l o g s p e e d   [ km / h ] s p e e d r e f   [ km / h ] over a wide range of speeds, the speed of reference was considered constant throughout the experiments for homogeneous data processing. The analysis could be improved by applying different values of the speed reference based on standard recommendations as is proposed by ISO [46] and CNOSSOS-EU for traffic noise emissions modelling.
  • The asphalt discrimination rate in urban scenarios could be improved with the application of this acoustics-based method using electric vehicles, due to the lower crossing speed (<35 km/h) as reported in [22].
  • The detection of road surface quality by unsupervised learning has been evaluated by comparisons with applications of supervised classification metrics (i.e., accuracy) [35]. On the other hand, the present work proposes the “Estimation of road- section discrimination” which is based on the actual length of the road.
  • Comparison of the signals acquired by two or more microphones (e.g., one for each tire) could be included in future research steps. This would improve the detection of wear of the pavement, such as potholes, cracks, and bumps. A shorter-time window could also be included for impulsive noise events processing.

5. Conclusions

The main contribution of this paper is to present the performance of asphalt surface characterization based on the TPIN signal produced along the wheel path. This unsupervised technique allows the automatic discrimination of a small number of clusters related to predominant asphalt superficial condition.
Moreover, this method takes advantage of some improvement opportunities detected in previous related works, such as the inclusion of driving condition influences and the detection of the differences over the asphalt length without a catalogue of damages.
Three kinds of scenarios were considered to check the performance of the methodology such as a reference route, one urban avenue, and one urban street circuit. All these pavements possessed same characteristics of regular aggregate compositions, i.e., no one of them were neither open-graded pavements, rubber asphalts or poroelastic surfaces. These measurements were limited to weather conditions without rain and dry asphalt. The conclusions of the study are as follows:
  • The superficial condition of the studied roads is closely related to the rolling sound footprint and TPIN amplitudes in the frequency and time domains. These relations allow the interpretation of the clustering results.
  • An advantage of the application of UL over supervised techniques is the possibility of detecting areas with homogeneous rolling noise footprint without knowledge of the current road status. These localized zones are related to the homogeneous condition of the road status (deteriorated or not). The results were compared throughout further conventional visual inspections.
  • The implemented methodology has allowed the automatic and continuous discrimination of the state of the asphalt surface along the wheel trajectory. From these results, the surface discrimination of the wheel path on single lane roads can reach 92 % (i.e., the reference road and the urban street circuit). Multiple observations allow to evaluate better the TPIN from a narrower wheel track area.
  • Whereas in the case of the urban scenarios of roads with more than one lane, the discrimination rate decreases up to 57%. This because of the discrimination system must deal with different variables such as the speed limit, traffic flow, and a wider inspected area. Especially, when the vehicle changes lanes during each trip as it could happened during naturalistic driving behaviour.
  • The present acoustic-based method allows the inspection of road facilities with nonstop traffic inspections, non-destructive approach, and opportunistic scenario.
  • The mapping report contributes to pavement management through visual information. The surveyed areas producing different TPIN footprints assist in road maintenance planning, traffic noise mitigation activities, road condition warning reports.
  • In the present research phase, only corrections due to driving characteristics (speed and acceleration) were included. Future developments are also expected to incorporate corrections due to the variability of other conditions during driving, such as vehicle load, driver behaviour, tire inflation pressure, tire tread pattern, temperature, humidity, vehicle engine, pavement materiality, etc.
  • The technique could be improved for the detection of punctual defects such as potholes or manholes through refined time windowing in the signal processing and spatial resolution in the geoprocessing.
  • However, the consideration of these new conditions will surely imply complexity in the clustering interpretation.

Author Contributions

C.R.-R.; writing—original draft preparation, methodology, in situ recordings, signal processing, data curation. R.M.; Writing—review, discussion, and conclusions C.A.; conceptualization, data analysis, writing—review and editing G.d.A.; Supervision and conceptualization. All authors have read and agreed to the published version of the manuscript.

Funding

The data collection for this research was performed during the period 2018–2021 within the framework of the grant “Convocatoria Abierta 2017” (SENESCYT) by Ecuadorian Government, received by the author C. Ramos-Romero.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Features based on the triangular filters and the coefficients applied to the linear correction of the influence of speed and acceleration on TPIN levels.
Table A1. Features based on the triangular filters and the coefficients applied to the linear correction of the influence of speed and acceleration on TPIN levels.
FilterTriangular Filter
T f
Coefficientsp-ValuesR2FilterTriangular Filter
T f
Coefficientsp-ValuesR2
f L f c f U B s
Speed
B a
Accel.
SpeedAccel. f L f c f U B s
Speed
B a
Accel.
SpeedAccel.
150.065.080.038.72.20.00.076.726553.9593.3635.533.31.70.00.087.7
265.080.095.037.41.80.00.078.927593.3635.5680.833.61.70.00.086.7
380.095.0110.032.91.60.00.074.328635.5680.8729.232.01.90.00.085.3
495.0110.0125.025.21.10.00.060.529680.8729.2781.133.61.80.00.085.6
5110.0125.0140.022.20.90.00.057.830729.2781.1836.737.31.70.00.087.5
6125.0140.0155.022.50.90.00.059.731781.1836.7896.236.61.80.00.086.5
7140.0155.0170.023.50.90.00.060.832836.7896.2960.036.11.70.00.087.0
8155.0170.0185.025.60.90.00.065.533896.2960.01028.434.71.70.00.085.8
9170.0185.0200.027.20.80.00.065.034960.01028.41101.535.51.70.00.084.5
10185.0200.0215.027.20.90.00.067.7351028.41101.51179.936.01.80.00.081.5
11200.0215.0230.025.60.80.00.065.8361101.51179.91263.935.01.60.00.079.2
12215.0230.0245.020.51.20.00.057.0371179.91263.91353.934.31.80.00.079.2
13230.0245.0260.018.01.10.00.053.1381263.91353.91450.234.41.80.00.082.5
14245.0260.0278.516.91.20.00.048.7391353.91450.21553.433.81.70.00.081.5
15260.0278.5298.317.71.50.00.051.2401450.21553.41664.031.91.60.00.080.9
16278.5298.3319.620.11.60.00.061.2411553.41664.01782.432.01.70.00.082.2
17298.3319.6342.325.01.60.00.071.8421664.01782.41909.328.41.90.00.078.8
18319.6342.3366.730.41.80.00.072.9431782.41909.32045.231.42.00.00.079.7
19342.3366.7392.829.51.70.00.073.7441909.32045.22190.727.71.90.00.078.4
20366.7392.8420.730.41.90.00.083.4452045.22190.72346.627.52.10.00.077.0
21392.8420.7450.730.92.00.00.084.1462190.72346.62513.627.62.00.00.082.2
22420.7450.7482.731.61.80.00.084.3472346.62513.62692.527.61.80.00.083.0
23450.7482.7517.133.61.80.00.084.2482513.62692.52884.228.11.90.00.081.3
24482.7517.1553.933.21.60.00.086.5492692.52884.23089.427.91.90.00.079.6
25517.1553.9593.332.71.70.00.087.5502884.23089.43309.328.11.90.00.079.8
In, f L , f c , and f U are the lower, central, and upper frequencies of the triangular filters, respectively. The coefficients for Equation (1) by each triangular filter are B s for the speed and B a for the acceleration. Additionally, the p-values and the coefficient of determination R 2 ( % ) are presented. The highlighted filters (20 to 50) correspond to the selected 31 frequency features.

References

  1. Ma, Q.; Yang, H.; Mayhue, A.; Sun, Y.; Huang, Z.; Ma, Y. E-Scooter safety: The riding risk analysis based on mobile sensing data. Accid. Anal. Prev. 2021, 151, 105954. [Google Scholar] [CrossRef] [PubMed]
  2. Praticò, F.G.; Fedele, R.; Naumov, V.; Sauer, T. Detection and monitoring of bottom-up cracks in road pavement using a machine-learning approach. Algorithms 2020, 13, 81. [Google Scholar] [CrossRef] [Green Version]
  3. Sandberg, U.; Ejsmont, J.A. Tyre/road noise sources and generation mechanisms. In Tire/Road Noise Reference Book; INFORMEX: Evere, Belgium, 2002. [Google Scholar]
  4. Li, T. Influencing parameters on tire–pavement interaction noise: Review, experiments, and design considerations. Designs 2018, 2, 38. [Google Scholar] [CrossRef] [Green Version]
  5. Alonso, J.; López, J.M.; Pavón, I.; Recuero, M.; Asensio, C.; Arcas, G.d.; Bravo, A.E. On-board wet road surface identification using tyre/road noise and Support Vector Machines. Appl. Acoust. 2014, 76, 407–415. [Google Scholar] [CrossRef]
  6. Freitas, E.F.; Martins, F.F.; Oliveira, A.; Segundo, I.R.; Torres, H. Traffic noise and pavement distresses: Modelling and assessment of input parameters influence through data mining techniques. Appl. Acoust. 2018, 138, 147–155. [Google Scholar] [CrossRef] [Green Version]
  7. Mioduszewski, P. Low noise pavements in Poland. In Proceedings of the INTER-NOISE and NOISE-CON Congress and Conference Proceedings, New York, NY, USA, 19–22 August 2012. [Google Scholar]
  8. Radopoulou, S.C.; Brilakis, I. Automated detection of multiple pavement defects. J. Comput. Civ. Eng. 2017, 31. [Google Scholar] [CrossRef] [Green Version]
  9. Del Pizzo, L.G.; Bianco, F.; Moro, A.; Schiaffino, G.; Licitra, G. Relationship between tyre cavity noise and road surface characteristics on low-noise pavements. Transp. Res. Part D Transp. Environ. 2021, 98, 102971. [Google Scholar] [CrossRef]
  10. Wang, Z.; Zhan, J.; Duan, C.; Guan, X.; Zhong, Z.; Cao, Z. Road Surface Recognition Based on Vision and Tire Noise. In Proceedings of the 2021 5th CAA International Conference on Vehicular Control and Intelligence (CVCI), Tianjin, China, 29–31 October 2021. [Google Scholar]
  11. Dong, Z.; Ye, S.; Gao, Y.; Fang, G.; Zhang, X.; Xue, Z.; Zhang, T. Rapid Detection Methods for Asphalt Pavement Thicknesses and Defects by a Vehicle-Mounted Ground Penetrating Radar (GPR) System. Sensors 2016, 16, 2067. [Google Scholar] [CrossRef] [Green Version]
  12. Wang, S.; Zhao, S.; Al-Qadi, I.L. Continuous real-time monitoring of flexible pavement layer density and thickness using ground penetrating radar. NDT E Int. 2018, 100, 48–54. [Google Scholar] [CrossRef]
  13. 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]
  14. Zakeri, H.; Nejad, F.M.; Fahimifar, A. Image based techniques for crack detection, classification and quantification in asphalt pavement: A review. Arch. Comput. Methods Eng. 2017, 24, 935–977. [Google Scholar] [CrossRef]
  15. Diaby, I.; Germain, M.; Goïta, K. Evidential Data Fusion for Characterization of Pavement Surface Conditions during Winter Using a Multi-Sensor Approach. Sensors 2021, 21, 8218. [Google Scholar] [CrossRef]
  16. Kyriakou, C.; Christodoulou, S.E.; Dimitriou, L. Smartphone-based pothole detection utilizing artificial neural networks. J. Infrastruct. Syst. 2019, 25, 04019019. [Google Scholar] [CrossRef]
  17. Gueta, L.B.; Sato, A. Classifying road surface conditions using vibration signals. In Proceedings of the 2017 Asia-Pacific Signal and Information Processing Association Annual Summit and Conference (APSIPA ASC), Kuala Lumpur, Malaysia, 12–15 December 2017. [Google Scholar]
  18. Chen, C.; Seo, H.; Zhao, Y.; Chen, B.; Kim, J.; Choi, Y.; Bang, M. Pavement damage detection system using big data analysis of multiple sensor. In Proceedings of the International Conference on Smart Infrastructure and Construction 2019 (ICSIC) Driving Data-Informed Decision-Making; Institution of Civil Engineers (ICE): London, UK, 2019. [Google Scholar]
  19. Ganji, M.R.; Ghelmani, A.; Golroo, A.; Sheikhzadeh, H. A brief review on the application of sound in pavement monitoring and comparison of tire/road noise processing methods for pavement macrotexture assessment. Arch. Comput. Methods Eng. 2021, 28, 2977–3000. [Google Scholar] [CrossRef]
  20. Shamsabadi, S.S. Design and Implementation of Pavemon: A Gis Web-Based Pavement Monitoring System Based on Large Amounts of Heterogeneous Sensors Data; Northeastern University: Boston, MA, USA, 2015. [Google Scholar]
  21. Li, T. A state-of-the-art review of measurement techniques on tire–pavement interaction noise. Measurement 2018, 128, 325–351. [Google Scholar] [CrossRef]
  22. Verheijen, E.; Jabben, J. Effect of Electric Cars on Traffic Noise and Safety; National Institute for Public Health and the Environment of Netherlands: Bilthoven, The Netherlands, 2010. [Google Scholar]
  23. Osman, M.; May, D. Relative Influence of Pavement Texture and Tire Type on Pavement/Tire Noise; 0148-7191; SAE Technical Paper; SAE International: Warrendale, PA, USA, 1980. [Google Scholar]
  24. Zhang, Y.; McDaniel, J.G.; Wang, M.L. Pavement macrotexture measurement using tire/road noise. J. Civ. Struct. Health Monit. 2015, 5, 253–261. [Google Scholar] [CrossRef]
  25. Kalliris, M.; Kanarachos, S.; Kotsakis, R.; Haas, O.; Blundell, M. Machine learning algorithms for wet road surface detection using acoustic measurements. In Proceedings of the 2019 IEEE International Conference on Mechatronics (ICM), Ilmenau, Germany, 18–20 March 2019. [Google Scholar]
  26. Ramos-Romero, C.; León-Ríos, P.; Al-Hadithi, B.M.; Sigcha, L.; De Arcas, G.; Asensio, C. Identification and mapping of asphalt surface deterioration by tyre-pavement interaction noise measurement. Measurement 2019, 146, 718–727. [Google Scholar] [CrossRef]
  27. Masino, J.; Foitzik, M.-J.; Frey, M.; Gauterin, F. Pavement type and wear condition classification from tire cavity acoustic measurements with artificial neural networks. J. Acoust. Soc. Am. 2017, 141, 4220–4229. [Google Scholar] [CrossRef]
  28. Ambrosini, L.; Gabrielli, L.; Vesperini, F.; Squartini, S.; Cattani, L. Deep neural networks for road surface roughness classification from acoustic signals. In Proceedings of the 144th Audio Engineering Society Convention, Milan, Italy, 23–26 May 2018. [Google Scholar]
  29. Doğan, D. Road-types classification using audio signal processing and SVM method. In Proceedings of the 2017 25th Signal Processing and Communications Applications Conference (SIU), Antalya, Turkey, 15–18 May 2017. [Google Scholar]
  30. Paulo, J.P.; Coelho, J.B.; Figueiredo, M.A. Statistical classification of road pavements using near field vehicle rolling noise measurements. J. Acoust. Soc. Am. 2010, 128, 1747–1754. [Google Scholar] [CrossRef]
  31. Kongrattanaprasert, W.; Nomura, H.; Kamakura, T.; Ueda, K. Automatic detection of road surface conditions using tire noise from vehicles. IEICE Tech. Rep. 2009, 108, 55–60. [Google Scholar]
  32. Abdić, I.; Fridman, L.; Brown, D.E.; Angell, W.; Reimer, B.; Marchi, E.; Schuller, B. Detecting road surface wetness from audio: A deep learning approach. In Proceedings of the 2016 23rd International Conference on Pattern Recognition (ICPR), Cancun, Mexico, 4–8 December 2016. [Google Scholar]
  33. Pepe, G.; Gabrielli, L.; Ambrosini, L.; Squartini, S.; Cattani, L. Detecting road surface wetness using microphones and convolutional neural networks. In Proceedings of the 146th Audio Engineering Society Convention, Dublin, Ireland, 20–23 March 2019. [Google Scholar]
  34. Safont, G.; Salazar, A.; Rodriguez, A.; Vergara, L. Multichannel Signal Processing for Road Surface Identification. In Proceedings of the ICASSP 2020–2020 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP), Barcelona, Spain, 4–8 May 2020. [Google Scholar]
  35. David, J.; De Pessemier, T.; Dekoninck, L.; De Coensel, B.; Joseph, W.; Botteldooren, D.; Martens, L. Detection of road pavement quality using statistical clustering methods. J. Intell. Inf. Syst. 2020, 54, 483–499. [Google Scholar] [CrossRef]
  36. Ganji, M.R.; Ghelmani, A.; Golroo, A.; Sheikhzadeh, H. Mean texture depth measurement with an acoustical-based apparatus using cepstral signal processing and support vector machine. Appl. Acoust. 2020, 161, 107168. [Google Scholar] [CrossRef]
  37. Zhang, Y.; McDaniel, J.G.; Wang, M.L. Pavement macrotexture estimation using principal component analysis of tire/road noise. In Proceedings of the Nondestructive Characterization for Composite Materials, Aerospace Engineering, Civil Infrastructure, and Homeland Security 2014, San Diego, CA, USA, 10–13 March 2014. [Google Scholar]
  38. Van Hauwermeiren, W.; Filipan, K.; Botteldooren, D.; De Coensel, B. Opportunistic monitoring of pavements for noise labeling and mitigation with machine learning. Transp. Res. Part D Transp. Environ. 2021, 90, 102636. [Google Scholar] [CrossRef]
  39. Li, F.F.; Cox, T.J. Digital Signal Processing in Audio and Acoustical Engineering; CRC Press: Boca Raton, FL, USA, 2019. [Google Scholar]
  40. Flintsch, G.W.; De León, E.; McGhee, K.K.; AI-Qadi, I.L. Pavement surface macrotexture measurement and applications. Transp. Res. Rec. 2003, 1860, 168–177. [Google Scholar] [CrossRef]
  41. Saykin, V.V.; Zhang, Y.; Cao, Y.; Wang, M.L.; McDaniel, J.G. Pavement macrotexture monitoring through sound generated by a tire-pavement interaction. J. Eng. Mech. 2013, 139, 264–271. [Google Scholar] [CrossRef]
  42. Sakhaeifar, M.; Banihashemrad, A.; Liao, G.; Waller, B. Tyre–pavement interaction noise levels related to pavement surface characteristics. Road Mater. Pavement Des. 2018, 19, 1044–1056. [Google Scholar] [CrossRef]
  43. Combrinck, H.; Botha, E. On the Mel-Scaled Cepstrum; Department of Electrical and Electronic Engineering, University of Pretoria: Pretoria, South Africa, 1996. [Google Scholar]
  44. Junoh, A.K.; Mohd Nopiah, Z.; Wan Muhamad, W.Z.A.; Nor, M.J.M. Multi objective optimization of noise and vibration in passenger car cabin by using goal programming approach. In Advanced Materials Research; Trans Tech Publications: Wollerau, Switzerland, 2012; Volume 383–390, pp. 976–983. [Google Scholar]
  45. Slaney, M. Auditory toolbox: A matlab toolbox for auditory modeling work. Interval Res. Corp. Tech. Rep. 1998, 10, 1998. [Google Scholar]
  46. ISO 11819-2:2017; Acoustics—Measurement of the Influence of Road Surfaces on Traffic Noise—Part 2: The Close-Proximity Method. ISO: Geneva, Switzerland, 2017.
  47. Patel, E.; Kushwaha, D.S. Clustering cloud workloads: K-means vs gaussian mixture model. Procedia Comput. Sci. 2020, 171, 158–167. [Google Scholar] [CrossRef]
  48. Van Der Maaten, L.; Postma, E.; Van den Herik, J. Dimensionality reduction: A comparative. J Mach Learn Res 2009, 10, 13. [Google Scholar]
  49. Van der Maaten, L.; Hinton, G. Visualizing data using t-SNE. J. Mach. Learn. Res. 2008, 9, 2579–2605. [Google Scholar]
  50. Cao, Y.; Wang, L. Automatic selection of t-SNE perplexity. arXiv 2017, arXiv:1708.03229, preprint. [Google Scholar]
  51. Linderman, G.C.; Steinerberger, S. Clustering with t-SNE, provably. SIAM J. Math. Data Sci. 2019, 1, 313–332. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Agis, D.; Pozo, F. A frequency-based approach for the detection and classification of structural changes using t-SNE. Sensors 2019, 19, 5097. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Akman, O.; Comar, T.; Hrozencik, D.; Gonzales, J. Data clustering and self-organizing maps in biology. In Algebraic and Combinatorial Computational Biology; Elsevier: Amsterdam, The Netherlands, 2019. [Google Scholar]
  54. Murtagh, F.; Legendre, P. Ward’s hierarchical agglomerative clustering method: Which algorithms implement Ward’s criterion? J. Classif. 2014, 31, 274–295. [Google Scholar] [CrossRef] [Green Version]
  55. Halkidi, M.; Batistakis, Y.; Vazirgiannis, M. On clustering validation techniques. J. Intell. Inf. Syst. 2001, 17, 107–145. [Google Scholar] [CrossRef]
  56. Corduneanu, A.; Bishop, C.M. Variational Bayesian model selection for mixture distributions. In Proceedings of the Artificial intelligence and Statistics, Key West, FL, USA, 4–7 January 2001. [Google Scholar]
  57. Wattenberg, M.; Viégas, F.; Johnson, I. How to use t-SNE effectively. Distill 2016, 1, e2. [Google Scholar] [CrossRef]
  58. Ramos-Romero, C.; Cermeño, J.M.; Asensio, C. Shifts detection in the road surface condition through tyre/road noise analysis and pattern recognition approach. In Proceedings of the Euronoise 2021, Madeira, Portugal, 25–27 October 2021. [Google Scholar]
Figure 1. Methodological flowchart.
Figure 1. Methodological flowchart.
Sensors 22 09686 g001
Figure 2. Sensor’s set-up.
Figure 2. Sensor’s set-up.
Sensors 22 09686 g002
Figure 3. Dataset design.
Figure 3. Dataset design.
Sensors 22 09686 g003
Figure 4. Selection of number of clusters k = 3 by Dendrogram (left) and Elbow plot (right).
Figure 4. Selection of number of clusters k = 3 by Dendrogram (left) and Elbow plot (right).
Sensors 22 09686 g004
Figure 5. Reference-route data clustering results. t-SNE clusters formed (left). Geopositioned clusters (middle-right).
Figure 5. Reference-route data clustering results. t-SNE clusters formed (left). Geopositioned clusters (middle-right).
Sensors 22 09686 g005
Figure 6. Reference-route data clustering results. Features and driving conditions by formed clusters: Frequency domain (up); Peak level, Equivalent continuous sound level, Crest Factor, Speed and Acceleration (down).
Figure 6. Reference-route data clustering results. Features and driving conditions by formed clusters: Frequency domain (up); Peak level, Equivalent continuous sound level, Crest Factor, Speed and Acceleration (down).
Sensors 22 09686 g006
Figure 7. Urban avenue data clustering results. t-SNE clusters formed (left). Geopositioned clusters (right).
Figure 7. Urban avenue data clustering results. t-SNE clusters formed (left). Geopositioned clusters (right).
Sensors 22 09686 g007
Figure 8. Urban avenue data clustering results. Features and driving conditions by formed clusters: Frequency domain (up); Peak level, Equivalent continuous sound level, Crest Factor, Speed and Acceleration (down).
Figure 8. Urban avenue data clustering results. Features and driving conditions by formed clusters: Frequency domain (up); Peak level, Equivalent continuous sound level, Crest Factor, Speed and Acceleration (down).
Sensors 22 09686 g008
Figure 9. Urban Street circuit data clustering results. t-SNE clusters formed (left). Geopositioned clusters (right).
Figure 9. Urban Street circuit data clustering results. t-SNE clusters formed (left). Geopositioned clusters (right).
Sensors 22 09686 g009
Figure 10. Urban Street circuit. Features and driving conditions by formed clusters: Frequency domain (up); Peak level, Equivalent continuous sound level, Crest Factor, Speed and Acceleration (down).
Figure 10. Urban Street circuit. Features and driving conditions by formed clusters: Frequency domain (up); Peak level, Equivalent continuous sound level, Crest Factor, Speed and Acceleration (down).
Sensors 22 09686 g010
Table 2. Elements of dataset and the length of the inspected routes.
Table 2. Elements of dataset and the length of the inspected routes.
Roadway
ID
Circuit Length
[km]
PassesTravel Length Approx. [km]Dataset Elements [n]
1-s Readings
(Speed ≥ 35 km/h)
Data for lineal model 39.004504
Reference road 2.74410.961430
Urban avenue5.82421.522266
Urban street circuit2.81411.241998
Table 3. Discrimination rate of the surface of wheel-path interaction zone by clustering.
Table 3. Discrimination rate of the surface of wheel-path interaction zone by clustering.
Roadway
ID
Circuit Length [km]Road-Section by Visual
Inspection.
“Actual
Condition”
Road-Section Length l R S
[km]
Clustered Sections Lengths
after Geoprocessing
l C [km]
“Not Assigned” Section
[km]
“No Data” Section [km]Dominant
Cluster
C k
on
Road
Sections
Estimation of the Road-Section
Discrimination [%]
Percentage of Route
without Clear
Discrimination
[%]
Percentage of
No-Processed Route [%]
C1C2C3 l n a l n d l C l R S l n d 100 l n a l R S 100 l n d l R S 100
University
campus
2.54renewed1.541.420.00-0.000.12C192.210.007.79
distressed1.000.160.76-0.000.08C276.000.008.00
Urban avenue5.39renewed1.640.180.940.220.200.10C257.3212.206.10
distressed3.742.270.460.520.300.20C160.538.025.33
Urban street
circuit
2.81renewed0.480.040.360.060.020.00C275.004.170.00
distressed2.351.700.340.180.070.06C172.342.982.55
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ramos-Romero, C.; Asensio, C.; Moreno, R.; de Arcas, G. Urban Road Surface Discrimination by Tire-Road Noise Analysis and Data Clustering. Sensors 2022, 22, 9686. https://doi.org/10.3390/s22249686

AMA Style

Ramos-Romero C, Asensio C, Moreno R, de Arcas G. Urban Road Surface Discrimination by Tire-Road Noise Analysis and Data Clustering. Sensors. 2022; 22(24):9686. https://doi.org/10.3390/s22249686

Chicago/Turabian Style

Ramos-Romero, Carlos, César Asensio, Ricardo Moreno, and Guillermo de Arcas. 2022. "Urban Road Surface Discrimination by Tire-Road Noise Analysis and Data Clustering" Sensors 22, no. 24: 9686. https://doi.org/10.3390/s22249686

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