6 January 2017 Mapping asphalt pavement aging and condition using multiple endmember spectral mixture analysis in Beijing, China
Author Affiliations +
J. of Applied Remote Sensing, 11(1), 016003 (2017). doi:10.1117/1.JRS.11.016003
Asphalt road reflectance spectra change as pavement ages. This provides the possibility for remote sensing to be used to monitor a change in asphalt pavement conditions. However, the relatively narrow geometry of roads and the relatively coarse spatial resolution of remotely sensed imagery result in mixtures between pavement and adjacent landcovers (e.g., vegetation, buildings, and soil), increasing uncertainties in spectral analysis. To overcome this problem, multiple endmember spectral mixture analysis (MESMA) was used to map the asphalt pavement condition using Worldview-2 satellite imagery in this study. Based on extensive field investigation and in situ measurements, aged asphalt pavements were categorized into four stages—preliminarily aged, moderately aged, heavily aged, and distressed. The spectral characteristics in the first three stages were further analyzed, and a MESMA unmixing analysis was conducted to map these three kinds of pavement conditions from the Worldview-2 image. The results showed that the road pavement conditions could be detected well and mapped with an overall accuracy of 81.71% and Kappa coefficient of 0.77. Finally, a quantitative assessment of the pavement conditions for each road segment in this study area was conducted to inform road maintenance management.
Pan, Zhang, Tian, Jin, Luo, and Yang: Mapping asphalt pavement aging and condition using multiple endmember spectral mixture analysis in Beijing, China



Asphalt roads play a crucial role in modern society and economy as an infrastructure for faster flow of people, goods, and services.1,2 Asphalt road surfaces are made of the mixture of bitumen and rocky components and have many benefits, including cost efficiency, comfort, durability, and up to 100% reusability.1 However, the quality of asphalt roads deteriorates gradually over time, which is primarily due to the effect of aging and deterioration of their surface.3 After being paved, new asphalt pavement will go through a series of aging processes caused by a variety of factors such as temperature, oxidation, loads, and water.4,5 Generally speaking, the aging process can be divided into four successive conditions: preliminarily aged (PA), moderately aged (MA), heavily aged (HA), and distressed. Furthermore, road maintenance engineers have found that the asphalt pavement with slight distresses tends to be destroyed completely in a short period if not timely and effectively maintained.2 Thus, it is imperative for road maintenance crews to assess the pavement condition and make timely maintenance plans to prolong road life.6,7 Some approaches have been previously developed to map road surface condition. For example, one common practice is extensive field inspections by engineers to evaluate several physical parameters, such as pavement condition index (PCI) and structural index (SI).5,8 In recent years, an advanced pavement management system (PMS) mounted on a vehicle has been used for road maintenance purposes. The PMS has a set of sophisticated subsystems, such as a global position system (GPS), geographic information system (GIS), and ground penetrating radar (GPR)910.11.12 and is capable of providing comprehensive information about road conditions. However, all of these approaches are labor intensive and quite time consuming.2,5 Therefore, only the major highways in some urban areas have been inspected, whereas the secondary highways have been either neglected or rarely inspected in any way.5

As an alternative, hyperspectral remote sensing has the advantages of wide spectral coverage and fine spectral resolution and has the potential of being used to map the asphalt pavement condition.7,1314.15 There has been evidence that pavement at different aging stages shows distinct spectral characteristics in reflectance and emissivity.5,16,17 Thus, it is possible to establish a quantitative relationship between road surface characteristics and their spectral signature.3,1819.20.21 However, due to the limitation of the spatial resolution of most remote sensing systems and the high spatial variability of the surface objects, mixed pixels are very common even in high-resolution images, especially for those captured in heterogeneous urban areas.22 Because roads are linear features imbedded in a background of other types of land covers, asphalt pavement is usually part of a mixed pixel, significantly increasing the uncertainty in image analysis of road condition.

To overcome this problem, we apply multiple endmember spectral mixture analysis (MESMA, Roberts23), a method to retrieve the fraction of all the materials in a mixed pixel, to mapping the asphalt pavement condition based on Worldview-2 satellite imagery in this study. Prior to the analysis, several in situ investigations and spectral measurements were conducted in the area. Based on the data collected in the field, a subpixel unmixing model was constructed to analyze the imagery and detect pavement conditions. Finally, a qualitative and quantitative assessment of the results was performed to validate the model.


Material and Methods


Aging Process of Asphalt Pavement

Based on previous studies and our field investigation, asphalt pavement shows a dramatic decline of bitumen proportion in the preliminarily aged condition due to the effects of volatilization, absorption, oxidation, and photochemical reaction on the road surface, resulting in a brittle and vulnerable pavement [Fig. 1(b)].5,24 Along with the further erosion of bitumen, moderately aged road surfaces expose their stone content [Fig. 1(c)]. The asphalt pavement then enters the heavily aged stage in which the rock pebbles are exposed at the road surface as relatively small granules due to continuous tearing as well as weathering effect [Fig. 1(d)]. During the aging processes, the deformability and structural strength of the road pavement are degraded. Eventually, a wide range of road pavement distresses occur in forms of linear cracks, alligator cracks, potholes, and rutting [Fig. 1(e)]. In order to quantify the degree of asphalt pavement aging, the Munsell Neutral Value Scale Card (MNVSC) was used to describe the color change of the asphalt pavement quantitatively during our field inspection.18,25 The MNVSC is a standard scale of neutral grays placed in sequence between black (0% reflectance) and white (100% reflectance) in 37 levels. Colors within each level are nondifferentiable to human eyes under standard viewing conditions. The standard color names and their corresponding values and spectral reflectance established by ISCC-NBS are illustrated in Table 1. According to the field investigation and visual interpretation, the aging conditions of asphalt can be categorized into three classes, preliminarily aged [N0.5/ to N4.25/], moderately aged (N4.25 to N6.75/], and heavily aged (N6.75/, N8.75/]. The range of (N8.75/, N9.5/] was not included in the classification to match the color of any kind of asphalt pavement. In addition, the pavement distresses were also excluded because the spectral features of distresses could not be appropriately or accurately characterized by MNVSC.

