Next Article in Journal
An Exploratory Study Investigating Gender Effects on Using 3D Maps for Spatial Orientation in Wayfinding
Previous Article in Journal
Effect of the Long-Term Mean and the Temporal Stability of Water-Energy Dynamics on China’s Terrestrial Species Richness
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Landsat-Based Land Cover Change in the Beijing-Tianjin-Tangshan Urban Agglomeration in 1990, 2000 and 2010

1
State Key Laboratory of Remote Sensing Science, Institute of Remote Sensing and Digital Earth, Chinese Academy of Sciences, Beijing 100101, China
2
Department of Geographical Sciences, University of Maryland, College Park, MD 20742, USA
*
Author to whom correspondence should be addressed.
ISPRS Int. J. Geo-Inf. 2017, 6(3), 59; https://doi.org/10.3390/ijgi6030059
Submission received: 29 November 2016 / Revised: 17 February 2017 / Accepted: 20 February 2017 / Published: 24 February 2017

Abstract

:
Rapid urbanization dramatically changes the local environment. A hybrid classification method is designed and applied to multi-temporal Landsat images and ancillary data to obtain land cover change datasets. A support vector machine (SVM) classifier is used to classify multi-temporal Landsat Enhanced Thematic Mapper Plus (ETM+) images that were collected in 2000 at the pixel level. These images are also segmented with the mean shift method. The impervious surface is refined based on a combination of the segmented objects and the SVM classification results. The changed areas in 1990 and 2010 are determined by comparing the Thematic Mapper (TM) and ETM+ images via the re-weighted multivariate alteration detection transformation method. The TM images that were masked as changed areas in 1990 and 2000 are input into the SVM classifier. Land cover maps for 1990 and 2010 are produced by combining the unchanged area in 2000 with the new classes of the changed areas in 1990 and 2010. Land cover change has continuously accelerated since 1990. Remarkably, arable land decreased, while the impervious surface area significantly increased.

1. Introduction

Land cover and land use provide information that improves our understanding of the interactions between humans and the environment [1]. Land cover influences the energy balance and the carbon and hydrological cycles and, thus, plays an important role in global change research [2,3,4]. Land cover directly affects the physical characteristics of the land surface, such as soil moisture, albedo, temperature and transpiration, so many scientific studies require information regarding the spatial distribution and dynamic changes of land cover [5,6,7]. Many projects, such as the International Geosphere Biosphere Programme (IGBP), the Global Observations of Forest Cover and Global Observations of Land Dynamics (GOFC-GOLD) and the Global Rain Forest Mapping (GRFM) project, have been proposed to understand the land cover and land cover changes [8,9,10].
During the past two decades, China’s economic development has resulted in the rapid expansion of urban areas [11,12,13], and some urban agglomerations have gradually formed, such as the Beijing-Tianjin-Tangshan (BTT), Wuhan, Yangtze River Delta and Pearl River Delta urban agglomerations [14]. The dramatic land cover changes around urban agglomerations cause a variety of problems, such as regional climate change, the urban heat island effect and increased greenhouse gas emissions [15,16,17,18,19]. Many climate models, such as the Regional Atmospheric Modeling System (RAMS), the Penn State University/National Center for Atmospheric Research (PSU/NCAR) mesoscale model (MM5) and the Weather Research and Forecasting Model (WRF), require land cover as input data to analyze the effects of land cover changes on the climate [20,21,22,23,24]. BTT is the largest urban agglomeration in North China [14]. Using climate models to simulate the effects of land cover changes on the regional climate in the BTT region requires long-term land cover products.
Currently, many land cover and land use datasets that cover the BTT region are freely available. Datasets with resolutions from 300 m–1000 m were primarily developed based on Medium Resolution Imaging Spectrometer (MERIS), Advanced Very High Resolution Radiometer (AVHRR), Moderate Resolution Imaging Spectroradiometer (MODIS) and Satellite Pour l’Observation de la Terre (SPOT) Vegetation [25,26,27,28]. These datasets represent the land surface characteristics in different years, but none of these datasets are continuous from 1990–2010. Furthermore, datasets with coarse resolution have low classification accuracy, especially in urban areas [21,29,30]. The information from classification products with 30-m spatial resolution is more detailed than that from coarser products, but the number of classes in classification products with 30-m resolution is generally lower than what is required by climate models [1,31,32]. In addition to the above global- and national-scale datasets, many land cover classifications of Beijing have been conducted, but the other regions in the BTT urban agglomeration have rarely been addressed [33,34,35,36].
Differences between land cover products are more susceptible to classification methods and classification systems [37]. The land covers must have identical semantics or extremely high accuracy to monitor any changes [38]. The above-mentioned land covers have obvious differences in thematic classes, classification methods and data. Comparing inconsistent maps creates unreliable conclusions [38]. Consistent land covers in 1990, 2000 and 2010 must be developed to understand the land cover changes from the expansion of the BTT urban agglomeration.
Images from remote sensing satellites have become established as the main basis for the development of land covers [26,28,39]. Classifications can be divided into pixel-based and object-based classifications based on pixels or pixel clusters [40]. The ‘salt and pepper’ problem often occurs in pixel-based classifications, so object-based classifications were developed to reduce the noise [35,41,42,43]. Segmentation is a critical step in object-based classifications. The segmentation of an image combines a group of spectrally-similar and spatially-adjacent pixels as an object [44,45,46]. Generally, some parameters must be set for image segmentation, such as the spatial scale, bandwidth and shapes; however, selecting one parameter that is suitable for all land cover classes is difficult [35]. Meanwhile, integrating expert knowledge and the effective attributes of the objects to be classified is time consuming [47]. Therefore, methods to combine pixel-based and object-based approaches have been used in some classifications [31,48,49,50].
Change detection is an important component of land cover studies. As a post-classification method, the comparison of classification results can be performed to detect changes [51]. Furthermore, pre-classification algorithms, such as image differencing, change vector analysis, image regression and image ratios, have been used to distinguish changed areas between images that are acquired on different dates [52,53,54]. Post-classification change detection could reduce the influence of radiometric calibration and sensor differences between different images, but the detection accuracy greatly depends on the classification accuracy of each land cover [55,56]. In contrast, pre-classification approaches can avoid the cumulative errors of the two classification results because only the changed area is updated. Generally, pre-classification approaches require images that are acquired on similar dates, but the iteratively re-weighted multivariate alteration detection transformation (IR-MAD) algorithm can be used to combine the relative radiometric matches and change detection [57]. The IR-MAD algorithm has no strict limit on acquiring images in the same phenological period. In this study, the IR-MAD algorithm is used to identify changed areas.
The aim of this study is to identify consistent land covers in 1990, 2000 and 2010 in the BTT urban agglomeration by using Landsat images. The classification system and accuracies of the land covers should satisfy the requirements of the RAMS climate model. To achieve this goal, we adopted a hybrid classification approach that integrated the advantages of both pixel-based and object-based approaches to classify multi-temporal Landsat images in 2000. The land covers in 1990 and 2010 were produced based on the land cover in 2000 by using the pre-classification change detection method. The characteristics of the urbanization in the BTT region were analyzed based on these three land covers.

2. Materials

2.1. Study Area

The BTT urban agglomeration (38°28′–41°3′ N, 115°23′–119°48′ E) includes 5 cities: Beijing, Langfang, Tianjin, Tangshan and Qinhuangdao. As the Chinese capital, Beijing promotes the economic development of its surrounding areas. Similarly, as a municipality, Tianjin represents the economic center of the Pohai Bay region. The influence of Beijing and Tianjin has served to accelerate the economic development of the adjacent cities. The location of the study area is shown in Figure 1. The study area is surrounded by the Taihang Mountains to the west, the Yanshan Mountains to the north and the Pohai Sea to the east. The main portion of the study area is arable land and is located in the northeastern area of the North China plain. Wheat and maize are widely planted in the plains areas, whereas forest and shrub are the major classes of vegetation in the mountainous areas.

2.2. Data and Pre-Processing

The Landsat satellite series has continually observed Earth since 1972 and has accumulated an enormous number of time series images [38]. Landsat images have been widely used in land cover classification because of their stable imaging quality [1,31,32,58,59,60]. In this study, Level 1 terrain-corrected (L1T) Landsat images were selected as the primary data source for the land cover classifications. The cloud coverage in the selected Landsat images was less than 5%. The Landsat images were atmospherically corrected to produce the surface reflectance by using the Landsat Ecosystem Disturbance Adaptive Processing System (LEDAPS). LEDAPS uses the 6SV (Second Simulation of a Satellite Signal in the Solar Spectrum vector code) radiative transfer code. The aerosol characterizations are derived independently from each Landsat acquisition, assume a fixed continental aerosol type and use ancillary water vapor [61,62,63]. The L1T Landsat images were processed to surface reflectance products by applying LEDAPS. The influence of the difference in the illumination geometry in the Thematic Mapper (TM)/ETM+ images was substantially eliminated.
Covering the entire study area required TM/ETM+ images from six scenes. The images that were used in this study are listed in Table 1. The images that were acquired in 1990 and 2010 were Landsat 5 TM images, and those that were acquired in 2000 were Landsat 7 ETM+ images. The data quality of Landsat 7 ETM+ images is superior to that of Landsat 5 TM images, so three Landsat 7 ETM+ images were used for every scene in 2000 [60,64]. Multi-temporal ETM+ images were utilized in the classification for 2000 to reduce spectral confusion. The scene selection was based on the vegetation dynamics of the target land cover types over a growing season. Crops in the study area could be harvested twice in one year. However, some arable lands were cultivated only once in one year. Multi-temporal images can compensate for cropland classification errors from harvesting or an absence of planting and capture agricultural changes during the crop’s life [65]. Multi-temporal images also make full use of phenological differences in different classes. The three acquisition dates of the ETM+ images covered the early, peak and late growing seasons [66]. The choice of appropriate acquisition dates would facilitate maximizing the differences between types. The major green vegetation in the BTT region was agricultural crop. Considering the phenology of the crops and the maximum separability among dry farmland, paddy field and forest, the three acquired dates were in March, August and October. We experienced difficulty selecting three temporal images in one year because of the influences of clouds and rain. Therefore, high-quality images that were acquired in the previous or following year were selected. One Landsat 5 TM image was selected for every scene in 1990 and 2010. An extra TM image was added for scene 121/33 in 1990 to compensate for the influence of cloud coverage.
In addition to the Landsat images, which constituted the primary data source, Shuttle Radar Topography Mission (SRTM) digital elevation model (DEM) data were included in the classification [67]. Forest, shrub and grass were easily confused with agricultural crops in the study area because of the spectral similarity in the Landsat images. However, these forest, shrub and grass were mainly distributed in mountainous areas, and agricultural crops were mainly distributed in plains in the study area. A large difference in slope was observed between the two terrains. In addition, the Conversion of Cropland to Forest Program (CCFP) has been implemented by the Chinese government since 1999 [68]. In the CCFP, cultivated land along steep slopes is converted into forest, shrub and grass [69]. Therefore, any forest, shrub and grass on mountains exhibits obvious slope differences between cultivated lands. Raw three-arc-second-resolution DEM data were mosaicked and re-projected to the same projection as the TM/ETM+ images and then cut to the size same as each scene. The DEM was resampled to 30-m horizontal resolution by using the nearest neighbor resampling method. The Landsat images and the slopes from the DEM were used as the input data for the SVM classifier. Only the Landsat images were used when using the SVM to train the samples, and the accuracy of 5-fold cross-validation was 87.98%. When the slope data were added, the accuracy reached 92.63%.

3. Methods

The study area included many rural settlements that belonged to the impervious surface class. During the growing season, these rural settlements were extensively sheltered by trees that were planted in streets or courtyards. During the dormant season, the spectral characteristics of the rural settlements were similar to those of the bare arable land. This complicated environment reduced the classification accuracy of rural settlements when using pixel-based classification methods, and rural settlements were typically subject to omission errors. The segment object was the cluster of pixels of various classes. The characteristics that described a single pixel were not suitable to the object. Using the SVM to classify the segment objects required more characteristics. Image segmentation was performed and combined with the SVM classification results to improve the classification accuracy of the impervious surface class [31,70].
The flow diagram of the classification process is shown in Figure 2. During classification, the multi-temporal ETM+ images from 2000 and the slope data were first classified scene-by-scene by the SVM classifier. The pixel-based and object-based approaches were integrated to improve the accuracy of the impervious surface. Therefore, the mean shift algorithm was used to segment the three multi-temporal Landsat ETM+ images to generate objects. The segmented objects were combined with the SVM classification result to revise the impervious surface class. The pre-classification change detection method was used to find the changed area; thus, one of the three images from 2000 was selected for comparison with images from 1990 and 2010 by using the IR-MAD approach [52,57]. Only the areas on the Landsat images from 1990 and 2010 that were masked by the changed area were classified with the SVM classifier. The land cover for 2000 was updated with the classification results of the changed area to develop the land cover for 1990 and 2010.

3.1. Classification System and Training Samples

Considering the spatial resolution of the TM/ETM+ images and the requirements of some climate models, the classification system comprised 9 classes that referenced some existing classification systems: forest, shrub, grass, dry farmland, paddy field, bare land, impervious surface, aquatic vegetation and water bodies.
The training samples were primarily collected by visual interpretation of the ETM+ images from 2000 with reference to Google Earth. The rest were collected from field observations in 2000. The sample points that were collected by visual interpretation were first randomly distributed in the six scenes. Then, a polygon that included pixels with the same class was manually drawn at a single point location. If the random sample was not in a pure pixel cluster, we selected a nearby polygon with the same class as the random points. Some land cover classes, such as water bodies, bare land and paddy fields, had lower area proportions in the study area. If only random samples were used, the sample counts of these classes would be significantly lower than those of the other classes. More samples were added by artificial selection for the classes with lower area proportions to mitigate the influence of these unbalanced sample counts. Finally, 327 polygons (26893 pixels) were identified.
The samples were identified by visual interpretation and drawing polygons, so the classes of samples may have included mistakes. All of the samples were input into the SVM classifier that was trained with the same samples to check for mistakes. If the predicted class was not the same as that input, the class of the sample was subjected to additional checks. If the class of the sample was correct according to Google Earth images, the class was retained, even if the SVM classifier output a different class. If the class of the sample was incorrect, the class was revised. This process was repeated until we were confident that the samples included no mistakes.

3.2. SVM Classification

Recently, the SVM classifier has been widely used in remote sensing classification with good results [71,72,73]. SVM uses the kernel function to map the spectral feature space to high-dimensional space and calculates the inner product in this high-dimensional space. SVM performs classification by using an optimal hyperplane, which is determined based on the limited supporting samples [74]. SVM can achieve high classification accuracy with relatively small numbers of training samples [73]. In this study, the Libsvm (Version 3.18) software package was used [75]. This software package can normalize the input data and search for the optimal parameters. Bands 1–5 and 7 of the Landsat images were set as the main input data for the SVM classifier. The slope data from the DEM were also included in the classification. Therefore, the input data for the year 2000 had 19 layers, including one slope layer and 18 bands (each scene corresponded to 6 bands) from the Landsat images of three scenes. Each scene was individually classified, and the results of six scenes were mosaicked together.
During classification, the 19 layers were normalized by using the scale function, and the radial basis function (RBF) was set as the kernel function. The Libsvm software package performed numerous classification tests to search for the optimal parameters by setting the searching step length of the parameters. Specifically, the samples of every scene were first scaled to the range [−1, +1]. Then, svm_type was set as c_svc; the 10-fold cross-validation was performed; and the step lengths of the penalty parameter and the gamma parameter in RBF were set to 1. The iteration termination threshold was 0.001. SVM classification was conducted once the optimal model parameters were determined.

3.3. Image Segmentation

