Skip to main content
ORNL DAAC HomeNASA Home

DAAC Home > Get Data > NASA Projects > Arctic-Boreal Vulnerability Experiment (ABoVE) > User guide

ABoVE: Landsat-derived Annual Dominant Land Cover Across ABoVE Core Domain, 1984-2014

Documentation Revision Date: 2020-11-09

Dataset Version: 1

Summary

This dataset provides two 30-m resolution time series products of annual land cover classifications over the Arctic Boreal Vulnerability Experiment (ABoVE) core domain for each year of the period 1984-2014. The data are the annual dominant plant functional type in a given 30-m pixel derived from Landsat surface reflectance, landcover training data mapped across the ABoVE domain (using Random Forests modeling, with clustering and interpretation of field photography) and very high resolution imagery to assign land cover classifications. One product has a 15-class land cover classification that breaks out forest and shrub types into several additional classes; the other product provides a simplified, 10-class approach. Classification accuracy assessment results are provided per year. Assessments were based on a probability-based random sample of reference data that supported statistically robust estimation of areas and uncertainties in mapped areas.

To facilitate analyses, the Landsat 5 Thematic Mapper and Landsat 7 Enhanced Thematic Mapper data were reprojected into the ABoVE standard Alber’s equal area projection (30-m grid) and were gridded into 180 x 180 km ABoVE standard B grid tiles (Loboda et al., 2019). This resulted in a total of 175 B grid tiles over the ABoVE core domain. Land cover classifications (10- and 15-class products) are provided for each year from 1984 to 2014 for each tile.

This dataset provides 350 files with landcover data in GeoTIFF (.tif) format which includes 175 files for the 15-class product and 175 files for the 10-class product. Each file is a 31-band GeoTIFF, with each band providing one year of data for the 1984-2014 period, where band 1 is year 1984, band 2 is 1985, and so on until band 31, which is 2014. The accuracy assessment data are provided in two comma-separated (.csv) data files.

Figure 1. Land cover (10-class product) for the year 1984 in the ABoVE domain for analysis tile BH03V04. The pixel resolution is 30 m and a tile is 180 x 180 km. Refer to Figure 3 for the landcover in the same tile for the year 2014. (from ABoVE_LandCover_BH03V04.tif)

Citation

Wang, J.A., D. Sulla-Menashe, C.E. Woodcock, O. Sonnentag, R.F. Keeling, and M.A. Friedl. 2019. ABoVE: Landsat-derived Annual Dominant Land Cover Across ABoVE Core Domain, 1984-2014. ORNL DAAC, Oak Ridge, Tennessee, USA. https://doi.org/10.3334/ORNLDAAC/1691

Table of Contents

  1. Dataset Overview
  2. Data Characteristics
  3. Application and Derivation
  4. Quality Assessment
  5. Data Acquisition, Materials, and Methods
  6. Data Access
  7. References

Dataset Overview

This dataset provides two 30-m resolution time series products of annual land cover classifications over the Arctic Boreal Vulnerability Experiment (ABoVE) core domain for each year of the period 1984-2014. The data are the annual dominant plant functional type in a given 30-meter pixel derived from time series of Landsat surface reflectance, landcover training data mapped across the ABoVE domain using Random Forests modeling, with clustering and interpretation of field photography and very high resolution imagery to assign land cover classifications. One product has a 15-class land cover classification that breaks out forest and shrub types into several additional classes; the other product provides a simplified, 10-class approach. Classification accuracy assessment results are provided per year. Assessments were based on a probability-based random sample of reference data that supported statistically robust estimation of areas and uncertainties in mapped areas.

To facilitate analyses, the Landsat 5 Thematic Mapper and Landsat 7 Enhanced Thematic Mapper data were reprojected into the ABoVE standard Alber’s equal area projection (30-m grid) and were gridded into 180 x 180 km ABoVE standard B grid tiles (Loboda et al., 2019). This resulted in a total of 175 B grid tiles over the ABoVE core domain. Land cover classifications (10- and 15-class products) are provided for each year from 1984 to 2014 for each tile.

Project: Arctic-Boreal Vulnerability Experiment