Fig. 1

The aging processes of asphalt pavement over time. (a) Fresh road pavement, (b) preliminarily aged pavement, (c) moderately aged pavement, (d) heavily aged pavement, and (e) pavement distresses (e.g. linear crack, alligator crack, pothole, etc.)


Table 1

ISCC-NBS color names, values, and spectral reflectance of MNVSC and the corresponding asphalt pavement conditions.

Color namesValueReflectanceAsphalt pavement condition
Black[N0.5/, N2.25/][0.6%, 3.8%]Preliminarily aged
Dark gray to black(N2.25/, N2.75/](3.8%, 5.5%]
Dark gray(N2.75/, N4.25/](5.5%, 13.7%]
Medium to dark gray(N4.25/, N4.75/](13.7%, 17.6%]Moderately aged
Medium gray(N4.75/, N6.25/](17.6%, 33.0%]
Medium to light gray(N6.25/, N6.75/](33.0%, 39.5%]
Light gray(N6.75/, N8.25/](39.5%, 63.6%]Heavily aged
White to light gray(N8.25/, N8.75/](63.6%, 73.4%]
White(N8.75/, N9.5/](73.4%, 90.0%]Not concerned


Study Area and Data Acquisition

In order to investigate the spectral features of the aged asphalt road pavement, and the possibility of mapping pavement conditions from remotely sensed data, several asphalt roads located in the southwest of Beijing, China, were selected (Fig. 2). The East Baiyang Road (marked in red) in the south of the image has three lanes in both directions. The Yangguang South Street (marked in yellow) and the East Ring Road (marked in green) are located in the west and the east of the image, respectively. The number of lanes of the Yangguang South Street is the same as that of the Baiyang East Road. The East Ring Road has four lanes in each direction. In the middle of the image are the two significantly narrower roads, the Liangxiang East Area No. 14 and No. 16 Roads, linking the Baiyang East Road and the Yangguang South Street. Both possess two or more lanes in their main roads.

Fig. 2

The study area and the road segments draped on the Worldview-2 image.


Although few distresses were found on these road pavements, it is obvious that each road was in a different aging condition. For instance, the Liangxiang East Area No. 14 Road was just repaved in early 2013 and is now preliminarily aged, whereas the Liangxiang East Area No. 16 Road and the Yangguang South Street appear to be in a worse condition because they were built in 2006 and have not been repaved since then. In contrast, the majority of the pavement in both the Baiyang East Road and the East Ring Road is moderately aged, showing the loss of bitumen and some outcrops of stone component.

In this study, Worldview-2 satellite images were chosen for mapping the asphalt pavement conditions (Fig. 2). The Worldview-2 sensor has eight bands spanning from visible to near-infrared wavelengths (from 0.45 to 1.04  μm) that are useful to identifying the information about the aging of asphalt road pavement. In principle, the WorldView-2 multispectral imagery has a resolution of 1.84 m. In this study, however, the WorldView-2 dataset provided by Kosmos Company of Remote Sensing was resampled by the data provider to 2-m spatial resolution and was captured on September 21, 2013. Subsequently, preprocessing was performed including radiometric calibration, and atmospheric and geometric correction to obtain the original land surface reflectance for further spectral analysis. First, based on the gain and offset data stored in the image metadata, WorldView-2 imagery was calibrated using ENVI 5.3 software package. The atmospheric correction was then conducted for the Worldview-2 image to retrieve land surface reflectance using the FLAASH module integrated with ENVI 5.3. In the module, the atmospheric model mid-latitude summer and aerosol model urban were applied. Additionally, given that accuracy assessment required an accurately geolocation match between the Worldview-2 image and field GPS points, image registration based on 11 GPS points was also conducted using a first-order polynomial and nearest neighbor resampling, resulting in an RMS error of 0.83. Finally, in order to avoid the impact of off-road areas, we delineated the road areas into five polygons based on information from field investigation and visual interpretation of the Worldview-2 image (Fig. 2).


Spectral Measurements with Analytical Spectral Devices Spectrometer

Field measurement of pavement spectra was conducted three times in the study area between April and July in 2015. An analytical spectral devices full range (ASD Field Spec-FR, ASD Corporation) was used to record spectral data. The instrument has three detectors covering the visible- and near-infrared (VNIR) and a short-wave infrared (SWIR1 and SWIR2) band, with a spectral sampling interval of 1.4 nm for the VNIR detector and 2.0 nm for the SWIR detectors. Field spectral measurement was conducted between Beijing local time 11:00 am and 13:00 pm under a clear sky, and the dark current was removed at every beginning of measurement. The optical fiber bundle collected the reflected radiation with a 25-deg conical field of view at 50 cm above the pavement, corresponding to a 22  cm×22  cm region. In addition to the ground targets, a white reference was also measured with the Spectralon Panel (Labsphere Inc., North Sutton, New Hampshire) to standardize and calculate the spectral reflectance of all the ground targets. The spectrum of every target including the white reference was recorded five times per measurement, and then the mean of five spectra was used in further processing. In addition to the asphalt pavement, the spectral information of other related objects including vegetation, cement sidewalks, white traffic lines, and bare soil was also collected, and ancillary information such as digital photos, pavement distress types, pavement health condition, and the asphalt color evaluated using the Munsell Neutral Value Scale (glossy collection) were collected and recorded through the field work. In total, a spectral library consisting 500 spectra was developed and used in the following spectral analysis.

The spectra were preprocessed including spectral smoothing and removal of the water absorption bands using the software package of ViewSpecPro 6.0.15, and a MATLAB Toolbox for spectral automatic postprocessing developed by NERC Field Spectroscopy Facility.26


Asphalt Pavement Condition Mapping Using MESMA

Some asphalt pavement pixels are covered by some other objects (e.g., vegetation, cement sidewalks, soil) even though the Worldview-2 image has a relatively fine spatial resolution of 2 m (Fig. 3). Spectral mixture analysis (SMA) was, therefore, performed to calculate the fraction of pavement from the mixed pixels. MESMA is based on the principle of linear spectral mixture analysis shown by the following equation:


where ρλ is the reflectance of the image acquired by remote sensor in band λ, and ρiλ is the reflectance of the ith endmember; fi represents the area fraction of the ith endmember in the mixed pixel; N is the number of endmembers, and ϵλ denotes the residuals produced by the model. The linear unmixing in Eq. (1) usually needs to meet the constraints in the following equation:



Fig. 3

A portion of the Worldview-2 false color image in the Baiyang East Road area.


The accuracy of the model can be evaluated using root mean square error (RMSE)


where M is the total number of the image bands. To account for the considerable within-class variability, MESMA allows the number and types of endmembers to vary from pixel to pixel. Generally speaking, the quality of the spectral mixture analysis results highly depends on the availability of representative endmembers.27 Endmembers can be derived from either the image directly (i.e., image endmembers28) by using an endmember extraction algorithm, or by referencing to a spectral library (i.e., reference endmembers) based on field spectral measurements.29 Considering the scale effect caused by the different influence of atmosphere and the bias of geographic location between the in situ measurements and image spectra, image endmembers with the same scale as the image are easier to associate with image features.30 For instance, vegetation spectra collected in the field do not easily “scale up” to match spectra measured by a remote sensor because field spectra are collected on the scale of individual leaves or branches and cannot adequately capture the multiple-scattering environment of a canopy.31 Therefore, image endmembers are extracted and utilized in the study and the ENVI-add on Viper tools were used to perform MESMA.27 The flowchart for the analysis is illustrated in Fig. 4.

Fig. 4

The flowchart of the MESMA with the Worldview-2 image.


Based on the GPS and other road attribute data taken in situ, a spectral library was first developed using a VIPER tool from several regions of interest (ROIs) on the Worldview-2 image. The library contains about 40 pixels for each representative surface type, covering the asphalt pavements in preliminarily aged, moderately aged, and heavily aged conditions, as well as other related objects including vegetation, cement sidewalk, white traffic line, bare soil, and shade. The library was optimized by analyzing the endmember average root mean square error (EAR), a metric proposed by Dennison and Roberts,32 to select the most representative endmembers with the lowest EAR value for each land cover class. EAR was calculated for each endmember by averaging the RMSE of the set of models that use that endmember to unmix the spectra belonging to the same land cover class


where A is the modeled spectra class, n is the number of spectra in class A, and Ai, Aj is the endmember in class A. The term n1 accounts for the endmember modeling itself, which produces a zero RMSE. EAR is equivalent to averaging the RMSE for all the models of a single endmember within the same land cover class.

Meanwhile, because too many or too few endmembers will increase model error,27 and the aging process is progressive in each pavement condition, five endmembers for each type of aged pavements and three nonpavement endmembers for every other class were extracted in this study. In addition, one shade endmember that we selected from the mean of the top three shade endmembers with the lowest EAR was used as well. As illustrated in Fig. 5, the spectra of the selected endmembers were used to build the spectral unmixing model. A coupled strategy is proposed in our study by combining the two-endmember model (EM2) with the three-endmember model (EM3). Namely, if the RSME by EM2 and that by EM3 is greater than σ, EM3 is adopted; otherwise, EM2 should be adopted. This strategy is simple but quite effective per our experiment.31,33

Fig. 5

The endmember spectra resample to Worldview bands and used in the spectral unmixing model.


