// 定义研究区域的坐标点
var cor = [
[44.44379054392027, 34.879626513916506],
[49.27777491892027, 34.879626513916506],
[49.27777491892027, 38.95343252297563],
[44.44379054392027, 38.95343252297563],
[44.44379054392027, 34.879626513916506]
];
// 创建多边形区域
var roi = ee.Geometry.Polygon(cor);
// 将地图中心设置为研究区域
Map.centerObject(roi);
// 定义时间范围
var time_start = '2010';
var time_end = '2015';
// 加载MODIS NDSI数据集
var modis = ee.ImageCollection("MODIS/061/MOD10A1")
.select('NDSI_Snow_Cover')
.filterDate(time_start, time_end);
// 加载MODIS土地覆盖数据集并提取水体
var lc = ee.ImageCollection("MODIS/061/MCD12Q1").select('LC_Type1').mode();
var wb = lc.eq(17).not();
// 将水体数据裁剪至研究区域并添加至地图
Map.addLayer(wb.clip(roi), [], 'water_body', false);
// 计算NDSI最大值并更新掩膜
var ndsi_max = modis.max().updateMask(wb);
// 将NDSI最大值裁剪至研究区域并添加至地图
Map.addLayer(ndsi_max.clip(roi), [], 'ndsi_2010', false);
// 设置阈值并创建掩膜
var thr = ndsi_max.gt(50);
Map.addLayer(thr.clip(roi), [], 'thr50', false);
// 更新掩膜
var mask = thr.updateMask(thr);
Map.addLayer(mask.clip(roi), [], 'mask50', false);
// 计算雪盖面积
var area = mask.multiply(ee.Image.pixelArea().divide(1e6));
Map.addLayer(area.clip(roi), [], 'area', false);
// 计算雪盖总面积
var snow_area = ee.Number(area.reduceRegion({
reducer: ee.Reducer.sum(),
geometry: roi,
scale: 500,
maxPixels: 1e13
}).get(0));
print('snow cover area km2', snow_area);
// 映射NDSI图像并计算雪盖面积时间序列
var ndsi = modis.map(function (img) {
var mask_water = img.updateMask(wb);
var ndsi_thr = mask_water.gt(20);
var ndsi_mask = ndsi_thr.updateMask(ndsi_thr);
var ndsi_area = ndsi_mask.multiply(ee.Image.pixelArea().divide(1e6));
return ndsi_area.copyProperties(img, img.propertyNames());
});
print(
ui.Chart.image.series(ndsi, roi, ee.Reducer.sum(),
500, 'system:time_start')
);
// 导出特定时间范围的NDSI图像至Google Drive
Export.image.toDrive({
image: ndsi.filterDate('2011', '2012')
.filter(ee.Filter.calendarRange(1, 1, 'month')).toBands().clip(roi),
description: 'modis_ndsi',
scale: 500,
region: roi,
crs: 'EPSG:4326'
});
- 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.
- 64.
- 65.
- 66.
- 67.
- 68.
- 69.
- 70.
- 71.
- 72.
- 73.
- 74.
- 75.
- 76.
- 77.
- 78.
- 79.
- 80.
- 81.
- 82.