Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Seasonal Variation of Colored Dissolved Organic Matter in Barataria Bay, Louisiana, Using Combined Landsat and Field Data
Previous Article in Journal
Mapping Two-Dimensional Deformation Field Time-Series of Large Slope by Coupling DInSAR-SBAS with MAI-SBAS
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mapping Impervious Surface Distribution with Integration of SNNP VIIRS-DNB and MODIS NDVI Data

1
School of Resource and Environmental Sciences, Wuhan University, Wuhan 430079, China
2
Zhejiang Provincial Key Laboratory of Carbon Cycling in Forest Ecosystems and Carbon Sequestration, School of Environmental & Resource Sciences, Zhejiang A & F University, Lin’an 311300, China
3
Institute of Photogrammetry and Remote Sensing, Chinese Academy of Surveying and Mapping, Beijing 100830, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2015, 7(9), 12459-12477; https://doi.org/10.3390/rs70912459
Submission received: 20 May 2015 / Revised: 13 September 2015 / Accepted: 15 September 2015 / Published: 22 September 2015

Abstract

:
Data from the U.S. Defense Meteorological Satellite Program’s Operational Line-scan System are often used to map impervious surface area (ISA) distribution at regional and global scales, but its coarse spatial resolution and data saturation produce high inaccuracy in ISA estimation. Suomi National Polar-orbiting Partnership (SNPP) Visible Infrared Imaging Radiometer Suite’s Day/Night Band (VIIRS-DNB) with its high spatial resolution and dynamic data range may provide new insights but has not been fully examined in mapping ISA distribution. In this paper, a new variable—Large-scale Impervious Surface Index (LISI)—is proposed to integrate VIIRS-DNB and Moderate Resolution Imaging Spectroradiometer (MODIS) normalized difference vegetation index (NDVI) data for mapping ISA distribution. A regression model was established, in which LISI was used as an independent variable and the reference ISA from Landsat images was a dependent variable. The results indicated a better estimation performance using LISI than using a single VIIRS-DNB or MODIS NDVI variable. The LISI-based approach provides accurate spatial patterns from high values in core urban areas to low values in rural areas, with an overall root mean squared error of 0.11. The LISI-based approach is recommended for fractional ISA estimation in a large area.

Graphical Abstract

1. Introduction

Impervious surface area (ISA) is defined as any man-made places such as buildings, highways, streets, and parking lots where water cannot penetrate [1,2,3,4]. Population increases and economic conditions have led to rapid ISA expansion in the past three decades, especially in China [5,6,7,8]. The rapid ISA increment has produced serious ecological and environmental problems [9,10] such as vegetation loss [11], water pollution [2], climate change [12,13,14], and urban heat islands [9,15] at local, regional, and even global scales. Therefore, increasing attention has focused on mapping ISA distribution in the past two decades [3,16,17].
Many studies on ISA mapping focused on individual cities using high or medium spatial resolution images such as QuickBird, IKONOS, Landsat, and ASTER [3]. Considering the labor intensity and cost in regional or global ISA mapping, coarse spatial resolution images such as the U.S. Defense Meteorological Satellite Program’s Operational Line-scan System (DMSP-OLS) and Moderate Resolution Imaging Spectroradiometer (MODIS) have been used in recent years [18,19,20]. Two types of variables—vegetation abundance and nighttime light data—are often used. Because of the inverse correlation between ISA and vegetation indices, MODIS normalized difference vegetation index (NDVI) has been used to map ISA in a large area [21,22,23,24]. However, vegetation distribution in an urban landscape is influenced by many factors such as terrain, climate, population, economic conditions, and cultures; thus, using these data alone for large-scale ISA mapping may generate high inaccuracy [25]. Since nighttime light data such as DMSP-OLS and the U.S. Suomi National Polar-orbiting Partnership (SNPP) Visible Infrared Imaging Radiometer Suite’s Day/Night Band (VIIRS-DNB) can effectively record human activities in a large area, they have been used for estimating population, energy consumption, economic data, and mapping ISA distribution [26,27,28,29,30,31]. Like vegetation data, the nighttime light data are also influenced by different factors such as economic conditions, mixed pixel problems, and data saturation (especially for DMSP-OLS) [32,33]. In previous research, DMSP-OLS was often used for mapping ISA at regional and global scales because these data have been available since 1992 at no cost [34,35,36,37]. However, the original DMSP-OLS data have a spatial resolution of 2.7 km and data range of 0–63. The mixed-pixel problem and data saturation in DMSP-OLS make it inaccurate in mapping ISA distribution in a large area. The SNNP VIIRS DNB, which was launched in October 2011, has many advantages over the DMSP-OLS data: improved spatial resolution (375 m and 750 m at nadir for VIIRS DNB, compared with 2.7 km for DMSP-OLS) and enlarged data range (14 bit for VIIRS DNB compared with 6 bit for DMSP-OLS), co-location with multispectral measurements on VIIRS and other NPOESS sensors, and elimination of cross-track pixel size variation [33,38,39]. The new generation of nighttime light data from VIIRS-DNB was first released by NOAA/NGDC in early 2013. Since then, this dataset has been employed to estimate electric power consumption and economic conditions [29,40,41,42]. It is believed that these data will be valuable for mapping ISA distribution in a large area and may provide more accurate ISA estimation than DMSP-OLS.
Much previous research explored the thresholding approach of the DMSP-OLS data for mapping ISA distribution in a large area [37,43,44,45]. Since Lu et al. used the combination of DMSP-OLS and MODIS NDVI data to estimate ISA in South and East China [25], scientists have explored the combined use of DMSP-OLS and MODIS NDVI or SPOT VGT and confirmed the promise of improving ISA mapping performance [6,46,47,48,49,50]. However, the coarse spatial resolution in MODIS and DMSP-OLS cannot be effectively used to produce details of spatial patterns. The improved spatial resolution in VIIRS-DNB and MODIS NDVI will provide new insights for ISA mapping in a large area. Therefore, the objective of this research was to propose a new variable through integrating VIIRS-DNB and MODIS NDVI data to improve ISA mapping performance in a large area.

2. Study Area and Datasets

2.1. Study Area

