GEE计算NDVI、NDTI、EVI等指数

GEE计算NDVI、NDTI、EVI等指数

GEE具有海量的遥感数据，所以利用GEE代码在谷歌云平台计算NDVI、NDTI、EVI等指数很方便

首先确定研究区。以山西省为例

代码如下：

var roi = ee.FeatureCollection("users/lilei655123/shanxi");

//归一化热差指数 (NDTI)，用于区分城市地区的人行道和屋顶。

function rmL8Cloud(image) {

var cloudShadowBitMask = (1 << 3);

var cloudsBitMask = (1 << 5);

var qa = image.select('pixel_qa');

var mask = qa.bitwiseAnd(cloudShadowBitMask).eq(0)

.and(qa.bitwiseAnd(cloudsBitMask).eq(0));

return image.updateMask(mask)

.copyProperties(image)

.copyProperties(image, ["system:time_start",'system:time_end']);

}

var star_date = '2020-07-01'//定义起始时间

var end_date = '2020-09-30'//定义终止时间

var L8_ = ee.ImageCollection("LANDSAT/LC08/C01/T1_SR").filterBounds(roi)

.filterDate(star_date, end_date)

//.filter(ee.Filter.lt('CLOUD_COVER',10))//云量设//按云量排序//云量过滤

.map(rmL8Cloud)

.median()

.clip(roi)

print(L8_)

// display the image

Map.addLayer(L8_,{min:0.0,max:0.4, bands:'B5,B4,B3'},'Image');

Map.centerObject(roi)

var img = L8_.select(['B2','B3', 'B4', 'B5','B6', 'B7'])

//NDVI

var ndvi = img.expression(

'((NIR - RED) / (NIR + RED))', {

'NIR': img.select('B5'),

'RED': img.select('B4')

}).rename('NDVI');

//NDTI

var ndti = img.expression(

'((SWIR1-SWIR2)/(SWIR1 + SWIR2))',{

'SWIR1': img.select('B6'),

'SWIR2': img.select('B7')

}).rename('NDTI');

//EVI

var evi = img.expression(

'2.5 * ((NIR - RED) / (NIR + 6 * RED - 7.5 * BLUE + 1))', {

'NIR': img.select('B5'),

'RED': img.select('B4'),

'BLUE': img.select('B2')

}).rename('EVI');

//NDBI

var ndbi = img.expression(

'(MIR - NIR)/(MIR + NIR)', {

'NIR': img.select('B5'),

'MIR': img.select('B6')

}).rename('NDBI');

//MNDWI

var mndwi = img.expression(

'(GREEN - MIR) / (GREEN + MIR)', {

'GREEN': img.select('B3'),

'MIR': img.select('B6')

}).rename('MNDWI')

Map.addLayer(ndvi, {'min':0,'max':1,'palette':['#A9A9A9','00FF00']}, 'NDVI')

//导出函数

Export.image.toDrive({

image: ndvi,//

description: 'ndvi2020',//

scale: 30,//

region:roi,//

maxPixels:1e13,

fileFormat: 'GeoTIFF',

});

GEE代码链接：​​https://code.earthengine.google.com/8e1e1e3fe4cbf8e6e8da8f987c8aaa33​​

版权声明：本文内容由网络用户投稿，版权归原作者所有，本站不拥有其著作权，亦不承担相应法律责任。如果您发现本站中有涉嫌抄袭或描述失实的内容，请联系我们jiasou666@gmail.com 处理，核实后本网站将在24小时内删除侵权内容。