Viper tools can provide a maximum of four spectral libraries including one for shade endmember. Each of these libraries refers to one type of endmember in a mixed pixel. Specifically, Viper tools define two shade options, photometric shade and nonphotometric shade. The former defaults to zero reflectance, whereas the latter enables a user to provide a shade spectrum from a spectral library. In this study, we chose the filtered shade endmember by EAR as the input. The EM2 models (Table 2) included 27 endmembers in the first library (namely, 15 of three kinds of aging pavement and 12 of four nonpavement objects) and the shade endmember in the second library. They were processed in a constrained mode with a set of recommended parameters by Roberts et al.34 (e.g., minimum and maximum fraction constraints of 0.05 and 1.05, shade constraint of 0.8, RMSE of 0.025, and residual constraint of 0.025). In addition, 180 EM3 models (Table 2) in the same constraint mode as above were applied to the image as well, in which the models were comprised of 15 asphalt pavement endmembers in the first library, 12 endmembers of other objects in the second library and one shade endmember in the third library. After generating the fraction image of each class of land cover (i.e., PA pavement, MA pavement, HA pavement, vegetation, sidewalks, traffic line, and soil), shade normalization following 31 and 33 was performed, that is, each nonshade fraction was divided by the sum of the nonshade fractions for that pixel. For those pixels modeled by EM2 models, the resulting shade normalized fraction would be 100%. Finally, for more effective visual interpretation for pavement in different aging degrees, a classification image for aging pavement distribution was generated based on the criteria: the value of a pixel would be set to 1 if the fraction of a pavement endmember in it is >0.5.

Table 2

The statistics of the EM2 and EM3 models in the experiment.

Model typeFirst endmemberSecond endmemberThird endmemberNumber of pixelsPercentage (%)
EM2PA pavementShadeNull30845.88
MA pavementShadeNull767814.63
HA pavementShadeNull1456827.76
Traffic linesShadeNull34686.61
EM3PA pavementVegetationShade10351.97
PA pavementTraffic linesShade3560.68
PA pavementSoilShade1520.29
PA pavementSidewalksShade2780.53
MA pavementVegetationShade18643.55
MA pavementTraffic linesShade3560.68
MA pavementSoilShade2940.56
MA pavementSidewalksShade5561.06
HA pavementVegetationShade27565.25
HA pavementTraffic linesShade7671.46
HA pavementSoilShade2650.51
HA pavementSidewalksShade8251.57


Accuracy Assessment Procedure

In order to evaluate the accuracy of the fraction image of every target, the Worldview-2 multispectral image was pan-sharpened with its corresponding 0.5-m panchromatic image, to generate a multispectral image of 0.5-m resolution for the study area. Then, the classification of pan-sharpened image was performed by visual interpretation, which was based on the priori knowledge from our three exhaustive field investigations in the study area. A classification image consisting of five classes (PA pavement, MA pavement, HA pavement, shade, and nonpavement land covers) was produced and visually inspected for accuracy, and errors were manually corrected. The visually classified image was then used to calculate the fraction of each class within each 4×4  pixel neighborhood, from which a 2-m fraction image was generated and used as reference in the following accuracy assessment. Afterward, a number of samples for each type of endmember was selected from the 2-m reference image and used to conduct linear regression analysis with the modeled fraction at the sample pixel location in the fraction image generated from the MESMA analysis. To further assess the accuracy of the classification result, a number of random sample points of the asphalt pavement in the study area was selected from the field measurements, regardless of the endmembers used in the unmixing analysis.


Results and Analysis


Spectral Analysis of the Aged Asphalt Pavement

The asphalt pavement shows different spectral patterns in different conditions of aging due to the change of ingredients and structure of the road pavement. Some representative spectral curves of the asphalt pavement for the first three conditions are illustrated in Fig. 6. The spectral curves labeled as P1, P2, and P3 were taken from the road pavement in the preliminarily aged condition. Specifically, P1 was from fresh asphalt, whereas P2 and P3 had been just lightly aged. The average reflectance of these three pavement samples was very low and the peak reflectance was <10% at the wavelength of 0.35 to 2.5  μm. Two hydrocarbon absorption zones were found around the wavelengths of 1.72 and 2.3  μm, which is consistent with the findings of Herold.5 The photos M1, M2, and M3 in Fig. 6 were taken of the pavements in a moderately aged condition, in which the lithologic content was exposed to the road surface. The overall spectral reflectance in this condition increased as expected, and the absorption of hydrocarbon observed for samples in the lightly aged condition could no longer be found. Instead, a strong absorption at the wavelength of 2.33  μm by the lithologic component of the pavement was found. Moreover, the spectral slope from visible- to near-infrared appears to be sharper. In the heavily aged condition (e.g., H1, H2, and H3 in Fig. 6), the milled rock pebbles and less asphalt in the pavement resulted in a visibly flatter pavement. Hence, the overall reflectance increased considerably (e.g., a peak of 24% observed for H3). The reflectance in the visible band continuously increases further through the aging process, and three weak absorption features related with iron oxide emerged at 0.52, 0.67, and 0.87  μm.16,35 After a turning point at 0.78  μm, the reflectance showed a steady decline with a strong absorption by mineral at 2.33  μm.

Fig. 6

The spectra and corresponding digital photos and color values of the asphalt pavement in preliminarily aging stage (P1, P2, and P3), moderately aging stage (I1, I2, and I3), and heavily aging stage (L1, L2, and L3).


In sum, it is clear that the asphalt pavement in different aging conditions presents different spectral patterns from 0.35 to 2.5  μm, especially in visible- and near-infrared bands. This makes it possible to map the asphalt pavement conditions based on Wordview-2 imagery.


