利用中低分辨率DEM提取坡耕地坡度信息的误差分析图文精.docx

1、利用中低分辨率DEM提取坡耕地坡度信息的误差分析图文精利用中低分辨率DEM 提取坡耕地坡度信息的误差分析王莉1, 2, 3, 贺秀斌1, 23, 鲍玉海1, 2, 南宏伟1, 2, 3, 刘艳锋1, 2, 3(1. 中国科学院成都山地灾害与环境研究所山地环境演变与调控重点实验室, 四川成都610041; 2. 中国科学院水利部成都山地灾害与环境研究所, 四川成都610041; 3. 中国科学院研究生院, 北京100039摘要利用数字高程模型(DE M 实现坡耕地数据的自动提取和侵蚀状况调查, 并与实地勘察及相关研究结果进行对比。结果表明, 三峡库区平行岭谷地貌、出露的侏罗纪紫色泥(页 岩和砂岩

2、互层结构及人为活动是小尺度坡耕地地形的关键控制因素。关键词坡耕地; DE M 分辨率; 土地利用; 水土流失中图分类号S127文献标识码A 文章编号0517-6611(2010 07-03644-03Error Ana lysis on Gradient I nforma tion of Slopi ng Land D er i ved from the L ow 2resoluti on D E MW ANG L i et a l (Key Laboratory ofMountain Environment and Regulation, I nstitute ofMountain Haza

3、rds and Environment, Chinese Acad 2emy of Sciences, Chengdu, Sichuan 610041Abstract By using DE M , sl ope inf or mati on was extracted automatically and s oil erosi on situation was investigated . Then compared with field survey and related research . The results showed that parallel ridge valley i

4、n Three Gorges Reservoir A rea, the interbedded structures f or med by Jurassic purp le clay (page rocks and human activities are the key controlling factors of s mall 2scale sl op ing terrain . Key words Slop ing land; Res olution of DE M; Land use; Soil er osion基金项目中国科学院西部行动计划项目(KZ CX22XB2207201 ;

5、 国家科技支撑计划项目(2008BAD98B04 。作者简介王莉(1985- , 女, 陕西西安人, 硕士研究生, 研究方向:土壤侵蚀与水土保持。3通讯作者, 研究员, 博士生导师, E 2mail:xiubinhi m de . ac . cn 。收稿日期2009211217坡耕地研究是水土流失治理及水土保持规划的核心内容1, 掌握坡耕地的数量、土流失具有重要意义。数字高程模型(D igital on 2el, DE M , S 三维向量的有限序列, , 数学表达式为z =f (x, y , 其中, z 表示高程, x, y 表示水平坐标, 也可用经纬度来描述。由于DE M 描述的是地面高程

6、信息, 在测绘、水文、气象、地貌、地质、土壤、工程建设、通讯、军事等国民经济和国防建设等领域有着广泛应用2。诸多生态环境建设工程项目, 尤其是水土流失监测与水土保持规划工作的开展都需要DE M 数据支持3, 利用DE M 为信息源自动提取地面坡度, 建立坡耕地资源数据库, 已成为最重要的技术方法, 且得到了广泛应用4。笔者利用数字高程模型实现对忠县坡耕地的自动提取, 结合该区土壤侵蚀图对坡耕地侵蚀程度进行分级, 并与实地勘察及相关研究结果进行对比, 分析了利用中低分辨率DE M 提取坡耕地坡度信息误差的来源及其对水土流失调查结果的影响。1研究区概况与研究方法1. 1研究区概况忠县位于重庆市中部

7、、三峡库区腹心地带, 30033035N, 107310814E, 东西长66. 45k m, 南北宽60. 15k m, 总面积2176k m 2。境内溪河纵横交错, 地形复杂, 地貌由金华山、方斗山、猫耳山3个背斜和其间的拔山、忠州2个向斜构成。地处暖湿亚热带东南季风区的忠县, 雨量充沛, 年降雨量约1200mm, 相对湿度达80%。坡耕地分布广泛, 土壤类型主要为紫色土。1. 2研究方法采用基于忠县150000地形图数据源的三角网法生成DE M, 即首先将数字化地形图生成拓扑关系完S (N TI N 转化为DE M (图1 。 25m 。图1忠县数字高程模型F i g . 1 D E M

8、 of Zhongxi a n County图2忠县坡度分级图F i g . 2M ap of slope cl a ssi f i ca ti on i n Zhongx i a n County在已经生成的DE M 基础上, 利用GI S 软件空间分析功能中的表面分析功能生成坡度分级(图2 。目前, 利用DE M 提取地面坡度的算法可归纳为5种3:四块法、空间矢量分析法、拟合平面法、拟合曲面法、直接解法。其中, 前3种方法是为解求地面平均坡度而设计的, 后2种方法是为解求地面最大坡度而设计的。在A rc Map 中, 坡度的计算采用拟合责任编辑陈秀晨责任校对张士敏安徽农业科学, Journ

9、al of Anhui Agri . Sci . 2010, 38 (7 :3644-3646曲面法, 即采用3像元3像元的窗口, 每个窗口中心为1个高程点, 点e 的坡度计算方法如下:坡度=A tan(dz /dx 2+(dz/dy 257. 29587式中, dz /dx =(c +2f +i -(a +2d +g /(8x _cell _size ,dz /dy =(g +2h +i -(a +2b +c /(8y _cell _size a b c d e f ghi按照国家通用标准, 将地面坡度划分为6级, 按此分级生成坡度分级图, 并把坡度分级图由栅格格式转为矢量格 式。在此基础上

10、, 应用GI S 强大的空间分析功能结合忠县土地利用图(图3 和土壤侵蚀图(图4 提取坡耕地信息及其水土流失信息, 分析计算机提取结果与实际结果出现偏差的原因。 忠县土地利用图和土壤侵蚀图利用遥感影像解译获得。图3忠县土地利用图F i g . 3M ap of l and use i n Zhongx i a n Coun ty图4忠县土壤侵蚀图F i g . 4Map of so il erosi on i n Zhongxi a n County2结果与分析2. 1坡耕地信息提取结果将矢量化的坡度分级与该区域数字化土地利用进行综合分析, 利用“查询”功能将忠县各坡度级的耕地面积提取出来并统

11、计, 得到各个坡度水田和旱地的面积和比例(表1 。由表1可知, 大于25的水田和旱地的面积分别为30. 50和104. 52k m 2, 比例分别为14. 71%和11. 56%。我国实施“退耕还林”政策后, 25以上坡耕地的面积已大量减少。2. 2坡耕地侵蚀分级结果综合分析忠县数字化土地利用和土壤侵蚀, 耕地中水田大多为微度侵蚀, 其面积占总水田面积的87. 03%; 而旱地的侵蚀程度则较严重, 中、强度侵蚀的旱地占总旱地面积的89. 59%(表2 。HE 等6通过137Cs核素示踪法对基于坡面的土壤侵蚀进行研究, 结果表明该区20的坡耕地土壤平均侵蚀速率为2415t/(k m 2a 5,

12、根据我国1997年5月实施的土壤侵蚀分类分级标准(表3 , 属于轻度侵蚀。而基于DE M 获取的忠县1525的坡耕地侵蚀结果显示, 中度以上侵蚀程度的耕地面积占了90%, 与上述研究结果存在差异。表1不同坡度水田、旱地面积及所占比例Table 1Area and proporti on of paddy f i eld and dry l and i n di fferen tslopes坡度Sl ope 水田Paddy field面积km2比例%Pr oporti on 旱地D ryland面积k m2比例%Pr oporti on 981. 12. 88 1. 422Erosi on le

13、vel of dry l and i n d i fferent slopes坡度Sl ope 面积k m 2微度侵蚀M icr o 2degree 轻度侵蚀L ight degree 中度侵蚀Moderate强度侵蚀I ntensity 极强度侵蚀Str ongest 350. 0800. 016. 006. 750. 028584表3土壤侵蚀强度分级标准Table 3C l assi f i cati on st andard of so il erosi on i n ten sity级别Level平均侵蚀模数t/(k m 2a Mean er osi on modulus 平均流失厚度

14、mm /aMean l oss thickness 微度M icr o 2degree200, 600, 10001500011. 12. 3误差来源分析利用中低分辨率DE M 自动提取的坡耕地信息, 与实际调查结果之间存在较大差距, 从而影响了其土壤侵蚀程度分级的结果, 该区特殊的地貌特征和人类活动对坡耕地的改造是造成误差的主要原因。2. 3. 1地貌特征对坡度提取的影响。忠县的地貌是由金华山、方斗山、猫耳山3个背斜和其间的拔山、忠州2个向斜构成, 属于典型的平行岭谷地貌, 土壤则以出露的侏罗纪紫色泥(页 岩与砂岩为主, 共同作用形成一种阶梯状的地貌特征。这一特殊的地貌结构使得该区的地形颇为

15、复杂, 由图5可知, 整个大的坡面坡度为18, 由于特殊的土壤和地形构造, 在这一大的坡面上形成了多个坡度较小的坡, 这些小坡将整个坡面划分成许多段, 相邻的坡面之间存在一个坎, 从而形成一层层阶梯状坡地, 而每段坡面的坡度一般都小于546338卷7期王莉等利用中低分辨率DE M 提取坡耕地坡度信息的误差分析8, 坡长也较短, 一般在5m 左右。图5忠县平行岭谷地貌特征F i g . 5Geomorphologi c character isti cs of the para llel r i dge va l 2ley i n Zhongx i a n Coun ty中低分辨率的DE M 只

16、能从宏观的角度描述地形的起伏, 而对于坡面上较小的坡度变化却很难体现, 进而影响了坡耕地信息提取的准确性, 坡度信息提取的误差则进一步影响了对坡耕地水土流失现状的评价。DE M 对坡度的提取有着较大影响, DE 率的降低, 坡度衰减, 6。提高DE M 的分辨率, 从微观的角度去描述这一地形特征, 对准确提取坡度信息及评价坡耕地土壤侵蚀状况有重要意义。此外, 人类活动对坡耕地微地形的改造也是影响坡度信息提取的重要因素。2. 3. 2坡度信息误差对坡耕地水土流失现状评价的影响。基于DE M 的坡耕地信息提取是在宏观尺度上进行的, 尤其是在其分辨率较低时, 对于起伏变化较大的坡地难以描述, 而人为

17、活动对坡耕地的改造则更加大了坡耕地自动提取工作的难度, 由于农民对坡耕地的改造, 将流失的土壤又返回到坡面上, 从而一定程度上减轻了水土流失程度。2种因素共同影响了坡耕地坡度信息的提取结果, 而坡度信息的误差则进一步影响了对坡耕地水土流失现状的评价结果。因此, 在评价坡耕地坡度对土壤侵蚀程度的影响时, 应考虑并研究人为因素对地形改造的影响。3结论与讨论应用DE M 对坡度信息的提取是区域宏观坡耕地数据和水土流失现状评价的基础, 不同比例与分辨率DE M 在地形信息载量上存在明显差异, 对坡度影响较大。该研究利用忠县150000地形图生成分辨率为25m 的DE M, 基于该分辨率DE M 自动提

18、取的坡耕地信息在评价其水土流失现状时与实地勘察和相关研究结果存在一定差异, 其原因主要为该地区地貌特殊, 平行岭谷地貌和出露的侏罗纪紫色泥(页 岩与砂岩互层结构共同作用形成的地形较为复杂, 造成该区土地利用多以小地块为主的镶嵌结构, 而分辨率较低的DE M 对地形起伏的描述误差较大。另外, 人为修建的梯田(水田 、坎田以及垄沟整地措施对坡耕地微观地形改变的影响较大。因此, 在坡耕地水土流失现状评价中, 单从中低分辨率DE M 获取坡度信息来判断水土流失程度, 将会产生较大误差。采用不同分辨率DE M 和不同类型DE M 对坡度信息的提取有重要影响。DE M 分辨率越高, 对与坡度有关的地貌和水

19、文参数的描述越精确, 中低分辨率坡度可以通过变换, , 实现了2, 对较低分辨率, 该方法仍不能分型变换是一个值得注意的研究方法, 但有待进一步的验证。针对DE M 生成方法, 目前一种专业化的DE M 生产软件ANUDE M 越来越广泛地被应用, 其差值生成的DE M 比基于TI N 方法建立的DE M 能更加精确地反映重要水文地貌特征和地面起伏7。此外, 加强实地勘察力度, 并与计算机自动提取技术相结合, 才是正确评价坡耕地水土流失的方法。参考文献1李晓燕, 李壁成, 安韶山. 基于GI S 的小流域坡耕地调查方法探讨J .水土保持通报, 2001, 21(5 :37-38.2郭鹏, 孙艳

20、玲, 白洁, 等. 数字高程模型在坡耕地调查中的应用J .国土资源遥感, 2003(4 :59-61.3汤国安, 杨勤科, 张勇, 等. 不同比例尺DE M 提取地面坡度的精度研究以在黄土丘陵沟壑区的试验为例J.水土保持通报, 2001, 21(1 :53-56.4汤国安, 赵牡丹, 李天文, 等. DE M 提取黄土高原地面坡度的不确定性J.地理学报, 2003, 58(6 :824-830.5HE X B, X U Y B, Z HANG X B. Traditi onal far ming syste m f or s oil conser 2vati on on sl ope far

21、mland in s outh western China J .S oil &Tillage Re 2search, 2007, 94:193-200.6杨勤科, 贾大韦, 李锐, 等. 基于DE M 的坡度研究现状与展望J .水土保持通报, 2007, 27(1 :146-149. 7杨勤科,MC V I A C AR TI M R, 李领涛, 等. ANUDE M -专业化数字高程模型插值算法及其特点J.干旱地区农业研究, 2006, 24(3 :36-40.(上接第3467页16米琴, 曹长年. 齐墩果酸对红细胞免疫和细胞毒性的影响J .天水师范学院学报, 2004, 24(5 :18-19.17青海省革命委员会卫生局. 青海省药品标准M.青海省革命委员会卫生局, 1977.18任延明, 袁明, 任世存. 湿生扁蕾提取物对小鼠实验性肝损伤的保护作用J.青海医学院学报, 2008, 29(4 :253-255.19曹长年, 米琴, 赵宙兴, 等. 藏药湿生扁蕾有效成分的研究J.青海农林科技, 2004(3 :12.20刘丽莎, 姬可平. 秦艽种子发芽特性的研究J .中草药, 2002, 33(3 :269-271.21高宏, 李府垠, 秦临喜, 等. 湿生扁蕾种子处理和田间覆盖试验J.甘肃农业大学学报, 2006(3 :40-43.6463安徽农业科学2010 年