//MAPAS COMBINADOS DE PENDIENTE Y ORIENTACION DE LADERAS WWW.GISANDBEERS.COM
//www.gisandbeers.com/analisis-morfologicos-combinados-pendiente-ladera-google-earth-engine/

//Acceso a coleccion de datos altitudinales
var DEM = ee.Image("USGS/3DEP/10m").select("elevation");

//Procesado de modelos digitales de terreno basados en morfologia de orientacion, pendiente y sombras
var Pendiente =  ee.Terrain.slope (DEM);
var Laderas =  ee.Terrain.aspect (DEM);
var Hillshade = ee.Terrain.hillshade(DEM, 300, 40);

// Puntuacion de valores morfologicos
var ClasesPendiente = Pendiente  
          .where(Pendiente.lte(5), 10)
          .where(Pendiente.gt(5).and(Pendiente.lte(20)), 20) 
          .where(Pendiente.gt(20).and(Pendiente.lte(40)), 30)
          .where(Pendiente.gt(40), 40); 
var ClasesOrientacion = Laderas  
          .where(Laderas.lte(0), 0)
          .where(Laderas.gt(0).and(Laderas.lte(22.5)), 1)
          .where(Laderas.gt(22.5).and(Laderas.lte(67.5)), 2) 
          .where(Laderas.gt(67.5).and(Laderas.lte(112.5)), 3)
          .where(Laderas.gt(112.5).and(Laderas.lte(157.5)), 4)
          .where(Laderas.gt(157.5).and(Laderas.lte(202.5)), 5)
          .where(Laderas.gt(202.5).and(Laderas.lte(247.5)), 6)
          .where(Laderas.gt(247.5).and(Laderas.lte(292.5)), 7)
          .where(Laderas.gt(292.5).and(Laderas.lte(337.5)), 8)
          .where(Laderas.gt(337.5), 1); 
var MDT = ClasesPendiente.add(ClasesOrientacion);

//Visualizacion de MDT morfologicos
Map.addLayer (Pendiente,{
  min: 0.0,
  max: 90.0,
  palette: [
    '3ae237', 'b5e22e', 'd6e21f', 'fff705', 'ffd611', 'ffb613', 'ff8b13',
    'ff6e08', 'ff500d', 'ff0000', 'de0101', 'c21301', '0602ff', '235cb1',
    '307ef3', '269db1', '30c8e2', '32d3ef', '3be285', '3ff38f', '86e26f'
  ]},
'Pendientes',false);

Map.addLayer (Laderas,{
  min: 0.0,
  max: 360.0,
  palette: [
    '3ae237', 'b5e22e', 'd6e21f', 'fff705', 'ffd611', 'ffb613', 'ff8b13',
    'ff6e08', 'ff500d', 'ff0000', 'de0101', 'c21301', '0602ff', '235cb1',
    '307ef3', '269db1', '30c8e2', '32d3ef', '3be285', '3ff38f', '86e26f'
  ]},
'Laderas',false);

Map.addLayer(Hillshade, {min:25,max:230}, "Hillshade");

Map.addLayer(MDT, {min: 19, max :48, opacity: 0.8,
palette: ['#b0b0b0','#a1a1a1','#98b581','#72a890','#7c8ead','#8c75a0','#b47ba1','#cb8b8f',
'#c5a58a','#bdbf89','#8dc458','#3dab71','#5078b6','#77479d','#c04d9c','#e76f7a','#e2a66c',
'#d6db5e','#84d600','#00ab44','#0068c0','#6c00a3','#ca009c','#ff5568','#ffab47','#f4fa00'
]}, 'Pendiente-Orientacion',true);

// Incorporacion de masas de agua sobre MDT combinado de pendiente y ladera
var MasasAgua = ee.Image('JRC/GSW1_3/GlobalSurfaceWater').select('occurrence');
Map.addLayer (MasasAgua, {min: 0.0, max: 100,
  palette: ['#8ba5ff']},
  'Masas de agua',true);

Map.setCenter (-114.39,36.13, 11);
print('GIS AND BEERS:  WWW.GISANDBEERS.COM');