植被信息提取

1、NDVI指數(shù)與植被覆蓋度的計(jì)算步驟收集整理資料如下NDVI:歸一化植被指數(shù)和植物的蒸騰作用、太陽(yáng)光的截取、光合作用以及地表凈初級(jí)生產(chǎn)力等密切相關(guān)。1、NDVI的應(yīng)用：檢測(cè)植被生長(zhǎng)狀態(tài)、植被覆蓋度和消除部分輻射誤差等；2、-1<=NDVI<=1,負(fù)值表示地面覆蓋為云、水、雪等，對(duì)可見(jiàn)光高反射；0表示有巖石或裸土等，NIR和R近似相等；正值，表示有植被覆蓋，且隨覆蓋度增大而增大；3、NDVI的局限性表現(xiàn)在，用非線性拉伸的方式增強(qiáng)了NIR和R的反射率的對(duì)比度。對(duì)于同一幅圖象，分別求RVI和NDVI時(shí)會(huì)發(fā)現(xiàn)，RVI值增加的速度高于NDVI增加速度，即NDVI對(duì)高植被區(qū)具有較低的靈敏度；4、

2、NDVI能反映出植物冠層的背景影響，如土壤、潮濕地面、學(xué)、枯葉、粗超度等，且與植被覆蓋有關(guān)；erdas里面利用NDVI提取植被指數(shù)的步驟：1、在Erdas的主工具中選擇Interpreter模塊，出現(xiàn)ImageInterpreter對(duì)話框2、然后選擇SpectralEnhancement,會(huì)彈出SpectralEnhancement對(duì)話3、再選擇Indices選項(xiàng)出現(xiàn)Indices對(duì)話框以TM數(shù)據(jù)為例進(jìn)行說(shuō)明，選擇InputFile,選擇OutputFile,在OutputOptions的Sensor中選擇LandstTM在SelectFunction里面選擇NDVI,DataType默認(rèn)為F

3、loat不用改變，可以發(fā)現(xiàn)最下面的Function顯示band4-band3/band4+band3,這個(gè)就是NDVI的計(jì)算公式。最后選擇OK即可完成，這里要注意的是沒(méi)有OutputFile的話Ok按鈕時(shí)不能使用的。如果NDVI計(jì)算的話在ENVI是最方便的在Transform菜單下就有，同時(shí)ENVI的波段計(jì)算功能也很方便完成NDVI計(jì)算。另外，ERDASMODELNDVI分類首先說(shuō)如何做NDVI,雖然ERDAS1有個(gè)現(xiàn)成專門可以做NDVI的的地方，但是我們注意到TM4+TM3T能為0,當(dāng)除數(shù)為0時(shí)系統(tǒng)會(huì)報(bào)錯(cuò)，所以應(yīng)該在分母上加0.001或0.0001都可以。這樣分母就不會(huì)為0了，同時(shí)注意輸出圖

4、象類型要是floatsingle,否則做出來(lái)的結(jié)果可能是空白圖象。給NDVI圖象進(jìn)行再分類：注意輸出圖象類型為thematic。其中中間一步要這樣設(shè)置：注意先選中左框里的原始NDVI波段，然后點(diǎn)ADDCLUMMN鈕，并定義你的范圍。ENVI中提取NDVI值1、打開(kāi)圖像2、主菜單中Transforms-NDVI,顯示如下窗口選擇影像-OK3、在以下窗口中輸入相應(yīng)的參數(shù)點(diǎn)擊OK即可。說(shuō)明：選擇浮點(diǎn)型，ENVI數(shù)值范圍保持為-11;選擇字節(jié)型，鍵入最小NDVI值，該值將被拉伸為0;鍵入最大NDVI值，該值將被拉伸為255,獲得的ENVI將被拉伸為0255的范圍。這樣得到的ndvi圖式HDRR#式的，

