Cropland Mapping with Sentinel-2 and Dynamic World

// Define region of interest (Rif Mountains and Middle Atlas, Morocco)
var roi = ee.Geometry.Rectangle([-6, 32, -2, 35]);

// Function to mask clouds using the QA60 band
var maskS2clouds = function(image) {
  var qa = image.select('QA60');
  // Bits 10 and 11 are clouds and cirrus, respectively.
  var cloudBitMask = 1 << 10;
  var cirrusBitMask = 1 << 11;
  // Both flags should be set to zero, indicating clear conditions.
  var mask = qa.bitwiseAnd(cloudBitMask).eq(0)
      .and(qa.bitwiseAnd(cirrusBitMask).eq(0));
  return image.updateMask(mask).divide(10000);
};

// Load Sentinel-2 SR data and filter by date and region
var s2 = ee.ImageCollection('COPERNICUS/S2_SR')
    .filterDate('2023-01-01', '2023-12-31')
    .filterBounds(roi)
    .map(maskS2clouds);

// Function to calculate NDVI
var addNDVI = function(image) {
  var ndvi = image.normalizedDifference(['B8', 'B4']).rename('NDVI');
  return image.addBands(ndvi);
};

// Function to calculate EVI
var addEVI = function(image) {
  var evi = image.expression(
      '2.5 * (NIR - RED) / (NIR + 6 * RED - 7.5 * BLUE + 1)', {
        'NIR': image.select('B8'),
        'RED': image.select('B4'),
        'BLUE': image.select('B2')
      }).rename('EVI');
  return image.addBands(evi);
};

// Apply NDVI and EVI calculations
var s2_ndvi_evi = s2.map(addNDVI).map(addEVI);

// Load Dynamic World data
var dw = ee.ImageCollection('GOOGLE/DYNAMIC_WORLD/V1')
    .filterDate('2023-01-01', '2023-12-31')
    .filterBounds(roi);

// Get the most frequent land cover class from Dynamic World
var dw_mode = dw.select('label').reduce(ee.Reducer.mode());

// Export the cropland map to Google Drive
Export.image.toDrive({
  image: dw_mode,
  description: 'Cropland_Map_Morocco_2023',
  folder: 'GEE_Exports',
  scale: 10,
  region: roi,
  maxPixels: 1e13
});