Next Article in Journal
The Impact of Local Acquisition Time on the Accuracy of Microwave Surface Soil Moisture Retrievals over the Contiguous United States
Previous Article in Journal
Scale Issues Related to the Accuracy Assessment of Land Use/Land Cover Maps Produced Using Multi-Resolution Data: Comments on “The Improvement of Land Cover Classification by Thermal Remote Sensing”. Remote Sens. 2015, 7(7), 8368–8390
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Reply

Response to Johnson B.A. Scale Issues Related to the Accuracy Assessment of Land Use/Land Cover Maps Produced Using Multi-Resolution Data: Comments on “The Improvement of Land Cover Classification by Thermal Remote Sensing”. Remote Sens. 2015, 7, 8368–8390

1
Department of Geography, Ludwig Maximilian University of Munich, Munich 80333, Germany
2
Institute for Water Management, Hydrology and Hydraulic Engineering (IWHW), University of Natural Resources and Life Sciences, Vienna 1180, Austria
*
Author to whom correspondence should be addressed.
Remote Sens. 2015, 7(10), 13440-13447; https://doi.org/10.3390/rs71013440
Submission received: 24 August 2015 / Revised: 23 September 2015 / Accepted: 10 October 2015 / Published: 15 October 2015

Abstract

:
Following the suggestion made by Johnson (Johnson B.A., 2015), a polygon-based cross validation (CV) method is compared to the pixel-based CV method to classify different levels of land cover categories using a single-date Landsat 8 image and time series of Landsat TM images. Also, different variants of band combinations, with and without the thermal bands, were considered. The results demonstrate that the inclusion of thermal information into the classification process will improve the classification performance, as was already shown in our original study (Sun and Schulz, 2015). However, it is also demonstrated that the polygon-based CV method produced lower overall accuracy values when compared to the pixel-based CV method. This confirms the argument made by Johnson that a correlation of calibration and validation data due to random sampling of multi-scale data will overestimate the performance of the classifier, and independent polygon-based CV methods have to be applied instead.

Graphical Abstract

1. Introduction

We very much appreciate Johnson’s comments [1] on our recently published results [2] on the improvement by thermal remote sensing for land cover classification. In our original paper, we made use of single-date Landsat 8 and time series of Landsat 4/5 images to investigate the potential of thermal information for an improved land cover classification in the Attert Catchment in the Grand Duchy of Luxembourg [2]. The classification results were assessed by a 10-fold pixel based cross-validation (CV) method where pixel were randomly selected and the overall accuracy (OA) taken as an evaluation measure. We found that the inclusion of thermal bands can improve the accuracy of the land cover classification when added to the multispectral bands. Based on the accuracy assessment, we also concluded that the time series of thermal images alone produced similar classification results when compared to all other VIS/NIR and TIR band combinations.
In his comment, Johnson pointed out that the high accuracy data produced by the thermal images were likely caused by the overestimation of the pixel-based CV method [1]. He also recommended performing the CV at the region-of-interest (ROI)/polygon level, thereby avoiding the spatial autocorrelation between training and validation pixels given multi-resolution data.
We acknowledge the various arguments around the accuracy assessment of the multi-scale remotely sensed images [3,4,5] and the numerous data fusion methods [6,7]. While our experiments focused on the evaluation of the thermal bands compared to the other multispectral variables, the assessment problems of the random sampling procedure related to the scale issues of the resampled thermal images were neglected. In this response, we have carefully taken Johnson’s suggestion into account and evaluated all our original land cover classification with the polygon-based CV method.
In the following sections, we briefly introduce the new assessment procedure and show the classification results for both the single-date and time series of image applications. While we will briefly introduce the data used, a detailed description of the images and the classification algorithm can be found in [2].

2. Polygon-Based Cross Validation Method

Based on the suggestion of [1], we added a polygon-based CV analysis to our study. Training and validation samples were collected using the area of interest tool of the drawing toolbox in the ERDAS software (ERDAS, Inc., Atlanta, GA, USA); the polygons here are the same as in the original paper [2]. However, in order to make sure the pixels from the polygon as pure as possible, here the pixels were kept only when their central points are inside the polygon; the polygons were then refined. This leads to the slightly smaller sample sizes compared to the original paper. As we are not aiming at evaluating the effects of the training size on the classification accuracy, here the comparison between the pixel-based CV method and the polygon-based CV method was conducted using the same size of samples (pixels with central point inside the polygon).
For all land cover categories (except for water bodies) more than 10 polygons of ground truth data exists and there is little difference among the size of each polygon. For the 10-fold polygon-based cross validation, all polygons were split into 10 smaller sets and reorganized into 10 groups of data. The classification model was trained using pixels from nine groups of polygons (training data); the resulting model is then validated on the remaining group (validation data), and this procedure is repeated 10-times so that each group is used for validation once. In this way the correlation of multi-resolutions calibration and validation is avoided.

