Exportar a Google Drive valores mensuales de precipitación obtenidos en GEE a partir de imágenes CHIRPS daily

En un post anterior, se consideró la exportación a Google Drive de una lista con los valores diarios de precipitación obtenidas para una reducción puntual de una colección de imágenes diarias de CHIRPS (Climate Hazards Group InfraRed Precipitation with Station Data). Esta lista puede ser manipulada en Python para obtener los promedios mensuales y anuales para los años considerados en la serie.

Con el fin de corroborar la validez de la síntesis mensual de datos producidos por GEE, se copió, adaptó y ejecutó un código obtenido de aquí. El código finalmente producido fue el siguiente:

var fg_points = ee.Geometry.Point (-66.00, 9.22); // Valle La Pascua

// Set years and month
var startYear = 2012;
var endYear = 2016;
var years = ee.List.sequence(startYear, endYear);
var months = ee.List.sequence(1,12);
// load the image collection
var Daily = ee.ImageCollection("UCSB-CHG/CHIRPS/DAILY");

// make monthly summed mosaics
// loop over the years and months to get summed monthly images
var byMonth = ee.ImageCollection(ee.FeatureCollection(years.map(function(y){
  var yearCollection = Daily.filter(ee.Filter.calendarRange(y, y, 'year'));
  var byYear = ee.ImageCollection.fromImages(
    months.map(function(m) {
      var summedImage = yearCollection.filter(ee.Filter.calendarRange(m, m, 'month'))
        .reduce(ee.Reducer.sum()); 
      var date = ee.Date.fromYMD(y, m, 1).format("MM_dd_YYYY");
      return summedImage.set('system:time_start', ee.Date.fromYMD(y, m, 1)).rename(date);
      //.set('month', m).set('year', y); // eventually set year and month 
  }));
  return byYear;
})).flatten());

// filter the empty one out
var outputMonthly = byMonth.filter(ee.Filter.listContains('system:band_names', 'constant').not())
  .sort('system:time_start').toBands();
print(outputMonthly);

var features = outputMonthly.reduceRegions(fg_points, ee.Reducer.first(), 30);
//print(features);

Export.table.toDrive(features,
"PrecipitationHTL",   //task name
"GEE_Folder",         //folder name
"PrecipitationHTL");  //file name

Después de bajar de Google Drive el archivo correspondiente, se observó que el orden de los meses estaba alterado (para series más voluminosas también ocurre lo mismo para los años). Después de corregido el orden, se computaron los promedios mensuales y totales y éste concordó con el obtenido con los valores diarios. El código anterior es fiable.

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