// Load the Hansen Global Forest Change dataset var gfc = ee.Image('UMD/hansen/global_forest_change_2023_v1_11'); // Select the layers for tree cover in 2000 and forest loss var treeCover2013 = gfc.select('treecover2000'); var forestLoss = gfc.select('loss'); // Clip the images to the AOI var treeCover2000_aoi = treeCover2013.clip(edo); var forestLoss_aoi = forestLoss.clip(edo); // Set visualization parameters var treeCoverVis = { min: 0, max: 100, palette: ['red', 'yellow', 'green'] }; var forestLossVis = { min: 0, max: 1, palette: ['red'] }; // Add the layers to the map Map.centerObject(edo, 6); Map.addLayer(treeCover2000_aoi, treeCoverVis, 'Tree Cover 2000'); Map.addLayer(forestLoss_aoi, forestLossVis, 'Forest Loss'); // Calculate the total area of forest loss var lossPixels = forestLoss_aoi.reduceRegion({ reducer: ee.Reducer.sum(), geometry: edo, scale: 30, maxPixels: 1e9 }); // Convert the result to hectares (1 pixel = 30m x 30m = 900 m² = 0.09 hectares) var lossHectares = ee.Number(lossPixels.get('loss')).multiply(0.09); print('Total forest loss in hectares:', lossHectares); // Optional: Export the forest loss layer to Google Drive Export.image.toDrive({ image: forestLoss_aoi, description: 'Forest_Loss', scale: 30, region: edo, fileFormat: 'GeoTIFF' });
What will be the unit of K factor if done by this method ?? 2. I am getting error in last step, I am not able to convert polygon into raster , can you tell me what should i do , i am getting an error like soil_clip.K_factor doesnt exist in table but i have exact same in table. Error 000622 and Error 000728. If you can tell then it willl be of great help
I did this for 2 locations but i got f_silt as 0.3, I am guessing its some problem and in one location i got f_oc as 1 for 4 type of soil in that location which again sounds like an error. because of this i am getting K factor value as in range of 0.005 which is very low wrt given values by other authors. Can you tell me where I am going wrong with this??
This was so helpful. I'm currently investigating land use and land cover changes in a catchment and I had to use IDRISI TERRSET. I know understand the subset method,,training stage and classification stage a lot better.
Thanks. I have added an ascii land cover map from arcgis to terrset. But terrset opens/displays the ascii file using a color ramp not as per original land classes in arcgis. Kindly could you know how treat the ascii file?
Qgis has different plugins for different analysis but idrisi terrset is the best for prediction. you can either use land change modeller or cellula automata.
