GEE中利用哨兵二号影像计算2018-2020每月NDVI代码实现
时间: 2023-06-19 19:09:09 浏览: 258
GEE上不同影像NDVI的计算公式.txt
5星 · 资源好评率100%
在GEE中,可以使用Sentinel-2影像计算每月的NDVI,并且可以使用哨兵二号影像的数据来实现这个目的。下面是一个计算2018-2020年每月NDVI的代码实现:
```javascript
// 定义计算NDVI函数
var addNDVI = function(image){
var ndvi = image.normalizedDifference(['B8', 'B4']).rename('NDVI');
return image.addBands(ndvi);
};
// 定义日期范围
var startDate = '2018-01-01';
var endDate = '2020-12-31';
var startYear = 2018;
var endYear = 2020;
// 加载哨兵二号影像
var sentinel2 = ee.ImageCollection('COPERNICUS/S2_SR')
.filterBounds(geometry)
.filterDate(startDate, endDate)
.map(addNDVI);
// 将每个月的影像合并成一个影像
var months = ee.List.sequence(1, 12);
var years = ee.List.sequence(startYear, endYear);
var byMonthYear = ee.ImageCollection.fromImages(
years.map(function (y) {
return months.map(function (m) {
var monthImage = sentinel2.filter(ee.Filter.calendarRange(y, y, 'year'))
.filter(ee.Filter.calendarRange(m, m, 'month'))
.reduce(ee.Reducer.mean());
return monthImage.set('year', y).set('month', m).set('system:time_start', ee.Date.fromYMD(y, m, 1));
});
}).flatten());
// 将每个月的NDVI作为一个图层添加到地图上
var ndviPalette = ['FFFFFF', 'CE7E45', 'DF923D', 'F1B555', 'FCD163', '99B718',
'74A901', '66A000', '529400', '3E8601', '207401', '056201',
'004C00', '023B01', '012E01', '011D01', '011301'];
var ndviParams = {min: -0.2, max: 1, palette: ndviPalette};
byMonthYear = byMonthYear.map(function(image){
return image.clip(geometry);
});
Map.addLayer(byMonthYear.select('NDVI'), ndviParams, 'NDVI');
```
这个代码实现首先定义了一个计算NDVI的函数`addNDVI`,然后加载了哨兵二号影像,将每个月的影像合并成一个影像,并将每个月的NDVI作为一个图层添加到地图上。其中,`ndviPalette`是一个NDVI的颜色调色板,用于将NDVI值转换为对应的颜色。
阅读全文