Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Analysis of land cover evolution within the built-up areas of provincial capital cities in northeastern China based on nighttime light data and Landsat data

  • Zhiwei Xie ,

    Contributed equally to this work with: Zhiwei Xie, Yaohui Han

    Roles Conceptualization, Supervision, Validation, Writing – review & editing

    Affiliation School of Transportation Engineering, Shenyang Jianzhu University, Hunnan District, Shenyang, China

  • Yaohui Han ,

    Contributed equally to this work with: Zhiwei Xie, Yaohui Han

    Roles Methodology, Writing – original draft, Writing – review & editing

    Current address: Liaoning Water Conservancy and Hydropower Survey and Design Research Institute CO., LTD, Heping District, Shenyang, China

    Affiliation School of Transportation Engineering, Shenyang Jianzhu University, Hunnan District, Shenyang, China

  • Lishuang Sun ,

    Roles Conceptualization, Formal analysis, Supervision

    645695280@qq.com

    Affiliation School of Transportation Engineering, Shenyang Jianzhu University, Hunnan District, Shenyang, China

  • Jiwei Ping

    Roles Data curation, Validation

    Affiliation School of Transportation Engineering, Shenyang Jianzhu University, Hunnan District, Shenyang, China

Abstract

Mastering the evolution of urban land cover is important for urban management and planning. In this paper, a method for analyzing land cover evolution within urban built-up areas based on nighttime light data and Landsat data is proposed. The method solves the problem of inaccurate descriptions of urban built-up area boundaries from the use of single-source diurnal or nocturnal remote sensing data and was able to achieve an effective analysis of land cover evolution within built-up areas. Four main procedures are involved: (1) The neighborhood extremum method and maximum likelihood method are used to extract nighttime light data and the urban built-up area boundaries from the Landsat data, respectively; (2) multisource urban boundaries are obtained using boundary pixel fusion of the nighttime light data and Landsat urban built-up area boundaries; (3) the maximum likelihood method is used to classify Landsat data within multisource urban boundaries into land cover classes, such as impervious surface, vegetation and water, and to calculate landscape indexes, such as overall landscape trends, degree of fragmentation and degree of aggregation; (4) the changes in the multisource urban boundaries and landscape indexes were obtained using the abovementioned methods, which were supported by multitemporal nighttime light data and Landsat data, to model the urban land cover evolution. Using the cities of Shenyang, Changchun and Harbin in northeastern China as experimental areas, the multitemporal landscape index showed that the integration and aggregation of land cover in the urban areas had an increasing trend, the natural environment of Shenyang and Harbin was improving, while Changchun laid more emphasis on the construction of artificial facilities. At the same time, the method proposed in this paper to extract built-up areas from multi-source city data showed that the user accuracy, production accuracy, overall accuracy and Kappa coefficient are at least 3%, 1%, 1% and 0.04 higher than the single-source data method.

Introduction

Urbanization is a complex phenomenon that has a profound relationship with land cover evolution (LCE) and social dynamics [1, 2]. As the urbanization process has progressed, it has changed from continuous expansion to compression within compact areas. Compact urbanization mainly reflects land cover evolution within urban built-up areas (UBAs), which is manifested by changes in surface landscape patterns, such as impervious ground, vegetation and water [35]. LCE can directly reflect the trend of urban master planning in a specific period, which has important practical significance for urban planning and management.

LCE is a description of the changing characteristics and patterns of land surface morphology, not only to analyze the shifts between land cover classes, but also to identify the evolutionary patterns implicit [6]. Land cover change detection and evolution pattern analysis are the main research components of LCE [7]. Along with the development of satellite technology and remote sensing image processing technology, remote sensing has played an important role in the analysis of LCE [8]. The main methods for detecting land cover change are the post-classification comparison method and the direct comparison method [9]. Subject to current automatic image change detection techniques, post-classification comparison is still the dominant technical approach, and supervised classification methods are mostly used [10]. Land cover evolution models are used to describe, interpret, predict, optimize, and support decision making on land cover change. The land cover evolution model analysis model mainly includes quantitative analysis model, ecological quality analysis model, land use degree model, spatial transfer change analysis model and change prediction model [11, 12]. Most of the above models use a single or composite index to express land change indicators and pay less attention to the evolution of the spatial distribution pattern.

The urban landscape pattern has a strong correlation with the spatial distribution pattern of impervious ground, vegetation and water, which can be applied to the study of the urban heat island effect, green space optimization, the mechanism of water logging disaster formation, and other issues [13, 14]. For the study of compact urbanization, researchers have analyzed the trends of urban construction by analyzing changes in landscape patterns within built-up areas to quantitatively reflect LCE [15, 16]. The landscape index provides morphological information about highly condensed landscapes that reflects the structural composition and spatial configuration of the landscape elements. Fragstats software is one of the landscape index analysis platforms that is widely used by researchers [17, 18].

