Google Earth Engine(GEE)对指定地点Sentinel-2 Level1C数据进行NDVI指数的计算
程序员文章站
2022-05-26 19:10:26
...
Google Earth Engine (GEE) 是由谷歌公司开发的众多应用之一。借助谷歌公司超强的服务器运算能力以及与NASA的合作关系,GEE平台将Landsat/Sentinel等可以公开获取的遥感图像数据存储在谷歌的磁盘阵列中,使得GEE用户可以方便的提取、调用和分析海量的遥感大数据资源。
本文基于GEE平台,实现对指定地点Sentinel-2 Level1C数据进行NDVI指数的计算,示例地点为湖南省长沙市。代码如下:
(JS)
var cs = China.filterBounds(point).geometry();
var roi = cs;
function NDVI_V3(img) {
var ndvi = img.clip(roi).normalizedDifference(["B8","B4"]);
return ndvi;
}
//使用QA波段去云
function maskS2clouds(image) {
var qa = image.select('QA60');
// Bits 10 and 11 are clouds and cirrus, respectively.
var cloudBitMask = 1 << 10;
var cirrusBitMask = 1 << 11;
// Both flags should be set to zero, indicating clear conditions.
var mask = qa.bitwiseAnd(cloudBitMask).eq(0)
.and(qa.bitwiseAnd(cirrusBitMask).eq(0));
return image.updateMask(mask).divide(10000);
}
//sentinel2 and roi
var dataset = ee.ImageCollection("COPERNICUS/S2")
.filterDate("2018-01-01", "2018-12-31")
.filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 20))
.map(maskS2clouds);
var img = dataset.median().clip(roi)
//var img = ee.Image(dataset.filterBounds(roi).first());
print(img)
var ndvi3 = NDVI_V3(img);
var visParam = {
min: -0.2,
max: 0.8,
palette: 'FFFFFF, CE7E45, DF923D, F1B555, FCD163, 99B718, 74A901, 66A000, 529400,' +
'3E8601, 207401, 056201, 004C00, 023B01, 012E01, 011D01, 011301'
};
var rgbVis = {
min: 0.0,
max: 3048,
bands: ['B8', 'B4', 'B3'],
};
Map.addLayer(img, rgbVis, "raw_img");
Map.addLayer(ndvi3, visParam, "ndvi_3");
Map.centerObject(roi, 9);
注:①使用前,需要import我国的行政区划文件,其次需要选取例如湖南省长沙市境内的一点(point),最后import导入S2图 像。
②关于我国行政区划文件,分享百度网盘:链接:https://pan.baidu.com/s/1vrT5AmRTX8KmMMbpGtTAcg
提取码:o5vt
效果:
参考文章: