Abstract
Landscape metrics (LM) play a crucial role in fields such as urban planning, ecology, and environmental research, providing insights into the ecological and functional dynamics of ecosystems. However, in dynamic systems, generating thematic maps for LM analysis poses challenges due to the substantial data volume required and issues such as cloud cover interruptions. The aim of this study was to compare the accuracy of land cover maps produced by three temporal aggregation methods: median reflectance, maximum normalised difference vegetation index (NDVI) and a two-date image stack using Sentinel-2 (S2), and then to analyse their implications for LM calculation. The Google Earth Engine platform facilitated data filtering, image selection, and aggregation while mitigating cloud cover effects. A random forest algorithm was employed to classify five land cover classes across ten sites, with classification accuracy assessed using global measurements and Kappa index. LM were then quantified. The analysis revealed that S2 data provided a high-quality, cloud-free dataset suitable for analysis, ensuring a minimum of 25 cloud-free pixels over the study period. The two-date and median methods exhibited superior land cover classification accuracy compared to the max NDVI method. In particular, the two-date method resulted in lower fragmentation-heterogeneity and complexity metrics in resulting maps compared to the median and max NDVI methods. Nevertheless, median method holds promise for integration into operational land cover mapping programs, particularly for larger study areas exceeding the width of S2 swath coverage. These results highlight the importance of appropriate temporal aggregation techniques when using satellite data for landscape analysis and monitoring.