Through sprawling urbanization, the rate of urbanization in the post-World War II developing world has been remarkable. However, the sprawling urbanization has resulted in the decay of the original urban areas, and the central urban area at night has become an empty city [19]. This has resulted in a huge waste of land and energy resources, and has caused ecological damage around the city. As a result, there is a trend towards developing compact urbanization and an urgent need to assess the land cover evolution within built-up [20] areas. Using remote sensing image data to determine urban built-up area boundaries (UBAB) is an important process for analysis of LCE and compact urbanization. Depending on the number of data sources used, the methods for extracting UBA can be divided into single-source data and multisource data methods. The single-source data method generally uses nighttime light data or Landsat data separately. In recent years, due to the strong stability of nighttime light data and the high integrity of UBA extraction results, many studies have been conducted, such as those on the neighborhood extreme value method, statistical data reference method and gradient spatial constraint method [2123]. However, nighttime light data can only capture the spatial extent of urban areas with obvious light at night, which does not include factories, development zones and other areas at the edges of cities that are not lit at night [24, 25]. Landsat data are also widely used to extract UBA, but the complexity of spectral features has resulted in the misclassification of some land cover types [26]. In addition, the completeness of the extraction results is insufficient, and the gap areas formed by the vegetation and water need to be filled manually [27]. The combination of nighttime light data and Landsat data can improve the accuracy of UBA extraction, and researchers have proposed methods, such as the fusion variable method, weighted pixel averaging method, and feature fusion method [2830].

China is the second-largest economy in the world and is urbanizing at a remarkable speed. Northeastern China, which is an old industrial base, plays an important role in China's economic structure. Since the implementation of the policy of revitalizing the old industrial base in northeastern China in 2003, the urbanization rate in this area has been the highest in the country. Due to the rapid expansion of urban areas, the quality of urbanization in northeastern China is poor [31, 32]. To develop a compact urbanization, it is of great practical significance to ascertain the landscape pattern changes in the capital cities of northeastern China and to gain insight into the LCE [33].

In this paper, we propose a method for analyzing the land cover evolution within built-up areas that can obtain the tendency of urban land construction and allow a deeper understanding of compact urbanization. First, the neighborhood extreme method and maximum likelihood method are used to extract the UBAB. Second, the single-source UBAB data described above is integrated with the weighted boundary pixel fusion algorithm to build the multisource UBAB. Then, Landsat data within the UBAB are classified into classes, such as impervious ground, vegetation and water, by the maximum likelihood method, and the corresponding landscape pattern indexes are calculated. Finally, the multitemporal landscape pattern indexes are used to obtain the landscape pattern changes and to analyze the land cover evolution trends. To verify the effectiveness of the method, Shenyang, Changchun and Harbin, capital cities in northeastern China, were used as experimental areas and compared with other methods. At the end of this paper, the main advantages and disadvantages of this method are described. The technical flowchart of this paper is shown in Fig 1.

thumbnail
Fig 1. Analysis of land cover evolution within built-up areas based on nighttime light data and Landsat data.

https://doi.org/10.1371/journal.pone.0239371.g001

Materials and methods

Study areas

This study selected three capital cities in northeastern China as the research areas: Shenyang city, Liaoning Province; Changchun city, Jilin Province; and Harbin city, Heilongjiang Province. As China's first heavy industrial and agricultural bases, the capital cities of these three northeastern provinces had earlier economic development and once had the highest level of urbanization in China. In addition, these cities are the political, economic, cultural, educational and transportation centers of the region. Shenyang is the largest city in northeastern China, with a total area of 13,000 square kilometers and a resident population of 8,316,000 in 2019 [34]; Changchun is located in the geographical center of northeastern China, with a total area of 20,500 square kilometers and a resident population of 7,667,000 in 2019 [35]; and Harbin is the northernmost provincial capital city in China, with a total area of 53,100 square kilometers and a resident population of 9,515,000 in 2019 [36].

Data

Nighttime light data.

The nighttime light data used in this paper were the DMSP/OLS (Defense Meteorological Satellite Program/Operational Linescan System) published by NOAA. These data have the advantages of long time series and high stability and are widely used to extract urban built-up areas [3739]. The DMSP/OLS data range from 1992 to 2013 and have a spatial resolution of 1000 m. In this study, the DMSP/OLS stable nighttime light data from 2000 to 2012 were selected, and the effects of oversaturation and discontinuity were reduced by relative radiometric correction.

Landsat data.

Landsat imagery courtesy of NASA Goddard Space Flight Center and U.S.