The Arctic-Boreal Vulnerability Experiment (ABoVE) is a NASA Terrestrial Ecology Program field campaign taking place in Alaska and western Canada between 2016 and 2021. Climate change in the Arctic and Boreal region is unfolding faster than anywhere else on Earth, resulting in reduced Arctic sea ice, thawing of permafrost soils, decomposition of long-frozen organic matter, widespread changes to lakes, rivers, coastlines, and alterations of ecosystem structure and function. ABoVE seeks a better understanding of the vulnerability and resilience of ecosystems and society to this changing environment.

Related Publication

Wang, J. A., Sulla-Menashe, D., Woodcock, C. E., Sonnentag, O., Keeling, R. F., & Friedl, M. A. 2020. Extensive Land Cover Change Across Arctic-Boreal Northwestern North America from Disturbance and Climate Forcing. Global Change Biology 26, 807–822. https://doi.org/10.1111/gcb.14804

Related Dataset

Loboda, T.V., E.E. Hoy, and M.L. Carroll. 2019. ABoVE: Study Domain and Standard Reference Grids, Version 2. ORNL DAAC, Oak Ridge, Tennessee, USA. https://doi.org/10.3334/ORNLDAAC/1527

Acknowledgements

This research was funded by the NASA Terrestrial Ecology- Arctic and Boreal Vulnerability Experiment, grant numbers NNX17AE74G and NNX15AU63A and by the National Science Foundation (NSF) under grant number DGE-1247312..

Data Characteristics

Spatial Coverage: Alaska and Canada

ABoVE Reference Locations

Domain: Core ABoVE

The Landsat 5 Thematic Mapper and Landsat 7 Enhanced Thematic Mapper data were reprojected into the ABoVE standard Alber’s equal area projection (30-m grid) and were gridded into 180 x 180 km ABoVE standard B grid tiles (Loboda et al., 2019). This resulted in a total of 175 B grid tiles over the ABoVE core domain.

The two landcover classifications are provided for each tile for a total of 350 data files. The respective B grid tiles are identified in the filenames.

Spatial Resolution: 30 m

Temporal Coverage: 1984-01-01 to 2014-12-31

Study Area: (all latitudes and longitudes given in decimal degrees)

Site (Region) Westernmost Longitude Easternmost Longitude Northernmost Latitude Southernmost Latitude
ABoVE domain: Alaska and Canada -170.0058 -98.97401389 76.225675 50.25900556

 

Data File Information

This dataset provides 350 files with landcover data in GeoTIFF (.tif) format which includes 175 files for the 15-class product and 175 files for the 10-class product. Each file is a 31-band GeoTIFF, with each band providing one year of data for the 1984-2014 period, where band 1 is year 1984, band 2 is 1985, and so on until band 31, which is 2014. The data values in each band represent the annual dominant plant functional type class (1-15 or 1-10; definitions in Table 2 & 3 below) in each 30-m pixel. The respective B grid tiles are identified in the filenames.

Land cover class accuracy assessment data are provided in two comma-separated (.csv) data files.

User Note:  Using the ABoVE Standard Reference Grid System

As described, the Landsat data were reprojected into the ABoVE standard Alber’s equal area projection (30-m grid) and were then gridded into 180 x 180 km ABoVE standard B grid tiles. Each of the resulting data filenames includes the respective B grid tile coordinates (BhXXvZZ).

If a user wants to select landcover data files (B grid tiles) based on a defined bounding box, they can identify the respective ABoVE standard B grid tiles based on location by using the shapefiles provided in the related dataset (Loboda et al., 2019; https://doi.org/10.3334/ORNLDAAC/1527).

Table 1. File names and descriptions

File names Descriptions
ABoVE_LandCover_BhXXvZZ.tif

Annual landcover classification of 15 plant functional types. The data provide the dominant plant functional type class for a given 30-m pixel in the identified ABoVE B grid tile (BhXXvZZ).

Where XX is the horizontal ABoVE tile location and ZZ is the vertical ABoVE tile location for the respective B grid.

ABoVE_LandCover_Simplified_BhXXvZZ.tif

Annual landcover classification of 10 plant functional types. The data provide the dominant plant functional type class for a given 30-m pixel in the identified ABoVE grid tile (BhXXvZZ).

Where XX is the horizontal ABoVE tile location and ZZ is the vertical ABoVE tile location for the respective B grid.

accuracy_assess_1984-2014.csv Classification accuracy assessment results (for the simplified 10-class version) for each year provided as error matrices. The columns indicate reference data and the rows indicate the classification. Refer to Table 4 for example data.
accuracy_summary_1984-2014.csv Provides the producer's accuracy (inverse of error rate of omission), the user's accuracy (inverse of error rate of commission) for each class, as well as the overall accuracy for each year.

GeoTIFF Data Descriptions

Two different land cover classifications are provided to describe the dominant plant functional type at a given pixel: 1) a classification with 15 different classes that breaks out forest and shrub types into several additional classes, and 2) a simplified approach with 10 different classes where forest and shrub types are aggregated.

