Next Article in Journal
Analysis and Experimental Validation of Efficient Coded OFDM for an Impulsive Noise Environment
Next Article in Special Issue
Initial Assessment of the LEO Based Navigation Signal Augmentation System from Luojia-1A Satellite
Previous Article in Journal
A Smart Pillow for Health Sensing System Based on Temperature and Humidity Sensors
Previous Article in Special Issue
Potentiality of Using Luojia 1-01 Nighttime Light Imagery to Investigate Artificial Light Pollution
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mapping Urban Extent Using Luojia 1-01 Nighttime Light Imagery

1
State Key Laboratory of Information Engineering in Surveying, Mapping and Remote Sensing, Wuhan University, Wuhan 430079, China
2
Collaborative Innovation Centre of Geospatial Technology, Wuhan 430079, China
3
School of Economics, Wuhan Donghu University, Wuhan 430212, China
4
Key Laboratory of the Ministry of Land and Resources for Law Evaluation Engineering, Wuhan 430074, China
*
Authors to whom correspondence should be addressed.
Sensors 2018, 18(11), 3665; https://doi.org/10.3390/s18113665
Submission received: 29 September 2018 / Revised: 22 October 2018 / Accepted: 24 October 2018 / Published: 29 October 2018
(This article belongs to the Special Issue The Design, Data Processing and Applications of Luojia 1-01 Satellite)

Abstract

:
Luojia 1-01 satellite, launched on 2 June 2018, provides a new data source of nighttime light at 130 m resolution and shows potential for mapping urban extent. In this paper, using Luojia 1-01 and VIIRS nighttime light imagery, we compared several methods for extracting urban areas, including Human Settlement Index (HSI), Simple Thresholding Segmentation (STS) and SVM supervised classification. According to the accuracy assessment, the HSI method using LJ1-01 data had the best performance in urban extent extraction, which presented the largest Kappa Coefficient value, 0.834, among all the results. For the urban areas extracted by VIIRS based HSI method, the largest Kappa Coefficient value was 0.772. In contrast, the largest Kappa Coefficient values obtained by STS method were 0.79 and 0.7512 respectively when using LJ1-01 and VIIRS data, while for SVM method the values were 0.7829 and 0.7486 when using Landsat-LJ and Landsat-VIIRS composite data respectively. The experimented results demonstrated that the utilization of nighttime light imagery can largely improve the accuracy of urban extent extraction and LJ1-01 data, with a higher resolution and more abundant spatial information, can lead to better identification results than its predecessors.

1. Introduction

Over the past half-century, most regions of the world, especially developing countries, have experienced a phase of high-speed urbanization [1,2,3]. In China, since the implementation of the reform and opening-up policy, the urbanization rate has increased by about 1% per year, from 17.9% in 1978 to 58.5% in 2017 [4]. The rapid development of urbanization has brought dramatic changes in urban areas, as well as a series of regional socio-economic and ecological problems caused by disorderly urban sprawl [5,6,7,8,9,10,11]. Therefore, accurate and timely measurement of urban areas, which is vital for analyzing urbanization dynamics and controlling problems mentioned above, has become an important topic in urban research [12,13,14,15,16].
As an objective and real-time data source, remote sensing satellite images have been widely used in urban area mapping and urbanization monitoring. Currently, coarse resolution images, such as the nighttime light (NTL) data have proven effective in urban extent extraction at regional and global scales [17,18,19,20,21]. At cloud-free nights, the visible light images of the earth’s surface detected by remote sensors are remotely sensed NTL images. Unlike daytime remote sensing images, the city lights recorded by the NTL images are closely related to human activities [22]. Therefore, it provides a unique perspective for socio-economic analysis and has been widely used in the fields of regional economy [23,24,25,26,27], urbanization [28,29,30], power consumption estimation [31], conflict assessment [32,33,34] and so forth. According to previous research, the Defense Meteorological Satellite Program’s Operational Linescan System (DMSP/OLS) data and the Day/Night Band (DNB) in the Visible Infrared Imaging Radiometer Suite (VIIRS) on the Suomi National Polar-orbiting Partnership Satellite are the two most commonly used NTL data sources [35,36]. Launched in the 1970s, the DMSP/OLS has a rich archive of data but it also has some shortcomings such as coarse spatial resolution (~2.7 km), blooming effects, saturation in urban cores and lack of on-board calibration, which may lead to misestimates of urban areas [31,37,38,39]. While the VIIRS DNB data, as a successor to the DMSP/OLS, has provided significant improvements such as finer spatial resolution (~740 m), on-board calibration and wider radiometric measurement range, largely reducing the saturation and blooming problems of the DMSP/OLS data [40,41,42,43].
Generally, the urban extent extraction methods using NTL data can be divided into three types, that is, thresholding-based, classification-based and index-based methods [44,45,46]. The thresholding-based method selects suitable values (local or global) to segment images [45,47,48]. As single-threshold has been proven to be problematic for urban areas at different levels of development, a multiple-threshold method is widely adopted in urban extent extraction. Shi et al. [49] utilized statistical data from the government as reference and when the urban areas extracted results and the reference data produced the minimum difference, the DNB values used for extraction were determined as the optimal thresholds. Zhou et al. [50] developed a cluster-based method to map urban areas, where a logistic regression model was used for potential urban clusters and optimal thresholds were estimated by cluster size and overall NTL magnitude. The classification-based method is to recognize NTL data as a kind of grayscale image and apply it to various classification algorithms. Jing et al. [51] utilized the integration of NTL and MODIS data and evaluated the accuracy of four classification algorithms for urban extent extraction. Xu et al. [52] and He et al. [53] mapped urban areas in China using the integrated NTL, normalized difference vegetation index (NDVI) and land surface temperature (LST) support vector machine (SVM) classification method. Nevertheless, the accuracy of classification methods is largely influenced by the selection of training samples, while solid a priori knowledge is not always available in the majority of cases. Besides, this process is labor-intensive and time-consuming and the accuracy of extraction results may suffer from the design of class schemes [54,55]. In the way of index-based methods, Lu et al. [56] proposed the human settlement index (HSI) by integrating DMSP/OLS NTL and MODIS NDVI data to extract urban areas in southeastern China. Zhang et al. [57] developed a vegetation adjusted NTL urban index (VANUI) to study urban structures on the global scale. Liu et al. [58] proposed a normalized urban area composite index (NUACI), involving the normalized difference water index (NDWI) and enhanced vegetation index (EVI) to increase the accuracy of urban extent extraction. However, limited by the spatial resolution of data sources, these methods are applied at regional or global scales in most cases.
In June 2018, Luojia 1-01 (LJ1-01), a new generation of NTL remote sensing satellite, was launched in China successfully. It has provided dramatic improvements over its predecessors, in terms of increased spatial resolution (~130 m), high radiometric quantization (14 bits) with a swath of 250 km. This new generation of NTL data source brings new possibilities and insights to the study of urban extent extraction, especially the improvement of image spatial resolution, which may not only increase the accuracy of urban extent mapping but also show more details of urban structure. Furthermore, urban extent extraction using LJ1-01 NTL data can to some extent fill the gaps in the field of NTL remote sensing in local-scale research.
The main purpose of this paper was to assess the accuracy of urban extent extraction at the city scale using LJ1-01 nighttime light imagery. For comparison, index-based, thresholding-based and classification-based methods were all adopted in the experiment. Besides, VIIRS DNB imagery were also utilized to reveal the difference between two kinds of NTL data. Finally, the Google Map high-resolution satellite images were used as reference and a visual interpretation method was applied for accuracy assessment.

