Next Article in Journal
The Effect of CSR Fit and CSR Authenticity on the Brand Attitude
Previous Article in Journal
Challenges Caused by Increased Use of E-Powered Personal Mobility Vehicles in European Cities
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Uncertainty Problems in Image Change Detection

1
School of Geomatics and Urban Spatial Information, Beijing University of Civil Engineering and Architecture, Beijing 100044, China
2
Department of Geography, University of Calgary, Calgary, AB T2N 1N4, Canada
3
Department of Geography, University of Wisconsin-Milwaukee, Milwaukee, WI 53211, USA
4
Key Laboratory of Environmental Change and Natural Disaster, Ministry of Education, Academy of Disaster Reduction and Emergency Management; Faculty of Geographical Science, Beijing Normal University, Beijing 100875, China
5
College of Mechanical Engineering and Automation, Fuzhou University, Fuzhou 350108, China
*
Author to whom correspondence should be addressed.
Sustainability 2020, 12(1), 274; https://0-doi-org.brum.beds.ac.uk/10.3390/su12010274
Submission received: 27 October 2019 / Revised: 19 December 2019 / Accepted: 20 December 2019 / Published: 30 December 2019

Abstract

:
Image Change Detection (ICD) methods are widely adopted to update large area land use/cover products. Uncertainty problems, however, are well known in such techniques, and a transparent assessment is necessary. In this study, a framework was proposed for evaluating binary land change utilizing remote sensing images. First, two widely adopted ICD methods were used to establish change maps. Second, binary decisions on Change (C) and Non-Change (NC) classes were reached through thresholding on change maps. Then, results were evaluated using two sampling designs: random sampling and stratified sampling. Analysis of results suggests that (1) for random sampling, with an increasing threshold on change variables, the overall accuracy increases and shows a large variance, which is highly correlated with the C omission error; and (2) comparatively, for stratified sampling, in which two strata (i.e., C and NC) were set, the overall accuracy shows a smaller variance and is highly associated with the NC commission error. The significant trends in accuracy assessments indicate the trade-offs between the C and NC classification errors in a binary decision and can present superficial or perfunctory accuracy evaluation in certain circumstances that the causes of error sources and uncertainty problems in ICD are not fully understood.

1. Introduction

Land Use/Cover Change (LUCC), also known as Land Change (LC), reflects a close interaction between the economic development and environmental biodiversity [1] and shapes our decisions on sustainability. Benefitting from the availability of remote sensing image archives, many regions have established LUCC databases [2,3,4]. To closely monitor LUCC, researchers accordingly have sped up the pace in renewing LUCC products [4,5,6,7]. For example, a new generation of National Land Cover Database (NLCD) products named NLCD 2016 contains LUCC products for nearly two decades [8]. All these products utilize Landsat images as their main data sources [5,8]. It was reported that the NLCD 2006 had the first land change product, employing image change detection (ICD) technology [2], followed by the NLCD 2011, adopting improved solutions [5]. Similarly, the Chinese national LUCC database has products for six periods (1980s, 1995, 2000, 2005, 2008, and 2010) [3]. To utilize LUCC products in different applications, it is necessary to understand the workflow and qualities of these products.
Techniques for Image Change Detection (ICD), that automatically correlate and compare sets of images taken from the same area at different times [9], are widely adopted to automate or semi-automate the updating of LUCC products. Over the past decades, numerous ICD methods, using geospatial data, have been developed [9,10,11,12,13,14,15,16]. In this study, we only focused on bi-temporal (before–after) and binary (Change (C)/Non-Change (NC)) land change detection. To facilitate our work, we divided ICD techniques into two categories: bi-temporal image analysis and image-and-map analysis, based on whether thematic land classes are used to facilitate the image analysis or not. Bi-temporal image analysis is favorable for its high efficiency and effectiveness, by implementing analysis at a pixel level without ground truth data [11]. Nevertheless, Canty [17] pointed out that if change detection is carried out at the pixel level (as opposed to using segments or objects), error (typically >5%) may corrupt or even dominate the true change signal depending on its strength. Instead of adopting bi-temporal images, “image-and-map analysis” takes the thematic maps as auxiliary data in ICD. As images are not directly compared to each other during the land change detection process, radiometric and phenology differences could be well resolved [18]. The “Cross-Correlation” algorithm, designed by Koeln and Bissonnette [19], can be used to evaluate the differences between an existing land cover map (Date 1) and a recent single-date multispectral image (Date 2) [11,16,20]. Unlike bi-temporal image analysis, pixel-based spectral analysis can contribute to spatial analysis at an object-level by segmenting heterogeneous land surface into homogenous classes first. In image-and-map analysis, ICD can be regarded as a posterior step after image classification, so the final accuracy depends on the quality of the classified image of each date [12]. Consistent with the aforesaid, it is clear that there is no single method that is accurate enough to handle all the challenges when it comes to ICD. Therefore, researchers can combine several methods to produce LUCC products. For example, the Comprehensive Change Detection Method (CCDM), a knowledge- and trajectory-based solution, was utilized to produce the 2011 NLCD. In this technical framework, various image indices had been designed and used to either compute the change magnitude of land cover or indicate the change direction of biomass. At the same time, the thematic maps in CCDM facilitate the change decisions. In this paper, we adopted ICD techniques to automate land change detection.
Uncertainty problems are inherent parts of computer modeling [21,22,23]. Studies [24,25,26,27,28] have identified common error sources in most of the LUCC studies produced thus far. The above was corroborated with a number of other groups of researchers [29,30] who further added that most of the challenges exist without well-established solutions for this task. Among these challenges in ICD, some significant ones need to be highlighted. The primary challenge exists in the simple dichotic form of land change decision. For binary land change detection, image pixels are divided into two categories, Change (C) and No Change (NC), adopting various statistical techniques such as thresholding or hypothesis testing. Such techniques are criticized for failing to consider the change features based on type, severity, and continuity [31]. Furthermore, to select a threshold is complicated as the proper value is scene-specific and experientially selecting a value is not robust [32,33]. In many cases, values were adopted arbitrarily. For example, researchers have begun their analysis by proposing a distribution of their observations and then set their thresholds on the point of one standard deviation away from expected values (mean). Another significant challenge is that “change is a rare feature over a given landscape” [30]. According to Jensen [11], in most regional projects, the amount of actual change over a 1–5 year period is probably not greater than 10% of the total area. Such land features add to our difficulties in ICD modeling and spatial sampling. Moreover, land changes are dynamic processes. Samples need to be collected in 3 dimensional space [34]. Not only the spatial heterogeneity, but also the temporal characteristic of land change scenarios should be considered for spatial sampling, as the uncertainty of the final estimation accumulates in the trinity of spatial sampling and statistical inference [34]. As the aforementioned violations of both ICD modeling and sampling can introduce uncertainties in the final evaluation, one needs to meet many criteria to avoid the pitfalls during the assessment process for land change detection [30]. In this study, to explore uncertainty problems in ICD, two categories of ICD methods (i.e., “bi-temporal image analysis” plus “image-and-map analysis”) were adopted and binary change decisions were achieved by simulating various values of thresholds. Such decisions were evaluated based on two common sampling strategies, through which the uncertainty problems in ICD were illustrated.

