朋友们好!本次我们对ERA5、FLDAS、CHIRPS和TerraClimate数据集的年均降水量进行了分析。大家只需根据自己的研究区域和时间段进行相应调整,即可获取所需的数据。希望这对大家的研究有所帮助!
一、GEE代码
var ET_region = table;
var collection1 = ee.ImageCollection("ECMWF/ERA5_LAND/MONTHLY_AGGR")
.filterDate('1985-01-01', '2023-12-31')
.select('total_precipitation_sum');
var collection3 = ee.ImageCollection("NASA/FLDAS/NOAH01/C/GL/M/V001")
.filterDate('1985-01-01', '2023-12-31')
.select('Rainf_f_tavg');
var collection4 = ee.ImageCollection("UCSB-CHG/CHIRPS/DAILY")
.filterDate('1985-01-01', '2023-12-31')
.select('precipitation');
var collection5 = ee.ImageCollection("IDAHO_EPSCOR/TERRACLIMATE")
.filterDate('1985-01-01', '2023-12-31')
.select('pr');
function getYearlyData(collection, datasetName, bandName, factor) {
var range = collection.reduceColumns(ee.Reducer.minMax(), ['system:time_start']);
var diff = ee.Date(range.get('max')).difference(ee.Date(range.get('min')), 'year');
var yearlyData = ee.List.sequence(0, diff).map(function(n) {
var start = ee.Date(range.get('min')).advance(n, 'year');
var end = start.advance(1, 'year');
var yearlyImage = collection
.filterDate(start, end)
.sum()
.multiply(factor)
.set('system:time_start', start);
var year = ee.Date(yearlyImage.get('system:time_start')).get('year');
var meanValue = yearlyImage.reduceRegion({
reducer: ee.Reducer.mean(),
geometry: ET_region,
scale: 50000
}).get(bandName);
return ee.Feature(null, {
'Year': year,
'Dataset': datasetName,
'Precipitation (mm)': meanValue
});
});
return ee.List(yearlyData);
}
var yearlyERA5 = getYearlyData(collection1, 'ERA5', 'total_precipitation_sum', 1000);
var yearlyFLDAS = getYearlyData(collection3, 'FLDAS', 'Rainf_f_tavg', 24 * 60 * 60 * 30);
var yearlyCHIRPS = getYearlyData(collection4, 'CHIRPS', 'precipitation', 1);
var yearlyTerraClimate = getYearlyData(collection5, 'TerraClimate', 'pr', 1);
var allYearlyData = ee.List(yearlyERA5)
.cat(yearlyFLDAS)
.cat(yearlyCHIRPS)
.cat(yearlyTerraClimate);
Export.table.toDrive({
collection: ee.FeatureCollection(allYearlyData),
description: 'YearlyPrecipitationData',
fileFormat: 'CSV',
});
二、成为会员