Pavement Condition Mapping

A synthetic unmixing analysis based on the coupled strategy proposed in Sec. 2.4 was performed with a threshold value (σ) of 0.01. This means that if the difference of RMSE by EM2 and EM3 is >0.01, the EM3 model is used to unmix the pixel; otherwise, the EM2 model is used to retrieve the fraction of the endmembers. Table 2 shows the number and percentage of the pixels processed by the different models. As can be seen, 81.5% of the pixels were unmixed by the EM2 model, whereas 18.11% of the pixels were unmixed by the EM3 model.

The fraction maps of the pavement, shade, and nonpavement endmembers are shown in Fig. 7. It can be seen that the Liang Xiang East Area No. 14 Road has a much bigger fraction of PA endmember, suggesting this road pavement was mostly in the PA condition [Fig. 7(a)]. The East Ring Road and the Baiyang Road had a larger fraction of MA endmember [Fig. 7(b)], and the Yangguang South Street had a larger fraction of HA endmember [Fig. 7(c)]. The fraction of the other nonpavement land covers was mapped in Fig. 7(d), and that of shade was illustrated in Fig. 7(e).

Fig. 7

The fraction maps for the (a) preliminarily aged, (b) moderately aged, (c) heavily aged pavements, (d) nonpavement land covers (i.e., vegetation, sidewalks, traffic line, and soil), and (e) shades based on the subpixel unmixing approach.


For more effective visual interpretation, the value of a pixel is set to 1 if the fraction of a pavement endmember in it is >0.5. The three thresholding fraction images of the heavily, moderately, and preliminarily aged pavements were used to generate an RGB false composite image (Fig. 8). The image classification showed a fairly good agreement between the modeled road surface condition and the actual measurement (Fig. 8).

Fig. 8

The fraction image for the heavily aging pavement (HA), the preliminarily aging pavement (PA), and the moderately aging pavement (MA) in RGB and the photos taken in A, B, C, and D sites.


Specifically, it can be seen that the pavement health conditions in the sites of A, B, C, and D are visually consistent with what is shown in the photos taken at these locations. For example, the surface of the pavement that was evaluated as heavily aged (in red) shows to be in a bright gray color (meaning less asphalt) on the corresponding photo. The pavement of the Liangxiang East Area No. 14 Road was classified as being in the preliminary aging stage.

The photos taken of site A (Fig. 8) suggest that the road had been repaved recently and therefore was in relatively the best condition. The Yangguang South Street had been in service for about 9 years without any maintenance and was mostly recognized as heavily aged (bright red pixels in Fig. 8). The intersection of the East Ring Road and the Baiyang East Road was moderately aged. In site D, the different lanes of the Baiyang East Road appear to be in different aging stages.


Accuracy Assessment

According to the accuracy assessment procedure described in Sec. 2.5, the pan-sharpened Worldview-2 multispectral image was classified into five classes (PA pavement, MA pavement, HA pavement, shade, and nonpavement land covers) by visual interpretation. Afterward, 20 samples for each type of endmember were selected from the 2-m reference image and used to conduct linear regression analysis with the modeled fraction at the sample pixel location in the fraction image generated from the MESMA analysis. The results are presented in Fig. 9 and show that the MESMA-derived fraction was well correlated with the reference fraction (R2>0.7). For the classification result, 246 random sample points of the asphalt pavement in the study area were selected from the field measurements (Fig. 10). The relationship between the Munsell color values and the aging stages collected in field investigation was used to match the classes of the 165 pavement samples with their color values (Fig. 10). Statistical analysis showed that 81.71% of the points could be correctly matched with a Kappa coefficient of 0.77.

Fig. 9

Correlation analysis between the modeled and reference fraction: (a) preliminary aged pavement, (b) moderately aged pavement, (c) heavily aged pavement, (d) nonpavement objects, and (e) shade.


Fig. 10

The matching result between the color values and the classification results. (The size of the circles indicates the number of the matched or nonmatched points. The color values from 0.5 to 8.75 in x-axis correspond to the Munsell color values in Table 1. Y-axis presents five types of land cover classification results from the experiment above, i.e., PA: preliminarily aging pavement; MA, moderately aging pavement; HA, heavily aging pavement; NP, nonpavement objects; SD, shade.)



Statistical Analysis of Pavement Conditions

To quantitatively evaluate the aging condition of the asphalt road pavement in the study area, the area and percentage of the pavement in the three aging conditions were calculated and listed in Table 3. The results suggest that the Liangxiang East Area No. 14 Road was in the best condition, with 59.40% and 25.45% of its pavement being preliminarily aged and moderately aged, respectively (Fig. 11). In contrast, the Yangguang South Street and the Liangxiang East Area No. 16 Road show a significantly different pattern. The percentage of heavily aged pavement on these two roads reached up to 88.17% and 80.02%, respectively, presenting an urgent need for maintenance. For the East Ring Road, about 53.03% was classified as being moderately aged while 44.53% as heavily aged. It is obvious that the information derived from remotely sensed data is valuable for the management department to make informed plans for effective and timely road maintenance.

Table 3

The area of the road pavements in the three aging conditions.