2. Materials and Methods

2.1. Study Area and Data

The city of Beijing, China, is located in the northern hemisphere with latitude from 39°26’ to 41°03’ and longitude from 115°25’ to 117°30’. The area of Beijing is about 16,800 km2, which is mainly composed of mountains (10,418 km2) and plains (6390 km2). The population of Beijing is approximately 21,516,000 [35]. As is the case for many metropolises in the world, Beijing has been facing fast urbanization trends. To explore the inter-annual land change patterns of this metropolis during the period 2010–2015, we downloaded Landsat images from the United States Geological Survey Website [36]. We considered several factors while selecting images, including (1) sensor selection; (2) anniversary date synchronization [37]; (3) band selection; (4) cloud cover < 20%; and (5) quality ≥ 7. Two Landsat satellites, namely, Landsat 5 and Landsat 7, were simultaneously in service in 2010. When Landsat 5 became obsolete in 2013, we decided to switch to Landsat 7 and Landsat 8. Unfortunately, because some images from Landsat 7 needed significant corrections due to the Scan Line Corrector (SLC) failure, we chose Landsat 8 with the Operational Land Imager (OLI). As a result, the bi-temporal images (2010-12-14 TM and 2014-12-25 OLI image) were selected for this pilot study. The spatial resolutions for these bands for the two sensors are both 30 meters. Besides Landsat images, a 2010 classification map was acquired from the GlobaLand30, “30 m Global Land Cover (GLC) data product” [38], which is a global dataset with level 1 land cover classification scheme based on the year 2010. The test area, a sub-area of 518 × 482 pixels from a Landsat image (Figure 1) in the eastern plain of Beijing undergoes rapid urban expansion [39] and was hence selected for a pilot study. In the test area, there are three land cover types: artificial surface, cultivated land, and water, respectively composing 49.6%, 49.6%, and 0.8% of the study area (Figure 2).

2.2. Methodology

In this study, a framework was designed to explore the uncertainty problems in ICD. It was composed of three main modules (Figure 3): (1) module 1: Input; (2) module 2: Change Detection; (3) module 3: Evaluation. In module 1, geospatial data were prepared. In module 2, two ICD methods (i.e., “bi-temporal image analysis” plus “image-and-map analysis”) were adopted to construct change maps, followed by a simulation on the values of thresholds to have different land change decisions. In module 3, accuracy assessments were implemented by first designing spatial sampling and computing error matrices.

2.2.1. Change Variables

Two change variables were established adopting “bi-temporal image analysis” and “image-and-map analysis”, respectively. Two pre-processing steps were taken to change the two raw images to “clear” images by removing scattering and absorption effects. Firstly, the digital numbers (DNs) of the two Landsat scenes were rescaled to the top of atmosphere (TOA) reflectance. Secondly, the surface (ground) reflectance free from the atmospheric effect was computed. Furthermore, spatial registration was checked visually. The original spatial accuracy of two original Landsat images was about one pixel within the test area. Then, the Normalized Difference Vegetation Index (NDVI), a standardized index of greenness (relative biomass), was chosen as the change variable for bi-temporal images. To be rendered into a color map easily, the NDVI values were scaled to a range of 0–200. By subtracting two NDVI images, one from the other (after image minus before image), the differencing NDVI (dNDVI) could be achieved.
Another change variable, the Z-score, was calculated according to the Cross-Correlation (CC) algorithm. First, the mean and standard deviation of the reflectance of each band on the Date 2 image (2014 OLI image) were calculated for each land class in Date 1 map (GlobaLand30). Then, the Z-score was calculated as the normalized distance from each pixel value to the “expected” value (mean) of that class (Equation (1)):
z   =   i ( ( B i μ ic ) / σ ic ) 2
where z is the Z-score, i is the band number, c is the land cover type, B i is the ith band of Landsat image, and μ ic and σ ic are the mean and standard deviation of the ith band of the image over the c land cover type area, respectively [19].
Once change maps are constructed, land change decisions can be achieved by adopting thresholding techniques based on the assumptions made from the change magnitudes: the higher the change magnitude is, the higher the chance of having change will be.

