基于遥感数据的流域土壤侵蚀强度快速估测方法_谭炳香

收稿日期:2004-07-12;修订日期:2005-01-18

基金项目:北京市自然科学基金重点项目“官厅库区土壤水分承载量及植被恢复机理研究”(6011003)。

作者简介:谭炳香(1966-),女,副研究员,研究方向为遥感在林业中的应用,现主要从事遥感和GIS 技术的森林资源监测与评估,森林参

数遥感估测等方面的研究。

基于遥感数据的流域土壤侵蚀强度快速估测方法

谭炳香1,李增元1,王彦辉2,于澎涛2,柳立兵3

(1.中国林业科学研究院资源信息研究所,北京　100091;2.中国林业科学研究院森林生态环境与保护

研究所,北京　100091;3.北京市延庆县水利局水土保持工作站,北京　102100)

摘要:以北京延庆县境内的妫水河流域为例,提出了一种基于遥感数据的土壤侵蚀强度快速估测方法。首先,利用遥感数据和植被指数模型提取流域内土地利用类型信息和植被覆盖度信息;其次,利用数字高程模型数据生成坡度图;然后,结合土壤侵蚀强度分级指标,将坡度图与土地类型图、植被覆盖度图空间叠加,判断和计算侵蚀强度等级,结果获得了流域土壤侵蚀强度等级图;最后,计算了流域的年平均侵蚀模数。结果表明,妫水河流域的土壤侵蚀以微度和轻度为主,所占面积比例为74.88%,极度和剧烈侵蚀很少,不到总面积的2%。整个流域的年侵蚀模数估计为1746.1t /km 2·a 。关　键　词:遥感;植被覆盖;土壤侵蚀;流域

中图分类号:TP79 文献标识码:A 文章编号:1004-0323(2005)02-0215-06

1　引　言

植被具有截持降雨、减缓径流、保土固土等生态功能,尤其是对控制土壤流失起着决定性的作用,植被盖度(特别是地面覆盖度)的大小直接影响着土壤流失的强弱。已有的大量对比观测和试验发现,植被防治土壤侵蚀的能力随植被覆盖度的增加而增强。在其它条件一定时,侵蚀量与植被覆盖度成反比关系〔1〕。因此,进行植被覆盖调查和监测,是预测水土流失和制定水土保持规划的最基础的工作。但是,常规的植被覆盖地面调查费时、费力、成本高,而且不利于大范围的植被监测调查。近些年来,遥感技术的快速发展和广泛应用为植被覆盖信息的提取提供了可靠的技术途径。相对于传统的地面调查,遥感方法具有耗资少、周期短、覆盖广等特点,遥感信息已成为无实测资料的广大地区植被覆盖估算的一种极其重要的信息源〔2～5〕。目前,有关植被覆盖的遥感研究主要集中在植被类型分类制图方面,而对植被覆盖度的应用研究相对不足。

应用遥感技术能够及时、快速地获取流域内的大范围区域植被覆盖和土壤侵蚀强度的空间分布状

况,为流域水土保持规划提供科学依据,加强水源涵养型植被建设,改善水源区生态环境,增强植被的水源涵养与土壤保护功能。因此,本文以北京市延庆县

的妫水河流域为例,试图提出一种基于遥感数据的土壤侵蚀强度快速估测方法,包括流域植被覆盖遥感监测、流域土壤侵蚀强度图绘制、流域年平均侵蚀模数计算等环节。

2　试验区及其数据获取

2.1　试验区概况

官厅水库是北京市仅次于密云水库的第二大水源,是首都水资源保护和水环境治理的重点地区。妫水河是官厅水库的一级流域和主要水源之一,主要分布在北京市延庆县境内,代表着华北半干旱区的水源地,具有相对丰富的前期调查研究基础,开展华北半干旱区水源涵养型植被研究较为理想。妫水河流域是官厅水库主要水源之一,位于北京延庆县境内。延庆县属延怀盆地的一部分,北、东、南三面环山,一般海拔600～800m ,西面是官厅水库。平原面积5.4万hm 2,占总面积的27.2%,山区面积14.5万hm 2,占总面积的72.8%,是华北平原到坝

第20卷　第2期2005年4月

遥　感　技　术　与　应　用

REM OT E S EN SING TEC HNOLOGY AND APPLICA TION

Vol .20　No .2Apr .2005