Road namePreliminary aging pavement (m2)Moderately aging pavement (m2)Heavily aging pavement (m2)
Liangxiang East Area No. 14 Road7367.973156.671879.65
Baiyang East Road6378.8721,451.6711,479.34
East Ring Road567.8912,351.6710,371.36
Liangxiang East Area No. 16 Road165.391023.454975.82
Yangguang South Street1056.461289.2917,489.87

Fig. 11

The percentage of the area of the aged pavements in the roads (%).




The application of the proposed method in mapping road pavement condition in the study area has been shown preliminarily to be effective. The coupling strategy of two- and three-endmember models allows for more accurate mapping of asphalt pavement conditions, especially near the road edges. For example, the EM2 models built for the study area are shown in Fig. 12. Different colors are used to indicate different two-endmember models used in the unmixing analysis. Specifically, 82% of the image pixels (i.e., the ratio of the number of modeled pixels to the number of entire image pixels without background) were unmixed successfully using the EM2 model except those near the border of roads and greenbelts (e.g., black pixels in the areas A and B in Fig. 12) because of the lack of appropriate endmembers in the models. In contrast, the EM3 model was more effective for those mixed pixels that are located between roads and greenbelts except those covered by vehicles (e.g., black pixels in the site A in Fig. 13) and 98.5% of valid pixels (i.e., the ratio of the number of modeled pixels to the number of entire image pixels without a background) with a low average RMSE were modeled successfully. It was found that the EM2 model was more appropriate for analyzing those pixels within the central road area, which was much less impacted by surrounding objects. The EM3 model was more suitable for analyzing those pixels near the road edges, where more nonpavement objects such as soil and vegetation tended to be captured in a pixel as well.

Fig. 12

The EM2 model image in the study area. (Color pixels indicate different EM2 models, and black pixels indicate unsuccessful unmixing.)


Fig. 13

The EM3 model image in the study area. (Color pixels indicate different EM3 models, and black pixels indicate unsuccessful unmixing.)


The analysis result is clearly valuable for the highway maintenance department to make timely and suitable maintenance plans. However, our approach also has some weakness. First of all, the shade endmember is simply treated as a low reflective object, despite that shading is more complicated in reality. The proposed method can be further improved by using multiseasonal images to remove or reduce the shades by the objects near roads prior to the subpixel unmixing analysis. Second, based on the in situ measured pavement spectral data, it can be seen in Fig. 3 that more spectral signatures can be observed in the wavelengths >1.0  μm, especially in those >2.0  μm. Unfortunately, the Worldview satellite does not acquire data to cover this spectral range. On the other hand, the distress condition of asphalt road pavement cannot be accurately detected well with our proposed approach due to the fact that its spectral signature is not very distinct. LiDAR might be a good data source for providing the 3-D geometric information of the distress pavement. In addition, bridging the generated pavement condition maps with the operational regulations of the maintenance department is yet to be studied for more practical applications.



This paper presents an experimental study of mapping road pavement conditions from the WorldView-2 multispectral imagery based on an MESMA. The analysis of the in situ measured spectra of asphalt road pavement showed that the asphalt pavements in different aging conditions present quite different spectral patterns over the range of 0.35 to 2.5  μm. The aging conditions of the road pavements were detected and assessed through the spectral unmixing analysis of the Worldview-2 image in our study area located in the southwest of Beijing. The assessment of the analysis results showed that the multiple endmember mixture analysis is capable of classifying the road pavements that are preliminarily aged, moderately aged, and heavily aged, with an overall accuracy of 81.71% and the Kappa coefficient of 0.77, when the in situ data were used as reference. It was found that the two-endmember model seems to be more appropriate for analyzing the pixels representing the central road areas over the three-endmember model, due to the fact that fewer endmembers exist in these areas. The three-endmember model, however, performs better for analyzing the road edge areas, where the surface is often covered by a mixture of natural and artificial objects. The proposed approach integrating the two- and three-endmember models has shown to be effective and has the great potential in supporting the planning and decision making of road pavement maintenance.

Important spectral characteristics such as the dramatic change of absorption intensity at the wavelength of 2.3  μm might be used to develop an alternative approach for detecting and mapping pavement condition. The use of hyperspectral imagery with higher spectral and spatial resolution and wider spectral coverage is expected to produce better results than what is from multispectral imagery. Moreover, a more accurate detection and mapping of pavement condition may be achieved by synthetically analyzing the spectral information from a remote sensor and the geometric information from a LiDAR sensor in future work.


This study was financially supported by two grants from the National Natural Science Foundation of China (No. 41571331) and the China Transport Telecommunications and Information Center (No. GFZX0404080102). The authors would like to give our great thanks to the anonymous reviewers for their valuable comments and suggestions.


1. M. C. Rubio et al., “Warm mix asphalt: an overview,” J. Cleaner Prod. 24, 76–84 (2012). http://dx.doi.org/10.1016/j.jclepro.2011.11.053 Google Scholar

2. S. S. Pulugurtha et al., “Preliminary results towards developing thresholds for pavement condition maintenance: safety perspective,” Proc. Soc. Behav. Sci. 104, 302–311 (2013). http://dx.doi.org/10.1016/j.sbspro.2013.11.123 Google Scholar

3. M. Herold et al., “Imaging spectrometry and asphalt road surveys,” Transp. Res. Part C: Emerging Technol. 16(2), 153–166 (2008). http://dx.doi.org/10.1016/j.trc.2007.07.001 Google Scholar