Geological Survey, downloaded from USGS (https://earthexplorer.usgs.gov/). Landsat data are image data published by NASA and are widely used to analyze urbanization due to their high spatial resolution. In this paper, Landsat TM data from 2000 to 2012 were selected for the extraction of built-up areas; data with interference from the sensors or the atmosphere were excluded. In addition, the reference built-up areas that were consistent with the experimental time period were identified by artificial interpretation [40].

Methodology

Data preprocessing

Nighttime light data.

The DMSP/OLS data were set to a Lambert isometric cone projection with a resampling accuracy of 1 km. To address the grayscale saturation of DMSP/OLS, we selected the 2017 F16 data as the baseline data and applied the following one-dimensional quadratic regression model to the other data for saturation correction [39]: (1) where DN and DNb represent the grayscale values of the image before and after correction, respectively; p and q are regression coefficients, and m is a constant.

To address the problem of fluctuating gray values in multitemporal nighttime light data, formula (2) was used to perform a relative radiation correction for the image. (2) where DN(n−1,i) and DN(n,i) represent the gray values of pixel i in years n−1 and n, respectively.

Landsat data.

Landsat imagery courtesy of NASA Goddard Space Flight Center and U.S.

Geological Survey, downloaded from USGS (http://landsat.visibleearth.nasa.gov/). To eliminate the effects of atmospheric water vapor, oxygen and methane on the radiation from the features [41], the FLAASH model was used for atmospheric correction of the images. Subsequently, the multitemporal remote sensing images were calibrated using a polynomial algorithm.

Urban built-up area extraction based on single source data

Urban built-up area extraction based on the neighborhood extremum method using DMSP/OLS data.

The neighborhood extremum method uses the change in grayscale of the nighttime light data from the built-up to the non-built-up areas as the basis, using the locations of gray mutation point as boundary points. The closed curve formed by combining all boundary points is the UBAB. First, the linear filtering method is used to construct the neighborhood difference image. The gray values of the DMSP/OLS data gradually decrease from built-up to non-built-up areas and show a gradient of changes. To enhance the grayscale gradient variation, an 8-neighborhood extremum template was constructed. The template calculates the absolute difference in the gray values between the central image and the rest of the images; the maximum difference is used as the grayscale value of the central image. The template calculates the absolute difference in the gray values between the central pixel and the rest of the pixels, taking the maximum difference as the gray value of the central pixel. The entire image is filtered using the 8-neighborhood extreme template, and a neighborhood difference image reflecting the grayscale variation feature was obtained, as shown in Fig 2. The black and white areas in the figure correspond to areas with small and large changes in gray values, respectively.

thumbnail
Fig 2. Neighborhood difference image.

(a)Shenyang;(b) Changchun;(c) Harbin.

https://doi.org/10.1371/journal.pone.0239371.g002

Then, the UBABs were extracted using an extreme value search method. The neighborhood difference pixels at the junction of the built-up and non-built-up areas have the greatest variation in gray values, which can be expressed as follows: (3) where DN(i,j) is the target pixel and represents the gray value of the pixel in row i, column j. We extracted the pixels that satisfy the above formula, set the remaining pixels to 0, and constructed the boundary extreme image, as shown in Fig 3(A), 3(B) and 3(C). The gray values of the cross section of boundary extreme images are extracted, and they are shown in Fig 3(A), 3(E) and 3(F). The collection of pixels with gray values greater than or equal to the lowest peak point was used as the boundary, resulting in the boundary image. Finally, object-oriented classification was used to extract the built-up areas of the boundary image, and the UBAB extraction from the DMSP/OLS data was achieved.

thumbnail
Fig 3. Boundary extreme images and diagram of cross section line gray values.

(a) Boundary extreme images of Shenyang; (b) Boundary extreme images of Changchun; (c) Boundary extreme images of Harbin; (d) Schematic of cross section line gray values of Shenyang; (e) Schematic of cross section line gray values of Changchun; (f) Schematic of cross section line gray values of Harbin.

https://doi.org/10.1371/journal.pone.0239371.g003

Urban built-up area extraction based on the maximum likelihood method using landsat data.

Remote sensing image classification methods are divided into supervised classification and unsupervised classification. Unsupervised classification only depends on the distribution of pixel gray value in spectral feature space, which makes it difficult to achieve accurate matching between pixel clusters and ground classes. Because of its simple calculation and strong stability, maximum likelihood method is still widely used in supervised classification algorithm [42]. The maximum likelihood method is an image classification method that combines the distances of samples to the centers of known classes and spectral distribution features. This method is widely used to extract urban built-up area from Landsat data [43]. The image classes were impervious surface, vegetation and water, and 100 training samples were selected for each category through artificial interpretation. The mean vector and covariance matrix of the training samples were calculated to obtain the corresponding probability discriminant function and Bayesian discriminant rule. Finally, the image data were classified, and the classification results were postprocessed with clumping, sieving, removal of outliers, etc. The pixels belonging to the impervious surface class were used as built-up areas, while the water and vegetation (urban green spaces) located within the construction land were also used as built-up areas and the rest as non-built-up areas, thus enabling the UBAB extraction.

Urban built-up areas extraction based on the multisource boundary pixels fusion method

Advantages and disadvantages of urban built-up area extraction using single-source data.

The results of the urban built-up area boundary extraction with DMSP/OLS data and Landsat data are shown in Fig 4, and the advantages and disadvantages of the above results were assessed as follows:

  1. The extraction of the urban built-up area boundaries using DMSP/OLS data had the advantage of stability and high integrity. The light sensitivity of the DMSP/OLS data may be missed in areas with weak lighting, as shown in Fig 4(B). In addition, the light spillover nature of DMSP/OLS data has led to a wide range of UBABs.
  2. The high spatial resolution of the Landsat data makes it possible to better reflect details. Landsat data have difficulty obtaining high-quality time series due to clouds and sensor errors. The spectral complexity of the data also affects the completeness and correctness of the results, as shown in Fig 4(A).
thumbnail
Fig 4. Advantages and disadvantages of the urban built-up area boundary extraction based on single source data.

(a) Built-up area extraction results of Shenyang in 2000 based on single source data; (b) Schematic of omissions in the extraction results with Landsat data; (c) Schematic of omissions in the extraction results with the DMSP/OLS data. (USGS/NASA Landsat).

https://doi.org/10.1371/journal.pone.0239371.g004

Based on the above analysis, to solve the problem of omission detection in built-up areas extracted from single source data, this paper proposes a multisource boundary pixel fusion method. The method combines the advantages of multisource data to optimize UBAB.

Multisource boundary pixel fusion method. Landsat data and DMSP/OLS built-up area boundaries were used as reference boundaries and target boundaries, as shown in Fig 5.

thumbnail
Fig 5. Schematic of the multisource data boundary fusion method.

https://doi.org/10.1371/journal.pone.0239371.g005

First, the tangent L is constructed for each pixel O(x0,y0) in the reference boundary. Next, the normal distribution of point O is constructed; the intersection with the target boundary is points A and B, and the nearest intersection A(x1,y1) is chosen. Finally, the fusion of points O and A is calculated using Eq (4) to obtain the fusion point C(x,y). Multisource data fusion boundaries are formed by the fusion points as follows: (4) where a and b are the fusion weights, a+b = 1 and in general a = b = 0.5.

Analysis of land cover evolution within built-up areas based on landscape index

The urban built-up area boundaries were obtained using the multisource boundary pixel fusion method, and land cover classes, such as impervious surface, water and vegetation, were obtained by classifying the Landsat data within the boundaries by the maximum likelihood method. We calculated the landscape index to reflect land cover evolution. The principle of selecting landscape indexes is that they can represent the characteristics of ecological state and landscape pattern, and have low redundancy [44]. Landscape index can be divided into three categories: overall trend index, fragmentation index and aggregation index [45]. In combination with the research focus of compact urbanization in this paper, Class area (CA), Percent of landscape (PL), Patch density (PD), Aggregation index (AI) were selected from above categories respectively.

(1) Overall trend index

① Class area (CA): CA is the size of the land cover and reflects the overall level of development of the land cover within UBA; it can be expressed as follows: (5) where ai is the area of class i.

② Percent of landscape (PL): PL is a measure of the landscape composition that can be used to assess the tendency of land cover change within an UBA. The formula is as follows: (6) where aij is the area of patch j that belongs to class i. A denotes the total landscape area and 0<PLAND<1.

(2) Fragmentation index

Patch density (PD): PD is the class number of land cover in a unit area, judging the fragmentation degree of land cover within UBA. The formula is as follows: (7) where bn is the number of patch n, and A is the total area of all patches.

(3) Aggregation index

Aggregation index (AI): AI reflects the aggregation degree of the patches belonging to the same class and the spatial configuration of the land cover. The formula is as follows: (8) where gi is the number of adjacent patches belonging to class i and max→gi is the largest number of adjacent patches.

Results

Urban built-up area extraction

The DMSP/OLS data and Landsat data from 2000, 2004, 2008 and 2012 were used as experimental data. The extraction results of neighborhood extremum method and maximum likelihood method are shown in Figs 6 and 7.

thumbnail
Fig 6. Extraction results of the urban built-up areas by the neighborhood extremum method.

https://doi.org/10.1371/journal.pone.0239371.g006

thumbnail
Fig 7. Extraction results of the urban built-up areas by the maximum likelihood method.

(USGS/NASA Landsat).

https://doi.org/10.1371/journal.pone.0239371.g007

The UBABs of Shenyang, Changchun and Harbin were extracted with the multisource boundary pixel fusion method, and the results are shown in Fig 8 and Table 1. The results showed that the UBA in all three cities had shown a clear growth trend. Shenyang, Changchun and Harbin had the fastest growth in the urban built-up areas from 2008 to 2012, 2008 to 2012 and 2000 to 2004, respectively.

thumbnail
Fig 8. Extraction results of the urban built-up areas by the multisource boundary pixel fusion method.

(USGS/NASA Landsat).

https://doi.org/10.1371/journal.pone.0239371.g008

Analysis of land cover evolution within urban built-up areas

Land cover classification.

The Landsat data were clipped using the urban built-up area boundaries and classified into impervious surfaces, water and vegetation using the maximum likelihood method, and the classification results are shown in Fig 7 and Table 2. As shown in Fig 9, the land cover classes within the three UBAs showed trends of expansion, with a higher degree of expansion in impervious surfaces and vegetation. The average growth rates of the Shenyang land cover classes were above 20%, with the highest growth in vegetation from 2000 to 2004 and in water and impervious surfaces from 2004 to 2008. The average growth rate of the impervious surfaces in Changchun was approximately 50%, with the fastest growth period from 2000 to 2004. The growth rate of land cover in Harbin was approximately 30%, with the fastest growth in vegetation and impervious surfaces occurring from 2000 to 2004 and in water from 2004 to 2008.

thumbnail
Table 2. Statistics of land cover classification from 2000 to 2012 (km2).

https://doi.org/10.1371/journal.pone.0239371.t002

Land cover evolution analysis based on landscape pattern index.

The landscape indexes of the impervious surfaces, vegetation and water were calculated separately, including the CA, PL, PD and AI, and the statistical results are shown in Table 3. The specific analysis was as follows:

  1. Shenyang: Landscape index trend chart for Shenyang is shown in Fig 10. CA: all land cover classes showed an increasing trend, with the fastest growth in impervious surfaces, which reflected the expansion of the built-up areas; PL: the impervious surfaces covered the largest percentage of the built-up areas and showed a decreasing trend, while vegetation and water occupancy showed increasing trends; PD: all land cover classes showed downward trends, with the largest rate of decline in vegetation, which reflected reduced fragmentation and increased homogeneity of the land cover; AI: all land cover classes showed downward trends, with the highest rate of decline in vegetation, which reflected the increased agglomeration of land cover classes. In general, impervious surface indexes had the highest values, but the occupancy, aggregation and integrity of vegetation increased, and the ecological environment of Shenyang tended to be optimized.
  2. Changchun: Landscape index trend chart for Shenyang is shown in Fig 11. CA: all land cover classes showed increasing trends, with the highest rate of growth in impervious surfaces, which reflected the expansion of built-up areas. PL: the impervious surface showed an upward trend, while vegetation and water showed downward trends, reflecting the fact that the area occupied by impervious surfaces was still increasing; PD: all land cover classes showed a downward trend, with the highest rate of decline in vegetation, reflecting the reduced fragmentation of the land cover; AI: the impervious surface and vegetation show a downward trend and the water showed an upward trend, reflecting increasing aggregation in the former classes and decreasing aggregation in the latter. Overall, the impervious surface of Changchun is being further strengthened, with greater emphasis on the construction of artificial features.
  3. Harbin: Landscape index trend chart for Shenyang is shown in Fig 12. CA: all land cover classes showed growing trends, with impervious surfaces showing the highest growth rate, which reflected their greater contribution to urban expansion. PL: the proportion of impervious surface was the highest, but the proportion of vegetation was increasing, reflecting the improving natural environment of the city; PD: the fragmentation of all land cover classes showed a decreasing trend, with vegetation fragmentation decreasing the fastest, reflecting the increasing homogeneity of land cover; AI: the areas of all land cover classes increased, but the trend towards impervious surfaces was more pronounced. Overall, ecological construction in Harbin was increasing, while the land cover was in an aggregated state.
thumbnail
Table 3. Landscape index within urban build-up areas from 2000 to 2012.

https://doi.org/10.1371/journal.pone.0239371.t003

In summary, the spatial extent of the land cover classes in Shenyang, Changchun and Harbin all showed growing trends, with increases in both integration and aggregation. In addition, the areas occupied by vegetation in Shenyang and Harbin were increasing, indicating that the natural environment was improving, while Changchun was focused on the construction of artificial facilities.

Discussion

Urban built-up area extraction accuracy

To verify the validity of this paper's method, the built-up areas of Shenyang, Changchun and Harbin were extracted using the artificial interpretation method. The built-up areas of the DMSP/OLS and Landsat data were extracted with the neighborhood extreme value method and the maximum likelihood method, respectively, and compared with the multisource boundary pixel fusion method.

As seen from Table 4 the user accuracy, production accuracy, overall accuracy and kappa coefficients of the boundary image fusion method were 88%, 90%, 99% and 0.88, respectively, which were at least 3%, 1%, 1% and 0.04 better than the other two methods. The method in this paper had higher accuracy than the built-up area extraction method based on single source data.

thumbnail
Table 4. Confusion matrix for the extraction of built-up areas.

https://doi.org/10.1371/journal.pone.0239371.t004

Land cover classification accuracy

To assess the accuracy of the land cover classification, the classification results of water, vegetation and impervious surfaces were assessed for accuracy. The 100 validation samples were randomly generated for each class, and their classes were interpreted visually. The interpretation results were compared with the computer classification results, and the average accuracy is shown in Table 5. As shown in Table 5, the overall accuracy and Kappa coefficient of the method used in this paper were above 95% and 0.90, respectively, which meets the accuracy requirements.

thumbnail
Table 5. Confusion matrix for the land cover classification.

https://doi.org/10.1371/journal.pone.0239371.t005

Boundary pixel fusion weight analysis

Different combinations of weights for the multisource boundary pixel fusion algorithm can affect the accuracy of built-up area extraction. We set a and b to different values and calculated the Kappa coefficients of the extraction results, which are shown in Fig 13. As seen from Fig 13, setting the weights to a = b = 0.5 give the highest accuracy.

Conclusions

Through analyses of provincial capitals in northeastern China, the validity of the proposed method for analyzing land cover evolution within urban built-up areas was demonstrated. We extracted urban built-up areas using DMSP/OLS data and Landsat data and analyzed land cover evolution within built-up areas using landscape indexes. To improve the extraction accuracy of the built-up areas based on single-source remote sensing data, the multisource boundary pixel fusion method was proposed, which effectively integrates the information advantages of multisource data. Using the urban built-up areas as the spatial extent, the land cover evolutions of Shenyang, Changchun and Harbin from 2000 to 2012 were analyzed through the landscape indexes, and the following results were obtained: the built-up areas exhibited clear expansion trends, with increased integration and agglomeration of land cover classes, such as vegetation, water and impervious surfaces; the construction of the ecological environments in Shenyang and Harbin were being strengthened, while Changchun was still focused on the construction of artificial facilities. The pillar industry of provincial capitals in northeastern China is industry. Compact urbanization needs to transfer the industrial facilities in the city center to the periphery, so as to improve the land use efficiency in the main urban area. The paper had shown that the land use efficiency in the main urban areas is improving, and Shenyang and Harbin have advantages in ecological environment. It has been said that further improvement of land use efficiency is very beneficial to the economic development of northeastern China.

Although the method in this paper improves the accuracy of the extraction of urban built-up areas based on remote sensing data, there are still some shortcomings to be studied. First, the analysis of compact urbanization in this paper is limited to the land cover evolution within built-up areas, without considering the impact of other socio-economic factors, and the proposals for the development of compact cities are less well covered. Second, this paper adopted the DMSP/OLS data as a nighttime light data source; these data are subject to certain constraints due to the acquisition time and spatial resolution, and the subsequent consideration would be to select VIIRS/NPP data, which have a higher spatial resolution. Third, the spatial resolution of the Landsat data still poses a disadvantage, and the selection of GF-2, WordView-2 or other high-resolution remote sensing images could improve the extraction accuracy, but the amount of data operations would increase; reasonable classification methods need further study. In addition, this paper only used landscape indexes to describe the land cover evolution and did not predict evolution trends; in future research, integration with cellular automata (CA), Markov, and other predictive models will be the focus of our approach.

Acknowledgments

We gratefully thank the data distribution agencies who provided the publicly released data used in this study. The DMSP/OLS data and the Landsat data are obtained from the National Oceanic and Atmospheric Administration (NOAA) and National Aeronautics and Space Administration (NASA), respectively.

References

  1. 1. Lu Y. Analysis and simulation of land use/Land cover change in Heifei city spatiotemporal data fusion of remote sensing. Doctor, Anhui University, Hefei, China, 2019. (in Chinese)
  2. 2. He J, Shi X, Gong J. Modeling the spatial of urban agglomeration considering their spatial interaction: A case study of Wuhan metropolitan area. Geomatics and Information Science of Wuhan University, 2016, 41(06):462–467. https://doi.org/10.13203/j.whugis20140250
  3. 3. Kim J, Li W, Newman G, Kil S H, Park S Y. The influence of urban landscape spatial patterns on single-family housing prices. Environment and Planning B: Urban Analytics and City Science, 2018, 745(1): 26–43. https://doi.org/10.1177/0265813516663932
  4. 4. Feng Y, Liu Y, Tong X. Spatiotemporal variation of landscape patterns and their spatial determinants in Shanghai, China. Ecological Indicators, 2018, 87: 22–32. https://doi.org/10.1016/j.ecolind.2017.12.034
  5. 5. Xiao R, Wang G, Zhang Q, Zhang Z. Multi-scale analysis of relationship between landscape pattern and urban river water quality in different seasons. Scientific Reports, 2016, 6(1): 1–10. pmid:28442746
  6. 6. Espinoza-Molina D, Bahmanyar R, Datcu M, et al. Land-cover evolution class analysis in Image Time Series of Landsat and Sentinel-2 based on Latent Dirichlet Allocation. 2017 9th International Workshop on the Analysis of Multitemporal Remote Sensing Images (MultiTemp). IEEE, 2017, 1–4. https://doi.org/10.1109/Multi-Temp.2017.8035261
  7. 7. Wang C, Lei S, Ekmore A J, et al. Integrating temporal evolution with cellular automata for simulating land cover change. Remote Sensing, 2019, 11(3): 301. https://doi.org/10.3390/rs11030301
  8. 8. Quader M A, Agrawal S, Kervyn M. Multi-decadal land cover evolution in the Sundarban, the largest mangrove forest in the world. Ocean & Coastal Management, 2017, 139: 113–124. https://doi.org/10.1016/j.ocecoaman.2017.02.008
  9. 9. Lu D, Mausel P, Brondizio E, et al. Change detection techniques. International Journal of Remote Sensing, 2004, 25(12): 2365–2401. https://doi.org/10.1080/0143116031000139863
  10. 10. Zhang C, Yue P, Tapete D, et al. A deeply supervised image fusion network for change detection in high resolution bi-temporal remote sensing images. ISPRS Journal of Photogrammetry and Remote Sensing, 2020, 166: 183–200. https://doi.org/10.1016/j.isprsjprs.2020.06.003
  11. 11. Veldkamp A, Lambin E F. Predicting land-use change. Agriculture, Ecosystems & Environment, 2001, 85(1–3): 1–6. https://doi.org/10.1016/S0167-8809(01)00199-2
  12. 12. Tao L. Land cover change and its ecological environment response based on"3S" technique in Dongting Lake region. Doctor, Hunan Agricultural University, Guangzhou, China, 2017. (in Chinese)
  13. 13. Pan J, Li X, Feng Z. Analysis of spatial and temporal patterns of impervious surfaces and vegetation covers in Lanzhou base on the V-I-AP model. Resources Science, 2010, 32(03): 130–137. (in Chinese)
  14. 14. Nong D H, Lepczyk C A, Miura T, et al. Quantifying urban growth patterns in Hanoi using landscape expansion modes and time series spatial metrics. PloS one, 2018, 13(5): e0196940. pmid:29734346
  15. 15. Liu P, Jia S, Han R, et al. Landscape pattern and ecological security assessment and prediction using remote sensing approach. Journal of Sensors, 2018, 2018. https://doi.org/10.1155/2018/1058513
  16. 16. Zhang W, CHEN F. Impacts of grading rule on urban thermal landscape pattern research[J]. Sustainability, 2018, 10(7): 2514. https://doi.org/10.3390/su10072514
  17. 17. Hong J, Xiao C, Zhi L, et al. The investigation of spatiotemporal patterns of landscape fragmentation during rapid urbanization in Changsha city. Remote Sensing for Land & Resources, 2011, 23(1): 133–137. https://doi.org/10.6046/gtzyyg.2011.01.27
  18. 18. Yu H, Zhou T, Li A, et al. Landscape changes and ecological effects on Dali city, Yunnan province. Geo-Informatics in Resource Management and Sustainable Ecosystem. Springer, Berlin, Heidelberg, 2015: 22–31. https://doi.org/10.1007/978-3-662-49155-3_3
  19. 19. Li Y. Thinking of Compact City. Urban Planning Forum 2007, 1: 87–90. (in Chinese)
  20. 20. Sun C, Lin T, Zhao Q, et al. Spatial pattern of urban green spaces in a long-term compact urbanization process—A case study in China. Ecological indicators, 2019, 96: 111–119. https://doi.org/10.1016/j.ecolind.2017.09.043
  21. 21. Sun L S, Han Y H, Xie Z W. Neighborhood extremum method of extracting urban built-up areas using night lighting Data. Geomatics and Information Science of Wuhan University,2019. https://doi.org/10.13203/j.whugis20190010.
  22. 22. Su Y, Chen X, Wang C, et al. A new method for extracting built-up urban areas using DMSP-OLS nighttime stable lights: A case study in the Pearl River Delta, southern China. GIScience & Remote Sensing, 2015, 52(2): 218–238. https://doi.org/10.1080/15481603.2015.1007778
  23. 23. Xie Z, Ye X, Zheng Z, Li D, Sun L, Li R, et al. Modeling Polycentric Urbanization Using Multisource Big Geospatial Data. Remote Sensing, 2019, 11(3): 310. https://doi.org/10.3390/rs11030310
  24. 24. Proville J, Zavala-Araiza D, Wagner G. Night-time lights: A global, long term look at links to socio-economic trends. PloS One, 2017, 12(3): e0174610. pmid:28346500
  25. 25. Liu J, Zhao G, Ni Y. Overview of researches based on urban extracting Via DMSP/OLS night-time light data: Taking Fuzhou as an example. Geomatics & Spatial Information Technology, 2017, 40(9), 175–179. (in Chinese)
  26. 26. Yan M C. Research and Contrast on Several Vegetation Classification Methods of Hight-Resolution Satellite Image Data. Journal of Remote Sensing, 2007,11(02): 93–98. (in Chinese)
  27. 27. Karanam H K, BabuNeela V. Study of normalized difference built-up (NDBI) index in automatically mapping urban areas from Landsat TM imagery. International Journal of Scientific Research and Review, 2017, 6: 239–248.
  28. 28. Wang R, Li J, Li X, Chen X. Extraction of built up areas using combination of Landsat and DMSP data—a case study of Jiangxi Province. Journal of Central China Normal University (Natural Science), 2018, 52(1): 130–136. https://doi.org/10.19603/j.cnki.1000-1190.2018.01.021
  29. 29. Ma X, Tong X, Liu S, et al. Optimized sample selection in SVM classification by combining with DMSP-OLS, Landsat NDVI and GlobeLand30 products for extracting urban built-up areas. Remote Sensing, 2017, 9(3): 236. https://doi.org/10.3390/rs9030236
  30. 30. Yao W. Pixel level and feature level remote sensing image fusion methods and applications. Doctor, Dalian University of Technology, Dalian, China, 2011. (in Chinese)
  31. 31. Chen J. Study on the similarity of the industrial structure between the northeastern provinces in China. Master's, Jilin University, Changchun, China, 2012. (in Chinese)
  32. 32. Zhen X. Strategy research on the development problems of the small and micro Industrial enterprises in Liaoning province. Master, Jilin University, Changchun, China, 2014. (in Chinese)
  33. 33. Li Y. Analysis on the difference of economic development level in Jilin Province. Master, Jilin University, Changchun, China, 2017. (in Chinese)
  34. 34. Li X, Hu X M, Ma Y, et al. Impact of planetary boundary layer structure on the formation and evolution of air-pollution episodes in Shenyang, Northeast China. Atmospheric Environment, 2019, 214: 116850. https://doi.org/10.1016/j.atmosenv.2019.116850
  35. 35. Wang L, Wang J, Tan X, et al. Analysis of NOx pollution characteristics in the atmospheric environment in Changchun city. Atmosphere, 2020, 11(1): 30. https://doi.org/10.3390/atmos11010030
  36. 36. Yuan Y, Yang M, Wu J, et al. Assessing bus transit service from the perspective of elderly passengers in Harbin, China[J]. International Journal of Sustainable Transportation, 2019, 13(10): 761–776. https://doi.org/10.1080/15568318.2018.1512691
  37. 37. Shi L, Huang C, Sun H. Research on spatial-temporal pattern of Yulin economic development based on DMSP/OLS nighttime imagery. Journal of Geomatics, 2018 (1): 89–92. https://doi.org/10.14188/j.2095-6045.2016094
  38. 38. Wu J S, Li S and Zhang X W. Research on saturation correction for long-time series of DMSP-OLS nighttime light dataset in China. Journal of Remote Sensing, 2018, 22(4): 621–632. (in Chinese)
  39. 39. Liu Z, He C, Zhang Q, Huang Q, Yang Y. Extracting the dynamics of urban expansion in China using DMSP/OLS nighttime light data from 1992 to 2008. Landscape and Urban Planning, 2012, 106(1): 62–72. https://doi.org/10.1016/j.landurbplan.2012.02.013
  40. 40. Zhou W, Newsam S, Li C, et al. Learning low dimensional convolutional neural networks for high-resolution remote sensing image retrieval. Remote Sensing, 2017, 9(5): 489. https://doi.org/10.3390/rs9050489
  41. 41. Jia K, Wei X, Gu X, et al. Land cover classification using Landsat 8 operational land imager data in Beijing, China. Geocarto international, 2014, 29(8): 941–951. https://doi.org/10.1080/10106049.2014.894586
  42. 42. Ali M Z, Qazi W, Aslam N. A comparative study of ALOS-2 PALSAR and landsat-8 imagery for land cover classification using maximum likelihood classifier. The Egyptian Journal of Remote Sensing and Space Science, 2018, 21: S29–S35. https://doi.org/10.1016/j.ejrs.2018.03.003
  43. 43. Henits L, Jürgens C, Mucsi L. Seasonal multitemporal land-cover classification and change detection analysis of Bochum, Germany, using multitemporal Landsat TM data. International Journal of Remote Sensing, 2016, 37(15): 3439–3454. https://doi.org/10.1080/01431161.2015.1125558
  44. 44. Yang X, Zheng X Q, Chen R. A land use change model: Integrating landscape pattern indexes and Markov-CA. Ecological Modelling, 2014, 283: 1–7. https://doi.org/10.1016/j.ecolmodel.2014.03.011
  45. 45. Gong Y. The Spatial-temporal pattern evaluation of wetland landscape and its driving mechanism in Changsha. Doctor, Central South University of Forestry and Technology, Changsha, China, 2013. (in Chinese)