3. Results and Discussion

In this section, the accuracy statistics for a land cover classification applying both the original pixel based CV method and the polygon-based CV method are compared. Input data (a single-date Landsat 8 image on 21 July 2013, as well as time series of images from Landsat 4/5) and classification methods (Random Forest and the k-NN algorithms) are kept as in the original paper [2].

3.1. Three-Level Classification Based on the Single-Date Landsat 8 Image (S1)

In a first part, three variants of spectral band combinations of the single Landsat 8 image are used as input data and include: (i) Bands4, only considering bands 2 to 5 in the VIS/NIR spectral region, without thermal bands; (ii) Bands6T, as Bands4, but adding the two thermal bands 10 and 11; and (iii) Bands10T, including all bands, except the panchromatic band.
Table 1 summarizes the accuracy evaluation results for both the polygon-based and the pixel-based CV methods. The OA values for the classification of Level 1 land cover category obtained by the polygon-based CV are almost identical to the pixel-based evaluation for the three variants (97% to 98%). The OA values of Level 2 and Level 3 categories decreased about 5% and 9%, respectively, when switching to the polygon-based CV method. However, the increase in performance when adding the thermal bands into the classification is still pronounced and our conclusions drawn in the original paper still hold.
Table 1. The mean values of overall accuracy (OA) calculated by a polygon-based 10-fold cross validation (CV) method for the three variants from Landsat 8 in 2013 classified by k-NN and Random Forest: Bands4, only considering bands 2 to 5 in the visible and near infrared spectral region without thermal bands; Bands6T, as Bands4, but adding the two thermal bands 10 and 11, and Bands10T, including all bands, except the panchromatic band (see [2] for a detailed channel/band description); k-NN5 and RF represents the nearest neighbor method with k = 5 and Random Forest, respectively.
Table 1. The mean values of overall accuracy (OA) calculated by a polygon-based 10-fold cross validation (CV) method for the three variants from Landsat 8 in 2013 classified by k-NN and Random Forest: Bands4, only considering bands 2 to 5 in the visible and near infrared spectral region without thermal bands; Bands6T, as Bands4, but adding the two thermal bands 10 and 11, and Bands10T, including all bands, except the panchromatic band (see [2] for a detailed channel/band description); k-NN5 and RF represents the nearest neighbor method with k = 5 and Random Forest, respectively.
Image Classification Accuracy in 2013Level 1 (%)Level 2 (%)Level 3 (%)
Assessed by Polygon-Based CVk-NN5RFk-NN5RFk-NN5RF
Bands497.097.583.382.567.567.7
Bands6T96.596.788.788.879.179.2
Bands10T96.997.688.087.580.878.0
Assessed by Pixel-Based CVk-NN5RFk-NN5RFk-NN5RF
Bands498.498.586.287.276.578.2
Bands6T98.298.492.192.987.788.1
Bands10T98.798.693.593.589.389.2

3.2. Two-Level Classification Based on Time Series of Images (TS1 and TS2)