2. Study Area and Data

2.1. Study Area

Wuhan, situated in the central part of China, is the capital city of Hubei Province. Wuhan covers an area of 8494.41 km2 and had a population of 10.89 million by the end of 2017 [59]. Due to the confluence of the Yangtze River and its greatest branch—the Han River, Wuhan is divided into three districts, namely Wuchang, Hankou and Hanyang. Wuhan has been known as the “City of Hundred Lakes”, where urban water area accounts for 1/4 of the total area. It has also been regarded as the “Thoroughfare to Nine Provinces”, serving as the transportation junction which links the east and west, as well as the north and the south [60]. Since the implementation of the reform and opening-up policy in 1978, the whole Chinese society has undergone dramatic socio-economic development, which accommodates urban expansion. Wuhan has also experienced an impressive social and economic transformation, such as industrialization and rapid urbanization and has become one of the fastest-growing cities in China during the past three decades.

2.2. Data

2.2.1. LJ1-01

LJ1-01, launched on 2 June 2018, is the first satellite in the Luojia-1 scientific experimental satellite series owned by Wuhan University. It is the first remote sensing satellite focusing on nighttime light in China and also the first low-orbit satellite with earth observation and satellite navigation enhancement functions. LJ1-01 is a 20 kg-level micro-nano satellite, equipped with a high sensitivity night-light camera that has a spectral bandwidth of 0.319 μm. It can obtain high-precision nighttime light imagery with a dynamic range of up to 14 bits at night, with a spatial resolution of 130 meters and a swath of 250 km. It will provide an objective basis for the research of socio-economic parameter estimation, eco-environmental disaster monitoring, major event assessment, public health and other fields, carrying out dynamic monitoring of macroeconomic operations in China and the world. The data is free to download at the High-Resolution Earth Observation System of the Hubei Data and Application Center (http://59.175.109.173:8888/, accessed in June 2018). Since the radiometric calibration for LJ1-01 imagery is still under improvement, we used digital number (DN) value in this experiment for analysis. The comparison of parameters for DMSP/OLS, NPP/VIIRS and LJ1-01 NTL data are shown in Table 1.

2.2.2. VIIRS DNB

VIIRS DNB Cloud-Free Composites (version 1) are a suite of average radiance composite images, where the data impacted by stray light, lightning, lunar illumination and cloud-cover have been excluded. The products are produced in 15 arc-second geographic grids, spanning the globe from 75 N latitude to 65 S. In this research, the VIIRS DNB data was utilized for comparison with LJ1-01 NTL data in urban extent extraction. Data were downloaded from the NOAA National Geophysical Data Center (https://ngdc.noaa.gov/eog/viirs/download_dnb_composites.html). Since the monthly composite images of June 2018 have not been released yet, the images of May 2018 were selected for the experiment.

2.2.3. Landsat 8 OLI

Landsat 8 Operational Land Imager (OLI) Level-2 data for the study area were obtained from the United States Geological Survey (USGS) website (https://earthexplorer.usgs.gov/). Level-2 data products, also called surface reflectance products, provide an estimate of the surface spectral reflectance at a 30-m spatial resolution in the absence of atmospheric scattering or absorption [61]. Since the Landsat image in June 2018 was partly covered by clouds, the closest cloud-free image in April 2018 was selected for this research. The data was utilized in both HSI and SVM methods to compare their performance in urban extent extraction.

2.2.4. Other Auxiliary Data

Google Map remote sensing satellite images, with a spatial resolution of 1 m, were utilized as reference data in the visual interpretation accuracy assessment. The study area under different data sources are shown in Figure 1.

3. Methods

3.1. Data Preprocessing

First of all, the data in the study area—Wuhan city, was clipped from the original datasets. Then, for consistency among different data sources, the LJ1-01 NTL, the VIIRS DNB, the Landsat 8 OLI products and the Google Map satellite images were all re-projected to Albers Conical Equal Area projection. After that, a bilinear interpolation algorithm was applied to resample the LJ1-01 NTL and VIIRS DNB images to the same resolution as that of the Landsat data (i.e., 30 m) and each image has a size of 3000 × 3000 pixels. In addition, since land-reflected moon light has not been excluded from the monthly composite products of VIIRS DNB data, it is necessary to remove background noise before use. In this research, all pixel values in VIIRS imagery were subtracted by 0.5 first, then pixels with negative DN values were assigned values of zero.
To apply the HSI method developed by Lu et al. [56], the NDVI image was generated first, as Equation (1) shows. However, this method has been proven to be difficult to separate urban areas and water bodies accurately [62]. Thus, to enhance the performance of the HSI, the NDWI image was also generated as a mask to remove water areas from the original HSI result. The NDWI proposed by Mcfeeters [63] was derived as Equation (2):
  N D V I = N I R R N I R + R  
  N D W I = G N I R G + N I R  
where N I R , R and G are Landsat 8 OLI surface reflectance images in near-infrared, red and green band respectively.

3.2. Urban Extent Extraction Based on Different Methods

3.2.1. Human Settlement Index

The Human Settlement Index, proposed by Lu et al. [56], is a method of extracting residential sites from the integration of nighttime light data and the NDVI index. In theory, NDVI images generated by Equation (1) have values ranging from −1 to 1. However, as the LJ1-01 NTL data has a dynamic radiometric range of 14 bits, the DN values in the study area ranging from 0 to 1879. Thus, it is necessary to match the NTL data with the range of NDVI index.
First, pixels with zero values were masked out and a natural logarithmic transformation was adopted to the LJ1-01 NTL images to maintain the rich details in the nighttime light. Then, the data was normalized as Equation (3):
  N T L n o r = N T L N T L m i n N T L m a x N T L m i n  
where N T L n o r is the normalized value of the LJ1-01 NTL image, N T L m a x and N T L m i n are the minimum and maximum values in the LJ1-01 NTL image respectively. Then the HSI was generated as Equation (4):
  H S I = ( 1 N D V I ) + N T L n o r ( 1 N T L n o r ) + N D V I + N T L n o r × N D V I  
After that, regions with DN values greater than 0 in the NDWI image were regarded as water bodies, the mask was generated and the corresponding regions in the HSI image were removed. Lastly, multiple thresholds (from 0.25 to 1.2 with an interval of 0.05) were applied to extract urban areas from the HSI image.
The processing of VIIRS DNB data was exactly the same as that of LJ1-01 data mentioned above. The logarithm transformation and normalization algorithm were adopted first, then the HSI was generated, with water bodies removed by NDWI. Lastly, multiple thresholds were applied to extract urban areas from the HSI image.

3.2.2. Simple Thresholding Segmentation

In order to compare with the result of the HSI method in Section 3.2.1, the Simple Thresholding Segmentation (STS) method was also adopted in this research using LJ1-01 NTL and VIIRS DNB data respectively [64]. Since it was difficult to determine the optimal threshold for urban extent extraction, multiple thresholds have been tested and finally, 1 to 20 with an interval of 1 have been selected as threshold values to extracted urban areas from the two kinds of NTL imagery.

3.2.3. SVM Supervised Classification

As mentioned in the Introduction section, thresholding-based, classification-based and index-based methods are the three most commonly used methods for urban extent mapping. Thus, in this research, the SVM supervised classification method was also adopted for comparison [65,66]. First, the normalized difference built-up index (NDBI) was generated. The NDBI index was proposed by Zha et al. [67] for urban built-up areas extraction using Landsat TM data. When applied to Landsat 8 OLI data, it was derived as [68]:
  N D B I = S W I R 1 N I R S W I R 1 + N I R  
where S W I R 1 and N I R are Landsat 8 OLI surface reflectance images in shortwave infrared 1 and near-infrared band respectively. Then, in addition to Landsat band 1~7, the generated NDBI image was also added as an additional band to enhance the information of built-up areas in Landsat data. The image was classified into five land-cover types presented in the study area, that is, bare soil, water bodies, vegetation, built-up land in high reflectance and built-up land in low reflectance. For each land-cover type, over 1000 pixels were selected as training samples based on interpretation of high-resolution images from Google Earth and separability for all the type-pairs were greater than 1.8. After performing the SVM classification, bare soil, water bodies and vegetation in the results were merged into non-urban areas, while the built-up land covers were recognized as urban areas.
In order to compare the results using different datasets, NTL data was also utilized in this experiment for classification-based urban extent extraction. The resampled LJ1-01 and VIIRS nighttime light imagery was added to Landsat band 1~7 as an additional band respectively and the Landsat-LJ and Landsat-VIIRS composite data were generated. Then the SVM classifier was applied using the same training samples as above. Lastly, the extraction results went for accuracy assessment, which will be presented in Section 3.3.

3.3. Accuracy Assessment

A visual interpretation method was utilized to quantify the performance of different urban extent extraction methods presented in Section 3.2. Google Map remote sensing satellite imagery was recognized as reference data, that is, “ground truth” data, with a spatial resolution of 1 m.
First, 800 randomly sampled points in the study area were generated. Then each point was labeled as “urban area” or “non-urban area” by visual interpretation using Google Map satellite images. After that, extracted urban areas obtained in Section 3.2 were recognized as “labeling results” and the confusion matrix was established for each method to perform the accuracy assessment.
The confusion matrix is a specific table that is often used to describe the accuracy of classification results [69]. Generally, rows in the confusion matrix correspond to classes, while columns correspond to predicted classes (or vice versa). Several metrics are often used as performance assessment measures based on confusion matrix, that is, Producer’s Accuracy (PA), User’s Accuracy (UA), Overall Accuracy (OA) and Kappa Coefficient (KC). PA measures the omission error, indicating the probability that an urban pixel is correctly identified. UA measures the commission error, indicating the probability that a labeled urban pixel is truly urban. OA indicates the proportion of pixels that are correctly identified and lastly, KC indicates the agreement between “labeling results” and “ground truth” by combining both types of errors, which can provide a more comprehensive assessment [70,71]. These metrics are computed as:
  P A k = x k k x k +  
  U A k = x k k x + k  
  O A = x k k N  
  K C = N x k k x k + x + k N 2 x k + x + k  
where x k k is the number of pixels that are correctly identified, x k + is the total number of pixels that belong to class k , x + k is the total number of pixels identified as class k and N is the total number of pixels in the dataset.

4. Results

Figure 2, Figure 3, Figure 4, Figure 5 and Figure 6 show some of the urban extent extraction results based on different methods. Figure 2 and Figure 3 show urban areas of Wuhan extracted by the HSI method using LJ1-01 and VIIRS data respectively, while Figure 4 and Figure 5 show urban areas extracted by STS method using LJ1-01 and VIIRS data respectively. And urban areas shown in Figure 6 were extracted by SVM method using Landsat only, Landsat-LJ composite and Landsat-VIIRS composite data. By comparing different results, it can be found that urban extent extracted by the HSI method using LJ1-01 and VIIRS data were generally similar. In the results derived from the STS method, a lot of spatial information within the city was lost, especially when using VIIRS data. While in contrast, the urban extent extracted by SVM method were overestimated in suburban areas especially when using Landsat data only.
The accuracy assessment for the performance of different methods are presented in Figure 7. It can be found that the urban extent extracted by the HSI method using LJ1-01 data with a threshold of 0.65 showed the largest KC value, 0.834, among all the results. Compared with other methods, the HSI method using LJ1-01 data had the best performance, with all KC values larger than 0.797 when the thresholds were between 0.5 and 0.75. In terms of the HSI method using VIIRS data, the KC values increased gradually from 0.2434 with the increase of the thresholds, reached the maximum value of 0.772 when the threshold was 0.95 and then began to slowly decline. In contrast, the KC value obtained by STS method showed different change trends from HSI method. When using LJ1-01 data, the KC values remained at a relatively higher level around 0.77 when the thresholds were between 3 and 12, then began to decrease with the increase of the thresholds. When using VIIRS data, the change trend was similar but the maximum KC value only reached 0.7512, which was lower than using LJ1-01 data. As for the SVM method, the KC values for using Landsat only, Landsat-LJ composite and Landsat-VIIRS composite data were 0.6427, 0.7829 and 0.7486 respectively.
In terms of PA, the HSI method using LJ1-01 and VIIRS data both maintained significantly high PA values when the thresholds were less than 0.7 and they also had similar change trends. In contrast, the PA values of the STS method using LJ1-01 and VIIRS data were relatively low and dropped sharply with increasing thresholds. In terms of UA, all the methods that used LJ1-01 data kept at a relatively high level and the SVM method using Landsat-LJ composite data had the best result. In the matter of OA, LJ1-01 data still had better performance than VIIRS data. The largest OA values for HSI, STS and SVM method when using LJ1-01 data were 0.9313, 0.915 and 0.9213 respectively, all larger than those of VIIRS data used ones.

5. Discussion

5.1. The Advantages of NTL Data in Urban Area Extraction

In this research, we compared several commonly used methods for extracting urban areas, including NTL based and non-NTL based methods. The experimental results demonstrated that the adding of NTL to Landsat data can improve the accuracy of urban extent extraction results.
Urban areas shown in Figure 6a were extracted by the SVM method using Landsat data only, while urban areas shown in Figure 6b,c were extracted by the same method, except that the data were replaced with Landsat-LJ composite and Landsat-VIIRS composite data respectively. As can be seen from the results, with the integration of NTL data, the overestimation in urban suburbs has been largely reduced and the number of small patches has also decreased. According to Figure 7, after adding VIIRS and LJ1-01 data, the KC values of SVM method increased from 0.6427 to 0.7486 and 0.7829 respectively, while the values of PA, UA and OA metrics all increased in varying degrees.
For a more intuitive comparison of NTL based and non-NTL based methods, urban extent extraction results with the largest KC values for HSI and STS method are presented in Figure 8a,b respectively, while Figure 8c presents the result of SVM method using Landsat data only, along with the Google Map satellite image in Figure 8d for comparison. It can be clearly seen that the HSI method (LJ1-01 data used, threshold = 0.65) had the best extraction result of urban extent, while the SVM method using only Landsat data had the worst. Specifically, the STS method (LJ1-01 data used, threshold = 5) shows a large overestimation in urban cores, while SVM method contrarily shows an underestimation in urban cores and overestimation in urban suburbs, resulting in a large number of redundant small patches. This is consistent with the accuracy assessment results in Figure 7. The values of OA and KC metrics both exhibited the same decreasing trends among these three results, that is, HIS > STS > SVM. This finding demonstrates that the SVM method using only Landsat data to extract urban extent may not be accurate enough. It is difficult to separate bare soils from urban land, which may cause overestimation in urban suburbs. The utilization of NTL data, which is able to represent luminous built-up areas at night, can greatly improve the accuracy of urban extent extraction.

5.2. The Advantages of LJ1-01 Data Compared with VIIRS

In this research, we also compared the performance of different NTL data in urban extent extraction. LJ1-01 NTL and VIIRS DNB images were utilized to contrast two main nighttime light data sources. Figure 9a,b displays the urban extent extraction results with the largest KC values for the HSI method using LJ1-01 and VIIRS data respectively and the results with the largest KC values for STS method are shown in Figure 9c,d. It can be seen that HSI methods using both LJ1-01 and VIIRS data were effective in urban extent extraction. In addition to human settlements, the luminous bridges and roads at night can also be identified accurately. However, the urban areas extracted based on VIIRS data were more fragmented with a large number of small patches and many urban pixels have not been identified. In contrast, the LJ1-01 based HSI method exhibited better performance, which extracted urban extent more accurately with less fragments and patches. This result is supported by the accuracy assessment according to Figure 7. The largest KC value for VIIRS based HSI method was 0.772, while it increased to 0.834 for the LJ1-01 based method. In terms of other metrics, the PA value increased from 0.871 to 0.9493, the UA value from 0.8043 to 0.824 and the OA value from 0.9075 to 0.9313 for VIIRS and LJ1-01 data. Moreover, it can be seen from Figure 9c,d that when applying the STS method, the extracted urban extent based on LJ1-01 data contained more urban spatial information than the VIIRS data. This finding can be obtained from the accuracy assessment in Figure 7 as well. The largest KC value for LJ1-01 based STS method, 0.79, was higher than that for the VIIRS based method. In terms of the SVM method, the PA, UA, OA and KC values using Landsat-LJ composite data were 0.7235, 0.9813, 0.9213 and 0.7829 respectively, all of which were larger than the corresponding values using Landsat-VIIRS composite data. The above results demonstrate that the LJ1-01 NTL data, which has a higher spatial resolution, can bring more abundant spatial details in urban extent extraction and lead to better identification results compared with VIIRS DNB data.

5.3. Prospects for the Future

According to the experiments and analysis above, we have concluded that LJ1-01 NTL data, with a spatial resolution of 130 m, contains more abundant spatial information of urban heterogeneity and can lead to a more accurate result when used in urban extent extraction. However, there are still some limitations in the application of LJ1-01 NTL data, which is worth more efforts in the future.
As a new nighttime light imagery source, the Luojia 1-01 satellite was launched in June 2018, which means that only images from June 2018 to now can be obtained for the time being. The lack of multi-temporal images makes LJ1-01 NTL data temporarily unable to be applied to long-term urban dynamic monitoring. However, if LJ1-01 NTL data can be integrated with other nighttime light data which has similar spatial resolution and long-time archived images, such as the International Space Station (ISS) imagery, this problem can be alleviated to a certain extent. The ISS imagery is a kind of multispectral photograph taken by astronauts in space using digital cameras, which has both daytime and nighttime archives [72]. All the photographs can be freely downloaded at “Gateway to Astronaut Photography of Earth” (https://eol.jsc.nasa.gov/), which is run by the National Aeronautics and Space Administration (NASA). Compared with other commonly used nighttime light data, the ISS imagery has several advantages, such as higher spatial resolution, three spectral bands in the visible range and variable overpass times [73]. It has an abundant data archive since last century and during the past decade, the quality of its imagery, especially nighttime light imagery has increased dramatically. However, for most global cities, its high-quality imagery archiving is not enough to do research on urban dynamic monitoring. With the release of LJ1-01 data, the integration of LJ1-01 and ISS nighttime light images may provide a solution to this problem. By integrating various types of high-resolution NTL data for a consistent observation, we can identify urban expansion areas between different periods in a fine scale, as well as better understand the dynamic patterns of urban land expansion. Moreover, further analysis in other geographical regions or at a global scale is to be performed. By combining high-resolution NTL data from different sources and other remote sensing data, an accurate extraction of urban extent in different periods and regions can be achieved. Therefore, a continuing mapping of urban dynamics at the global scale will also be possible in the future.

6. Conclusions

In this research, we compared different methods of urban extent extraction in Wuhan city, including a Human Settlement Index using both LJ1-01 and VIIRS images, Simple Thresholding Segmentation using both LJ1-01 and VIIRS images and SVM supervised classification using Landsat and NTL composite data. The performance of different methods was assessed by reference samples from visual interpretation of Google Map satellite images. The results demonstrate that the HSI method using LJ1-01 NTL data had the best performance in urban extent extraction, which can not only accurately extract human settlements but also identify luminous bridges and roads at night. In addition, it presented the largest Kappa Coefficient value, 0.834, among all the results when the threshold was 0.65, while other metrics remained at a high level. The relevant findings of this research are as follows:
(1) Compared with non-NTL methods, the addition of NTL data can largely improve the accuracy of the extraction results. For SVM classification method, the values of all accuracy assessment metrics have increased in varying degrees after adding LJ1-01 and VIIRS data respectively. By comparing the urban extent extraction results which have the largest KC values for different methods, we can find out that non-NTL based SVM classification method is likely to cause overestimation in urban suburbs, while NTL based methods exhibited better performance according to the accuracy assessment.
(2) LJ1-01 NTL data, which contains more spatial details, can lead to better identification results than VIIRS DNB data in urban extent extraction. The urban areas extracted by the VIIRS based HSI method were fragmented with a large number of small patches, while the LJ1-01 based method identified urban land more accurately. The same finding can be obtained from the accuracy assessment. The largest KC value for the VIIRS based method was 0.772, while it increased to 0.834 for the LJ1-01 based method. As for the STS and SVM methods, the LJ1-01 NTL data also presented better performance compared with VIIRS DNB data.
(3) LJ1-01, as a new generation of nighttime light imagery, has provided higher spatial resolution, wider radiometric measurement range and richer urban dynamic information over its predecessors. It displays a great potential for urban extent extraction at city scale and is bound to be more widely used in urban mapping. However, the current lack of multi-temporal images temporarily limits its application in urban dynamic monitoring. By integrating LJ1-01 with other high-resolution nighttime light imagery, we may find a solution to this problem. Moreover, the integration of high-resolution nighttime light imagery and other remote sensing data may provide a possibility to map urban extent in different periods and regions and then monitor urban dynamics continuously at the global scale in the future.

Author Contributions

X.L., H.X. and D.L. conceived and designed the experiments; L.Z. performed the experiments and wrote the paper; and all authors edited the paper.

Funding

This research was supported by Youth Fund of Wuhan Donghu University under grant no. 2018dhsk004, the Open Fund of Key Laboratory of the Ministry of Land and Resources for Law Evaluation Engineering under grant no. CUGFP-1804, National Natural Science Foundation of China under grant no. 41771386 and 41671518 and the Fundamental Research Funds for the Central Universities, China University of Geosciences (Wuhan).

Acknowledgments

We give thanks to the editors and anonymous reviewers for their valuable comments to improve our manuscript.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Wu, J. Urban ecology and sustainability: The state-of-the-science and future directions. Landsc. Urban Plan. 2014, 125, 209–221. [Google Scholar] [CrossRef]
  2. Seto, K.C.; Michail, F.; Burak, G.; Reilly, M.K. A meta-analysis of global urban land expansion. PLoS ONE 2011, 6, e23777. [Google Scholar] [CrossRef] [PubMed]
  3. Grimm, N.B.; Faeth, S.H.; Golubiewski, N.E.; Redman, C.L.; Wu, J.; Bai, X.; Briggs, J.M. Global change and the ecology of cities. Science 2008, 319, 756–760. [Google Scholar] [CrossRef] [PubMed]
  4. Statistics, C.N.B.O. Statistical Communiqué of the People’s Republic of China on the 2017 National Economic and Social Development. Available online: http://www.stats.gov.cn/tjsj/zxfb/201802/t20180228_1585631.html (accessed on 21 August 2018).
  5. Yuan, F.; Bauer, M.E. Comparison of impervious surface area and normalized difference vegetation index as indicators of surface urban heat island effects in landsat imagery. Remote Sens. Environ. 2007, 106, 375–386. [Google Scholar] [CrossRef]
  6. Superczynski, S.D.; Christopher, S.A. Exploring land use and land cover effects on air quality in central alabama using gis and remote sensing. Remote Sens. 2011, 3, 2552–2567. [Google Scholar] [CrossRef]
  7. Shahbaz, M.; Lean, H.H. Does financial development increase energy consumption? The role of industrialization and urbanization in tunisia. Energy Policy 2012, 40, 473–479. [Google Scholar] [CrossRef] [Green Version]
  8. Salazar, A.; Baldi, G.; Hirota, M.; Syktus, J.; Mcalpine, C. Land use and land cover change impacts on the regional climate of non-amazonian south america: A review. Glob. Planet. Chang. 2015, 128, 103–119. [Google Scholar] [CrossRef]
  9. Lu, Y.; Coops, N.C.; Hermosilla, T. Regional assessment of pan-pacific urban environments over 25 years using annual gap free landsat data. Int. J. App. Earth Obs. Geoinf. 2016, 50, 198–210. [Google Scholar] [CrossRef]
  10. Kalnay, E.; Cai, M. Impact of urbanization and land-use change on climate. Nature 2003, 423, 528–531. [Google Scholar] [CrossRef] [PubMed]
  11. Foley, J.A.; Defries, R.; Asner, G.P.; Barford, C.; Bonan, G.; Carpenter, S.R.; Chapin, F.S.; Coe, M.T.; Daily, G.C.; Gibbs, H.K. Global consequences of land use. Science 2005, 309, 570–574. [Google Scholar] [CrossRef] [PubMed]
  12. Zhang, L.; Peng, J.; Liu, Y.; Wu, J. Coupling ecosystem services supply and human ecological demand to identify landscape ecological security pattern: A case study in beijing–tianjin–hebei region, china. Urban Ecosyst. 2017, 20, 1–14. [Google Scholar] [CrossRef]
  13. Liu, X.; Liang, X.; Li, X.; Xu, X.; Ou, J.; Chen, Y.; Li, S.; Wang, S.; Pei, F. A future land use simulation model (flus) for simulating multiple land use scenarios by coupling human and natural effects. Landsc. Urban Plan. 2017, 168, 94–116. [Google Scholar] [CrossRef]
  14. Inkoom, J.N.; Nyarko, B.K.; Antwi, K.B. Explicit modeling of spatial growth patterns in shama, ghana: An agent-based approach. J. Geovis. Spat. Anal. 2017, 1, 7. [Google Scholar] [CrossRef]
  15. Chen, Y.; Li, X.; Liu, X.; Ai, B.; Li, S. Capturing the varying effects of driving forces over time for the simulation of urban growth by using survival analysis and cellular automata. Landsc. Urban Plan. 2016, 152, 59–71. [Google Scholar] [CrossRef]
  16. Chen, Y.; Li, X.; Liu, X.; Ai, B. Modeling urban land-use dynamics in a fast developing city using the modified logistic cellular automaton with a patch-based simulation strategy. Int. J. Geogr. Inf. Sci. 2014, 28, 234–255. [Google Scholar] [CrossRef]
  17. Yu, B.L.; Tang, M.; Wu, Q.S.; Yang, C.S.; Deng, S.Q.; Shi, K.F.; Peng, C.; Wu, J.P.; Chen, Z.Q. Urban built-up area extraction from log-transformed npp-viirs nighttime light composite data. IEEE Geosci. Remote Sens. Lett. 2018, 15, 1279–1283. [Google Scholar] [CrossRef]
  18. Xue, X.Y.; Yu, Z.L.; Zhu, S.C.; Zheng, Q.M.; Weston, M.; Wang, K.; Gan, M.Y.; Xu, H.W. Delineating urban boundaries using landsat 8 multispectral data and viirs nighttime light data. Remote Sens. 2018, 10, 799. [Google Scholar] [CrossRef]
  19. Sharma, R.C.; Tateishi, R.; Hara, K.; Gharechelou, S.; Iizuka, K. Global mapping of urban built-up areas of year 2014 by combining modis multispectral data with viirs nighttime light data. Int. J. Digit. Earth 2016, 9, 1004–1020. [Google Scholar] [CrossRef]
  20. Guo, W.; Li, G.Y.; Ni, W.J.; Zhang, Y.H.; Lu, D.S. Exploring improvement of impervious surface estimation at national scale through integration of nighttime light and proba-v data. Gisci. Remote Sens. 2018, 55, 699–717. [Google Scholar] [CrossRef]
  21. Kotarba, A.Z.; Aleksandrowicz, S. Impervious surface detection with nighttime photography from the international space station. Remote Sens. Environ. 2016, 176, 295–307. [Google Scholar] [CrossRef]
  22. Li, D.; Li, X. An overview on data mining of nighttime light remote sensing. Acta Geodaetica et Cartographica Sinica 2015, 44, 591–601. [Google Scholar]
  23. Chen, X.; Nordhaus, W.D. Using luminosity data as a proxy for economic statistics. Proc. Natl. Acad. Sci. USA 2011, 108, 8589–8594. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Henderson, J.V.; Storeygard, A.; Weil, D.N. Measuring economic growth from outer space. Am. Econ. Rev. 2012, 102, 994–1028. [Google Scholar] [CrossRef] [PubMed]
  25. Li, X.; Chen, X.; Zhao, Y.; Xu, J.; Chen, F.; Li, H. Automatic intercalibration of night-time light imagery using robust regression. Remote Sens. Lett. 2013, 4, 45–54. [Google Scholar] [CrossRef]
  26. 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]
  27. Xu, H.; Yang, H.; Li, X.; Jin, H.; Li, D. Multi-scale measurement of regional inequality in mainland china during 2005–2010 using dmsp/ols night light imagery and population density grid data. Sustainability 2015, 7, 13469–13499. [Google Scholar] [CrossRef]
  28. Ma, T. Quantitative responses of satellite-derived nighttime lighting signals to anthropogenic land-use and land-cover changes across china. Remote Sens. 2018, 10, 1447. [Google Scholar] [CrossRef]
  29. Yu, B.; Shu, S.; Liu, H.; Song, W.; Wu, J.; Wang, L.; Chen, Z. Object-based spatial cluster analysis of urban landscape pattern using nighttime light satellite images: A case study of china. Int. J. Geogr. Inf. Sci. 2014, 28, 2328–2355. [Google Scholar] [CrossRef]
  30. Zhang, Q.; Seto, K. 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]
  31. Letu, H.; Hara, M.; Yagi, H.; Naoki, K.; Tana, G.; Nishio, F.; Shuhei, O. Estimating energy consumption from night-time dmps/ols imagery after correcting for saturation effects. Int. J. Remote Sens. 2010, 31, 4443–4458. [Google Scholar] [CrossRef]
  32. Li, X.; Liu, S.; Jendryke, M.; Li, D.; Wu, C. Night-time light dynamics during the iraqi civil war. Remote Sens. 2018, 10, 858. [Google Scholar] [CrossRef]
  33. Li, X.; Li, D. Can night-time light images play a role in evaluating the syrian crisis? Int. J. Remote Sens. 2014, 35, 6648–6661. [Google Scholar] [CrossRef]
  34. Li, X.; Li, D.; Xu, H.; Wu, C. Intercalibration between dmsp/ols and viirs night-time light images to evaluate city light dynamics of syria’s major human settlement during syrian civil war. Int. J. Remote Sens. 2017, 38, 1–18. [Google Scholar] [CrossRef]
  35. Li, X.C.; Zhou, Y.Y. Urban mapping using dmsp/ols stable night-time light: A review. Int. J. Remote Sens. 2017, 38, 6030–6046. [Google Scholar] [CrossRef]
  36. Li, X.; Elvidge, C.; Zhou, Y.Y.; Cao, C.Y.; Warner, T. Remote sensing of night-time light. Int. J. Remote Sens. 2017, 38, 5855–5859. [Google Scholar] [CrossRef]
  37. Sutton, P.C. A scale-adjusted measure of “urban sprawl” using nighttime satellite imagery. Remote Sens. Environ. 2003, 86, 353–369. [Google Scholar] [CrossRef]
  38. Letu, H.; Hara, M.; Tana, G.; Nishio, F. A saturated light correction method for dmsp/ols nighttime satellite imagery. IEEE Trans. Geosci. Remote Sens. 2012, 50, 389–396. [Google Scholar] [CrossRef]
  39. Elvidge, C.D.; Keith, D.M.; Tuttle, B.T.; Baugh, K.E. Spectral identification of lighting type and character. Sensors 2010, 10, 3961–3988. [Google Scholar] [CrossRef] [PubMed]
  40. Small, C.; Elvidge, C.D.; Baugh, K. Mapping urban structure and spatial connectivity with viirs and ols night light imagery. Urban Remote Sens. Event 2013, 230–233. [Google Scholar] [CrossRef]
  41. 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]
  42. Elvidge, C.D.; Baugh, K.E.; Zhizhin, M.; Hsu, F.-C. Why viirs data are superior to dmsp for mapping nighttime lights. Proc. Asia-Pac. Adv. Netw. 2013, 35, 62–69. [Google Scholar] [CrossRef]
  43. Cao, C.; Xiong, X.; Blonski, S.; Liu, Q.; Guenther, B.; Weng, F. Suomi npp viirs on-orbit performance, data quality, and new applications. SPIE Asia-Pac. Remote Sens. 2012, 8528, 85280D. [Google Scholar]
  44. Dou, Y.Y.; Liu, Z.F.; He, C.Y.; Yue, H.B. Urban land extraction using viirs nighttime light data: An evaluation of three popular methods. Remote Sens. 2017, 9, 175. [Google Scholar] [CrossRef]
  45. Ma, W.T.; Li, P.J. An object similarity-based thresholding method for urban area mapping from visible infrared imaging radiometer suite day/night band (viirs dnb) data. Remote Sens. 2018, 10, 263. [Google Scholar] [CrossRef]
  46. Zhang, X.Y.; Li, P.J. A temperature and vegetation adjusted ntl urban index for urban area mapping and analysis. ISPRS J. Photogramm. Remote Sens. 2018, 135, 93–111. [Google Scholar] [CrossRef]
  47. Tong, L.; Hu, S.; Frazier, A.E. Mixed accuracy of nighttime lights (ntl)-based urban land identification using thresholds: Evidence from a hierarchical analysis in wuhan metropolis, china. Appl. Geogr. 2018, 98, 201–214. [Google Scholar] [CrossRef]
  48. Li, K.N.; Chen, Y.H. A genetic algorithm-based urban cluster automatic threshold method by combining viirs dnb, ndvi, and ndbi to monitor urbanization. Remote Sens. 2018, 10, 277. [Google Scholar] [CrossRef]
  49. Shi, K.; Huang, C.; Yu, B.; Yin, B.; Huang, Y.; Wu, J. Evaluation of npp-viirs night-time light composite data for extracting built-up urban areas. Remote Sens. Lett. 2014, 5, 358–366. [Google Scholar] [CrossRef]
  50. 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]
  51. Jing, W.; Yang, Y.; Yue, X.; Zhao, X. Mapping urban areas with integration of dmsp/ols nighttime light and modis data using machine learning techniques. Remote Sens. 2015, 7, 12419–12439. [Google Scholar] [CrossRef]
  52. Xu, M.; He, C.; Liu, Z.; Dou, Y. How did urban land expand in china between 1992 and 2015? A multi-scale landscape analysis. PLoS ONE 2016, 11, e0154839. [Google Scholar] [CrossRef] [PubMed]
  53. He, C.; Liu, Z.; Tian, J.; Ma, Q. Urban expansion dynamics and natural habitat loss in china: A multiscale landscape perspective. Glob. Chang. Biol. 2014, 20, 2886–2902. [Google Scholar] [CrossRef] [PubMed]
  54. Zhang, Q.; Wang, P.; Chen, H.; Huang, Q.L.; Jiang, H.B.; Zhang, Z.J.; Zhang, Y.M.; Luo, X.; Sun, S.J. A novel method for urban area extraction from viirs dnb and modis ndvi data: A case study of chinese cities. Int. J. Remote Sens. 2017, 38, 6094–6109. [Google Scholar] [CrossRef]
  55. Wang, R.; Wan, B.; Guo, Q.H.; Hu, M.S.; Zhou, S.P. Mapping regional urban extent using npp-viirs dnb and modis ndvi data. Remote Sens. 2017, 9, 862. [Google Scholar] [CrossRef]
  56. 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]
  57. Zhang, Q.L.; 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]
  58. Liu, X.; Hu, G.; Ai, B.; Li, X.; Shi, Q. A normalized urban areas composite index (nuaci) based on combination of dmsp-ols and modis for mapping impervious surface area. Remote Sens. 2015, 7, 17168–17189. [Google Scholar] [CrossRef]
  59. Bureau, W.S. Statistical Communiqué for Wuhan’s National Economic and Social Development in 2017. Available online: http://www.whtj.gov.cn/details.aspx?id=3957 (accessed on 22 August 2018).
  60. Government, W.M. Wuhan’s Advantages. Available online: http://english.wh.gov.cn/whgk_3581/whys/201204/t20120419_125685.html (accessed on 22 August 2018).
  61. Survey, U.S.G. Landsat 8 Oli/Tirs Level-2 Data Products–Surface Reflectance. Available online: https://lta.cr.usgs.gov/L8Level2SR (accessed on 8 September 2018).
  62. Ma, T.; Xu, T.; Huang, L.; Zhou, A. A human settlement composite index (hsci) derived from nighttime luminosity associated with imperviousness and vegetation indexes. Remote Sens. 2018, 10, 455. [Google Scholar] [CrossRef]
  63. 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]
  64. Otsu, N. A threshold selection method from gray-level histograms. IEEE Trans. Syst. Man. Cybern. 2007, 9, 62–66. [Google Scholar] [CrossRef]
  65. Cristianini, N.; An, S.T.J. Introduction to Support Vector Machines; China Machine Press: Beijing, China, 2005; pp. 658–664. [Google Scholar]
  66. Chapelle, O.; Vapnik, V.; Bousquet, O.; Mukherjee, S. Choosing multiple parameters for support vector machines. Mach. Learn. 2002, 46, 131–159. [Google Scholar] [CrossRef]
  67. Zha, Y.; Gao, J.; Ni, S. Use of normalized difference built-up index in automatically mapping urban areas from tm imagery. Int. J. Remote Sens. 2003, 24, 583–594. [Google Scholar] [CrossRef]
  68. Macarof, P.; Statescu, F. Comparasion of ndbi and ndvi as indicators of surface urban heat island effect in landsat 8 imagery: A case study of iasi. Present Environ. Sustain. Dev. 2017, 11, 141–150. [Google Scholar] [CrossRef]
  69. Congalton, R.G. A review of assessing the accuracy of classification of remotely sensed data. Remote Sens. Environ. 1991, 37, 270–279. [Google Scholar] [CrossRef]
  70. Lewis, H.G.; Brown, M. A generalized confusion matrix for assessing area estimates from remotely sensed data. Int. J. Remote Sens. 2001, 22, 3223–3235. [Google Scholar] [CrossRef] [Green Version]
  71. Foody, G.M. Status of land cover classification accuracy assessment. Remote Sens. Environ. 2002, 80, 185–201. [Google Scholar] [CrossRef]
  72. Kyba, C.; Garz, S.; Kuechly, H.; De Miguel, A.; Zamorano, J.; Fischer, J.; Hölker, F. High-resolution imagery of earth at night: New sources, opportunities and challenges. Remote Sens. 2014, 7, 1–23. [Google Scholar] [CrossRef] [Green Version]
  73. Miguel, A.S.D.; Castaño, J.G.; Zamorano, J.; Pascual, S.; Ángeles, M.; Cayuela, L.; Martinez, G.M.; Challupner, P.; Kyba, C.C.M. Atlas of astronaut photos of earth at night. Astron. Geophys. 2014, 55, 4–36. [Google Scholar] [CrossRef]
