Paddy rice cropping patterns (PRCPs) play important roles in both agroecosystem modeling and food security. Although paddy rice maps have been generated over several regions using satellite observations, few studies have focused on mapping diverse smallholder PRCPs, which include crop rotation and are dominant cropping structures in South China. Here, an approach called the feature selection and hierarchical classification (FSHC) method was proposed to effectively identify paddy rice and its rotation types. Considering the cloudy and rainy weather in South China, a harmonized Landsat and Sentinel-2 (HLS) surface reflectance product was employed to increase high-quality observations. The FSHC method consists of three processes: cropping intensity mapping, feature selection, and decision tree (DT) model development. The FSHC performance was carefully evaluated using crop field samples obtained in 2018 and 2019. Results suggested that the derived cropping intensity map based on the Savitzky–Golay (S-G) filtered normalized difference vegetation index (NDVI) time series was reliable, with an overall accuracy greater than 93%. Additionally, the optimal spectral (i.e., normalized difference water index (NDWI) and land surface water index (LSWI)) and temporal (start-of-season (SOS) date) features for distinguishing different PRCPs were successfully identified, and these features are highly related to the critical growth stage of paddy rice. The developed DT model with three hierarchical levels based on optimal features performed satisfactorily, and the identification accuracy of each PRCP can be achieved approximately 85%. Furthermore, the FSHC method exhibited similar performances when mapping PRCPs in adjacent years. These results demonstrate that the proposed FSHC approach with HLS data can accurately extract diverse PRCPs over fragmented croplands; thus, this approach represents a promising opportunity for generating refined crop type maps.

https://www.mdpi.com/2072-4292/15/4/1034