图13 林地植被覆盖率结果
2.4 城镇用地
城镇用地使用的是等密度模型，城镇用地是林地和建筑的混合 1) NDVI最大应取林地的NDVI最大值0.990759
2) 选取纯净的城镇样本，计算NDVI 并统计，得到NDVI最小即NDVI0值为-0.168374
图14 城镇用地NDVI统计结果
3) 将数据带入上面的公式①，利用公式fg=float(b1+0.168374)/float(0.990759+0.168374)进行
波段运算得到结果，其中b1为林地计算得到的NDVI。 4) 计算后的影像进行掩膜只留下城镇用地的部分。 结果如下：
11
图15 城镇用地植被覆盖率结果
2.5 农田
农田是农田和裸土的混合，适合非等密度模型。非密度模型与等密度模型相似，该模型也假设像元中植被类型较为单一；但此时的植被垂直密度较小，即LAI<<∞。因而植被覆盖部分的NDVI值( NDVIg) 需由Bear定律确定。
·······································②
其中：
·······················③
对于LAI的估算可用如下方法：
其中，i = 1 ,2 ，分别对应于TM3 ,TM4 波段，s 指土壤，DN∞i可从遥感影像中获得；r1∞ ,r2∞据Price 研究可分别取0.05 和0.7。 土壤线方程DNs2?DNs1?a'?b'
参数a′, b′可从可见光和近红外波段的散点图中通过回归法求取。
c1, c2取0.6 和0.21 ，是Price根据Landsat TM 卫星资料推算出的，适用于垂直密度不高的农田和草地。 1) 求土壤线方程
12
图16 提取土壤线上的点
土壤线方程4000030000y = 1.158x -1739NIR2000010000005000100001500020000250003000035000Red 图17 拟合土壤线方程
2) 通过程序计算LAI
图18 LAI计算结果
3) 选择纯净的农田样本，计算其NDVI并统计最大值为0.8725，NDVI0为裸土的NDVI，在
13
2.3中已经算出，应取0.036493
4) 将求得的数据带入公式③，通过波段运算求NDVIg，公式为：
0.8725-(0.8725-0.036493)*exp(-b1)，其中b1未步骤2中求得的LAI图。
5) 根据公式②求fg=float(b1-0.036493)/float(b2-0.036493)，其中b1为农田的NDVI，b2为计
算得到的NDVIg。
6) 计算后的影像进行掩膜只留下农田的部分。 结果如下：
图19 农田植被覆盖率结果
2.6试验区相加
利用波段运算将每个试验区计算求得的植被覆盖度图相加即得到最终的植被覆盖率图。结果如下：
图20 植被覆盖率结果
14
【实验结果】
图21 植被覆盖率专题图
15