Map forests for the period 2000 – 2017
step 1: add your aoi
// set the area of interest var chl = ee.FeatureCollection("users/servirmekong/Vietnam/CHL_Boundary").geometry();
step 2: import the tree canopy cover and tree height
// import tree height var th = ee.ImageCollection("projects/servir-mekong/yearly_primitives_smoothed/tree_height"); // import tree canopy cover var cc = ee.ImageCollection("projects/servir-mekong/yearly_primitives_smoothed/tree_canopy");
step 3: Set forest definitions
// Forest definitions // Tree canopy greater than (%) var ccThreshold = 10; // Tree height greater than (meters) var thThreshold = 5; // minimum mapping unit var mmu = 5;
step 4: Set time period
// set year of interest var year = 2017; var startDate = ee.Date.fromYMD(year,1,1); var endDate = ee.Date.fromYMD(year,12,31);
step 5: filter for time
// get image for the year var thImage = ee.Image(th.filterDate(startDate,endDate).first()); var ccImage = ee.Image(cc.filterDate(startDate,endDate).first());
step 6: Create boolean maps for treeheight and canopy cover using the thresholds
// Set the thresholds of the image var treeheight = thImage.gt(thThreshold); var canopy = ccImage.gt(ccThreshold);
step 7: Create boolean maps for treeheight and canopy cover using the thresholds
// create the forest layer var forest = treeheight.add(canopy).eq(2) forest = forest.mask(forest) // add the forest layer to the map Map.addLayer(forest.clip(chl),{palette:"darkgreen"},"Forests without minimum mapping unit");
step 8: Add the minimum mapping unit for the final map
// apply the minimum mapping unit var map = forest.connectedPixelCount(mmu+2).gte(mmu); Map.addLayer(map.updateMask(map).clip(chl),{palette:"darkgreen"},"Forests with minimum mapping unit");
Find the example here