// VISOR DE IMAGENES SATELITE DE LA ISLA DE LA PALMA
// http://www.gisandbeers.com/visor-imagenes-satelite-la-palma-google-earth-engine/

// MAPA BASE Y AOI LA PALMA
var AOI = ee.Geometry.Rectangle([-18.03, 28.867, -17.685, 28.43]);
var AOILava = ee.Geometry.Rectangle([-17.94, 28.59, -17.85, 28.64]);
var Todoque = ee.Geometry.Point([-17.92, 28.62]);
var SRTM = ee.Image ('CGIAR/SRTM90_V4');
var HillShade = ee.Terrain.hillshade (SRTM);
Map.addLayer (ee.Image(0).byte(), {}, "Background");
Map.addLayer (HillShade, {min: 100, max:255}, 'Hillshade',false);
var Base = ee.Image ('COPERNICUS/S2_SR/20210826T120319_20210826T120321_T28RBS').clip(AOI);
Map.addLayer (Base, {max: 4000.0, min: 200.0, gamma: 1.0, 
     bands: ['B4', 'B3', 'B2']}, 'La Palma',false);
var TituloMapa = ui.Label({
  value: 'ERUPCIÓN DEL VOLCÁN CUMBRE VIEJA, LA PALMA', 
  style: {fontWeight: 'bold', fontSize: '20px', position: 'bottom-right'}}); 
Map.add(TituloMapa);
Map.setCenter (-17.92, 28.65, 11); 


//COLECCIONES DATOS RADAR SENTINEL 1
var Radar = ee.ImageCollection('COPERNICUS/S1_GRD')
.filter(ee.Filter.eq('instrumentMode', 'IW')).select('VH')
.filter(ee.Filter.eq('orbitProperties_pass', 'ASCENDING')) 
.filterMetadata('resolution_meters', 'equals' , 10) 
.filterBounds(AOI)
.map(function(ImagenRadar) {var BordeTile = ImagenRadar.lt(-30.0);
var Mascara = ImagenRadar.mask().and(BordeTile.not());
return ImagenRadar.updateMask(Mascara);});
var RadarT0 = ee.Filter.date('2021-09-01', '2021-09-16');
var T0VH = ee.Image.cat(Radar.filter(RadarT0).mean()).clip(AOI);
Map.addLayer(T0VH, {min:-25,max:5, gamma: 1}, 'Radar pre',false);
var RadarT1 = ee.Filter.date('2021-09-16', '2021-10-01');
var T1VH = ee.Image.cat(Radar.filter(RadarT1).mean()).clip(AOI);
var RadarT2 = ee.Filter.date('2021-10-02', '2021-10-14');
var T2VH = ee.Image.cat(Radar.filter(RadarT2).mean()).clip(AOI);
Map.addLayer(T2VH, {min:-25,max:5, gamma: 1}, 'Radar post',false);
var Coladas = T0VH.addBands(T1VH).addBands(T2VH);
Map.addLayer(Coladas, {min:-25,max:5,gamma: 1}, 'Coladas',false);


//COLECCIONES OPTICAS SENTINEL 2 Y LANDSAT 8
var ColeccionSentinel = ee.ImageCollection ('COPERNICUS/S2')
  .filterDate ('2021-09-16' ,'2021-10-16')
  .filterBounds (Todoque);
//print (ColeccionSentinel);

var S20Septiembre = ee.Image ('COPERNICUS/S2_SR/20210920T120331_20210920T120326_T27RYM').clip(AOI);
Map.addLayer (S20Septiembre, {max: 5000.0, min: 0.0, gamma: 0.8, 
     bands: ['B12', 'B11', 'B4']}, 'S2 20 sept',false);
var S25Septiembre = ee.Image ('COPERNICUS/S2_SR/20210925T120319_20210925T120321_T27RYM').clip(AOI);
Map.addLayer (S25Septiembre, {max: 5000.0, min: 0.0, gamma: 0.8, 
     bands: ['B12', 'B11', 'B4']}, 'S2 25 sept',false);
