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

Twcat_tree

发布于 2024-07-09 08:08:43

1060

发布于 2024-07-09 08:08:43

文章被收录于专栏：二猫の家

简介

在本篇博客中，我们将使用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'
});