利用
TM
计算
NDVI
问题——
NDVI
值不连续
NDVI
(
Normalized Difference Vegetation Index,
归一化植被指数,标准差异植被指数
)
表达式:
NDVI=(p(nir)-p(red))/(p(nir)+p(red))
其中
p
(
nir
)表示近红外波段,
p
(
red
)表示
可见光红波段,
它和植物的蒸腾作用、
太阳光的截取、
光合作用以及地表净初生产力
(
NPP
)
等密切相关。
1
、
NDVI
的应用:检测植被生长状态、植被覆盖度和消除部分辐射误差等;
2
、
-1<=NDVI<=1
,负值表示地面覆盖为云、水、雪等,对可见光高反射;
0
表示有岩石或
裸土等,
NIR
和
R
近似相等;正值,表示有植被覆盖,且随覆盖度增大而增大;
3
、
NDVI
的局限性表现在,用非线性拉伸的方式增强了
NIR
和
R
的反射率的对比度。对于
同一幅图象,分别求
RVI
和
NDVI
时会发现,
RVI
值增加的速度高于
NDVI
增加速度,即
NDVI
对高植被区具有较低的灵敏度;
4
、
NDVI
能反映出植物冠层的背景影响,如土壤、潮湿地面、雪、枯叶、粗超度等,且与
植被覆盖有关。
以
TM
为例:
NDVI=
(
B4-B3
)
/(B4+B3)
1.arcgis
打开
TM
,
B4
,
B3
两波段数据
2.
空间分析里面有个栅格计算器的工具(
spatial analyst-raster calculator
)
输入(
B4-B3
)
/(B4+B3)
3.
点击
evaluate
(执行)
4.
查看结果
这个时候你会发现
NDVI
的值范围是正确的(
-1,1
),但是它是不连续变化的。根据研究区
以及经验判断这个结果是错的。
分析原因,寻找解决的办法。
——
数据类型的为整形
需要修改栅格数据的类型,改为
float pointing
浮点型!
修改办法:
栅格计算器里面有
float
工具如下图:
转换
4,3
两个波段的数据
type
为
float point
再次运算。