Mapping asphalt pavement aging and condition using multiple endmember spectral mixture analysis in Beijing, China

Abstract. 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.


Introduction
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) [9][10][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,[13][14][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,[18][19][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, Roberts 23 ), 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.

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.

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.
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 nearinfrared 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: E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 1 ; 1 1 6 ; 4 4 7 where ρ 0 λ is the reflectance of the image acquired by remote sensor in band λ, and ρ iλ is the reflectance of the ith endmember; f i 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: E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 2 ; 1 1 6 ; 3 6 2 The accuracy of the model can be evaluated using root mean square error (RMSE) E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 3 ; 1 1 6 ; 3 0 8 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 endmembers 28 ) 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. 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 E Q -T A R G E T ; t e m p : i n t r a l i n k -; e 0 0 4 ; 1 1 6 ; 4 4 8 Fig. 4 The flowchart of the MESMA with the Worldview-2 image.
where A is the modeled spectra class, n is the number of spectra in class A, and A i , A j is the endmember in class A. The term n − 1 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 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.

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 pansharpened 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.

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.
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 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.
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).
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).
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 (R 2 > 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.

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  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.) 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.

Discussion
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. 11 The percentage of the area of the aged pavements in the roads (%).  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. 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.

Conclusion
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.