Improving Fractional Impervious Surface Mapping Performance through Combination of DMSP-OLS and MODIS NDVI Data

Impervious surface area (ISA) is an important parameter for many studies such as urban climate, urban environmental change, and air pollution; however, mapping ISA at the regional or global scale is still challenging due to the complexity of impervious surface features. The Defense Meteorological Satellite Program’s Operational Linescan System (DMSP-OLS) data have been used for ISA mapping, but high uncertainty existed due to mixed-pixel and data-saturation problems. This paper presents a new index called normalized impervious surface index (NISI), which is an integration of DMSP-OLS and Moderate Resolution Imaging Spectroradiometer (MODIS) normalized difference vegetation index (NDVI) data, in order to reduce these problems. Meanwhile, this newly developed index is compared with previously used indices—Human Settlement Index (HSI) and Vegetation Adjusted Nighttime light Urban Index (VANUI)—in ISA mapping performance. We selected China as an example to map fractional ISA distribution through a support vector regression approach based on the relationship between the index and Landsat-derived ISA data. The results indicate that the proposed NISI provided better ISA estimation accuracy than HSI and VANUI, especially when the fractional ISA in a pixel is relatively large (i.e., >0.6) or very small (i.e., <0.2). This approach can be used to rapidly update ISA datasets at regional and global scales.


