很短的几句代码,可是我却花了很长的时间才写出来,因为array那里的除法运算结果老是不对,正常在-1-1之间的。从别的资料摘来处理NDVI计算的array代码处,出现了很多问题,可能它用了什么优化计算的函数,但是结果不对,果断放弃了。直接硬算,影像波段是整数,转成浮点型数字就行,然后再参与运算得出了正确的结果范围。这个效率还是不行,用c++&&GDAL计算速率快得多了。

Python&&GDAL实现NDVI的计算方式