The mean shift algorithm has been widely used in image segmentation and identification tracking [76,77]. As a kernel density estimate method, the mean shift does not require prior knowledge of the number of segmentations and sets no limits on the segmented shape [78]. In this study, the EDISON V1.1 software was used. The 18 reflective bands of the three scenes of Landsat images were transformed based on principal component analysis (PCA) to reduce the calculation iterations of the mean shift. The first three principal components were selected as the segmented layers. The Gaussian function was taken as the kernel function, and the bandwidths of the geometry space and color space were set to 15 and 9.5, respectively.
In this study, the segmented object was not set as the input data for the SVM classifier. A set of appropriate thresholds was used according to the area percentages of the classes within an object to determine whether the object belonged to the impervious surface class. Forest and dry farmland were the main factors because the interference factors to the rural settlements mainly included the trees and dry farmland around or inside the rural settlements. The segmented objects were overlaid with the SVM classification to refine the impervious surface. Generally, if the area fractions of impervious surface (Ai), forest (Af) and dry farmland (Ad) within an object met the judgment conditions, the polygon was designated as an impervious surface. These thresholds were selected based on multiple tests. The following rules were designed through the stratified extraction of the impervious surface:
If Ai ≥ 50%, class = impervious surface,
If 40% ≤ Ai < 50% and 40% ≤ Af < 50%, class = impervious surface,
If 30% ≤ Ai < 40% and 30% ≤ Af < 40% and 10% ≤ Ad < 30%, class = impervious surface
where Ai, Af and Ad are the area fractions of impervious surface, forest and dry farmland, respectively.
The rural settlements were mainly subjected to omission errors, so objects were classified as impervious surface when the Ai pixels were greater than 50% in a segmented object. If an object had more than 40% impervious surface pixels and forest pixels, the total pixels of these two classes were more than 80%. Basically, the rural settlements matched this situation. The majority of impervious surface objects were identified through these two equations. The proposed Equation (3) considered that the proportions of rural settlements, trees and farmland were all relatively low.
Furthermore, the area percentage of the small impervious surface patch may have been less than 30% if a small impervious surface patch was located in a large object. Therefore, the impervious surface patch would be deleted, even if the impervious surface patch was correctly classified. An extra rule was set to avoid this situation: if one impervious surface patch comprised more than 10 adjacent pixels, the impervious surface patch would be reserved regardless of the percentages of the classes in one object. If an impervious surface patch was deleted, the adjacent class with the largest area would be assigned to the deleted location.

3.4. Land cover Change

The IR-MAD algorithm was used to detect the changed area, and then, the TM images of the changed area were classified by the SVM classifier. IR-MAD judged the unchanged pixels by iteratively calculating the chi-square distribution of the difference in the bi-temporal images [52,57]. The IR-MAD linearly transformed the original images R and T to new images U and V based on the canonical correlation analysis (CCA) in Equation (4):
U = a T R ,   V = b T T ,
where U and V are the canonical correlation images, a T and b T are the transform matrixes and R and T are the bi-temporal images.
The bands of the canonical correlation images were referred to as the canonical variates. The canonical variates were arranged according to the correlations. If the correlation was relatively low, this pair of canonical variates had more change information. The differences in each pair of canonical variables were mutually uncorrelated.
According to the central limit theorem, the differences in the canonical variates approximately corresponded to a Gaussian distribution, and the sum of squares of the differences in the canonical variates corresponded to a chi-square distribution [52,57]. The differences in the canonical variates and the weights of unchanged pixels were calculated by using Equations (5)–(7):
M i = U N i + 1 V N i + 1 ,   i = 1 , 2 , ... , N ,
Z = i = 1 N ( M i σ M i ) 2 ,
Pr   ( n o   c h a n g e ) = 1 P χ 2 ; N ( Z ) ,
where i is the band of the image, M i is the difference in the variates, σ M i is the standard deviation of M i , N is the total number of bands, P χ 2 ; N ( Z ) is the χ 2 distributed with N degrees of freedom and Pr   ( n o   c h a n g e ) is the weight.
Pr   ( n o   c h a n g e ) was used to weigh each pixel. Then, Equations (4)–(7) were iterated until no significant changes were observed in the canonical correlations. These pixels could then be determined as unchanged pixels if their chi-square distribution probability was lower than 0.9 [52].
In this study, we calculated the chi-square images from 2010 and 1990 based on the images from 2000. The threshold was slightly lower than 0.95, which was set by the developer, so some unchanged pixels were selected in the identified area. Although some extra unchanged pixels were selected, mostly changed pixels were chosen, so the omission error was reduced.
The samples were collected from the unchanged area to classify the changed pixels with the SVM classifier. The land covers from 1990 and 2010 were produced by updating the land cover from 2000 with the classification of the changed areas for 1990 and 2010. Figure 3 shows this process.

4. Results

4.1. Accuracy Assessment

The assessment samples were selected through stratified random sampling. Bare land and aquatic vegetation had the smallest area proportions, so we supplemented some points to these two classes. The class of the assessment sample was identified by visual interpretation. A total of 526 points from these nine classes were collected to assess the classification accuracy. The raw assessment points and error matrixes were shown in supplementary materials. The error matrix was used to calculate the overall accuracy, Kappa coefficient, user’s accuracy and producer’s accuracy for the three periods [79]. The overall accuracy is the proportion of the area that was correctly mapped. The user’s accuracy is the proportion of correctly-classified pixels with regard to all of the pixels that are classified as this class in the classified image. The producer’s accuracy is the proportion of correctly-classified pixels with regard to all of the pixels of that ground truth class [80]. The overall accuracy, user’s accuracy and producer’s accuracy were adjusted in terms of the unbiased estimator of the area proportions [80]. The specific calculation process is shown in Equations (8)–(11) [81]. The results of the accuracy assessment are shown in Table 2.
p ^ i j = ( N i . N ) ( n i j n i . )
where N i . is the number of pixels of class i in the entire map, N is the number of pixels in the entire map, n i . is the number of pixels from a sample of class i and n i j is the sample count in row i (map category) and column j (reference category) in the error matrix of the sample counts.
p ^ PA = p ^ i j p . j
p ^ U A = p ^ i j p i .
p ^ O A = p ^ k k
where p ^ PA is the corrected producer’s accuracy, p ^ UA is the corrected user’s accuracy and p ^ OA is the corrected overall accuracy. p . j is the sum of the p ^ ij in column j of the estimated error matrix. p i . is the sum of the p ^ ij in row i of the estimated error matrix. p ^ kk is the diagonal value of the estimated error matrix.
The land cover in 2000 had the highest overall accuracy and Kappa coefficient, which were 89.13% and 0.87 respectively. The overall accuracies and Kappa coefficients were 87.94% and 87.76% and 0.85 and 0.85 for 1990 and 2010, respectively. Overall, the accuracy for 2000 was higher than that for 1990, and the accuracy for 1990 was higher than that for 2010; however, the difference was minor. We were more concerned about the dry farmland, forest and impervious surface in the nine classes. Indeed, these three classes comprised approximately 70% of the study area and were more seriously affected by urbanization. The user’s accuracy and producer’s accuracy of these three classes all exceeded 80%.
A sub-area was selected to compare the land cover in 2000 with China’s Land Use/cover Dataset (CLUDs) for the year 2000 [32]. This sub-area includes mountains and plains in Landsat Path/Row 123/032. The land covers that were classified in this study and the CLUDs are shown in Figure 4. Figure 4a–c shows the Landsat images in 1990, 2000 and 2010. Figure 4d–f shows the land covers that were classified in this study for Figure 4a–c. Figure 4g is the CLUD for the year 2000. This CLUD was derived by visual interpretation [32]. When comparing the land cover in 2000 with the CLUD for the year 2000, the main differences were the distribution of grass and shrub. The interpreted boundaries of grass and shrub in this CLUD were artificial drawings. These polygons, which were drawn by humans, could not adequately present the interleaving and transitions among vegetation. However, if no water was present in a riverbed, this area would be classified as bare land or grass in the land cover for the year 2000. In the CLUD for the year 2000, this riverbed was interpreted as a river. This difference shows the flexibility in visual interpretation.

4.2. Impervious Surface Refinement