Introduction
Rapid population increase and economic conditions have resulted in unprecedented urbanization in China in the past four decades [1][2][3].For example, China's population increased from 1.0 billion in 1981 to 1.27 billion in 2000 and to 1.37 billion in 2015 [4], while the percent of urban population increased from 20.9% in 1982 to 36.2% in 2000 and to 56.1% in 2015 [5].The rapid urbanization has produced important impacts on urban climate, air pollution, energy consumption, and ecosystem services, requiring mapping and monitoring of urban expansion in timely ways [6][7][8][9][10][11]. Remote sensing, due to its unique characteristics in data collection and presentation, has been extensively applied for mapping urban distribution and dynamic change [3,8,[12][13][14][15][16][17].However, the complex urban landscapes and confusion of different land cover types make it difficult to directly use remote sensing data for mapping urban land cover distribution and detecting its change.Impervious surface area (ISA) has been regarded as an important attribute to represent urban distribution and expansion [18].Therefore, much research has been conducted to accurately map ISA distribution using different remote sensing data, especially Landsat imagery, due to its long-term data availability at no cost and suitable spectral and spatial resolutions [19][20][21][22].
At national and global scales, the Defense Meteorological Satellite Program's Operational Linescan System (DMSP-OLS) data are often used to map ISA distribution [23][24][25].The threshold-based approach has been used to extract ISA [26][27][28], but the results may have high uncertainties, depending on the selection of thresholds, and no thresholds can effectively represent the ISA amount due to the impacts of different economic conditions on nighttime light data in a large area [29].The DMSP-OLS data originally have spatial resolution of approximately 2.7 km, but the product was resampled to a 1 km cell size with data range of 0-63, thus the data have serious mixed-pixel and data-saturation problems [30][31][32][33].The blooming effect of DMSP-OLS data especially overestimates ISA in urban-suburban regions.Therefore, DMSP-OLS data alone cannot provide accurate ISA mapping due to the limitation of the data themselves.
Moderate Resolution Imaging Spectroradiometer (MODIS) is another source of coarse spatial resolution data for ISA mapping at regional and global scales [13,[34][35][36][37][38].Schneider et al. [36] used time-series MODIS data to map urban extent at the global scale using supervised decision tree classification.Similar research was also conducted at the continental scale using MODIS data to detect urban area change [13].One critical step in using MODIS data is to distinguish ISA from other non-vegetation covers such as bare lands and water, in addition to reducing the mixed-pixel problem causing by coarse spatial resolution.Therefore, ISA mapping using MODIS data alone is difficult because of the spectral confusion between ISA and other land covers such as bare soils and water and mixed-pixel problems in urban landscapes.
DMSP-OLS and MODIS data have different characteristics in representing land covers.DMSP-OLS represents the nighttime light that is closely related to ISA distribution, but the coarse spatial resolution (1 km), data saturation (data range is 0-63), and blooming effect make it inaccurate in ISA mapping.In contrast, the MODIS normalized difference vegetation index (NDVI) has relatively higher spatial resolution than DMSP-OLS (i.e., 250 m for MODIS NDVI vs. 1 km for the DMSP-OLS product), and has better representation for different land covers such as vegetation and non-vegetation types (e.g., ISA, water, bare soils).The 16-bit data (0-65,535) in MODIS data significantly reduce the data-saturation problem, thus providing much detailed information for separating land covers.These different characteristics of DMSP-OLS and MODIS datasets provide the opportunity to make full use of both advantages to obtain abundant urban land cover information.
Previous studies have proven that the integration of DMSP-OLS and MODIS NDVI can make effective use of both advantages to reduce the blooming effects, data-saturation, and mixed-pixel problems in DMSP-OLS, thus considerably improving the ISA mapping performance at regional and global scales [29,[39][40][41].Since Lu et al. [29] proposed a Human Settlement Index (HSI) based on the combination of DMSP-OLS and MODIS NDVI data to estimate settlement areas in south and east China, much research has been conducted to combine DMSP-OLS data and a vegetation index (e.g., MODIS NDVI or EVI, SPOT VEGETATION) for improving ISA mapping performance [40][41][42][43].The Vegetation Adjusted nighttime light Urban Index (VANUI) based on the normalized DMSP-OLS and the mean value of the time-series MODIS NDVI data has been regarded as a better index than HSI in mapping ISA distribution [40].Based on a similar idea, Zhuo et al. [41] proposed the enhanced vegetation index (EVI)-adjusted nighttime light index (EANTLI)-to reduce the data saturation problem in DMSP-OLS data.A recent study using the integration of the Suomi National Polar-orbiting Partnership Visible Infrared Imaging Radiometer Suite's Day/Night Band (VIIRS DNB) and MODIS NDVI further improved ISA estimation by using higher spatial resolution data (VIIRS DNB with 750 m vs. DMSP-OLS with 1 km spatial resolution) and reducing the data-saturation problems [44].
In ISA mapping research, linear regression analysis is commonly used to establish the relationship between ISA and a DMSP-OLS-based or NDVI-based index [29,44], but very little research had explored the nonparametric approaches to map urban distribution [42].Although previous research has made good progress in integrating DMSP-OLS and MODIS data to improve ISA mapping performance, there are some shortcomings in previously used indices.For example, the HSI [29] cannot be effectively used to estimate ISA in the low vegetation cover regions such as the western part of China due to dry weather; the VANUI [33] resulted in inaccurate ISA estimation for some regions such as the moist tropical regions because the mean value of MODIS NDVI time series data cannot remove the impacts of clouds/shadows and their contamination in the 8-day or 16-day composite.Previous research indicated that the ISA mapping accuracy in low or high ISA density regions are still poor [29,40,44], and a new approach is required to improve ISA estimation.Therefore, this research aims at developing a new index called Normalized Impervious Surface Index (NISI) by integrating DMSP-OLS and MODIS NDVI data to improve the ISA estimation accuracy, and to explore a nonparametric approach-a support vector regression (SVR)-to build the relationship between the proposed index and reference ISA data for estimating fractional ISA.

Study Area and Datasets
China has experienced unprecedented urbanization rates since the early 1980s due to the rural-to-urban population migration and economic conditions [2,3,45,46].The urbanization patterns and rates vary greatly from coastal regions to central and western parts of the country and from southeastern to northeastern parts [45][46][47].In order to develop ISA reference data for use in ISA modeling, ten metropolitan areas in different geographic locations and climate zones were selected (Figure 1).They have various population sizes, economic conditions, and administrative areas as summarized in Table 1.The complex landscapes, economic conditions, and topographic conditions in China provide an ideal example to explore approaches to rapidly map ISA spatial distribution and dynamic change.
In ISA mapping research, linear regression analysis is commonly used to establish the relationship between ISA and a DMSP-OLS-based or NDVI-based index [29,44], but very little research had explored the nonparametric approaches to map urban distribution [42].Although previous research has made good progress in integrating DMSP-OLS and MODIS data to improve ISA mapping performance, there are some shortcomings in previously used indices.For example, the HSI [29] cannot be effectively used to estimate ISA in the low vegetation cover regions such as the western part of China due to dry weather; the VANUI [33] resulted in inaccurate ISA estimation for some regions such as the moist tropical regions because the mean value of MODIS NDVI time series data cannot remove the impacts of clouds/shadows and their contamination in the 8-day or 16-day composite.Previous research indicated that the ISA mapping accuracy in low or high ISA density regions are still poor [29,40,44], and a new approach is required to improve ISA estimation.Therefore, this research aims at developing a new index called Normalized Impervious Surface Index (NISI) by integrating DMSP-OLS and MODIS NDVI data to improve the ISA estimation accuracy, and to explore a nonparametric approach-a support vector regression (SVR)-to build the relationship between the proposed index and reference ISA data for estimating fractional ISA.

Study Area and Datasets
China has experienced unprecedented urbanization rates since the early 1980s due to the rural-to-urban population migration and economic conditions [2,3,45,46].The urbanization patterns and rates vary greatly from coastal regions to central and western parts of the country and from southeastern to northeastern parts [45][46][47].In order to develop ISA reference data for use in ISA modeling, ten metropolitan areas in different geographic locations and climate zones were selected (Figure 1).They have various population sizes, economic conditions, and administrative areas as summarized in Table 1.The complex landscapes, economic conditions, and topographic conditions in China provide an ideal example to explore approaches to rapidly map ISA spatial distribution and dynamic change.Different data sets as summarized in Table 2 were used in this research.The DMSP-OLS nighttime light data products (Version 4) with 6-bit data range and cell size of 1 × 1 km have the capability to acquire city lights, even low intensity from small residences and traffic lights, while sunlight, moonlight, glare, clouds, fires, and lighting features from auroras were removed through one-year data [30,48].The MODIS NDVI 16-day composite (MOD13Q1) was generated from the daily MODIS Level-2G (L2G) surface reflectance.The constrained view angle-maximum value composite (CV-MVC) method was used to minimize the cloud contamination.The maximum NDVI value with the smallest view angle or the closest-to-nadir view for each pixel per 16-day composite interval was generated using MODIS Vegetation Indices (MOD13).The water mask product (MOD44W) with the cell size of 250 m × 250 m utilizes the SWBD (SRTM Water Body Data) and complements it with information from 250 m MODIS data to create a complete representation of global surface water (see detailed description in [49]).Landsat 8 Operational Land Imager (OLI) reflective bands were used to map ISA distribution for the selected sites.

Development of Impervious Surface Area Data from Landsat OLI Images
The Landsat 8 OLI data were re-projected from the Universal Transverse Mercator coordinate system to Albers Conical Equal Area projection, and the nearest-neighbor algorithm was used to resample the Landsat 8 OLI imagery with 30 m spatial resolution into a pixel size of 25 m × 25 m.A combination of thresholding and cluster analysis was used to extract ISA data from Landsat imagery [44].We first produced NDVI and normalized difference water index (NDWI) to extract vegetation and water bodies.Based on analysis of vegetation and water sample plots, we selected a threshold of 0.3 in NDVI data to extract vegetation pixels (i.e., if a pixel has NDVI value of greater than 0.3, the pixel was assigned as vegetation) and a threshold of 0.1 in NDWI to extract water pixels (i.e., if a pixel has NDWI value of greater than 0.1, the pixel was assigned as water).After masking vegetation and water, the remaining pixels were mainly ISA and bare soils.Their spectral signatures were then extracted from the OLI imagery, and cluster analysis was conducted to classify the remaining pixels into 100 clusters.The analyst merged the clusters into ISA and others based on analysis of their spectral signatures.The ISA imagery was visually examined by overlaying it on the Landsat 8 OLI color composite to make sure the ISA pixels were properly extracted with high accuracy.Finally, the ISA imagery was aggregated from a cell size of 25 m × 25 m to a cell size of 250 m × 250 m to generate fractional ISA imagery based on the calculation of mean value within the window size of 10 by 10 pixels.

Development of Impervious Surface Area Data from Landsat OLI Images
The Landsat 8 OLI data were re-projected from the Universal Transverse Mercator coordinate system to Albers Conical Equal Area projection, and the nearest-neighbor algorithm was used to resample the Landsat 8 OLI imagery with 30 m spatial resolution into a pixel size of 25 m × 25 m.A combination of thresholding and cluster analysis was used to extract ISA data from Landsat imagery [44].We first produced NDVI and normalized difference water index (NDWI) to extract vegetation and water bodies.Based on analysis of vegetation and water sample plots, we selected a threshold of 0.3 in NDVI data to extract vegetation pixels (i.e., if a pixel has NDVI value of greater than 0.3, the pixel was assigned as vegetation) and a threshold of 0.1 in NDWI to extract water pixels (i.e., if a pixel has NDWI value of greater than 0.1, the pixel was assigned as water).After masking vegetation and water, the remaining pixels were mainly ISA and bare soils.Their spectral signatures were then extracted from the OLI imagery, and cluster analysis was conducted to classify the remaining pixels into 100 clusters.The analyst merged the clusters into ISA and others based on analysis of their spectral signatures.The ISA imagery was visually examined by overlaying it on the Landsat 8 OLI color composite to make sure the ISA pixels were properly extracted with high accuracy.Finally, the ISA imagery was aggregated from a cell size of 25 m × 25 m to a cell size of 250 m × 250 m to generate fractional ISA imagery based on the calculation of mean value within the window size of 10 by 10 pixels.

Preprocessing and Integration of DMSP-OLS and MODIS NDVI Data
Both DMSP-OLS and MODIS NDVI datasets were re-projected into Albers Conical Equal Area coordinate system.The nearest-neighbor resampling algorithm was used to resample the DMSP-OLS data to a cell size of 1 km × 1 km and resample the time-series MODIS NDVI data to a cell size of 250 m × 250 m.The DMSP-OLS data were normalized using Equation (1): where OLS nor is the normalized DMSP-OLS image.OLS min and OLS max represent the minimum and maximum values in the DMSP-OLS imagery, here they are 0 and 63.The cloud and shadow contamination problem is common in the time-series MODIS NDVI data, especially in subtropical and tropical regions.In order to eliminate this problem and to reduce the confusion of bare soils in agricultural lands and ISA, a maximum algorithm using Equation ( 2) was used to produce a new composite: where NDVI 1 , NDVI 2 , . . ., NDVI n are the MOD13Q1 NDVI time series images.The global water data (MOD44W) were used to mask out water bodies from both OLS nor and NDVI max datasets, thus they have a data range between 0 and 1.
In DMSP-OLS data, coarse spatial resolution and small data range (0-63) result in data saturation and blooming effects [29,48], while the relatively high spatial resolution and high radiometric resolution in MODIS data make it possible to better distinguish different land cover types [29,36,37].Because of the sensitivity of red and near-infrared wavelengths in separation of vegetation and non-vegetation types (e.g., ISA and water), NDVI has almost no data saturation problems in urban landscapes [29].The different characteristics of DMSP-OLS and MODIS NDVI data result in the integrated use of both data sets, and this integration has been proven to produce better ISA extraction performance than using individual data alone [29,[39][40][41].It is critical to take advantage of both data sources in reducing data saturation problems in urban landscapes and reducing the blooming effects in suburban regions.
The common indices used in previous research are HSI [29,43] and VANUI [40], but they cannot effectively extract ISA when the proportion of ISA in a pixel is very small (e.g., less than 0.2) or very high (e.g., greater than 0.8).Here, we propose a new index called Normalized Impervious Surface Index (NISI) (Equation ( 3)) based on OLSnor (1 km × 1 km) and NDVImax (250 m × 250 m) for improving ISA mapping performance: Here, we used NDVImax instead of NDVImean to eliminate cloud and shadow problems in tropical and subtropical regions and reduce the impacts of bare soils on the extraction of ISA.The inclusion of 2(1 − 0.5OLS nor )(1 + NDVI max ) is to enhance the separation of vegetation and non-vegetation types in the urban landscapes when vegetation covers account for a small portion in a unit.The use of 0.5 before OLSnor is to reduce the influence of OLSnor values, and the use of 2 is to make sure the NISI values fall within the range of 0 and 1.
(1−NDVI max ) (1+NDVI max ) is to enhance the difference between vegetation and non-vegetation, and OLS nor (1−0.5OLSnor ) is to enhance the difference of nighttime light values within the urban landscape.Therefore, NISI is to enlarge the difference between ISA and other land covers for better separation of ISA from other land covers in the complex urban landscape.The proposed index was compared to OLSnor and NDVImax using a typical line going through the urban landscape to examine the effectiveness of the new index in reducing data saturation and the blooming effects.In addition to the proposed NISI, another two indices-HSI and VANUI-were also calculated [29,40].

Mapping Fractional Impervious Surface Area Distribution with Support Vector Regression
The support vector machine is a supervised nonparametric statistical learning technique and has attracted great attention in recent decades due to its good performance in remote sensing-based classification [14,42,50,51].Support vector machine is an integrated software package for supporting category classification and regression estimation.In this research, we used the support vector regression (SVR) to map fractional ISA distribution.The Gaussian radial basis function (RBF) kernel was used in the SVR modeling.There are three major steps: (1) input the reference samples from Landsat OLI and corresponding samples from the NISI imagery; (2) use the RBF kernel of the SVR modeling and optimize parameters through a cross-validation method; and (3) input NISI imagery and map fractional ISA distribution based on the optimized parameters.The same procedure is also used for ISA prediction using HSI and VANUI imagery.
A random sampling technique was used to collect samples from the fractional ISA imagery, NISI, HSI and VANUI images for each selected city.The fractional ISA data from the Landsat 8 OLI images were used as a dependent variable and the individual index (e.g., NISI, HSI and VANUI) was used separately as an independent variable to establish the relationship between fractional ISA reference data and NISI, HSI and VANUI, respectively, using the SVR approach.A total of 1500 ISA samples were extracted from ten reference images with 150 samples for each, as well as from the NISI, HSI and VANUI images.Of these samples, 1000 samples were randomly selected for modeling, and the remaining 500 samples were used for evaluation of ISA estimates.The optimized SVR model was used to predict ISA distribution for the entire study area.

Evaluation of the Impervious Surface Area Estimation Results
Accuracy assessment is often conducted at pixel level using the error matrix approach [52], but this approach is not suitable for evaluation of fractional ISA estimates in this research.Therefore, correlation coefficient (R), root mean squared error (RMSE), relative RMSE (RMSEr), and residual analysis were used to assess the fitness between the ISA estimates and reference data [29,44].Meanwhile, the accuracy assessment was conducted at different ISA groups (e.g., <0.2, [0.2-0.4),[0.4-0.6), [0.6-0.8), and ≥0.8) to understand potential under/overestimation.In addition to the traditional approaches using R and RMSE for accuracy assessment, we also examined the accuracy of spatial prediction of fractional ISA maps by calculating means ( μ) using Equation (4) [53], where ŷi and y i represent the predicted ISA and reference data, and N and n represent the number of total predicted samples and reference samples:

A Comparative Analysis of the Proposed Index and Individual OLSnor and NDVImax Data
The effectiveness of the NISI data comparing with individual OLSnor and NDVImax data alone is illustrated in Figure 3.It indicates that the data-saturation problem in urban regions and blooming effects in the urban-rural interfaces are obvious in OLSnor data (Figure 3a,b), but not for NDVImax and NISI (Figure 3c,d).The NDVImax values vary, depending on the composition of vegetation and ISA in the urban landscape.The high spatial resolution in NDVImax data reduced the blooming effects (Figure 3c), and NISI provided much detailed information compared to OLSnor data alone, especially in the urban regions (Figure 3d).This implies that NISI has much more abundant information about urban land cover structures than OLSnor, thus the spatial patterns of ISA distribution may be much improved when NISI is used for ISA estimation.

Analysis of Impervious Surface Area Estimates
The accuracy assessment results in Table 3 indicate that the newly proposed NISI has slightly better performance than previously used two indices (i.e., HSI and VANUI).Overall, the NISI-based ISA prediction has the highest R and lowest RMSE and RMSEr values, followed by the VANUI-based predictions.HSI has relatively poor performance.The spatial patterns of ISA distributions using these indices are illustrated in Figure 4.It indicates that the VANUI-derived ISA distribution (Figure 4b) has obviously the highest number of ISA pixels, followed by the HSI-(Figure 4a) and NISI-derived ISA results (Figure 4c)).Considering one city, for example, the VANUI-derived ISA distribution in Beijing (Figure 4(b1)) has the highest number of yellow pixels (ISA in 0.5-0.7),and NISI-derived results have the highest number of red pixels (Figure 4(c1), ISA > 0.8), and HSI-derived results have the lowest red and yellow pixels (Figure 4(a1), ISA > 0.5).In contrast, the VANUI-derived results have the highest number of blue pixels (i.e., ISA < 0.2), followed by the HSI-derived results, while the NISI-derived results have the least.This may imply that the NISI can better reduce the data saturation problem in the urban landscape, especially when ISA accounts for high proportions in a unit, and VANUI may have higher chances to overestimate ISA distribution than NISI when ISA is very small.
The scatterplots between estimates and reference data and residual analysis results (see Figure 5) can be better used to explain the prediction performance.It is obvious that when ISA is high (i.e., greater than 0.8), the NISI (Figure 5(a3)) has better performance in ISA prediction than HSI (Figure 5(a1)) and VANUI (Figure 5(a2)).The residual distributions indicate that when ISA is small (e.g., less than 0.4), NISI and VANUI have similar performance (Figure 5(b2,b3)) and have better performance than HSI (Figure 5(b1)).For the three indices, the majority of the test samples have residual values within the range of −0.3 and 0.3.However, Figure 5(b1,b2,b3) also shows clearly the underestimation and overestimation problems.For example, when ISA is greater than 0.8, ISA pixels are commonly underestimated; and when ISA is less than 0.1, the ISA pixels are overestimated.
The under-or overestimation problems vary, depending on various cities, as shown in Figure 6.Overall, the ISA distribution using three indices for different cities present clear spatial patterns from the highest values in core urban areas to the lowest values in rural areas.However, the ISA estimation performances using the three indices are different; for example, NISI predicted the highest number of pixels with red (ISA > 0.9), VANUI predicted the highest number of pixels with yellow (ISA in 0.5-0.7)and blue (ISA < 0.1), while HSI predicted the highest number of pixels with cyan (ISA in 0.2-0.4) in these cities.All three indices overestimated ISA in Urumqi, and underestimated in Chengdu and Beijing, in comparison with the Landsat-derived ISA distribution.This situation may imply that ISA estimates using three indices have different uncertainties within the large area due to the unbalanced economic conditions and different climate conditions.has better performance when ISA is very small or high.When ISA is in medium range in a pixel (i.e., 0.4-0.6), the three indices have similar performance.When ISA is in the range of 0.2 and 0.4, NISI has slightly poor performance than HSI and VANUI.

Discussion
DMSP-OLS data have been widely used for ISA mapping at regional and global scales [24,54].However, the problems in mixed pixels (1 km spatial resolution), data saturation (0-63 data range), blooming effects, and the impacts of different socioeconomic conditions on DMSP-OLS radiances [29,30,40,41] have been recognized as major factors resulting in poor performance in extracting ISA distribution.Therefore, much research has explored the approaches to reduce these problems.One solution is to combine MODIS and DMSP-OLS data to produce a new product with improved spatial resolution and data ranges by taking the advantages of MODIS data in spatial resolution (250 m), radiometric resolution (12 bits), and good separation capacity in vegetation and non-vegetation (e.g., ISA, water) types.For example, previous research has proposed HSI [29], VANUI [40] and EANTLI [41] to reduce the data saturation problem, and has been proven effective in improving ISA estimation accuracy and spatial patterns.This research further confirmed the necessity of integrating DMSP-OLS and MODIS data in improving ISA mapping performance, and the newly developed NISI can further improve ISA mapping performance.This research indicates that the newly proposed NISI has better capability in ISA mapping than HSI and VANUI when ISA is very small (ISA is less than 0.2) and high (ISA is higher than 0.6).VANUI has the potential to overestimate ISA pixels with small values.This is because the mean values used in NDVI time series data sets cannot eliminate the impacts of clouds and bare soils, while the maximum values of NDVI time series data sets can solve this problem.Another advantage of using NISI is its capability in further improving the spatial structure in urban landscapes, thus NISI has better ISA prediction performance than HSI and VANUI when ISA is high.As global Landsat data are available, the integration of DMSP-OLS and Landsat NDVI will be another choice for improving ISA mapping [55].The advance of computing (e.g., Google Earth Engine, Google Cloud) and satellite technologies will make the effective use of multi-source data an important research topic in improving ISA mapping performance at regional and global scales [18].
The ISA estimation uncertainty varies in different cities (see Figure 6).Some cities in arid/semiarid regions such as Urumqi mainly suffered from an overestimation problem, while Dongguan-Shenzhen in south China and Chengdu in west China appeared to suffer from an underestimation problem.In addition to the economic conditions that directly affect DMSP-OLS radiance values, climate condition may be another factor indirectly affecting DMSP-OLS radiances.Previous research has shown that vegetation indices or abundance has highly negative correlation with ISA in urban landscapes, thus ISA can be estimated from the vegetation index [11].However, this relationship in arid and semiarid regions becomes poor because of limited vegetation covers in and around urban regions.In this situation, integration of vegetation indices and DMSP-OLS data cannot reduce the blooming effects and data saturation problems.To date, there is little research to explore the potential to accurately extract ISA in arid/semiarid regions.One solution may be use of temperature during the winter season when buildings in winter have higher temperature than surrounding bare soils due to human activities in urban areas.Use of MODIS temperature products may be valuable, but their coarse spatial resolution is problematic.An alternative is to use Landsat thermal bands [20].More research is needed to explore the incorporation of temperature data into DMSP-OLS or/and MODIS data to improve ISA mapping performance.ISA modeling based on different climate zones and/or economic conditions may be needed in the future for further improving estimation accuracy.

Conclusions
Regional and global ISA mapping has attracted much attention in the past two decades due to its importance in different applications such as urban environmental change, population and economic studies, and climate change.This research proposes a new index called NISI and compares it to previously used indices-HSI and VANUI-in ISA estimation in China using the SVR approach.This research has proven the promise of accurately mapping ISA distribution in a large area.Comparing with the two previously used indices, HSI and VANUI, the newly developed index can further improve ISA mapping performance, especially when the ISA proportion in a pixel is very small (e.g., <0.2) and relatively high (e.g., >0.6).This research confirmed the necessity to integrate DMSP-OLS and MODIS NDVI data to improve ISA estimation.More research should be explored to improve ISA estimation in arid and semiarid regions.

Figure 1 .
Figure 1.The study area-China and selected ten cities in red squares.Figure 1.The study area-China and selected ten cities in red squares.

Figure 1 .
Figure 1.The study area-China and selected ten cities in red squares.Figure 1.The study area-China and selected ten cities in red squares.

Figure 2 Figure 2
Figure 2 illustrates the framework of mapping fractional ISA distribution using the combination of DMSP-OLS, MODIS NDVI, and a limited number of Landsat 8 OLI data.The major steps include: (1) developing ISA data from Landsat 8 OLI data for the typical sites and aggregating the pixel ISA data at cell size of 25 m × 25 m to fractional ISA data at cell size of 250 m × 250 m; (2) preprocessing data for coarse spatial resolution images and implementing data integration using the indices; (3) developing ISA estimation models based on different indices using SVR; and (4) comparing and evaluating ISA estimation results.

Figure 2 .
Figure 2. Framework of mapping fractional ISA distribution through combination of multisource remote sensing data (ISA-impervious surface area; NISI-Normalized Impervious Surface Index; HSI-Human Settlement Index; VANUI-Vegetation Adjusted Nighttime light Urban Index; NDVI-Normalized Difference Vegetation Index).

Figure 2 .
Figure 2. Framework of mapping fractional ISA distribution through combination of multisource remote sensing data (ISA-impervious surface area; NISI-Normalized Impervious Surface Index; HSI-Human Settlement Index; VANUI-Vegetation Adjusted Nighttime light Urban Index; NDVI-Normalized Difference Vegetation Index).

Figure 3 .
Figure 3.A false color composite based on (a) OLSnor, NDVImax, and NISI as red, green, and blue; and corresponding digital values of (b) OLSnor, (c) NDVImax, and (d) NISI based on the line in the urban landscape of Beijing.(Note: OLSnor is the normalized DMSP-OLS image, and NDVImax is a composite of NDVI time series data using the maximum algorithm).

Figure 3 .
Figure 3.A false color composite based on (a) OLSnor, NDVImax, and NISI as red, green, and blue; and corresponding digital values of (b) OLSnor, (c) NDVImax, and (d) NISI based on the line in the urban landscape of Beijing.(Note: OLSnor is the normalized DMSP-OLS image, and NDVImax is a composite of NDVI time series data using the maximum algorithm).

Figure 4 .
Figure 4. Predicted impervious surface area distributions with the support vector regression method based on (a) HSI, (b) VANUI, and (c) NISI.Here (a1), (b1), and (c1) represent a typical location of impervious surface distribution corresponding to a, b, and c index respectively.

Figure 4 .
Figure 4. Predicted impervious surface area distributions with the support vector regression method based on (a) HSI, (b) VANUI, and (c) NISI.Here (a1), (b1), and (c1) represent a typical location of impervious surface distribution corresponding to a, b, and c index respectively.

Table 1 .
A summary of population, gross domestic product (GDP), and administrative area for the selected ten cities in 2012.

Table 2 .
Datasets used in research.

Table 3 .
Accuracy assessment results from Human Settlement Index (HSI), Vegetation Adjusted Nighttime light Urban Index (VANUI), and Normalized Impervious Surface Index (NISI) using Support Vector Regression (SVR).

Table 4 .
Evaluation of map applications by estimating the map means of fractional impervious surface area distributions.Note: Diff is the difference between estimates using different indices and mean from the samples.The smaller Diff values imply better estimation performance.

Table 5 .
Comparison of root mean squared error (RMSE) and relative RMSE (RMSEr) results among three variables based on impervious surface area (ISA) ranges.