Dominant tree species map of Switzerland
We created a tree species map of Switzerland for the dominant tree species in the forested areas. The spatial resolution of the map is 10 m and the coordinate system is ETRS89-extended / LAEA Europe (EPSG 3035). The map comprises Sentinel-2 index time series from the year 2020, a digital elevation model and species reference data from the Swiss National Forest Inventory. The map is available as raster (.tif) or vector dataset (.gpkg).
Access will be granted upon request.
In total, the following 15 species were mapped:
Abies alba, Acer pseudoplatanus, Alnus glutinosa, Alnus incana, Betula pendula, Castanea sativa, Fagus sylvatica, Fraxinus excelsior, Picea abies, Pinus cembra, Pinus mugo arborea, Pinus sylvestris, Quercus petraea, Quercus robur, Sorbus aucuparia.
Approach
Data
- Swiss National Forest Inventory Data (stand species with > 60 % dominance in upper canopy; on at least more than 9 plots dominant)
- Sentinel-2 time series (2020, Indices: CCI, CIRE, NDMI, EVI, NDVI)
- Digital elevation model (DEM) (swissalti3d, 5 m)
- Biogeographical regions (Federal Office for the Environment FOEN)
- Forest mask 2017 (Approach: Waser et al., 2015)
Modeling approach
We identified the most meaningful variables that led to separation of the respective groups by using random forest models with a forward feature selection (Meyer et al., 2018; Ververidis & Kotropoulos, 2005). In this approach, the final random forest model is solely built from the selected meaningful variables. By identifying meaningful variables, we can determine which variables might influence the grouping. Further, to avoid overfitting and overly optimistic results, we applied 10-fold spatial cross-validation and put all pixels from a plot in the same spatial fold.
The modeling was realized using the CAST package in R (Meyer et al., 2022), based on the well-known caret package (Kuhn, 2022). We used the ranger package in R (Wright & Ziegler, 2017) to implement the random forest models, due to its short computation time.
Training data for modeling
- 295 Sentinel-2, DEM & Biogeographical variables
- 10525 tree species pixels
Selected variables for final model
- EVI of 2020.05.16
- NDMI of 2020.03.12
- CIRE of 2020.04.16
- NDMI of 2020.07.05
- CCI of 2020.05.11
- dem
- CCI of 2020.08.14
- NDMI of 2020.08.24
- CCI of 2020.12.22
- NDMI of 2020.04.21
- NDMI of 2020.11.17
- NDMI of 2020.08.09
- CIRE of 2020.03.22
- CIRE of 2020.08.09
- CCI of 2020.11.02
- CIRE of 2020.06.10
Overall Accuracy of final model
Nationwide prediction
- Predicted throughout forest mask 2017 (Approach: Waser et al., 2015)
- Not applied on incomplete Sentinel-2 time series (own category in final map: incomplete_ts)
- Applied the Area of Applicability (Meyer 2022) to sort out pixels outside of the feature space; basically where the model had not the same values for pixels as in the available training data
Be aware that the map is only validated with the training data itself, an independent validation with other data sources remains missing
References
-
Kuhn, M. (2022). Classification and Regression Training. 6.0-93.
-
Meyer, H., Reudenbach, C., Hengl, T., Katurji, M., & Nauss, T. (2018). Improving performance of spatio-temporal machine learning models using forward feature selection and target-oriented validation. Environmental Modelling and Software, 101, 1-9. https://doi.org/10.1016/j.envsoft.2017.12.001
-
Meyer, H., Milà, C., & Ludwig, M. (2022). CAST: 'caret' Applications for Spatial-Temporal Models. 0.7.0.
-
Ververidis, D., & Kotropoulos, C. (2005). Sequential forward feature selection with low computational cost. 2005 13th European Signal Processing Conference.
-
Waser, L., Fischer, C.,Wang, Z., & Ginzler, C. (2015). Wall-to-Wall Forest Mapping Based on Digital Surface Models from Image-Based Point Clouds and a NFI Forest Definition. Forests, 6, 12, 4510–4528.
-
Wright, M. N., & Ziegler, A. (2017). ranger: A Fast Implementation of Random Forests for High Dimensional Data in C++ and R. Journal of Statistical Software, 77(1), 1-17. https://doi.org/doi:10.18637/jss.v077.i01
This work was supported by:
- Swiss National Science Foundation
(Grant/Award: 200021_184605)