Mapping Polylepis Forest Using Sentinel, PlanetScope Images, and Topographical Features with Machine Learning

Globally, there is a significant trend in the loss of native forests, including those of the Polylepis genus, which are essential for soil conservation across the Andes Mountain range. These forests play a critical role in regulating water flow, promoting soil regeneration, and retaining essential nutrients and sediments, thereby contributing to the soil conservation of the region. In Ecuador, these forests are often fragmented and isolated in areas of high cloud cover, making it difficult to use remote sensing and spectral vegetation indices to detect this forest species. This study developed twelve scenarios using medium- and high-resolution satellite data, integrating datasets such as Sentinel-2 and PlanetScope (optical), Sentinel-1 (radar), and the Sigtierras project topographic data. The scenarios were categorized into two groups: SC1-SC6, combining 5 m resolution data, and SC7-SC12, combining 10 m resolution data. Additionally, each scenario was tested with two target types: multiclass (distinguishing Polylepis stands, native forest, Pine, Shrub vegetation, and other classes) and binary (distinguishing Polylepis from non-Polylepis). The Recursive Feature Elimination technique was employed to identify the most effective variables for each scenario. This process reduced the number of variables by selecting those with high importance according to a Random Forest model, using accuracy and Kappa values as criteria. Finally, the scenario that presented the highest reliability was SC10 (Sentinel-2 and Topography) with a pixel size of 10 m in a multiclass target, achieving an accuracy of 0.91 and a Kappa coefficient of 0.80. For the Polylepis class, the User Accuracy and Producer Accuracy were 0.90 and 0.89, respectively. The findings confirm that, despite the limited area of the Polylepis stands, integrating topographic and spectral variables at a 10 m pixel resolution improves detection accuracy.