The impervious surface class from the SVM included many incorrect trivial speckles because this class is affected by the complexity of the environment. However, integration with segmented objects effectively reduced these incorrect trivial speckles.
The impervious surface of a sub-image was visually interpreted with reference to the high-resolution images to evaluate the effect of the impervious surface refinement. Detailed information is shown in Table 3. In this table, ‘no change’ means that no changes occurred; ‘add correct’ indicates that the impervious surface was correctly added to the SVM classification; ‘add error’ indicates that the impervious surface was incorrectly added; and ‘remove correct’ and ‘remove error’ refer to correctly or incorrectly deleted impervious surface, respectively. The percentage was calculated based on the area of impervious land that was identified by the SVM. Overall, the percentage of correct refinement far exceeded the percentage of incorrect refinement. In this sub-area, the percentages of impervious surface pixels that were added and removed correctly were 18.99% and 13.48%, respectively. The corresponding percentages for incorrect addition and removal were 2.89% and 1.99%, respectively. Therefore, the accuracy of the impervious surface classification was improved by 27.59%.
Figure 5 shows the local impervious surface that was classified by the SVM, the objects that were segmented by the mean shift method and the refined impervious surface after synthesizing the SVM and mean shift results. Figure 5a is a false-color Landsat image. Figure 5b shows the vector boundary of the segmented objects when overlaid onto the sub-image in Figure 5a. The outlines of the impervious surface are clearly delineated. In Figure 5c, the impervious surface from the SVM contains many trivial speckles. Figure 5d shows the refined impervious surface. Detailed information regarding the integration is presented in Figure 5e.

4.3. Change Mask

A threshold of 0.9 was manually chosen for the chi-square image based on the IR-MAD algorithm to extract the changed area. Figure 6 shows the sub-images, which contain complex changes in various classes. Figure 6a,b is the false-color Landsat sub-image that were acquired on 2 August 1999 and 22 September 2009, respectively. The main classes in the sub-images were an airport, towns, rural settlements, dry farmland and rivers. The airport site at the center grew almost two-fold from 1999–2009. The bright cyan areas correspond to buildings and roads. Forests are shown in dark red in the upper corner along the river. A substantial amount of dry farmland (bright red or light cyan) became impervious surface. Figure 6c is the chi-square image; the grayscale values reflect the intensity of the changes. The chi-square image exhibited a clear response to these class changes. Figure 6d is the mask of the changed area in yellow, and the background is the sub-image in Figure 6a. Clearly, dramatic changes occurred around the airport.

4.4. Land Cover and Change

The land cover in the BTT region has dramatically changed over the past two decades. The total areas of the study region in 1990, 2000 and 2010 were 55,508 km2, 55,660 km2 and 55,855 km2, respectively. Thus, the area increased by 152 km2 and 195 km2 from 1990–2000 and from 2000–2010, respectively. The increased areas corresponded to sea reclamation and thus were concentrated in the coastal region. The increased land from the sea mainly included impervious surface and bare land. The land cover maps for the three periods are shown in Figure 7, and the local details are shown in Figure 4.
The areas and proportions of the classes were adjusted to correct for the estimated error based on the assessment samples. The adjusted areas and proportions are shown in Table 4. The proportions of each class for the three periods are also shown in Figure 8.
The sum of the areas of dry farmland, forest and impervious surface always exceeded 70% of the total study area during the two analyzed decades. The most remarkable changes were the reduction in arable land and the increase in impervious surface between 1990 and 2010. Dry farmland changed the most significantly from 1990–2000, shrinking by 1260 km2, with its proportion of the total area decreasing from 47.1% down to 44.7%. Subsequently, dry farmland decreased further to 42.09% from 2000–2010. The area of paddy fields was approximately 5.76% in 1990 and decreased to 4.92% and 2.95% in 2000 and 2010, respectively. In contrast, the area of impervious surface increased the most during the studied period. The area of impervious surface increased from 8.72%–11.45% from 1990–2000 and then to 16.94% by 2010, nearly doubling in only two decades. The increased area of impervious surface was mostly attributable to the conversion of arable land. Indeed, approximately 90% of the increased area of impervious surface between 1990 and 2000 was originally arable land; from 2000–2010, this value was approximately 80%.
Compared to the changes seen in the areas of dry farmland and impervious surface, the changes in the areas of other classes were relatively small. The forest area exhibited a slight decrease from 1990–2000, followed by an increase of 0.72% by 2010. The decreased forest area was mainly converted to dry farmland and impervious surface from 1990–2000. Beginning in 1999, China began to convert cropland to forest, so some dry farmland with steep slopes and some shrub land were converted to forest land. Although water bodies had some increments (i.e., the enclosed sea), the total area that was covered by water bodies continuously decreased. Clearly, the water levels of the reservoirs in this area declined.

5. Discussion

5.1. Integration of Pixel-Based and Object-Based Classification Methods

The combination of pixel-based classification and object-based image segmentation was critical in this study. The land surface characteristics of the study area are very complex. For example, some agricultural land was not planted, so its spectral characteristics were similar to those of bare land and buildings. Many greenhouses were located on agricultural land and appeared similar to rural settlements in the multi-spectral images. Furthermore, buildings in some rural settlements were sheltered by trees. This situation is not as prevalent in urban areas because buildings are generally higher than trees. During classification, some rural settlements were classified as several disconnected components when using the SVM classifier. However, the buildings and trees in these rural settlements could be segmented as one object. These objects could comprise one or several types of land cover pixels, so the spectral characteristics of these objects were more complex than those of single pixels. Therefore, the classification of these objects required more information, such as the shape, perimeter, area and location relative to other objects. When we attempted to use some object features, such as the mean, median, maximum and minimum, as the input data for the SVM classifier, the classification results were not improved relative to those that were only based on pixels. Indeed, the mean, maximum and minimum values of the objects widely varied, especially for objects with pixels of several classes. Thus, the identification of objects requires more knowledge, such as the relationship between objects, the distribution pattern of pixels in an object and the shape of the object. Substantial additional work is required to explore the knowledge for every class to classify objects with the SVM classifier.
Image segmentation generally requires segmentation parameters to be set. However, large differences in spatial scale existed between the land cover classes, so properly segmenting the images for all of the land cover classes with one spatial scale parameter was difficult. In this study, the segmentation parameters were set to a scale that was suitable for rural settlements to improve the classification accuracy of the impervious surface. We mainly considered the proportions of the buildings, trees and agricultural land in a single object to combine the advantages of both pixel-based and object-based approaches. In this study, the error that was associated with impervious surface based on the SVM classifier was mainly omission error. Therefore, the thresholds were set to supplement the impervious surface when combining the SVM classification results and segmented objects. Finally, the gaps in impervious surface were filled, and the small misclassified patches were removed. The integrating rules could be easily understood, but the rules were constructed by many experiments. These rules were empirical rules and should be redesigned for different regions or different types. In addition, types with large-scale differences will not be adequately segmented because of the limitations of the segmentation-scale parameters. Self-adaptive segmentation could be a potential solution.

5.2. Change Detection

Using a consistent classification method is critical to ensure that the classification results for different periods are comparable. The spectral characteristics of different land covers usually have some confusion or overlap and are not always separated by a clear-cut boundary. The transitions between some classes, such as between bare land and sparse grassland and between wetlands and water bodies, are continuous. Classification errors often appear near the boundaries between these classes because of their spectral similarity. Each satellite image has unique Sun-illumination geometry and weather characteristics. Despite using the same classifier to classify images that are acquired at different dates, ensuring that the boundaries of the unchanged land cover in the transition region are similarly identified is difficult. This type of discrepancy in boundary-adjacent areas can result in errors when comparing land covers between different periods. Generally, only a portion of the land cover in a region changes over a given period. The error in land cover changes is lower if only the changed areas are updated and the unchanged areas are retained. However, updating the changed area particularly relies on the ability to accurately identify the changed area. If the identified changed area is substantially greater than the actual changed area, the effect of updating is the same as classifying the entire image. In contrast, the omission error increases if the identified changed area is less than the actual changed area. If the changed area can be appropriately identified, however, updating only the changed area ensures that the new land cover product is consistent with the base product.
In this study, the threshold of changed pixels was set slightly lower than the recommended threshold by the software developer, so the identified changed area was slightly greater than the actual changed area. This choice was made to reduce the omission error. Overall, the method of updating the changed area reduced the required computation for the classification and improved the consistency of the land covers between different periods.

6. Conclusions

China’s urbanization process has continuously progressed, and some urban agglomerations have formed. This urbanization has led to more frequent and intense land cover changes, thus affecting the local climate. Multi-temporal Landsat images and ancillary data were classified to extract LCLUC information for the BTT urban agglomeration.
In this study, the SVM classifier was used to classify multi-temporal Landsat images that were acquired in 2000. The objects that were identified via image segmentation were integrated with the SVM land cover classification results to improve the classification accuracy of impervious surface. The change detection method IR-MAD was used to update the land cover in 1990 and 2010 based on the land cover in 2000. Finally, consistent land covers were developed for 1990, 2000 and 2010. Remarkably, the area of arable land in the BTT urban agglomeration decreased over the past two decades, while the impervious surface area increased. The area of impervious surface increased nearly two-fold from 1990–2010. Meanwhile, the decreased area of arable land, including dry farmland and paddy fields, was slightly less than the increased area of impervious surface. The method that was used in this study was proven to be effective after validating the classification results. Thus, this method could be used for similar applications.