File properties:

Bands: 31

No Data Value: 255

EPSG: 102001

Resolution: 30 m

Number rows x columns: 6,000 x 6,000

Missing values primarily occur as part of the ocean mask or for screened permanent snow/ice pixels.

Table 2. Land cover classes in the 15-class approach

ID Label Description
1 Evergreen Forest Area dominated by tall woody vegetation (> 3m tall) and over 60% canopy coverage with primarily (>75%) evergreen phenological habit (canopy maintains green foliage year-round).
2 Deciduous Forest Area dominated by tall woody vegetation (> 3m tall) and over 60% canopy coverage with primarily (>75%) deciduous phenological habit (annual cycle of leaf-on and leaf-off periods).
3 Mixed Forest Area dominated by tall woody vegetation (> 3m tall) and over 60% canopy coverage with neither forest type (deciduous or evergreen) exceeding over 60% of the area.
4 Woodland Area dominated by tall woody vegetation (> 3m tall) with between 30-60% canopy coverage. Frequently co-exists with peatlands and typically, but not always, evergreen in phenological habit.
5 Low Shrub Area dominated by dense hemi-prostrate to low-erect shrubs (5-30cm in height) with >60% area coverage. Analogous to "prostrate dwarf-shrub" and primarily occurring in tundra areas.
6 Tall Shrub Area dominated by woody vegetation between 50cm and 3m tall and shrub canopy coverage >60% coverage. Typically, but not always, deciduous phenological habit.
7 Open Shrubs Area with woody vegetation less than 3m tall and between 30-60% canopy coverage. Shrubs typically underlain by herbaceous or barren land cover.
8 Herbaceous Area dominated by herbaceous land cover greater than 60% land cover and tree/shrub cover less than 10%. 
9 Tussock Tundra Tundra-specific herbaceous land dominated by Eriophorum vaginatum and other tussock-forming herbaceous species, coverage over 60%.
10 Sparsely Vegetated 10-30% canopy coverage, any vegetation but typically herbaceous/bryophyte, with rock underneath
11 Fen Hydrologically connected, sedge/grass dominated wetland
12 Bog Ombrotrophic, peat and shrub dominated wetland
13 Shallows/littoral Lakes <1m deep with some vegetation or shoreline mixed with water/land
14 Barren <10% vegetation, mostly rock
15 Water Oceans, lakes, and rivers, either salt-water or freshwater. 

 

Table 3. Land cover classes in the 10-class approach

ID Label Combines the following classes from the 15-class approach Combines the following labels:
1 Evergreen Forest 1,4 Evergreen Forest, Woodland
2 Deciduous Forest 2,3 Deciduous Forest, Mixed Forest
3 Shrubland 5,6,7 Low Shrub, Tall Shrub, Open Shrub
4 Herbaceous 8,9 Herbaceous, Tussock Tundra
5 Sparsely Vegetated 10 Sparsely Vegetated
6 Barren 14 Barren
7 Fen 11 Fen
8 Bog 12 Bog
9 Shallows/Littoral 13 Shallows/Littoral
10 Water 15 Water

 

Accuracy assessment

File: accuracy_assess_1984-2014.csv

This file provides error matrices for classification accuracy assessment results by year, for the years 1984-2014 for the 10-class products. The results for the year 2014 are provided as an example. Reference data are indicated in columns; mapped classification data are in indicated rows. In the error matrix, correct classifications are shown along the diagonal and erroneous classifications are shown in the other cells. Refer to Table 3 for the landcover definitions.

Table 4. 10-class classification accuracy assessment for the year 2014.