The second part of the analysis focus on the differences between polygon-based and pixel-based CV methods when using time series of Landsat 4/5 images as input in the land cover classification system (Table 2 and Table 3). The time series of images consist of two groups: TS1, including 7 images between 1984 and 1990, and TS2, including 6 images from 2006 to 2011. Five variants of the times series were analyzed including the following different band combinations: B3B4, the combinations of band 3 and band 4 of Landsat TM; PC3, the first three principal components of all VIS/NIR bands; 6Bands, all bands except the thermal band; Thermal only the single thermal band; and 7Bands, the combination of all 7 Landsat 4/5 bands.
Table 2. Overall Accuracy of Level 1 classification assessed by the polygon-based CV method using five variants of time series of images: B3B4, the combination of band 3 and band 4; 3PC, the first three principal components of all bands; 6Bands, the combination of six bands; Thermal, the thermal band; and 7Bands, the combination of all seven bands.
Table 2. Overall Accuracy of Level 1 classification assessed by the polygon-based CV method using five variants of time series of images: B3B4, the combination of band 3 and band 4; 3PC, the first three principal components of all bands; 6Bands, the combination of six bands; Thermal, the thermal band; and 7Bands, the combination of all seven bands.
Image NumberB3B4 (%)3PC (%)6Bands (%)Thermal (%)7Bands (%)
k-NNRFk-NNRFk-NNRFk-NNRFk-NNRF
TS1-186.687.791.691.286.687.787.890.596.295.5
TS1-291.990.994.894.495.696.288.488.896.597.0
TS1-394.294.995.995.896.496.188.490.996.896.3
TS1-495.696.296.896.796.497.188.787.096.896.7
TS1-596.196.296.896.297.397.091.791.396.897.3
TS1-696.195.997.296.396.897.491.789.396.997.1
TS1-796.096.697.497.297.797.693.294.697.797.0
TS2-191.491.595.694.295.895.580.086.095.995.1
TS2-295.195.997.296.997.597.084.589.397.897.0
TS2-397.197.298.097.698.497.892.692.798.398.2
TS2-497.097.897.897.997.898.295.092.698.798.0
TS2-598.197.998.798.398.099.094.294.498.398.3
TS2-698.598.298.698.198.698.495.495.498.398.4
Table 3. Overall Accuracy of Level 2 classification assessed by the polygon-based CV method using five variants of time series of images: B3B4, the combination of band 3 and band 4; 3PC, the first three principal components of all bands; 6Bands, the combination of six bands; Thermal, the thermal band; and 7Bands, the combination of all seven bands.
Table 3. Overall Accuracy of Level 2 classification assessed by the polygon-based CV method using five variants of time series of images: B3B4, the combination of band 3 and band 4; 3PC, the first three principal components of all bands; 6Bands, the combination of six bands; Thermal, the thermal band; and 7Bands, the combination of all seven bands.
Image NumberB3B4 (%)3PC (%)6Bands (%)Thermal (%)7Bands (%)
k-NNRFk-NNRFk-NNRFk-NNRFk-NNRF
TS1-163.362.073.372.279.480.339.251.679.779.4
TS1-272.673.378.777.981.582.050.452.980.481.8
TS1-377.976.979.381.983.482.855.457.981.983.3
TS1-480.881.585.085.586.184.262.260.781.783.2
TS1-582.783.784.184.584.686.166.366.383.385.2
TS1-683.283.486.085.985.284.766.868.885.186.2
TS1-782.983.084.485.085.186.369.972.085.486.6
TS2-168.869.683.782.485.786.042.048.783.684.5
TS2-279.480.388.489.890.190.957.358.590.689.8
TS2-386.886.289.992.488.791.060.963.991.291.0
TS2-487.889.690.691.691.891.565.764.992.791.6
TS2-587.790.192.092.691.192.370.371.592.293.0
TS2-689.191.392.192.390.291.974.570.492.093.3
The best OA value from the pixel-based CV is almost the same as in the original paper, the details can be found in the Supplementary materials (Tables S1 and S2). For simplification, we here only list the OA values of the polygon-based CV for Level 1 and Level 2 land cover categories classification. Table 2 shows OA values of the polygon-based CV for the classification of Level 1 land cover categories, whereby images within each group (TS1, TS2) are added subsequently as input. The classification performance increases with increasing numbers of images for all variants, a behavior that is already shown in our original study using a random pixel-based CV. Also, the 7Bands variant including the thermal band still achieved the best overall performance, especially for smaller number of images included. Using the full set of available images all variants performed almost equally well, with OA values of 96%–98.5%. When only using the thermal band, classification performance is reduced by 4.5% compared to the pixel-based CV in our original paper, nicely demonstrating the overestimation of performance when correlation of multi-resolution calibration and validation data are existent, as commented by Johnson [1].
The differences between the pixel-based based and a polygon-based 10-fold CV method are summarized in Figure 1. Here, OA values of all five variants from TS1 are summarized in a single box-plot for each time-step in the left image of Figure 1. It is clearly seen that on average the polygon-based CV method produced significant lower OA for both time series (TS1, TS2), again supporting the issue raised by Johnson’s comment. Two variants of 6Bands and 7Bands were selected to show the detailed variation for the two methods in the right image of Figure 1. Besides the lower OA in comparison to the pixel-based CV method, the polygon-based CV still produced higher average OA for the 7Bands compared to the 6Bands without thermal band (the right image in Figure 1).
Repeating this analysis for the classification of Level 2 land cover categories, the differences in the performance measure (OA) between both CV methods is even more pronounced. Table 3 summarizes the classification results. The average OA values are in general lower, as we analyze more specific land cover categories. The best OA for TS1 and TS2 are 86.6% and 93.3% when including the full set of images. The 7Bands variant including the thermal band still achieved the best OA value of 86.6% for TS1, which is 10.9% lower than 97.5% from the pixel-based CV method. The best OA of TS2 from polygon-based CV method is about 5% lower than the corresponding value for the pixel-based CV method. Again, the Thermal variant including only the single thermal band, showed only a relatively weak performance with OA values of 72% and 74.5% for the both time series, compared to 96.8% 95.8% for the pixel-based CV method. Figure 2, taking TS2 as an example, summarizes the differences between the pixel-based and the polygon-based CV methods in the left image and displays the variation of the selected 6Bands and 7Bands variants, again supporting the issue raised by [1].
Figure 1. The distribution of OA values for the Level 1 land cover category classification using times series TS1, the polygon-based and pixel-based 10-fold cross validation methods and the Random Forest methods. (Left): all variants (Table 2) are summarized in a single box-plot; (Right): OA comparison of selected variants of 6Bands without thermal band and 7Bands with thermal band.
Figure 1. The distribution of OA values for the Level 1 land cover category classification using times series TS1, the polygon-based and pixel-based 10-fold cross validation methods and the Random Forest methods. (Left): all variants (Table 2) are summarized in a single box-plot; (Right): OA comparison of selected variants of 6Bands without thermal band and 7Bands with thermal band.
Remotesensing 07 13440 g001
Figure 2. The distribution of OA values for the Level 2 land cover category classification using times series TS2, the polygon-based and pixel-based 10-fold cross validation methods and the Random Forest methods. (Left): all variants (Table 3) are summarized in a single box-plot; (Right): OA comparison of selected variants of 6Bands without thermal band and 7Bands with thermal band.
Figure 2. The distribution of OA values for the Level 2 land cover category classification using times series TS2, the polygon-based and pixel-based 10-fold cross validation methods and the Random Forest methods. (Left): all variants (Table 3) are summarized in a single box-plot; (Right): OA comparison of selected variants of 6Bands without thermal band and 7Bands with thermal band.
Remotesensing 07 13440 g002