Supplementary Materials

The following are available online at www.mdpi.com/2220-9964/6/3/59/s1, Table S1: Error matrix for the year 1990 (unit: point), Table S2: Error matrix for the year 2000 (unit: point), Table S3: Error matrix for the year 2010 (unit: point).

Acknowledgments

This research was funded by the National Natural Science Foundation of China (Grant Nos. 41201444 and 41201442).

Author Contributions

Aqiang Yang designed the framework of this research and finished the final version of this paper. Guoqing Sun supervised the research efforts.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Gong, P.; Wang, J.; Yu, L.; Zhao, Y.; Zhao, Y.; Liang, L.; Niu, Z.; Huang, X.; Fu, H.; Liu, S.; et al. Finer resolution observation and monitoring of global land cover: First mapping results with Landsat TM and ETM+ data. Int. J. Remote Sens. 2013, 34, 2607–2654. [Google Scholar] [CrossRef]
  2. Deng, X.; Shi, Q.; Zhang, Q.; Shi, C.; Yin, F. Impacts of land use and land cover changes on surface energy and water balance in the Heihe river basin of China, 2000–2010. Phys. Chem. Earth Parts A/B/C 2015, 79–82, 2–10. [Google Scholar] [CrossRef]
  3. Chen, Y.; Li, X.; Liu, X.; Ai, B. Analyzing land cover change and corresponding impacts on carbon budget in a fast developing sub-tropical region by integrating MODIS and Landsat TM/ETM+ images. Appl. Geogr. 2013, 45, 10–21. [Google Scholar] [CrossRef]
  4. Yira, Y.; Diekkrüger, B.; Steup, G.; Bossa, A.Y. Modeling land use change impacts on water resources in a tropical west African Catchment (Dano, Burkina Faso). J. Hydrol. 2016, 537, 187–199. [Google Scholar] [CrossRef]
  5. Pielke, R.A. Land use and climate change. Science 2005, 310, 1625–1626. [Google Scholar] [CrossRef] [PubMed]
  6. Buchanan, G.M.; Nelson, A.; Mayaux, P.; Hartley, A.; Donald, P.F. Delivering a global, terrestrial, biodiversity observation system through remote sensing. Conserv. Biol. 2009, 23, 499–502. [Google Scholar] [CrossRef] [PubMed]
  7. Bontemps, S.; Herold, M.; Kooistra, L.; van Groenestijn, A.; Hartley, A.; Arino, O.; Moreau, I.; Defourny, P. Revisiting land cover observation to address the needs of the climate modeling community. Biogeosciences 2012, 9, 2145–2157. [Google Scholar] [CrossRef] [Green Version]
  8. Turner Ii, B.L.; Skole, D.; Sanderson, S.; Fischer, G.; Fresco, L.; Leemans, R. Land-Use and Land Cover Change; Science/Research Plan; IGBP: Stockholm, Sweden, 1995; p. 132. [Google Scholar]
  9. Justice, C.O.; Smith, R.; Gill, A.M.; Csiszar, I. A review of current space-based fire monitoring in Australia and the gofc/gold program for international coordination. Int. J. Wildland Fire 2003, 12, 247–258. [Google Scholar] [CrossRef]
  10. Rosenqvist, A.; Shimada, M.; Chapman, B.; Freeman, A.; De Grandi, G.; Saatchi, S.; Rauste, Y. The global rain forest mapping project—A review. Int. J. Remote Sens. 2000, 21, 1375–1387. [Google Scholar] [CrossRef]
  11. Cai, H.; Yang, X.; Xu, X. Spatiotemporal patterns of urban encroachment on cropland and its impacts on potential agricultural productivity in China. Remote Sens. 2013, 5, 6443–6460. [Google Scholar] [CrossRef]
  12. Su, S.; Jiang, Z.; Zhang, Q.; Zhang, Y. Transformation of agricultural landscapes under rapid urbanization: A threat to sustainability in Hang-Jia-Hu Region, China. Appl. Geogr. 2011, 31, 439–449. [Google Scholar] [CrossRef]
  13. Li, B.; Chen, D.; Wu, S.; Zhou, S.; Wang, T.; Chen, H. Spatio-temporal assessment of urbanization impacts on ecosystem services: Case study of Nanjing city, China. Ecological Indicators 2016, 71, 416–427. [Google Scholar] [CrossRef]
  14. Zhang, Q.; Su, S. Determinants of urban expansion and their relative importance: A comparative analysis of 30 major metropolitans in China. Habitat Int. 2016, 58, 89–107. [Google Scholar] [CrossRef]
  15. Li, X.; Zhou, W.; Ouyang, Z.; Xu, W.; Zheng, H. Spatial pattern of greenspace affects land surface temperature: Evidence from the heavily urbanized Beijing metropolitan area, China. Landsc. Ecol. 2012, 27, 887–898. [Google Scholar] [CrossRef]
  16. Odindi, J.O.; Bangamwabo, V.; Mutanga, O. Assessing the value of urban green spaces in mitigating multi-seasonal urban heat using MODIS land surface temperature (LST) and Landsat 8 data. Int. J. Environ. Res. 2015, 9, 9–18. [Google Scholar]
  17. Chen, X.L.; Zhao, H.M.; Li, P.X.; Yin, Z.Y. Remote sensing image-based analysis of the relationship between urban heat island and land use/cover changes. Remote Sens. Environ. 2006, 104, 133–146. [Google Scholar] [CrossRef]
  18. Du, H.; Wang, D.; Wang, Y.; Zhao, X.; Qin, F.; Jiang, H.; Cai, Y. Influences of land cover types, meteorological conditions, anthropogenic heat and urban area on surface urban heat island in the Yangtze River delta urban agglomeration. Sci. Total Environ. 2016, 571, 461–470. [Google Scholar] [CrossRef] [PubMed]
  19. Sato, Y.; Higuchi, A.; Takami, A.; Murakami, A.; Masutomi, Y.; Tsuchiya, K.; Goto, D.; Nakajima, T. Regional variability in the impacts of future land use on summertime temperatures in Kanto region, the Japanese megacity. Urban For. Urban Green. 2016, 20, 43–55. [Google Scholar] [CrossRef]
  20. Cotton, W.R.; Pielke, R.A.; Walko, R.L.; Liston, G.E.; Tremback, C.J.; Jiang, H.; McAnelly, R.L.; Harrington, J.Y.; Nicholls, M.E.; Carrio, G.G.; et al. RAMS 2001: Current status and future directions. Meteorol. Atmos. Phys. 2003, 82, 5–29. [Google Scholar] [CrossRef]
  21. Loveland, T.R.; Reed, B.C.; Brown, J.F.; Ohlen, D.O.; Zhu, Z.; Yang, L.; Merchant, J.W. Development of a global land cover characteristics database and IGBP discover from 1 km AVHRR data. Int. J. Remote Sens. 2000, 21, 1303–1330. [Google Scholar] [CrossRef]
  22. McQueen, J.T.; Valigura, R.A.; Stunder, B.J.B. Evaluation of the RAMS model for estimating turbulent fluxes over the Chesapeake bay. Atmos. Environ. 1997, 31, 3803–3819. [Google Scholar] [CrossRef]
  23. Santos-Alamillos, F.J.; Pozo-Vázquez, D.; Ruiz-Arias, J.A.; Tovar-Pescador, J. Influence of land-use misrepresentation on the accuracy of WRF wind estimates: Evaluation of GLCC and corine land-use maps in Southern Spain. Atmos. Res. 2015, 157, 17–28. [Google Scholar] [CrossRef]
  24. Mahmood, R.; Leeper, R.; Quintanar, A.I. Sensitivity of planetary boundary layer atmosphere to historical and future changes of land use/land cover, vegetation fraction, and soil moisture in Western Kentucky, USA. Glob. Planet. Chang. 2011, 78, 36–53. [Google Scholar] [CrossRef]
  25. Hansen, M.C.; Defries, R.S.; Townshend, J.R.G.; Sohlberg, R. Global land cover classification at 1 km spatial resolution using a classification tree approach. Int. J. Remote Sens. 2000, 21, 1331–1364. [Google Scholar] [CrossRef]
  26. Friedl, M.A.; Sulla-Menashe, D.; Tan, B.; Schneider, A.; Ramankutty, N.; Sibley, A.; Huang, X. MODIS collection 5 global land cover: Algorithm refinements and characterization of new datasets. Remote Sens. Environ. 2010, 114, 168–182. [Google Scholar] [CrossRef]
  27. An, Y.; Zhao, W.; Zhang, Y. Accuracy assessments of the globcover dataset using global statistical inventories and fluxnet site data. Acta Ecol. Sin. 2012, 32, 314–320. [Google Scholar] [CrossRef]
  28. Bartholome, E.; Belward, A.S. GLC2000: A new approach to global land cover mapping from earth observation data. Int. J. Remote Sens. 2005, 26, 1959–1977. [Google Scholar] [CrossRef]
  29. Sertel, E.; Robock, A.; Ormeci, C. Impacts of land cover data quality on regional climate simulations. Int. J. Climatol. 2010, 30, 1942–1953. [Google Scholar] [CrossRef]
  30. Ge, J.; Qi, J.; Lofgren, B.M.; Moore, N.; Torbick, N.; Olson, J.M. Impacts of land use/cover classification accuracy on regional climate simulations. J. Geophys. Res. Atmos. 2007, 112. [Google Scholar] [CrossRef]
  31. Chen, J.; Chen, J.; Liao, A.; Cao, X.; Chen, L.; Chen, X.; He, C.; Han, G.; Peng, S.; Lu, M.; et al. Global land cover mapping at 30 m resolution: A pok-based operational approach. ISPRS J. Photogram. Remote Sens. 2015, 103, 7–27. [Google Scholar] [CrossRef]
  32. Liu, J.; Kuang, W.; Zhang, Z.; Xu, X.; Qin, Y.; Ning, J.; Zhou, W.; Zhang, S.; Li, R.; Yan, C.; et al. Spatiotemporal characteristics, patterns, and causes of land-use changes in China since the late 1980s. J. Geogr. Sci. 2014, 24, 195–210. [Google Scholar] [CrossRef]
  33. Tian, Y.; Yin, K.; Lu, D.; Hua, L.; Zhao, Q.; Wen, M. Examining land use and land cover spatiotemporal change and driving forces in Beijing from 1978 to 2010. Remote Sens. 2014, 6, 10593–10611. [Google Scholar] [CrossRef]
  34. Jia, K.; Liang, S.; Wei, X.; Yao, Y.; Su, Y.; Jiang, B.; Wang, X. Land cover classification of Landsat data with phenological features extracted from time series MODIS NDVI data. Remote Sens. 2014, 6, 11518–11532. [Google Scholar] [CrossRef]
  35. Yu, W.; Zhou, W.; Qian, Y.; Yan, J. A new approach for land cover classification and change analysis: Integrating backdating and an object-based method. Remote Sens. Environ. 2016, 177, 37–47. [Google Scholar] [CrossRef]
  36. Wu, W.; Zhao, S.; Zhu, C.; Jiang, J. A comparative study of urban expansion in Beijing, Tianjin and Shijiazhuang over the past three decades. Landsc. Urban Plan. 2015, 134, 93–106. [Google Scholar] [CrossRef]
  37. McCallum, I.; Obersteiner, M.; Nilsson, S.; Shvidenko, A. A spatial comparison of four satellite derived 1 km global land cover datasets. Int. J. Appl. Earth Obs. Geoinf. 2006, 8, 246–255. [Google Scholar] [CrossRef]
  38. Sexton, J.O.; Urban, D.L.; Donohue, M.J.; Song, C. Long-term land cover dynamics by multi-temporal classification across the Landsat-5 record. Remote Sens. Environ. 2013, 128, 246–258. [Google Scholar] [CrossRef]
  39. Zhong, B.; Ma, P.; Nie, A.; Yang, A.; Yao, Y.; Lü, W.; Zhang, H.; Liu, Q. Land cover mapping using time series HJ-1/CCD data. Sci. China Earth Sci. 2014, 57, 1790–1799. [Google Scholar] [CrossRef]
  40. Lu, D.; Weng, Q. A survey of image classification methods and techniques for improving classification performance. Int. J. Remote Sens. 2007, 28, 823–870. [Google Scholar] [CrossRef]
  41. Li, M.; Ma, L.; Blaschke, T.; Cheng, L.; Tiede, D. A systematic comparison of different object-based classification techniques using high spatial resolution imagery in agricultural environments. Int. J. Appl. Earth Obs. Geoinf. 2016, 49, 87–98. [Google Scholar] [CrossRef]
  42. Yan, G.; Mas, J.F.; Maathuis, B.H.P.; Xiangmin, Z.; Van Dijk, P.M. Comparison of pixel-based and object-oriented image classification approaches—A case study in a coal fire area, wuda, Inner Mongolia, China. Int. J. Remote Sens. 2006, 27, 4039–4055. [Google Scholar] [CrossRef]
  43. Shiraishi, T.; Motohka, T.; Thapa, R.B.; Watanabe, M.; Shimada, M. Comparative assessment of supervised classifiers for land use–land cover classification in a tropical region using time-series PalSAR mosaic data. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2014, 7, 1186–1199. [Google Scholar] [CrossRef]
  44. Walter, V. Object-based classification of remote sensing data for change detection. ISPRS J. Photogram. Remote Sens. 2004, 58, 225–238. [Google Scholar] [CrossRef]
  45. Benz, U.C.; Hofmann, P.; Willhauck, G.; Lingenfelder, I.; Heynen, M. Multi-resolution, object-oriented fuzzy analysis of remote sensing data for GIS-ready information. ISPRS J. Photogram. Remote Sens. 2004, 58, 239–258. [Google Scholar] [CrossRef]
  46. Blaschke, T. Object based image analysis for remote sensing. ISPRS J. Photogram. Remote Sens. 2010, 65, 2–16. [Google Scholar] [CrossRef]
  47. Vieira, M.A.; Formaggio, A.R.; Rennó, C.D.; Atzberger, C.; Aguiar, D.A.; Mello, M.P. Object based image analysis and data mining applied to a remotely sensed Landsat time-series to map sugarcane over large areas. Remote Sens. Environ. 2012, 123, 553–562. [Google Scholar] [CrossRef]
  48. Guindon, B.; Zhang, Y.; Dillabaugh, C. Landsat urban mapping based on a combined spectral–spatial methodology. Remote Sens. Environ. 2004, 92, 218–232. [Google Scholar] [CrossRef]
  49. Geneletti, D.; Gorte, B.G.H. A method for object-oriented land cover classification combining Landsat TM data and aerial photographs. Int. J. Remote Sens. 2003, 24, 1273–1286. [Google Scholar] [CrossRef]
  50. Myint, S.W.; Gober, P.; Brazel, A.; Grossman-Clarke, S.; Weng, Q. Per-pixel vs. Object-based classification of urban land cover extraction using high spatial resolution imagery. Remote Sens. Environ. 2011, 115, 1145–1161. [Google Scholar] [CrossRef]
  51. Abd El-Kawy, O.R.; Rød, J.K.; Ismail, H.A.; Suliman, A.S. Land use and land cover change detection in the western nile delta of Egypt using remote sensing data. Appl. Geogr. 2011, 31, 483–494. [Google Scholar] [CrossRef]
  52. Canty, M.J.; Nielsen, A.A. Linear and kernel methods for multivariate change detection. Comput. Geosci. 2012, 38, 107–114. [Google Scholar]
  53. Xian, G.; Homer, C.; Fry, J. Updating the 2001 national land cover database land cover classification to 2006 by using Landsat imagery change detection methods. Remote Sens. Environ. 2009, 113, 1133–1147. [Google Scholar] [CrossRef]
  54. Lu, D.; Mausel, P.; Brondízio, E.; Moran, E. Change detection techniques. Int. J. Remote Sens. 2004, 25, 2365–2401. [Google Scholar] [CrossRef]
  55. Chen, Z.; Wang, J. Land use and land cover change detection using satellite remote sensing techniques in the mountainous Three Gorges area, China. Int. J. Remote Sens. 2010, 31, 1519–1542. [Google Scholar] [CrossRef]
  56. Chen, X.; Chen, J.; Shi, Y.; Yamaguchi, Y. An automated approach for updating land cover maps based on integrated change detection and classification methods. ISPRS J. Photogram. Remote Sens. 2012, 71, 86–95. [Google Scholar] [CrossRef]
  57. Canty, M.J.; Nielsen, A.A. Automatic radiometric normalization of multitemporal satellite imagery with the iteratively re-weighted mad transformation. Remote Sens. Environ. 2008, 112, 1025–1036. [Google Scholar] [CrossRef]
  58. Kim, D.-H.; Sexton, J.O.; Noojipady, P.; Huang, C.; Anand, A.; Channan, S.; Feng, M.; Townshend, J.R. Global, Landsat-based forest-cover change from 1990 to 2000. Remote Sens. Environ. 2014, 155, 178–193. [Google Scholar] [CrossRef]
  59. Hansen, M.C.; Egorov, A.; Potapov, P.V.; Stehman, S.V.; Tyukavina, A.; Turubanova, S.A.; Roy, D.P.; Goetz, S.J.; Loveland, T.R.; Ju, J.; et al. Monitoring conterminous United States (CONUS) land cover change with web-enabled Landsat data (WELD). Remote Sens. Environ. 2014, 140, 466–484. [Google Scholar] [CrossRef]
  60. Yang, L.; Huang, C.; Homer, C.G.; Wylie, B.K.; Coan, M.J. An approach for mapping large-area impervious surfaces: Synergistic use of Landsat-7 ETM+ and high spatial resolution imagery. Can. J. Remote Sens. 2003, 29, 230–240. [Google Scholar] [CrossRef]
  61. Ju, J.; Roy, D.P.; Vermote, E.; Masek, J.; Kovalskyy, V. Continental-scale validation of MODIS-based and LEDAPS Landsat ETM+ atmospheric correction methods. Remote Sens. Environ. 2012, 122, 175–184. [Google Scholar] [CrossRef]
  62. Maiersperger, T.K.; Scaramuzza, P.L.; Leigh, L.; Shrestha, S.; Gallo, K.P.; Jenkerson, C.B.; Dwyer, J.L. Characterizing LEDAPS surface reflectance products by comparisons with aeronet, field spectrometer, and MODIS data. Remote Sens. Environ. 2013, 136, 1–13. [Google Scholar] [CrossRef]
  63. Masek, J.G.; Vermote, E.F.; Saleous, N.E.; Wolfe, R.; Hall, F.G.; Huemmrich, K.F.; Gao, F.; Kutler, J.; Lim, T.-K. A Landsat surface reflectance dataset for North America, 1990–2000. IEEE Geosci. Remote Sens. Lett. 2006, 3, 68–72. [Google Scholar] [CrossRef]
  64. Tucker, C.J.; Grant, D.M.; Dykstra, J.D. Nasa’s global orthorectified Landsat data set. Photogram. Eng. Remote Sens. 2004, 70, 313–322. [Google Scholar] [CrossRef]
  65. Zheng, B.; Campbell, J.B.; de Beurs, K.M. Remote sensing of crop residue cover using multi-temporal Landsat imagery. Remote Sens. Environ. 2012, 117, 177–183. [Google Scholar] [CrossRef]
  66. Limin, Y.; Homer, C.; Hegge, K.; Chengquan, H.; Wylie, B.; Reed, B. A Landsat 7 scene selection strategy for a national land cover database. In Proceedings of the IEEE 2001 International Geoscience and Remote Sensing Symposium (IGARSS 2001), Sydney, Ausralia, 9–13 July 2001; Volume 1123, pp. 1123–1125.
  67. Satge, F.; Denezine, M.; Pillco, R.; Timouk, F.; Pinel, S.; Molina, J.; Garnier, J.; Seyler, F.; Bonnet, M.-P. Absolute and relative height-pixel accuracy of srtm-gl1 over the south american andean plateau. ISPRS J. Photogramm. Remote Sens. 2016, 121, 157–166. [Google Scholar] [CrossRef]
  68. Bennett, M.T. China’s sloping land conversion program: Institutional innovation or business as usual? Ecol. Econ. 2008, 65, 699–711. [Google Scholar] [CrossRef]
  69. Xu, J.; Yin, R.; Li, Z.; Liu, C. China’s ecological rehabilitation: Unprecedented efforts, dramatic impacts, and requisite policies. Ecol. Econ. 2006, 57, 595–607. [Google Scholar] [CrossRef]
  70. Desclee, B.; Bogaert, P.; Defourny, P. Forest change detection by statistical object-based method. Remote Sens. Environ. 2006, 102, 1–11. [Google Scholar] [CrossRef]
  71. Huang, C.; Davis, L.S.; Townshend, J.R.G. An assessment of support vector machines for land cover classification. Int. J. Remote Sens. 2002, 23, 725–749. [Google Scholar] [CrossRef]
  72. Melgani, F.; Bruzzone, L. Classification of hyperspectral remote sensing images with support vector machines. IEEE Trans. Geosci. Remote Sens. 2004, 42, 1778–1790. [Google Scholar] [CrossRef]
  73. Shao, Y.; Lunetta, R.S. Comparison of support vector machine, neural network, and cart algorithms for the land cover classification using limited training data points. ISPRS J. Photogram. Remote Sens. 2012, 70, 78–87. [Google Scholar] [CrossRef]
  74. Srivastava, P.K.; Han, D.; Rico-Ramirez, M.A.; Bray, M.; Islam, T. Selection of classification techniques for land use/land cover change investigation. Adv. Space Res. 2012, 50, 1250–1265. [Google Scholar] [CrossRef]
  75. Chang, C.-C.; Lin, C.-J. LIBSVM: A library for support vector machines. ACM Trans. Intell. Syst. Technol. 2011, 2, 27. [Google Scholar] [CrossRef]
  76. Fukunaga, K.; Hostetler, L. The estimation of the gradient of a density function, with applications in pattern recognition. IEEE Trans. Inf. Theory 1975, 21, 32–40. [Google Scholar] [CrossRef]
  77. Cheng, Y. Mean shift, mode seeking, and clustering. IEEE Trans. Pattern Anal. Mach. Intell. 1995, 17, 790–799. [Google Scholar] [CrossRef]
  78. Comaniciu, D.; Meer, P. Mean shift: A robust approach toward feature space analysis. IEEE Trans. Pattern Anal. Mach. Intell. 2002, 24, 603–619. [Google Scholar] [CrossRef]
  79. Congalton, R.G. A review of assessing the accuracy of classifications of remotely sensed data. Remote Sens. Environ. 1991, 37, 35–46. [Google Scholar] [CrossRef]
  80. Olofsson, P.; Foody, G.M.; Stehman, S.V.; Woodcock, C.E. Making better use of accuracy data in land change studies: Estimating accuracy and area and quantifying uncertainty using stratified estimation. Remote Sens. Environ. 2013, 129, 122–131. [Google Scholar] [CrossRef]
  81. Zhang, J.; Pham, T.-T.-H.; Kalacska, M.; Turner, S. Using Landsat Thematic Mapper records to map land cover change and the impacts of reforestation programmes in the borderlands of southeast Yunnan, China: 1990–2010. Int. J. Appl. Earth Obs. Geoinf. 2014, 31, 25–36. [Google Scholar] [CrossRef]
