Agregación de series de tiempo de precipitación en GEE

La agregación temporal de series de tiempo en Google Earth Engine (GEE) se hace con base en la clase ‘ee.Reducer’. La estadística simple que se usa en la agregación puede ser del tipo minimum, maximum, mean, median, standard deviation, sum, histogram, linear regression, list, etc. En el caso de la precipitación, se va a disponer de series diarias adquiridas de Climate Hazards Group InfraRed Precipitation with Station (CHIRPS) con 0.05° de resolution. Para producir un sólo mapa anual (2017 en este caso) la serie de 365 rásters anuales se agregará con sum.

En el código siguiente se accede a la ImageCollection mediante el id “UCSB-CHG/CHIRPS/DAILY” y se extrae para el área de Venezuela (continental e insular) mediante la geometría de su FeatureCollection filtrado de “USDOS/LSIB/2013”. La bounding box de Venezuela también se emplea para exportar el mapa resultante a Google Drive. Se usó la escala por defecto de Google Drive (asignó el valor de 0.00898315°; aproximadamente 1 km) porque la original (0.05°) producía errores de exportación (probablemente me falta asignar un parámetro para indicarle que la exportación no es en metros).

//Selecting International Boundary Polygons
var countries = ee.FeatureCollection("USDOS/LSIB/2013");

//Selecting Boundary Polygon of Venezuela
var venezuela = countries.filter(ee.Filter.inList('cc',['VE']));

var bounds = venezuela.geometry().bounds();

var geometry = ee.Geometry.Rectangle([-74.769, 0.084, -59.52, 14.389]);

print(bounds, geometry);

////Selecting Image Collection of dayly precipitation
var precipCollection = ee.ImageCollection("UCSB-CHG/CHIRPS/DAILY")
  .filterDate('2017-01-01', '2017-12-31');

print(precipCollection);

// reduce the image collection to one image by summing the 365 daily rasters
var annualPrecip = precipCollection.reduce(ee.Reducer.sum());

var colorizedVis = {
  min: 0.0,
  max: 4000.0,
  palette: ['dee8e9', 'b8fff5', '77ebda','22ab98', '0047e6'],
  
};

var annualPrecip = annualPrecip.clip(venezuela.geometry())

Map.addLayer(annualPrecip, colorizedVis, 'annualPrecip');

// Export the image, specifying scale and region.
Export.image.toDrive({
  image: annualPrecip,
  region:bounds,
  description: 'precipitation_yearly'
});

La ejecución del código anterior produjo la visualización siguiente en GEE:

reducers1

El mismo mapa visualizado en QGIS:

reducers2

Anuncios
Esta entrada fue publicada en Google Earth Engine. Guarda el enlace permanente.

Responder

Por favor, inicia sesión con uno de estos métodos para publicar tu comentario:

Logo de WordPress.com

Estás comentando usando tu cuenta de WordPress.com. Cerrar sesión /  Cambiar )

Google+ photo

Estás comentando usando tu cuenta de Google+. Cerrar sesión /  Cambiar )

Imagen de Twitter

Estás comentando usando tu cuenta de Twitter. Cerrar sesión /  Cambiar )

Foto de Facebook

Estás comentando usando tu cuenta de Facebook. Cerrar sesión /  Cambiar )

Conectando a %s