一次解决Landsat 、哨兵Sentinel 21个光谱指数

     大部分光谱指数如NDVI植被指数计算起来是比较简单的，用ArcGIS或者其他GIS软件中栅格计算器即可搞定。但是有一些指数计算起来并没那么轻松，如土壤调节植被指数SAVI。

      推荐一个开源软件Orfeo Tool Box (OTB)，可以一次性计算21个光谱指数。下载地址：https://www.orfeo-toolbox.org/download/。 盘地址在文末。

21个光谱指数如下:

植被指数：NDVI、TNDVI、RVI、SAVI、TSAVI、MSAVI、MSAVI2、GEMI、IPVI等。

水体指数：NDWI、NDWI2、MNDWI、NDTI

土壤指数：RI、CI、BI、BI2

建筑指数：ISU

下载后解压，打开批处理bat文件

    打开已经合成的波段。

     合成波段可在ArcToolbox中数据管理工具——>栅格——>栅格处理——>波段合成完成，注意选择Landsat 11个波段的顺序，哨兵Sentinel同样。

合成后加入，我们可以在RGB中选择4、3、2验证一下，没有问题。

打开view——>OTB-Application browser,选择辐射指数。你会发现这里还有很多影像处理工具。

     打开后，在输入img处输入合成的影像，输出位置设置好，在下面波段选择处一一对应好各个波段。

    注意中红外MIR  对于Landsat可以不选，对应需要该波段参数的指数计算可以忽略。

然后shift所有光谱指数，执行，耐心等待。

      计算好后，你会发现“结果tif”里有21个光谱指数结果，我们在QGIS（或者ArcMap）中打开“结果tif”中波段1（代表NDVI）。

       再利用QGIS中Raster——>Raster Calculator计算一下源Landsat中4、5波段计算的NDVI比较一下：

你会发现两者NDVI数值一样。当然你还可以利用这21个波段继续组合。

软件下载地址： 

链接：https://pan.baidu.com/s/1m7co1CptqhPHnDJyZ973xg

提取码：1111