Figure 1. The map shows the location of the study area in China. The magnified map is a false-color mosaic (red-green-blue (RGB): Bands 4, 3 and 2, respectively) of Landsat Enhanced Thematic Mapper (ETM) images. The cities of Beijing, Langfang, Tianjin, Tangshan and Qinhuangdao compose the urban agglomeration.
Figure 1. The map shows the location of the study area in China. The magnified map is a false-color mosaic (red-green-blue (RGB): Bands 4, 3 and 2, respectively) of Landsat Enhanced Thematic Mapper (ETM) images. The cities of Beijing, Langfang, Tianjin, Tangshan and Qinhuangdao compose the urban agglomeration.
Ijgi 06 00059 g001
Figure 2. Flow diagram of the land cover and land use change (LCLUC) classification. IR-MAD, iteratively re-weighted multivariate alteration detection.
Figure 2. Flow diagram of the land cover and land use change (LCLUC) classification. IR-MAD, iteratively re-weighted multivariate alteration detection.
Ijgi 06 00059 g002
Figure 3. Flow diagram of the land cover change classification.
Figure 3. Flow diagram of the land cover change classification.
Ijgi 06 00059 g003
Figure 4. Details of the classifications of the sub-area in the JJT: (a) Landsat 5 TM sub-image at Path/Row 123/032 that was acquired on 7 September 1992; (b) Landsat 7 ETM+ sub-image that was acquired on 2 August 1999; (c) Landsat 5 TM sub-image that was acquired on 22 September 2009; (df) classification results for (a–c), respectively; and (g) China’s Land-Use/cover Datasets (CLUDs) in 2000.
Figure 4. Details of the classifications of the sub-area in the JJT: (a) Landsat 5 TM sub-image at Path/Row 123/032 that was acquired on 7 September 1992; (b) Landsat 7 ETM+ sub-image that was acquired on 2 August 1999; (c) Landsat 5 TM sub-image that was acquired on 22 September 2009; (df) classification results for (a–c), respectively; and (g) China’s Land-Use/cover Datasets (CLUDs) in 2000.
Ijgi 06 00059 g004
Figure 5. Objects that were segmented by the mean shift method and the refined impervious surface after synthesizing the object and SVM results: (a) ETM+ sub-image at Path 123/Row 032 that was acquired on 21 October 1999; (b) boundaries of objects that were segmented by the mean shift method (yellow lines); (c) impervious surface that was classified by the SVM classifier; (d) refined impervious surface based on the integration of the segmented objects and the SVM classification result; and (e) detailed information regarding the integration.
Figure 5. Objects that were segmented by the mean shift method and the refined impervious surface after synthesizing the object and SVM results: (a) ETM+ sub-image at Path 123/Row 032 that was acquired on 21 October 1999; (b) boundaries of objects that were segmented by the mean shift method (yellow lines); (c) impervious surface that was classified by the SVM classifier; (d) refined impervious surface based on the integration of the segmented objects and the SVM classification result; and (e) detailed information regarding the integration.
Ijgi 06 00059 g005aIjgi 06 00059 g005b
Figure 6. Changed area that was determined by comparing two images: (a) ETM+ sub-image at Path 123/Row 032 that was acquired on 2 August 1999; (b) Landsat 5 TM sub-image that was acquired on 22 September 2009; (c) chi-square image from the IR-MAD algorithm; and (d) mask of the changed area, which is shown in yellow against the background image (i.e., Figure 6a).
Figure 6. Changed area that was determined by comparing two images: (a) ETM+ sub-image at Path 123/Row 032 that was acquired on 2 August 1999; (b) Landsat 5 TM sub-image that was acquired on 22 September 2009; (c) chi-square image from the IR-MAD algorithm; and (d) mask of the changed area, which is shown in yellow against the background image (i.e., Figure 6a).
Ijgi 06 00059 g006
Figure 7. Classifications of the Beijing-Tianjin-Tangshan (BTT) region for 1990, 2000 and 2010: (a) land cover in 1990; (b) land cover in 2000; and (c) land cover in 2010. The area that is indicated by the black rectangle is shown in detail in Figure 4.
Figure 7. Classifications of the Beijing-Tianjin-Tangshan (BTT) region for 1990, 2000 and 2010: (a) land cover in 1990; (b) land cover in 2000; and (c) land cover in 2010. The area that is indicated by the black rectangle is shown in detail in Figure 4.
Ijgi 06 00059 g007
Figure 8. Percentages of different land cover areas for the three periods.
Figure 8. Percentages of different land cover areas for the three periods.
Ijgi 06 00059 g008
Table 1. TM/ETM+ images that were used in this study. The acquisition dates comprise the year, month and day. All of the images for 2000 were ETM+ images, whereas those for 1990 and 2010 were TM images. The images that were acquired on 27 August 1993 for 121/33 were used to complement areas that were influenced by cloud cover.
Table 1. TM/ETM+ images that were used in this study. The acquisition dates comprise the year, month and day. All of the images for 2000 were ETM+ images, whereas those for 1990 and 2010 were TM images. The images that were acquired on 27 August 1993 for 121/33 were used to complement areas that were influenced by cloud cover.
Path/Row121/32121/33122/32122/33123/32123/33
Acquired time09/23/199108/19/199007/30/199207/30/199209/07/199205/26/1989
08/27/1993
05/02/200006/19/200005/25/200005/12/200105/16/200005/16/2000
06/19/200009/07/200008/11/199909/01/200108/02/199907/01/1999
09/07/200010/12/200111/01/200011/01/200010/21/200012/10/2000
09/27/201009/11/201008/30/200908/30/200909/22/200909/22/2009
Table 2. Adjusted accuracy assessment of the classifications for the three periods. “Overall” is the overall accuracy, “Kappa” is the Kappa coefficient, “PA” is the producer’s accuracy and “UA” is the user’s accuracy.
Table 2. Adjusted accuracy assessment of the classifications for the three periods. “Overall” is the overall accuracy, “Kappa” is the Kappa coefficient, “PA” is the producer’s accuracy and “UA” is the user’s accuracy.
Class199020002010
PAUAPAUAPAUA
Forest87.2789.3389.3289.3385.6787.01
Grass61.2981.8262.3985.1964.7183.33
Shrub72.8784.4473.3882.9872.3684.78
Dry Farmland95.9388.8996.1090.4095.3289.60
Paddy Field71.3782.9372.8990.0069.5585.00
Bare Land69.0691.4961.4595.6569.3689.36
Impervious surface82.4383.0286.9283.9384.1780.36
Water Bodies90.0393.1891.3590.9194.0292.68
Aquatic Vegetation76.2687.5075.7587.1875.8187.50
Overall87.9489.1387.76
Kappa0.850.870.85
Table 3. Changes in impervious surface after integration based on the sub-image.
Table 3. Changes in impervious surface after integration based on the sub-image.
Identified by SVMNo ChangeAdd CorrectAdd ErrorRemove CorrectRemove Error
Pixel count59,32350,14411,263171479971182
Area (km2)53.3945.1310.141.547.201.06
Percentage (%)100.0084.5318.992.8913.481.99
Table 4. Error-adjusted areas and proportions of the classes for the three periods.
Table 4. Error-adjusted areas and proportions of the classes for the three periods.
199020002010
Area (km2)Percent (%)Area (km2)Percent (%)Area (km2)Percent (%)
Forest10,90819.6510,83219.4611,18520.03
Grass18283.2915172.7313282.38
Shrub45048.1152949.514,4357.94
Dry Farmland26,14247.1024,88244.7023,50942.09
Paddy Field31995.7627414.9217173.07
Bare Land4410.795200.937921.42
Impervious Surface48428.72637311.45946016.94
Water Bodies33576.0531915.7331315.61
Aquatic Vegetation2870.523100.562980.53
Sum55,508100.0055,660100.0055,855100.00

