像元二分模型

1、ENVI下植被覆盖度的遥感估算植被覆盖度是指植被(包括叶、茎、枝)在地面的垂直投影面积占统计区总面积的百分比。容易与植 被覆盖度混淆的概念是植被盖度，植被盖度是指植被冠层或叶面在地面的垂直投影面积占植被区总面积的 比例。两个概念主要区别就是分母不一样。植被覆盖度常用于植被变化、生态环境研究、水土保持、气候 等方面。植被覆盖度的测量可分为地面测量和遥感估算两种方法。地面测量常用于田间尺度，遥感估算常用于 区域尺度。估算模型目前已经发展了很多利用遥感测量植被覆盖度的方法，较为实用的方法是利用植被指数近似估算植被覆盖度，常用的植被指数为NDVI。下面是李苗苗等在像元二分模型的基础上研究的模型：VFC

2、 = (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) /( VFCmax- VFCmin) (3)利用这个模型计算植被覆盖度的关键是计算NDVIsoil和NDV

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

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

5、，利用研究区地区的矢量数据生成的ROI建立一个掩膜文件。)CompuStatistics Input Fil选择统计文件及掩膜文件 Compute Statistics ParametersJ c Stits J HiEtersjiisCovari anceSamples Resi ze Factor 1|门门门门Lin«s Remire Factor 1-00000V Output to the ScreenOutput to a Statisties FileOutput to 筑 Text Rep cart File计算统计参数(3) 得到研究区的统计结果。在统计结果中，最后一

6、列表示对应NDVI值的累积概率分布。我们分别取累积概率为5%和90%的NDVI值作为NDVImin和NDVImax。这里得到：NDVImax=0.522991NDVImin=0.031766 Statistics Results: NDVI.img E 回Select Plot*Clear PlotMin/MciH/Macin：yDW*innq .THie Options0.5>00占2X-0.5kOBand Number-0.053379186292011410.10241 1055单-0.0463302268B2238290.12471 2301-0 0402802902325285

7、20.15951.3097-0,033730320242848760.17601.5657-0.0271814577B3306540.25161 8172-0.0206313B7923694460.21322.0304-0 0140817027S4397210.38622.4167-0.00753241S54439060.02302 4397-0.0009S21044255483310.57393 013E.0 0055607B6336269640.43223.44570 012117542896812530.29943.74410 01G667313347626370.44734.19140

8、 025217755338331700.41514 6065Ll0 031766820039209730.4551n in： 0.03B3168656310075360.47575 53730 0443663333610513720.46035 99310.0514159341411E47860.51346.51150 0579650900512737910.48927.00070.06451F9697G13707690.53307 53360 071064105010147S7790.5771S.11080.0776149802015737990.5307S.64950 0841648657

9、116603700 . 47589.1253T1统计结果（4）根据公式（4）,我们可以将整个地区分为三个部分：当NDVI小于0.031766, VFC取值为0; NDVI大于0.522991, VFC取值为1;介于两者之间的像元使用公式（4）计算。利用ENVI主菜单->Basic Tools->BandMath，在公式输入栏中输入：(b1 It 0.031766)*0+(b1 gt 0.522991)*1+(b1 ge 0.031766 and bl le 0.522991)* (b1-0.031766)/(0.522991-0.031766)b1:选择NDVI图像(5) 得到一个单波段的植被覆盖度图像文件，像元值表示这个像元内的平均植被覆盖度。在Display显示(6) 选择 Tools->Color Mapping-&