2.2.2. Sampling and Accuracy Assessments

To evaluate the change detection, two common sample designs (random and stratified sampling) were adopted by selecting 100 samples, respectively. The referenced data were gathered through Google Earth, and the available high-resolution scenes in the Beijing area (2010-11-11, 2010-11-19, 2014-11-22, and 2014-11-24) were chosen. Moreover, the decision on C/NC mainly was based on the usage of lands, for example, from barren soil to factories or from buildings to the desolated area. Samples were gathered by checking the approximate extent of 5 × 5 matrix areas on the Google Earth images, taking the discrepancy of the spatial resolutions between image sources into consideration. In random sampling, 12 samples were referenced as C and the left 88 were as NC. The disproportion of C/NC samples was due to the fact that most of the area remains unchanged and only a small amount of area had changed. In stratified sampling, two strata (i.e., C and NC) were set and 50 samples were randomly drawn from each stratum. Such a strategy was adopted in an earlier study by Jin et al. [5] to produce the NLCD products. Based on the above two sample designs, an accuracy assessment was carried out by calculating the binary change/no change confusion matrix [40]. Three types of accuracy could be calculated: overall accuracy, producer’s accuracy, and user’s accuracy.

3. Results

Both the spatial and frequency distributions of the two change detectors (i.e., dNDVI and Z-score/CC) are illustrated in Figure 4. Looking at the spatial distribution of dNDVI (Figure 4a), the red color denotes high dNDVI values, and the green color denotes lower dNDVI values. Based on our personal observation, some red color areas could be visually identified as the vegetation. Figure 4b shows the histogram of dNDVI. It was a two-sided bell-shaped distribution, not only denoting the change magnitude but also showing the change directions (positive and negative values of change values). A series of thresholds were derived using 12 decrement numbers, and thresholds were calculated by setting different deviations from the mean of dNDVI (Equation (2)).
T = μ + 0.25 × σ × N
where T denotes for thresholds, μ is the expected mean of dNDVI, σ is the standard deviation of dNDVI, and N ranges from 1 to 12.
Figure 4c shows the Z-scores, which are all positive. Pixels with higher Z-scores (red color) were not randomly distributed over the image area, and obvious spatial aggregation can be found in many artificial surface areas with high Z-scores. Figure 4d shows the histogram of the Z-scores. Judging from the curve of the histogram, it is one-sided with all positive values, only denoting the change magnitude and a majority of pixels falls close to zero. The values of thresholds were simulated using 12 decrement numbers: 1, 2, 3, 4, 5, 6, 10, 20, 30, 40, 50, 100.
Simulations of thresholds on two different change interpreters were carried out on two different sampling designs based on the error matrix. Five indices were calculated, namely, the overall accuracy, omission error of C and NC, commission error of C as well as NC. Judging from the evaluation results in the random sampling (Figure 5a,b), we can find that for both thee two change variables, the trends in accuracy assessments were consistent: (1) with the increase of the thresholds, the overall accuracy increased dramatically. After certain thresholds, the overall accuracy reached the highest level, which was the proportion of the NC samples in the whole samples; (2) from the perspective of the producers, the C omission error decreased significantly, while the NC omission error increased dramatically as the threshold went higher; (3) from the perspective of users, the C commission error remained quite higher than the NC type; (5) the trend of overall accuracy had a corresponding tendency with the C omission error. Looking at the number of C/NC samples (Figure 5c,d), it is clear that the higher the threshold values were, the fewer the C candidates and the more NC ones.
Looking at the accuracy assessments in stratified sampling (Figure 6a,b), for both the two change variables, the trends in accuracy assessments were also apparent: (1) with the increase of the thresholds, the overall accuracy kept a steady low level below 0.6; (2) similar to the previous experiments, the C omission error increased significantly, while the NC omission error decreased correspondingly; (3) unlike the pattern in the previous stage, the trend that the C commission error higher than the NC type was lost, instead they fluctuated at the same level; (4) the trend of overall accuracy had an excellent corresponding tendency with the NC commission error. Judging from the number of C/NC samples (Figure 6c,d), again, the higher the thresholds, the fewer C candidates and more NC ones. The sum of the numbers of the C/NC samples remained 100. In binary change detection, land change was put into the dichotic forms, so C and NC classes were complementary not only in image space but also in sample space.
Again, such significant trends in accuracy assessments can be viewed from two perspectives. From the users’ perspective, which is the prediction side or the map side, when we increased the threshold on change magnitude, the predicted C candidates decreased, and the unqualified candidate will be the NC ones, as the pool of the samples remains the same. From the producers’ perspective, the C and NC samples were always the same. With the increase of the threshold, the readjusted predicted C and NC samples were further referenced to C/NC classes.
The above evaluation helps visualize a transparent accuracy assessment on ICD. To be mentioned, in Figure 5 and Figure 6, a convergent point can be identified. The converged point for these four error (C OE/CE, NC OE/CE) lines seemed to denote a trade-off balance among different errors. For binary land change detection, the gain and loss of the C/NC sides had an effect on the accuracy assessment.

