Accurate crop mapping is vital for monitoring agricultural resources, food security, and ecosystem sustainability. Advances in high-resolution sensing technologies now enable precise, large-scale crop mapping, improving agricultural management and decision-making. However, in scenarios where balancing precision and computational resources is important, obtaining the optimal feature combination (especially newly proposed features) and strategies from the rich feature sets contained in multi-source remote sensing imagery remains one of the challenges. In this paper, we propose a hierarchical feature optimization method, incorporating a newly reported vegetation feature, for mapping crop types by combining the Sentinel-1 Synthetic Aperture Radar (SAR) and Sentinel-2 optical imagery within the Google Earth Engine (GEE) platform. The method first calculates spectral features, texture features, polarization features, vegetation index features, and crop phenological features, with a particular focus on infrared band features and the newly developed Kernel Normalized Difference Vegetation Index (kNDVI). These 126 features are then selected to construct 15 crop type mapping models based on different feature combinations and a random forest (RF) classifier. Feature selection was performed using the feature correlation analysis and random forest recursive feature elimination (RF-RFE) to identify the optimal subset. The experiment was conducted in the Linhe region, covering an area of 2333 km2. The resulting 10 m crop map, generated by the optimal model (Model 15) with 34 key features, demonstrated that integrating multi-source features significantly enhances mapping accuracy. The model achieved an overall accuracy of 90.10% across five crop types (corn, wheat, sunflower, soybean, and beet), outperforming other representative feature optimization methods, Relief-F (87.50%) and CFS (89.60%). The study underscores the importance of feature optimization and reduction of redundant features while also showcasing the effectiveness of red edge and infrared features, as well as the kNDVI, in mapping crop type.
Liu et al. (Fri,) studied this question.