Wildfire Risk Mapping with MODIS and DeepSeek

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

// Filter MODIS active fire data for 2023
var fires = ee.ImageCollection('MODIS/061/MCD14ML/C6_1')
  .filterDate('2023-01-01', '2023-12-31')
  .filterBounds(roi);

// Function to mask out pixels with low confidence
var maskLowConfidence = function(image) {
  return image.updateMask(image.select('confidence').gt(50));
};

// Apply the mask
vAR firesFiltered = fires.map(maskLowConfidence);

//Mosaic the fire collection and select the FRP band
var fire_frp = firesFiltered.mosaic().select('frp');

// Function to calculate NDVI
var calculateNDVI = function(image) {
  return image.normalizedDifference(['Nir', 'Red']).rename('NDVI');
};

// Function to calculate EVI
var calculateEVI = function(image) {
  return image.expression(
    '2.5 * ((NIR - RED) / (NIR + 6 * RED - 7.5 * BLUE + 1))', {
      'NIR': image.select('Nir'),
      'RED': image.select('Red'),
      'BLUE': image.select('Blue')
    }).rename('EVI');
};

// Load MODIS surface reflectance data
var modisSR = ee.ImageCollection('MODIS/061/MOD09GA/006')
  .filterDate('2023-01-01', '2023-12-31')
  .filterBounds(roi);

// Apply scaling factors and rename bands
var modisSRScaled = modisSR.map(function(image) {
  return image.multiply(0.0001).select('sur_refl_b01', 'sur_refl_b02', 'sur_refl_b04').rename('Red', 'Nir', 'Blue').copyProperties(image, ['system:time_start']);
});

// Calculate NDVI and EVI
var ndvi = modisSRScaled.map(calculateNDVI).mosaic();
var evi = modisSRScaled.map(calculateEVI).mosaic();

//Export the risk map
Export.image.toDrive({
  image: riskMap,
  description: 'Wildfire_Risk_Map_Morocco_2023',
  folder: 'GEE_Exports',
  region: roi,
  scale: 30, // Adjust scale as needed
  maxPixels: 1e13
});