Figure 1. Study area under different data sources: (a) LJ1-01 NTL image in June 2018; (b) VIIRS DNB image in May 2018; (c) Landsat 8 OLI false color composite (R/G/B = 5/4/3) in April 2018; (d) Google Map remote sensing satellite image in 2018.
Figure 1. Study area under different data sources: (a) LJ1-01 NTL image in June 2018; (b) VIIRS DNB image in May 2018; (c) Landsat 8 OLI false color composite (R/G/B = 5/4/3) in April 2018; (d) Google Map remote sensing satellite image in 2018.
Sensors 18 03665 g001
Figure 2. Extracted urban areas of Wuhan by HSI method using LJ1-01 NTL data: (a) threshold = 0.4; (b) threshold = 0.8; (c) threshold = 1.2.
Figure 2. Extracted urban areas of Wuhan by HSI method using LJ1-01 NTL data: (a) threshold = 0.4; (b) threshold = 0.8; (c) threshold = 1.2.
Sensors 18 03665 g002
Figure 3. Extracted urban areas of Wuhan by HSI method using VIIRS DNB data: (a) threshold = 0.6; (b) threshold = 0.9; (c) threshold = 1.2.
Figure 3. Extracted urban areas of Wuhan by HSI method using VIIRS DNB data: (a) threshold = 0.6; (b) threshold = 0.9; (c) threshold = 1.2.
Sensors 18 03665 g003
Figure 4. Extracted urban areas of Wuhan by STS method using LJ1-01 NTL data: (a) threshold = 2; (b) threshold = 10; (c) threshold = 18.
Figure 4. Extracted urban areas of Wuhan by STS method using LJ1-01 NTL data: (a) threshold = 2; (b) threshold = 10; (c) threshold = 18.
Sensors 18 03665 g004
Figure 5. Extracted urban areas of Wuhan by STS method using VIIRS DNB data: (a) threshold = 2; (b) threshold = 10; (c) threshold = 18.
Figure 5. Extracted urban areas of Wuhan by STS method using VIIRS DNB data: (a) threshold = 2; (b) threshold = 10; (c) threshold = 18.
Sensors 18 03665 g005
Figure 6. Extracted urban areas of Wuhan by SVM method using: (a) Landsat data only; (b) Landsat and LJ1-01 composite data; (c) Landsat and VIIRS composite data.
Figure 6. Extracted urban areas of Wuhan by SVM method using: (a) Landsat data only; (b) Landsat and LJ1-01 composite data; (c) Landsat and VIIRS composite data.
Sensors 18 03665 g006
Figure 7. Accuracy assessment of urban extent extraction for: (a) HSI method using LJ1-01 data; (b) HSI method using VIIRS data; (c) STS method using LJ1-01 data; (d) STS method using VIIRS data; (e) SVM method.
Figure 7. Accuracy assessment of urban extent extraction for: (a) HSI method using LJ1-01 data; (b) HSI method using VIIRS data; (c) STS method using LJ1-01 data; (d) STS method using VIIRS data; (e) SVM method.
Sensors 18 03665 g007aSensors 18 03665 g007bSensors 18 03665 g007c
Figure 8. Urban extent extraction results with the largest Kappa Coefficient for: (a) HSI method (LJ1-01 data used and threshold = 0.65); (b) STS method (LJ1-01 data used and threshold = 5); (c) SVM method using Landsat data only; (d) Google Map satellite image for comparison.
Figure 8. Urban extent extraction results with the largest Kappa Coefficient for: (a) HSI method (LJ1-01 data used and threshold = 0.65); (b) STS method (LJ1-01 data used and threshold = 5); (c) SVM method using Landsat data only; (d) Google Map satellite image for comparison.
Sensors 18 03665 g008
Figure 9. Urban extent extraction results with the largest Kappa Coefficient for: (a) HSI method using LJ1-01 data (threshold = 0.65); (b) HSI method using VIIRS data (threshold = 0.95); (c) STS method using LJ1-01 data (threshold = 5); (d) STS method using VIIRS data (threshold = 9).
Figure 9. Urban extent extraction results with the largest Kappa Coefficient for: (a) HSI method using LJ1-01 data (threshold = 0.65); (b) HSI method using VIIRS data (threshold = 0.95); (c) STS method using LJ1-01 data (threshold = 5); (d) STS method using VIIRS data (threshold = 9).
Sensors 18 03665 g009
Table 1. Comparison of parameters for different NTL data.
Table 1. Comparison of parameters for different NTL data.
ParametersDMSP/OLSNPP/VIIRSLJ1-01
Available Period1992–2013November 2011–presentJune 2018–present
CountryThe U.S.The U.S.China
Spatial Resolution2.7 km740 m130 m
Swath3000 km3000 km250 km
Spectrum Range0.5–0.9 μm0.5–0.9 μm0.46–0.98 μm
Radiometric Resolution6 bits14 bits14 bits
SaturationSaturatedNot saturatedNot saturated

Share and Cite

MDPI and ACS Style

Li, X.; Zhao, L.; Li, D.; Xu, H. Mapping Urban Extent Using Luojia 1-01 Nighttime Light Imagery. Sensors 2018, 18, 3665. https://doi.org/10.3390/s18113665

AMA Style

Li X, Zhao L, Li D, Xu H. Mapping Urban Extent Using Luojia 1-01 Nighttime Light Imagery. Sensors. 2018; 18(11):3665. https://doi.org/10.3390/s18113665

Chicago/Turabian Style

Li, Xi, Lixian Zhao, Deren Li, and Huimin Xu. 2018. "Mapping Urban Extent Using Luojia 1-01 Nighttime Light Imagery" Sensors 18, no. 11: 3665. https://doi.org/10.3390/s18113665

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop