基于 GEE 下载 1958—2023 年全球 0.04 度气候数据(蒸散、气温、干旱、降水、土壤湿度、辐射、饱和水汽压差、风速)

目录

1 完整代码

2 运行结果


1 完整代码

根据需要自行修改“起始年份”“截止年份”即可。

var geometry = table;
Map.centerObject(geometry, 6);

var terraclimate = ee.ImageCollection("IDAHO_EPSCOR/TERRACLIMATE")
  .filterBounds(geometry);

var startYear = 1958;
var endYear = 1960;

var sumBands = ['aet', 'pet', 'pr', 'srad'];
var meanBands = ['pdsi', 'soil', 'vpd', 'vs'];

for (var year = startYear; year <= endYear; year++) {
  var yearCollection = terraclimate.filterDate(year + '-01-01', year + '-12-31');

  var maxTempMean = yearCollection.select('tmmx').mean().multiply(0.1);
  var minTempMean = yearCollection.select('tmmn').mean().multiply(0.1);
  var annualMeanTemp = maxTempMean.add(minTempMean).divide(2);

  var sumImage = yearCollection.select(sumBands).sum()
    .multiply(ee.Image([0.1, 0.1, 1, 0.1]))
    .toDouble()
    .rename(sumBands);

  var meanImage = yearCollection.select(meanBands).mean()
    .multiply(ee.Image([0.01, 0.1, 0.01, 0.01]))
    .toDouble()
    .rename(meanBands);

  var annualImage = sumImage.addBands(meanImage).addBands(annualMeanTemp.rename('tmean'))
    .clip(geometry)
    .set('year', year);
    
    

  Export.image.toDrive({
    image: annualImage.select('aet'),
    description: 'AET_' + year,
    fileNamePrefix: 'AET_' + year,
    folder: 'TerraClimate_Annual',
    scale: 4638.3,
    region: geometry,
    crs: "EPSG:4326",
    maxPixels: 1e13
  });
  
  Export.image.toDrive({
    image: annualImage.select('pet'),
    description: 'PET_' + year,
    fileNamePrefix: 'PET_' + year,
    folder: 'TerraClimate_Annual',
    scale: 4638.3,
    region: geometry,
    crs: "EPSG:4326",
    maxPixels: 1e13
  });


  Export.image.toDrive({
    image: annualImage.select('pr'),
    description: 'Precipitation_' + year,
    fileNamePrefix: 'Precipitation_' + year,
    folder: 'TerraClimate_Annual',
    scale: 4638.3,
    region: geometry,
    crs: "EPSG:4326",
    maxPixels: 1e13
  });

  Export.image.toDrive({
    image: annualImage.select('tmean'),
    description: 'Avg_Temperature_' + year,
    fileNamePrefix: 'Avg_Temperature_' + year,
    folder: 'TerraClimate_Annual',
    scale: 4638.3,
    region: geometry,
    crs: "EPSG:4326",
    maxPixels: 1e13
  });

  Export.image.toDrive({
    image: annualImage.select('pdsi'),
    description: 'PDSI_' + year,
    fileNamePrefix: 'PDSI_' + year,
    folder: 'TerraClimate_Annual',
    scale: 4638.3,
    region: geometry,
    crs: "EPSG:4326",
    maxPixels: 1e13
  });

  Export.image.toDrive({
    image: annualImage.select('soil'),
    description: 'Soil_Moisture_' + year,
    fileNamePrefix: 'Soil_Moisture_' + year,
    folder: 'TerraClimate_Annual',
    scale: 4638.3,
    region: geometry,
    crs: "EPSG:4326",
    maxPixels: 1e13
  });

  Export.image.toDrive({
    image: annualImage.select('srad'),
    description: 'Radiation_' + year,
    fileNamePrefix: 'Radiation_' + year,
    folder: 'TerraClimate_Annual',
    scale: 4638.3,
    region: geometry,
    crs: "EPSG:4326",
    maxPixels: 1e13
  });

  Export.image.toDrive({
    image: annualImage.select('vpd'),
    description: 'VPD_' + year,
    fileNamePrefix: 'VPD_' + year,
    folder: 'TerraClimate_Annual',
    scale: 4638.3,
    region: geometry,
    crs: "EPSG:4326",
    maxPixels: 1e13
  });

  Export.image.toDrive({
    image: annualImage.select('vs'),
    description: 'Wind_Speed_' + year,
    fileNamePrefix: 'Wind_Speed_' + year,
    folder: 'TerraClimate_Annual',
    scale: 4638.3,
    region: geometry,
    crs: "EPSG:4326",
    maxPixels: 1e13
  });
}

2 运行结果

点击RUN即可下载数据
帕默尔干旱指数PDSI_2020
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值