GEE | 多源数据集年降水量对比分析

职场   2025-01-26 17:07   广东  

朋友们好!本次我们对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',});

二、成为会员

GIS遥感数据处理应用
会员:数据处理,ArcGIS/Python/MATLAB/R/GEE教学,指导作图和论文。
 最新文章