//VISUALIZACION DE MAPAS DE CONTAMINACION NO2 TROPOSFERICO
//http://www.gisandbeers.com/montar-mapas-de-contaminacion-ante-coronavirus/

//Llamamos a la coleccion de datos Sentinel 5P y seleccionamos la banda de NO2 troposferico
var TROPOMISensor = ee.ImageCollection('COPERNICUS/S5P/NRTI/L3_NO2')
  .select('tropospheric_NO2_column_number_density');

//Definimos los dos momentos de analisis de concentracion de contaminacion
var NO2Previo = TROPOMISensor.filterDate('2020-03-09', '2020-03-15').mean(); //Primer momento
var NO2Posterior = TROPOMISensor.filterDate('2020-03-16', '2020-03-20').mean(); //Segundo momento

//Dibujamos en el visor la zona AOI con ayuda de las herramientas de edición geométrica

//Simbolizamos la rampa de valores de NO2 troposferico con valores maximos-minimos y transparencia
var SimbologiaNO2  = {min:0, max: 0.0002, opacity: 0.70,
  palette: ["black", "blue", "purple", "cyan", "green", "yellow", "red"]};

//Definimos un nuevo mapa base del territorio, por ejemplo, usando el mosaico SRTM
//A partir del DEM SRTM definimos una nueva variable para generar un mapa de sombras
var SRTM = ee.Image('CGIAR/SRTM90_V4');
var HillShade = ee.Terrain.hillshade(SRTM);

//Añadimos un fondo de mapa en color negro para generar contraste de colores
Map.addLayer (ee.Image(0).byte(), {}, "Background");

//Simbolizamos el mapa de sombras hillshade con un degradado de grises y lo añadimos a la vista
Map.addLayer (HillShade, {min: 100, max:255}, 'Hillshade');

//Añadimos el mapa de contaminacion previo al visor de Earth Engine y aplicamos su simbologia 
Map.addLayer (NO2Previo, SimbologiaNO2, 'Emisiones Previas de NO2 troposferico');

//Iniciamos la comparacion panelando la vista con swipe
//Añadimos el mapa de contaminacion del momento posterior y los mapas base que deseamos visualizar
var Comparador = ui.Map();
Comparador.addLayer (ee.Image(0).byte(), {}, "Background");
Comparador.addLayer (HillShade, {min: 100, max:255}, 'Hillshade');
Comparador.addLayer (NO2Posterior, SimbologiaNO2, 'Emisiones Posteriores de NO2 troposferico');
var SWIPE = ui.Map.Linker([ui.root.widgets().get(0), Comparador]);
var SWIPE2 = ui.SplitPanel({
  firstPanel: SWIPE.get(0),
  secondPanel: SWIPE.get(1),
  orientation: 'horizontal', 
  wipe: true,
  style: {stretch: 'both'}});
ui.root.widgets().reset([SWIPE2]);

//Descargamos las imagenes de ambos momentos en Google Drive a resolucion 0,01 grado (1100m)
//Clipeamos y recortamos la zona territorial por los limites geometricos de la zona AOI
var NO2PrevioAOI= NO2Previo.clip (geometry);
var NO2PosteriorAOI= NO2Posterior.clip (geometry);

Export.image.toDrive({
  image: NO2PrevioAOI.select("tropospheric_NO2_column_number_density"),
  description: 'NO2Previo',
  scale: 1100,
  region: geometry});
Export.image.toDrive({
  image: NO2PosteriorAOI.select("tropospheric_NO2_column_number_density"),
  description: 'NO2Posterior',
  scale: 1100,
  region: geometry});

//Centramos la vista en las coordenadas y zoom de la zona que deseemos visualizar
Map.setCenter (-3.72,39.70, 6);