4. D. M. L. Marcelo et al., “Characterization of aging processes on the asphalt mixture surface,” in 2nd Int. Symp. on, Asphalt Pavements et Environment, Transportation Research Board, of The National Academies, France, p. 10 (2012). Google Scholar

5. M. Herold and D. Roberts, “Spectral characteristics of asphalt road aging and deterioration: implications for remote-sensing applications,” Appl. Opt. 44(20), 4327–4334 (2005).APOPAI0003-6935 http://dx.doi.org/10.1364/AO.44.004327 Google Scholar

6. N. Bandara and M. Gunaratne, “Current and future pavement maintenance prioritization based on rapid visual condition evaluation,” J. Transp. Eng. 127(2), 116–123 (2001).JTPEDI0733-947X http://dx.doi.org/10.1061/(ASCE)0733-947X(2001)127:2(116) Google Scholar

7. E. Schnebele et al., “Review of remote sensing methodologies for pavement management and assessment,” Eur. Transp. Res. Rev. 7(2), 1–19 (2015). http://dx.doi.org/10.1007/s12544-015-0156-6 Google Scholar

8. Y. U. Shah et al., “Development of overall pavement condition index for urban road network,” Proc. Soc. Behav. Sci. 104, 332–341 (2013). http://dx.doi.org/10.1016/j.sbspro.2013.11.126 Google Scholar

9. K. Themistocleous et al., “Damage assessment using advanced non-intrusive inspection methods: integration of space, UAV, GPR, and field spectroscopy,” Proc. SPIE 9229, 92291O (2014).PSISDG0277-786X http://dx.doi.org/10.1117/12.2069507 Google Scholar

10. R. Haas, W. R. Hudson and J. P. Zaniewski, Modern Pavement Management, Krieger Publishing Company, Malabar, Florida (1994). Google Scholar

11. M. Ahmed, C. Haas and R. Haas, “Toward low-cost 3D automatic pavement distress surveying: the close range photogrammetry approach,” Can. J. Civ. Eng. 38(12), 1301–1313 (2011).CJCEB81208-6029 http://dx.doi.org/10.1139/l11-088 Google Scholar

12. S. Chambon and J.-M. Moliard, “Automatic road pavement assessment with image processing: review and comparison,” Int. J. Geophys. 2011, 1–20 (2011). http://dx.doi.org/10.1155/2011/989354 Google Scholar

13. R. B. Gomez, “Hyperspectral imaging: a useful technology for transportation analysis,” Opt. Eng. 41(9), 2137–2143 (2002). http://dx.doi.org/10.1117/1.1497985 Google Scholar

14. M. R. Resende, L. L. B. Bernucci and J. A. Quintanilha, “Monitoring the condition of roads pavement surfaces: proposal of methodology using hyperspectral images,” J. Transp. Lit. 8(2), 201–220 (2014). http://dx.doi.org/10.1590/S2238-10312014000200009 Google Scholar

15. E. P. Paska, State-of-the-Art Remote Sensing Geospatial Technologies in Support of Transportation Monitoring and Management, The Ohio State University, Columbus, Ohio (2009). Google Scholar

16. A. Mei et al., “Spectroradiometric laboratory measures on asphalt concrete: preliminary results,” Proc. Soc. Behav. Sci. 53, 514–523 (2012). http://dx.doi.org/10.1016/j.sbspro.2012.09.902 Google Scholar

17. J. Marcondes et al., “Spectral analysis of highway pavement roughness,” J. Transp. Eng. 117(5), 540–549 (1991).JTPEDI0733-947X http://dx.doi.org/10.1061/(ASCE)0733-947X(1991)117:5(540) Google Scholar

18. A. Mei et al., “Integration of field and laboratory spectral data with multi-resolution remote sensed imagery for asphalt surface differentiation,” Remote Sens. 6(4), 2765–2781 (2014).RSEND3 http://dx.doi.org/10.3390/rs6042765 Google Scholar

19. C. Manzo et al., “Spectral modelling used to identify the aggregates index of asphalted surfaces and sensitivity analysis,” Constr. Build. Mater. 61, 147–155 (2014).CBUMEZ http://dx.doi.org/10.1016/j.conbuildmat.2014.02.056 Google Scholar

20. S. Pascucci et al., “Road asphalt pavements analyzed by airborne thermal remote sensing: preliminary results of the Venice Highway,” Sensors 8(2), 1278–1296 (2008).SNSRES0746-9462 http://dx.doi.org/10.3390/s8021278 Google Scholar

21. A. Moropoulou et al., “An application of thermography for detection of delaminations in airport pavements,” NDT&E Int. 34(5), 329–335 (2001). http://dx.doi.org/10.1016/S0963-8695(00)00047-5 Google Scholar

22. C. Small, “High spatial resolution spectral mixture analysis of urban reflectance,” Remote Sens. Environ. 88(1), 170–186 (2003).RSEEA70034-4257 http://dx.doi.org/10.1016/j.rse.2003.04.008 Google Scholar

23. D. A. Roberts et al., “Mapping chaparral in the Santa Monica Mountains using multiple endmember spectral mixture models,” Remote Sens. Environ. 65(3), 267–279 (1998).RSEEA70034-4257 http://dx.doi.org/10.1016/S0034-4257(98)00037-6 Google Scholar

