Forests in Southeast Asia are experiencing some of the highest rates of deforestation and degradation in the world, with natural forest species being replaced by cropland and plantation monoculture. In this work, we have developed an innovative method to accurately map rubber and palm oil plantations using fusion of Landsat-8, Sentinel 1 and 2. We applied cloud and shadow masking, bidirectional reflectance distribution function (BRDF), atmospheric and topographic corrections to the optical imagery and a speckle filter and harmonics for Synthetic Aperture Radar (SAR) data. In this workflow, we created yearly composites for all sensors and combined the data into a single composite. A series of covariates were calculated from optical bands and sampled using reference data of the land cover classes including surface water, forest, urban and built-up, cropland, rubber, palm oil and mangrove. This training dataset was used to create biophysical probability layers (primitives) for each class. These primitives were then used to create land cover and probability maps in a decision tree logic and Monte-Carlo simulations. Validation showed good overall accuracy (84%) for the years 2017 and 2018. Filtering for validation points with high error estimates improved the accuracy up to 91%. We demonstrated and concluded that error quantification is an essential step in land cover classification and land cover change detection. Our overall analysis supports and presents a path for improving present assessments for sustainable supply chain analyses and associated recommendations.