China, as the world’s biggest soybean importer and fourth-largest producer, needs accurate mapping of its planting areas for global food supply stability. The challenge lies in gathering and collating ground survey data for different crops. We proposed a spatiotemporal migration method leveraging vegetation indices’ temporal characteristics. This method uses a feature space of six integrals from the crops’ phenological curves and a concavity-convexity index to distinguish soybean and non-soybean samples in cropland. Using a limited number of actual samples and our method, we extracted features from optical time-series images throughout the soybean growing season. The cloud and rain-affected data were supplemented with SAR data. We then used the random forest algorithm for classification. Consequently, we developed the 10-meter resolution ChinaSoybean10 maps for the ten primary soybean-producing provinces from 2019 to 2022. The map showed an overall accuracy of about 93%, aligning significantly with the statistical yearbook data, confirming its reliability. This research aids soybean growth monitoring, yield estimation, strategy development, resource management, and food scarcity mitigation, and promotes sustainable agriculture.