envi统计ndvi变化的区域面积并制图输入

喝醉酒的猫儿

怎样用envi计算不同时间ndvi，给出变化的图像统计变化的区域面积，并对变化的区域进行制图输入

jimi21

请论坛各位高手回答下楼主的提问，会奖励相应的铜板和威望的哦！

地质少侠

帮顶下了

地质少侠

转载的希望对楼主有用！

本人初次学习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))可行。

以上是个人多次尝试总结，如果有不对请指定出来。继续改进。

宁静的大西北

两幅ndvi图用波段运算相减生成图1。然后build mask进行掩膜图2（对图1像元为0，也就是没变化的区域）。利用掩膜图2 apply mask到图1上，统计区域就行了