Google Earth Engine - подсчитывает дни самого длинного сухого периода за период - PullRequest
Новая
       8

Google Earth Engine - подсчитывает дни самого длинного сухого периода за период

0 голосов
/ 24 октября 2018

Я пытаюсь сопоставить максимальное количество дней подряд с дождем <1 мм в Google Earth Engine.</p>

Это ссылка на код

https://code.earthengine.google.com/22b5c20d2700a2ffb5989f892838ac58

Сначала я переклассифицирую коллекцию с 0, если дождь <= 1, и 1, если> 1.Затем я запускаю код, который должен считать дни самого длинного сухого периода, но он может сделать это только в том случае, если сухой период достигнет конца периода времени.

Например, если я ищусамый длинный сухой период за 4 дня я получаю следующие серии: 

rain days 1 2 3 4   output
          0,0,1,1 = 0 dry days
          0,1,0,0 = 2 dry days
0 = rain<=1 and 
1 = rain>1 (as per the first step)

Кто-нибудь может помочь в решении этой проблемы?Спасибо

1 Ответ

0 голосов
/ 02 ноября 2018

Я не думаю, что вы были далеко в вашем коде, который вы предоставили.Чтобы отслеживать сухие заклинания, вы должны использовать .iterate().Я попробовал ваше приложение немного по-другому, и вместо того, чтобы классифицировать данные перед итерацией, я вычисляю, какие пиксели являются сухими каждый день, и переносит накопленные дни, когда сухой пиксель, в противном случае он устанавливается равным нулю: 

// DATA
var collection = ee.ImageCollection("UCSB-CHG/CHIRPS/DAILY");

// Define time range
var startyear = 2000;
var endyear = 2017;

var startmonth = 1;
var endmonth = 12;

// Set date in ee date format
var startdate = ee.Date.fromYMD(startyear,startmonth,1);
var enddate = ee.Date.fromYMD(endyear,endmonth,31);

// Filter data
var datain_t = collection.filterDate(startdate, enddate)
  .filter(ee.Filter.calendarRange(startmonth,endmonth, 'month'))
  .select("precipitation").map(function(img){
     return img.addBands(ee.Image.constant(0).uint8().rename('counter'));
  })
  .sort('system:time_start');

// // START 
var dataset = datain_t
.filterDate("2016-08-01","2016-08-30")
.sort('system:time_start:');
print(dataset,"dataset");

var precipThresh = 1; // mm

function drySpells(img, list){
  // get previous image
  var prev = ee.Image(ee.List(list).get(-1));
  // find areas gt precipitation threshold (gt==0, lt==1)
  var dry = img.select('precipitation').lt(precipThresh);
  // add previous day counter to today's counter
  var accum = prev.select('counter').add(dry).rename('counter');
  // create a result image for iteration
  // precip < thresh will equal the accumulation of counters
  // otherwise it will equal zero
  var out = img.select('precipitation').addBands(
        img.select('counter').where(dry.eq(1),accum)
      ).uint8();
  return ee.List(list).add(out);
}

// create first image for iteration
var first = ee.List([ee.Image(dataset.first())]);

// apply dry speall iteration function
var maxDrySpell = ee.ImageCollection.fromImages(
    dataset.iterate(drySpells,first)
).max(); // get the max value

// display results
Map.addLayer(maxDrySpell.select('counter'),{min:0,max:30,palette:'#9ecae1,#ffffff,#ffeda0,#feb24c,#f03b20'},'Max Dry Spells');

Вот ссылка на код: https://code.earthengine.google.com/80b4c0f7e82a5f0da316af1d2a55dd59

Не пытайтесь запустить этот анализ в течение слишком длительного периода времени, иначе Earth Engine выдаст ошибку.Надеюсь, это поможет!

...