谷歌地球引擎GEE网格代数和波段计算NDVI计算

  本文主要对GEE中的栅格代数与波段计算操作加以介绍。

  本文是谷歌地球引擎(Google Earth Engine,GEE)系列教学文章的第六篇,更多GEE文章请参考专栏:GEE学习与应用_疯狂学习GIS的博客-CSDN博客。

  首先,依据第二篇GEE教学博客(Google Earth Engine谷歌地球引擎GEE数据导入后筛选并显示_疯狂学习GIS的博客-CSDN博客)中内容,我们将Landsat 8 Collection 1 Tier 1的大气表观反射率TOA Reflectance产品导入GEE,并对其做好重命名。接下来,同样依据第二篇GEE教学博客(Google Earth Engine谷歌地球引擎GEE数据导入后筛选并显示_疯狂学习GIS的博客-CSDN博客)中内容,将Landsat 8 Collection 1 Tier 1的大气表观反射率TOA Reflectance产品按照时间进行选取(具体时间范围大家任意选择即可),并求取对应时间范围内的平均值;同时,将经过上述操作后的图层在地图中加以显示。

var landsat_202006=landsat_8.filterDate('2020-06-01','2020-06-30').mean();
Map.addLayer(landsat_202006,{},"Landsat");

  接下来,我们就开始栅格计算。在第五篇GEE教学博客(Google Earth Engine谷歌地球引擎GEE多图层像元条件筛选与掩膜处理_疯狂学习GIS的博客-CSDN博客)中,我们利用normalizedDifference()函数计算了NDVI图层;而在本文中,我们则尝试手动进行波段运算的方式计算NDVI图层,具体代码如下。其中,.select()表示从遥感影像中选择对应的波段,.subtract()函数表示波段之间进行相减的操作;同理,.add()函数表示相加操作,.divide()表示相除操作。

var band_4=landsat_202006.select('B4');
var band_5=landsat_202006.select('B5');
var numerator=band_5.subtract(band_4);
var denominator=band_5.add(band_4);
var NDVI=numerator.divide(denominator);
Map.addLayer(NDVI,{},"NDVI");

  运行上述代码,可以看到NDVI图层已经生成。

  这里还有一点需要注意:在对最原始的Landsat 8 Collection 1 Tier 1大气表观反射率TOA Reflectance产品(数据类型为ee.ImageCollection)按照时间进行选取并求取对应时间范围内的平均值后(数据类型为ee.Image),遥感影像各波段才可以参与计算。也就是说,.subtract()等这些波段运算函数仅对于ee.Image数据有效;而若是ee.ImageCollection等数据类型需要参与波段计算,还需要映射等其它操作。

  此外,我们还可以按照第四篇GEE教学博客(Google Earth Engine谷歌地球引擎GEE栅格数据图层可视化设置代码嵌入_疯狂学习GIS的博客-CSDN博客_栅格数据可视化)中方法,对NDVI图层加以手动可视化配置,并将配置嵌入到代码中,从而实现每一次代码运行得到的NDVI图层可视化设置都是一致的。

  具体代码如下:

Map.addLayer(NDVI,{min:[-0.2111],max:[0.6778],palette:["ff2a08","fff810","56ff10","1889ff","e510ff"]},"NDVI");

  运行后可以看到NDVI图层已经具有默认的可视化配置了~

发表评论

您的电子邮箱地址不会被公开。

此站点使用Akismet来减少垃圾评论。了解我们如何处理您的评论数据