4. Discussion

Errors and uncertainties are build-in imperfections in ICD and make it hard to make accurate decision in a chaotic situation (conflicting information from conflicting sources) [41]. As Mowrer [42] said, perhaps the worst nightmare of a natural resource manager is to appear “uncertain” to the public or to admit that there is an “error” in the decision process being presented. However, noise can give information on the state of the error source [43], and biases are as crucial as the correctness in land change detection. As too much emphasis had been put on the overall accuracy of a map, without considering the causes of the results [25], in this discussion, we shift our focus from the accuracy assessment to the causes of it by answering the following three questions: (1) How was the performance of the ICD methods? (2) Was the sampling biased or not? (3) Why are trends so significant in the accuracy analysis of binary land change detection?

4.1. Bi-Temporal Image Analysis Versus Image-and-Map Analysis

Two common binary change detection methods were adopted in this study, and their performance was expected to be judged from the accuracy analysis. If we simply judge from the overall accuracy, it seems that the performance of ICD under certain thresholds can be very good. Yet, after a careful check, reporting the overall accuracy with the omission and commission errors, it would show us that both of the two ICD methods in this study do not have a satisfying performance for land change detection. For bi-temporal image analysis, assumptions of identical image situations are critical. Yet, many assumptions were not carefully checked. For example, at the stage of image prepossessing, although photometric calibrations were adopted to minimize the effects of different sensors, Landsat solar geometry was not provided per pixel and TOA adjustment within a given scene did not vary per pixel [44], so TOA correction did not enhance within-scene fidelity. For image-and-map analysis, the classification map was regarded as the baseline in ICD so that the spectral center of land classes can be calculated. Yet, classification errors exist. Meanwhile, we assume that the majority of each class remains unchanged between paired images so that the spectral statistical features of land classes remain unchanged between two dates, which is also hard to guarantee.
When remote sensing images (time shots), which provide a snapshot of the land surface, are utilized to represent land dynamics, which are more appropriately viewed as fuzzy phenomena [23], uncertainty problems definitely exist. In binary change detection, land change dynamics is interpreted either as abrupt change (the complete replacement of one cover type by another, also called “real change”) or as subtle change (the modification of land cover without changing its overall classification, also called “pseudo change”) [5,31]. Subtle change has the spectral variance caused by the phenology problem (cyclical changes in the condition of the ground cover that do not relate to class change) [27]. Detection of image differences may be confused with problems in phenology and cropping, and such problems may be exacerbated by limited image availability, poor quality in temperate zones (especially clouds during critical times in spring of rapid phonological change), and difficulties in calibrating poor images [45]. Regions with high change magnitude are visualized in Figure 4a (vegetation area). It hints that the signal of phenology change highly interweaves with the signal of real land cover change. Without breaking these two types of signals apart, it is hard to present a satisfactory work. To minimize phenology effect, it is practical to adopt two pairs of images, one pair within leaf-on season and the other within leaf-off season to mask the area of subtle change [5].

4.2. Random Sampling Versus Stratified Sampling

Samples are regarded as a subset of individuals from within a population to estimate characteristics of the whole population [34]. To construct unbiased estimates, samples should concern the selection probabilities [46]. As land change is a rare feature of a given landscape, the C samples are scarce in the random sample pool. When thresholding is adopted, these samples have a greater chance of being detected as NC than C class. By setting high thresholds, the disproportion of C/NC samples can result in high overall accuracy, yet with high C omission and commission error and low NC omission and commission error. Although random sampling is a probabilistic sampling, the spatial and temporal heterogeneity of the study area is not considered. In this respect, random sampling cannot represent land change dynamics.
The optimal selection of spatial sampling locations depends in part on the spatial structure [47]. Many studies claimed that spatial stratified homogeneity can achieve more efficient spatial sampling and inference [34]. Theoretically, to achieve unbiased sampling, the evenness in the spatial-temporal domain should be considered. From the perspective of the spatial domain, the landscape is heterogeneous with various land classes; from the perspective of the temporal domain, the proportions of C/NC area are not equal. As our interest is in C/NC rather than from–to changes, two strata (i.e., C/NC) were chosen in stratified sampling. With the increase of the value of threshold, the trend of the trade-offs between C/NC errors still remains. Or in other words, the result is still biased. In this case, it is due to the abundant C samples. But stratified sampling has better performance than random sampling, as smaller variance of overall accuracy has been identified. In this respect, stratified sampling is superior to random sampling.

4.3. Thresholding and Sensitivity Analysis