var L26Septiembre = ee.Image ('LANDSAT/LC08/C01/T1_RT/LC08_208040_20210926').clip(AOI);
Map.addLayer (L26Septiembre, {max: 35000, min: 5000.0, gamma: 1.0,
     bands: ['B7', 'B6', 'B4']}, 'L8 26 sept',false);
var S30Septiembre = ee.Image ('COPERNICUS/S2/20210930T120331_20210930T120327_T27RYM').clip(AOI);
Map.addLayer (S30Septiembre, {max: 5000.0, min: 0.0, gamma: 0.8, 
     bands: ['B12', 'B11', 'B4']}, 'S2 30 sept',false);
var S5Octubre = ee.Image ('COPERNICUS/S2/20211005T120329_20211005T120323_T28RBS').clip(AOI);
Map.addLayer (S5Octubre, {max: 5000.0, min: 0.0, gamma: 0.8, 
     bands: ['B12', 'B11', 'B4']}, 'S2 5 Oct',false);
var S10Octubre = ee.Image ('COPERNICUS/S2/20211010T120331_20211010T120328_T28RBS').clip(AOI);
Map.addLayer (S10Octubre, {max: 5000.0, min: 0.0, gamma: 0.8, 
     bands: ['B12', 'B11', 'B4']}, 'S2 10 Oct',false);
var S15Octubre = ee.Image ('COPERNICUS/S2/20211015T120329_20211015T120324_T28RBS').clip(AOI);
Map.addLayer (S15Octubre, {max: 5000.0, min: 0.0, gamma: 0.8, 
     bands: ['B12', 'B11', 'B4']}, 'S2 15 Oct',true);


// COLECCION S02 SENTINEL 5P
var SO2 = ee.ImageCollection('COPERNICUS/S5P/NRTI/L3_SO2')
  .select ('SO2_column_number_density').filterDate ('2021-09-16', '2021-10-13').mean(); 
Map.addLayer (SO2, {min:0, max: 0.009,
  palette: ['042333', '2c3395', '744992', 'b15f82', 'eb7958', 'fbb43d', 'e8fa5b']}, 'Emisiones SO2', false);


// DESCARGA DE IMAGENES
var IMG15Octubre = S15Octubre.getThumbURL({
  'min': 0, 'max': 5000, 'dimensions': 3000, 'crs': 'EPSG:3857',
  'bands': ['B12', 'B11', 'B4'],
  'region': ee.Geometry.Rectangle([-18.03, 28.867, -17.685, 28.43])});
print('Dia 15 octubre (S2)', IMG15Octubre);
var IMG10Octubre = S10Octubre.getThumbURL({
  'min': 0, 'max': 5000, 'dimensions': 3000, 'crs': 'EPSG:3857',
  'bands': ['B12', 'B11', 'B4'],
  'region': ee.Geometry.Rectangle([-18.03, 28.867, -17.685, 28.43])});
print('Dia 10 octubre (S2)', IMG10Octubre);
var IMG5Octubre = S5Octubre.getThumbURL({
  'min': 0, 'max': 5000, 'dimensions': 3000, 'crs': 'EPSG:3857',
  'bands': ['B12', 'B11', 'B4'],
  'region': ee.Geometry.Rectangle([-18.03, 28.867, -17.685, 28.43])});
print('Dia 5 octubre (S2)', IMG5Octubre);
var IMG30Septiembre = S30Septiembre.getThumbURL({
  'min': 0, 'max': 5000, 'dimensions': 3000, 'crs': 'EPSG:3857',
  'bands': ['B12', 'B11', 'B4'],
  'region': ee.Geometry.Rectangle([-18.03, 28.867, -17.685, 28.43])});
print('Dia 30 septiembre (S2)', IMG30Septiembre);
var IMG26Septiembre = L26Septiembre.getThumbURL({
  'min': 5000, 'max': 35000, 'dimensions': 3000, 'crs': 'EPSG:3857',
  'bands': ['B7', 'B6', 'B4'],
  'region': ee.Geometry.Rectangle([-18.03, 28.867, -17.685, 28.43])});
