基于NDVI指数的微山湖水生植被分级研究

凌娟

（枣庄学院旅游与资源环境学院，山东 枣 庄 2 77100）

归一化植被指数(简称NDVI)是反映地表植被覆盖的重要指标，可通过对遥感影像的处理而提取，进而通过数学模型计算植被覆盖度［1］。在利用NDVI对植被覆盖进行监测评价方面，国内外已进行大量的研究与试验,是遥感估算植被覆盖度研究中最常用的植被指数［2-3］。

1 研究区概况

微山湖系中国北方最大的淡水湖，位于山东鲁南地区，南北长230km，总面积2100km2。京杭大运河纵贯南北。湖区水生植物生长茂盛，种类繁多。主要植物群落有挺水植被带和沉水植被带，水生植被盖度达89.9%，有“水下森林”的美誉。

2 研究方法

2.1 遥感数据选取

Landsat陆地卫星是全球生态观测研究最系统、数据连续性最好的遥感卫星，在地学、环境科学领域广泛应用。

2.2 遥感数据预处理

影像数据经辐射定标、去除条带、图像剪裁、大气校正、获取地表反射率等处理。几何纠正采用高斯克里格投影，基准面采用WGS-84。通过波段分析法，获取湖区543标准假彩色图像。

2.3 NDVI分级

NDVI指数被定义为近红外波段（DNNIR）与可见光红光波段（DNR）数值之差与这2个波段数值之和的比值，即 NDVI=（DNNIR－DNR）/（DNNIR＋DNR）和植物的蒸腾作用、太阳光的截取、光合作用以及地表净初级生产力等密切相关［4］。-1＜=NDVI＜=1，负值表示地面覆盖为云、水、雪等，对可见光高反射；0表示有岩石或裸土等，NIR和R近似相等；正值表示有植被覆盖，且随覆盖度增大而增大。

分析表明NDVI值最小值-0.16378，最大值0.51312，数据合理，满足密度分割的需要。以标准差分级法为基础，通过在NDVI平均值中加、减标准差产生类型间距，结合NDVI指数值定义及现场查验，将分级阈值定为：-0.16378、-0.13998、0.047752、0.097991、0.298902和0.51312，共分为疑似采砂造成的裸露沙地、可通航的开放水域、沉水植物生长区、挺水植物生长区、陆生植物生长区共5个区域。影像数据像元边长为30m×30m，每个像元代表约900m2的面积。通过像元数对整个湖区植被信息提取计算，形成微山湖NDVI指数分级图（见图1）。

图1 研究区NDVI指数分级图Fig.1 NDVI index classification map of the study area

3 结果与讨论

目前常用的多光谱遥感数据很难精确识别水华和水草。利用NDVI指数分级，可及时有效的反映水生植被分布状况。对遥感影像数据进行处理和水生植被分布信息提取，可以作为水生植被监测的信息源。

在研究区NDVI指数分级图中深绿色代表生长最旺盛的陆生植被，多为沿湖稻田，主要分布于研究区北半部；挺水植物生长范围最广；沉水植物沿湖分布，随着水深的增大而递减。灰色水体区域较春季遥感影像明显增大，植被繁茂，航道堵塞严重。湖区北部、中部发现2处小型红色区域，位于浅水区，有交通痕迹，少植被分布，疑似偷采或废弃的采沙场。

根据NDVI分级信息和水生植被分级统计，结合水文条件、人类活动，研究区植被分布特点如下：①开放水域受京杭大运河航运等人类活动影响，水生植被分布极不均匀，是NDVI指数值波动最大的区域；②沉水植被区和挺水植被区是湖区水生植被主要分布区域，像元、NDVI指数值较为平均，无明显波动；③陆生植被区有着丰富的植被资源，NDVI指数值稳定，基本无区域性变化；④沉水植被区和挺水植被区处于养殖区和护堤附近，约占湖区总面积的40%，是开展湿地旅游、渔业养殖主要区域。大量人工构建物破坏了水面的完整性，渔网、水箱构成一个个独立区域，阻碍航运，富营养化水体，可导致藻类暴发。

4 结语

利用这种分级方法处理的图像能够及时、准确的反映出水生植被区域分布和等级分区，为湖区生态环境监测治理提供依据。研究发现渔业养殖活动对湖区植被分布有着较大影响。在NDVI密度分割影像中发现疑似采砂区域，对于违法采砂活动，应采取有效措施，消除采砂活动对湖区生态环境的影响。

［1］ 蔡博峰，刘春兰，陈操操，等.露天煤矿生态环境遥感监测与评价方法研究——以霍林河～号露天矿为例［J］.中国矿业，2009，18（6）：61-64．

［2］ 李忠峰，李雪梅，蔡运龙，等.基于SPOT VEGETATION数据的榆林地区土地覆盖变化研究［J］.干早区资源与环境，2007,21（2）:56-59.

［3］ 霍艾迪，张广军，武苏里，等.利用MODIS-NDVI进行沙化土地评价研究——以 陕西省北部地区为例［J］.干早地区农业研究，2008，26（2）：154-158．

［4］ 赵英时.遥感应用分析原理与方法［M］.北京：科学出版社，2003:374-375.