Thresholding is widely adopted in ICD, yet also widely criticized for its “hard/crisp” solutions [33]. First, it requires an arbitrary value of threshold, and there are few instances where this value can be rationally chosen. Secondly, our decisions are dichotomous. There are only two outcomes (i.e., C/NC). Because such techniques cannot quantify our degree of uncertainty, many researchers resorted to inferential statistical tools, such as hypothesis analysis, to achieve “soft” decisions [48,49,50,51]. In hypothesis analysis, a distribution for the observations is presented and then thresholds set on the point bearing certain a probability value showing the decision strength [51,52]. No matter what types of techniques are adopted for decision making, we want to be confident about our decisions. Unfortunately, it is impossible because the probability of two types of errors (i.e., Type I error and Type II error) cannot be simultaneously minimized [52]. Faced with this problem, statisticians adopt the strategy of controlling only α (level of significance), that is, to select a very low value for α, the probability of making a Type I error [52] which is tantamount to saying “better to accept something false than reject something true”. This is the same posture we had when we set the value of the thresholds. We started with a low value for the threshold, so that the probability of missing C pixels was low, as we were interested in change, and it was better for us not to set the value of the threshold high. To make a decision, costs are usually evaluated to compare the risk (change) of decisions [53]. As it is impossible for us to control C and NC errors simultaneously, not only in image space, but also sample space, it is a great challenge for us to determine the cost and risk.
Another inferential technique, sensitivity analysis (SA), a measure of the variations for a given input on a given output [54], has been emphasized as an essential part of modelling [55,56,57]. The advantage of such a technique is prominent because no assumptions are needed, compared to the technique of hypothesis analysis [58]. Sensitivity analysis allows modelers to improve and communicate the quality of their subjective beliefs about the merits of different strategies [59] and has been used in land surface change detection [4]. A practical way of SA on land change can be implemented by designing an unbiased sampling and calculating error matrix [30]. Similarly, in this study, the sensitivity of thresholds on land change decisions was evaluated, adopting the method of traditional accuracy assessments. Moreover, it is claimed that SA can identify the critical control point, because different options can be identified in such dynamic analysis, so that a decision maker can evaluate these options [59]. In this study, we set arbitrary values of thresholds to make land change decisions. As land change is rare, we intended to increase the values of the thresholds. If we have unbiased sampling, our experiments can help us to locate a critical point of setting the arbitrary value of the threshold in Figure 5 and Figure 6, and it may be useful for decision making.
To be mentioned, both the careful sensitivity analysis of model-based inference and the paradox existing in its applications were emphasized in the literature [60]. In this study, significant trends were found in accuracy assessment. Yet, after further exploration on the mechanism of binary land change detection, we concluded that the failure to represent land change dynamics, especially for spatial sampling, would lead to dysfunctional accuracy analysis. To avoid perfunctory evaluation in binary change detection, it is necessary to find the causes of the results.

Author Contributions

Conceptualization and analysis, W.W.; funding acquisition, W.F.; writing—original draft, W.W. and M.H.-B.; writing—review & editing, W.W., C.W. and W.N.; supervision, M.H.-B., C.W. and W.F. All authors have read and agreed to the published version of the manuscript.

Funding

This work is mainly supported by the National Key Research and Development Program of China (No. 2018YFC1508803), and jointly supported by the National Key Research and Development Program of China (No. 2017YFA0604903). This work is also funded by Open Fund of State Laboratory of Information Engineering in Surveying, Mapping and Remote Sensing, Wuhan University (Grant No.19E01).

Acknowledgments

The authors would like to thank anonymous reviewers for their valuable comments and language editing which have greatly improved the quality of this manuscript.

Conflicts of Interest

No potential conflict of interest was reported by the authors.