4. Conclusions

In this response to [1], a polygon-based CV method, as suggested by [1], was applied to evaluate a land cover classification for three different levels of land cover categories. The classification was based on (i) a single-date Landsat 8 image; and (ii) time series of Landsat 4/5 images. The performance of classification results using the polygon-based CV were compared to a pixel-based CV method as applied in our original application [2].
For the single-date Landsat 8 image, the polygon-based method achieved almost similar accuracy values when compared to the pixel-based method, for all three-levels of land cover categories and for both classification methods used. When using time series of images, five different variants of band combinations (see Table 2) with and without thermal information have been considered.
The accuracy of the Level 1 classification decreased but to a very acceptable and still useful level when compared to the commonly recommended standard of 85% [8] (the best OA of Thermal is 94.6% and the best OA of 7Bands is 98.4%). The most obvious decline in performance is observed in the classification results for the Level 2 category, of which the best OA among the five variants is 93.3% and only 74.5% when using only the thermal bands.
Consistent with our former findings, the inclusion of the thermal bands still improved the land cover classification in comparison to only using the VIS/NIR bands, also when assessing the classification results with a polygon type of CV approach. This has also been shown by other researcher, Eisavi et al. [9] applied the random forest classifier to the multi-temporal spectral and thermal features in land cover classification and found that the contribution of multi-temporal thermal information led to a considerable increase in the accuracy data. When using time series of thermal images to classify land cover at the Level 2 category, the performance and OA values were significantly lower for the polygon based CV when compared to the pixel-based evaluation for all band combinations considered. Again, the inclusion of thermal information improved the classification results on various levels.
In summary, a clear effect of correlation in the samples for calibration and validation due to multi-resolution data could be observed here. Classification accuracy (OA) was highly overestimated when ignoring correlation effects in the selection of calibration and validation data using time series of images as input. We therefore strongly support the comments made by Johnson and can only support his recommendation concerning the appropriate choice of CV methods. However, the choice of the CV method did not change our original conclusions, in that the inclusion of thermal data into the classification process, can significantly improve classification results.

Supplementary Files

Supplementary File 1

Acknowledgments

