GEE代码实例教程详解：长时间序列风速分析

简介

在本篇博客中，我们将使用Google Earth Engine (GEE) 对长时间序列的风速数据进行分析。通过ERA5再分析数据集，我们可以计算2010年至2024年间的平均风速，并与1980年至2020年的风速数据进行比较。

背景知识

ERA5数据集

ERA5是ECMWF（European Centre for Medium-Range Weather Forecasts，欧洲中期天气预报中心）提供的全球气候再分析数据集，具有高空间和时间分辨率。

风速计算

风速是风向的两个分量（u和v分量）的矢量合成，通常使用平方和的平方根来计算。

完整代码

// 定义研究区域的坐标点
var cor = [[52.24091120672926, 28.94538993989301],[59.00848933172926, 28.94538993989301],[59.00848933172926, 33.89295282767506],[52.24091120672926, 33.89295282767506],[52.24091120672926, 28.94538993989301]
];// 创建多边形区域
var roi = ee.Geometry.Polygon(cor);// 加载ERA5数据集并选择风速分量
var era5 = ee.ImageCollection("ECMWF/ERA5/DAILY").select('u_10m', 'v_10m').filterDate('2010', '2024');// 计算风速
var windspeed = era5.map(function (img) {var u = img.select('u_10m');var v = img.select('v_10m');var wind_speed = u.pow(2).add(v.pow(2)).sqrt().rename('wind_speed');return wind_speed.copyProperties(img, img.propertyNames());
});// 绘制风速时间序列图表
print(ui.Chart.image.series(windspeed, roi, ee.Reducer.mean(),27000, 'system:time_start')
);// 加载1980年至2020年的ERA5数据集
var wind40 = ee.ImageCollection("ECMWF/ERA5/DAILY").select('u_10m', 'v_10m').filterDate('1980', '2020');// 计算风速
var windspeed40 = wind40.map(function (img) {var u = img.select('u_10m');var v = img.select('v_10m');var wind_speed = u.pow(2).add(v.pow(2)).sqrt().rename('wind_speed');return wind_speed.copyProperties(img, img.propertyNames());
});// 定义函数提取风速平均值和日期
function get_value(img) {var date = ee.Date(img.get('system:time_start')).format().slice(0, 10);var value = img.reduceRegion({reducer: ee.Reducer.mean(),geometry: roi,scale: 27000}).get('wind_speed');return [date, value];
}// 提取风速数据列表
var size = windspeed40.size();
var windlist = windspeed40.toList(size).map(get_value);// 将风速数据转换为特征集合
var wind_feature = ee.FeatureCollection(windlist.map(function (el) {el = ee.List(el);return ee.Feature(null, {'date': el.get(0),'windspeed': el.get(1).getInfo()});
}));// 导出风速数据表至Google Drive
Export.table.toDrive({collection: wind_feature,description: 'windspeed40',fileFormat: 'CSV'
});