||
实际应用中我们常常会遇见需要对研究区域的长时序多幅遥感影像分别进行波段运算得到NDVI或者进行其他处理,本质上还是批量进行波段运算(或者其他数学运算)。
参考GEE example code:
// Map an expression over a collection.
// Computes the mean NDVI and SAVI by mapping an expression over a collection
// and taking the mean.
var 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 NDVI = function(image) {
return image.expression('float(nir - red) / (nir + red)',
{
'nir': image.select('B4'),
'red': image.select('B3'),
});
};
//或者这样写 注意到函数名的位置
function add_NDVI (image){
var NDVI = image.normalizedDifference(['B4','B3'])
return image.addBands(NDVI)//波段名为'nd'
}
//或者
var NDVI = function(image) {
return image.addBands(image.normalizedDifference(["B4", "B3"]).rename("NDVI"));//对波段命名
};
// 计算SAVI(Soil Adjusted Vegetation Index)的函数
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 the functions over the collection, reduce to mean and display.
Map.addLayer(collection.map(NDVI).mean(), vis, 'Mean NDVI');
Map.addLayer(collection.map(SAVI).mean(), vis, 'Mean SAVI');
多日平均SAVI
多日平均NDVI
【参考】
https://zhuanlan.zhihu.com/p/129385465
点滴分享,福泽你我!Add oil!
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-9-1 14:35
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社