李京忠１，２，刘咏梅１，冯非凡３，高 原１，常 伟１，韩亚磊２
（１．西北大学城市与环境学院，西安 ７１０１２７； ２．许昌学院城乡规划与园林学院，河南 许昌 ４６１０００；
３．中国地质大学地球物理与空间信息学院，武汉 ４３００７４）
摘要：利用１９９０年和２０００年两期河南省ＴＭ遥感影像数据，通过对数据进行大气校正、辐射校正、图像镶嵌、图像分割等数据预处理工作，然后计算出每期影像的ＮＤＶＩ（归一化植被指数），利用像元二分模型把植被指数转换为植被覆盖度，并对这两个时期的植被覆盖度进行对比分析。结果表明，从植被覆盖度的时间变化上看，１０年间极低与低植被覆盖度类型的面积减少，但是中、中高和高植被覆盖度类型的面积明显增加，这说明河南省的植被覆盖总体上是好转的；从植被覆盖度的空间变化上看，豫北濮阳、鹤壁，豫东商丘，豫南驻马店等区域植被覆盖度增长明显，而豫西平顶山、三门峡及商丘南部植被覆盖度减少变化明显。
教育期刊网 http://www.jyqkw.com
关键词 ：ＴＭ影像；河南省；植被覆盖度；ＮＤＶＩ
中图分类号：Ｘ８２６ 文献标识码：A 文章编号：0439－８114（２０15）16－3908-03
DOI:10.14088/j.cnki.issn0439-8114.2015.16.019
收稿日期：２０１４－０９－１６
基金项目：国家自然科学基金项目（４１１７１２２５）；２０１３年河南省科技计划项目（１３２４００４１１０２１）
作者简介：李京忠（１９７８－），男，山东临沂人，讲师，博士，主要从事遥感信息提取研究，（电话）１８６３９７２２３８０（电子信箱）ｚｈｏｎｇ＿ｌｉｊ＠１６３．ｃｏｍ；
通信作者，刘咏梅（１９７０－），女，陕西西安人，副教授，博士，主要从事区域水土流失遥感监测，（电子信箱）ｌｉｕｙｍ＠ｎｗｕ．ｅｄｕ．ｃｎ。
植被覆盖度是指示地表植物生长态势和变化的重要参数，可以通过该参数研究特定区域内的植被变化情况，对于了解该区域的生态环境现状和发展趋势具有重要的现实意义［１－３］。植被覆盖度的传统量算方法主要就是地面测量，限制较大而且精度不高，难以实现大面积的信息提取［４］。基于遥感影像的提取方法具有同步获取大区域范围内的数据信息并且受主观影响因素较小、精度较高等特点，成为当今地面植被覆盖度信息提取和植被监测的主要手段［５］。
河南省是中国重要的农业粮食基地、农副产品加工基地，同时也是国家东西部地区的连接枢纽和中原经济区的主要承载地。因此河南省的粮食产量和农副产品的加工量对于国家的粮食安全、经济发展有着不可忽视的重要作用［６］。而植被覆盖度的优劣会引起生态环境的变化，这对粮食的产量和农副产品的加工量有着决定性的影响［７］。本研究以多时相的遥感影像为数据源，利用像元二分模型估算区域植被覆盖度，对比不同时期的植被覆盖状况从而得到研究区植被覆盖度的空间和时间变化特征，并分析其与相应区域的经济发展、城市化进程的关系，为河南省经济发展过程中的生态环境评价提供数据支撑。
１ 研究区概况
河南省位于中国中部，北纬３１°２３′－３６°２２′，东经１１０°２１′－１１６°３９′之间。国土面积１６．７万ｋｍ２，居全国省（区、市）第十七位，约占全国总面积的１．７３％。河南省地处沿海与中西部结合部，地势西高东低，北、西、南三面依次由太行山脉、伏牛山脉、桐柏山脉和大别山脉沿省界呈半环形分布；中、东部属华北平原南部区域；西南部为南阳盆地，跨越黄河、淮河、海河、长江四大水系，山水相连。在全省面积中，山地丘陵面积７．４万ｋｍ２，占全省总面积的４４．３％；平原和盆地面积９．３万ｋｍ２，占总面积的５５．７％。复杂多样的土地类型为农、林、牧、渔业的综合发展提供了有利的条件。
２ 数据与方法
２．１ 数据源与预处理
数据主要来自国际科学数据服务平台、美国地质调查局下属的地球资源观测和科学中心（ＥＲＯＳ）、马里兰大学地球数据中心的Ｌａｎｄｓａｔ ＴＭ影像数据及河南省的行政区域矢量图。
本研究中收集河南省在１９９０年和２０００年共２８景的ＴＭ数据，要求该数据尽量在植被覆盖较好的夏季７、８月份，并对遥感数据进行辐射校正、地形校正、几何精校正等预处理工作。然后对所有的影像数据进行匀色镶嵌处理，并根据河南省的行政区划图进行裁剪得到研究区遥感影像集。
２．２ 植被覆盖度的提取
植被指数是利用红光和红外波段的不同组合进行植被信息指示的参数，可以定量地表明植被的活力，相对于单个波段来说该参数对植被生物量的探测具有更好的灵敏性［８］。目前，全球关于植被指数的模型超过１５０种，但是这些植被指数模型中只有极少数经过系统验证并应用到实践中［９］。其中归一化植被指数（ＮＤＶＩ）由于检测灵敏度高、检测范围较宽，可消除地形和群落结构的阴影和辐射干扰，削弱太阳高度角和大气所带来的噪音等优点，取得了广泛的应用。ＮＤＶＩ对植被状况的评估非常有效，已被广泛应用到地表植被生产力评估、植被覆盖估算及变化趋势分析的研究中［１０，１１］。本研究采用常见的归一化植被指数，其计算公式为：
ＮＤＶＩ＝（ＮＩＲ－Ｒ）/（ＮＩＲ＋Ｒ） （１）
其中ＮＩＲ为近红外波段，Ｒ为红波段［１２］。
本研究以ＮＤＶＩ作为主要参数应用于像元二分模型计算河南省植被覆盖度，其基本原理是一个像元的ＮＤＶＩ值可以表达为由绿色植被部分所贡献的信息与由无植被覆盖（裸土）部分所贡献的信息两部分组成，像元二分模型植被覆盖度（ＶＦＣ）的计算公式为：
ＶＦＣ＝（ＮＤＶＩ－ＮＤＶＩｍｉｎ）÷（ＮＤＶＩｍａｘ－ＮＤＶＩｍｉｎ） （２）
ＮＤＶＩｍａｘ和ＮＤＶＩｍｉｎ分别为区域内最大和最小的ＮＤＶＩ值。由于不可避免存在噪声，ＮＤＶＩｍａｘ和ＮＤＶＩｍｉｎ一般取一定置信度范围内的最大值与最小值，也就是植被纯像元和土壤纯像元。由于植被分布、近邻地物辐射等因素影响，置信度的取值主要根据研究区影像的实际情况来定，在没有大量实测数据作为参考的情况下，通常采用累积百分数５％和９５％作为置信区间［１３］。本研究中采用累积百分数５％作为裸土的纯像元，累积百分数９５％作为植被纯像元，根据式（２）计算出１９９０年和２０００年的植被覆盖度图像。
２．３ 植被覆盖度分级
在植被覆盖度结果图的基础上，参考马志勇等［１４］对植被覆盖度等级分类范围的分类标准，将植被覆盖度图像分为５级（表１）。按照分级标准，进行密度分割，相同植被级别赋予同一颜色，利用专家系统决策树方法生成植被覆盖度分级图。
３ 结果与分析
３．１ 植被覆盖度的时间变化
依据河南省１９９０、２０００年的植被覆盖度等级图，可以计算１９９０－２０００年不同植被覆盖等级的面积及１０年间变化。河南省１０年间植被覆盖面积有所增加，由１９９０年的１５７ ４８３．９５ ｋｍ２增加到２０００年的１５７ ６２７．３４ ｋｍ２，１０年间的植被覆盖的增加面积为１４３．３９ ｋｍ２，其年增长率为１４．3 ｋｍ２／年。并统计各个植被覆盖等级的面积及１０年间面积变化情况见表２。
对表２进行分析发现，１９９０年低覆盖度类型所占面积最大，为３７ ８３４．７４ ｋｍ２；其次是高覆盖度类型，为３７ ４３３．２９ ｋｍ２；第三是中覆盖度的类型，面积为３３ ３１４．４４ ｋｍ２；极低覆盖度类型的最少，面积为２３ ３６７．６６ ｋｍ２。２０００年时植被覆盖度面积最大的是高覆盖度类型，面积为３９ ０１８．４７ ｋｍ２；其次是中覆盖度类型，面积为３６ ０１０．３８ ｋｍ２；排在第三位的是低覆盖度类型，面积为３２ ５７１．２７ ｋｍ２；面积最少的是极低覆盖度类型，面积仅为２０ ６１４．２７ ｋｍ２。从植被覆盖分等级的变化角度上分析，１０年间极低与低植被覆盖度类型的面积减少，这和城市发展的扩展和人类活动的影响密不可分，造成土地利用属性的改变，使得极低与低覆盖度的面积减少，其中低覆盖度类型减少得最多，为５ ２６３．４７ ｋｍ２。但是中、中高和高植被覆盖度类型的面积都明显增加，这说明河南省的植被覆盖总体上是好转的。特别是森林植被的覆盖上，中、高植被覆盖度的增加量明显，其中中高覆盖度类型增加的面积最多，为３ ８７９．１３ ｋｍ２。
３．２ 植被覆盖度的空间变化
由图１可以看出，河南省１９９０年植被覆盖度大于４０％的区域主要分布在中部的漯河、周口、许昌地区，南部的驻马店、大别山信阳境内及西部的伏牛山、熊耳山、崤山等地区；２０００年植被覆盖度大于４０％的区域主要分布在北部濮阳、新乡、鹤壁地区，东部的商丘、开封及周口部分区域，南部的大别山信阳境内及西部的三门峡、洛阳、南阳部分地区。
为了分析植被覆盖度的空间变化情况，在ＥＮＶＩ（完整的遥感图像处理平台）中，对研究区域内两期植被覆盖图像进行变化检测分析，以１９９０年的植被覆盖度图像作为变化前的图像，２０００年的植被覆盖度图像作为变化后的图像，通过变化检测，得到河南省１９９０－２０００年植被覆盖动态变化图（图２）。从图２可以看出豫北濮阳、鹤壁，豫东商丘，豫南驻马店等区域植被覆盖度增长明显，而豫西平顶山、三门峡及商丘南部植被覆盖度减少明显。
４ 小结与讨论
利用遥感影像和ＧＩＳ相结合的手段，可以快速、准确地反映研究区域内地表植被覆盖的动态变化。将提取的植被覆盖度分级图以及动态变化图与研究河南省行政区划图以及河南省经济社会发展变化的特征相结合，可以看出河南省在１９９０－２０００年间植被覆盖面积有所增加，植被覆盖变化的主要地物类型包括林地、草地和人工园林等，在植被覆盖减少区域主要是由于地方经济建设用地的占用、农业结构的调整、水土流失、森林滥砍滥伐等因素。随着城市化进程的日益发展，土地利用类型发生了很大变化，最重要的原因是农业用地向非农业用地的转变。同时，植被覆盖度的空间分布特点也与经济发展有着密切联系，其分布特点受到经济建设与城市化发展的影响。河南省植被覆盖度变化的空间特征明显，即植被覆盖增加的区域多是加大了生态环境的保护力度，并增加了植树造林的面积，使区域保持生态平衡；而植被覆盖减少的区域，一方面是由于城市化发展过程不可避免的对林地、草地等植被覆盖区域的占用，以进行城市建筑、道路、公共设施的建设，另一方面则是人们不重视保护环境，大量砍伐森林，加之对山区的开发，使具有高植被覆盖度的林区减少，部分区域甚至消失。
教育期刊网 http://www.jyqkw.com
参考文献：
［１］ 王冬梅，孟兴民，邢 钊，等．基于ＲＳ的武都区植被覆盖度动态变化及其驱动力分析［Ｊ］．干旱区资源与环境，２０１２，２６（１１）：９２－９７．
［２］ ＬＥＰＲＩＥＵＲ Ｃ， ＫＥＲＲ Ｙ Ｈ， ＭＡＳＴＯＲＣＨＩＯ Ｓ， ｅｔ ａｌ． Ｍｏｎｉｔｏｒｉｎｇ ｖｅｇｅｔａｔｉｏｎ ｃｏｖｅｒ ａｃｒｏｓｓ ｓｅｍｉ－ａｒｉｄ ｒｅｇｉｏｎｓ ｃｏｍｐａｒｉｓｏｎ ｏｆ ｒｅｍｏｔｅ ｏｂｓｅｒｖａｔｉｏｎｓ ｆｒｏｍ ｖａｒｉｏｕｓ ｓｃａｌｅｓ［Ｊ］． Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｊｏｕｒｎａｌ ｏｆ Ｒｅｍｏｔｅ Ｓｅｎｓｉｎｇ， ２０００， ２１（２）：２８１－３００．
［３］ ＰＵＲＥＶＤＯＲＪ Ｔ， ＴＡＴＥＩＳＨＩ Ｒ， ＩＳＨＩＹＡＭＡ Ｔ． Ｒｅｌａｔｉｏｎｓｈｉｐｓ ｂｅｔｗｅｅｎ ｐｅｒｃｅｎｔ ｖｅｇｅｔａｔｉｏｎ ｃｏｖｅｒ ａｎｄ ｖｅｇｅｔａｔｉｏｎ ｉｎｄｉｃｅｓ［Ｊ］． Ｉｎｔｅｒｎａｔｉｏｎａｌ Ｊｏｕｒｎａｌ ｏｆ Ｒｅｍｏｔｅ Ｓｅｎｓｉｎｇ，１９９８，１９（１８）：３５１９－３５３５．
［４］ 李苗苗，吴炳方，颜长珍，等．密云水库上游植被覆盖度的遥感估算［Ｊ］．资源科学，２００４，２６（４）：１５７－１６４．
［５］ 王晓江，胡尔查，李爱平，等．基于ＭＯＤＩＳ ＮＤＶＩ的内蒙古大青山自然保护区植被覆盖度的动态变化特征［Ｊ］．干旱区资源与环境，２０１４，２８（８）：６１－６５．
［６］ 李 峰．我国中部农业循环经济发展战略研究［Ｄ］．武汉：武汉大学，２０１３．
［７］ 杨 强．基于遥感的榆林地区生态脆弱性研究［Ｄ］．南京：南京大学，２０１２．
［８］ 田庆久，闵祥军．植被指数研究进展［Ｊ］．地球科学进展，１９９８，１３（４）：３２７－３３３
［９］ 邓书斌，陈秋锦．植被光谱特征与植被指数综述［Ａ］．杭州师范大学遥感与地球科学研究院．第十七届中国遥感大会摘要集［Ｃ］．杭州：杭州师范大学，２０１０．
［１０］ ＰＡＲＵＥＬＯ Ｊ Ｍ， ＥＰＳＴＥＩＮ Ｈ Ｅ， ＬＡＵＥＮＲＯＴＨ Ｗ Ｋ， ｅｔ ａｌ． Ａ ＮＰＰ ｅｓｔｉｍａｔｅｓ ｆｒｏｍ ＮＤＶＩ ｆｏｒ ｔｈｅ ｃｅｎｔｒａｌ ｇｒａｓｓｌａｎｄ ｒｅｇｉｏｎ ｏｆ ｔｈｅ Ｕｎｉｔｅｄ Ｓｔａｔｅｓ［Ｊ］． Ｅｃｏｌｏｇｙ， １９９７， ７８（３）：９５３－９５８．
［１１］ 左丽君，张增祥，董婷婷，等．ＭＯＤＩＳ／ＮＤＶＩ和ＭＯＤＩＳ／ＥＶＩ在耕地信息提取中的应用及对比分析［Ｊ］．农业工程学报，２００８，２４（３）：１６７－１７２．
［１２］ ＬＹＯＮ Ｊ Ｇ， ＹＵＡＮ Ｄ， ＬＵＮＥＴＴＡ Ｒ Ｓ， ｅｔ ａｌ． Ａ ｃｈａｎｇｅ ｄｅｔｅｃｔｉｏｎ ｅｘｐｅｒｉｍｅｎｔ ｕｓｉｎｇ ｖｅｇｅｔａｔｉｏｎ ｉｎｄｉｃｅｓ［Ｊ］． Ｐｈｏｔｏｇｒａｍｍｅｔｒｉｃ Ｅｎｇｉｎｅｅｒｉｎｇ ａｎｄ Ｒｅｍｏｔｅ Ｓｅｎｓｉｎｇ， １９９８， ６４（２）：１４３－１５０．
［１３］ 刘 静，银 山，张国盛，等．毛乌素沙地１７年间植被覆盖度变化的遥感监测［Ｊ］．干旱区资源与环境，２００９，２３（７）：１６２－１６７．
［１４］ 马志勇，沈 涛，张军海，等．基于植被覆盖度的植被变化分析［Ｊ］．测绘通报，２００７（３）：４５－４８．