目标:获取某区域的NDVI的影像
结果:在某一时段范围内,有多少期影像,最终就返回多少幅影像,并以tif格式输出
注意事项:同上一篇
参考:
-
(2条消息) GEE学习笔记二 GEE获取Sentine-2与Landsat8 NDVI数据并进行时序变化分析_ndvi时间序列曲线怎么得到_努力的噗噗的博客-CSDN博客,编程 | gee逐期导出研究区NDVI(基于Landsat) - 知乎 (zhihu.com)
代码
// 矢量范围
var geometry = ee.FeatureCollection("users/yuanmengjia9991/dashalong")
// 去云函数
function maskL8sr(image){
// Bits 4 and 3 are cloud shadow and cloud, respectively.
var cloudShadowBitMask = (1 << 4);
var cloudsBitMask = (1 << 3);
// Get the pixel QA band.
var qa = image.select('QA_PIXEL');
// Both flags should be set to zero, indicating clear conditions.
var mask = qa.bitwiseAnd(cloudShadowBitMask).eq(0)
.and(qa.bitwiseAnd(cloudsBitMask).eq(0));
return image.updateMask(mask);
}
// 定义NDVI函数
function NDVI_LS(image){
var ndvi = image.expression(
"(nir-red)/(nir+red)",
{
"nir": image.select("SR_B5").multiply(0.0000275).add(-0.2),
"red": image.select("SR_B4").multiply(0.0000275).add(-0.2)
}).rename('NDVI');
return image.addBands(ndvi)
}
var col = ee.ImageCollection("LANDSAT/LC08/C02/T1_L2")
.filterDate('2016-3-1','2021-10-31')
.filterBounds(geometry)
.filter(ee.Filter.lte('CLOUD_COVER',10))
.map(maskL8sr)
.map(NDVI_LS)
.select('NDVI');
//引入批量导出函数
var batch = require('users/fitoprincipe/geetools:batch')
//导出影像
batch.Download.ImageCollection.toDrive(col, 'yakou_NDVI', {
scale: 30,
region: geometry
})