Since 2011, over half of China’s population has lived in cities (National Bureau of Statistics of China, 2012). China has been experiencing rapid urbanization during the past three decades [5,51,52]. It is necessary to update ISA distribution frequently because of its rapid dynamic change and the requirements of urban planning and management. We chose all of China as a study area because of its large area with wide variations in population densities and economic conditions, as well as different landscape patterns and cultural customs. We chose six typical cities—Beijing, Chengdu, Kunming, Shanghai, Wuhan, and Urumqi—as sample sites to verify variables’ versatility and accuracy (see Figure 1). Beijing and Shanghai are the two biggest megacities in China with the largest populations and gross domestic products (GDPs); Wuhan and Chengdu are megacities located in Central and West China, respectively; Kunming and Urumqi are located in Southwest and Northwest China, respectively. Table 1 provides a summary of populations, GDPs, and areas of the six cities.
Figure 1. The study area—all of China and six selected cities.
Figure 1. The study area—all of China and six selected cities.
Remotesensing 07 12459 g001
Table 1. Summary of populations, gross domestic products (GDPs), and total areas for the six cities.
Table 1. Summary of populations, gross domestic products (GDPs), and total areas for the six cities.
CitiesPopulation (Million)GDP (Billion RMB)Area (km2)
Beijing in North China20.691780.116,800.00
Shanghai in East China23.802010.16340.50
Wuhan in Central China10.12800.38494.41
Chengdu in Central China11.73813.812,390.00
Kunming in Southwest China6.53301.121,001.28
Urumqi in Northwest China3.35206.015,173.13

2.2. Datasets

Three types of remote sensing data—VIIRS-DNB (two-month composite product), MODIS 16-day NDVI composite (MOD13Q1), and Landsat 8 OLI (Operational Land Imager) data—were used in this research (see Table 2). The VIIRS-DNB and MODIS NDVI data were acquired in 2012, and Landsat 8 OLI images were acquired in 2013 because no Landsat data were available in 2012.
Table 2. Remote sensing data used in research.
Table 2. Remote sensing data used in research.
DataAcquisition DateDescription
VIIRS-DNBTwo-month composite product in April and October 2012A spectral range of 500–900 nm; highly sensitive to very low levels of visible light at night with zero moonlight; spatial resolution of 743 m.
MODIS NDVI (MOD13Q1)16-day MODIS NDVI composite between April and October 2012 (h23v04-h23v05, h24v04-h24v05, h25v03-h25v06, h26v03-h26v06, h27v04-h27v06, h28v05-h28v07, h29v06); total number of scenes: 247Gridded level-3 product with 250 m spatial resolution.
Landsat 8 OLI imagerypath/row: acquisition dateSix multispectral bands with 30 m and one panchromatic band with 15 m spatial resolution were used. Two thermal bands with 100 m spatial resolution were not used due to their relatively coarse spatial resolution.
123/32: 1 September 2013
118/38: 29 August 2013
123/39: 12 May 2013
129/39: 20 April 2013
129/43: 20 April 2013
143/29: 28 August 2013
Note: VIIRS-DNB data were downloaded from National Geophysical Data Center [53]; MODIS NDVI time-series data were downloaded from the NASA Goddard Space Flight Center (GSFC) [54]; and Landsat 8 OLI imagery was downloaded from the United States Geological Survey (USGS) (http://earthexplorer.usgs.gov/).

3. Methods

Figure 2 illustrates the framework for mapping ISA in a large area using a combined use of VIIRS-DNB and MODIS NDVI data. As a comparison, single VIIRS-DNB and MODIS NDVI variables were also used with the same samples. The major steps include: (1) produce ISA reference data from Landsat 8 OLI data for the selected cities; (2) develop a new index based on the combination of VIIRS-DNB and MODIS NDVI data; (3) establish regression models for ISA estimation; and (4) evaluate ISA estimates.
Figure 2. Framework of mapping ISA distribution using the Large-scale Impervious Surface Index (LISI).
Figure 2. Framework of mapping ISA distribution using the Large-scale Impervious Surface Index (LISI).
Remotesensing 07 12459 g002

3.1. Produce ISA Reference Data from Landsat 8 OLI Imagery

Landsat 8 OLI imagery covers 11 bands, including eight reflective bands (e.g., visible, near infrared, and shortwave infrared) with 30 m spatial resolution, one panchromatic band with 15 m spatial resolution, and two thermal infrared bands with 100 m spatial resolution [55]. Band 1 (violet-deep blue), band 9 (Cirrus), and bands 10 and 11 (thermal infrared) were not used in this research. The Universal Transverse Mercator (UTM) coordinate system in the Landsat 8 OLI data was re-projected to Albers Conical Equal Area projection. A combination of thresholding and cluster analysis was used to extract ISA data, as illustrated in Figure 3. The major steps include (1) conduct data fusion to produce a new dataset with improved spatial resolution; (2) produce vegetation indices to mask out vegetation and water; (3) extract spectral signatures for the remaining pixels and conduct cluster analysis; and (4) merge the clusters into ISA and others and evaluate ISA results.
The spatial resolution of remote sensing data is an important factor in effectively extracting ISA data because of the complex land-cover composition in urban landscapes [3]. Previous research has indicated that improved spatial resolution in Landsat data is valuable in ISA mapping and the wavelet-merging technique is an effective tool to integrate multispectral and panchromatic data into a new dataset [56] and thus this technique is used in this research.
Figure 3. Framework of extracting ISA data from Landsat 8 OLI imagery.
Figure 3. Framework of extracting ISA data from Landsat 8 OLI imagery.
Remotesensing 07 12459 g003
Many vegetation indices have been developed [21,57,58]. NDVI may be the most commonly used index for vegetation-related studies and NDWI (normalized difference water index) for water extraction [59,60]. Thus, they are used in this research to extract vegetation and water from the fused multispectral imagery. Based on analysis of vegetation samples, a threshold of 0.3 in NDVI data was selected to extract the vegetation pixels. Meanwhile, a threshold of 0.1 in NDWI was used to extract water pixels. After masking the vegetation and water pixels in the fused multispectral imagery, the remaining pixels are mainly ISA and bare soils. The spectral signatures from the fused imagery were extracted for the remaining pixels, and a cluster analysis (ISODATA in this research) was conducted to classify the remaining pixels into 100 clusters. The analyst was then to merge the clusters into ISA and others based on analysis of their spectral signature and QuickBird images in Google Earth. Finally, ISA results with 15 m spatial resolution were produced. Although no quantitative accuracy assessment for the ISA results was conducted, visual examination of the final ISA images by overlaying them on the corresponding Landsat 8 OLI color composite indicates the reliability and good quality; thus, the extracted ISA data from Landsat were used in modeling and accuracy assessment for the ISA estimation in a large area.

3.2. Develop Large-Scale Impervious Surface Index Data through a Combination of VIIRS-DNB and MODIS NDVI Data

The VIIRS-DNB data were re-projected from geographic (Lat/Lon) system into Albers Conical Equal Area projection and resampled to a cell size of 750 m by 750 m using the nearest-neighbor algorithm. Since the VIIRS-DNB data contain fires, gas flares, volcanoes, and background noises, they must be removed before using the data for ISA mapping. In this research, we used a threshold value of 0.5 to remove them. As the majority of data values were less than 65 (some had extremely high values) we set the threshold of 65 as the maximum, which means that all the pixel values greater than 65 were assigned to 65. In order to keep all data sources in the same range, between 0 and 1, the VIIRS-DNB image was normalized with Equation (1):
DNB nor = DNB DNB min DNB max DNB min
where DNB nor is a fractional VIIRS-DNB image having data ranges between 0 and 1,  DNB min and DNB max are the minimum and maximum values, respectively.
The MODIS NDVI data (MOD13Q1 product here: 16-day composite with 250 m spatial resolution) were re-projected from sinusoidal projection to Albers Conical Equal Area projection, and the nearest-neighbor resampling algorithm was used during the reprojection procedure. In theory, NDVI imagery has values ranging from −1 to +1. Due to the coarse spatial resolution in MODIS NDVI, the land surface covers, not including water bodies, in a large area have data ranges between 0 and 1 during the growing season.
Since non-vegetation land covers such as farmlands (e.g., bare soils after harvest), water, and ISA have similar NDVI values, it is necessary to reduce the confusion among them. One effective approach is to produce a maximum NDVI image from multi-temporal NDVI images [25]:
NDVI max = MAX [ NDVI 1 , NDVI 2 , , NDVI n ]
where NDVI 1 ,  NDVI 2 , ..., NDVI n are the multitemporal MOD13Q1 NDVI images acquired in 2012. Another important role for Equation (2) is to remove the impact of cloud contamination. Therefore, the final NDVImax imagery is cloud-free and has a data range between 0 and 1.
In urban landscapes, vegetation abundance is closely related to the patterns of settlements [61]; thus, vegetation indices have been used to estimate ISA. However, non-vegetation land covers, such as bare soils and water bodies within or outside the urban areas have NDVI values similar to ISA. Therefore, individual NDVI data are difficult to use directly for separating ISA from water and bare soils. Oppositely, the nighttime light data represent the urban area and have features considerably different from data for other regions, but these data are also influenced by different economic conditions [44,62]. Therefore, a combination of nighttime light and NDVI data has complementary features, as previous research has confirmed [25,46].
Two common combination variables—human settlement index [25] and vegetation adjusted normalized urban index [6,46]—have been proposed for mapping ISA distribution in a large area. Both variables are based on DMSP-OLS and MODIS NDVI with the same spatial resolution data (1000 m). Here we present a new combination variable called Large-scale Impervious Surface Index (LISI) based on VIIRS-DNB data with spatial resolution of 750 m and MODIS NDVImax with spatial resolution of 250 m:
LISI = ( 1 NDVI max ) × DNB nor
Since NDVI max is negatively related to ISA, 1- NDVI max is used to keep the values between 0 and 1 and is positively related to ISA. Compared to DMSP-OLS data, VIIRS-DNB data have high spatial resolution and a richly dynamic data range. Since some objects, such as airports and new, tall buildings [41] have very high DNB values, the expression DNB nor can smooth this effect. The combination of DNB nor and ( 1 NDVI max ) not only highlights urban ISA but also improves spatial resolution. As a comparison of these datasets, Figure 4 illustrates the improvement by combining both features into the new dataset.
Figure 4. A comparison of four datasets in Beijing City; (a) VIIRS-DNB with 750 m spatial resolution; (b) 1-NDVImax with 250 m spatial resolution; (c) LISI with improved spatial resolution (250 m); and (d) Landsat 8 OLI color composite with 15 m spatial resolution.
Figure 4. A comparison of four datasets in Beijing City; (a) VIIRS-DNB with 750 m spatial resolution; (b) 1-NDVImax with 250 m spatial resolution; (c) LISI with improved spatial resolution (250 m); and (d) Landsat 8 OLI color composite with 15 m spatial resolution.
Remotesensing 07 12459 g004

3.3. Map ISA Distribution with Regression Models

Regression analysis was used to develop ISA estimation models. The dependent variable—fractional ISA data—was obtained from the Landsat 8 OLI images, which were aggregated from a cell size of 15 m to 250 m and 750 m, respectively, using a mean algorithm to match the cell size of LISI and 1-NDVImax (250 m) and DNBnor (750 m). The independent variable is DNBnor, 1-NDVImax, and LISI, respectively. In order to develop the regression models, a random sampling technique was used to collect samples for each selected city. A total of 4800 samples were extracted from six reference images with 800 samples for each. The values for the same sample locations were extracted from DNBnor, 1-NDVImax, and LISI. Of these samples, 3600 samples were used for modeling and the remaining 1200 samples were used to evaluate ISA estimates. The coefficient of determination (R2) was used to evaluate the fitness of regression models. Meanwhile, three typical cities—Beijing, Wuhan, and Urumqi—were selected for a comparative analysis of their ISA spatial patterns, which were developed from the three regression models based on DNBnor, 1-NDVImax, and LISI.

3.4. Conduct Evaluation of ISA Estimates

Overall accuracy, kappa coefficient, and producer’s and user’s accuracies are often used for pixel-level classification evaluation [63]. However, these approaches are not suitable for the evaluation of fractional ISA estimates [3]. As in previous research [64], we used correlation coefficient (R) and root mean squared error (RMSE) to evaluate the results. In addition to the evaluation of ISA estimates in overall China, the accuracy assessment was also conducted at five ISA groups—very low, low, medium, high, and very high—based on 0.2 intervals of the ISA reference data (between 0 and 1), and at individual cities. The objectives of different accuracy assessment methods are to understand the error sources, whether they are from different levels of ISA values or from different spatial locations due to various population densities and economic conditions.

4. Results

4.1. Analysis of ISA Spatial Distribution

A comparison of R2 values among three regression models (see Figure 5) indicate that the LISI-based model has the highest performance and the 1-NDVImax-based model has the lowest performance. This is reasonable because NDVI is influenced by environmental and geographic factors such as bare soils, moisture, and compositions of different land-cover types. This study indicates that individual MODIS NDVI data are not suitable for ISA mapping in China because of the considerable difference of vegetation conditions in Eastern and Western China. Although using VIIRS DNB provides better estimation performance than NDVI, its estimation variation is wide due to the impacts of different economic conditions on nighttime light data. Figure 6 illustrates the ISA distribution using the LISI-based model, indicating that large amounts of ISA are distributed along the coastal regions and in the central metropolis of major cities, and much smaller amounts of ISA in Western and Northwestern China. This figure also clearly shows the spatial pattern differences between a high ISA proportion in core urban areas and a low ISA proportion in rural regions (see Figure 6a–c).
Figure 5. The relationships between reference (a) ISA and DNBnor, (b) ISA and 1-NDVImax, and (c) ISA and LISI.
Figure 5. The relationships between reference (a) ISA and DNBnor, (b) ISA and 1-NDVImax, and (c) ISA and LISI.
Remotesensing 07 12459 g005
Figure 6. ISA distribution in China using the LISI-based model, highlighting the ISA spatial patterns from high ISA proportions in core urban areas to low proportions in rural regions; (a) Urumqi; (b) Beijing; and (c) Wuhan.
Figure 6. ISA distribution in China using the LISI-based model, highlighting the ISA spatial patterns from high ISA proportions in core urban areas to low proportions in rural regions; (a) Urumqi; (b) Beijing; and (c) Wuhan.
Remotesensing 07 12459 g006
In order to better explain the ISA estimation performances, three cities—Beijing, Wuhan, and Urumqi—having different economic conditions and locations (see Table 1) were selected as examples. Beijing and Wuhan are located in North and Central China, respectively, with much better vegetation distributions than Urumqi in the west of China due to the latter’s dry weather. A comparison of ISA distribution (see Figure 7) indicates that the LISI-based model highlights the ISA spatial patterns from the highest values in core urban areas to the lowest values in rural areas, and the NDVI-based approach cannot effectively extract ISA distribution in urban regions, resulting in considerable underestimation when ISA density is high in an urban area. Figure 7 shows that the LISI-based method provides better ISA estimation performance than DNB- and NDVI-based approaches for each city.

4.2. Comparative Analysis of ISA Estimates

The overall accuracy assessment results in Table 3 also confirm that the LISI-based approach provides the best accuracy with R value of 0.81 and RMSE value of 0.11, much improved compared to the 1-NDVImax-based approach. The scatterplots between estimates and reference data (Figure 8) indicate a reasonably good estimation performance using the DNBnor-based approach and best estimation performance using the LISI-based approach. In contrast, the 1-NDVImax-based approach has very poor estimation performance because it cannot effectively estimate ISA values when the ISA proportion in a pixel is higher than 0.5.
Figure 7. A comparison of the ISA distributions from three datasets for three cities, highlighting the better spatial patterns from the LISI-based model than the other two results; (a1a3) represent ISA distribution using the DNBnor in Beijing, Wuhan, and Urumqi; (b1b3) represent ISA distribution using 1-NDVImax in Beijing, Wuhan, and Urumqi; and (c1c3) represent ISA distribution using LISI in Beijing, Wuhan, and Urumqi.
Figure 7. A comparison of the ISA distributions from three datasets for three cities, highlighting the better spatial patterns from the LISI-based model than the other two results; (a1a3) represent ISA distribution using the DNBnor in Beijing, Wuhan, and Urumqi; (b1b3) represent ISA distribution using 1-NDVImax in Beijing, Wuhan, and Urumqi; and (c1c3) represent ISA distribution using LISI in Beijing, Wuhan, and Urumqi.
Remotesensing 07 12459 g007
Table 3. A comparison of overall accuracy assessments among three datasets.
Table 3. A comparison of overall accuracy assessments among three datasets.
VariableRRMSE
DNBnor0.7290.132
1-NDVImax0.5630.160
LISI0.8120.113
Figure 8. The relationships between ISA estimates and corresponding reference data among three datasets; (a) ISA estimates from DNBnor; (b) ISA estimates from 1-NDVImax; and (c) ISA estimates from LISI.
Figure 8. The relationships between ISA estimates and corresponding reference data among three datasets; (a) ISA estimates from DNBnor; (b) ISA estimates from 1-NDVImax; and (c) ISA estimates from LISI.
Remotesensing 07 12459 g008
The analysis above is based on overall estimation performance, but may not reveal some details caused by different environmental or economic conditions. Table 4 summarizes RMSE results, which are grouped into five ISA levels—very low, low, medium, high, and very high—based on reference data ranges: <0.2, [0.2–0.4), [0.4–0.6), [0.6–0.8), and ≥0.8. These results confirmed that the LISI-based approach produces better performance than the other two. Comparing the RMSE values at different ISA levels indicates that the LISI-based approach improves ISA estimation performance when ISA is relatively low or high, but when ISA is at medium level (i.e., between 0.4 and 0.6), the DNBnor-based approach provides better estimation performance than the LISI-based approach. Table 4 indicates that a single 1-NDVImax variable is not suitable for ISA estimation in China, DNBnor is valuable when ISA is at medium level, and LISI is recommended for ISA estimation in a large area.
Table 4. Comparison of RMSE results among three estimation approaches at five ISA levels.
Table 4. Comparison of RMSE results among three estimation approaches at five ISA levels.
GroupData RangeRMSE
DNBnor1-NDVImaxLISI
Very low<0.20.1000.1210.076
Low0.2–0.40.2170.1880.181
Medium0.4–0.60.1820.4610.195
High0.6–0.80.226Null0.179
Very high≥0.80.310Null0.215
Overall 0.1320.1600.113
Note: Null represents there are no data in ISA estimates.
R and RMSE were also used to evaluate the estimation performance at individual cities, and the results are summarized in Table 5. Again, the LISI-based approach provided better estimation performance than the other two variables in all cities, implying the robustness of the LISI variable in ISA estimation. The 1-NDVImax variable has good performance in Beijing and Chengdu, but very poor performance in Urumqi, implying that the 1-NDVImax variable is not suitable for ISA estimation in cities where vegetation accounts for a limited proportion and bare soils/desert have important impacts on the NDVI values. The DNBnor-based approach provides reliable estimation in different cities, but relatively poorer performance than the LISI-based approach.
Table 5. Comparison of accuracy analysis results from three variables at different cities.
Table 5. Comparison of accuracy analysis results from three variables at different cities.
CityDNBnor1-NDVImaxLISI
RRMSERRMSERRMSE
Beijing0.7440.1700.8450.1360.8830.119
Shanghai0.8040.1050.6610.1330.8050.105
Wuhan0.7930.0960.7400.1060.8400.086
Chengdu0.8010.1040.8010.1040.8440.093
Kunming0.7390.1090.7640.1040.8220.092
Urumqi0.8140.0760.2330.1270.8520.068

5. Discussions

DMSP-OLS data are often used for ISA estimation in a large area using a thresholding-based approach [44,45,65], but a pixel-based approach produces high inaccuracy in ISA estimation due to the mixed-pixel problem (2.7 km spatial resolution for the original data) and data saturation (6 bits) [33]. Since the DMSP-OLS data are seriously influenced by different economic conditions, no suitable thresholds can be used to accurately extract ISA data, thus, the spatial patterns of extracted ISA distribution are often poor [25]. The VIIRS-DNB data considerably reduce the problems in DMSP-OLS because of the improvement in spatial resolution and data ranges; thus, the ISA estimation using VIIRS-DNB for ISA estimation in China has reasonably good estimation performance, with an overall RMSE of 0.13, as shown in this research.
Due to the problems inherited in DMSP-OLS data, previous research has explored the combination of DMSP-OLS and MODIS NDVI data for ISA mapping. Lu et al. (2008) first proposed the human settlement index based on the combination of DMSP-OLS and MODIS NDVI data for mapping human settlements in China and obtained much better estimation performance (including the improved spatial patterns and estimation accuracy) than using individual DMSP-OLS data [25]. Zhang et al. (2013) further proposed the vegetation adjusted normalized urban index to improve the ISA estimation [46]. The combined indices are based on coarse spatial resolution images (e.g., 1 km spatial resolution) and small data ranges (0–63 in DMSP-OLS data), thus, the improvement is influenced by the data limitations themselves. This research proposed the LISI variable—an integration of VIIRS-DNB and MODIS NDVI—which has provided much improvement in ISA estimation because of the integration of the improved spatial resolution and data ranges in VIIRS-DNB data and MODIS NDVI with 250 m spatial resolution. The LISI variable combined the advantages of both VIIRS-DNB and MODIS NDVI data; that is, VIIRS-DNB can effectively reflect the difference between urban and non-urban regions, and MODIS NDVI can better capture the differences within the urban inner structure (e.g., the different compositions of ISA, vegetation, and water). Additionally, this research indicates that the LISI-based approach can provide robust ISA estimation in different cities, and it is especially valuable in improving the ISA estimation when the ISA proportion in a pixel is relatively low or high.
Previous research using DMSP-OLS for ISA mapping has indicated that unbalanced economic conditions in a large area affect the nighttime light values; therefore, multiple thresholds are often used to extract ISA data for different cities (or locations) [44]. However, this research shows that the VIIRS-DNB reduced this problem. Even when using the same DNB-based model for ISA mapping in China, the different economic conditions such as in Shanghai, Chengdu, and Urumqi did not considerably affect the ISA estimation performance, as shown in Table 5. When combining the RMSE (Table 5) and economic conditions (Table 1) for the same city, a scatterplot showing the relationships between RMSE and GDP indicates that economic condition is indeed an important factor affecting the estimation errors (see Figure 9). For example, Beijing and Shanghai have the largest populations and GDPs in China, but their different characteristics in spatial patterns of urban landscapes and different land-cover compositions affect the ISA estimation performance. Use of individual VIIRS-DNB data has relatively high estimation errors when economic condition is very high or relatively low, as shown in Figure 9. However, the LISI can considerably reduce this problem and is, thus, recommended for ISA mapping in a large area. We encourage researchers to explore the use of LISI in different countries. Furthermore, this research only explores the ISA estimation using a linear regression model, and more research is needed in the future to explore nonparametric algorithms such as a support vector machine or neural network for developing ISA estimation models.
Figure 9. The relationship between RMSE and economic conditions, showing the effects of different economic conditions on ISA estimation performance through the comparison of different variables.
Figure 9. The relationship between RMSE and economic conditions, showing the effects of different economic conditions on ISA estimation performance through the comparison of different variables.
Remotesensing 07 12459 g009

6. Conclusions

Through comparative analysis of VIIRS-DNB, MODIS NDVI, and the proposed LISI variable for ISA mapping in a large area, we obtained the following conclusions:
(1)
VIIRS-DNB can be used for ISA mapping in a large area with an overall RMSE of 0.13. However, the areas having higher ISA proportion produced higher errors. Additionally, very high or very low economic conditions influenced ISA estimation performance. This implies that individual VIIRS-DNB data may produce inaccurate spatial patterns of ISA distribution if the study area covers urban landscapes having considerably different economic conditions;
(2)
Individual MODIS NDVI is not a good variable for ISA mapping in a large area, especially in areas with very low vegetation covers, such as Western China. However, in some large cities such as Chengdu and Kunming in this research, NDVI can produce ISA estimates with similar to or even better performance than VIIRS-DNB. This implies that MODIS NDVI is valuable, but it is critical to properly use it in ISA estimation;
(3)
The proposed LISI variable combined advantages of both VIIRS-DNB and MODIS NDVI features and provided much improved ISA estimation performance, especially the improved spatial patterns. Overall, the LISI-based approach has an RMSE of 0.11 and has much-improved estimation performance when ISA proportion is high, compared to the other two datasets. Therefore, LISI is recommended for ISA estimation in a large area.

Acknowledgments

The study was financially supported by the Zhejiang A & F University’s Research and Development Fund for the talent startup project (2013FR052) and facility support from the Center for Global Change and Earth Observations at Michigan State University. The authors would like to thank Changming Zhu for his assistance in data processing.

Author Contributions

Wei Guo is responsible for the data collection and processing, resultant analysis, and preparation of manuscript. Dengsheng Lu is responsible for the design of the study, technical support, resultant analysis, and preparation of the manuscript. Yanlan Wu and Jixian Zhang contributed to the manuscript.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Brabec, E.; Schulte, S.; Richards, P.L. Impervious surfaces and water quality: A review of current literature and its implications for watershed planning. J. Plan. Lit. 2002, 16, 499–514. [Google Scholar] [CrossRef]
  2. Arnold, C.L.; Gibbons, C.J. Impervious surface coverage: The emergence of a key environmental indicator. J. Am. Plan. Assoc. 1996, 62, 243–258. [Google Scholar] [CrossRef]
  3. Lu, D.; Li, G.; Kuang, W.; Moran, E. Methods to extract impervious surface areas from satellite images. Int. J. Digit. Earth 2014, 7, 93–112. [Google Scholar] [CrossRef]
  4. Wu, C.; Murray, A.T. Estimating impervious surface distribution by spectral mixture analysis. Remote Sens. Environ. 2003, 84, 493–505. [Google Scholar] [CrossRef]
  5. Kuang, W.; Liu, J.; Zhang, Z.; Lu, D.; Xiang, B. Spatiotemporal dynamics of impervious surface areas across china during the early 21st century. Chin. Sci. Bull. 2012, 58, 1691–1701. [Google Scholar] [CrossRef]
  6. Ma, Q.; He, C.; Wu, J.; Liu, Z.; Zhang, Q.; Sun, Z. Quantifying spatiotemporal patterns of urban impervious surfaces in China: An improved assessment using nighttime light data. Landsc. Urban Plan. 2014, 130, 36–49. [Google Scholar] [CrossRef]
  7. Schneider, A.; Mertes, C.M.; Tatem, A.J.; Tan, B.; Sulla-Menashe, D.; Graves, S.J.; Patel, N.N.; Horton, J.A.; Gaughan, A.E.; Rollo, J.T.; et al. A new urban landscape in East-Southeast Asia, 2000–2010. Environ. Res. Lett. 2015, 10. [Google Scholar] [CrossRef]
  8. Schneider, A.; Chang, C.; Paulsen, K. The changing spatial form of cities in Western China. Landsc. Urban Plan. 2015, 135, 40–61. [Google Scholar] [CrossRef]
  9. Imhoff, M.L.; Zhang, P.; Wolfe, R.E.; Bounoua, L. Remote sensing of the urban heat island effect across biomes in the continental USA. Remote Sens. Environ. 2010, 114, 504–513. [Google Scholar] [CrossRef]
  10. Kaufmann, R.K.; Seto, K.C.; Schneider, A.; Liu, Z.; Zhou, L.; Wang, W. Climate response to rapid urban growth: Evidence of a human-induced precipitation deficit. J. Clim. 2007, 20, 2299–2306. [Google Scholar] [CrossRef]
  11. Wu, J. Urban sustainability: An inevitable goal of landscape research. Landsc. Ecol. 2009, 25, 1–4. [Google Scholar] [CrossRef]
  12. Carlson, T.N.; Traci Arthur, S. The impact of land use—Land cover changes due to urbanization on surface microclimate and hydrology: A satellite perspective. Glob. Planet. Chang. 2000, 25, 49–65. [Google Scholar] [CrossRef]
  13. Jennings, D.; Taylor Jarnagin, S. Changes in anthropogenic impervious surfaces, precipitation and daily streamflow discharge: A historical perspective in a mid-Atlantic subwatershed. Landsc. Ecol. 2002, 17, 471–489. [Google Scholar] [CrossRef]
  14. Xian, G. Satellite remotely-sensed land surface parameters and their climatic effects for three metropolitan regions. Adv. Space Res. 2008, 41, 1861–1869. [Google Scholar] [CrossRef]
  15. Buyantuyev, A.; Wu, J. Urban heat islands and landscape heterogeneity: Linking spatiotemporal variations in surface temperatures to land-cover and socioeconomic patterns. Landsc. Ecol. 2009, 25, 17–33. [Google Scholar] [CrossRef]
  16. Schneider, A. Monitoring land cover change in urban and peri-urban areas using dense time stacks of Landsat satellite data and a data mining approach. Remote Sens. Environ. 2012, 124, 689–704. [Google Scholar] [CrossRef]
  17. Weng, Q. Remote sensing of impervious surfaces in the urban areas: Requirements, methods, and trends. Remote Sens. Environ. 2012, 117, 34–49. [Google Scholar] [CrossRef]
  18. Schneider, A.; Friedl, M.A.; Potere, D. Mapping global urban areas using MODIS 500-m data: New methods and datasets based on “urban ecoregions”. Remote Sens. Environ. 2010, 114, 1733–1746. [Google Scholar] [CrossRef]
  19. Schneider, A.; Friedl, M.A.; Potere, D. A new map of global urban extent from MODIS satellite data. Environ. Res. Lett. 2009, 4. [Google Scholar] [CrossRef]
  20. Mertes, C.M.; Schneider, A.; Sulla-Menashe, D.; Tatem, A.J.; Tan, B. Detecting change in urban areas at continental scales with MODIS data. Remote Sens. Environ. 2015, 158, 331–347. [Google Scholar] [CrossRef]
  21. Huete, A.; Didan, K.; Miura, T.; Rodriguez, E.P.; Gao, X.; Ferreira, L.G. Overview of the radiometric and biophysical performance of the MODIS vegetation indices. Remote Sens. Environ. 2002, 83, 195–213. [Google Scholar] [CrossRef]
  22. Knight, J.; Voth, M. Mapping impervious cover using multi-temporal MODIS NDVI data. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2011, 4, 303–309. [Google Scholar] [CrossRef]
  23. Yang, F.; Matsushita, B.; Fukushima, T.; Yang, W. Temporal mixture analysis for estimating impervious surface area from multi-temporal MODIS NDVI data in Japan. ISPRS J. Photogramm. Remote Sens. 2012, 72, 90–98. [Google Scholar] [CrossRef] [Green Version]
  24. Yang, F.; Matsushita, B.; Yang, W.; Fukushima, T. Mapping the human footprint from satellite measurements in Japan. ISPRS J. Photogramm. Remote Sens. 2014, 88, 80–90. [Google Scholar] [CrossRef] [Green Version]
  25. Lu, D.; Tian, H.; Zhou, G.; Ge, H. Regional mapping of human settlements in Southeastern China with multisensor remotely sensed data. Remote Sens. Environ. 2008, 112, 3668–3679. [Google Scholar] [CrossRef]
  26. Sutton, P.; Roberts, D.; Elvidge, C.; Baugh, K. Census from heaven: An estimate of the global human population using night-time satellite imagery. Int. J. Remote Sens. 2010, 22, 3061–3076. [Google Scholar] [CrossRef]
  27. Cao, C.; Bai, Y. Quantitative analysis of VIIRS DNB nightlight point source for light power estimation and stability monitoring. Remote Sens. 2014, 6, 11915–11935. [Google Scholar] [CrossRef]
  28. Cao, X.; Wang, J.; Chen, J.; Shi, F. Spatialization of electricity consumption of China using saturation-corrected DMSP-OLS data. Int. J. Appl. Earth Obs. Geoinf. 2014, 28, 193–200. [Google Scholar] [CrossRef]
  29. Shi, K.; Yu, B.; Huang, Y.; Hu, Y.; Yin, B.; Chen, Z.; Chen, L.; Wu, J. Evaluating the ability of NPP-VIIRS nighttime light data to estimate the gross domestic product and the electric power consumption of china at multiple scales: A comparison with DMSP-OLS data. Remote Sens. 2014, 6, 1705–1724. [Google Scholar] [CrossRef]
  30. Elvidge, C.D.; Baugh, K.E.; Anderson, S.J.; Sutton, P.C.; Ghosh, T. The night light development index (NLDI): A spatially explicit measure of human development from satellite data. Soc. Geogr. 2012, 7, 23–35. [Google Scholar] [CrossRef]
  31. Elvidge, C.; Safran, J.; Tuttle, B.; Sutton, P.; Cinzano, P.; Pettit, D.; Arvesen, J.; Small, C. Potential for global mapping of development via a nightsat mission. GeoJournal 2007, 69, 45–53. [Google Scholar] [CrossRef]
  32. Elvidge, C.D.; Erwin, E.H.; Baugh, K.E.; Ziskin, D.; Tuttle, B.T.; Ghosh, T.; Sutton, P.C. Overview of DMSP nightime lights and future possibilities. In Proceedings of Urban Remote Sensing Event, 2009 Joint, Shanghai, China, 20–22 May 2009; pp. 1–5.
  33. Elvidge, C.D.; Baugh, K.E.; Zhizhin, M.; Hsu, F.C. Why VIIRS data are superior to DMSP for mapping nighttime lights. In Proceedings of the Asia-Pacific Advanced Network, Boulder, CO, USA; 2013; Volume 35, p. 62. [Google Scholar]
  34. Sutton, P.C. A scale-adjusted measure of “urban sprawl” using nighttime satellite imagery. Remote Sens. Environ. 2003, 86, 353–369. [Google Scholar] [CrossRef]
  35. Elvidge, C.; Tuttle, B.; Sutton, P.; Baugh, K.; Howard, A.; Milesi, C.; Bhaduri, B.; Nemani, R. Global distribution and density of constructed impervious surfaces. Sensors 2007, 7, 1962–1979. [Google Scholar] [CrossRef]
  36. Zhang, Q.; Seto, K.C. Mapping urbanization dynamics at regional and global scales using multi-temporal DMSP/OLS nighttime light data. Remote Sens. Environ. 2011, 115, 2320–2329. [Google Scholar] [CrossRef]
  37. Zhou, Y.; Smith, S.J.; Elvidge, C.D.; Zhao, K.; Thomson, A.; Imhoff, M. A cluster-based method to map urban area from DMSP/OLS nightlights. Remote Sens. Environ. 2014, 147, 173–185. [Google Scholar] [CrossRef]
  38. Lee, T.E.; Miller, S.D.; Turk, F.J.; Schueler, C.; Julian, R.; Deyo, S.; Dills, P.; Wang, S. The NPOESS VIIRS day/night visible sensor. Bull. Am. Meteorol. Soc. 2006, 87, 191–199. [Google Scholar] [CrossRef]
  39. Miller, S.D.; Mills, S.P.; Elvidge, C.D.; Lindsey, D.T.; Lee, T.F.; Hawkins, J.D. Suomi satellite brings to light a unique frontier of nighttime environmental sensing capabilities. Proc. Natl. Acad. Sci. USA 2012, 109, 15706–15711. [Google Scholar] [CrossRef] [PubMed]
  40. Li, X.; Xu, H.; Chen, X.; Li, C. Potential of NPP-VIIRS nighttime light imagery for modeling the regional economy of China. Remote Sens. 2013, 5, 3057–3081. [Google Scholar] [CrossRef]
  41. Ma, T.; Zhou, C.; Pei, T.; Haynie, S.; Fan, J. Responses of Suomi-NPP VIIRS-derived nighttime lights to socioeconomic activity in China’s cities. Remote Sens. Lett. 2014, 5, 165–174. [Google Scholar] [CrossRef]
  42. Ma, T.; Zhou, Y.; Wang, Y.; Zhou, C.; Haynie, S.; Xu, T. Diverse relationships between Suomi-NPP VIIRS night-time light and multi-scale socioeconomic activity. Remote Sens. Lett. 2014, 5, 652–661. [Google Scholar] [CrossRef]
  43. Imhoff, M.L.; Lawrence, W.T.; Elvidge, C.; Paul, T.; Levine, E.; Prevalsky, M.; Brown, V. Using nighttime DMSP/OLS images of city lights to estimate the impact of urban land use on soil resources in the U.S. Remote Sens. Environ. 1997, 59, 105–117. [Google Scholar] [CrossRef]
  44. Henderson, M.; Yeh, E.T.; Gong, P.; Elvidge, C.; Baugh, K. Validation of urban boundaries derived from global night-time satellite imagery. Int. J. Remote Sens. 2003, 24, 595–609. [Google Scholar] [CrossRef]
  45. Lawrence, W.T.; Imhoff, M.L.; Kerle, N.; Stutzer, D. Quantifying urban land use and impact on soils in Egypt using diurnal satellite imagery of the earth surface. Int. J. Remote Sens. 2010, 23, 3921–3937. [Google Scholar] [CrossRef]
  46. Zhang, Q.; Schaaf, C.; Seto, K.C. The vegetation adjusted NTL urban index: A new approach to reduce saturation and increase variation in nighttime luminosity. Remote Sens. Environ. 2013, 129, 32–41. [Google Scholar] [CrossRef]
  47. Cao, X.; Chen, J.; Imura, H.; Higashi, O. A SVM-based method to extract urban areas from DMSP-OLS and spot VGT data. Remote Sens. Environ. 2009, 113, 2205–2209. [Google Scholar] [CrossRef]
  48. Pandey, B.; Joshi, P.K.; Seto, K.C. Monitoring urbanization dynamics in India using DMSP/OLS night time lights and spot-VGT data. Int. J. Appl. Earth Obs. Geoinf. 2013, 23, 49–61. [Google Scholar] [CrossRef]
  49. Shao, Z.; Liu, C. The integrated use of DMSP-OLS nighttime light and MODIS data for monitoring large-scale impervious surface dynamics: A case study in the Yangtze River Delta. Remote Sens. 2014, 6, 9359–9378. [Google Scholar] [CrossRef]
  50. Zhang, X.; Li, P.; Cai, C. Regional urban extent extraction using multi-sensor data and one-class classification. Remote Sens. 2015, 7, 7671–7694. [Google Scholar] [CrossRef]
  51. Kuang, W.; Chi, W.; Lu, D.; Dou, Y. A comparative analysis of megacity expansions in China and the U.S.: Patterns, rates and driving forces. Landsc. Urban Plan. 2014, 132, 121–135. [Google Scholar] [CrossRef]
  52. Ma, T.; Zhou, Y.; Zhou, C.; Haynie, S.; Pei, T.; Xu, T. Night-time light derived estimation of spatio-temporal characteristics of urbanization dynamics using DMSP/OLS satellite data. Remote Sens. Environ. 2015, 158, 453–464. [Google Scholar] [CrossRef]
  53. VIIRS-DNB. Two-Month Composite Product. Available online: http://ngdc.noaa.gov/eog/viirs/download_monthly.hOLIl (accessed on 9 February 2014).
  54. MODIS NDVI (MOD13Q1). Available online: http://ladsweb.nascom.nasa.gov/data/search.hOLIl (accessed on 2 December 2013).
  55. Landsat 8 OLI Description. Available online: http://landsat.usgs.gov/landsat8.php (accessed on 19 April 2015).
  56. Lu, D.; Li, G.; Moran, E.; Batistella, M.; Freitas, C.C. Mapping impervious surfaces with the integrated use of Landsat thematic mapper and radar data: A case study in an urban-rural landscape in the Brazilian Amazon. ISPRS J. Photogramm. Remote Sens. 2011, 66, 798–808. [Google Scholar] [CrossRef] [Green Version]
  57. Bannari, A.; Morin, D.; Bonn, F.; Huete, A.R. A review of vegetation indices. Remote Sens. Rev. 1995, 13, 95–120. [Google Scholar] [CrossRef]
  58. McDonald, A.J.; Gemmell, F.M.; Lewis, P.E. Investigation of the utility of spectral vegetation indices for determining information on coniferous forests. Remote Sens. Environ. 1998, 66, 250–272. [Google Scholar] [CrossRef]
  59. McFeeters, S.K. The use of the normalized difference water index (NDWI) in the delineation of open water features. Int. J. Remote Sens. 1996, 17, 1425–1432. [Google Scholar] [CrossRef]
  60. Xu, H. Modification of normalised difference water index (NDWI) to enhance open water features in remotely sensed imagery. Int. J. Remote Sens. 2006, 27, 3025–3033. [Google Scholar] [CrossRef]
  61. Weng, Q.; Lu, D.; Schubring, J. Estimation of land surface temperature—Vegetation abundance relationship for urban heat island studies. Remote Sens. Environ. 2004, 89, 467–483. [Google Scholar] [CrossRef]
  62. Small, C.; Pozzi, F.; Elvidge, C.D. Spatial analysis of global urban extent from DMSP-OLS night lights. Remote Sens. Environ. 2005, 96, 277–291. [Google Scholar] [CrossRef]
  63. Congalton, R.G.; Green, K. Assessing the Accuracy of Remotely Sensed Data: Principles and Practices, 2nd ed.; CRC Press: Boca Raton, FL, USA, 2009. [Google Scholar]
  64. Lu, D.; Weng, Q. Spectral mixture analysis of aster images for examining the relationship between urban thermal features and biophysical descriptors in Indianapolis, Indiana, USA. Remote Sens. Environ. 2006, 104, 157–167. [Google Scholar] [CrossRef]
  65. Imhoff, M.; Lawrence, W.T.; Stutzer, D.C.; Elvidge, C.D. A technique for using composite DMSP/OLS “city lights” satellite data to map urban area. Remote Sens. Environ. 1997, 61, 361–370. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Guo, W.; Lu, D.; Wu, Y.; Zhang, J. Mapping Impervious Surface Distribution with Integration of SNNP VIIRS-DNB and MODIS NDVI Data. Remote Sens. 2015, 7, 12459-12477. https://doi.org/10.3390/rs70912459

AMA Style

Guo W, Lu D, Wu Y, Zhang J. Mapping Impervious Surface Distribution with Integration of SNNP VIIRS-DNB and MODIS NDVI Data. Remote Sensing. 2015; 7(9):12459-12477. https://doi.org/10.3390/rs70912459

Chicago/Turabian Style

Guo, Wei, Dengsheng Lu, Yanlan Wu, and Jixian Zhang. 2015. "Mapping Impervious Surface Distribution with Integration of SNNP VIIRS-DNB and MODIS NDVI Data" Remote Sensing 7, no. 9: 12459-12477. https://doi.org/10.3390/rs70912459

Article Metrics

Back to TopTop