print('Dia 26 septiembre (L8)', IMG26Septiembre);
var IMG25Septiembre = S25Septiembre.getThumbURL({
  'min': 0, 'max': 5000, 'dimensions': 3000, 'crs': 'EPSG:3857',
  'bands': ['B12', 'B11', 'B4'],
  'region': ee.Geometry.Rectangle([-18.03, 28.867, -17.685, 28.43])});
print('Dia 25 septiembre (S2)', IMG25Septiembre);
var IMG20Septiembre = S20Septiembre.getThumbURL({
  'min': 0, 'max': 5000, 'dimensions': 3000, 'crs': 'EPSG:3857',
  'bands': ['B12', 'B11', 'B4'],
  'region': ee.Geometry.Rectangle([-18.03, 28.867, -17.685, 28.43])});
print('Dia 20 septiembre (S2)', IMG20Septiembre);


// TIMELAPSE
var Timelapse = {crs: 'EPSG:3857', dimensions: '600', region: AOILava, bands: ['B12', 'B11', 'B4'], min: 0.0,  max: 5000, framesPerSecond: 2,};
var Animacion = ui.Thumbnail({image: ColeccionSentinel,  params: Timelapse,  style: {position: 'bottom-right', width: '400px'}});
Map.add(Animacion);


// COMPARATIVA SWIPE
var Comparador = ui.Map();
Comparador.addLayer (ee.Image(0).byte(), {}, "Background");
Comparador.addLayer (Base, {max: 4000.0, min: 200.0, gamma: 1.0, bands: ['B4', 'B3', 'B2']}, 'La Palma',true);
Comparador.addLayer(T0VH, {min:-25,max:0, gamma: 0.9}, 'Radar pre',false);
Comparador.addLayer(T2VH, {min:-25,max:0, gamma: 0.9}, 'Radar post',false);
Comparador.addLayer(Coladas, {min:-25,max:0, gamma: 0.9}, 'Coladas',false);
Comparador.addLayer (S20Septiembre, {max: 5000.0, min: 0.0, gamma: 0.8, 
     bands: ['B12', 'B11', 'B4']}, 'S2 20 sept',false);
Comparador.addLayer (S25Septiembre, {max: 5000.0, min: 0.0, gamma: 0.8, 
     bands: ['B12', 'B11', 'B4']}, 'S2 25 sept',false);
Comparador.addLayer (L26Septiembre, {max: 35000, min: 5000.0, gamma: 1.0,
     bands: ['B7', 'B6', 'B4']}, 'L8 26 sept',false);
Comparador.addLayer (S30Septiembre, {max: 5000.0, min: 0.0, gamma: 0.8, 
     bands: ['B12', 'B11', 'B4']}, 'S2 30 sept',false);
Comparador.addLayer (S5Octubre, {max: 5000.0, min: 0.0, gamma: 0.8, 
     bands: ['B12', 'B11', 'B4']}, 'S2 5 Oct',false);
Comparador.addLayer (S10Octubre, {max: 5000.0, min: 0.0, gamma: 0.8, 
     bands: ['B12', 'B11', 'B4']}, 'S2 10 Oct',false);
Comparador.addLayer (S15Octubre, {max: 5000.0, min: 0.0, gamma: 0.8, 
     bands: ['B12', 'B11', 'B4']}, 'S2 15 Oct',false);
Comparador.addLayer (SO2, {min:0, max: 0.009,
  palette: ['042333', '2c3395', '744992', 'b15f82', 'eb7958', 'fbb43d', 'e8fa5b']}, 'Emisiones SO2', false);
var SWIPE = ui.Map.Linker([ui.root.widgets().get(0), Comparador,]);
var SWIPE2 = ui.SplitPanel({
  firstPanel: SWIPE.get(1),
  secondPanel: SWIPE.get(0),
  orientation: 'horizontal', 
  wipe: true,
  style: {stretch: 'both',  width: '100%'}});
ui.root.widgets().reset([SWIPE2]);
Comparador.setCenter(-17.92, 28.65, 11);