Google Earth Engine谷歌地球引擎GEE栅格代数与NDVI波段计算手动求取

  本文主要对GEE中的栅格代数与波段计算操作加以介绍。本文是谷歌地球引擎（Google Earth Engine，GEE）系列教学文章的第六篇，更多GEE文章请参考专栏：GEE学习与应用（https://blog.csdn.net/zhebushibiaoshifu/category_11081040.html）。

  首先，依据第二篇GEE教学博客（https://blog.csdn.net/zhebushibiaoshifu/article/details/117296956）中内容，我们将Landsat 8 Collection 1 Tier 1的大气表观反射率TOA Reflectance产品导入GEE，并对其做好重命名。接下来，同样依据第二篇GEE教学博客（https://blog.csdn.net/zhebushibiaoshifu/article/details/117296956）中内容，将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教学博客（https://blog.csdn.net/zhebushibiaoshifu/article/details/118909637）中，我们利用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教学博客（https://blog.csdn.net/zhebushibiaoshifu/article/details/117573181）中方法，对NDVI图层加以手动可视化配置，并将配置嵌入到代码中，从而实现每一次代码运行得到的NDVI图层可视化设置都是一致的。

  具体代码如下：

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

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

欢迎关注公众号/CSDN/知乎/微博：疯狂学习GIS