目录
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 运行结果

