【GEE】全球气象数据ERA5风速—波段处理及下载

xtigao

已于 2024-09-10 14:57:03 修改

阅读量904

点赞数 4

分类专栏： GEE 文章标签： arcgis gis google earth

于 2024-09-10 14:49:46 首次发布

本文链接：https://blog.youkuaiyun.com/qq_39546299/article/details/142098875

版权

GEE 专栏收录该内容

4 篇文章

订阅专栏

10m的风速均值是由两个波段u方向风速风量(u_component_of_wind_10m)和v方向风速风量 (v_component_of_wind_10m)开平方求得。
计算era5中的每年每8天的均值并输出。

Alt

var nmg =**********;
  
  // 定义一个函数来计算给定日期范围内的8天风速均值  
function calculateWindSpeed(startDate, endDate) {  
  var era5Subset = ee.ImageCollection("ECMWF/ERA5_LAND/DAILY_AGGR")  
    .select(['u_component_of_wind_10m', 'v_component_of_wind_10m'])  
    .filterDate(startDate, endDate)  
    .filterBounds(nmg);  
  
  var windSpeedImage = era5Subset.map(function(image) {  
    return image.expression(  
      'sqrt(u*u + v*v)',  
      {'u': image.select('u_component_of_wind_10m'), 'v': image.select('v_component_of_wind_10m')}  
    ).rename('wind_speed_10m');  
  }).reduce(ee.Reducer.mean());  
  
  var firstImageDate = ee.Date(era5Subset.first().get('system:time_start'));  
  var dateStr = firstImageDate.format('YYYYMMDD');  
  
  return windSpeedImage.set('eightDayStart', dateStr);  
}  
// 客户端循环来计算每个8天窗口的均值并导出  
for (var i = 2000; i <= 2020; i++) {  

// 定义一个日期数组（这里我们手动设置几个窗口作为示例）  
  var startDateStrings = [i + '-01-01',i + '-01-09',i + '-01-17',i + '-01-25',
i + '-02-02',i + '-02-10',i + '-02-18',i + '-02-26',i + '-03-06',i + '-03-14',
i + '-03-22',i + '-03-30',i + '-04-07',i + '-04-15',i +'-04-23',i + '-05-01',
i + '-05-09',i + '-05-17',i + '-05-25',i + '-06-02',i + '-06-10',i + '-06-18',
i + '-06-26',i + '-07-04',i + '-07-12',i + '-07-20',i + '-07-28',i + '-08-05',
i + '-08-13',i +'-08-21',i + '-08-29',i +'-09-06',i + '-09-14',i + '-09-22',
i + '-09-30',i + '-10-08',i + '-10-16',i + '-10-24',i + '-11-01',i + '-11-09',
i + '-11-17',i + '-11-25',i +'-12-03',i + '-12-11',i +'-12-19',i + '-12-27']; // 可以根据需要扩展  
  

  startDateStrings.forEach(function(startDateStr, index) {  
    var startDate = ee.Date(startDateStr);  
    var endDate = startDate.advance(7, 'day');  
    // 计算8天风速均值  
    var meanWindSpeedImage = calculateWindSpeed(startDate, endDate);  
    // ...（根据需要导出meanRadiationImage和meanWindSpeedImage）  
    var imageId = 'WD_' + i + ('0' + (index + 1)).slice(-2); // 使用年份和三位数的顺序号 
   
    // 导出风速图像到 Google Drive  
    Export.image.toDrive({  
      image: meanWindSpeedImage,  
      description: imageId,  
      fileNamePrefix: imageId,  
      scale: 1000,  
      region: nmg,  
      maxPixels: 1e13,  
      folder: 'GEE_ERA5_LAND_WIND' // 或者与太阳辐射相同的文件夹，根据需要调整  
    });  
  
    // 重置计数器（如果需要在同一年份内重置）  
    // 注意：这里的计数器可能不需要重置，因为它依赖于年份的外部循环  
  });  
}