year landcover evergreen_forest deciduous_forest shrub herb sparse_veg barren fen bog shallows water
2014 Everg. F 223 9 19 9 5 9 10 2 7 3
2014 Decid. F 5 119 7 1 0 0 4 0 2 0
2014 Shrub 9 14 161 20 2 2 11 0 0 2
2014 Herb 0 1 11 93 4 1 3 0 1 3
2014 Sparse Veg 2 0 10 16 108 1 2 2 1 0
2014 Barren 0 2 3 6 8 52 2 4 1 4
2014 Fen 4 9 18 6 2 0 61 1 2 1
2014 Bog 5 0 0 0 0 0 2 39 0 0
2014 Shallows 2 0 0 4 3 0 1 0 30 2
2014 Water 1 1 4 1 4 2 0 0 9 69

 

File: accuracy_summary_1984-2014.csv

This file provides the producer’s classification accuracy (inverse of error rate of omission), user’s classification accuracy (inverse of error rate of commission), and overall classification accuracy for each class in the 10-class products for each year.

Table 5. Example data for the years 2014 and 2013.

year accuracy evergreen_forest deciduous_forest shrub herb sparse_veg barren fen bog shallows water
2014 prodAcc 0.8884 0.7677 0.691 0.5962 0.7941 0.7761 0.6354 0.8125 0.566 0.8214
2014 userAcc 0.7534 0.8623 0.7285 0.7949 0.7606 0.6341 0.5865 0.8478 0.7143 0.7582
2014 overAcc 0.7467 0.7467 0.7467 0.7467 0.7467 0.7467 0.7467 0.7467 0.7467 0.7467
2013 prodAcc 0.8889 0.8125 0.7161 0.6562 0.8519 0.9155 0.6771 0.8776 0.6481 0.9059
2013 userAcc 0.8116 0.8904 0.7445 0.8268 0.7468 0.8125 0.619 0.86 0.7778 0.875
2013 overAcc 0.792 0.792 0.792 0.792 0.792 0.792 0.792 0.792 0.792 0.792

Application and Derivation

The data provide information on plant functional types in Arctic and Boreal biomes over a 31-year period which could help to quantify the loss, regrowth, and transitions between evergreen and deciduous tree species in the boreal forest biome at large spatial scales. See Wang et al (2019) for analysis details.

Quality Assessment

Two files are included that characterize the year by year accuracy for the simplified land cover classes across all years.