5、要想在ERDA阱打開(kāi)需要轉(zhuǎn)換成IMG格式的：打開(kāi)ERDA歌件，點(diǎn)擊import,顯示如下窗口：ENVI下植被覆蓋度的遙感估算植被覆蓋度是指植被（包括葉、莖、枝）在地面的垂直投影面積占統(tǒng)計(jì)區(qū)總面積的百分比。容易與植被覆蓋度混淆的概念是植被蓋度，植被蓋度是指植被冠層或葉面在地面的垂直投影面積占植被區(qū)總面積的比例。兩個(gè)概念主要區(qū)別就是分母不一樣。植被覆蓋度常用于植被變化、生態(tài)環(huán)境研究、水土保持、氣候等方面。植被覆蓋度的測(cè)量可分為地面測(cè)量和遙感估算兩種方法。地面測(cè)量常用于田間尺度，遙感估算常用于區(qū)域尺度。估算模型目前已經(jīng)發(fā)展了很多利用遙感測(cè)量植被覆蓋度的方法，較為實(shí)用的方法是利用植被指數(shù)近似估算植被覆

6、蓋度，常用的植被指數(shù)為NDVL下面是李苗苗等在像元二分模型的基礎(chǔ)上研究的模型：VFC=(NDVI-NDVIsoil)/(NDVIveg-NDVIsoil)(1)其中，NDVIsoil為完全是裸土或無(wú)植被覆蓋區(qū)域的NDVI值，NDVIveg則代表完全被植被所覆蓋的像元的NDVI值，即純植被像元的NDVI值。兩個(gè)值的計(jì)算公式為：NDVIsoil=(VFCmax*NDVImin-VFCmin*NDVImax/(VFCmax-VFCmin)(2)NDVIveg=(1-VFCmin)*NDVImax-(1-VFCmax)*NDVImin)/(VFCmaxVFCmin)(3)利用這個(gè)模型計(jì)算植被覆蓋度的關(guān)

7、鍵是計(jì)算NDVIsoil和NDVIveg這里有兩種假設(shè)：1) 當(dāng)區(qū)域內(nèi)可以近似取VFCmax=100%VFCmin=0%公式(1)可變?yōu)椋篤FC=(NDVI-NDVImin)/(NDVImax-NDVImin)(4)NDVImax和NDVImin分別為區(qū)域內(nèi)最大和最小的NDVI值。由于不可避免存在噪聲，NDVImax和NDVImin一般取一定置信度范圍內(nèi)的最大值與最小值，置信度的取值主要根據(jù)圖像實(shí)際情況來(lái)定。2) 當(dāng)區(qū)域內(nèi)不能近似取VFCmax=100%VFCmin=0%當(dāng)有實(shí)測(cè)數(shù)據(jù)的情況下，取實(shí)測(cè)數(shù)據(jù)中的植被覆蓋度的最大值和最小值作為VFCmafflVFCmin,這兩個(gè)實(shí)測(cè)數(shù)據(jù)對(duì)應(yīng)圖像的ND

8、VI作為NDVImax和NDVImin當(dāng)沒(méi)有實(shí)測(cè)數(shù)據(jù)的情況下，取一定置信度范圍內(nèi)的NDVImax和NDVIminVFCmax和VFCmin根據(jù)經(jīng)驗(yàn)估算。實(shí)現(xiàn)流程下面我們以“當(dāng)區(qū)域內(nèi)可以近似取VFCmax=100%VFCmin=0%情況下，整個(gè)影像中NDVIsoil和NDVIveg取固定值，介紹在ENVI中實(shí)現(xiàn)植被覆蓋度的計(jì)算方法。使用的數(shù)據(jù)是經(jīng)過(guò)幾何校正、大氣校正的TM影像。(1) 選擇Transform->NDVI,利用TM影像計(jì)算NDVI(2) 選擇BasicTools->Statistics->ComputeStatistics在文件選擇對(duì)話框中，利用研究區(qū)地區(qū)的矢量數(shù)

9、據(jù)生成的ROI建立一個(gè)掩膜文件。選擇統(tǒng)計(jì)文件及掩膜文件計(jì)算統(tǒng)計(jì)參數(shù)(3)得到研究區(qū)的統(tǒng)計(jì)結(jié)果。在統(tǒng)計(jì)結(jié)果中，最后一列表示對(duì)應(yīng)NDVI值的累積概率分布。我們分別取累積概率為5%和90%的NDVI值作為NDVImin和NDVImax這里得到：NDVImax=0.522991NDVImin=0.031766OptionsE A t統(tǒng)計(jì)結(jié)果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)得到一個(gè)單波段的植被覆蓋度圖像文件，像元值表示這個(gè)像元內(nèi)的平均植被覆蓋度。在Dis