Identification of soybean based on Sentinel-1/2 SAR and MSI imagery under a complex planting structure

Abstract

Soybeans, as the major global oil seed crops, are among the large-scale agricultural products imported into China. In this way, the accurate identification of their growing areas is the basis for agricultural decision-making and planting structure adjustment, which are of great significance to the national food security. Currently, there is insufficient research on the identification of soybean, using remote sensing, under complex agricultural planting structures. In this paper, Guoyang County, a typical soybean-producing region in the north of Anhui Province, China, was selected as the study area, and then the multi-temporal Sentinel-1/2 (S-1/2) microwave and optical multispectral data were integrated to obtain the spatial distribution of the planting areas, via a stepwise hierarchical extraction strategy. Considering the Jeffries-Matusita (JM) distance, August 18 (i.e., the early pod setting stage) was correspondingly determined as more appropriate for soybean extraction, compared with other time phases. For this purpose, a set of rules were established to eliminate the non-cropland pixels, and the total distribution of the field vegetation was derived accordingly. In total, 30 candidate features (viz. 10 spectral bands, 11 vegetation indices, 4 texture features, and 5 microwave polarization features) were selected, and their importance was evaluated based on the ReliefF algorithm. With the support of the typical cover type samples, the weight assessment method for the ReliefF feature was further bundled with three machine-learning (ML) algorithms, namely, the random forest (RF), the back-propagation neural network (BPNN), and the support vector machine (SVM), to single out the optimum subset of the features for soybean identification. The results finally revealed that the ReliefF-RF model outperformed the other two. With nine optimum features, Kappa coefficient also reached up to 0.77–0.85. Furthermore, the extraction accuracy produced by the optimum subset of the features was significantly greater than the original 10 spectral bands (Kappa coefficient of 0.75–0.84), albeit slightly lower than the effect of the total 30 features (Kappa coefficient of 0.78–0.85).