上高原的过渡地带。县内山脉大致走向北东与东西向,由中部北起佛爷顶,经九里梁至凤坨梁,形成一自然分水岭,分水岭以西为山前平原区,以东为山后区。平原区的高程均在510m 以下,由东北至西南倾斜,妫水河从东向西流入官厅水库,总流域面积1073.64km 2

,地形上构成一个典型的开启式盆地。除平原谷地耕作区外,延庆县大部分地区植被属落叶阔叶林向针阔混交林过渡的植被带,除落叶阔叶林外,还有暖性针叶林,次生的落叶灌丛和草本群落,其中以次生的落叶灌丛占最大面积。其分布特点为:盆地边缘有喜温耐旱的灌草植物,南部有少量喜温的果、柿树种。浅山丘陵大部分阳坡植被稀疏,覆盖度小,35°以上陡坡岩石裸露,少灌无草,水土流失严重,沟壑纵横,人为破坏严重,植被残破,为稀疏植被。

2.2　遥感数据获取

获取遥感图像时,时相的选择至关重要,因为最佳时相的选择是提高侵蚀调查与监测质量的基础和保障〔5〕。因此,我们选用了2001年5月19日的Landsa t5TM 多波段遥感图像(空间分辨率30m )和全色波段图像(15m)。5月中下旬的北京延庆地区乔灌木树种已展叶并开始旺盛生长,而农作物尚未覆盖地表,坡耕地在此季节较易从图像上判读出来,故土地类型的分类效果会最佳。

以最新版的1∶5万的地形图为基准,采用人机交互方式选择地面控制点,利用二次多项式和最近邻内插法对图像进行几何精度纠正,经重新选点检验,将误差控制在一个像元内。利用妫水河流域的边界线对遥感影像进行掩膜,得妫水河流域的TM 图像。将TM 多波段图像与全色波段图像融合,获得15m 空间分辨率的多波段图像。利用计算机自动分类和目视判读相结合的方法从遥感图像中提取土地利用类型信息,利用植被指数模型求算植被覆盖度。2.3　流域数字高程模型

从延庆县水利局获得了数字化的妫水河流域1∶1万等高线数据,利用遥感图像处理软件中的地形分析工具生成了妫水河流域的数字高程模型(DEM ),空间分辨率采样为与TM 图像一样大小的15m;再由DEM 数据生成坡度图。

3　技术路线与方法

3.1　技术路线

以辐射纠正和几何纠正的妫水河流域遥感图像

为信息源,提取土地利用类型信息和植被覆盖度信

息;利用遥感影像处理软件由1∶1万数字高程模型(DEM )数据生成坡度图。最后,结合土壤侵蚀强度分级指标,将坡度等级图与土地类型图和植被覆盖度等级图分别进行空间叠加,最后生成土壤侵蚀强度等级图,统计各侵蚀强度等级的面积,并计算年平均侵蚀模数(如图1所示)。

3.2　植被覆盖度的估算

植被覆盖度是指植被冠层的垂直投影面积与土地面积之比,它是描述森林植被覆盖特征的一个重要基础数据,广泛应用于森林水文、森林生态等方面的研究,也是许多生态模型中的基本参数。准确、快速地获取地表植被覆盖度及其变化信息,对揭示地表特征时空变化规律、探讨变化的驱动因子、分析评价区域生态环境等都具有重要的现实意义。植被覆盖度的监测方法分为地面实测和遥感监测两类。由于植被覆盖度具有显著的时空分异特性,大范围内森林植被覆盖的高密度地面动态实测的方法是不现实的,因而遥感监测已成为区域植被覆盖动态监测的主要手段〔6,7〕。遥感估算植被覆盖度的方法大致可归纳为两种:经验模型法和植被指数转换法。其中,经验模型法主要是通过建立实测植被覆盖度数据与植被指数之间的经验模型来求取大面积植被覆盖。由于经验模型法基于特定区域的实测数据,在局地使用效果可能较好,但在大面积推广应用时则受到限制。植被指数转换法则是在对光谱信号进行分析的基础上,通过建立植被指数和植被覆盖度的转换关系来直接估算植被覆盖度。由于植被指数转换法不依赖于植被覆盖度的实测数据,因而实际工作中使用较多。目前比较广泛使用的是Gutman

216 遥　感　技　术　与　应　用 第20卷