Following established practices for the assessment of remotely sensed land cover maps, the accuracy of the simplified land cover was assessed over each year by drawing a stratified random sample and creating a reference time series of land cover data for each site by visually interpreting very high-resolution imagery, field photography (geolocated photos provided by Kevin Turner (for the Old Crow Flats), Isla Myers-Smith (for Qikiqtaruk-Herschel Island), and Brendan Rogers, and time series of Landsat reflectance data (Olofsson et al., 2014). In reality, these data were stratified across a design intended to assess the accuracy of change detection. These reference data were adapted to year by year analyses and re-estimated year by year variances (Stehman, 2014). Table 6 provides an accuracy summary across all 31 years.

Table 6. Accuracy summary means

Producer’s accuracy refers to the inverse of error rate of omission, the user’s accuracy refers to the inverse of error rate of commission for each class. The overall accuracy for each year is also provided.

accuracy evergreen_forest deciduous_forest shrub herb sparse_veg barren fen bog shallows water
Producer's Acc 0.8891 0.7506 0.6949 0.6485 0.8096 0.8558 0.6418 0.8854 0.6567 0.8908
User's Acc 0.8266 0.881 0.6893 0.7882 0.7026 0.7902 0.6128 0.8586 0.7575 0.8883
Overall Acc 0.7735 0.7735 0.7735 0.7735 0.7735 0.7735 0.7735 0.7735 0.7735 0.7735

The land cover data does not correct for short-term changes in climate that may impact the spectral appearance of certain areas. For example, droughts may result in wetlands appearing drier, and subsequently being briefly classified as shrubs. Caution is advised when utilizing the wetland data for this reason (Wang et al., 2019).

Data Acquisition, Materials, and Methods

Study area

The study region included the core domain of NASA’s ABoVE project which includes the majority of Alaska and northwestern Canada and covers 4.67 x 106 km2. To facilitate analyses, the Landsat 5 Thematic Mapper and Landsat 7 Enhanced Thematic Mapper data were reprojected into the ABoVE standard Alber’s equal area projection (30-m grid) and were gridded into 180 x 180 km ABoVE standard B grid tiles (Loboda et al., 2019). This resulted in a total of 175 B grid tiles over the ABoVE core domain.

The area mapped was defined by an ocean-land mask that was adopted from the 10-m ocean mask from Natural Earth data repository and manually edited to closely fit the Arctic coastline in our region (Wang et al., 2019).

Data processing

The primary source of data for land cover and disturbance mapping was time series of surface reflectance and brightness temperature from Landsat 5 and Landsat 7 retrieved from the ABoVE Science Cloud, covering the years 1984 to 2014. The Fmask algorithm was used to pre-process the Landsat image time series at each site to remove ephemeral contamination from clouds, snow, ice, cloud shadow, and other atmospheric effects that were not screened by the USGS (Zhu and Woodcock, 2012).

Disturbance detection and data modeling using CCDC

The minimum mapping unit and unit of analysis were 30-meter Landsat pixels within the ABoVE tiling system. For each pixel, the Change Detection and Classification (CCDC) algorithm (Zhu and Woodcock, 2014) was used to detect statistical breaks in time series of surface reflectance data. These breaks are hypothesized to represent land cover changes and the time segments intervening the breaks are assumed to be of stable land cover. The CCDC synthetic data was used to calculate a suite of surface reflectance metrics across bands and seasons to characterize the spectral-seasonal behavior of each time segment (e.g. spring versus summer values or annual amplitude for each band). In addition, several spectral indices were calculated as input features and the seasonal metrics were applied to them as well.

Figure 2. Example of time series and disturbance detection for one Landsat band and one pixel using the Continuous Change Detection and Classification (CCDC) algorithm. In this example, surface reflectance for Landsat band 5 (shortwave-infrared) is shown. Top panel: points are time series of observations reflectance, grey lines are fitted harmonic models, dashed vertical lines are the approximate timing of breaks in land cover, and colored horizontal lines and arrows indicate between-break time segments of stable land cover. CCDC detects and removes outliers in its model fits. Bottom panel: characteristic seasonal cycle of the color-coded time segments indicated in the top panel. Colored points indicate observations falling within the correspondingly colored time segments and lines indicate the harmonic model fit across those values, averaged across all years within the time segment.

Training and validating the land cover mapping algorithms

Training sites were generated across all ABoVE tiles via stratified random sampling (up to 10 per class in each tile) from existing land cover maps for Alaska and Canada circa 2000 (Jin et al., 2017; Wulder et al., 2008), ensuring coverage across a range of climatic conditions, spectral feature space, and ecosystem types. In total, the training data set included 15,628 sampled from the Canadian and Alaskan land cover datasets, 1,531 sites from the 10-m Scotty Creek land cover map (Chasmer et al., 2014), 970 sites from the Toolik Lake tundra vegetation map (Greaves et al., 2017, unpublished dissertation-data provided by personal communication), and 1,553 sites from field photographs for a total of 19,682 training sites (Wang et al., 2019).

map of study area and points

Figure 3. Points represent locations that were included in the clustering algorithm that was used to define the land cover legend and explore the spectral-temporal space within the ABR. Darker blue color indicates the ABoVE Core Study Domain (used for this study), while the lighter blue indicates the ABoVE Extended Study Domain (outside the scope of this study). Red points indicate sites were manually labelled by visual interpretation of very high-resolution imagery or field photographs. Black points were not manually labelled but were included in the clustering and the training data for the Random Forest algorithm (Wang et al., 2019).

Landcover class development

The training site data were clustered using a tree-based clustering method combining the R package treeClust (Buttrey and Whittaker, 2015). Each time segment was defined by a specific statistical cluster and mapped across the domain using Random Forest (Breiman, 2001).  Each time segment was defined by a specific statistical cluster and mapped. The clusters were assigned meaningful land cover labels by visually interpreting a large sample (n = 10,136) of very high-resolution imagery (copyright DigitalGlobe), derived from IKONOS, Geo-Eye, QuickBird, Worldview-2, and Worldview-3, field photographs, and existing higher-resolution land cover maps (Chasmer et al., 2014; Greaves et al., 2017, unpublished dissertation-data provided by personal communication), as well as the Landsat surface reflectance time series. The clusters were cross-walked to land cover based on the interpretation of a set of land cover characteristics observed in these sources.

Landcover in 2014 for same pixel as Figure 1

Figure 4. Land cover (10-class product) for the year 2014 in the ABoVE domain for analysis tile BH03V04. The pixel resolution is 30 m and a tile is 180 x 180 km. Refer to Figure 1 for the landcover in the same tile for the year 1984. (ABoVE_LandCover_BH03V04.tif)

Several post-processing steps resolved certain artifacts, including the removal of spurious land cover changes and the reduction of certain misclassifications. Once clusters have been reassigned to meaningful land covers, per-pixel series of time segments were analyzed and consecutive time segments with the same land cover label (but possibly differing clusters) were merged into single time segments of the same land cover. Additionally, time segments were extended forward or backward in time to cover time periods that had no CCDC output within a pixel. Therefore, missing data should only exist where Fmask identified 100% snow and ice or outside the ocean-land mask. Spurious classifications of water classes in mountainous regions were identified by comparing them with the slope of the pixel, derived from ASTER elevation data, and existing annual Landsat-based maps of global surface water status (Pekel et al., 2016). Those pixels that were mapped as water but had non-zero slope and were not classified as water by Pekel et al. 2016, were set to their second most likely class after water, using the Random Forest probabilities (Wang et al., 2019).

 

Data Access

These data are available through the Oak Ridge National Laboratory (ORNL) Distributed Active Archive Center (DAAC).

ABoVE: Landsat-derived Annual Dominant Land Cover Across ABoVE Core Domain, 1984-2014

Contact for Data Center Access Information:

References

Breiman, L. Machine Learning (2001) 45:5. https://doi.org/10.1023/A:1010933404324

Buttrey, S.E., Whitaker, L.R., 2015. treeClust: An R package for tree-based clustering dissimilarities. The R Journal (2015) 7:2, pages 227-236

Chasmer, L., Hopkinson, C., Veness, T., Quinton, W., Baltzer, J., 2014. A decision-tree classification for low-lying complex land cover types within the zone of discontinuous permafrost. Remote Sensing of Environment 143, 73–84. https://doi.org/10.1016/j.rse.2013.12.016

Greaves, H. E. (2017). Applying Lidar and High-Resolution Multispectral Imagery for Improved Quantification and Mapping of Tundra Vegetation Structure and Distribution in the Alaskan Arctic (Doctoral dissertation, University of Idaho). https://ui.adsabs.harvard.edu/abs/2017PhDT........51G/abstract

Jin, S., Yang, L., Zhu, Z., Homer, C., 2017. A land cover change detection and classification protocol for updating Alaska NLCD 2001 to 2011. Remote Sensing of Environment 195, 44–55. https://doi.org/10.1016/j.rse.2017.04.021

Olofsson, P., Foody, G.M., Herold, M., Stehman, S.V., Woodcock, C.E., Wulder, M.A., 2014. Good practices for estimating area and assessing accuracy of land change. Remote Sensing of Environment 148, 42–57. https://doi.org/10.1016/j.rse.2014.02.015

Pekel, J.-F., Cottam, A., Gorelick, N., Belward, A.S., 2016. High-resolution mapping of global surface water and its long-term changes. Nature 540, 418. https://doi.org/10.1038/nature20584

Stehman, S.V., 2014. Estimating area and map accuracy for stratified random sampling when the strata are different from the map classes. International Journal of Remote Sensing 35, 4923–4939. https://doi.org/10.1080/01431161.2014.930207

Wang, J. A., Sulla-Menashe, D., Woodcock, C. E., Sonnentag, O., Keeling, R. F., & Friedl, M. A. 2020. Extensive Land Cover Change Across Arctic-Boreal Northwestern North America from Disturbance and Climate Forcing. Global Change Biology 26, 807–822. https://doi.org/10.1111/gcb.14804

Wulder, M.A., White, J.C., Cranny, M., Hall, R.J., Luther, J.E., Beaudoin, A., Goodenough, D.G., Dechka, J.A., 2008. Monitoring Canada’s forests. Part 1: Completion of the EOSD land cover project. Canadian Journal of Remote Sensing 34, 549–562. https://doi.org/10.5589/m08-066

Zhu, Z., Woodcock, C.E., 2014. Continuous change detection and classification of land cover using all available Landsat data. Remote sensing of Environment 144, 152–171. https://doi.org/10.1016/j.rse.2014.01.011 

Zhu, Z., Woodcock, C.E., 2012. Object-based cloud and cloud shadow detection in Landsat imagery. Remote Sensing of Environment 118, 83–94. https://doi.org/10.1016/j.rse.2011.10.028