植被信息提取

1、NDVI指数与植被覆盖度的计算步骤收集整理资料如下NDVI:归一化植被指数和植物的蒸腾作用、太阳光的截取、光合作用以及地表净初级生产力等密切相关。1、NDVI的应用：检测植被生长状态、植被覆盖度和消除部分辐射误差等；2、-1<=NDVI<=1,负值表示地面覆盖为云、水、雪等，对可见光高反射；0表示有岩石或裸土等，NIR和R近似相等；正值，表示有植被覆盖，且随覆盖度增大而增大；3、NDVI的局限性表现在，用非线性拉伸的方式增强了NIR和R的反射率的对比度。对于同一幅图象，分别求RVI和NDVI时会发现，RVI值增加的速度高于NDVI增加速度，即NDVI对高植被区具有较低的灵敏度；4、

2、NDVI能反映出植物冠层的背景影响，如土壤、潮湿地面、学、枯叶、粗超度等，且与植被覆盖有关；erdas里面利用NDVI提取植被指数的步骤：1、在Erdas的主工具中选择Interpreter模块，出现ImageInterpreter对话框2、然后选择SpectralEnhancement,会弹出SpectralEnhancement对话3、再选择Indices选项出现Indices对话框以TM数据为例进行说明，选择InputFile,选择OutputFile,在OutputOptions的Sensor中选择LandstTM在SelectFunction里面选择NDVI,DataType默认为F

3、loat不用改变，可以发现最下面的Function显示band4-band3/band4+band3,这个就是NDVI的计算公式。最后选择OK即可完成，这里要注意的是没有OutputFile的话Ok按钮时不能使用的。如果NDVI计算的话在ENVI是最方便的在Transform菜单下就有，同时ENVI的波段计算功能也很方便完成NDVI计算。另外，ERDASMODELNDVI分类首先说如何做NDVI,虽然ERDAS1有个现成专门可以做NDVI的的地方，但是我们注意到TM4+TM3T能为0,当除数为0时系统会报错，所以应该在分母上加0.001或0.0001都可以。这样分母就不会为0了，同时注意输出图

4、象类型要是floatsingle,否则做出来的结果可能是空白图象。给NDVI图象进行再分类：注意输出图象类型为thematic。其中中间一步要这样设置：注意先选中左框里的原始NDVI波段，然后点ADDCLUMMN钮，并定义你的范围。ENVI中提取NDVI值1、打开图像2、主菜单中Transforms-NDVI,显示如下窗口选择影像-OK3、在以下窗口中输入相应的参数点击OK即可。说明：选择浮点型，ENVI数值范围保持为-11;选择字节型，键入最小NDVI值，该值将被拉伸为0;键入最大NDVI值，该值将被拉伸为255,获得的ENVI将被拉伸为0255的范围。这样得到的ndvi图式HDRR#式的，

5、要想在ERDA阱打开需要转换成IMG格式的：打开ERDA歌件，点击import,显示如下窗口：ENVI下植被覆盖度的遥感估算植被覆盖度是指植被（包括叶、茎、枝）在地面的垂直投影面积占统计区总面积的百分比。容易与植被覆盖度混淆的概念是植被盖度，植被盖度是指植被冠层或叶面在地面的垂直投影面积占植被区总面积的比例。两个概念主要区别就是分母不一样。植被覆盖度常用于植被变化、生态环境研究、水土保持、气候等方面。植被覆盖度的测量可分为地面测量和遥感估算两种方法。地面测量常用于田间尺度，遥感估算常用于区域尺度。估算模型目前已经发展了很多利用遥感测量植被覆盖度的方法，较为实用的方法是利用植被指数近似估算植被覆

6、盖度，常用的植被指数为NDVL下面是李苗苗等在像元二分模型的基础上研究的模型：VFC=(NDVI-NDVIsoil)/(NDVIveg-NDVIsoil)(1)其中，NDVIsoil为完全是裸土或无植被覆盖区域的NDVI值，NDVIveg则代表完全被植被所覆盖的像元的NDVI值，即纯植被像元的NDVI值。两个值的计算公式为：NDVIsoil=(VFCmax*NDVImin-VFCmin*NDVImax/(VFCmax-VFCmin)(2)NDVIveg=(1-VFCmin)*NDVImax-(1-VFCmax)*NDVImin)/(VFCmaxVFCmin)(3)利用这个模型计算植被覆盖度的关

7、键是计算NDVIsoil和NDVIveg这里有两种假设：1) 当区域内可以近似取VFCmax=100%VFCmin=0%公式(1)可变为：VFC=(NDVI-NDVImin)/(NDVImax-NDVImin)(4)NDVImax和NDVImin分别为区域内最大和最小的NDVI值。由于不可避免存在噪声，NDVImax和NDVImin一般取一定置信度范围内的最大值与最小值，置信度的取值主要根据图像实际情况来定。2) 当区域内不能近似取VFCmax=100%VFCmin=0%当有实测数据的情况下，取实测数据中的植被覆盖度的最大值和最小值作为VFCmafflVFCmin,这两个实测数据对应图像的ND

8、VI作为NDVImax和NDVImin当没有实测数据的情况下，取一定置信度范围内的NDVImax和NDVIminVFCmax和VFCmin根据经验估算。实现流程下面我们以“当区域内可以近似取VFCmax=100%VFCmin=0%情况下，整个影像中NDVIsoil和NDVIveg取固定值，介绍在ENVI中实现植被覆盖度的计算方法。使用的数据是经过几何校正、大气校正的TM影像。(1) 选择Transform->NDVI,利用TM影像计算NDVI(2) 选择BasicTools->Statistics->ComputeStatistics在文件选择对话框中，利用研究区地区的矢量数

9、据生成的ROI建立一个掩膜文件。选择统计文件及掩膜文件计算统计参数(3)得到研究区的统计结果。在统计结果中，最后一列表示对应NDVI值的累积概率分布。我们分别取累积概率为5%和90%的NDVI值作为NDVImin和NDVImax这里得到：NDVImax=0.522991NDVImin=0.031766OptionsE A t统计结果NDVI0S9379 046830 0402£0 03373002031 014091 007S32 OOD9S2 。口”制 012117 019tG7 025217O3G3L6酰1415 &E75S5 O64S15 0i7106 4 O77bl4

10、。旧41641007536 1091372 118478( 127371137076 imm 1S73?9 1660370201141 223829 252852284S76330654 3b944b 431721 4430540J31 b?$64 6B12S3 726378381704473415174411914606595B15115Lj：Z74757 4皿5114 4m S330 E771 538747S81105643S12S3B6S&3 93G3£ 33414 S9DG5 96178 1DEO1O9602Q965711S£29 22686 29023 3

11、2C24 4” 7g 3879270275418 二 工044之5 70b33 54孙 S13G4 755331 23011.2C97弓iriri3 ：F5uit5 NT J：，n、？0 1247 0 15950 17&00,25160.21S20.38&20 02301.91722 03042.-41672.43573 0136禾I用 ENVI主菜单->Basic Tools->Band MathVFC取值为0 用公式（4）亍栏中输入：(b1 lt 0.031766 )*0+(b1 gt 0.522991 )*1+(b1 ge 0.031766 and bl le 0.522991)* (b1-0.031766)/ (0.522991-0.031766)0.031766,0.522991 , VFC取值为 1;NDVIn:WIjrnqNumber03176石020392073bl:选择NDVI图像(5)得到一个单波段的植被覆盖度图像文件，像元值表示这个像元内的平均植被覆盖度。在Dis