// 定义研究点坐标
var cor = [46.00719637247596, 37.199483887704204];
var loc = ee.Geometry.Point(cor);
// 加载WWF HydroSHEDS流域数据集并筛选研究区域
var roi = ee.FeatureCollection("WWF/HydroSHEDS/v1/Basins/hybas_5").filterBounds(loc);
Map.centerObject(roi);
Map.addLayer(roi);
// 定义时间范围
var time_start = '2001', time_end = '2024';
// 加载MODIS ET和PET数据集
var modis = ee.ImageCollection("MODIS/061/MOD16A2GF")
.select('ET', 'PET')
.filterDate(time_start, time_end);
// 计算平均ET
var et_mean = modis.select('ET').mean().multiply(0.1);
Map.addLayer(et_mean.clip(roi), [], 'et_mean', false);
// 计算夏季平均ET
var et_mean_summer = modis.select('ET')
.filter(ee.Filter.calendarRange(6, 8, 'month')).mean().multiply(0.1);
Map.addLayer(et_mean_summer.clip(roi), [], 'et_summer', false);
// 导出平均ET图像
Export.image.toDrive({
image: et_mean.clip(roi),
description: 'et_mean',
region: roi,
scale: 500,
crs: et_mean.getInfo().crs,
maxPixels: 1e13,
folder: 'evapotranspiration'
});
// 计算CWSI
var cwsi = modis.map(function (img) {
var et = img.select('ET').multiply(0.1);
var pet = img.select('PET').multiply(0.1);
var cwsi = img.expression('1 - (et / pet)', {
'et': et,
'pet': pet
}).rename('cwsi');
return cwsi.copyProperties(img, ['system:time_start', 'system:time_end']);
});
var cwsi_mean = cwsi.mean();
Map.addLayer(cwsi_mean.clip(roi), [], 'cwsi_mean', false);
// 加载MODIS土地覆盖数据集
var lc = ee.ImageCollection("MODIS/061/MCD12Q1")
.select('LC_Type1')
.mode();
// 筛选作物区域
var crop = lc.eq(12);
Map.addLayer(crop.clip(roi), [], 'crop', false);
// 应用掩膜到CWSI
var cwsi_crop = cwsi_mean.updateMask(crop);
Map.addLayer(cwsi_crop.clip(roi), [], 'cwsi_mean_crop', false);
- 1.
- 2.
- 3.
- 4.
- 5.
- 6.
- 7.
- 8.
- 9.
- 10.
- 11.
- 12.
- 13.
- 14.
- 15.
- 16.
- 17.
- 18.
- 19.
- 20.
- 21.
- 22.
- 23.
- 24.
- 25.
- 26.
- 27.
- 28.
- 29.
- 30.
- 31.
- 32.
- 33.
- 34.
- 35.
- 36.
- 37.
- 38.
- 39.
- 40.
- 41.
- 42.
- 43.
- 44.
- 45.
- 46.
- 47.
- 48.
- 49.
- 50.
- 51.
- 52.
- 53.
- 54.
- 55.
- 56.
- 57.
- 58.
- 59.
- 60.
- 61.
- 62.
- 63.