References

  1. Liu, J.; Zhang, Z.; Zhuang, D.; Wang, Y.; Zhou, W.; Zhang, S.; Li, R.; Jiang, N.; Wu, S. A study on the spatial-temporal dynamic changes of land-useand driving forces analyses of China in the 1990s. Geogr. Res. 2003, 22, 1–12. [Google Scholar] [CrossRef]
  2. Fry, J.; Xian, G.Z.; Jin, S.; Dewitz, J.; Homer, C.G.; Yang, L.; Barnes, C.A.; Herold, N.D.; Wickham, J.D. Completion of the 2006 national land cover database for the conterminous united states. Photogramm. Eng. Remote Sens. 2011, 77, 858–864. [Google Scholar]
  3. Zhang, Z.; Wang, X.; Zhao, X.; Liu, B.; Yi, L.; Zuo, L.; Wen, Q.; Liu, F.; Xu, J.; Hu, S. A 2010 update of National Land Use/Cover Database of China at 1: 100000 scale using medium spatial resolution satellite images. Remote Sens. Environ. 2014, 149, 142–154. [Google Scholar] [CrossRef]
  4. Xian, G.; Homer, C.G.; 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] [Green Version]
  5. Jin, S.; Yang, L.; Danielson, P.; Homer, C.; Fry, J.; Xian, G. A comprehensive change detection method for updating the National Land Cover Database to circa 2011. Remote Sens. Environ. 2013, 132, 159–175. [Google Scholar] [CrossRef] [Green Version]
  6. Xian, G.; Homer, C. Updating the 2001 National Land Cover Database Impervious Surface Products to 2006 using Landsat Imagery Change Detection Methods. Remote Sens. Environ. 2010, 114, 1676–1686. [Google Scholar] [CrossRef]
  7. Song, X.; Hansen, M.C.; Stehman, S.V.; Potapov, P.V.; Tyukavina, A.; Vermote, E.F.; Townshend, J.R. Global land change from 1982 to 2016. Nature 2018, 560, 639–643. [Google Scholar] [CrossRef]
  8. MRLC. Available online: http://www.mrlc.gov/nlcd (accessed on 3 December 2019).
  9. Singh, A. Review article digital change detection techniques using remotely-sensed data. Int. J. Remote Sens. 1989, 10, 989–1003. [Google Scholar] [CrossRef] [Green Version]
  10. Blaschke, T. Towards a framework for change detection based on image objects. Göttinger Geogr. Abh. 2005, 113, 1–9. [Google Scholar]
  11. Jensen, J.R. Introductory Digital Image Processing: A Remote-Sensing Perspective; Pearson Education, Inc.: Glenview, IL, USA, 2016. [Google Scholar] [CrossRef]
  12. Lu, D.; Mausel, P.; Brondizio, E.; Moran, E. Change detection techniques. Int. J. Remote Sens. 2004, 25, 2365–2401. [Google Scholar] [CrossRef]
  13. Radke, R.J.; Andra, S.; Al-Kofahi, O.; Roysam, B. Image change detection algorithms: A systematic survey. IEEE Trans. Image Process. 2005, 14, 294–307. [Google Scholar] [CrossRef] [PubMed]
  14. Ridd, M.K.; Liu, J. A Comparison of Four Algorithms for Change Detection in an Urban Environment. Remote Sens. Environ. 1998, 63, 95–100. [Google Scholar] [CrossRef]
  15. Tewkesbury, A.P.; Comber, A.J.; Tate, N.J.; Lamb, A.; Fisher, P.F. A critical synthesis of remotely sensed optical image change detection techniques. Remote Sens. Environ. 2015, 160, 1–14. [Google Scholar] [CrossRef] [Green Version]
  16. Civco, D.L.; Hurd, J.D.; Wilson, E.H.; Song, M.; Zhang, Z. A comparison of land use and land cover change detection methods. In Proceedings of the ASPRS-ACSM Annual Conference, Washington, DC, USA, 22–26 April 2002. [Google Scholar]
  17. Canty, M.J. Image Analysis, Classification and Change Detection in Remote Sensing: With algorithms for ENVI/IDL and Python; CRC Press: Florida, FL, USA, 2014. [Google Scholar]
  18. Hussain, M.; Chen, D.; Cheng, A.; Wei, H.; Stanley, D. Change detection from remotely sensed images: From pixel-based to object-based approaches. ISPRS J. Photogramm. Remote Sens. 2013, 80, 91–106. [Google Scholar] [CrossRef]
  19. Koeln, G.; Bissonnette, J. Cross-correlation analysis: Mapping landcover change with a historic landcover database and a recent, single-date multispectral image. In Proceedings of the 2000 ASPRS Annual Convention, Washington, DC, USA, 1–5 May 2000. [Google Scholar]
  20. Tarantino, C.; Adamo, M.; Lucas, R.; Blonda, P. Detection of changes in semi-natural grasslands by cross correlation analysis with WorldView-2 images and new Landsat 8 data. Remote Sens. Environ. 2016, 175, 65–72. [Google Scholar] [CrossRef] [PubMed]
  21. Hamby, D. A review of techniques for parameter sensitivity analysis of environmental models. Environ. Monit. Assess. 1994, 32, 135–154. [Google Scholar] [CrossRef]
  22. VanderWeele, T.J.; Tchetgen, E.J.T.; Halloran, M.E. Interference and sensitivity analysis. Stat. Sci. 2014, 29, 687. [Google Scholar] [CrossRef]
  23. Plewe, B. The Nature of Uncertainty in Historical Geographic Information. Trans. GIS 2002, 6, 431–456. [Google Scholar] [CrossRef]
  24. Hao, M.; Shi, W.; Deng, K.; Zhang, H.; He, P. An object-based change detection approach using uncertainty analysis for VHR images. J. Sens. 2016, 2016. [Google Scholar] [CrossRef]
  25. Patil, G.; Taillie, C. Modeling and interpreting the accuracy assessment error matrix for a doubly classified map. Environ. Ecol. Stat. 2003, 10, 357–373. [Google Scholar] [CrossRef]
  26. Zhu, Z.; Woodcock, C.E. Continuous change detection and classification of land cover using all available Landsat data. Remote Sens. Environ. 2014, 144, 152–171. [Google Scholar] [CrossRef] [Green Version]
  27. Dronova, I.; Gong, P.; Wang, L. Object-based analysis and change detection of major wetland cover types and their classification uncertainty during the low water period at Poyang Lake, China. Remote Sens. Environ. 2011, 115, 3220–3236. [Google Scholar] [CrossRef]
  28. Song, C.; Woodcock, C.E. Monitoring forest succession with multitemporal Landsat images: Factors of uncertainty. IEEE Trans. Geosci. Remote Sens. 2003, 41, 2557–2567. [Google Scholar] [CrossRef]
  29. Stehman, S.V. Estimating area from an accuracy assessment error matrix. Remote Sens. Environ. 2013, 132, 202–211. [Google Scholar] [CrossRef]
  30. Olofsson, P.; Foody, G.M.; Herold, M.; Stehman, S.V.; Woodcock, C.E.; Wulder, M.A. Good practices for estimating area and assessing accuracy of land change. Remote Sens. Environ. 2014, 148, 42–57. [Google Scholar] [CrossRef]
  31. Coppin, P.; Lambin, E.; Jonckheere, I.; Muys, B. Digital change detection methods in natural ecosystem monitoring: A review. In Analysis of Multi-Temporal Remote Sensing Images; World Scientific: Singapore, 2002; pp. 3–36. [Google Scholar]
  32. Rosin, P.L.; Ioannidis, E. Evaluation of global image thresholding for change detection. Pattern Recognit. Lett. 2003, 24, 2345–2356. [Google Scholar] [CrossRef] [Green Version]
  33. Rosin, P. Thresholding for change detection. In Proceedings of the Sixth International Conference on Computer Vision, Bombay, India, 4–7 January 1998; pp. 274–279. [Google Scholar]
  34. Wang, J.; Stein, A.; Gao, B.; Ge, Y. A review of spatial sampling. Spat. Stat. 2012, 2, 1–14. [Google Scholar] [CrossRef]
  35. Stats. Available online: http://www.stats.gov.cn/ (accessed on 5 October 2016).
  36. USGS. Available online: http://glovis.usgs.gov/ (accessed on 20 October 2016).
  37. Siewe, S.S. Change Detection Analysis of the Landuse and Landcover of the Fort Cobb Reservoir Watershed; Oklahoma State University: Stillwater, OK, USA, 2007. [Google Scholar]
  38. GlobaLand30. Available online: http://www.globallandcover.com/ (accessed on 10 June 2016).
  39. Li, H.; Song, W. Expansion of Rural Settlements on High-Quality Arable Land in Tongzhou District in Beijing, China. Sustainability 2019, 11, 5153. [Google Scholar] [CrossRef] [Green Version]
  40. Van Oort, P. Interpreting the change detection error matrix. Remote Sens. Environ. 2007, 108, 1–8. [Google Scholar] [CrossRef]
  41. Couclelis, H. The Certainty of Uncertainty: GIS and the Limits of Geographic Knowledge. Trans. GIS 2003, 7, 165–175. [Google Scholar] [CrossRef]
  42. Mowrer, H.T. Accuracy (re)assurance: Selling uncertainty assessment to the uncertain. In Spatial Accuracy Assessment: Land Information Uncertainty in Natural Resources; Jaton, K.L.A.A., Ed.; Ann Arbor Press: Chelsea, MI, USA, 1999. [Google Scholar]
  43. Thompson, S.K. Sampling, 3rd ed.; Wiley: Hoboken, HN, USA, 2012. [Google Scholar]
  44. Hansen, M.C.; Loveland, T.R. A review of large area monitoring of land cover change using Landsat data. Remote Sens. Environ. 2012, 122, 66–74. [Google Scholar] [CrossRef]
  45. Boori, M.S.; Amaro, V.E. Land use change detection for environmental management: Using multi-temporal, satellite data in the Apodi Valley of northeastern Brazil. Appl. GIS 2010, 6, 1–15. [Google Scholar] [CrossRef]
  46. Heckman, J.J. Sample selection bias as a specification error. Econom. J. Econom. Soc. 1979, 153–161. [Google Scholar] [CrossRef]
  47. Rogerson, P.A.; Delmelle, E.; Batta, R.; Akella, M.; Blatt, A.; Wilson, G. Optimal Sampling Design for Variables with Varying Spatial Importance. Geogr. Anal. 2004, 36, 177–194. [Google Scholar] [CrossRef]
  48. Biagini, F.; Campanino, M. Elements of Probability and Statistics: An. Introduction to Probability with de Finetti’s Approach and to Bayesian Statistics; Springer: New York, NY, USA, 2016. [Google Scholar]
  49. Harris, R.; Jarvis, C. Statistics for Geography and Environmental Science, 1st ed.; Routledge: London, UK, 2014. [Google Scholar] [CrossRef]
  50. Bruzzone, L.; Prieto, D.F. Automatic analysis of the difference image for unsupervised change detection. IEEE Trans. Geosci. Remote Sens. 2000, 38, 1171–1182. [Google Scholar] [CrossRef] [Green Version]
  51. Teng, S.; Chen, Y.; Cheng, K.; Lo, H. Hypothesis-test-based landcover change detection using multi-temporal satellite images–A comparative study. Adv. Space Res. 2008, 41, 1744–1754. [Google Scholar] [CrossRef]
  52. Burt, J.E.; Barber, G.M.; Rigby, D.L. Elementary Statistics for Geographers; Guilford Press: New York, NY, USA, 2009. [Google Scholar]
  53. Hwang, C.L.; Lin, M.J. Group Decision Making Under Multiple Criteria: Methods and Applications; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2012. [Google Scholar]
  54. Saltelli, A.; Tarantola, S.; Campolongo, F.; Ratto, M. Sensitivity Analysis in Practice: A Guide to Assessing Scientific Models; John Wiley & Sons, Ltd.: Hoboken, NJ, USA, 2004. [Google Scholar] [CrossRef]
  55. Seebach, L.; Strobl, P.; Vogt, P.; Mehl, W.; San-Miguel-Ayanz, J. Enhancing post-classification change detection through morphological post-processing–a sensitivity analysis. Int. J. Remote Sens. 2013, 34, 7145–7162. [Google Scholar] [CrossRef]
  56. Verburg, P.H.; Tabeau, A.; Hatna, E. Assessing spatial uncertainties of land allocation using a scenario approach and sensitivity analysis: A study for land use in Europe. J. Environ. Manag. 2013, 127, S132–S144. [Google Scholar] [CrossRef]
  57. Chen, Y.; Yu, J.; Khan, S. Spatial sensitivity analysis of multi-criteria weights in GIS-based land suitability evaluation. Environ. Model. Softw. 2010, 25, 1582–1591. [Google Scholar] [CrossRef]
  58. Iman, R.L.; Helton, J.C. An investigation of uncertainty and sensitivity analysis techniques for computer models. Risk Anal. 1988, 8, 71–90. [Google Scholar] [CrossRef]
  59. Christopher, F.H.; Patil, S.R. Identification and review of sensitivity analysis methods. Risk Anal. 2002, 22, 553–578. [Google Scholar] [CrossRef]
  60. Saltelli, A.; Annoni, P. How to avoid a perfunctory sensitivity analysis. Environ. Model. Softw. 2010, 25, 1508–1517. [Google Scholar] [CrossRef]