24. B. Ayalew et al., “Pavement management using hyperspectral imagery,” Proc. SPIE 5097, 207 (2003).PSISDG0277-786X http://dx.doi.org/10.1117/12.502588 Google Scholar

25. A. Mei et al., “Bitumen removal determination on asphalt pavement using digital imaging processing and spectral analysis,” Open J. Appl. Sci. 4(6), 46087 (2014). http://dx.doi.org/10.4236/ojapps.2014.46034 Google Scholar

26. NREC, “NERC field spectroscopy facility,” (2016),  http://fsf.nerc.ac.uk/, (12 June 2016). Google Scholar

27. J. Franke et al., “Hierarchical multiple endmember spectral mixture analysis (MESMA) of hyperspectral imagery for urban environments,” Remote Sens. Environ. 113(8), 1712–1723 (2009).RSEEA70034-4257 http://dx.doi.org/10.1016/j.rse.2009.03.018 Google Scholar

28. J. B. Adams, M. O. Smith, A. R. Gillespie, “Imaging spectroscopy: interpretation based on spectral mixture analysis,” in Remote Geochemical Analysis: Elemental and Mineralogical Composition, , C. M. Pieters and P. A. J. Englert, Eds., pp. 145–166, Cambridge University Press, Cambridge, United Kingdom (1993). Google Scholar

29. F. Fan and Y. Deng, “Enhancing endmember selection in multiple endmember spectral mixture analysis (MESMA) for urban impervious surface area mapping using spectral angle and spectral distance parameters,” Int. J. Appl. Earth Obs. Geoinf. 33, 290–301 (2014). http://dx.doi.org/10.1016/j.jag.2014.06.011 Google Scholar

30. T. Rashed et al., “Measuring the physical composition of urban morphology using multiple endmember spectral mixture models,” Photogramm. Eng. Remote Sens. 69(9), 1011–1020 (2003). http://dx.doi.org/10.14358/PERS.69.9.1011 Google Scholar

31. R. L. Powell et al., “Sub-pixel mapping of urban land cover using multiple endmember spectral mixture analysis: Manaus, Brazil,” Remote Sens. Environ. 106(2), 253–267 (2007).RSEEA70034-4257 http://dx.doi.org/10.1016/j.rse.2006.09.005 Google Scholar

32. P. E. Dennison and D. A. Roberts, “Endmember selection for multiple endmember spectral mixture analysis using endmember average RMSE,” Remote Sens. Environ. 87(2), 123–135 (2003).RSEEA70034-4257 http://dx.doi.org/10.1016/S0034-4257(03)00135-4 Google Scholar

33. D. A. Roberts et al., “Synergies between VSWIR and TIR data for the urban environment: an evaluation of the potential for the hyperspectral infrared imager (HyspIRI) decadal survey mission,” Remote Sens. Environ. 117, 83–101 (2012).RSEEA70034-4257 http://dx.doi.org/10.1016/j.rse.2011.07.021 Google Scholar

34. D. Roberts, K. Halligan and P. Dennison, VIPER Tools User Manual, UC Santa Barbara, Department of Geography, Visualization and Image Processing for Environmental Research Laboratory, Santa Barbara, California (2007). Google Scholar

35. M. Herold et al., “Spectrometry for urban area remote sensing—development and analysis of a spectral library from 350 to 2400 nm,” Remote Sens. Environ. 91(3), 304–319 (2004).RSEEA70034-4257 http://dx.doi.org/10.1016/j.rse.2004.02.013 Google Scholar


Yifan Pan received his BE degree in remote sensing from Shandong University of Science and Technology, Qingdao, China, in 2013. He is currently a PhD student in the Institute of Remote Sensing and GIS, Peking University, Beijing, China. His research focuses on the application of hyperspectral remote sensing.

Xianfeng Zhang received his BS degree from Beijing Normal University, his MSc degree from Peking University, Beijing, China, in 1990 and 1993, respectively, and his PhD from the University of Western Ontario, London, Canada, in 2005. He is currently an associate professor with the Institute of Remote Sensing and GIS, Peking University. His research interests include methods for information extraction from multisensor imagery and its applications in urban, wetland, and agricultural cropland environments.

Jie Tian is currently an assistant professor with Clark University, USA. He received his MS degree from University of Western Ontario and his PhD from Queen's University, Canada in 2004 and 2009, respectively.

Xu Jin is currently a master student at Peking University.

Lun Luo is currently a senior researcher.

Ke Yang is an associate researcher at China Transport Telecommunications and Information Center, Ministry of Transport of China.

© The Authors. Published by SPIE under a Creative Commons Attribution 3.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.
Yifan Pan, Xianfeng Zhang, Jie Tian, Xu Jin, Lun Luo, Ke Yang, "Mapping asphalt pavement aging and condition using multiple endmember spectral mixture analysis in Beijing, China," Journal of Applied Remote Sensing 11(1), 016003 (6 January 2017). https://doi.org/10.1117/1.JRS.11.016003 Submission: Received 12 July 2016; Accepted 16 December 2016
Submission: Received 12 July 2016; Accepted 16 December 2016




Remote sensing

RGB color model

Statistical analysis


Back to Top