Abstract
The continuous monitoring and real-time feedback on the uniformity of the mixture after paving and compaction is important to guide the construction and improve the life of the asphalt pavement. Thus, in this study, the latest digital image technology was used to collect images of asphalt pavement after paving and compaction according to the set parameters. The images collected were preprocessed in steps such as binary grayscale, filtering, histogram, and small particle filtering, which established the method for calculating the construction depth of the digital images of pavement. Moreover, a method was proposed to evaluate the uniformity of asphalt pavement using the construction depth coefficient of variation. The coefficients of variation were computed based on 100 images of the dependent project and compared with the actual images. The results indicated that the construction depth coefficient of variation of the digital images obtained by the developed calculation procedure was in agreement with the actual situation and combined the correlation degree of the calculation results. Small areas of 42 are recommended as the division of the asphalt pavement compaction uniformity evaluation method based on construction depth variation, and the use of as the evaluation standard for the segregation of asphalt pavement is recommended. The method proposed is a new and rapid test to evaluate and detect the asphalt pavement uniformity.
1. Introduction
The segregation of asphalt mixtures is a common issue to asphalt concrete pavement surfaces. This segregation is accelerated by the dynamic load of vehicles, resulting in the development of potholes and pileups, which affect both the aesthetics and performance of the pavement in serious cases [1, 2]. Thus, the segregation of asphalt mixtures after paving is one of the root causes of local damage of asphalt concrete pavement [3].
Many studies have been conducted by domestic and international scholars to determine the segregation of asphalt mixture. The uniformity of asphalt mixtures mainly refers to the uniformity of distribution of coarse and fine aggregates [4]. Brock [5] stated that pavement segregation can lead to changes in mixture gradation and asphalt content, resulting in poor pavement structural performance. Kennedy et al. [6] analysed the phenomenon and mechanism of asphalt mixture segregation and proposed that it is related to the nonuniform dispersion of coarse and fine aggregates. Nondestructive technologies have been developed for the identification and evaluation of asphalt mixture homogeneity, such as thermal infrared imager, ground penetrating radar, and laser surface profiler [7–9]. Gilbert [10] used a thermal infrared imager to test the temperature distribution of the asphalt mixture during the construction process to reflect the degree of mixture segregation. The author found that the temperature of the asphalt mixture did not necessarily lead to aggregate segregation. Schmitt et al. [11] and Liu [12] used pavement ground-penetrating radar, which can determine the distribution of density of asphalt concrete pavement.
As a new interdisciplinary technology, the birth of digital image processing can be dated back to the 1960s and has become increasingly mature along with the development of computer technology and digital signal processing technology. In recent years, digital image processing technology has developed rapidly in the field of road engineering and has been applied to the quantitative evaluation of asphalt mixture uniformity research [13, 14]. Azari [15] verified the characteristics of aggregate distribution in computed tomography scan images of asphalt mixture specimens using the normal distribution method in statistics. If the distribution characteristics satisfied the standard normal distribution, the distribution was found to be uniform. Bruno et al. [16] analysed the gradation of aggregates by digital image acquisition of core samples after drilling and core samples of on-site asphalt pavement and extracted the aggregate particles in the core samples using a digital image method.
Considering the abovementioned studies, the application of digital image technology in asphalt mixtures has been mainly used to evaluate and analyse aggregate particle uniformity and asphalt mixture paving uniformity, whereas few studies have focused on the uniformity of asphalt mixture after paving and compaction and its correlation analysis with traditional methods. Therefore, in this study, for the asphalt mixture compaction stage, digital image processing technology was used to construct the relationship between the construction depth in the traditional sand paving method and the construction depth using digital image processing. The evaluation method of asphalt pavement compaction uniformity based on the construction depth variation was proposed to evaluate the compaction uniformity of asphalt pavement with the variation coefficient as the evaluation index.
2. Image Acquisition and Processing Technology
2.1. Image Acquisition
Acquisition of images is the first step in evaluating the uniformity of aggregate distribution in asphalt mixtures. The sketch of asphalt mixture image acquisition and calculation is shown in Figure 1. In this study, the image acquisition was tested on a sunny day, the metal oxide semiconductor COMS industrial camera in Figure 2 height from the asphalt mixture paving surface was 60 cm, and the front is vertical downward. The area of paving asphalt mixture captured by each image was measured to be . The image pixel resolution was according to the performance of the COMS industrial camera used for image acquisition at the construction site, shooting height, etc.
2.2. Image Preprocessing
To calculate the statistical analysis of the collected asphalt mixture images, the photos need to be processed. We collected photos of the middle pavement surface layer AC-20 asphalt mixture (The pavement asphalt layer has three layers, the lower layer is ATB-30, the middle layer is AC-20, and the upper layer is AC-16.) at Shiyan highway after paving and compaction. A total of 586 photos (shooting without overlap) were preprocessed. Figure 3 shows an example of the preprocessing, taking photo SY-001 as an example, where SY denotes the name of Shiyan highway, and 001 denotes the photo number.
(a) SY-001 original image
(b) Histogram of SY-001 after equalization
(c) SY-001 binarized postprocessed image
The acquired raw images were preprocessed by MATLAB software, and the processing steps were as follows: (1) greying the colour images using the rgb2gray function; (2) Use function medfilt2 for spatial filtering and denoising removal to improve the image quality; (3) histogram equalization to enhance the image contrast and improve the image clarity; (4) morphological processing and watershed segmentation of the histogram equalized images to identify and segment the adherent particles. After this process, the final binary image was obtained (Figure 3(b)).
3. Calculation Method of Construction Depth of Asphalt Pavement
3.1. Computational Model
As the surface of asphalt pavement is rough and uneven, light irradiation on the pavement will be diffusely reflected, and the light intensity of diffuse reflection will vary with the roughness of the middle pavement surface. Therefore, when the camera is used to capture the image of asphalt mixture, different light intensities are displayed differently in the image.
In this study, asphalt mixture compaction was the research object and the shooting occurred in natural light on a clear day. Studies reported that there is a linear relationship between the grayscale value and the construction depth (i.e., the greater the pavement construction depth, the smaller the grayscale value of the image) [17]. Therefore, digital image processing technology was used to calculate the depth of asphalt pavement construction. The average construction depth of the pavement image was calculated by reflecting the convexity of the asphalt pavement surface by the brightness and darkness of the image pixels, and the difference in height between the raised and concave pavement was reflected by the grayscale difference to construct a three-dimensional (3D) surface [18].
The pavement images acquired through the camera were red-green-blue images, which could be represented after preprocessing by MATLAB toolbox functions as an of a two-dimensional (2D) matrix. If , then a 3D surface could be constructed (i.e., the surface texture construction map of the asphalt pavement). Figures 4 and 5 show the 2D matrix of the asphalt pavement image and the 3D surface of the asphalt pavement surface image, respectively.
Based on the data from the 2D matrix, the pavement construction depth was calculated as where is the pixel value; are the horizontal and vertical coordinates corresponding to the pixel, respectively; is the volume of space enclosed by the pixel value of the surface construction; denotes the largest grey value in the grey matrix; is the grey value of the pixel point whose coordinates are ; the integration region is equivalent to the sand spreading range in the sand spreading method; is the volume enclosed by the spatial surface and plane of the pavement surface construction, which is equivalent to the sand spreading volume in the sand spreading method.
Digital image processing technology was used to calculate the digital construction depth of the pavement to characterize the roughness. The higher the value, the rougher the pavement. The formula for calculating the digital construction depth of the pavement is where is the integral region, calculated as the product of the number of rows and columns of the matrix.
The pavement surface construction depth needs to be corrected by the scale correction factor to obtain the final digital construction depth of the asphalt pavement surface. The calculation formulas are where and are the maximum and minimum values of all elements in the matrix , respectively; is the empirical coefficient and takes the value of 1.36.
3.2. Calculation Procedure for the Coefficient of Variation of the Construction Depth
3.2.1. Calculation Method
The asphalt pavement images can be divided into equal parts by 22, 32, 42, 52, and 102 as shown in Figure 6, and the mean depth of construction, variance between regions () and coefficient of variation (), is calculated for each equal part of the area. The calculation formulas are where is the average digital construction depth of each aliquot (mm); is the average digital construction depth of the image (mm); is the number of small regions divided, , 9, 16, 25, 100.
(a) 22 division method
(b) 102 division method
3.2.2. Procedure for Calculating Coefficients of Variation for Digital Image Construction Depth
In this study, MATLAB software was used to write a fast procedure for calculating the coefficient of variation of the construction depth from asphalt pavement.
MATLAB software was run, and the pavement digital construction depth discrepancy coefficient of variation calculation program, whose main interface is shown in Figure 7, was loaded. The main functions of this calculation program include create task, load image, automate, draw pavement texture construction map, and calculate and export, which are controlled by five buttons on the main interface.
4. Correlation Analysis of the Coefficient of Variation () and Aggregate Segregation
4.1. Correlation Analysis of Coefficient of Variation of Construction Depth from Digital Images
The acquired images of the medium surface AC-20 asphalt mixture compacted pavement were imported into the MATLAB calculation program. The images were segmented according to five different aliquot methods of 22, 32, 42, 52 and 102, and the coefficient of variation of the pavement digital image construction were calculated separately. The calculation results are shown in Table 1.
The calculated results of the values were compared to the corresponding asphalt pavement images by visual observation and were overall consistent with the expected assumption that the higher the value for the digital image construction of the asphalt pavement, the higher the likelihood of aggregate segregation. This tendency can be observed in Figures 8(a) and 8(b) for uniformity of mix compaction and inhomogeneity, respectively. The results of the computation of the construction depth coefficient of variation in Table 1 were as follows: Figure 8(a), , , , , ; Figure 8(b), , , , , . Thus, the values were corroborated with the visually observed uniformity of the asphalt pavement aggregate distribution.
(a) TXYS-11 undissociated image
(b) TXYS-95 dissociated image
4.2. Correlation Analysis of Depth Coefficient of Variation for Digital Image Construction with Different Equal Area Division Methods
The digital image construction depth coefficient of variation can be calculated according to five image segmentation methods of 22, 32, 42, 52, and 102. From Table 1, the values of the five segmentation methods of each image were different, but from the overall analysis of the 100 acquired asphalt pavement images, they showed a certain correlation. Therefore, a correlation analysis between two values was conducted, and a more stable and representative partitioning method was selected as the standard for calculating the coefficient of variation for asphalt pavements based on the results of the correlation analysis.
Taking the image divided by 22 and 32 as an example, the distribution of the digital image construction depth coefficients of variation and is shown in Figure 9. As observed, the same image divided by 22 and 32 division methods correspond to different construction depth coefficients of variation and , respectively. However, from the overall analysis of 100 asphalt pavement images, the tendency of the construction depth coefficients of variation and is consistent. A linear regression with as the independent variable and as the variable is shown in Figure 10. The linear regression equation is , and the squared correlation coefficient is . Thus, there is a good linear correlation between and .
Using the same method, the linear regression equations and square of the correlation coefficients were obtained by constructing a correlation analysis of the coefficients of variation , , , , and for the digital images, as shown in Table 2.
From Table 2, the square of the correlation coefficient between and , , , ranged between 0.83 and 0.91, whereas the squares of the correlation coefficients between , , , and the other four division methods ranged between 0.91–0.96, 0.89–0.98, 0.88–0.98, and 0.83–0.98, respectively. Thus, the correlations between the division methods of 22 and 102 and the other four division methods were relatively weak, whereas there was a strong correlation between the coefficients of variation from the construction depth of digital image calculated by aliquoting images by 32, 42, and 52 division methods. Considering that the more regions are divided, the slower their operation speed, and combined with the correlation degree of the calculation results, we determined that the best evaluation method of asphalt pavement compaction uniformity based on the construction depth variation consists of dividing the image into 42 small areas.
4.3. Evaluation Criteria Based on the Coefficient of Variation of the Construction Depth
According to the results of the coefficient of variation of the construction depth based on asphalt pavement images calculated using MATLAB, combined with the visual observation of the asphalt pavement aggregate distribution uniformity, it was found that when took a value close, equal, or greater than 8, there was a relatively evident segregation of the asphalt pavement aggregate distribution based on visual observation. As a result, the preliminary evaluation criterion of asphalt pavement segregation was set to .
For , the corresponding , , , and values were calculated according to the correlation relationship equation in Table 2 as 7.875, 7.907, 8.071, and 8.334, respectively. From the calculated 100 asphalt pavement images, the images with values of , , , and greater than or equal to 8 were selected. After screening, those with a coefficient of variation of are presented in Table 3.
Using as the boundary value of aggregate segregation for the following layer of AC-20 asphalt mixture pavement, six asphalt pavement images were screened, and the asphalt pavement corresponding to these six images was considered to exhibit aggregate segregation, as shown in Figure 11. By observation, these asphalt pavements exhibited a more evident coarse and fine aggregate separation, indicating that the use of for aggregate inhomogeneity determination was in accordance with the actual situation.
(a) TXYS-53
(b) TXYS-95
(c) TXYS-15
(d) TXYS-30
(e) TXYS-20
(f) TXYS-39
5. Conclusions
In this study, the compaction uniformity of asphalt pavement was evaluated based on the construction depth segregation for the compacted pavement of AC-20 asphalt mixture using digital image processing technology and MATLAB software programming. The main findings include (1)A calculation method for calculating the construction depth of digital images of pavements was established by characterizing the height difference between raised and depressed pavements in terms of the image grayscale difference, and a method for evaluating the uniformity of asphalt pavements using the construction depth coefficient of variation was proposed(2)A program was written in MATLAB to quickly calculate the digital construction depth coefficients of variation for pavements, and the coefficient of variation was also calculated for 100 images of the dependent project and compared with the actual images; it was found that the coefficients of variation calculated by the developed program were in accordance with the actual situation(3)By analysing the correlation between different division regions two by two, it was found that the correlation between 22 and 102 division methods was weak, whereas the division methods of 32, 42, and 52 exhibited strong correlations. Combined with the operation speed, a small area of 42 is recommended as the division method of asphalt pavement compaction uniformity evaluation based on the construction depth variation(4)By visually comparing and analysing the observed and calculated values, it is recommended to use as the evaluation criterion for asphalt pavement segregation
Data Availability
The data used and/or analysed during the current study are available from the corresponding author on reasonable request.
Conflicts of Interest
The author declares that there is no conflict of interest to report regarding the present study.
Acknowledgments
This work was supported by the Chongqing Natural Science Foundation Project (Grant no. cstc2020jcyj-msxmX0320) and the Science and Technology Research Program of Chongqing Municipal Education Commission (Grant No. KJQN202104001).