Figure 1. Beijing 2014-12-25 Landsat-8 image (composited using false color) and pilot study area.
Figure 1. Beijing 2014-12-25 Landsat-8 image (composited using false color) and pilot study area.
Sustainability 12 00274 g001
Figure 2. Classification map and image for the study area: (a): 30 m Global Land Cover map; (b) 2010-12-14 Landsat-5 image (composited using false color).
Figure 2. Classification map and image for the study area: (a): 30 m Global Land Cover map; (b) 2010-12-14 Landsat-5 image (composited using false color).
Sustainability 12 00274 g002
Figure 3. Framework of image-based land change detection and accuracy assessment (NDVI: the Normalized Difference Vegetation Index).
Figure 3. Framework of image-based land change detection and accuracy assessment (NDVI: the Normalized Difference Vegetation Index).
Sustainability 12 00274 g003
Figure 4. Spatial and frequency distributions of two change maps. (a) map of dNDVI; (b) histogram of dNDVI; (c) map of Z-score; (d) histogram of Z-score.
Figure 4. Spatial and frequency distributions of two change maps. (a) map of dNDVI; (b) histogram of dNDVI; (c) map of Z-score; (d) histogram of Z-score.
Sustainability 12 00274 g004
Figure 5. Trends in sensitivity analysis based on random sampling. (a) Accuracy assessments with the thresholds on dNDVI; (b) accuracy assessments with the thresholds on Z-score; (c) trends of C/NC samples (from users’ side) with the thresholds on dNDVI; (d) trends of C/NC samples (from users’ side) with the thresholds on Z-score. OA: Overall accuracy; C OE: omission error of change; C CE: commission error of change; NC OE: omission error of no change; NC CE: commission error of no change).
Figure 5. Trends in sensitivity analysis based on random sampling. (a) Accuracy assessments with the thresholds on dNDVI; (b) accuracy assessments with the thresholds on Z-score; (c) trends of C/NC samples (from users’ side) with the thresholds on dNDVI; (d) trends of C/NC samples (from users’ side) with the thresholds on Z-score. OA: Overall accuracy; C OE: omission error of change; C CE: commission error of change; NC OE: omission error of no change; NC CE: commission error of no change).
Sustainability 12 00274 g005
Figure 6. Trends in sensitivity analysis based on stratified sampling. (a) Accuracy assessments with the thresholds on dNDVI; (b) accuracy assessments with the threshold values established for the Z-score; (c) trends of C/NC samples (from users’ side) with the threshold values established for the dNDVI; (d) trends of C/NC samples (from users’ side) with the threshold values established for the Z-score.
Figure 6. Trends in sensitivity analysis based on stratified sampling. (a) Accuracy assessments with the thresholds on dNDVI; (b) accuracy assessments with the threshold values established for the Z-score; (c) trends of C/NC samples (from users’ side) with the threshold values established for the dNDVI; (d) trends of C/NC samples (from users’ side) with the threshold values established for the Z-score.
Sustainability 12 00274 g006

Share and Cite

MDPI and ACS Style

Wang, W.; Hall-Beyer, M.; Wu, C.; Fang, W.; Nsengiyumva, W. Uncertainty Problems in Image Change Detection. Sustainability 2020, 12, 274. https://0-doi-org.brum.beds.ac.uk/10.3390/su12010274

AMA Style

Wang W, Hall-Beyer M, Wu C, Fang W, Nsengiyumva W. Uncertainty Problems in Image Change Detection. Sustainability. 2020; 12(1):274. https://0-doi-org.brum.beds.ac.uk/10.3390/su12010274

Chicago/Turabian Style

Wang, Wenyu, Mryka Hall-Beyer, Changshan Wu, Weihua Fang, and Walter Nsengiyumva. 2020. "Uncertainty Problems in Image Change Detection" Sustainability 12, no. 1: 274. https://0-doi-org.brum.beds.ac.uk/10.3390/su12010274

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