Share and Cite

MDPI and ACS Style

Yang, A.; Sun, G. Landsat-Based Land Cover Change in the Beijing-Tianjin-Tangshan Urban Agglomeration in 1990, 2000 and 2010. ISPRS Int. J. Geo-Inf. 2017, 6, 59. https://doi.org/10.3390/ijgi6030059

AMA Style

Yang A, Sun G. Landsat-Based Land Cover Change in the Beijing-Tianjin-Tangshan Urban Agglomeration in 1990, 2000 and 2010. ISPRS International Journal of Geo-Information. 2017; 6(3):59. https://doi.org/10.3390/ijgi6030059

Chicago/Turabian Style

Yang, Aqiang, and Guoqing Sun. 2017. "Landsat-Based Land Cover Change in the Beijing-Tianjin-Tangshan Urban Agglomeration in 1990, 2000 and 2010" ISPRS International Journal of Geo-Information 6, no. 3: 59. https://doi.org/10.3390/ijgi6030059

APA Style

Yang, A., & Sun, G. (2017). Landsat-Based Land Cover Change in the Beijing-Tianjin-Tangshan Urban Agglomeration in 1990, 2000 and 2010. ISPRS International Journal of Geo-Information, 6(3), 59. https://doi.org/10.3390/ijgi6030059

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