“…The downloaded Landsat 7 (ETM+) and Landsat 8 (OLI) MS bands, R, G, B, NIR, SWIR1, and SWIR2, including Sentinel-2 (MSI) Level-1C MS bands, B, G, R, VRE, VRE2, VRE3, NIR, NIR2, SWIR2, and SWIR3, were automatically processed using the SCP plugin by converting them from DN [Landsat] and scaled top of atmosphere (TOA) reflectance [Sentinel] into the TOA reflectance to reduce the inter-scene variability through a normalization for solar irradiance. Atmospheric correction of all images was carried out using an image-based technique called Dark Object Subtraction (DOS1) [77], as cited in [76]. Ordinary least squares regression (OLS) equations from Roy et al [78] were used to normalize the reflectance of one Landsat sensor to the other (ETM+ to OLI).…”