Solar Potential Mapping

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

// Import ERA5-Land hourly data for 2022
var era5 = ee.ImageCollection('ECMWF/ERA5_LAND/HOURLY')
  .filterDate('2022-01-01', '2022-12-31')
  .filterBounds(roi);

// Print the number of images in the collection
print('Number of ERA5 images:', era5.size());

// Import SRTM DEM data
var dem = ee.Image('USGS/SRTMGL1_003').clip(roi);

// Import land cover data
var landcover = ee.Image('COPERNICUS/Landcover/100m/Proba-V-C3/Global').select('discrete_classification').clip(roi);

// Calculate daily solar irradiance (example, needs adjustment based on ERA5 bands)
var dailyIrradiance = era5.map(function(image) {
  // Assuming 'surface_solar_radiation_downwards' is the relevant band and converting from J/m² to kWh/m²
  var irradiance = image.select('surface_solar_radiation_downwards').multiply(2.777e-7); // J/m² to kWh/m²
  return irradiance.set('system:time_start', image.get('system:time_start'));
});

// Calculate annual solar irradiance
var annualIrradiance = dailyIrradiance.sum();

// Calculate slope and aspect from DEM
var slope = ee.Terrain.slope(dem);
var aspect = ee.Terrain.aspect(dem);

// Forest cover mask (assuming forest is classified as 5 or 6 in landcover)
var forestMask = landcover.eq(5).or(landcover.eq(6));

// Export the suitability map to Google Drive
Export.image.toDrive({
  image: annualIrradiance,
  description: 'Solar_Irradiance_North_Morocco',
  folder: 'GEE_Exports',
  scale: 30, // Adjust scale as needed
  region: roi,
  maxPixels: 1e13
});