The authors would like to thank Brian A. Johnson for his comment on our original paper which initialized a very intensive and still ongoing discussion about the effects of correlated multi-scale data within calibration-validation exercises.
We would also like to sincerely thank the German Research Foundation (DFG) for funding this research through the CAOS (Catchments as Organized Systems) Research Unit (FOR 1598, SCHU1271/5-1), and the China Scholarship Committee for the support of the research and USGS for providing all the Landsat images.

Author Contributions

Liya Sun and Karsten Schulz designed this experiment and analyzed the data together. The results were compared and summarized in this manuscript by Liya Sun. Karsten Schulz has spent the most effort in revising the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Johnson, B.A. Scale issues related to the accuracy assessment of land use/land cover maps produced using multi-resolution data: Comments on “the improvement of land cover classification by thermal remote sensing”. Remote Sens. 2015, 7, 8368–8390. Remote Sens. 2015, 7, 13436–13439. [Google Scholar]
  2. Sun, L.; Schulz, K. The improvement of land cover classification by thermal remote sensing. Remote Sens. 2015, 7, 8368–8390. [Google Scholar] [CrossRef]
  3. Congalton, R.G.; Green, K. Assessing the Accuracy of Remotely Sensed Data: Principles and Practices; CRC Press/Taylor & Francis: Boca Raton, FL, USA, 2009. [Google Scholar]
  4. 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]
  5. Foody, G.M. Status of land cover classification accuracy assessment. Remote Sens. Environ. 2002, 80, 185–201. [Google Scholar] [CrossRef]
  6. Colditz, R.R.; Wehrmann, T.; Bachmann, M.; Steinnocher, K.; Schmidt, M.; Strunz, G.; Dech, S. Influence of image fusion approaches on classification accuracy: A case study. Int. J. Remote Sens. 2006, 27, 3311–3335. [Google Scholar] [CrossRef]
  7. Zhu, X.; Chen, J.; Gao, F.; Chen, X.; Masek, J.G. An enhanced spatial and temporal adaptive reflectance fusion model for complex heterogeneous regions. Remote Sens. Environ. 2010, 114, 2610–2623. [Google Scholar] [CrossRef]
  8. Anderson, J.R. A Land Use and Land Cover Classification System for Use with Remote Sensor Data; U.S. Government Printing Office (GPO): Washington, DC, USA, 1976. [Google Scholar]
  9. Eisavi, V.; Homayouni, S.; Yazdi, A.; Alimohammadi, A. Land cover mapping based on random forest classification of multitemporal spectral and thermal images. Environ. Monit. Assess. 2015, 187, 1–14. [Google Scholar] [CrossRef] [PubMed]

Share and Cite

MDPI and ACS Style

Sun, L.; Schulz, K. Response to Johnson B.A. Scale Issues Related to the Accuracy Assessment of Land Use/Land Cover Maps Produced Using Multi-Resolution Data: Comments on “The Improvement of Land Cover Classification by Thermal Remote Sensing”. Remote Sens. 2015, 7, 8368–8390. Remote Sens. 2015, 7, 13440-13447. https://doi.org/10.3390/rs71013440

AMA Style

Sun L, Schulz K. Response to Johnson B.A. Scale Issues Related to the Accuracy Assessment of Land Use/Land Cover Maps Produced Using Multi-Resolution Data: Comments on “The Improvement of Land Cover Classification by Thermal Remote Sensing”. Remote Sens. 2015, 7, 8368–8390. Remote Sensing. 2015; 7(10):13440-13447. https://doi.org/10.3390/rs71013440

Chicago/Turabian Style

Sun, Liya, and Karsten Schulz. 2015. "Response to Johnson B.A. Scale Issues Related to the Accuracy Assessment of Land Use/Land Cover Maps Produced Using Multi-Resolution Data: Comments on “The Improvement of Land Cover Classification by Thermal Remote Sensing”. Remote Sens. 2015, 7, 8368–8390" Remote Sensing 7, no. 10: 13440-13447. https://doi.org/10.3390/rs71013440

APA Style

Sun, L., & Schulz, K. (2015). Response to Johnson B.A. Scale Issues Related to the Accuracy Assessment of Land Use/Land Cover Maps Produced Using Multi-Resolution Data: Comments on “The Improvement of Land Cover Classification by Thermal Remote Sensing”. Remote Sens. 2015, 7, 8368–8390. Remote Sensing, 7(10), 13440-13447. https://doi.org/10.3390/rs71013440

Article Metrics

Back to TopTop