Агрегировать почасовые данные GPM к ежедневным в GEE - PullRequest
0 голосов
/ 25 мая 2019

Мне нужно агрегировать сумму 48 получасовых изображений в день из коллекции GPM, получая коллекцию imageCollection с полосой «осадков» и дневными изображениями

Я пытаюсь заполнить и повторить пустую коллекцию FeatureCollection, но я получаю пустую коллекцию без изображений

var dataset = ee.ImageCollection('NASA/GPM_L3/IMERG_V05')
var startdate = ee.Date.fromYMD(2014,3,1)
var enddate = ee.Date.fromYMD(2014,4,1)
var precipitation = dataset.filter(ee.Filter.date(startdate,enddate)).select('precipitationCal')
print(precipitation)

var difdate = enddate.difference(startdate, 'day')

// Time lapse
var lapse = ee.List.sequence(0, difdate.subtract(1))
var startdate = ee.Date('2014-01-01')
var listdates = lapse.map(function(day){
  return startdate.advance(day, 'day')
})

var pts = ee.FeatureCollection(ee.List([]))
var newft = ee.FeatureCollection(listdates.iterate(function(img, ft) {
  // Cast
  ft = ee.FeatureCollection(ft)
  var day = ee.Date(img)
  // Filter the collection in one day
  var day_collection = precipitation.filterDate(day, day.advance(1, 'day'))
  // Get the sum of all 24 images into one Image
  var sum = ee.Image(day_collection.sum())
  // Return the FeatureCollection with the new properties set
  return sum
}, listdates))

1 Ответ

0 голосов
/ 16 июля 2019

Пожалуйста, попробуйте мой пакет pkg_trend .aggregate_prob функция в нем, работает так же, как aggregate на языке R .

var imgcol_all = ee.ImageCollection('NASA/GPM_L3/IMERG_V05');

function add_date(img){
    var date  = ee.Date(img.get('system:time_start'));
    var date_daily = date.format('YYYY-MM-dd');
    return img.set('date_daily', date_daily);
}

var startdate = ee.Date.fromYMD(2014,3,1);
var enddate   = ee.Date.fromYMD(2014,4,1);
var imgcol = imgcol_all
    .filter(ee.Filter.date(startdate,enddate)).select('precipitationCal')
    .map(add_date);
// imgcol = pkg_trend.imgcol_addSeasonProb(imgcol); 
print(imgcol.limit(3), imgcol.size());

var pkg_trend = require('users/kongdd/public:Math/pkg_trend.js');
var imgcol_daily = pkg_trend.aggregate_prop(imgcol, "date_daily", 'sum');
print(imgcol_daily);

Map.addLayer(imgcol_daily, {}, 'precp daily');

Ссылка GEE https://code.earthengine.google.com/2e04ad4a4bee6789af23bfac42f63025

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...