本文分享如何利用自己创建表达式计算影像集,计算表达式为NDVI和SAVI
//筛选影像数据集
ar collection = ee.ImageCollection('LANDSAT/LE07/C01/T1_TOA')
.filterDate('2002-11-01', '2002-12-01');
//创建一个计算影像集的表达式函数_NDVI,选取波段名必须与影像波段名匹配
var NDVI = function(image) {
return image.expression('float(b("B4") - b("B3")) / (b("B4") + b("B3"))');
};
// 创建另一个表达式函数
var SAVI = function(image) {
return image.expression(
'(1 + L) * float(nir - red)/ (nir + red + L)',
{
'nir': image.select('B4'),
'red': image.select('B3'),
'L': 0.2
});
};
// 创建一个用于影像显示的调色板
var vis = {
min: 0,
max: 1,
palette: [
'FFFFFF', 'CE7E45', 'DF923D', 'F1B555', 'FCD163', '99B718',
'74A901', '66A000', '529400', '3E8601', '207401', '056201',
'004C00', '023B01', '012E01', '011D01', '011301'
]
};
Map.setCenter(-93.7848, 30.3252, 11);
// 图像显示
Map.addLayer(collection.map(NDVI).mean(), vis, 'Mean NDVI');
Map.addLayer(collection.map(SAVI).mean(), vis, 'Mean SAVI');
代码运行结果
SAVI显示
NDVI显示