A Phenology-Adaptive Rubber Plantation Mapping (PARM) Framework Coupling Sentinel-1 SAR and Optimally Selected Spectral Indices Across Heterogeneous Tropical Regions

Accurate mapping of rubber plantations is essential for sustainable land management and forest conservation in tropical regions. However, existing methods face two major challenges: persistent cloud cover limits the effectiveness of optical remote sensing in tropical areas, and regional phenological heterogeneity hinders the transferability of fixed-parameter approaches. This study proposes a Phenology-Adaptive Rubber Plantation Mapping (PARM) framework that integrates Sentinel-1 SAR time-series data with optimally selected spectral indices through a cascading constraint architecture. The framework operates as a structurally coherent system wherein SAR-derived phenological anchors explicitly govern downstream optical analysis across three internally dependent modules. First, three key phenological nodes—leaf-off start (LOS), fastest greening point (FGP), and full canopy point (FCP)—are extracted directly from SAR VH-polarization backscatter time series, enabling cloud-independent extraction of phenological temporal anchors. Second, the Jeffries–Matusita (JM) distance, evaluated within SAR-constrained phenological windows, is employed to identify the optimal vegetation and water indices for each region from six candidate spectral indices. Third, a time-weighted Rubber Plantation Discrimination Index (RPDI) is constructed using the selected indices and locally extracted phenological nodes, thereby amplifying the coupled signals of canopy greenness and moisture dynamics during critical phenological transitions. The framework was validated in Hainan Island and Vietnam, two regions with contrasting phenological regimes, using a spatial-block partitioning protocol (leave-one-subregion-out combined with DBSCAN-based clustering) designed to prevent samples from the same plantation from occurring in both training and test subsets. Within the Dynamic World forest mask, PARM achieved overall accuracies of 92.04% and 91.24%, respectively (93.57% and 91.00% on the fully held-out Qionghai City and Gia Lai province subregions), with Kappa coefficients exceeding 0.81 in both regions, consistently outperforming schemes based on raw spectral bands, individual spectral indices or direct multi-source time-series stacking. Error structure analysis revealed that residual classification failures are primarily associated with landscape fragmentation, stand immaturity, and residual cloud contamination, delineating the generalizability boundaries of the framework. These results demonstrate that tightly coupling SAR-based phenological characterization with adaptive optical index selection through a cascading constraint architecture provides a reliable foundation for rubber plantation mapping in cloud-prone tropical regions.