[Script] Créez une Carte de Perte en Sol (RUSLE) sur Google Earth Engine

// 1️⃣ Définir la limite du Maroc
var limites_maroc = ee.Geometry.Polygon([
  [[-17.0, 36.0], [-1.0, 36.0], [-1.0, 21.0], [-17.0, 21.0], [-17.0, 36.0]]
]);

// Charger le MNT SRTM et le clipper
var dem = ee.Image("USGS/SRTMGL1_003").clip(limites_maroc);

// Calcul de la pente et du facteur LS
var slope = ee.Terrain.slope(dem);
var LS = slope.expression(
  "pow((sin(slope * 3.1416 / 180) / 0.0896), 1.3)", 
  {slope: slope}
).clip(limites_maroc);

// Charger les précipitations CHIRPS et les clipper
var precip = ee.ImageCollection("UCSB-CHG/CHIRPS/DAILY")
  .filterDate('2020-01-01', '2020-12-31')
  .sum()
  .clip(limites_maroc);

// Facteur R basé sur la précipitation annuelle
var R = precip.multiply(0.5).clip(limites_maroc);

// Charger le facteur K (érodibilité du sol)
var K = ee.Image("OpenLandMap/SOL/SOL_TEXTURE-CLASS_USDA-TT_M/v02")
  .select('b0')
  .multiply(0.01)
  .clip(limites_maroc);

// Charger la couverture du sol (MODIS) et clipper
var landcover = ee.ImageCollection("MODIS/006/MCD12Q1")
  .filterDate('2020-01-01', '2020-12-31')
  .first()
  .select('LC_Type1')
  .clip(limites_maroc);

// Facteur C basé sur l'occupation des sols
var C = landcover.expression(
  "(b1 == 1) ? 0.01 : (b1 == 2) ? 0.05 : (b1 == 3) ? 0.2 : (b1 == 4) ? 0.4 : 0.5",
  {b1: landcover}
).clip(limites_maroc);

// Facteur P (fixé à 1)
var P = ee.Image.constant(1).clip(limites_maroc);

// Calcul final de l'érosion (RUSLE)
var erosion = R.multiply(K).multiply(LS).multiply(C).multiply(P).clip(limites_maroc);

// Affichage des données
Map.centerObject(limites_maroc, 6);
Map.addLayer(erosion, {min: 0, max: 200, palette: ['blue', 'green', 'yellow', 'orange', 'red', 'brown']}, "Érosion (RUSLE) - Maroc");