|
转载的希望对楼主有用!
本人初次学习ENVI软件,就碰到了一个大的bug,实验了一天终于搞定了!记录如下
NDVI的作用就不多说了,计算公式大家都知道是利用(近红外波段-红波段)/(近红外波段+红波段)的公式,在envi中使用transform下的NDVI公式时需要时组合波段(即多个波段组合一起的文件,如果各个波段分开的则无法使用)。用此方法计算的NDVI没问题。
但是如果你要是使用basictools下的bandmath的话就要注意了。
错误一:书写公式时如果用(b4-b3)/(b4+b3),那么你会发现结果一片黑或一片白,主要是计算值均取了整数,例如(50-10)/(50+10)=40/60=0因为整数运算,结果ndvi值为0。
错误二:查询资料发现必须指定波段的取值为float,计算公式改为float(b4-b3)/float(b4+b3),结果发现NDVI对正值均计算正常,比如植被,裸地等。但是对于水体并不算大家所谓的小于0的值,而是全部大于0,而且是大于1的值。郁闷之极,询问了很多人,解释说bandmath是针对反射率的,不能直接针对DN值计算,但是课本说用DN值后者反射率都可以啊。
几经尝试,发现下面的公式时可行的:NDVI=(float(b4)-float(b3))/(float(b4)+float(b3))可行。
以上是个人多次尝试总结,如果有不对请指定出来。继续改进。 |
|