基于 GEE 分析降水量趋势

目录

1 完整代码

2 运行结果

1 完整代码

var roi = table,
    chirps = ee.ImageCollection("UCSB-CHG/CHIRPS/PENTAD");
Map.centerObject(roi,7)
var styling = {color:"red",fillColor:"00000000"};
Map.addLayer(roi.style(styling),{},"geometry")
// 设置起始时间
var startyear = 2015;
var endyear = 2020;
var startdate = ee.Date.fromYMD(startyear,1,1);
var enddate = ee.Date.fromYMD(endyear,12,31)

// 设置年份list
var years = ee.List.sequence(startyear,endyear);

var months = ee.List.sequence(1,12);

// 筛选数据
var Pchirps = chirps.filterDate(startdate, enddate)
  // 按时间顺序降序排列
  .sort('system:time_start', false)
  .filterBounds(roi)
  .select("precipitation");
// 计算每个月的降水量
var MonthlyChirps =  ee.ImageCollection.fromImages(
  years.map(function (y) {
  return months.map(function(m){
    var w = Pchirps.filter(ee.Filter.calendarRange(y, y, 'year'))
             .filter(ee.Filter.calendarRange(m, m, 'month'))
             .sum();
    return w.set('year', y)
             .set('month', m)
             .set('system:time_start',ee.Date.fromYMD(y,m,1)) 
             .set('date', ee.Date.fromYMD(y,m,1))
});
}).flatten());

var mySample = ee.FeatureCollection(MonthlyChirps.map(function(img){
  var myimg = img.reproject('EPSG:32648',null,5000); 
  var time = myimg.get('date');
  var ft =  ee.Image(myimg).sampleRegions
  ({
     collection: roi,
     properties : ['Stations'],
     scale: 5000
  }).map(function(f) {return f.set('date', img.get('date'))})
  return ft;
})).flatten();
// 定义表
var title = {
  title: 'Rainfall',
  hAxis: {title: '时间'},
  vAxis: {title: '降水量(mm)'},
};
var chart = ui.Chart.image.seriesByRegion(
  MonthlyChirps, roi, ee.Reducer.mean(), 'precipitation', 5000, 'system:time_start', 'station')
    .setOptions(title);
// 输出结果
print(chart)

2 运行结果