“…Following the cloud mask, we calculated spectral indices, namely the Normalized Difference Vegetation Index (NDVI) [44,45] and the Normalized Difference Water Index (NDWI) [46], as well as Tasseled cap wetness, greenness and brightness (Table 1) [47]. We implemented a pixel-based compositing approach before the classification procedure and calculated temporal metrics [48][49][50][51] (e.g., 10th, 25th, 50th, 75th, 80th and 90th percentiles of annual time series, as well as a minima and maxima for both annual and growing season time series), which resulted in composites for each target year. We extracted eleven types of image features from each compositing result, including six surface reflectance (green, blue, red, NIR, SWIR1 and SWIR2) bands and the spectral indices (NDVI, NDWI, TC greenness, wetness, brightness).…”