谷歌地球引擎GEE：波段代数与波段计算

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

本文是谷歌地球引擎（Google Earth Engine，GEE）系列教学文章的第六篇，更多GEE 文章请参考专栏：# GEE学习与应用。

首先，依据# 谷歌地球引擎GEE：栅格数据的导入、筛选与可视化中内容，我们将Landsat 8 Collection 1 Tier 1 的大气表观反射率TOA Reflectance 产品导入GEE，并对其做好重命名。接下来，同样依据# 谷歌地球引擎GEE：栅格数据的导入、筛选与可视化中内容，将Landsat 8 Collection 1 Tier 1 的大气表观反射率TOA Reflectance产品按照时间进行选取（具体时间范围大家任意选择即可），并求取对应时间范围内的平均值；同时，将经过上述操作后的图层在地图中加以显示。

js 复制代码

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

接下来，我们就开始栅格计算。在# 谷歌地球引擎GEE：基于像素加以筛选、掩膜、提取分析中，我们利用normalizedDifference()函数计算了NDVI 图层；而在本文中，我们则尝试手动进行波段运算的方式计算NDVI 图层，具体代码如下。其中，.select()表示从遥感影像中选择对应的波段，.subtract()函数表示波段之间进行相减的操作；同理，.add()函数表示相加操作，.divide()表示相除操作。

js 复制代码

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：遥感影像可视化修改的代码方法中方法，对NDVI 图层加以手动可视化配置，并将配置嵌入到代码中，从而实现每一次代码运行得到的NDVI图层可视化设置都是一致的。

具体代码如下：

js 复制代码

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

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

至此，大功告成。