Evaluation and Selection of Multi-Spectral Indices to Classify Vegetation Using Multivariate Functional Principal Component Analysis

The identification, classification and mapping of different plant communities and habitats is of fundamental importance for defining biodiversity monitoring and conservation strategies. Today, the availability of high temporal, spatial and spectral data from remote sensing platforms provides dense time series over different spectral bands. In the case of supervised mapping, time series based on classical vegetation indices (e.g., NDVI, GNDVI, …) are usually input characteristics, but the selection of the best index or set of indices (which guarantees the best performance) is still based on human experience and is also influenced by the study area. In this work, several different time series, based on Sentinel-2 images, were created exploring new combinations of bands that extend the classic basic formulas as the normalized difference index. Multivariate Functional Principal Component Analysis (MFPCA) was used to contemporarily decompose the multiple time series. The principal multivariate seasonal spectral variations identified (MFPCA scores) were classified by using a Random Forest (RF) model. The MFPCA and RF classifications were nested into a forward selection strategy to identify the proper and minimum set of indices’ (dense) time series that produced the most accurate supervised classification of plant communities and habitat. The results we obtained can be summarized as follows: (i) the selection of the best set of time series is specific to the study area and the habitats involved; (ii) well-known and widely used indices such as the NDVI are not selected as the indices with the best performance; instead, time series based on original indices (in terms of formula or combination of bands) or underused indices (such as those derivable with the visible bands) are selected; (iii) MFPCA efficiently reduces the dimensionality of the data (multiple dense time series) providing ecologically interpretable results representing an important tool for habitat modelling outperforming conventional approaches that consider only discrete time series.