站长

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


<pre class="language-js line-numbers"><code>var landsat_202006=landsat_8.filterDate('2020-06-01','2020-06-30').mean();
Map.addLayer(landsat_202006,{},"Landsat");
</code></pre>
<img src="https://img.blogweb.cn/article/cdf797a6c2fe48c4a3694dc2e0219f0c.webp" loading="lazy" alt="谷歌地球引擎GEE：波段代数与波段计算">

<pre class="language-js line-numbers"><code>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");
</code></pre>
  运行上述代码，可以看到NDVI图层已经生成。
<img src="https://img.blogweb.cn/article/b8d55276284a4e12ae44d637b6cc0b88.webp" loading="lazy" alt="谷歌地球引擎GEE：波段代数与波段计算">
  这里还有一点需要注意：在对最原始的Landsat 8 Collection 1 Tier 1大气表观反射率TOA Reflectance产品（数据类型为<code>ee.ImageCollection</code>）按照时间进行选取并求取对应时间范围内的平均值后（数据类型为<code>ee.Image</code>），遥感影像各波段才可以参与计算。也就是说，<code>.subtract()</code>等这些波段运算函数仅对于<code>ee.Image</code>数据有效；而若是<code>ee.ImageCollection</code>等数据类型需要参与波段计算，还需要映射等其它操作。

<img src="https://img.blogweb.cn/article/f91ccde139454302917e58e781369983.webp" loading="lazy" alt="谷歌地球引擎GEE：波段代数与波段计算">
  具体代码如下：
<pre class="language-js line-numbers"><code>Map.addLayer(NDVI,{min:[-0.2111],max:[0.6778],palette:["ff2a08","fff810","56ff10","1889ff","e510ff"]},"NDVI");
</code></pre>
  运行后可以看到NDVI图层已经具有默认的可视化配置了~
<img src="https://img.blogweb.cn/article/073583e5c4ec4dd49eb9ffd9e833cd64.webp" loading="lazy" alt="谷歌地球引擎GEE：波段代数与波段计算">
  至此，大功告成。

本文主要对GEE中的栅格代数与波段计算操作加以介绍。运行上述代码，可以看到NDVI图层已经生成。这里还有一点需要注意：在对最原始的Landsat8Collection1Tier1大气表观反射率TOAReflectance产品（数据类型为）按照时间进行选取并求取对应时间范围内的平均值后（数据类型为），遥感影像各波段才可以参与计算。也就是说，等这些波段运算函数仅对于数据有效；而若是等数据类型需要参与波